We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit,K-symplectic in the extended phase space with long time energy conservation properties. They are based on ex...We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit,K-symplectic in the extended phase space with long time energy conservation properties. They are based on extending the original phase space to several copies of the phase space and imposing a mechanical restraint on the copies of the phase space. Explicit K-symplectic methods are constructed for two non-canonical Hamiltonian systems. Numerical tests show that the proposed methods exhibit good numerical performance in preserving the phase orbit and the energy of the system over long time, whereas higher order Runge–Kutta methods do not preserve these properties. Numerical tests also show that the K-symplectic methods exhibit better efficiency than that of the same order implicit symplectic, explicit and implicit symplectic methods for the original nonseparable non-canonical systems. On the other hand, the fourth order K-symplectic method is more efficient than the fourth order Yoshida’s method, the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om explicit K-symplectic methods for the extended phase space Hamiltonians, but less efficient than the the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om extended phase space symplectic-like methods with the midpoint permutation.展开更多
Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are as...An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are assumed to occur at the two ends of a linear-elastic beam element.The auxiliary differential equations governing the plastic rotational displacements and their corresponding hysteretic displacements are replaced with linearized differential equations.Then,the two sets of equations of motion for the original nonlinear system can be reduced to an expanded-order equivalent linearized equation of motion for equivalent linear systems.To solve the equation of motion for equivalent linear systems,the nonstationary random vibration analysis is carried out based on the explicit time-domain method with high efficiency.Finally,the proposed treatment method for initial values of equivalent parameters is investigated in conjunction with parallel computing technology,which provides a new way of obtaining the equivalent linear systems at different time instants.Based on the explicit time-domain method,the key responses of interest of the converged equivalent linear system can be calculated through dimension reduction analysis with high efficiency.Numerical examples indicate that the proposed approach has high computational efficiency,and shows good applicability to weak nonlinear and medium-intensity nonlinear systems.展开更多
Taking CPU time cost and analysis accuracy into account, dynamic explicit finite ele- ment method is adopted to optimize the forming process of autobody panels that often have large sizes and complex geometry. In this...Taking CPU time cost and analysis accuracy into account, dynamic explicit finite ele- ment method is adopted to optimize the forming process of autobody panels that often have large sizes and complex geometry. In this paper, for the sake of illustrating in detail how dynamic explicit finite element method is applied to the numerical simulation of the autobody panel forming process,an example of optimization of stamping process pain meters of an inner door panel is presented. Using dynamic explicit finite element code Ls-DYNA3D, the inner door panel has been optimized by adapting pa- rameters such as the initial blank geometry and position, blank-holder forces and the location of drawbeads, and satisfied results are obtained.展开更多
A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a...A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a structure does not need to be considered. This is because this subfamily is unconditionally stable for any instantaneous stiffness softening system, linear elastic system and instantaneous stiffness hardening system that might occur in the pseudodynamic testing of a real structure. In addition, it also offers good accuracy when compared to a general second-order accurate method for both linear elastic and nonlinear systems.展开更多
A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditiona...A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.展开更多
Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow pheno...Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.展开更多
It has been proven that the implicit method used to solve the vibration equation can be transformed into an explicit method,which is called the concomitant explicit method.The constant acceleration method's concom...It has been proven that the implicit method used to solve the vibration equation can be transformed into an explicit method,which is called the concomitant explicit method.The constant acceleration method's concomitant explicit method was used as an example and is described in detail in this paper.The relationship between the implicit method and explicit method is defined,which provides some guidance about how to create a new explicit method that has high precision and computational efficiency.展开更多
We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusi...We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.展开更多
In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on ...In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium, elastic single-phase medium and ideal fluid medium. This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling, and which does not need to solve sys-tem of linear equations. The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water, dam, sediment and basement rock. The vertical displacement at the top point of the dam is calculated and some conclusions are given.展开更多
The paper starts with a brief overview to the necessity of sheet metal forming simulation and the complexity of automobile panel forming, then leads to finite element analysis (FEA) which is a powerful simulation too...The paper starts with a brief overview to the necessity of sheet metal forming simulation and the complexity of automobile panel forming, then leads to finite element analysis (FEA) which is a powerful simulation tool for analyzing complex three-dimensional sheet metal forming problems. The theory and features of the dynamic explicit finite element methods are introduced and the available various commercial finite element method codes used for sheet metal forming simulation in the world are discussed,and the civil and international status quo of automobile panel simulation as well. The front door outer panel of one certain new automobile is regarded as one example that the dynamic explicit FEM code Dynaform is used for the simulation of the front door outer panel forming process. Process defects such as ruptures are predicted. The improving methods can be given according to the simulation results. Foreground of sheet metal forming simulation is outlined.展开更多
Time domain dynamic analysis of inclined dam-reservoir-foundation interaction was conducted using finite difference method （FDM）. The Timoshenko beam theory and the Euler-Bemoulli beam theory were implemented to dra...Time domain dynamic analysis of inclined dam-reservoir-foundation interaction was conducted using finite difference method （FDM）. The Timoshenko beam theory and the Euler-Bemoulli beam theory were implemented to draw out governing equation of beam. The interactions between the dam and the soil were modeled by using a translational spring and a rotational spring. A Sommerfeld＇s radiation condition at the infinity boundary of the fluid domain was adopted. The effects of the reservoir bottom absorption and surface waves on the dam-reservoir-foundation interaction due to the earthquake were studied. To avoid the instability of solution, a semi-implicit scheme was used for the discretization of the governing equation of dam and an explicit scheme was used for the discretization of the governing equation of fluid. The results show that as the slope of upstream dam increases, the hydrodynamic pressure on the dam is reduced. Moreover, when the Timoshenko beam theory is used, the system response increases.展开更多
By performing a particular spatial discretization to the nonlinear Schrodinger equation(NLSE),we obtain a non-integrable Hamiltonian system which can be decomposed into three integrable parts(L-L-N splitting).We integ...By performing a particular spatial discretization to the nonlinear Schrodinger equation(NLSE),we obtain a non-integrable Hamiltonian system which can be decomposed into three integrable parts(L-L-N splitting).We integrate each part by calculating its phase flow,and develop explicit symplectic integrators of different orders for the original Hamiltonian by composing the phase flows.A 2nd-order reversible constructed symplectic scheme is employed to simulate solitons motion and invariants behavior of the NLSE.The simulation results are compared with a 3rd-order non-symplectic implicit Runge-Kutta method,and the convergence of the formal energy of this symplectic integrator is also verified.The numerical results indicate that the explicit symplectic scheme obtained via L-L-N splitting is an effective numerical tool for solving the NLSE.展开更多
In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complic...In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.展开更多
In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum alo...In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum along energy are analyzed with suitable transformations. For numerical calculation, an implicit finite difference method is applied to solve a set of nonlinear dimensionless partial differential equations. Dimensionless velocity and temperature profile are also investigated due to the effects of assumed parameters in the concerned problem. An explicit finite difference technique is used to compute velocity and temperature profiles. The stability conditions are also examined.展开更多
In this paper,we analyze the explicit Runge-Kutta discontinuous Galerkin(RKDG)methods for the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology.The RKDG...In this paper,we analyze the explicit Runge-Kutta discontinuous Galerkin(RKDG)methods for the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology.The RKDG methods use a third order explicit total-variation-diminishing Runge-Kutta(TVDRK3)time discretization and upwinding numerical fluxes.By using the energy method,under a standard CourantFriedrichs-Lewy(CFL)condition,we obtain L2stability for general solutions and a priori error estimates when the solutions are smooth enough.The theoretical results are proved for piecewise polynomials with any degree k 1.Finally,since the solutions to this system are non-negative,we discuss a positivity-preserving limiter to preserve positivity without compromising accuracy.Numerical results are provided to demonstrate these RKDG methods.展开更多
In this paper,we propose two new explicit multi-symplectic splitting methods for the nonlinear Dirac(NLD)equation.Based on its multi-symplectic formulation,the NLD equation is split into one linear multi-symplectic sy...In this paper,we propose two new explicit multi-symplectic splitting methods for the nonlinear Dirac(NLD)equation.Based on its multi-symplectic formulation,the NLD equation is split into one linear multi-symplectic system and one nonlinear infinite Hamiltonian system.Then multi-symplectic Fourier pseudospectral method and multi-symplectic Preissmann scheme are employed to discretize the linear subproblem,respectively.And the nonlinear subsystem is solved by a symplectic scheme.Finally,a composition method is applied to obtain the final schemes for the NLD equation.We find that the two proposed schemes preserve the total symplecticity and can be solved explicitly.Numerical experiments are presented to show the effectiveness of the proposed methods.展开更多
Broken gap is an extremely dangerous state in the service of high-speed rails,and the violent wheel–rail impact forces will be intensified when a vehicle passes the gap at high speeds,which may cause a secondary frac...Broken gap is an extremely dangerous state in the service of high-speed rails,and the violent wheel–rail impact forces will be intensified when a vehicle passes the gap at high speeds,which may cause a secondary fracture to rail and threaten the running safety of the vehicle.To recognize the damage tolerance of rail fracture length,the implicit–explicit sequential approach is adopted to simulate the wheel–rail high-frequency impact,which considers the factors such as the coupling effect between frictional contact and structural vibration,nonlinear material and real geometric profile.The results demonstrate that the plastic deformation and stress are distributed in crescent shape during the impact at the back rail end,increasing with the rail fracture length.The axle box acceleration in the frequency domain displays two characteristic modes with frequencies around 1,637 and 404 Hz.The limit of the rail fracture length is 60 mm for high-speed railway at a speed of 250 km/h.展开更多
By using the Feynman-Kac formula and combining with Itˆo-Taylor expansion and finite difference approximation,we first develop an explicit third order onestep method for solving decoupled forward backward stochastic d...By using the Feynman-Kac formula and combining with Itˆo-Taylor expansion and finite difference approximation,we first develop an explicit third order onestep method for solving decoupled forward backward stochastic differential equations.Then based on the third order one,an explicit fourth order method is further proposed.Several numerical tests are also presented to illustrate the stability and high order accuracy of the proposed methods.展开更多
A model for the amount of pollution in lakes connected with some rivers is introduced.In this model,it is supposed the density of pollution in a lake has memory.The model leads to a system of fractional differential e...A model for the amount of pollution in lakes connected with some rivers is introduced.In this model,it is supposed the density of pollution in a lake has memory.The model leads to a system of fractional differential equations.This system is transformed into a system of Volterra integral equations with memory kernels.The existence and regularity of the solutions are investigated.A high-order numerical method is introduced and analyzed and compared with an explicit method based on the regularity of the solution.Validation examples are supported,and some models are simulated and discussed.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 11901564 and 12171466)。
文摘We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit,K-symplectic in the extended phase space with long time energy conservation properties. They are based on extending the original phase space to several copies of the phase space and imposing a mechanical restraint on the copies of the phase space. Explicit K-symplectic methods are constructed for two non-canonical Hamiltonian systems. Numerical tests show that the proposed methods exhibit good numerical performance in preserving the phase orbit and the energy of the system over long time, whereas higher order Runge–Kutta methods do not preserve these properties. Numerical tests also show that the K-symplectic methods exhibit better efficiency than that of the same order implicit symplectic, explicit and implicit symplectic methods for the original nonseparable non-canonical systems. On the other hand, the fourth order K-symplectic method is more efficient than the fourth order Yoshida’s method, the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om explicit K-symplectic methods for the extended phase space Hamiltonians, but less efficient than the the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om extended phase space symplectic-like methods with the midpoint permutation.
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
基金Fundamental Research Funds for the Central Universities under Grant No.2682022CX072the Research and Development Plan in Key Areas of Guangdong Province under Grant No.2020B0202010008。
文摘An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are assumed to occur at the two ends of a linear-elastic beam element.The auxiliary differential equations governing the plastic rotational displacements and their corresponding hysteretic displacements are replaced with linearized differential equations.Then,the two sets of equations of motion for the original nonlinear system can be reduced to an expanded-order equivalent linearized equation of motion for equivalent linear systems.To solve the equation of motion for equivalent linear systems,the nonstationary random vibration analysis is carried out based on the explicit time-domain method with high efficiency.Finally,the proposed treatment method for initial values of equivalent parameters is investigated in conjunction with parallel computing technology,which provides a new way of obtaining the equivalent linear systems at different time instants.Based on the explicit time-domain method,the key responses of interest of the converged equivalent linear system can be calculated through dimension reduction analysis with high efficiency.Numerical examples indicate that the proposed approach has high computational efficiency,and shows good applicability to weak nonlinear and medium-intensity nonlinear systems.
文摘Taking CPU time cost and analysis accuracy into account, dynamic explicit finite ele- ment method is adopted to optimize the forming process of autobody panels that often have large sizes and complex geometry. In this paper, for the sake of illustrating in detail how dynamic explicit finite element method is applied to the numerical simulation of the autobody panel forming process,an example of optimization of stamping process pain meters of an inner door panel is presented. Using dynamic explicit finite element code Ls-DYNA3D, the inner door panel has been optimized by adapting pa- rameters such as the initial blank geometry and position, blank-holder forces and the location of drawbeads, and satisfied results are obtained.
基金Science Council, Chinese Taipei Under Grant No. NSC-95-2221-E-027-099
文摘A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a structure does not need to be considered. This is because this subfamily is unconditionally stable for any instantaneous stiffness softening system, linear elastic system and instantaneous stiffness hardening system that might occur in the pseudodynamic testing of a real structure. In addition, it also offers good accuracy when compared to a general second-order accurate method for both linear elastic and nonlinear systems.
基金Project(Z110803)supported by the State Key Laboratory of Geomechanics and Geotechnical Engineering,ChinaProject(2008AA062303)supported by the National High Technology Research and Development Program of China
文摘A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.
基金King Mongkut’s University of Technology North Bangkok (KMUTNB)the Office of the Higher Education Commission (OHEC)the National Metal and Materials Technology Center (MTEC) for supporting this research work
文摘Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.
基金Fundamental Research Funds for the Central Universities
文摘It has been proven that the implicit method used to solve the vibration equation can be transformed into an explicit method,which is called the concomitant explicit method.The constant acceleration method's concomitant explicit method was used as an example and is described in detail in this paper.The relationship between the implicit method and explicit method is defined,which provides some guidance about how to create a new explicit method that has high precision and computational efficiency.
基金Open Access funding enabled and organized by Projekt DEAL.
文摘We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.
基金National Natural Scienccs Foundation of China (50178005).
文摘In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium, elastic single-phase medium and ideal fluid medium. This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling, and which does not need to solve sys-tem of linear equations. The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water, dam, sediment and basement rock. The vertical displacement at the top point of the dam is calculated and some conclusions are given.
文摘The paper starts with a brief overview to the necessity of sheet metal forming simulation and the complexity of automobile panel forming, then leads to finite element analysis (FEA) which is a powerful simulation tool for analyzing complex three-dimensional sheet metal forming problems. The theory and features of the dynamic explicit finite element methods are introduced and the available various commercial finite element method codes used for sheet metal forming simulation in the world are discussed,and the civil and international status quo of automobile panel simulation as well. The front door outer panel of one certain new automobile is regarded as one example that the dynamic explicit FEM code Dynaform is used for the simulation of the front door outer panel forming process. Process defects such as ruptures are predicted. The improving methods can be given according to the simulation results. Foreground of sheet metal forming simulation is outlined.
文摘Time domain dynamic analysis of inclined dam-reservoir-foundation interaction was conducted using finite difference method （FDM）. The Timoshenko beam theory and the Euler-Bemoulli beam theory were implemented to draw out governing equation of beam. The interactions between the dam and the soil were modeled by using a translational spring and a rotational spring. A Sommerfeld＇s radiation condition at the infinity boundary of the fluid domain was adopted. The effects of the reservoir bottom absorption and surface waves on the dam-reservoir-foundation interaction due to the earthquake were studied. To avoid the instability of solution, a semi-implicit scheme was used for the discretization of the governing equation of dam and an explicit scheme was used for the discretization of the governing equation of fluid. The results show that as the slope of upstream dam increases, the hydrodynamic pressure on the dam is reduced. Moreover, when the Timoshenko beam theory is used, the system response increases.
基金This research is partially supported by the Informatization Construction of Knowledge Innovation Projects of the Chinese Academy of Sciences“Supercomputing En-vironment Construction and Application”(INF105-SCE)National Natural Science Foundation of China(Grant Nos.10471145 and 10672143).
文摘By performing a particular spatial discretization to the nonlinear Schrodinger equation(NLSE),we obtain a non-integrable Hamiltonian system which can be decomposed into three integrable parts(L-L-N splitting).We integrate each part by calculating its phase flow,and develop explicit symplectic integrators of different orders for the original Hamiltonian by composing the phase flows.A 2nd-order reversible constructed symplectic scheme is employed to simulate solitons motion and invariants behavior of the NLSE.The simulation results are compared with a 3rd-order non-symplectic implicit Runge-Kutta method,and the convergence of the formal energy of this symplectic integrator is also verified.The numerical results indicate that the explicit symplectic scheme obtained via L-L-N splitting is an effective numerical tool for solving the NLSE.
基金National Natural Science Foundation of China Under Grant No. 50178065
文摘In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.
文摘In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum along energy are analyzed with suitable transformations. For numerical calculation, an implicit finite difference method is applied to solve a set of nonlinear dimensionless partial differential equations. Dimensionless velocity and temperature profile are also investigated due to the effects of assumed parameters in the concerned problem. An explicit finite difference technique is used to compute velocity and temperature profiles. The stability conditions are also examined.
基金supported by the University of Science and Technology of China Special Grant for Postgraduate ResearchInnovation and Practice+5 种基金the Chinese Academy of Science Special Grant for Postgraduate ResearchInnovation and PracticeDepartment of Energy of USA(Grant No.DE-FG02-08ER25863)National Science Foundation of USA(Grant No.DMS-1112700)National Natural Science Foundation of China(Grant Nos.1107123491130016 and 91024025)
文摘In this paper,we analyze the explicit Runge-Kutta discontinuous Galerkin(RKDG)methods for the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology.The RKDG methods use a third order explicit total-variation-diminishing Runge-Kutta(TVDRK3)time discretization and upwinding numerical fluxes.By using the energy method,under a standard CourantFriedrichs-Lewy(CFL)condition,we obtain L2stability for general solutions and a priori error estimates when the solutions are smooth enough.The theoretical results are proved for piecewise polynomials with any degree k 1.Finally,since the solutions to this system are non-negative,we discuss a positivity-preserving limiter to preserve positivity without compromising accuracy.Numerical results are provided to demonstrate these RKDG methods.
基金the open foundations of State Key Laboratory of High Performance Computing and State Key Laboratory of Aerodynamics.Y.C.gratefully acknowledges support from NUDT’s Innovation Foundation(Grant No.B110205)H.Z.was supported by the Natural Science Foundation of China(Grant No.11301525).
文摘In this paper,we propose two new explicit multi-symplectic splitting methods for the nonlinear Dirac(NLD)equation.Based on its multi-symplectic formulation,the NLD equation is split into one linear multi-symplectic system and one nonlinear infinite Hamiltonian system.Then multi-symplectic Fourier pseudospectral method and multi-symplectic Preissmann scheme are employed to discretize the linear subproblem,respectively.And the nonlinear subsystem is solved by a symplectic scheme.Finally,a composition method is applied to obtain the final schemes for the NLD equation.We find that the two proposed schemes preserve the total symplecticity and can be solved explicitly.Numerical experiments are presented to show the effectiveness of the proposed methods.
基金The work is supported by the National Natural Science Foundation of China(Nos.51608459,51778542 and U1734207)Fundamental Research Funds for the Central Universities(No.2682018CX01)Cultivation Program for the Excellent Doctoral Dissertation of Southwest Jiaotong University.
文摘Broken gap is an extremely dangerous state in the service of high-speed rails,and the violent wheel–rail impact forces will be intensified when a vehicle passes the gap at high speeds,which may cause a secondary fracture to rail and threaten the running safety of the vehicle.To recognize the damage tolerance of rail fracture length,the implicit–explicit sequential approach is adopted to simulate the wheel–rail high-frequency impact,which considers the factors such as the coupling effect between frictional contact and structural vibration,nonlinear material and real geometric profile.The results demonstrate that the plastic deformation and stress are distributed in crescent shape during the impact at the back rail end,increasing with the rail fracture length.The axle box acceleration in the frequency domain displays two characteristic modes with frequencies around 1,637 and 404 Hz.The limit of the rail fracture length is 60 mm for high-speed railway at a speed of 250 km/h.
基金supported by the NSF of China(No.12001539)the NSF of Hunan Province(No.2020JJ5647)China Postdoctoral Science Foundation(No.2019TQ0073).
文摘By using the Feynman-Kac formula and combining with Itˆo-Taylor expansion and finite difference approximation,we first develop an explicit third order onestep method for solving decoupled forward backward stochastic differential equations.Then based on the third order one,an explicit fourth order method is further proposed.Several numerical tests are also presented to illustrate the stability and high order accuracy of the proposed methods.
文摘A model for the amount of pollution in lakes connected with some rivers is introduced.In this model,it is supposed the density of pollution in a lake has memory.The model leads to a system of fractional differential equations.This system is transformed into a system of Volterra integral equations with memory kernels.The existence and regularity of the solutions are investigated.A high-order numerical method is introduced and analyzed and compared with an explicit method based on the regularity of the solution.Validation examples are supported,and some models are simulated and discussed.
