The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficu...The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.展开更多
A new variable time step method,which is called the backwards calculating time step method,is presented in this paper.It allows numerical simulation of soil freezing and thawing while avoiding "phase change missi...A new variable time step method,which is called the backwards calculating time step method,is presented in this paper.It allows numerical simulation of soil freezing and thawing while avoiding "phase change missing and overflowing".A sensitive heat capacity model is introduced through which the calculation errors are analyzed.Then the equation using the self-adjusted time step is presented and solved using finite differences.Through this equation,the time needed for a space cell to reach the phase change point temperature is calculated.Using this time allows the time step to be adjusted so that errors caused by "phase change missing and overflowing" are successfully eliminated.Above all,the obvious features of this method are an accelerated rate for adjusting the time step and simplifing the computations.An actual example proves that this method can accurately calculate the temperature fields during soil freezing and thawing.It is an improvement over traditional methods and can be widely used on complicated multi-dimensional phase change problems.展开更多
The physical model based on heat transfer theory and virtual boundary method for analyzing unsteady thermal field of rotor plate for eddy current retarder used in automobile is established and boundary conditions are ...The physical model based on heat transfer theory and virtual boundary method for analyzing unsteady thermal field of rotor plate for eddy current retarder used in automobile is established and boundary conditions are also defined. The finite element governing equation is derived by Galerkin method. The time differential item is discrete based on Galerkin format that is stable at any condition. And a new style of varying time step method is used in iteration process. The thermal field on the rotor plate at the radial and axle directions is analyzed and varying temperature at appointed points on two side-surfaces is measured. The testing and analytical data are uniform approximately. Finite element method can be used for estimating thermal field of the rotor plate at initial design stage of eddy current retarder.展开更多
Efforts have been made to solve the Dirac equation with axially deformed scalar and vector WoodsSaxon potentials in the coordinate space with the imaginary time step method. The results of the singleparticle energies ...Efforts have been made to solve the Dirac equation with axially deformed scalar and vector WoodsSaxon potentials in the coordinate space with the imaginary time step method. The results of the singleparticle energies thus obtained are consistent with those calculated with the basis expansion method, which demonstrates the feasibility of the imaginary time step method for the relativistic static problems.展开更多
To reduce energy consumption on summer air conditioning,a novel seasonal soil cold storage mode using natural energy is presented and two-dimensional transient heat transfer model of U-tube is developed. The three pro...To reduce energy consumption on summer air conditioning,a novel seasonal soil cold storage mode using natural energy is presented and two-dimensional transient heat transfer model of U-tube is developed. The three processes of cold storage in winter,shut-down in transition season and cold extraction in summer are simulated by using sensitive heat capacity method with variable time step. The changing of U-tube outlet water temperature in different periods,daily cold storage and cold extraction are estimated. The temperature field of the U-tube and soil around the tube is investigated. Simulations show that seasonal soil cold storage using natural cold energy is feasible in the north to Changchun,which provides theoretical support for seasonal soil cold storage application.展开更多
The convergence for the Imaginary Time Step(ITS)evolution with time step is investigated by performing the ITS evolution for the Schrdinger-like equation and the charge-conjugate Schrdinger-like equation deduced from ...The convergence for the Imaginary Time Step(ITS)evolution with time step is investigated by performing the ITS evolution for the Schrdinger-like equation and the charge-conjugate Schrdinger-like equation deduced from Dirac equation for the single proton levels of 12C in both the Fermi and Dirac seas.For the guaranteed convergence of the ITS evolution to the"exact"results,the time step should be smaller than a"critical"time stepΔtc for a given single-particle level.The"critical"time stepΔtc is more sensitive to the quantum numbers|κ|than to the energy of the single-particle level.For the single-particle levels with the sameκ,their"critical"time steps are in the same order.For the single-particle levels with similar energy,a relatively small(large)"critical"time step for larger(smaller)|κ|is needed.These conclusions can be used in the future self-consistent calculation to optimize the evolution procedure.展开更多
Taking the single neutron levels of 12C in the Fermi sea as examples,the optimization of the imaginary time step(ITS)evolution with the box size and mesh size for the Dirac equation is investigated.For the weakly boun...Taking the single neutron levels of 12C in the Fermi sea as examples,the optimization of the imaginary time step(ITS)evolution with the box size and mesh size for the Dirac equation is investigated.For the weakly bound states,in order to reproduce the exact single-particle energies and wave functions,a relatively large box size is required.As long as the exact results can be reproduced,the ITS evolution with a smaller box size converges faster,while for both the weakly and deeply bound states,the ITS evolutions are less sensitive to the mesh size.Moreover,one can find a parabola relationship between the mesh size and the corresponding critical time step,i.e.,the largest time step to guarantee the convergence,which suggests that the ITS evolution with a larger mesh size allows larger critical time step,and thus can converge faster to the exact result.These conclusions are very helpful for optimizing the evolution procedure in the future self-consistent calculations.展开更多
Material point method（MPM）was originally introduced for large deformation problems in solid mechanics applications.Later,it has been successfully applied to solve a wide range of material behaviors.However,previous ...Material point method（MPM）was originally introduced for large deformation problems in solid mechanics applications.Later,it has been successfully applied to solve a wide range of material behaviors.However,previous research has indicated that MPM exhibits numerical instabilities when resolving incompressible flow problems.We study Chorin＇s projection method in MPM algorithm to simulate material incompressibility.Two projection-type schemes,non-incremental projection and incremental projection,are investigated for their accuracy and stability within MPM.Numerical examples show that the non-incremental projection scheme provides stable results in single phase MPM framework.Further,it avoids artificial pressure oscillations and small time steps that are present in the explicit MPM approach.展开更多
The Euler-Lagrange approach combined with a discrete element method has frequently been applied to elucidate the hydrodynamic behavior of dense fluid-solid flows in fluidized beds. In this work, the efficiency and acc...The Euler-Lagrange approach combined with a discrete element method has frequently been applied to elucidate the hydrodynamic behavior of dense fluid-solid flows in fluidized beds. In this work, the efficiency and accuracy of this model are investigated. Parameter studies are performed; in these studies, the stiffness coefficient, the fluid time step and the processor number are varied under conditions with different numbers of particles and different particle diameters. The obtained results are compared with measurements to derive the optimum parameters for CFD/DEM simulations. The results suggest that the application of higher stiffness coefficients slightly improves the simulation accuracy. However, the average computing time increases exponentially. At larger fluid time steps, the results show that the average computation time is independent of the applied fluid time step whereas the simulation accuracy decreases greatly with increasing the fluid time step. The use of smaller time steps leads to negligible improvements in the simulation accuracy but results in an exponential rise in the average computing time. The parallelization accelerates the DEM simulations if the critical number for the domain decomposition is not reached. Above this number, the performance is no longer proportional to the number of processors. The critical number for the domain decomposition depends on the number of particles. An increase in solid contents results in a shift of the critical decomposition number to higher numbers of CPUs.展开更多
基金financial support from Hunan Provincial Natura1 Science Foundation of China,Grant Number:02JJY2085,for this study
文摘The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.
基金Project 2006G1662-00 supported by the Key Science and Technology Project of Heilongjiang Province
文摘A new variable time step method,which is called the backwards calculating time step method,is presented in this paper.It allows numerical simulation of soil freezing and thawing while avoiding "phase change missing and overflowing".A sensitive heat capacity model is introduced through which the calculation errors are analyzed.Then the equation using the self-adjusted time step is presented and solved using finite differences.Through this equation,the time needed for a space cell to reach the phase change point temperature is calculated.Using this time allows the time step to be adjusted so that errors caused by "phase change missing and overflowing" are successfully eliminated.Above all,the obvious features of this method are an accelerated rate for adjusting the time step and simplifing the computations.An actual example proves that this method can accurately calculate the temperature fields during soil freezing and thawing.It is an improvement over traditional methods and can be widely used on complicated multi-dimensional phase change problems.
基金Department of Science and Technology of Jiangsu Province,China(No. BE2003-46).
文摘The physical model based on heat transfer theory and virtual boundary method for analyzing unsteady thermal field of rotor plate for eddy current retarder used in automobile is established and boundary conditions are also defined. The finite element governing equation is derived by Galerkin method. The time differential item is discrete based on Galerkin format that is stable at any condition. And a new style of varying time step method is used in iteration process. The thermal field on the rotor plate at the radial and axle directions is analyzed and varying temperature at appointed points on two side-surfaces is measured. The testing and analytical data are uniform approximately. Finite element method can be used for estimating thermal field of the rotor plate at initial design stage of eddy current retarder.
基金Supported by National Natural Science Foundation of China (10435010, 10775004, 10221003)Major State Basic Research Development Program (2007CB815000)
文摘Efforts have been made to solve the Dirac equation with axially deformed scalar and vector WoodsSaxon potentials in the coordinate space with the imaginary time step method. The results of the singleparticle energies thus obtained are consistent with those calculated with the basis expansion method, which demonstrates the feasibility of the imaginary time step method for the relativistic static problems.
基金Sponsored by Heilongjiang Province Emphasis Science and Technology Project ( Grant No 2006G1662-00)
文摘To reduce energy consumption on summer air conditioning,a novel seasonal soil cold storage mode using natural energy is presented and two-dimensional transient heat transfer model of U-tube is developed. The three processes of cold storage in winter,shut-down in transition season and cold extraction in summer are simulated by using sensitive heat capacity method with variable time step. The changing of U-tube outlet water temperature in different periods,daily cold storage and cold extraction are estimated. The temperature field of the U-tube and soil around the tube is investigated. Simulations show that seasonal soil cold storage using natural cold energy is feasible in the north to Changchun,which provides theoretical support for seasonal soil cold storage application.
基金partly supported by the Major State 973 Program(Grant No.2007CB815000)the National Natural Science Foundation of China(Grant Nos.10775004 and 10975008)
文摘The convergence for the Imaginary Time Step(ITS)evolution with time step is investigated by performing the ITS evolution for the Schrdinger-like equation and the charge-conjugate Schrdinger-like equation deduced from Dirac equation for the single proton levels of 12C in both the Fermi and Dirac seas.For the guaranteed convergence of the ITS evolution to the"exact"results,the time step should be smaller than a"critical"time stepΔtc for a given single-particle level.The"critical"time stepΔtc is more sensitive to the quantum numbers|κ|than to the energy of the single-particle level.For the single-particle levels with the sameκ,their"critical"time steps are in the same order.For the single-particle levels with similar energy,a relatively small(large)"critical"time step for larger(smaller)|κ|is needed.These conclusions can be used in the future self-consistent calculation to optimize the evolution procedure.
基金supported partially by Guizhou Science and Technology Foundation(Grant No J[2010]2135)the National Basic Research Program of China(Grant No 2007CB815000)the National Natural Science Foundation of China(Grant Nos 10775004,10947013,and 10975008)
文摘Taking the single neutron levels of 12C in the Fermi sea as examples,the optimization of the imaginary time step(ITS)evolution with the box size and mesh size for the Dirac equation is investigated.For the weakly bound states,in order to reproduce the exact single-particle energies and wave functions,a relatively large box size is required.As long as the exact results can be reproduced,the ITS evolution with a smaller box size converges faster,while for both the weakly and deeply bound states,the ITS evolutions are less sensitive to the mesh size.Moreover,one can find a parabola relationship between the mesh size and the corresponding critical time step,i.e.,the largest time step to guarantee the convergence,which suggests that the ITS evolution with a larger mesh size allows larger critical time step,and thus can converge faster to the exact result.These conclusions are very helpful for optimizing the evolution procedure in the future self-consistent calculations.
基金financially supported by the Cambridge Commonwealth Trust and the European Union’s Seventh Framework Programme 662 for research,Technological Development and Demonstration under Grant Agreement No.PIAP-GA-663 2012-324522 (MPM Dredge)
文摘Material point method（MPM）was originally introduced for large deformation problems in solid mechanics applications.Later,it has been successfully applied to solve a wide range of material behaviors.However,previous research has indicated that MPM exhibits numerical instabilities when resolving incompressible flow problems.We study Chorin＇s projection method in MPM algorithm to simulate material incompressibility.Two projection-type schemes,non-incremental projection and incremental projection,are investigated for their accuracy and stability within MPM.Numerical examples show that the non-incremental projection scheme provides stable results in single phase MPM framework.Further,it avoids artificial pressure oscillations and small time steps that are present in the explicit MPM approach.
文摘The Euler-Lagrange approach combined with a discrete element method has frequently been applied to elucidate the hydrodynamic behavior of dense fluid-solid flows in fluidized beds. In this work, the efficiency and accuracy of this model are investigated. Parameter studies are performed; in these studies, the stiffness coefficient, the fluid time step and the processor number are varied under conditions with different numbers of particles and different particle diameters. The obtained results are compared with measurements to derive the optimum parameters for CFD/DEM simulations. The results suggest that the application of higher stiffness coefficients slightly improves the simulation accuracy. However, the average computing time increases exponentially. At larger fluid time steps, the results show that the average computation time is independent of the applied fluid time step whereas the simulation accuracy decreases greatly with increasing the fluid time step. The use of smaller time steps leads to negligible improvements in the simulation accuracy but results in an exponential rise in the average computing time. The parallelization accelerates the DEM simulations if the critical number for the domain decomposition is not reached. Above this number, the performance is no longer proportional to the number of processors. The critical number for the domain decomposition depends on the number of particles. An increase in solid contents results in a shift of the critical decomposition number to higher numbers of CPUs.
