Loading [MathJax]/jax/output/SVG/fonts/TeX/fontdata.js
Research article Special Issues

Three-Dimensional Heterogeneous Material Microstructure Reconstruction from Limited Morphological Information via Stochastic Optimization

  • We present a general computational framework that enables one to generate realistic 3D microstructure models of heterogeneous materials from limited morphological information via stochastic optimization. In our framework, the 3D material microstructure is represented as a 3D array, whose entries indicate the local state of that voxel. The limited structural data obtained in various experiments correspond to different mathematical transformations of the 3D array. Reconstructing the 3D material structure from such limited data is formulated as an inverse problem, originally proposed by Yeong and Torquato [Phys. Rev. E 57, 495 (1998)], which is solved using the simulated annealing procedure. The utility, versatility and robustness of our general framework are illustrated by reconstructing a polycrystalline microstructure from 2D EBSD micrographs and a binary metallic alloy from limited angle projections. Our framework can be also applied in the reconstructions based on small-angle x-ray scattering (SAXS) data and has ramifications in 4D materials science (e.g., charactering structural evolution over time).

    Citation: Jiao Yang. Three-Dimensional Heterogeneous Material Microstructure Reconstruction from Limited Morphological Information via Stochastic Optimization[J]. AIMS Materials Science, 2014, 1(1): 28-40. doi: 10.3934/matersci.2014.1.28

    Related Papers:

    [1] Jianxin Yang, Zhipeng Qiu, Xue-Zhi Li . Global stability of an age-structured cholera model. Mathematical Biosciences and Engineering, 2014, 11(3): 641-665. doi: 10.3934/mbe.2014.11.641
    [2] Mingtao Li, Xin Pei, Juan Zhang, Li Li . Asymptotic analysis of endemic equilibrium to a brucellosis model. Mathematical Biosciences and Engineering, 2019, 16(5): 5836-5850. doi: 10.3934/mbe.2019291
    [3] Linhua Zhou, Meng Fan, Qiang Hou, Zhen Jin, Xiangdong Sun . Transmission dynamics and optimal control of brucellosis in Inner Mongolia of China. Mathematical Biosciences and Engineering, 2018, 15(2): 543-567. doi: 10.3934/mbe.2018025
    [4] Zongmin Yue, Yuanhua Mu, Kekui Yu . Dynamic analysis of sheep Brucellosis model with environmental infection pathways. Mathematical Biosciences and Engineering, 2023, 20(7): 11688-11712. doi: 10.3934/mbe.2023520
    [5] Mingtao Li, Guiquan Sun, Juan Zhang, Zhen Jin, Xiangdong Sun, Youming Wang, Baoxu Huang, Yaohui Zheng . Transmission dynamics and control for a brucellosis model in Hinggan League of Inner Mongolia, China. Mathematical Biosciences and Engineering, 2014, 11(5): 1115-1137. doi: 10.3934/mbe.2014.11.1115
    [6] Zhisheng Shuai, P. van den Driessche . Impact of heterogeneity on the dynamics of an SEIR epidemic model. Mathematical Biosciences and Engineering, 2012, 9(2): 393-411. doi: 10.3934/mbe.2012.9.393
    [7] Yaoyao Qin, Xin Pei, Mingtao Li, Yuzhen Chai . Transmission dynamics of brucellosis with patch model: Shanxi and Hebei Provinces as cases. Mathematical Biosciences and Engineering, 2022, 19(6): 6396-6414. doi: 10.3934/mbe.2022300
    [8] Qiuyi Su, Jianhong Wu . Impact of variability of reproductive ageing and rate on childhood infectious disease prevention and control: insights from stage-structured population models. Mathematical Biosciences and Engineering, 2020, 17(6): 7671-7691. doi: 10.3934/mbe.2020390
    [9] Pengyan Liu, Hong-Xu Li . Global behavior of a multi-group SEIR epidemic model with age structure and spatial diffusion. Mathematical Biosciences and Engineering, 2020, 17(6): 7248-7273. doi: 10.3934/mbe.2020372
    [10] Tsanou Berge, Samuel Bowong, Jean Lubuma, Martin Luther Mann Manyombe . Modeling Ebola Virus Disease transmissions with reservoir in a complex virus life ecology. Mathematical Biosciences and Engineering, 2018, 15(1): 21-56. doi: 10.3934/mbe.2018002
  • We present a general computational framework that enables one to generate realistic 3D microstructure models of heterogeneous materials from limited morphological information via stochastic optimization. In our framework, the 3D material microstructure is represented as a 3D array, whose entries indicate the local state of that voxel. The limited structural data obtained in various experiments correspond to different mathematical transformations of the 3D array. Reconstructing the 3D material structure from such limited data is formulated as an inverse problem, originally proposed by Yeong and Torquato [Phys. Rev. E 57, 495 (1998)], which is solved using the simulated annealing procedure. The utility, versatility and robustness of our general framework are illustrated by reconstructing a polycrystalline microstructure from 2D EBSD micrographs and a binary metallic alloy from limited angle projections. Our framework can be also applied in the reconstructions based on small-angle x-ray scattering (SAXS) data and has ramifications in 4D materials science (e.g., charactering structural evolution over time).



    Brucellosis, a zoonotic disease, is a natural epidemic disease that is not only prevalent among livestock and humans but also widely spread in wild animals [1,2]. There are human and animal brucellosis in most countries of the world, and high-risk areas of epidemics are mainly distributed in developing countries, such as Syria, Jordan, Zambia, Mongolia [3,4,5,6]. It has brought huge losses to the livestock industry worldwide and serious health problems to livestock-related practitioners, there are hundreds of thousands of new cases reported annually [6,7,8]. Infected domestic and wild animals and their excreta are the main source of infection, it is contact and pathogen infection which are the main modes of brucellosis transmission [9]. Human brucellosis is rarely transmitted to susceptible animals, and there is no infection between people reported [10,11]. Therefore, it is the eradication of animal brucellosis that is the only way to solving human health problem, and understanding the mechanisms and risk factors of the spread of brucellosis is one of the first problems that must be solved.

    Brucellosis is transmitted to susceptible individuals mainly by contact with infectious individuals or by sucking pathogens from the environment [8]. In animals, the transmission of brucellosis mainly occurs between sexually mature animals, young animals may also be infected, but generally do not have any clinical manifestations and serological tests are usually non-positive, and the infection between young animals is so little as to be almost invisible [12]. In other words, sexually mature animals are very susceptible to brucellosis, latent infections may be found in young animals, but they are generally resistant [9]. It is important that many infected animals have a longer incubation period and that these animals remain serologically non-positive during this period. That is to say, latent animals may not be infectious and it is almost impossible to be found through detection, which is an important risk factor in the elimination of brucellosis [9]. Therefore, the impact of these mechanisms on brucellosis is worthy to be studied using mathematical models.

    Statistical methods such as descriptive statistics, correlation analysis and time series analysis have been widely used for quantitative assessment of risk management measures for brucellosis, among which is very worthy of concern about the study of bison, elk and livestock brucellosis (see [13,14,15,16,17]). Theoretical studies on the impacts of the transmission mechanisms of brucellosis and applied studies on the assessments of risk management measures have also been studied using the kinetic model (see [18,19,20,21,22]). Especially in recent years, many kinetic models with indirect transmission have been established to analyze brucellosis transmission (see [23,24,25,26,27]). Although there have been many studies on brucellosis, there are still many transmission mechanisms and risk factors that are not considered in existing models. For example, infected individuals may have no infectivity in the early stage, during which time is different; and animals that are not sexually mature are hardly infected with direct and indirect modes [9]. Therefore, in the present work, a multi-stage dynamic model with distributed time delay is proposed involving the above risk factors and general nonlinear incidences. The existence and uniqueness of the endemic equilibrium is analyzed, and the local and global asymptotic stability of equilibria is proved.

    The rest of this manuscript is constructed as follows. The dynamic model with distributed time delay and some preliminary results are given in Section 2. In Section 3, the content of research is the global dynamics of the disease-free equilibrium. The global stability of the equilibrium of persistent infection are analyzed in Section 4. In Section 5, the stability results are further explained by numerical simulation. A summary and further discussion is proposed in Section 6.

    According to the transmission mechanisms of animal brucellosis, we classify animal population into three compartments: The young susceptible compartment $ S_1(t) $, the adult susceptible compartment $ S_2(t) $ (sexually mature), the infected compartment $ I(t) $. Similar to the definition in literature [25], $ B(t) $ can be defined as the concentration or number of pathogens in the environment. Some explanations and assumptions about the kinetic model are listed as follows. (Ⅰ) Since brucellosis can cause abortion and significantly reduce the survival rate of young animals, we suppose that infected animals have no birth rate. (Ⅱ) The supplementary rate of young susceptible animal population $ S_1(t) $ is mainly derived from birth and import, it is assumed that the supplementary rate is $ A+bS_2(t) $. (Ⅲ) Usually, there are two different stages for infected animals, which are no infection force in the early stage and are infectious carriers when infected animals begin to shed brucella. Therefore, we make use of a delay $ \tau $ to describe the time from a susceptible individual to a infectious individual, and since the length of this time is varied, thus the delay $ \tau $ is a distributed parameter in the interval $ [0, h] $, where $ h $ is the maximum value of the delay. As a result, the susceptible class $ S_2(t) $ is reduced due to infection at the rate $ S_2 g(B)+S_2\int_0^{h} \rho(\tau)f(I(t-\tau))e^{-(\mu_2+c)\tau}d\tau $. (Ⅳ) Since the infected animal sheds the pathogen after the time $ \tau $, then the pathogen shedding rate of the infected animal is $ \int_0^{h}\rho(\tau)h(I(t-\tau))e^{-(\mu_2+c)\tau}d\tau $. (Ⅴ) In animal breeding, only the basic ewes breed newborns, and adult animals are widely used for slaughter or trade, then the birth rate $ b $ of adult animals is assumed to be less than the elimination rate $ \mu_2 $ of adult animals. Therefore, the modelling for animal brucellosis is given through the following distributed time-delay system:

    $ {dS1dt=A+bS2μ1S1σ1S1,dS2dt=σ1S1S2g(B)S2h0ρ(τ)f(I(tτ))e(μ2+c)τdτμ2S2,dIdt=S2h0ρ(τ)f(I(tτ))e(μ2+c)τdτ+S2g(B)(μ2+c)I,dBdt=h0ρ(τ)h(I(tτ))e(μ2+c)τdτdB. $ (2.1)

    Here, $ b $ is the birth rate of adult animals. $ \mu_1 $ and $ \mu_2 $ are the elimination rates of young and adult animals, respectively. $ \sigma_1 $ is the transfer rate from young susceptible individuals to sexually mature individuals. The culling rate is defined by $ c $ and $ d $ is the decay rate of the pathogen. $ f(I) $ and $ g(B) $ are contact and indirect infection rates, respectively. $ \rho(\tau) $ is a distributed function, it is non-negative, continuous and $ \int_0^{h} \rho(\tau)d\tau = 1 $.

    The initial conditions for system (2.1) are given as follows:

    $ {S1(x)=ϕ1(x),S2(x)=ϕ2(x),I(x)=ϕ3(x),B(x)=ϕ4(x),x[h,0],h>0,ϕ=(ϕ1,ϕ2,ϕ3,ϕ4)C+C. $ (2.2)

    Here, $ C $ denotes the Banach space $ C([-h, 0], \Re^4) $ of continuous functions mapping the interval $ [-h, 0] $ into $ \Re^4 $ with the sup-norm $ \|\phi\| = \sup_{x\in[-h, 0]}|\phi_i(x)|, i = 1, 2, 3, 4 $ for $ \phi\in C $. The nonnegative cone of $ C $ is defined as $ C_+ = C([-h, 0], \Re_+^4) $.

    In order to make epidemiological significance for system (2.1), $ f, g $ and $ h $ are assumed to be second-order continuous differentiable functions and satisfy the following hypotheses:

    $ (H_1) $ $ f(0) = g(0) = h(0) = 0 $ and $ f(I), g(B), h(I) > 0 $ for $ I, B > 0 $;

    $ (H_2) $ $ f'(I), g'(B) > 0 $ and $ f''(I), g''(B)\leq0 $ for $ B, I\geq0 $;

    $ (H_3) $ $ h'(I) > 0 $ and $ h''(I)\leq0 $ for $ I\geq0 $.

    The function $ f $ may be $ \beta I^p $ or saturation incidences $ k\ln(1+\frac{\lambda I}{k}) $ and $ \frac{\beta I^p}{1+kI^p} $ with constants $ \beta, p, \lambda, k > 0 $ [28,29]. The pathogen infection rate $ g $ may be $ \frac{\lambda B}{1+TB} $ or $ \eta(1-e^{-\alpha B}) $ with constants $ \lambda, \eta, \alpha > 0 $ and $ T\geq0 $ [30]. The function $ h $ may be $ kI $ with constants $ k > 0 $ [31].

    It is easy to verify that the system (2.1) always has a disease-free equilibrium $ E_0 = (S_1^0, S_2^0, 0, 0) $, where

    $ S_1^0 = \frac{A\mu_2}{\mu_2(\mu_1+\sigma_1)-b\sigma_1}, \hskip 0.5cm S_2^0 = \frac{A\sigma_1}{\mu_2(\mu_1+\sigma_1)-b\sigma_1}. $

    The reproduction number of sytem (2.1) is given by the following expression:

    $ R0=nS02fI(0)μ2+c+nS02gB(0)hI(0)d(μ2+c)=R01+R02. $ (2.3)

    where $ n = \int_0^{h}\rho(\tau)e^{-(\mu_2+c)\tau}d\tau $. $ e^{-(\mu_2+c)\tau} $ denotes the probability of survival from a newly infected susceptible individual to an infectious individual. $ n $ is the total survival rate of infected individuals. $ f_I(0) $ is the infection rate by an infectious individual. $ g_B(0) $ is the infection rate of brucella. $ \frac{1}{\mu_2+c} $ is the average life span of infected individual. $ \frac{h_I(0)}{d} $ represents the total number of brucella shed by infectious individuals. Therefore, based on the explanation in the literature [26], $ R_0 $ relies on direct and indirect infection and can be divided into two parts $ R_{01} $ and $ R_{02} $.

    In the next section, we establish that the solution of system (2.1) is non-negative and bounded and analyze the uniqueness of the equilibrium point of persistent infection. The following qualitative results indicate the nonnegativity and boundedness of the solution of system (2.1) with initial conditions (2.2).

    Theorem 2.1. $ (S_1(t), S_2(t), I(t), B(t)) $ is the solution of system (2.1) with the initial conditions (2.2), then $ S_1(t), S_2(t), I(t), B(t) $ are nonnegative and ultimately bounded.

    Proof. Since $ (\dot{S_1}(t)+\dot{S_2}(t)) = A > 0 $ for $ t\in[0, +\infty) $ when $ S_1(t) = S_2(t) = 0 $, it implies that $ S_1(t)+S_2(t)\geq0 $ for all $ t\in[0, +\infty) $. Therefore, if $ S_1(t) = 0 $, then $ S_2(t)\geq0 $ and $ \dot{S_1}(t) = A+bS_2 > 0 $, so that $ S_1(t)\geq0 $. We denote

    $ a1(t)=h0ρ(τ)f(I(tτ))e(μ2+c)τdτ+g(B)+μ2,a2(t)=S2(t)(h0ρ(τ)f(I(tτ))e(μ2+c)τdτ+g(B)), $

    and

    $ a_3(t) = \int_0^{h}\rho(\tau)h(I(t-\tau))e^{-(\mu_2+c)\tau}d\tau. \quad\quad\quad\quad\quad\hskip 0.3cm $

    From the last three equations in (2.1), it can conclude that

    $ S2(t)=S2(0)et0a1(η)dη+σ1t0etηa1(ξ)dξS1(η)dη0,I(t)=I(0)e(μ2+c)t+t0e(μ2+c)(tη)a2(η)dη0, $

    and

    $ B(t) = B(0)e^{-dt}+\int_0^te^{-d(t-\eta)}a_3(\eta)d\eta\geq0\quad\quad\quad\quad\quad\hskip 0.1cm $

    for all $ t\geq0 $. Thus, $ S_1(t), S_2(t), I(t), B(t)\geq0 $ for $ t\geq0 $.

    We now analyze the boundedness of the solution of system (2.1). In fact,

    $ (\dot{S_1}(t)+\dot{S_2}(t)+\dot{I}(t))\leq A-\mu_1S_1-(\mu_2-b)S_2-(\mu_2+c)I \quad\quad\quad\quad\quad\quad\hskip 0.2cm $
    $ \leq A-\mu(S_1(t)+S_2(t)+I(t)), $

    where $ \mu = \min\{\mu_1, \mu_2-b\} $.

    Hence, $ \limsup_{t\rightarrow\infty}(S_1(t)+S_2(t)+I(t))\leq \frac{A}{\mu} $, it implies that $ \limsup_{t\rightarrow\infty}I(t)\leq \frac{A}{\mu} $ and $ \limsup_{t\rightarrow\infty}B(t)\leq \frac{n}{d}h(\frac{A}{\mu}) $. Therefore, $ S_1(t), S_2(t), I(t), B(t) $ are ultimately bounded.

    Since the solution of system (2.1) is nonnegative and ultimately bounded, the set

    $ \Omega = \{(S_1(\cdot), S_2(\cdot), I(\cdot), B(\cdot))\in C_+: \|S_1+S_2+I\|\leq \frac{A}{\mu}, \|B\|\leq \frac{n}{d}h(\frac{A}{\mu})\}, $

    is positively invariant for system (2.1).

    In order to analyze the uniqueness of the positive solution of system (2.1), the following lemma is given:

    Lemma 2.1. Assume that conditions $ (H_1)-(H_{3}) $ are satisfied, the functions $ \frac{f(I)}{I} $, $ \frac{h(I)}{I} $, $ \frac{g(B)}{B} $ and $ \frac{g(h(I))}{I} $ are monotonic decreasing for $ I, B > 0 $.

    Proof. Since $ f''(I)\leq0 $, it shows that $ f'(I) $ is monotonic decreasing, it follows that

    $ \frac{f(I)}{I} = \frac{f(I)-f(0)}{I-0} = f'(\xi_1)\geq f'(I),\hskip 0.5cm \xi_1\in(0, I), $

    and

    $ (\frac{f(I)}{I})' = \frac{f'(I)I-f(I)}{I^2}\leq0. $

    That is to say, $ \frac{f(I)}{I} $ is a monotonically decreasing function. On the basis of the above method, it can show that $ \frac{h(I)}{I} $ and $ \frac{g(B)}{B} $ are also monotonic decreasing. Noting that

    $ g(h(I))I=g(h(I))g(h(0))h(I)h(0)h(I)I=g(h(ξ2))h(ξ3)g(h(I))h(I),ξ2,ξ3(0,I), $

    It can deduce that

    $ (\frac{g(h(I))}{I})' = \frac{g'(h(I))h'(I)I-g(h(I))}{I^2}\leq0. $

    Therefore, $ \frac{g(h(I))}{I} $ is also monotonic decreasing.

    For system (2.1), the endemic equilibrium $ E_* = (S_1^*, S_2^*, I^*, B^*) $ can be derived from the following algebraic equations:

    $ {A+bS2=μ1S1+σ1S1,σ1S1=S2(nf(I)+g(B))+μ2S2,S2(nf(I)+g(B))=(μ2+c)I,nh(I)=dB. $ (2.4)

    By direct calculation, it can be written as

    $ {S1=A+bS2μ1+σ1,m(S2)=S2(nf(I)+g(B)),S2(nf(I)+g(B))=(μ2+c)I,B=H(I). $

    where $ m(S_2) = \frac{A\sigma_1-(\mu_2(\mu_1+\sigma_1)-b\sigma_1)S_2}{\mu_1+\sigma_1} $ and $ H(I) = \frac{n}{d}h(I) $.

    Let us define

    $ F1(S2,I)m(S2)S2(nf(I)+g(H(I))),F2(S2,I)S2(nf(I)+g(H(I)))(μ2+c)I. $

    Similar to the method of analysis in literature [26], using Lemma 2.1, the following result can be summarized:

    Theorem 2.2. Assume that conditions $ (H_1)-(H_3) $ hold. Then there is a unique positive solution $ E_* = (S_1^*, S_2^*, I^*, B^*) $ of system (2.1) if and only if $ R_0 > 1 $.

    In this section, we show that the global stability of the equilibrium point $ E_0 $ of system (2.1) is independent of the initial value. The following conclusions are first obtained.

    Lemma 3.1. Assume that conditions $ (H_1)-(H_{3}) $ hold. The disease-free equilibrium $ E_0 $ is locally asymptotically stable if $ R_0\leq1 $ and is unstable if $ R_0 > 1 $.

    Proof. The characteristic equation at $ E_0 $ is

    $ |λ+μ1+σ1b00σ1λ+μ2S02fI(0)Γ(λ)S02gB(0)00λ+μ2+cΓ(λ)S02fI(0)S02gB(0)00hI(0)Γ(λ)λ+d|=0, $ (3.1)

    where

    $ \Gamma(\lambda) = \int_0^{h}\rho(\tau)e^{-(\mu_2+c+\lambda)\tau}d\tau. $

    It follows from (3.1) that

    $ (λ+μ1+σ1)(λ+μ2)((λ+d)(λ+μ2+cΓ(λ)S02fI(0))Γ(λ)S02gB(0)hI(0))=bσ1((λ+d)(λ+μ2+cΓ(λ)S02fI(0))Γ(λ)S02gB(0)hI(0)). $

    That is,

    $ (\lambda+\mu_1+\sigma_1)(\lambda+\mu_2)(\lambda+d)(\lambda+\mu_2+c) -b\sigma_1(\lambda+d)(\lambda+\mu_2+c) $
    $ \quad\hskip 0.1cm = ((\lambda+\mu_1+\sigma_1)(\lambda+\mu_2)-b\sigma_1) ((\lambda+d)\Gamma(\lambda)S_2^0f_I(0)+\Gamma(\lambda)S_2^0g_B(0)h_I(0)). $

    It follows that

    $ H_1(\lambda) ((\lambda+d)(\lambda+\mu_2+c)-(\lambda+d)\Gamma(\lambda)S_2^0f_I(0)-\Gamma(\lambda)S_2^0g_B(0)h_I(0)) = 0, $

    where $ H_1(\lambda) = (\lambda+\mu_1+\sigma_1)(\lambda+\mu_2)-b\sigma_1 $. Since $ \mu_1+\sigma_1+\mu_2 > 0 $ and $ \mu_2(\mu_1+\sigma_1)-b\sigma_1 > 0 $, then $ H_1(\lambda) $ consists of two roots which are negative real parts. Therefore, we only analyze the distribution of the roots of the following equation:

    $ H2(λ)=(λ+d)(λ+μ2+c)(λ+d)Γ(λ)S02fI(0)Γ(λ)S02gB(0)hI(0)=0. $ (3.2)

    Assume now that $ R_0 > 1 $, then

    $ H_2(0) = d(\mu_2+c)-dnS_2^0f_I(0)-nS_2^0g_B(0)h_I(0)) \lt 0, \hskip 0.5cm H_2(+\infty) = +\infty. $

    Hence $ H_2(\lambda) $ has at least one positive root in $ [0, +\infty) $, then $ E_0 $ is unstable if $ R_0 > 1 $.

    From (3.2), we have

    $ (λ+d)(λ+μ2+c)=Γ(λ)(λS02fI(0)+dS02fI(0)+S02gB(0)hI(0))=(μ2+c)F(λ)n(λR01+dR0), $

    or

    $ (λ+d)(λμ2+c+1)=R0F(λ)n(λR01R0+d). $ (3.3)

    Next, considering the case $ R_0\leq1 $. If $ \lambda = x+yi $ is a solution of (3.3), one shows that $ x < 0 $. Otherwise, $ x\geq 0 $ implies

    $ \left| \lambda+d\right| \gt \left| \lambda\frac{R_{01}}{R_0}+d\right| , \hskip 0.2cm \left| \frac{\lambda}{\mu_2+c}+1\right| \gt 1, \hskip 0.2cm \left| \frac{R_0F(\lambda)}{n}\right| \leq 1, $

    and thus

    $ \left|(\lambda+d)(\frac{\lambda}{\mu_2+c}+1)\right| \gt \left|\frac{R_0F(\lambda)}{n}(\lambda \frac{R_{01}}{R_0}+d)\right|, $

    this is a contradiction to (3.3). Therefore, all roots of Eq (3.3) have no zero and positive real parts when $ R_0\leq1 $, this shows that $ E_0 $ is locally asymptotically stable.

    Theorem 3.1. Assume that conditions $ (H_1)-(H_{3}) $ are established. If $ R_0\leq1 $, the disease-free equilibrium $ E_0 = (S_1^0, S_2^0, 0, 0) $ of system (2.1) is globally asymptotically stable.

    Proof. Since the functions $ \frac{f(I)}{I} $, $ \frac{g(B)}{B} $ and $ \frac{h(I)}{I} $ are decreasing, then we have

    $ nS2f(I)(μ2+c)IlimI0+nS02f(I)(μ2+c)I=nS02fI(0)μ2+cb1,S2g(B)dBlimB0+S02g(B)dB=S02gB(0)db2,nh(I)(μ2+c)IlimI0+nh(I)(μ2+c)I=nhI(0)μ2+c. $

    Define

    $ J = \left(10nhI(0)μ2+c1 \right), \hskip 0.5cm (a_1, a_2) = (b_1, b_2)J^{-1}. $

    We find $ a_1 = R_0 $ and define a Lyapunov functional $ L $ as follows:

    $ L(t) = L_1(t)+L_2(t)+L_3(t), $

    where

    $ L1=R0(S1S01S01lnS1S01+S2S02S02lnS2S02+I)+a2B,L2=R0S02h0Ψ(τ)f(I(tτ))dτ,L3=a2h0Ψ(τ)h(I(tτ))dτ, $

    and

    $ \Psi(\tau) = \int_\tau^{h}\varphi(s)ds, \hskip 0.5cm \varphi(s) = \rho(s)e^{-(\mu_2+c)s}. $

    Then the derivative of $ L_1 $ along the positive solutions of system (2.1) is

    $ dL1dt=R0(1S01S1)dS1dt+R0(1S02S2)dS2dt+R0dIdt+a2dBdt=R0(2A+bS02+σ1S01μ1S1(μ2b)S2S01AS1bS01S2S1σ1S02S1S2)+R0(S02g(B)+S02h0φ(τ)f(I(tτ))dτ(μ2+c)I)+a2(h0φ(τ)h(I(tτ))dτdB). $ (3.4)

    Calculating the derivative of $ L_2(t) $ along the solutions of system (2.1), one obtains

    $ dL2dt=R0S02h0Ψ(τ)df(I(tτ))dtdτ=R0S02h0Ψ(τ)df(I(tτ))dτdτ=R0S02(Ψ(τ)f(I(tτ)))|h0+R0S02h0dΨ(τ)dτf(I(tτ))dτ=R0S02nf(I(t))R0S02h0φ(τ)f(I(tτ))dτ. $ (3.5)

    Similar to the above-used method, it can obtain that

    $ \frac{dL_3}{dt} = a_2nh(I(t))-a_2\int_0^{h}\varphi(\tau)h(I(t-\tau))d\tau. $ (3.6)

    Combining the Eqs (3.4), (3.5) and (3.6), it follows that

    $ dLdt=dL1dt+dL2dt+dL3dt=R0μ1S01(2S1S01S01S1)+R0bS02(2S1S02S01S2S01S2S1S02)+R0(μ2b)S02(3S2S02S1S02S01S2S01S1)+R0(nS02f(I)(μ2+c)I,S02g(B)dB)((μ2+c)I,dB)T(R0,a2)(10nh(I)(μ2+c)I1)((μ2+c)I,dB)TR0(nS02fI(0)(μ2+c),S02gB(0)d)((μ2+c)I,dB)T(R0,a2)(10nhI(0)μ2+c1)((μ2+c)I,dB)T=(R01)(b1,b2)((μ2+c)I,dB)T0. $

    Similar to the analysis of Theorem 1 in the literature [32], it follows from Lemma 3.1 that the disease-free equilibrium $ E_0 $ is globally asymptotically stable by LaSalle's Invariance Principle [33].

    According to Lemma 3.1, the disease-free steady state solution $ E_0 $ is unstable when $ R_0 > 1 $. Using Theorem 4.2 in [34], the uniform permanence of system (2.1) can be proven, the process is ignored here. In this following section, by constructing a Lyapunov functional, the global stability of the equilibrium of persistent infection is proved. We first analyze its local stability.

    Lemma 4.1. Assume that conditions $ (H_1)-(H_{3}) $ are ture, If $ R_0 > 1 $, the endemic equilibrium $ E_* $ of system (2.1) exists and is locally asymptotically stable.

    Proof. For system (2.1), the characteristic equation at $ E_* $ is

    $ \left|λ+μ1+σ1b00σ1H3(λ)Γ(λ)S2fI(I)S2gB(B)0nf(I)g(B)H4(λ)S2gB(B)00Γ(λ)hI(I)λ+d\right| = 0, $

    where $ H_3(\lambda) = \lambda+nf(I^*)+g(B^*)+\mu_2 $, $ H_4(\lambda) = \lambda+\mu_2+c-\Gamma(\lambda)S_2^*f_I(I^*) $. By simple calculation, one obtains

    $ (λ+d)(λ+μ2+c)Φ1(λ)=Φ2(λ)Γ(λ)((λ+d)S2fI(I)+hI(I)S2gB(B))=Φ2(λ)Γ(λ)nΛn(λS2fI(I)Λ+d), $

    or

    $ (λ+d)(λμ2+c+1)Φ1(λ)=Φ2(λ)Γ(λ)nΛnμ2+c(λS2fI(I)Λ+d), $ (4.1)

    where

    $ Φ1(λ)=(λ+μ1+σ1)(λ+μ2+nf(I)+g(B))bσ1,Φ2(λ)=(λ+μ1+σ1)(λ+μ2)bσ1,Λ=S2fI(I)+1dS2hI(I)gB(B) $

    Noting that from Lemma 1

    $ Λnμ2+c=S2fI(I)n+1dS2hI(I)gB(B)nμ2+cS2(f(I)nI+ndh(I)Ig(B)B)μ2+c=S2(f(I)n+g(B))I(μ2+c)=1. $

    Assume $ \lambda = x+yi $ is a solution of (4.1). If $ x\geq0 $, we then have

    $ \left| \lambda+d\right| \gt \left| \lambda\frac{S_2^*f_I(I^*)}{\Lambda}+d\right| , \hskip 0.1cm \left| \frac{\lambda}{\mu_2+c}+1\right| \gt 1, \hskip 0.1cm \left| \frac{F(\lambda)}{n}\right| \leq 1, $

    and

    $ |Φ1(λ)Φ2(λ)|=|(λ+μ1+σ1)(λ+μ2+nf(I)+g(B))bσ1(λ+μ1+σ1)(λ+μ2)bσ1|=|λ+μ2+nf(I)+g(B)bσ1λ+μ1+σ1λ+μ2bσ1λ+μ1+σ1|=|M1+nf(I)+g(B)+M2iM1+M2i|>1, $

    where $ M_1 = x+\mu_2-\frac{b\sigma_1(x+\mu_1+\sigma_1)}{(x+\mu_1+\sigma_1)^2+y^2} > 0 $, $ M_2 = y+\frac{b\sigma_1y}{(x+\mu_1+\sigma_1)^2+y^2} $. So it concludes that

    $ \left|(\lambda+d)(\frac{\lambda}{\mu_2+c}+1)\Phi_1(\lambda)\right| \gt \left| \Phi_2(\lambda)\frac{F(\lambda)}{n}\frac{\Lambda n}{\mu_2+c} (\lambda \frac{S_2^*f_I(I^*)}{\Lambda}+d)\right|, $

    which is a contradiction to (4.1). Therefore, the Eq (4.1) can not have any roots with a nonnegative real part, this implies that $ E_* $ is locally asymptotically stable if $ R_0 > 1 $.

    Theorem 4.1. Assume that conditions $ (H_1)-(H_{3}) $ hold. If $ R_0 > 1 $, the endemic equilibrium $ E_* = (S_1^*, S_2^*, I^*, B^*) $ of system (2.1) is globally asymptotically stable.

    Proof. Define

    $ L1=S1S1S1lnS1S1+S2S2S2lnS2S2+IIIlnII+S2g(B)nh(I)(BBBlnBB). $

    Finding the time derivative of $ L_1 $ along the positive solutions of system (2.1) gives

    $ dL1dt=(1S1S1)dS1dt+(1S2S2)dS2dt+(1II)dIdt+S2g(B)nh(I)(1BB)dBdt=2A+bS2+σ1S1μ1S1(μ2b)S2S1AS1bS1S2S1σ1S2S1S2+S2g(B)+S2h0φ(τ)f(I(tτ))dτ(μ2+c)IIIS2g(B)IIS2h0φ(τ)f(I(tτ))dτ+S2g(B)nh(I)(h0φ(τ)h(I(tτ))dτdB+dB)S2g(B)nh(I)BBh0φ(τ)h(I(tτ))dτ.$

    By the Eq (2.4), we get

    $ dL1dt=μ1S1(2S1S1S1S1)+(μ2b)S2(3S2S2S1S2S1S2S1S1)+bS2(2S1S2S1S2S1S2S1S2)+nS2f(I)(3S1S1S1S2S2S1II)+S2g(B)(3+g(B)g(B)S1S1S1S2S2S1IIS2g(B)IS2g(B)I)+S2h0φ(τ)f(I(tτ))dτIIS2h0φ(τ)f(I(tτ))dτ+S2g(B)(h(I)h(I)BBBh(I)Bh(I)+1)+S2g(B)nh(I)(h0φ(τ)h(I(tτ))dτnh(I)+BBnh(I))S2g(B)nh(I)BBh0φ(τ)h(I(tτ))dτ. $ (4.2)

    Since the function $ \nu(x) = 1-x+\ln x $ is nonpositive for $ x > 0 $ and $ \nu(x) = 0 $ if and only if $ x = 1 $, so we define

    $ L_2 = -S_2^*f(I^*)\int_0^{h}\Psi(\tau)\nu(\frac{f(I(t-\tau))}{f(I^*)})d\tau. $

    A direct calculation shows that

    $ dL2dt=S2f(I)h0Ψ(τ)ddtν(f(I(tτ))f(I))dτ=S2f(I)h0Ψ(τ)ddτν(f(I(tτ))f(I))dτ=S2f(I)(Ψ(τ)ν(f(I(tτ))f(I))h0+h0φ(τ)ν(f(I(tτ))f(I))dτ)=S2f(I)h0φ(τ)ν(f(I(t))f(I))dτ+S2f(I)h0φ(τ)ν(f(I(tτ))f(I))dτ=S2f(I)h0φ(τ)(f(I)f(I)f(I(tτ))f(I)+lnf(I(tτ))f(I))dτ=nS2f(I)f(I)f(I)+S2f(I)h0φ(τ)(f(I(tτ))f(I)+lnf(I(tτ))f(I))dτ. $ (4.3)

    Define

    $ L_3 = -\frac{S_2^*g(B^*)}{n}\int_0^{h}\Psi(\tau)\nu(\frac{h(I(t-\tau))}{h(I^*)})d\tau. $

    Calculating the time derivative of $ L_3(t) $, one obtains

    $ dL3dt=S2g(B)nh0Ψ(τ)ddtν(h(I(tτ))h(I))dτ=S2g(B)nh0Ψ(τ)ddτν(h(I(tτ))h(I))dτ=S2g(B)nh0φ(τ)(h(I)h(I)h(I(tτ))h(I)+lnh(I(tτ))h(I))dτ=S2g(B)h(I)h(I)+S2g(B)nh0φ(τ)(h(I(tτ))h(I)+lnh(I(tτ))h(I))dτ.$ (4.4)

    For system (2.1), the following Lyapunov functional is considered:

    $ L = L_1+L_2+L_3. $

    From (4.2), (4.3) and (4.4), we can get

    $ dLdt=μ1S1(2S1S1S1S1)+(μ2b)S2(3S2S2S1S2S1S2S1S1)+bS2(2S1S2S1S2S1S2S1S2)+nS2f(I)(3+f(I)f(I)S1S1S1S2S2S1IIS2f(I)IS2f(I)I)+S2g(B)(3+g(B)g(B)S1S1S1S2S2S1IIS2g(B)IS2g(B)I)+S2g(B)(h(I)h(I)BBBh(I)Bh(I)+1)+S2f(I)h0φ(τ)F(τ)dτ+S2g(B)nh0φ(τ)H(τ)dτ, $ (4.5)

    where

    $ h0φ(τ)F(τ)dτ=h0φ(τ)(S2f(I)IS2f(I)IS2If(I(tτ))S2f(I)I+lnf(I(tτ))f(I))dτ=h0φ(τ)(ν(S2If(I(tτ))S2f(I)I)ν(S2f(I)IS2f(I)I))dτ=h0φ(τ)ν(S2If(I(tτ))S2f(I)I)dτnν(S2f(I)IS2f(I)I), $ (4.6)

    and

    $ h0φ(τ)H(τ)dτ=h0φ(τ)(h(I)Bh(I)BBh(I(tτ))h(I)B+lnh(I(tτ))h(I))dτ=h0φ(τ)ν(Bh(I(tτ))h(I)B)dτnν(h(I)Bh(I)B). $ (4.7)

    it follows from (4.5), (4.6) and (4.7) that

    $ dLdt=μ1S1(2S1S1S1S1)+(μ2b)S2(3S2S2S1S2S1S2S1S1)+bS2(2S1S2S1S2S1S2S1S2)+nS2f(I)(f(I)f(I)1)(1f(I)If(I)I)+nS2f(I)(ν(S1S1)+ν(S1S2S1S2)+ν(S2f(I)IS2f(I)I)+ν(f(I)If(I)I))+S2g(B)((g(B)g(B)1)(1g(B)Bg(B)B)+(h(I)h(I)1)(1h(I)Ih(I)I))+S2g(B)(ν(S1S1)+ν(S1S2S1S2)+ν(S2g(B)IS2g(B)I)+ν(g(B)Bg(B)B))+S2g(B)(ν(h(I)Bh(I)B)+ν(h(I)Ih(I)I))+S2f(I)(h0φ(τ)ν(S2If(I(tτ))S2f(I)I)dτnν(S2f(I)IS2f(I)I))+S2g(B)n(h0φ(τ)ν(Bh(I(tτ))h(I)B)dτnν(h(I)Bh(I)B))μ1S1(2S1S1S1S1)+(μ2b)S2(3S2S2S1S2S1S2S1S1)+bS2(2S1S2S1S2S1S2S1S2)+nS2f(I)(f(I)f(I)1)(1f(I)If(I)I)+S2g(B)((g(B)g(B)1)(1g(B)Bg(B)B)+(h(I)h(I)1)(1h(I)Ih(I)I))$

    By Lemma 2.1, it can conclude that

    $ \frac{dL}{dt} = \frac{dL_1}{dt}+\frac{dL_2}{dt}+\frac{dL_3}{dt}\leq 0. $

    The equality $ \frac{dL}{dt} = 0 $ suggests that $ \frac{S_1^*}{S_1} = 1 $, $ \frac{S_1S_2^*}{S_1^*S_2} = 1 $, $ \frac{f(I)}{f(I^*)} = 1 $ and $ \frac{g(B)}{g(B^*)} = 1 $, it implies that $ E_* $ is the maximum invariant set of system (2.1) in the set $ \{\frac{dL}{dt} = 0\} $. Using Lemma 4.1, the endemic steady state $ E_* $ is globally asymptotically stable.

    By choosing a specific kernel function and some infection functions, system (2.1) can be evolved into different dynamic models with time delays. In this section, some such examples are used to further illustrate theoretical results.

    Example: Consider the kernel function $ \rho(\tau) = \delta(\tau-\tau_0) $, where $ \delta $ is the Dirac delta function. Then system (2.1) can be rewritten as

    $ {dS1dt=A+bS2μ1S1σ1S1,dS2dt=σ1S1αS2BβS2f(I(tτ0))e(μ2+c)τ0μ2S2,dIdt=αS2B+βS2f(I(tτ0))e(μ2+c)τ0(μ2+c)I,dBdt=kI(tτ0)e(μ2+c)τ0dB. $ (5.1)

    Where $ A, b, \mu_1, \mu_2, \sigma_1, \alpha, \beta, c, k $ and $ d > 0 $, $ g(B) = \alpha B $ and $ h(I) = kI $. If $ f(I) = \beta I $, it is easy to verify that $ f(I), g(B) $ and $ h(I) $ satisfy the assumptions $ (H_1)-(H_3) $. Therefore, the basic reproduction number of system (5.1) can be defined as

    $ R0=βS02e(μ2+c)τ0μ2+c+kαS02e(μ2+c)τ0d(μ2+c). $

    In this case, if $ R_0\leq1 $, the disease-free equilibrium is globally asymptotically stable and the endemic equilibrium is also globally asymptotically stable if $ R_0 > 1 $ (see Figures 1 and 2). As is shown in Figure 2, initial values and time delays can not change dynamic properties of system (2.1) when $ R_0 $ is greater than one or less than or equal to one. However, from Figure 2b, it is easy to find that time delay $ \tau_0 $ has a significant impact on the positive equilibrium of system (5.1).

    Figure 1.  The time series diagrams of system (5.1) illustrate the effect of the initial values. The values of the parameters are $ A = 160, b = 0.04, \sigma_1 = 0.2, \mu_1 = 0.01, \beta = 2.1\times 10^{-5} $, $ \alpha = 8.158\times 10^{-6}, \mu_2 = 0.05 $, $ \tau_0 = 1, k = 2, d = 4 $. Figure 1a illustrates the time series diagrams when $ R_0 = 0.6463 $ $ (c = 0.3) $; Figure 1b shows the time series diagrams when $ R_0 = 1 (c = 0.2) $.
    Figure 2.  The time series diagrams of system (5.1) illustrate the effect of the initial values and the latency delay $ \tau_0 $. The values of the parameters are $ A = 160, b = 0.04, \sigma_1 = 0.2 $, $ \mu_1 = 0.01, \beta = 2.1\times 10^{-5} $, $ \alpha = 8.158\times 10^{-6}, \mu_2 = 0.05, c = 0.1 $, $ k = 2, d = 4 $. Figure 2a illustrates the time series diagrams when $ R_0 = 1.8420 $ $ (\tau_0 = 1) $; Figure 2b shows the time series diagrams when $ R_0 = 1.8420(\tau_0 = 1), 1.5854(\tau_0 = 2) $ and $ 1.3646(\tau_0 = 3) $.

    If $ f(I) = \beta I^2 $ [28], then $ f(I) $ does not satisfy hypothesis $ (H_2) $, that is, $ \frac{f(I)}{I} $ is a monotonically increasing function. From Figure 3, it is easy to see that the system (5.1) appears periodic oscillation behavior under certain conditions.

    Figure 3.  The time series diagrams and phase diagrams of system (5.1) illustrate the effect of the incidence rate $ f(I) $ with the initial value $ (2000, 5000, 2000, 10) $. The values of the parameters are $ A = 330, b = 0.04, \sigma_1 = 0.2 $, $ \mu_1 = 0.01, \beta = 1\times 10^{-9}, \alpha = 8.158\times 10^{-6} $, $ \mu_2 = 0.05, c = 0.02 $, $ k = 2, d = 4, \tau_0 = 1 $.

    In this article, a general $ S_1S_2IB $ dynamics model with distributed time delay for animal brucellosis is formulated. Under the assumptions of general biological significance, the non-negative and boundedness of the solution of system (2.1) is first proved. And then the global dynamics of the steady-state solution of system (2.1) are analyzed by constructing Lyapunov functional, it is found that the dynamic properties of equilibria depend on the basic reproduction number $ R_0 $: If $ R_0\leq 1 $, animal brucellosis will eventually die out regardless of the initial value; and if $ R_0 > 1 $, the spread of animal brucellosis is persistent and it eventually reaches the endemic steady state. These results imply that distributed time delay does not change the dynamic properties of system (2.1) when $ R_0 $ is greater than one or less than or equal to one. Finally, similar to numerical methods in [35,36], the stability results and other dynamic behaviors are further illustrated through numerical simulation, it turns out that the system experiences periodic oscillations if the assumption $ (H_2) $ is not satisfied. In other words, the system (2.1) may exhibit more complex dynamic behaviors if the function $ \frac{f}{I}, \frac{g}{B} $ or $ \frac{h}{I} $ is monotonically increasing. In these cases, the impact of distributed time delay on other dynamical behaviors of system (2.1) is not completely clear, then we leave these for future work.

    This research is partially supported by the National Youth Science Foundation of China (11501528), the National Sciences Foundation of China(11571324) and the Fund for Shanxi "1331KIRT".

    The authors declare that they have no competing interests.

    [1] Christensen RM (1979) Mechanics of Composite Materials.New York: Wiley .
    [2] Nemat-Nasser S, Hori M (1999) Micromechanics: Overall Properties of Heterogeneous Solids.Elsevier Science .
    [3] Torquato S (2002) Random Heterogeneous Materials: Microstructure and Macroscopic Properties.New York: Springer .
    [4] Sahimi M (2003) Heterogeneous Materials I: Linear Transport and Optical Properties, and II: Nonlinear and Breakdown Properties and Atomistic Modeling.New York: Springer .
    [5] Millar DIA (2012) Energetic Materials at Extreme Conditions.Springer Ph.D. Thesis .
    [6] Haymes RC (1971) Introduction to Space Science.New York: John Wiley and Sons Inc .
    [7] Thornton K, Poulsen HF (2008) Three-dimensional materials science: An intersection of three-dimensional reconstructions and simulations.MRS Bull 33: 587.
    [8] Brandon D, Kaplan WD (1999) Microstructural Characterization of Materials.New York: John Wiley and Sons .
    [9] Baruchel J, Bleuet P, Bravin A, et al. (2008) Advances in synchrotron hard x-ray based imaging.C R Physique 9: 624.
    [10] Kinney JH, Nichols MC (1992) X-ray tomographic microscopy (XTM) using synchrotron radiation.Annu Rev Mater Sci 22: 121.
    [11] Kak A, Slaney M (1988) Principles of Computerized Tomographic Imaging.IEEE Press .
    [12] Babout L, Maire E, Buffière JY (2001) Characterisationby X-ray computed tomography of decohesion, porosity growth and coalescence in model metal matrix composites.Acta Mater 49: 2055.
    [13] Borbély A, Csikor FF, Zabler S, et al. (2004) Three-dimensional characterization of the microstructure of a metal-matrix composite by holotomography.Mater Sci Engng A 367: 40.
    [14] Kenesei P, Biermann H, Borbély A (2005) Structure-property relationship in particle reinforced metal-matrix composites based on holotomography.Scr Mater 53: 787.
    [15] Weck A, Wilkinson DS, Maire E (2008) Visualization by x-ray tomography of void growth and coalescence leading to fracture in model materials.Acta Mater 56: 2919.
    [16] Toda H, Yamamoto S, Kobayashi M (2008) Direct measurement procedure for three-dimensional local crack driving force using synchrotron X-ray microtomography.Acta Mater 56: 6027.
    [17] Williams JJ, Flom Z, Amell AA (2010) Damage evolution in SiC particle reinforced Al alloy matrix composites by X-ray synchrotron tomography.Acta Mater 58: 6194.
    [18] Williams JJ, Yazzie KE, Phillips NC (2011) On the correlation between fatigue striation spacing and crack crowth rate: A three-dimensional (3-D) X-ray synchrotron tomography study.Metal Mater Trans 42: 3845.
    [19] Williams JJ, Yazzie KE, Phillips NC, et al.Understanding fatigue crack growth in aluminum alloys by in situ x-ray synchrotron tomography.Int J Fatigue in press .
    [20] Groeber M, Ghosh S, Uchic MD (2008) A framework for automated analysis and simulation of 3D polycrystalline microstructures. I: Statistical characterization.Acta Mater 56: 1257.
    [21] Groeber M, Ghosh S, Uchic MD (2008) A framework for automated analysis and simulation of 3D polycrystalline microstructures. II: Synthetic structure generation.Acta Mater 56: 1274.
    [22] Niezgoda SR, Fullwood DT, Kalidindi SR (2008) Delineation of the space of 2-point correlations in a composite material system.Acta Mater 56: 5286-5292.
    [23] Debye P, Anderson HR, Brumberger H (1957) Scattering by an inhomogeneous solid. II. The correlation function and its applications.J. Appl Phys 28: 679-683.
    [24] Yeong CLY, Torquato S (1998) Reconstructing random media.Phys Rev E 57: 495.
    [25] Yeong CLY, Torquato S (1998) Reconstructing random media: II. Three-dimensional reconstruction from two-dimensional cuts.Phys Rev E 58: 224.
    [26] Jiao Y, Stillinger FH, Torquato S (2007) Modeling heterogeneous materials via two-point correlation functions: Basic principles.Phys Rev E 76: 031110.
    [27] Jiao Y, Stillinger FH, Torquato S (2008) Modeling heterogeneous materials via two-point correlation functions: II. Algorithmic details and applications.Phys Rev E 77: 031135.
    [28] Jiao Y, Stillinger FH, Torquato S (2008) A superior descriptor of random textures and its predictive capacity.Proc Natl Acad Sci USA 106: 17634.
    [29] Roberts AP (1997) Statistical reconstruction of three-dimensional porous media from twodimensional images.Phys Rev E 56: 3203-3212.
    [30] Fullwood DT, Niezgoda SR, Kalidindi SR (2008) Microstructure reconstructions from 2-point statistics using phase-recovery algorithms.Acta Mater 56: 942-948.
    [31] Alireza H, Aliakbar S, Farhad A. (2011) Farhadpour A multiple-point statistics algorithm for 3D pore space reconstruction from 2D images.Advances in Water Resources 34: 1256-1267.
    [32] Tahmasebi P, Sahimi M (2013) Cross-correlation function for accurate reconstruction of heterogeneous media.Phys Rev Lett 110: 078002.
    [33] Torquato S (1986) Interfacial surface statistics arising in diffusion and flow problems in porous media.J Chem Phys 85: 4622.
    [34] Lu B, Torquato S (1992) Lineal path function for random heterogeneous materials.Phys Rev A 45: 922.
    [35] Torquato S, Avellaneda M (1991) Diffusion and reaction in heterogeneous media: Pore-size distribution, relaxation rimes, and mean survival time.J Chem Phys 95: 6477.
    [36] Prager S (1963) Interphase transfer in stationary two-phase media.Chem Eng Sci 18: 228.
    [37] Torquato S, Beasley JD, Chiew YC (1988) Two-point cluster function for continuum percolation.J Chem Phys 88: 6540.
    [38] Singh SS, Williams JJ, Jiao Y, et al. (2012) Modeling anisotropic multiphase heterogeneous materials via directional correlation functions: Simulations and experimental verification.Metall Mater Trans 43 A: 4470-4474.
    [39] Jiao Y, Pallia E, Chawla N (2013) Modeling and predicting microstructure evolution in lead/tin alloy via correlation functions and stochastic material reconstruction.Acta Mater 61: 3370.
    [40] Liu Y, Greene MS, Chen W, et al. (2013) Computational microstructure characterization and reconstruction for stochastic multiscale material design.Computer-Aided Design 45: 65.
    [41] Mikdam A, Belouettar R, Fiorelli D, et al. (2013) A tool for design of heterogeneous materials with desired physical properties using statistical continuum theory.Mater Sci Eng A 564: 493.
    [42] Kirkpatrick S, Gelatt CD, Vecchi MP (1983) Optimization by simulated annealing.Science 220: 671-680.
    [43] Sheehan N, Torquato S (2001) Generating Microstructures with Specified Correlation Functions.J Appl Phys 89: 53.
    [44] Rozman MG, Utz M (2002) Uniqueness of reconstruction of multiphase morphologies from two-point correlation functions.Physl Rev Lett 89: 135501.
    [45] Manko HH (2001) Solders and Soldering: Materials, Design, Production, and Analysis for Reliable Bonding.New York: McGraw-Hill .
    [46] Gommes CJ (2013) Three-dimensional reconstruction of liquid phases in disordered mesopores using in situ small-angle scattering.J Appl Cryst 46: 493-504.
    [47] Zhang H, Srolovitz DJ, Douglas JF (2009) Grain boundaries exhibit the dynamics of glass-forming liquids.Proc Natl Acad Sci USA 106: 7729-7734.
    [48] Chen Z, Chu KT, Srolovitz DJ (2010) Dislocation climb strengthening in systems with immobile obstacles: Three-dimensional level-set simulation study.Phys Rev B 81: 054104.
  • Reader Comments
  • © 2014 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5830) PDF downloads(972) Cited by(2)

Figures and Tables

Figures(5)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog