The new independent solutions of the nonlinear differential equation with time-dependent coefficients (NDE-TC) are discussed, for the first time, by employing experimental device called a drinking bird whose simple ba...The new independent solutions of the nonlinear differential equation with time-dependent coefficients (NDE-TC) are discussed, for the first time, by employing experimental device called a drinking bird whose simple back-and-forth motion develops into water drinking motion. The solution to a drinking bird equation of motion manifests itself the transition from thermodynamic equilibrium to nonequilibrium irreversible states. The independent solution signifying a nonequilibrium thermal state seems to be constructed as if two independent bifurcation solutions are synthesized, and so, the solution is tentatively termed as the bifurcation-integration solution. The bifurcation-integration solution expresses the transition from mechanical and thermodynamic equilibrium to a nonequilibrium irreversible state, which is explicitly shown by the nonlinear differential equation with time-dependent coefficients (NDE-TC). The analysis established a new theoretical approach to nonequilibrium irreversible states, thermomechanical dynamics (TMD). The TMD method enables one to obtain thermodynamically consistent and time-dependent progresses of thermodynamic quantities, by employing the bifurcation-integration solutions of NDE-TC. We hope that the basic properties of bifurcation-integration solutions will be studied and investigated further in mathematics, physics, chemistry and nonlinear sciences in general.展开更多
In this paper, the Auto-B?cklund transformation connected with the homogeneous balance method (HB) and the extended tanh-function method are used to construct new exact solutions for the time-dependent coefficients Ca...In this paper, the Auto-B?cklund transformation connected with the homogeneous balance method (HB) and the extended tanh-function method are used to construct new exact solutions for the time-dependent coefficients Calogero-Degasperis (VCCD) equation. New soliton and periodic solutions of many types are obtained. Furthermore, the soliton propagation is discussed under the effect of the variable coefficients.展开更多
To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main compon...To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main components:a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state,a quantitative relationship between critical irreversible deformation and complex stress state,and evolution characteristics of strength parameters.The proposed model is implemented in a self-developed numerical code,i.e.CASRock.The reliability of the model is validated through experiments.It is indicated that the time-dependent fracturing potential index(xTFPI)at a given time during the attenuation creep stage shows a negative correlation with the extent of excavationinduced damage.The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress,thereby highlighting the 3D stress-dependent characteristic of the model.Finally,the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed.展开更多
4In this paper,we study a nonlinear Petrovsky type equation with nonlinear weak damping,a superlinear source and time-dependent coefficients utt+△^2u-ki(t)|ut|m^2ut=k2(t)|u|^p-2u,x∈Ω,t>0,whereΩis a bounded doma...4In this paper,we study a nonlinear Petrovsky type equation with nonlinear weak damping,a superlinear source and time-dependent coefficients utt+△^2u-ki(t)|ut|m^2ut=k2(t)|u|^p-2u,x∈Ω,t>0,whereΩis a bounded domain in R^n.Under certain conditions on k1(t),k2(t)and the initial-boundary data,the upper bound for blow-up time of the solution with negative initial energy function is given by means of an auxiliary functional and an energy estimate met hod if p>m.Also,a lower bound of blow-up time are obtained by using a Sobolev-type inequality and a first order differential inequality technique for n=2,3 and n>4.展开更多
We investigate the impact of pairwise and group interactions on the spread of epidemics through an activity-driven model based on time-dependent networks.The effects of pairwise/group interaction proportion and pairwi...We investigate the impact of pairwise and group interactions on the spread of epidemics through an activity-driven model based on time-dependent networks.The effects of pairwise/group interaction proportion and pairwise/group interaction intensity are explored by extensive simulation and theoretical analysis.It is demonstrated that altering the group interaction proportion can either hinder or enhance the spread of epidemics,depending on the relative social intensity of group and pairwise interactions.As the group interaction proportion decreases,the impact of reducing group social intensity diminishes.The ratio of group and pairwise social intensity can affect the effect of group interaction proportion on the scale of infection.A weak heterogeneous activity distribution can raise the epidemic threshold,and reduce the scale of infection.These results benefit the design of epidemic control strategy.展开更多
In this paper, we apply Malliavin calculus to discuss when the solutions of stochastic differen-tial equations (SDE's) with time-dependent coefficients have smooth density. Under Hormander'scondition,we conclu...In this paper, we apply Malliavin calculus to discuss when the solutions of stochastic differen-tial equations (SDE's) with time-dependent coefficients have smooth density. Under Hormander'scondition,we conclude that the solutions of the SDE's have smooth density. As a consequence,we get the hypoellipticity for inhomogeneous differential operators.展开更多
Geomorphological features are commonly used to identify potential landslides.Nevertheless,overemphasis on these features could lead to misjudgment.This research proposes a process-oriented approach for potential lands...Geomorphological features are commonly used to identify potential landslides.Nevertheless,overemphasis on these features could lead to misjudgment.This research proposes a process-oriented approach for potential landslide identification that considers time-dependent behaviors.The method integrates comprehensive remote sensing and geological analysis to qualitatively assess slope stability,and employs numerical analysis to quantitatively calculate aging stability.Specifically,a time-dependent stability calculation method for anticlinal slopes is developed and implemented in discrete element software,incorporating time-dependent mechanical and strength reduction calculations.By considering the time-dependent evolution of slopes,this method highlights the importance of both geomorphological features and time-dependent behaviors in landslide identification.This method has been applied to the Jiarishan slope(JRS)on the Qinghai-Tibet Plateau as a case study.The results show that the JRS,despite having landslide geomorphology,is a stable slope,highlighting the risk of misjudgment when relying solely on geomorphological features.This work provides insights into the geomorphological characterization and evolution history of the JRS and offers valuable guidance for studying slopes with similar landslide geomorphology.Furthermore,the process-oriented method incorporating timedependent evolution provides a means to evaluate potential landslides,reducing misjudgment due to excessive reliance on geomorphological features.展开更多
We study equations in divergence form with piecewise Cαcoefficients.The domains contain corners and the discontinuity surfaces are attached to the edges of the corners.We obtain piecewise C^(1,α) estimates across th...We study equations in divergence form with piecewise Cαcoefficients.The domains contain corners and the discontinuity surfaces are attached to the edges of the corners.We obtain piecewise C^(1,α) estimates across the discontinuity surfaces and provide an example to illustrate the issue regarding the regularity at the corners.展开更多
This paper develops a generalized scalar auxiliary variable(SAV)method for the time-dependent Ginzburg-Landau equations.The backward Euler method is used for discretizing the temporal derivative of the time-dependent ...This paper develops a generalized scalar auxiliary variable(SAV)method for the time-dependent Ginzburg-Landau equations.The backward Euler method is used for discretizing the temporal derivative of the time-dependent Ginzburg-Landau equations.In this method,the system is decoupled and linearized to avoid solving the non-linear equation at each step.The theoretical analysis proves that the generalized SAV method can preserve the maximum bound principle and energy stability,and this is confirmed by the numerical result,and also shows that the numerical algorithm is stable.展开更多
Many experiments have demonstrated that resonant magnetic perturbation(RMP) can affect the turbulent transport at the edge of the tokamak. Through the Experimental Advanced Superconducting Tokamak(EAST) density modula...Many experiments have demonstrated that resonant magnetic perturbation(RMP) can affect the turbulent transport at the edge of the tokamak. Through the Experimental Advanced Superconducting Tokamak(EAST) density modulation experiment, the particle transport coefficients were calculated using the experimental data, and the result shows that the particle transport coefficients increase with RMP. In this study, the six-field two-fluid model in BOUT++ is used to simulate the transport before and after density pump-out induced by RMP,respectively referred as the case without RMP and the case with RMP. In the linear simulations,the instabilities generally decreases for cases with RMP. In the nonlinear simulation, ELM only appears in the case without RMP. Additionally, the particle transport coefficient was analyzed,and the result shows that the particle transport coefficient becomes larger for the case with RMP,which is consistent with the experimental conclusion. Moreover, its magnitude is comparable to the results calculated from experimental data.展开更多
This work is a simulation modelling with the LAMMPS calculation code of an electrode based on alkali metals (lithium, sodium and potassium) using the MEAM potential. For different multiplicities, two models were studi...This work is a simulation modelling with the LAMMPS calculation code of an electrode based on alkali metals (lithium, sodium and potassium) using the MEAM potential. For different multiplicities, two models were studied;with and without gap. In this work, we present the structural, physical and chemical properties of the lithium, sodium and potassium electrodes. For the structural properties, the cohesive energy and the mesh parameters were calculated, revealing that, whatever the chemical element selected, the compact hexagonal hcp structure is the most stable, followed by the face-centred cubic CFC structure, and finally the BCC structure. The most stable structure is lithium, with a cohesion energy of -6570 eV, and the lowest bcc-hcp transition energy of -0.553 eV/atom, followed by sodium. For physical properties, kinetic and potential energies were calculated for each of the sectioned chemical elements, with lithium achieving the highest value. Finally, for the chemical properties, we studied the diffusion coefficient and the activation energy. Only potassium followed an opposite order to the other two, with the quantities with lacunae being greater than those without lacunae, whatever the multiplicity. The order of magnitude of the diffusion coefficients is given by the relationship D<sub>Li</sub> > D<sub>Na</sub> > D<sub>k</sub> for the multiplicity 6*6*6, while for the activation energy the order is reversed.展开更多
The rise in breast cancer diagnoses among Chinese women has necessitated the use of X-ray breast screening,which carries a radiation risk.This study aimed to provide a dosimetry protocol for the Chinese female populat...The rise in breast cancer diagnoses among Chinese women has necessitated the use of X-ray breast screening,which carries a radiation risk.This study aimed to provide a dosimetry protocol for the Chinese female population to replace the traditional standard that utilizes simplified breast models,for the accurate estimation of the mean glandular dose of a patient undergoing digital breast tomosynthesis(DBT).The first set of detailed Chinese female breast models and representative breast parameters was constructed.Considering backscatter radiation and computational efficiency,we improved the combination of these models and the Chinese reference adult female whole-body voxel phantom.Image acquisition for four commercial DBT systems that are widely employed in China was simulated using the Monte Carlo method to obtain the normalized glandular dose coefficients of DBT(D_(gN)^(DBT))and the glandular depth dose(D_(g)^(dep)(z))for different breast characteristics and X-ray spectra.We calculated a series of D_(gN)^(DBT) values for breasts with different percentage mass glandularities(5%,25%,50%,75%,and 100%)and compressed breast thicknesses(2,3,4,5,6,and 7 cm)at various tube potentials(25,28,30,32,35,and 49 kV)and target/filter combinations(W/Rh,W/Al,Mo/Mo,Rh/Rh,and Rh/Ag).The parameter dependence of the breast characteristics and beam conditions on D_(gN)^(DBT) in detailed breast models was investigated.The D_(gN)^(DBT) results were 14.6-51.0%lower than those of the traditional dosimetry standard in China.The difference in D_(gN)^(DBT) was mainly due to a decrease in the depth of the main energy deposition area caused by the glandular distribution along the depth direction.The results obtained in this study may be used to improve breast dosimetry in China and provide more detailed information on risk assessment during DBT.展开更多
Adaptive fractional polynomial modeling of general correlated outcomes is formulated to address nonlinearity in means, variances/dispersions, and correlations. Means and variances/dispersions are modeled using general...Adaptive fractional polynomial modeling of general correlated outcomes is formulated to address nonlinearity in means, variances/dispersions, and correlations. Means and variances/dispersions are modeled using generalized linear models in fixed effects/coefficients. Correlations are modeled using random effects/coefficients. Nonlinearity is addressed using power transforms of primary (untransformed) predictors. Parameter estimation is based on extended linear mixed modeling generalizing both generalized estimating equations and linear mixed modeling. Models are evaluated using likelihood cross-validation (LCV) scores and are generated adaptively using a heuristic search controlled by LCV scores. Cases covered include linear, Poisson, logistic, exponential, and discrete regression of correlated continuous, count/rate, dichotomous, positive continuous, and discrete numeric outcomes treated as normally, Poisson, Bernoulli, exponentially, and discrete numerically distributed, respectively. Example analyses are also generated for these five cases to compare adaptive random effects/coefficients modeling of correlated outcomes to previously developed adaptive modeling based on directly specified covariance structures. Adaptive random effects/coefficients modeling substantially outperforms direct covariance modeling in the linear, exponential, and discrete regression example analyses. It generates equivalent results in the logistic regression example analyses and it is substantially outperformed in the Poisson regression case. Random effects/coefficients modeling of correlated outcomes can provide substantial improvements in model selection compared to directly specified covariance modeling. However, directly specified covariance modeling can generate competitive or substantially better results in some cases while usually requiring less computation time.展开更多
[Introduction] Accurate calculation of the hydrodynamic coefficients for floating structures and the investigation of the flow field distribution around floating bodies on the marine free surface are essential for imp...[Introduction] Accurate calculation of the hydrodynamic coefficients for floating structures and the investigation of the flow field distribution around floating bodies on the marine free surface are essential for improving the engineering design and application of marine structures.[Method] This study utilized the computational fluid dynamics(CFD) approach and the Reynolds Averaged NavierStokes(RANS) method and considered the effects of viscosity and free surface interactions on the hydrodynamic behavior of floating structures.By employing the dynamic mesh technique,this study simulated the periodic movements of simplified three-dimensional(3D)shapes:spheres,cylinders,and cubes,which were representative of complex marine structures.The volume of fluid(VOF) method was leveraged to accurately track the nonlinear behavior of the free surface.In this analysis,the added mass and damping coefficients for the fundamental modes of motion(surge,heave,and roll) were calculated across a spectrum of frequencies,facilitating the fast determination of hydrodynamic forces and moments exerted on floating structures.[Result] The results of this study are not only consistent with the results of the 3D potential flow theory but also further reflect the role of viscosity.This method can be used for precise calculation of the hydrodynamic coefficients of floating structures and for describing the flow field of such structures in motion on a free surface.[Conclusion] The methodology presented goes beyond the traditional potential flow approach.展开更多
This work is a simulation model with the LAMMPS calculation code of an electrode based on alkali metal oxides (lithium, sodium and potassium) using the Lennard Jones potential. For a multiplicity of 8*8*8, we studied ...This work is a simulation model with the LAMMPS calculation code of an electrode based on alkali metal oxides (lithium, sodium and potassium) using the Lennard Jones potential. For a multiplicity of 8*8*8, we studied a gap-free model using molecular dynamics. Physical quantities such as volume and pressure of the Na-O and Li-O systems exhibit similar behaviors around the thermodynamic ensembles NPT and NVE. However, for the Na2O system, at a minimum temperature value, we observe a range of total energy values;in contrast, for the Li2O system, a minimum energy corresponds to a range of temperatures. Finally, for physicochemical properties, we studied the diffusion coefficient and activation energy of lithium and potassium oxides around their melting temperatures. The order of magnitude of the diffusion coefficients is given by the relation Dli-O >DNa-O for the multiplicity 8*8*8, while for the activation energy, the order is well reversed EaNa-O > EaLi-O.展开更多
The populations of urban centers in Congo-Brazzaville have decided to develop various methods of water storage (concrete or masonry underground tanks) for domestic use, due to shortages in the distribution of water th...The populations of urban centers in Congo-Brazzaville have decided to develop various methods of water storage (concrete or masonry underground tanks) for domestic use, due to shortages in the distribution of water through the public network. However, questions remain as to the physico-chemical quality of the water stored in these tanks, when these structures are built in wet and relatively polluted areas. This paper presents a model of pollutant diffusion through the cementitious matrix (concrete) of tank walls simulated at a buried reservoir. The results of the experimental and numerical simulations show that certain concrete parameters, such as porosity, permeability and diffusivity, have a significant influence on the transfer of pollutants through the concrete walls, thus altering the physico-chemical quality of the stored water. The numerical models (1D) used to predict pollutant transfer and the quality of the stored water are consistent with those of the optimal control for identifying the diffusion coefficient. Major ion concentrations appear to be correlated with system porosity and diffusion coefficient. Nevertheless, the identification of the diffusion coefficient from the optimal control method, based on an explicit numerical resolution of a finite volume PDE for the approximation of the experiment, is not consistent with that of the optimal control method.展开更多
In this paper, some exact solutions of the (3 + 1)-dimensional variable-coefficient Yu-Toda-Sasa-Fukuyama equation are investigated. By using Hirota’s direct method and symbolic computation, we obtained N-soliton sol...In this paper, some exact solutions of the (3 + 1)-dimensional variable-coefficient Yu-Toda-Sasa-Fukuyama equation are investigated. By using Hirota’s direct method and symbolic computation, we obtained N-soliton solution. By using the long wave limit method, the N-order rational solution can be obtained from N-order soliton solution. Then, through the paired complexification of parameters, the lump solution is obtained from N-order rational solution. Meanwhile, we obtained a hybrid solution between 1-lump solution and N-soliton (N=1,2) by using the long wave limit method and parameter complex. Furthermore, four different sets of three-dimensional graphs of solitons, lump solutions and hybrid solutions are drawn by selecting four different sets of coefficient functions which include one set of constant coefficient function and three sets of variable coefficient functions.展开更多
Given a sample of regression data from (Y, Z), a new diagnostic plotting method is proposed for checking the hypothesis H0: the data are from a given Cox model with the time-dependent covariates Z. It compares two est...Given a sample of regression data from (Y, Z), a new diagnostic plotting method is proposed for checking the hypothesis H0: the data are from a given Cox model with the time-dependent covariates Z. It compares two estimates of the marginal distribution FY of Y. One is an estimate of the modified expression of FY under H0, based on a consistent estimate of the parameter under H0, and based on the baseline distribution of the data. The other is the Kaplan-Meier-estimator of FY, together with its confidence band. The new plot, called the marginal distribution plot, can be viewed as a test for testing H0. The main advantage of the test over the existing residual tests is in the case that the data do not satisfy any Cox model or the Cox model is mis-specified. Then the new test is still valid, but not the residual tests and the residual tests often make type II error with a very large probability.展开更多
When multiphysics coupling calculations contain time-dependent Monte Carlo particle transport simulations, these simulations often account for the largest part of the calculation time, which is insufferable in certain...When multiphysics coupling calculations contain time-dependent Monte Carlo particle transport simulations, these simulations often account for the largest part of the calculation time, which is insufferable in certain important cases. This study proposes an adaptive strategy for automatically adjusting the sample size to fulfil more reasonable simulations. This is realized based on an extension of the Shannon entropy concept and is essentially different from the popular methods in timeindependent Monte Carlo particle transport simulations, such as controlling the sample size according to the relative error of a target tally or by experience. The results of the two models show that this strategy can yield almost similar results while significantly reducing the calculation time. Considering the efficiency, the sample size should not be increased blindly if the efficiency cannot be enhanced further. The strategy proposed herein satisfies this requirement.展开更多
Time-dependent diffusion coefficient and conventional diffusion constant are calculated and analyzed to study diffusion of nanoparticles in polymer melts. A generalized Langevin equa- tion is adopted to describe the d...Time-dependent diffusion coefficient and conventional diffusion constant are calculated and analyzed to study diffusion of nanoparticles in polymer melts. A generalized Langevin equa- tion is adopted to describe the diffusion dynamics. Mode-coupling theory is employed to calculate the memory kernel of friction. For simplicity, only microscopic terms arising from binary collision and coupling to the solvent density fluctuation are included in the formalism. The equilibrium structural information functions of the polymer nanocomposites required by mode-coupling theory are calculated on the basis of polymer reference interaction site model with Percus-Yevick closure. The effect of nanoparticle size and that of the polymer size are clarified explicitly. The structural functions, the friction kernel, as well as the diffusion coefficient show a rich variety with varying nanoparticle radius and polymer chain length. We find that for small nanoparticles or short chain polymers, the characteristic short time non-Markov diffusion dynamics becomes more prominent, and the diffusion coefficient takes longer time to approach asymptotically the conventional diffusion constant. This constant due to the microscopic contributions will decrease with the increase of nanoparticle size, while increase with polymer size. Furthermore, our result of diffusion constant from mode- coupling theory is compared with the value predicted from the Stokes-Einstein relation. It shows that the microscopic contributions to the diffusion constant are dominant for small nanoparticles or long chain polymers. Inversely, when nanonparticle is big, or polymer chain is short, the hydrodynamic contribution might play a significant role.展开更多
文摘The new independent solutions of the nonlinear differential equation with time-dependent coefficients (NDE-TC) are discussed, for the first time, by employing experimental device called a drinking bird whose simple back-and-forth motion develops into water drinking motion. The solution to a drinking bird equation of motion manifests itself the transition from thermodynamic equilibrium to nonequilibrium irreversible states. The independent solution signifying a nonequilibrium thermal state seems to be constructed as if two independent bifurcation solutions are synthesized, and so, the solution is tentatively termed as the bifurcation-integration solution. The bifurcation-integration solution expresses the transition from mechanical and thermodynamic equilibrium to a nonequilibrium irreversible state, which is explicitly shown by the nonlinear differential equation with time-dependent coefficients (NDE-TC). The analysis established a new theoretical approach to nonequilibrium irreversible states, thermomechanical dynamics (TMD). The TMD method enables one to obtain thermodynamically consistent and time-dependent progresses of thermodynamic quantities, by employing the bifurcation-integration solutions of NDE-TC. We hope that the basic properties of bifurcation-integration solutions will be studied and investigated further in mathematics, physics, chemistry and nonlinear sciences in general.
文摘In this paper, the Auto-B?cklund transformation connected with the homogeneous balance method (HB) and the extended tanh-function method are used to construct new exact solutions for the time-dependent coefficients Calogero-Degasperis (VCCD) equation. New soliton and periodic solutions of many types are obtained. Furthermore, the soliton propagation is discussed under the effect of the variable coefficients.
基金supported by the National Natural Science Foundation of China(Grant No.52125903)the China Postdoctoral Science Foundation(Grant No.2023M730367)the Fundamental Research Funds for Central Public Welfare Research Institutes of China(Grant No.CKSF2023323/YT).
文摘To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main components:a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state,a quantitative relationship between critical irreversible deformation and complex stress state,and evolution characteristics of strength parameters.The proposed model is implemented in a self-developed numerical code,i.e.CASRock.The reliability of the model is validated through experiments.It is indicated that the time-dependent fracturing potential index(xTFPI)at a given time during the attenuation creep stage shows a negative correlation with the extent of excavationinduced damage.The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress,thereby highlighting the 3D stress-dependent characteristic of the model.Finally,the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed.
基金This paper is supported by the Natural Science Foundation of Shandong Province(No.ZR2018BA016).
文摘4In this paper,we study a nonlinear Petrovsky type equation with nonlinear weak damping,a superlinear source and time-dependent coefficients utt+△^2u-ki(t)|ut|m^2ut=k2(t)|u|^p-2u,x∈Ω,t>0,whereΩis a bounded domain in R^n.Under certain conditions on k1(t),k2(t)and the initial-boundary data,the upper bound for blow-up time of the solution with negative initial energy function is given by means of an auxiliary functional and an energy estimate met hod if p>m.Also,a lower bound of blow-up time are obtained by using a Sobolev-type inequality and a first order differential inequality technique for n=2,3 and n>4.
基金This work was supported by the National Natural Science Foundation of China(Grant No.12072340)the China Postdoctoral Science Foundation(Grant No.2022M720727)the Jiangsu Funding Program for Excellent Postdoctoral Talent(Grant No.2022ZB130).
文摘We investigate the impact of pairwise and group interactions on the spread of epidemics through an activity-driven model based on time-dependent networks.The effects of pairwise/group interaction proportion and pairwise/group interaction intensity are explored by extensive simulation and theoretical analysis.It is demonstrated that altering the group interaction proportion can either hinder or enhance the spread of epidemics,depending on the relative social intensity of group and pairwise interactions.As the group interaction proportion decreases,the impact of reducing group social intensity diminishes.The ratio of group and pairwise social intensity can affect the effect of group interaction proportion on the scale of infection.A weak heterogeneous activity distribution can raise the epidemic threshold,and reduce the scale of infection.These results benefit the design of epidemic control strategy.
基金The project supported by National Natural Science Foundation of China Crant 18971061
文摘In this paper, we apply Malliavin calculus to discuss when the solutions of stochastic differen-tial equations (SDE's) with time-dependent coefficients have smooth density. Under Hormander'scondition,we conclude that the solutions of the SDE's have smooth density. As a consequence,we get the hypoellipticity for inhomogeneous differential operators.
基金This research was supported by the National Natural Science Foundation of China(Grant Nos.41972284 and 42090054)This work was also supported by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project(Grant No.SKLGP2020Z005).
文摘Geomorphological features are commonly used to identify potential landslides.Nevertheless,overemphasis on these features could lead to misjudgment.This research proposes a process-oriented approach for potential landslide identification that considers time-dependent behaviors.The method integrates comprehensive remote sensing and geological analysis to qualitatively assess slope stability,and employs numerical analysis to quantitatively calculate aging stability.Specifically,a time-dependent stability calculation method for anticlinal slopes is developed and implemented in discrete element software,incorporating time-dependent mechanical and strength reduction calculations.By considering the time-dependent evolution of slopes,this method highlights the importance of both geomorphological features and time-dependent behaviors in landslide identification.This method has been applied to the Jiarishan slope(JRS)on the Qinghai-Tibet Plateau as a case study.The results show that the JRS,despite having landslide geomorphology,is a stable slope,highlighting the risk of misjudgment when relying solely on geomorphological features.This work provides insights into the geomorphological characterization and evolution history of the JRS and offers valuable guidance for studying slopes with similar landslide geomorphology.Furthermore,the process-oriented method incorporating timedependent evolution provides a means to evaluate potential landslides,reducing misjudgment due to excessive reliance on geomorphological features.
基金supported by National Natural Science Foundation of China(12061080,12161087 and 12261093)the Science and Technology Project of the Education Department of Jiangxi Province(GJJ211601)supported by National Natural Science Foundation of China(11871305).
文摘We study equations in divergence form with piecewise Cαcoefficients.The domains contain corners and the discontinuity surfaces are attached to the edges of the corners.We obtain piecewise C^(1,α) estimates across the discontinuity surfaces and provide an example to illustrate the issue regarding the regularity at the corners.
基金supported by the National Natural Science Foundation of China(12126318,12126302).
文摘This paper develops a generalized scalar auxiliary variable(SAV)method for the time-dependent Ginzburg-Landau equations.The backward Euler method is used for discretizing the temporal derivative of the time-dependent Ginzburg-Landau equations.In this method,the system is decoupled and linearized to avoid solving the non-linear equation at each step.The theoretical analysis proves that the generalized SAV method can preserve the maximum bound principle and energy stability,and this is confirmed by the numerical result,and also shows that the numerical algorithm is stable.
基金supported by the National Magnetic Confinement Fusion Program of China(No.2019YFE03090200)by National Natural Science Foundation of China(Nos.11975231,12175277 and 12305249).
文摘Many experiments have demonstrated that resonant magnetic perturbation(RMP) can affect the turbulent transport at the edge of the tokamak. Through the Experimental Advanced Superconducting Tokamak(EAST) density modulation experiment, the particle transport coefficients were calculated using the experimental data, and the result shows that the particle transport coefficients increase with RMP. In this study, the six-field two-fluid model in BOUT++ is used to simulate the transport before and after density pump-out induced by RMP,respectively referred as the case without RMP and the case with RMP. In the linear simulations,the instabilities generally decreases for cases with RMP. In the nonlinear simulation, ELM only appears in the case without RMP. Additionally, the particle transport coefficient was analyzed,and the result shows that the particle transport coefficient becomes larger for the case with RMP,which is consistent with the experimental conclusion. Moreover, its magnitude is comparable to the results calculated from experimental data.
文摘This work is a simulation modelling with the LAMMPS calculation code of an electrode based on alkali metals (lithium, sodium and potassium) using the MEAM potential. For different multiplicities, two models were studied;with and without gap. In this work, we present the structural, physical and chemical properties of the lithium, sodium and potassium electrodes. For the structural properties, the cohesive energy and the mesh parameters were calculated, revealing that, whatever the chemical element selected, the compact hexagonal hcp structure is the most stable, followed by the face-centred cubic CFC structure, and finally the BCC structure. The most stable structure is lithium, with a cohesion energy of -6570 eV, and the lowest bcc-hcp transition energy of -0.553 eV/atom, followed by sodium. For physical properties, kinetic and potential energies were calculated for each of the sectioned chemical elements, with lithium achieving the highest value. Finally, for the chemical properties, we studied the diffusion coefficient and the activation energy. Only potassium followed an opposite order to the other two, with the quantities with lacunae being greater than those without lacunae, whatever the multiplicity. The order of magnitude of the diffusion coefficients is given by the relationship D<sub>Li</sub> > D<sub>Na</sub> > D<sub>k</sub> for the multiplicity 6*6*6, while for the activation energy the order is reversed.
基金supported by the National Natural Science Foundation of China(Nos.U2167209 and 12175114)the National Key R&D Program of China(No.2021YFF0603600).
文摘The rise in breast cancer diagnoses among Chinese women has necessitated the use of X-ray breast screening,which carries a radiation risk.This study aimed to provide a dosimetry protocol for the Chinese female population to replace the traditional standard that utilizes simplified breast models,for the accurate estimation of the mean glandular dose of a patient undergoing digital breast tomosynthesis(DBT).The first set of detailed Chinese female breast models and representative breast parameters was constructed.Considering backscatter radiation and computational efficiency,we improved the combination of these models and the Chinese reference adult female whole-body voxel phantom.Image acquisition for four commercial DBT systems that are widely employed in China was simulated using the Monte Carlo method to obtain the normalized glandular dose coefficients of DBT(D_(gN)^(DBT))and the glandular depth dose(D_(g)^(dep)(z))for different breast characteristics and X-ray spectra.We calculated a series of D_(gN)^(DBT) values for breasts with different percentage mass glandularities(5%,25%,50%,75%,and 100%)and compressed breast thicknesses(2,3,4,5,6,and 7 cm)at various tube potentials(25,28,30,32,35,and 49 kV)and target/filter combinations(W/Rh,W/Al,Mo/Mo,Rh/Rh,and Rh/Ag).The parameter dependence of the breast characteristics and beam conditions on D_(gN)^(DBT) in detailed breast models was investigated.The D_(gN)^(DBT) results were 14.6-51.0%lower than those of the traditional dosimetry standard in China.The difference in D_(gN)^(DBT) was mainly due to a decrease in the depth of the main energy deposition area caused by the glandular distribution along the depth direction.The results obtained in this study may be used to improve breast dosimetry in China and provide more detailed information on risk assessment during DBT.
文摘Adaptive fractional polynomial modeling of general correlated outcomes is formulated to address nonlinearity in means, variances/dispersions, and correlations. Means and variances/dispersions are modeled using generalized linear models in fixed effects/coefficients. Correlations are modeled using random effects/coefficients. Nonlinearity is addressed using power transforms of primary (untransformed) predictors. Parameter estimation is based on extended linear mixed modeling generalizing both generalized estimating equations and linear mixed modeling. Models are evaluated using likelihood cross-validation (LCV) scores and are generated adaptively using a heuristic search controlled by LCV scores. Cases covered include linear, Poisson, logistic, exponential, and discrete regression of correlated continuous, count/rate, dichotomous, positive continuous, and discrete numeric outcomes treated as normally, Poisson, Bernoulli, exponentially, and discrete numerically distributed, respectively. Example analyses are also generated for these five cases to compare adaptive random effects/coefficients modeling of correlated outcomes to previously developed adaptive modeling based on directly specified covariance structures. Adaptive random effects/coefficients modeling substantially outperforms direct covariance modeling in the linear, exponential, and discrete regression example analyses. It generates equivalent results in the logistic regression example analyses and it is substantially outperformed in the Poisson regression case. Random effects/coefficients modeling of correlated outcomes can provide substantial improvements in model selection compared to directly specified covariance modeling. However, directly specified covariance modeling can generate competitive or substantially better results in some cases while usually requiring less computation time.
文摘[Introduction] Accurate calculation of the hydrodynamic coefficients for floating structures and the investigation of the flow field distribution around floating bodies on the marine free surface are essential for improving the engineering design and application of marine structures.[Method] This study utilized the computational fluid dynamics(CFD) approach and the Reynolds Averaged NavierStokes(RANS) method and considered the effects of viscosity and free surface interactions on the hydrodynamic behavior of floating structures.By employing the dynamic mesh technique,this study simulated the periodic movements of simplified three-dimensional(3D)shapes:spheres,cylinders,and cubes,which were representative of complex marine structures.The volume of fluid(VOF) method was leveraged to accurately track the nonlinear behavior of the free surface.In this analysis,the added mass and damping coefficients for the fundamental modes of motion(surge,heave,and roll) were calculated across a spectrum of frequencies,facilitating the fast determination of hydrodynamic forces and moments exerted on floating structures.[Result] The results of this study are not only consistent with the results of the 3D potential flow theory but also further reflect the role of viscosity.This method can be used for precise calculation of the hydrodynamic coefficients of floating structures and for describing the flow field of such structures in motion on a free surface.[Conclusion] The methodology presented goes beyond the traditional potential flow approach.
文摘This work is a simulation model with the LAMMPS calculation code of an electrode based on alkali metal oxides (lithium, sodium and potassium) using the Lennard Jones potential. For a multiplicity of 8*8*8, we studied a gap-free model using molecular dynamics. Physical quantities such as volume and pressure of the Na-O and Li-O systems exhibit similar behaviors around the thermodynamic ensembles NPT and NVE. However, for the Na2O system, at a minimum temperature value, we observe a range of total energy values;in contrast, for the Li2O system, a minimum energy corresponds to a range of temperatures. Finally, for physicochemical properties, we studied the diffusion coefficient and activation energy of lithium and potassium oxides around their melting temperatures. The order of magnitude of the diffusion coefficients is given by the relation Dli-O >DNa-O for the multiplicity 8*8*8, while for the activation energy, the order is well reversed EaNa-O > EaLi-O.
文摘The populations of urban centers in Congo-Brazzaville have decided to develop various methods of water storage (concrete or masonry underground tanks) for domestic use, due to shortages in the distribution of water through the public network. However, questions remain as to the physico-chemical quality of the water stored in these tanks, when these structures are built in wet and relatively polluted areas. This paper presents a model of pollutant diffusion through the cementitious matrix (concrete) of tank walls simulated at a buried reservoir. The results of the experimental and numerical simulations show that certain concrete parameters, such as porosity, permeability and diffusivity, have a significant influence on the transfer of pollutants through the concrete walls, thus altering the physico-chemical quality of the stored water. The numerical models (1D) used to predict pollutant transfer and the quality of the stored water are consistent with those of the optimal control for identifying the diffusion coefficient. Major ion concentrations appear to be correlated with system porosity and diffusion coefficient. Nevertheless, the identification of the diffusion coefficient from the optimal control method, based on an explicit numerical resolution of a finite volume PDE for the approximation of the experiment, is not consistent with that of the optimal control method.
文摘In this paper, some exact solutions of the (3 + 1)-dimensional variable-coefficient Yu-Toda-Sasa-Fukuyama equation are investigated. By using Hirota’s direct method and symbolic computation, we obtained N-soliton solution. By using the long wave limit method, the N-order rational solution can be obtained from N-order soliton solution. Then, through the paired complexification of parameters, the lump solution is obtained from N-order rational solution. Meanwhile, we obtained a hybrid solution between 1-lump solution and N-soliton (N=1,2) by using the long wave limit method and parameter complex. Furthermore, four different sets of three-dimensional graphs of solitons, lump solutions and hybrid solutions are drawn by selecting four different sets of coefficient functions which include one set of constant coefficient function and three sets of variable coefficient functions.
文摘Given a sample of regression data from (Y, Z), a new diagnostic plotting method is proposed for checking the hypothesis H0: the data are from a given Cox model with the time-dependent covariates Z. It compares two estimates of the marginal distribution FY of Y. One is an estimate of the modified expression of FY under H0, based on a consistent estimate of the parameter under H0, and based on the baseline distribution of the data. The other is the Kaplan-Meier-estimator of FY, together with its confidence band. The new plot, called the marginal distribution plot, can be viewed as a test for testing H0. The main advantage of the test over the existing residual tests is in the case that the data do not satisfy any Cox model or the Cox model is mis-specified. Then the new test is still valid, but not the residual tests and the residual tests often make type II error with a very large probability.
基金supported by the CAEP Found (No.CX20200028)Youth Program of National Natural Science Foundation of China (No.11705011).
文摘When multiphysics coupling calculations contain time-dependent Monte Carlo particle transport simulations, these simulations often account for the largest part of the calculation time, which is insufferable in certain important cases. This study proposes an adaptive strategy for automatically adjusting the sample size to fulfil more reasonable simulations. This is realized based on an extension of the Shannon entropy concept and is essentially different from the popular methods in timeindependent Monte Carlo particle transport simulations, such as controlling the sample size according to the relative error of a target tally or by experience. The results of the two models show that this strategy can yield almost similar results while significantly reducing the calculation time. Considering the efficiency, the sample size should not be increased blindly if the efficiency cannot be enhanced further. The strategy proposed herein satisfies this requirement.
基金This work was supported by the National Natural Science Foundation of China (No.21173152), the Ministry of Education of China (No.NCET-11-0359 and No.2011SCU04B31), and the Science and Technology Department of Sichuan Province (No.2011HH0005).
文摘Time-dependent diffusion coefficient and conventional diffusion constant are calculated and analyzed to study diffusion of nanoparticles in polymer melts. A generalized Langevin equa- tion is adopted to describe the diffusion dynamics. Mode-coupling theory is employed to calculate the memory kernel of friction. For simplicity, only microscopic terms arising from binary collision and coupling to the solvent density fluctuation are included in the formalism. The equilibrium structural information functions of the polymer nanocomposites required by mode-coupling theory are calculated on the basis of polymer reference interaction site model with Percus-Yevick closure. The effect of nanoparticle size and that of the polymer size are clarified explicitly. The structural functions, the friction kernel, as well as the diffusion coefficient show a rich variety with varying nanoparticle radius and polymer chain length. We find that for small nanoparticles or short chain polymers, the characteristic short time non-Markov diffusion dynamics becomes more prominent, and the diffusion coefficient takes longer time to approach asymptotically the conventional diffusion constant. This constant due to the microscopic contributions will decrease with the increase of nanoparticle size, while increase with polymer size. Furthermore, our result of diffusion constant from mode- coupling theory is compared with the value predicted from the Stokes-Einstein relation. It shows that the microscopic contributions to the diffusion constant are dominant for small nanoparticles or long chain polymers. Inversely, when nanonparticle is big, or polymer chain is short, the hydrodynamic contribution might play a significant role.