A numerical study is presented for the fully developed two-dimensional laminar flow of viscous incompressible fluid through a curved square duct for the constant curvature δ = 0.1. In this paper, a spectral-based com...A numerical study is presented for the fully developed two-dimensional laminar flow of viscous incompressible fluid through a curved square duct for the constant curvature δ = 0.1. In this paper, a spectral-based computational algorithm is employed as the principal tool for the simulations, while a Chebyshev polynomial and collocation method as secondary tools. Numerical calculations are carried out over a wide range of the pressure gradient parameter, the Dean number, 100 ≤ Dn ≤ 3000 for the Grashof number, Gr, ranging from 100 to 2000. The outer wall of the duct is treated heated while the inner wall cooled, the top and bottom walls being adiabatic. The main concern of the present study is to find out the unsteady flow behavior i.e. whether the unsteady flow is steady-state, periodic, multi-periodic or chaotic, if Dn or Gr is increased. It is found that the unsteady flow is periodic for Dn = 1000 at Gr = 100 and 500 and at Dn = 2000, Gr = 2000 but steady-state otherwise. It is also found that for large values of Dn, for example Dn = 3000, the unsteady flow undergoes in the scenario “periodic→chaotic→periodic”, if Gr is increased. Typical contours of secondary flow patterns and temperature profiles are also obtained, and it is found that the unsteady flow consists of single-, two-, three- and four-vortex solutions. The present study also shows that there is a strong interaction between the heating-induced buoyancy force and the centrifugal force in a curved square passage that stimulates fluid mixing and consequently enhance heat transfer in the fluid.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Diffusions of multiple components have numerous applications such as underground water flow, pollutant movement, stratospheric warming, and food processing. Particularly, liquid hydrogen is used in the cooling process...Diffusions of multiple components have numerous applications such as underground water flow, pollutant movement, stratospheric warming, and food processing. Particularly, liquid hydrogen is used in the cooling process of the aeroplane. Further, liquid nitrogen can find applications in cooling equipment or electronic devices, i.e., high temperature superconducting(HTS) cables. So, herein, we have analysed the entropy generation(EG), nonlinear thermal radiation and unsteady(time-dependent) nature of the flow on quadratic combined convective flow over a permeable slender cylinder with diffusions of liquid hydrogen and nitrogen. The governing equations for flow and heat transfer characteristics are expressed in terms of nonlinear coupled partial differential equations. The solutions of these equations are attempted numerically by employing the quasilinearization technique with the implicit finite difference approximation. It is found that EG is minimum for double diffusion(liquid hydrogen and heat diffusion)than triple diffusion(diffusion of liquid hydrogen, nitrogen and heat). The enhancing values of the radiation parameter R_(d) and temperature ratio θ_(w) augment the fluid temperature for steady and unsteady cases as well as the local Nusselt number. Because, the fluid absorbs the heat energy released due to radiation, and in turn releases the heat energy from the cylinder to the surrounding surface.展开更多
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.展开更多
The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resour...The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation(DES) method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.展开更多
We present an efficient approach to solve multi-dimensional time-dependent Schr?dinger equation(TDSE)in an intense laser field.In this approach,each spatial degree of freedom is treated as a distinguishable quasi-part...We present an efficient approach to solve multi-dimensional time-dependent Schr?dinger equation(TDSE)in an intense laser field.In this approach,each spatial degree of freedom is treated as a distinguishable quasi-particle.The non-separable Coulomb potential is regarded as a two-body operator between different quasi-particles.The time-dependent variational principle is used to derive the equations of motion.Then the high-order multi-dimensional problem is broken down into several lower-order coupled equations,which can be efficiently solved.As a demonstration,we apply this method to solve the two-dimensional TDSE.The accuracy is tested by comparing the direct solutions of TDSE using several examples such as the strong-field ionization and the high harmonic generation.The results show that the present method is much more computationally efficient than the conventional one without sacrificing accuracy.The present method can be straightforwardly extended to three-dimensional problems.Our study provides a flexible method to investigate the laser-atom interaction in the nonperturbative regime.展开更多
In this paper,the non-harmonic resonance of Bernoulli viscoelastic beams,Kirchhoff viscoelastic plates,Timoshenko viscoelastic beams,and Mindlin viscoelastic plates subjected to time-dependent exponentially decreasing...In this paper,the non-harmonic resonance of Bernoulli viscoelastic beams,Kirchhoff viscoelastic plates,Timoshenko viscoelastic beams,and Mindlin viscoelastic plates subjected to time-dependent exponentially decreasing transverse distributed load is investigated for the first time.The constitutive equations are expressed utilizing Boltzmann integral law with a constant bulk modulus.The displacement vector is approximated by employing the separation of variables method.The Laplace transformation is used to transfer equations from the time domain to the Laplace domain and vice versa.The novel point of the proposed method is to express,prove and calculate the critical time in which the displacement will be several times the displacement at time zero.In addition,this new method calculates the maximum deflection at the critical time,explicitly and exactly,without any need to follow the time-displacement curve with a low computational cost.Additionally,the proposed method introduces the critical range of time so that the responses are greater than the responses at time zero.展开更多
This paper mainly studies the blowup phenomenon of solutions to the compressible Euler equations with general time-dependent damping for non-isentropic fluids in two and three space dimensions. When the initial data i...This paper mainly studies the blowup phenomenon of solutions to the compressible Euler equations with general time-dependent damping for non-isentropic fluids in two and three space dimensions. When the initial data is assumed to be radially symmetric and the initial density contains vacuum, we obtain that classical solution, especially the density, will blow up on finite time. The results also reveal that damping can really delay the singularity formation.展开更多
On the basis of the quasi-geostrophic vorticity equation,theoretical research has been down upon the evolution of the amplitude of solitary Rossby waves employing the perturbation method,and come to the conclusion tha...On the basis of the quasi-geostrophic vorticity equation,theoretical research has been down upon the evolution of the amplitude of solitary Rossby waves employing the perturbation method,and come to the conclusion that the evolution of the amplitude satisfies the variable coefficient Korteweg-de Vries(KdV) equation.展开更多
文摘A numerical study is presented for the fully developed two-dimensional laminar flow of viscous incompressible fluid through a curved square duct for the constant curvature δ = 0.1. In this paper, a spectral-based computational algorithm is employed as the principal tool for the simulations, while a Chebyshev polynomial and collocation method as secondary tools. Numerical calculations are carried out over a wide range of the pressure gradient parameter, the Dean number, 100 ≤ Dn ≤ 3000 for the Grashof number, Gr, ranging from 100 to 2000. The outer wall of the duct is treated heated while the inner wall cooled, the top and bottom walls being adiabatic. The main concern of the present study is to find out the unsteady flow behavior i.e. whether the unsteady flow is steady-state, periodic, multi-periodic or chaotic, if Dn or Gr is increased. It is found that the unsteady flow is periodic for Dn = 1000 at Gr = 100 and 500 and at Dn = 2000, Gr = 2000 but steady-state otherwise. It is also found that for large values of Dn, for example Dn = 3000, the unsteady flow undergoes in the scenario “periodic→chaotic→periodic”, if Gr is increased. Typical contours of secondary flow patterns and temperature profiles are also obtained, and it is found that the unsteady flow consists of single-, two-, three- and four-vortex solutions. The present study also shows that there is a strong interaction between the heating-induced buoyancy force and the centrifugal force in a curved square passage that stimulates fluid mixing and consequently enhance heat transfer in the fluid.
基金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.
基金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 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.
文摘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.
文摘Diffusions of multiple components have numerous applications such as underground water flow, pollutant movement, stratospheric warming, and food processing. Particularly, liquid hydrogen is used in the cooling process of the aeroplane. Further, liquid nitrogen can find applications in cooling equipment or electronic devices, i.e., high temperature superconducting(HTS) cables. So, herein, we have analysed the entropy generation(EG), nonlinear thermal radiation and unsteady(time-dependent) nature of the flow on quadratic combined convective flow over a permeable slender cylinder with diffusions of liquid hydrogen and nitrogen. The governing equations for flow and heat transfer characteristics are expressed in terms of nonlinear coupled partial differential equations. The solutions of these equations are attempted numerically by employing the quasilinearization technique with the implicit finite difference approximation. It is found that EG is minimum for double diffusion(liquid hydrogen and heat diffusion)than triple diffusion(diffusion of liquid hydrogen, nitrogen and heat). The enhancing values of the radiation parameter R_(d) and temperature ratio θ_(w) augment the fluid temperature for steady and unsteady cases as well as the local Nusselt number. Because, the fluid absorbs the heat energy released due to radiation, and in turn releases the heat energy from the cylinder to the surrounding surface.
基金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.
基金supported by National Natural Science Foundation of China (Grant No.11172040)
文摘The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation(DES) method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.12204545 and 12274294)the Program for NUE independent research and development。
文摘We present an efficient approach to solve multi-dimensional time-dependent Schr?dinger equation(TDSE)in an intense laser field.In this approach,each spatial degree of freedom is treated as a distinguishable quasi-particle.The non-separable Coulomb potential is regarded as a two-body operator between different quasi-particles.The time-dependent variational principle is used to derive the equations of motion.Then the high-order multi-dimensional problem is broken down into several lower-order coupled equations,which can be efficiently solved.As a demonstration,we apply this method to solve the two-dimensional TDSE.The accuracy is tested by comparing the direct solutions of TDSE using several examples such as the strong-field ionization and the high harmonic generation.The results show that the present method is much more computationally efficient than the conventional one without sacrificing accuracy.The present method can be straightforwardly extended to three-dimensional problems.Our study provides a flexible method to investigate the laser-atom interaction in the nonperturbative regime.
文摘In this paper,the non-harmonic resonance of Bernoulli viscoelastic beams,Kirchhoff viscoelastic plates,Timoshenko viscoelastic beams,and Mindlin viscoelastic plates subjected to time-dependent exponentially decreasing transverse distributed load is investigated for the first time.The constitutive equations are expressed utilizing Boltzmann integral law with a constant bulk modulus.The displacement vector is approximated by employing the separation of variables method.The Laplace transformation is used to transfer equations from the time domain to the Laplace domain and vice versa.The novel point of the proposed method is to express,prove and calculate the critical time in which the displacement will be several times the displacement at time zero.In addition,this new method calculates the maximum deflection at the critical time,explicitly and exactly,without any need to follow the time-displacement curve with a low computational cost.Additionally,the proposed method introduces the critical range of time so that the responses are greater than the responses at time zero.
文摘This paper mainly studies the blowup phenomenon of solutions to the compressible Euler equations with general time-dependent damping for non-isentropic fluids in two and three space dimensions. When the initial data is assumed to be radially symmetric and the initial density contains vacuum, we obtain that classical solution, especially the density, will blow up on finite time. The results also reveal that damping can really delay the singularity formation.
基金supported by the Meteorological Special Project of China(GYHY200806005)the National Natural Sciences Foundation of China(40805028,40675039,40575036)the Key Technologies R&D Program of China(2009BAC51B04)
文摘On the basis of the quasi-geostrophic vorticity equation,theoretical research has been down upon the evolution of the amplitude of solitary Rossby waves employing the perturbation method,and come to the conclusion that the evolution of the amplitude satisfies the variable coefficient Korteweg-de Vries(KdV) equation.