A decoupling method with different subdomain time steps for the non-stationary Navier- Stokes/Darcy model is formulated and analyzed. The method has asynchronous time steps, which adopt small time steps in the fluid r...A decoupling method with different subdomain time steps for the non-stationary Navier- Stokes/Darcy model is formulated and analyzed. The method has asynchronous time steps, which adopt small time steps in the fluid region and large time steps in the porous region. It saves relatively large amount of CPU time. Stability and convergence of the method are proved. The numerical results are presented to illustrate the features of the proposed method.展开更多
In the fast-evolving landscape of digital networks,the incidence of network intrusions has escalated alarmingly.Simultaneously,the crucial role of time series data in intrusion detection remains largely underappreciat...In the fast-evolving landscape of digital networks,the incidence of network intrusions has escalated alarmingly.Simultaneously,the crucial role of time series data in intrusion detection remains largely underappreciated,with most systems failing to capture the time-bound nuances of network traffic.This leads to compromised detection accuracy and overlooked temporal patterns.Addressing this gap,we introduce a novel SSAE-TCN-BiLSTM(STL)model that integrates time series analysis,significantly enhancing detection capabilities.Our approach reduces feature dimensionalitywith a Stacked Sparse Autoencoder(SSAE)and extracts temporally relevant features through a Temporal Convolutional Network(TCN)and Bidirectional Long Short-term Memory Network(Bi-LSTM).By meticulously adjusting time steps,we underscore the significance of temporal data in bolstering detection accuracy.On the UNSW-NB15 dataset,ourmodel achieved an F1-score of 99.49%,Accuracy of 99.43%,Precision of 99.38%,Recall of 99.60%,and an inference time of 4.24 s.For the CICDS2017 dataset,we recorded an F1-score of 99.53%,Accuracy of 99.62%,Precision of 99.27%,Recall of 99.79%,and an inference time of 5.72 s.These findings not only confirm the STL model’s superior performance but also its operational efficiency,underpinning its significance in real-world cybersecurity scenarios where rapid response is paramount.Our contribution represents a significant advance in cybersecurity,proposing a model that excels in accuracy and adaptability to the dynamic nature of network traffic,setting a new benchmark for intrusion detection systems.展开更多
The matrix version of Symmetric Successive Over Relaxation(matrix-SSOR)scheme has been proved to be more efficient than the standard Lower-Upper Symmetric Gauss-Seidel(LUSGS),but less robust for high-speed flows.In or...The matrix version of Symmetric Successive Over Relaxation(matrix-SSOR)scheme has been proved to be more efficient than the standard Lower-Upper Symmetric Gauss-Seidel(LUSGS),but less robust for high-speed flows.In order to ulteriorly improve the convergence rate as well as numerical stability of matrix-SSOR,two improvements regarding entropy fix and local time step have been proposed and validated.Firstly,an augmented entropy fix method is imposed on the inviscid Jacobian matrix and proved to be effective in two high-speed flows,in which the key parameter in entropy fix is discussed and found to be insensitive within appropriate range of values.Since the time step also has great effects on the numerical stability and convergence rate,a modified cell residual adapted local time step method with consideration of the residual history is developed,which is found to be effective for increasing the convergence rate when the matrix-SSOR is applied,but invalid when the LU-SGS is used.The proposed modified local time step method is also insensitive to the key parameter within appropriate range of values.The two modifications can be conveniently implanted into analogous matrix-type implicit schemes to improve the numerical performance.展开更多
In order to further achieve the balance between the calculation accuracy and efficiency of the transient analysis of the aero-engine disc cavity system,an Optimized Time-adaptive Aerother-mal Coupling calculation(OTAC...In order to further achieve the balance between the calculation accuracy and efficiency of the transient analysis of the aero-engine disc cavity system,an Optimized Time-adaptive Aerother-mal Coupling calculation(OTAC)method has been proposed.It combines one-dimensional tran-sient calculation of air system,Conventional Sequence Staggered(CSS)method,Time-adaptive Aerothermal Coupling calculation(TAC)method and differential evolution optimization algorithm to obtain an efficient and high-precision aerothermal coupling calculation method of air system.Considering both the heat conduction in the solid domain and the flow in the fluid domain as unsteady states in the OTAC,the interaction of fluid-solid information within a single coupling time step size was implemented based on the CSS method.Furthermore,the coupling time step size was automatically adjusted with the number of iterations by using the Proportional-Integral-Deri vative(PID)controller.Results show that when compared with the traditional loosely coupling method with a fixed time step size,the computational accuracy and efficiency of the OTAC method are improved by 8.9%and 30%,respectively.Compared with the tight coupling calculation,the OTAC method can achieve a speedup of 1 to 2 orders of magnitude,while the calculation error is maintained within 6.1%.展开更多
A growing interest has been devoted to the contra-rotating propellers （CRPs） due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibrati...A growing interest has been devoted to the contra-rotating propellers （CRPs） due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibration. Compared with the single-screw system, it is more difficult for the open water performance prediction because forward and aft propellers interact with each other and generate a more complicated flow field around the CRPs system. The current work focuses on the open water performance prediction of contra-rotating propellers by RANS and sliding mesh method considering the effect of computational time step size and turbulence model. The validation study has been performed on two sets of contra-rotating propellers developed by David W Taylor Naval Ship R ＆ D center. Compared with the experimental data, it shows that RANS with sliding mesh method and SST k-ω turbulence model has a good precision in the open water performance prediction of contra-rotating propellers, and small time step size can improve the level of accuracy for CRPs with the same blade number of forward and aft propellers, while a relatively large time step size is a better choice for CRPs with different blade numbers.展开更多
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.展开更多
There are two models in use today to analyze structural responses when subjected to earthquake ground motions, the Displacement Input Model （DIM） and the Acceleration Input Model （AIM）. The time steps used in dire...There are two models in use today to analyze structural responses when subjected to earthquake ground motions, the Displacement Input Model （DIM） and the Acceleration Input Model （AIM）. The time steps used in direct integration methods for these models are analyzed to examine the suitability of DIM. Numerical results are presented and show that the time-step for DIM is about the same as for AIM, and achieves the same accuracy. This is contrary to previous research that reported that there are several sources of numerical errors associated with the direct application of earthquake displacement loading, and a very small time step is required to define the displacement record and to integrate the dynamic equilibrium equation. It is shown in this paper that DIM is as accurate and suitable as, if not more than, AIM for analyzing the response of a structure to uniformly distributed and spatially varying ground motions.展开更多
Errors due to split time stepping are discussed for an explicit free–surface ocean model. In commonly used split time stepping, the way of time integration for the barotropic momentum equation is not compatible with ...Errors due to split time stepping are discussed for an explicit free–surface ocean model. In commonly used split time stepping, the way of time integration for the barotropic momentum equation is not compatible with that of the baroclinic one. The baroclinic equation has three–time–level structure because of leapfrog scheme. The barotropic one, however, has two–time–level structure when represented in terms of the baroclinic time level, on which the baroclinic one is integrated. This incompatibility results in the splitting errors as shown in this paper. The proper split time stepping is therefore proposed in such a way that the compatibility is kept between the barotropic and baroclinic equations. Its splitting errors are shown extremely small, so that it is particularly relevant to long–term integration for climate studies. It is applied to a free–surface model for the North Pacific Ocean.展开更多
This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditiona...This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.展开更多
A multiple time step algorithm, called reversible reference system propagator algorithm, is introduced for the long time molecular dynamics simulation. In contrast to the conventional algorithms, the multiple time met...A multiple time step algorithm, called reversible reference system propagator algorithm, is introduced for the long time molecular dynamics simulation. In contrast to the conventional algorithms, the multiple time method has better convergence, stability and efficiency. The method is validated by simulating free relaxation and the hypervelocity impact of nano-clusters. The time efficiency of the multiple time step method enables us to investigate the long time interaction between lattice dislocations and low-angle grain boundaries.展开更多
The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integ...The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.展开更多
With the cell vertex finite volume discretization in space and second order backward implicit discretization in time, 2D unsteady Navier Stokes equations are solved by a dual time stepping method to simulate compr...With the cell vertex finite volume discretization in space and second order backward implicit discretization in time, 2D unsteady Navier Stokes equations are solved by a dual time stepping method to simulate compressible viscous flow around rigid airfoils in arbitrary unsteady motion. The selection of physical time step is not restricted by stability condition any more, and most of the successful acceleration techniques used in steady calculations can be implemented to increase the computation efficiency.展开更多
The heat transfer during the casting solidification process includes the heat radiation of the high temperature casting and the mold,the heat convection between the casting and the mold,and the heat conduction inside ...The heat transfer during the casting solidification process includes the heat radiation of the high temperature casting and the mold,the heat convection between the casting and the mold,and the heat conduction inside the casting and from the casting to the mold. In this paper,a formula of time step in simulation of solidification is derived,considering the heat radiation,convection and conduction based on the conservation of energy. The different heat transfer conditions between the conventional sand casting and the permanent mold casting are taken into account in this formula. The characteristics of heat transfer in the interior and surface of the casting are also considered. The numerical experiments show that this formula can avoid computational dispersion,and improve the computational efficiency by about 20% in the simulation of solidification process.展开更多
This paper focuses on the study of the stability of explicit time integration algorithm for dynamic problem by the Extended Finite Element Method(XFEM).A new enrichment scheme of crack tip is proposed within the frame...This paper focuses on the study of the stability of explicit time integration algorithm for dynamic problem by the Extended Finite Element Method(XFEM).A new enrichment scheme of crack tip is proposed within the framework of XFEM.Then the governing equations are derived and evolved into the discretized form.For dynamic problem,the lumped mass and the explicit time algorithm are applied.With different grid densities and different forms of Newmark scheme,the Dynamic Stress Intensity Factor(DSIF)is computed by using interaction integral approach to reflect the dynamic response.The effectiveness of the proposed scheme is demonstrated through the numerical examples,and the critical time stepping in different situations are listed and analyzed to illustrate the factors that affect the numerical stability.展开更多
A self-adaptive precise algorithm in the time domain was employed to solve 2-D nonlinear coupled heat and moisture transfer problems. By expanding variables at a discretized time interval, the variations of variables ...A self-adaptive precise algorithm in the time domain was employed to solve 2-D nonlinear coupled heat and moisture transfer problems. By expanding variables at a discretized time interval, the variations of variables can be described more precisely,and a nonlinear coupled initial and boundary value problem was converted into a series of recurrent linear boundary value problems which are solved by FE technique. In the computation, no additional assumption and the nonlinear iteration are required, and a criterion for self-adaptive computation is proposed to maintain sufficient computing accuracy for the change sizes of time steps. In the numerical comparison, the variations of material properties with temperature, moisture content, and both temperature and moisture content are taken into account, respectively. Satisfactory results have been obtained, indicating that the proposed approach is capable of dealing with complex nonlinear problems.展开更多
In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entrop...In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entropy solution in the limit if Courant number is less than or equal to 1.展开更多
A natural generalization of random choice finite difference scheme of Harten and Lax for Courant number larger than 1 is obtained. We handle interactions between neighboring Riemann solvers by linear superposition of ...A natural generalization of random choice finite difference scheme of Harten and Lax for Courant number larger than 1 is obtained. We handle interactions between neighboring Riemann solvers by linear superposition of their conserved quantities. We show consistency of the scheme for arbitrarily large Courant numbers. For scalar problems the scheme is total variation diminishing.A brief discussion is given for entropy condition.展开更多
A positive leap second will be introduced in the UTC time scale UTC(JATC)、UTC (CSAO) and UTC time signals of BPL、BPM transmittings at the end of June 1997. The sequence of dates of the UTC second markers will be:
This paper presents two comparisons or tests for a Lagrangian model of zooplankton dispersion:numerical schemes and time steps.Firstly,we compared three numerical schemes using idealized circulations.Results show that...This paper presents two comparisons or tests for a Lagrangian model of zooplankton dispersion:numerical schemes and time steps.Firstly,we compared three numerical schemes using idealized circulations.Results show that the precisions of the advanced Adams-Bashfold-Moulton(ABM) method and the Runge-Kutta(RK) method were in the same order and both were much higher than that of the Euler method.Furthermore,the advanced ABM method is more efficient than the RK method in computational memory requirements and time consumption.We therefore chose the advanced ABM method as the Lagrangian particle-tracking algorithm.Secondly,we performed a sensitivity test for time steps,using outputs of the hydrodynamic model,Symphonie.Results show that the time step choices depend on the fluid response time that is related to the spatial resolution of velocity fields.The method introduced by Oliveira et al.in 2002 is suitable for choosing time steps of Lagrangian particle-tracking models,at least when only considering advection.展开更多
The long-term stability of the roof is particularly important in designing underground rock structures.To estimate the durability of roof strata in underground excavation,a computation scheme of subcritical crack grow...The long-term stability of the roof is particularly important in designing underground rock structures.To estimate the durability of roof strata in underground excavation,a computation scheme of subcritical crack growth is proposed in this study.By adopting the proposed method,the potential collapse location of strata is derivable in accordance with a static model,the durability of roof strata can be estimated,a dynamic time step control strategy is achieved to balance the accuracy and speed of computing,and the initial crack size of rock can be estimated.In addition to the above,a mechanical model of underground excavation with non-uniformly distributed loads and partially yielded foundation is presented as the prototypical case.A set of case studies is carried out that showcase a power correlation between applied stress and roof durability.The allowable applied tensile stress for a 100-year life cycle is about 76%of the tensile strength.By using the proposed subcritical crack growth computation scheme,the roof stability in an underground excavation can be identified not only from the spatial view but also from the temporal perspective.展开更多
基金The authors would like to thank the referees for helpful comments and suggestions, which lead to substantial improvements of the presentation.The work was Supported by the National Nature Foundation of China （NO.11401422）, the Provincial Natural Science Foundation of Shanxi （NO.2015011001, 2014011005-4）,China Postdoctoral Science Foundation funded project （No.2015M582336）
文摘A decoupling method with different subdomain time steps for the non-stationary Navier- Stokes/Darcy model is formulated and analyzed. The method has asynchronous time steps, which adopt small time steps in the fluid region and large time steps in the porous region. It saves relatively large amount of CPU time. Stability and convergence of the method are proved. The numerical results are presented to illustrate the features of the proposed method.
基金supported in part by the Gansu Province Higher Education Institutions Industrial Support Program:Security Situational Awareness with Artificial Intelligence and Blockchain Technology.Project Number(2020C-29).
文摘In the fast-evolving landscape of digital networks,the incidence of network intrusions has escalated alarmingly.Simultaneously,the crucial role of time series data in intrusion detection remains largely underappreciated,with most systems failing to capture the time-bound nuances of network traffic.This leads to compromised detection accuracy and overlooked temporal patterns.Addressing this gap,we introduce a novel SSAE-TCN-BiLSTM(STL)model that integrates time series analysis,significantly enhancing detection capabilities.Our approach reduces feature dimensionalitywith a Stacked Sparse Autoencoder(SSAE)and extracts temporally relevant features through a Temporal Convolutional Network(TCN)and Bidirectional Long Short-term Memory Network(Bi-LSTM).By meticulously adjusting time steps,we underscore the significance of temporal data in bolstering detection accuracy.On the UNSW-NB15 dataset,ourmodel achieved an F1-score of 99.49%,Accuracy of 99.43%,Precision of 99.38%,Recall of 99.60%,and an inference time of 4.24 s.For the CICDS2017 dataset,we recorded an F1-score of 99.53%,Accuracy of 99.62%,Precision of 99.27%,Recall of 99.79%,and an inference time of 5.72 s.These findings not only confirm the STL model’s superior performance but also its operational efficiency,underpinning its significance in real-world cybersecurity scenarios where rapid response is paramount.Our contribution represents a significant advance in cybersecurity,proposing a model that excels in accuracy and adaptability to the dynamic nature of network traffic,setting a new benchmark for intrusion detection systems.
基金supported by the National Natural Science Foundation of China(Nos.12272397 and 11902334),the National Numerical Wind Tunnel Project,China。
文摘The matrix version of Symmetric Successive Over Relaxation(matrix-SSOR)scheme has been proved to be more efficient than the standard Lower-Upper Symmetric Gauss-Seidel(LUSGS),but less robust for high-speed flows.In order to ulteriorly improve the convergence rate as well as numerical stability of matrix-SSOR,two improvements regarding entropy fix and local time step have been proposed and validated.Firstly,an augmented entropy fix method is imposed on the inviscid Jacobian matrix and proved to be effective in two high-speed flows,in which the key parameter in entropy fix is discussed and found to be insensitive within appropriate range of values.Since the time step also has great effects on the numerical stability and convergence rate,a modified cell residual adapted local time step method with consideration of the residual history is developed,which is found to be effective for increasing the convergence rate when the matrix-SSOR is applied,but invalid when the LU-SGS is used.The proposed modified local time step method is also insensitive to the key parameter within appropriate range of values.The two modifications can be conveniently implanted into analogous matrix-type implicit schemes to improve the numerical performance.
基金support of the National Natural Science Foundation of China (No.52007002)the Science Center for Gas Turbine Project,China (No.P2022-A-II-007-001)the Fundamental Research Funds for the Central Universities,China (No.NS2023010).
文摘In order to further achieve the balance between the calculation accuracy and efficiency of the transient analysis of the aero-engine disc cavity system,an Optimized Time-adaptive Aerother-mal Coupling calculation(OTAC)method has been proposed.It combines one-dimensional tran-sient calculation of air system,Conventional Sequence Staggered(CSS)method,Time-adaptive Aerothermal Coupling calculation(TAC)method and differential evolution optimization algorithm to obtain an efficient and high-precision aerothermal coupling calculation method of air system.Considering both the heat conduction in the solid domain and the flow in the fluid domain as unsteady states in the OTAC,the interaction of fluid-solid information within a single coupling time step size was implemented based on the CSS method.Furthermore,the coupling time step size was automatically adjusted with the number of iterations by using the Proportional-Integral-Deri vative(PID)controller.Results show that when compared with the traditional loosely coupling method with a fixed time step size,the computational accuracy and efficiency of the OTAC method are improved by 8.9%and 30%,respectively.Compared with the tight coupling calculation,the OTAC method can achieve a speedup of 1 to 2 orders of magnitude,while the calculation error is maintained within 6.1%.
基金supported by the National Natural Science Foundation of China(Grant No.51079157)
文摘A growing interest has been devoted to the contra-rotating propellers （CRPs） due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibration. Compared with the single-screw system, it is more difficult for the open water performance prediction because forward and aft propellers interact with each other and generate a more complicated flow field around the CRPs system. The current work focuses on the open water performance prediction of contra-rotating propellers by RANS and sliding mesh method considering the effect of computational time step size and turbulence model. The validation study has been performed on two sets of contra-rotating propellers developed by David W Taylor Naval Ship R ＆ D center. Compared with the experimental data, it shows that RANS with sliding mesh method and SST k-ω turbulence model has a good precision in the open water performance prediction of contra-rotating propellers, and small time step size can improve the level of accuracy for CRPs with the same blade number of forward and aft propellers, while a relatively large time step size is a better choice for CRPs with different blade numbers.
基金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.
文摘There are two models in use today to analyze structural responses when subjected to earthquake ground motions, the Displacement Input Model （DIM） and the Acceleration Input Model （AIM）. The time steps used in direct integration methods for these models are analyzed to examine the suitability of DIM. Numerical results are presented and show that the time-step for DIM is about the same as for AIM, and achieves the same accuracy. This is contrary to previous research that reported that there are several sources of numerical errors associated with the direct application of earthquake displacement loading, and a very small time step is required to define the displacement record and to integrate the dynamic equilibrium equation. It is shown in this paper that DIM is as accurate and suitable as, if not more than, AIM for analyzing the response of a structure to uniformly distributed and spatially varying ground motions.
基金Hundred Talent Program of Chinese Academy of Sciences under Grant No. 0300YQ000101. Partly supported by the National Natural Sci
文摘Errors due to split time stepping are discussed for an explicit free–surface ocean model. In commonly used split time stepping, the way of time integration for the barotropic momentum equation is not compatible with that of the baroclinic one. The baroclinic equation has three–time–level structure because of leapfrog scheme. The barotropic one, however, has two–time–level structure when represented in terms of the baroclinic time level, on which the baroclinic one is integrated. This incompatibility results in the splitting errors as shown in this paper. The proper split time stepping is therefore proposed in such a way that the compatibility is kept between the barotropic and baroclinic equations. Its splitting errors are shown extremely small, so that it is particularly relevant to long–term integration for climate studies. It is applied to a free–surface model for the North Pacific Ocean.
基金The project supported by the National Key Basic Research and Development Foundation of the Ministry of Science and Technology of China （G2000048702, 2003CB716707）the National Science Fund for Distinguished Young Scholars （10025208）+1 种基金 the National Natural Science Foundation of China （Key Program） （10532040） the Research Fund for 0versea Chinese （10228028）.
文摘This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.
基金The project supported by the National Natural Science Foundation of China(the 973 Project 2004CB619304).
文摘A multiple time step algorithm, called reversible reference system propagator algorithm, is introduced for the long time molecular dynamics simulation. In contrast to the conventional algorithms, the multiple time method has better convergence, stability and efficiency. The method is validated by simulating free relaxation and the hypervelocity impact of nano-clusters. The time efficiency of the multiple time step method enables us to investigate the long time interaction between lattice dislocations and low-angle grain boundaries.
文摘The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.
文摘With the cell vertex finite volume discretization in space and second order backward implicit discretization in time, 2D unsteady Navier Stokes equations are solved by a dual time stepping method to simulate compressible viscous flow around rigid airfoils in arbitrary unsteady motion. The selection of physical time step is not restricted by stability condition any more, and most of the successful acceleration techniques used in steady calculations can be implemented to increase the computation efficiency.
基金The project is supported by the National Natural Science Foundation of China. (Grant No. 50605024).
文摘The heat transfer during the casting solidification process includes the heat radiation of the high temperature casting and the mold,the heat convection between the casting and the mold,and the heat conduction inside the casting and from the casting to the mold. In this paper,a formula of time step in simulation of solidification is derived,considering the heat radiation,convection and conduction based on the conservation of energy. The different heat transfer conditions between the conventional sand casting and the permanent mold casting are taken into account in this formula. The characteristics of heat transfer in the interior and surface of the casting are also considered. The numerical experiments show that this formula can avoid computational dispersion,and improve the computational efficiency by about 20% in the simulation of solidification process.
基金the National Natural Science Foundation of China(No.11672101,No.11372099)the 12th Five-Year Supporting Plan Issue(No.2015 BAB07B10)+1 种基金Jiangsu Province Natural Science Fund Project(No.BK 20151493)the Postgraduate Research and Innovation Projects in Jiangsu Province(No.2014B 31614)for the financial support.
文摘This paper focuses on the study of the stability of explicit time integration algorithm for dynamic problem by the Extended Finite Element Method(XFEM).A new enrichment scheme of crack tip is proposed within the framework of XFEM.Then the governing equations are derived and evolved into the discretized form.For dynamic problem,the lumped mass and the explicit time algorithm are applied.With different grid densities and different forms of Newmark scheme,the Dynamic Stress Intensity Factor(DSIF)is computed by using interaction integral approach to reflect the dynamic response.The effectiveness of the proposed scheme is demonstrated through the numerical examples,and the critical time stepping in different situations are listed and analyzed to illustrate the factors that affect the numerical stability.
文摘A self-adaptive precise algorithm in the time domain was employed to solve 2-D nonlinear coupled heat and moisture transfer problems. By expanding variables at a discretized time interval, the variations of variables can be described more precisely,and a nonlinear coupled initial and boundary value problem was converted into a series of recurrent linear boundary value problems which are solved by FE technique. In the computation, no additional assumption and the nonlinear iteration are required, and a criterion for self-adaptive computation is proposed to maintain sufficient computing accuracy for the change sizes of time steps. In the numerical comparison, the variations of material properties with temperature, moisture content, and both temperature and moisture content are taken into account, respectively. Satisfactory results have been obtained, indicating that the proposed approach is capable of dealing with complex nonlinear problems.
基金Supported in part by the National Natural Science of China, NSF Grant No. DMS-8657319.
文摘In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entropy solution in the limit if Courant number is less than or equal to 1.
基金The Project Supported by National Natural Science Foundation of China.
文摘A natural generalization of random choice finite difference scheme of Harten and Lax for Courant number larger than 1 is obtained. We handle interactions between neighboring Riemann solvers by linear superposition of their conserved quantities. We show consistency of the scheme for arbitrarily large Courant numbers. For scalar problems the scheme is total variation diminishing.A brief discussion is given for entropy condition.
文摘A positive leap second will be introduced in the UTC time scale UTC(JATC)、UTC (CSAO) and UTC time signals of BPL、BPM transmittings at the end of June 1997. The sequence of dates of the UTC second markers will be:
基金Supported by the Project Lagrangian Zooplankton Computation and Experiment (CNRS Programme EC2CO)the National Natural Science Fundation of China (Nos. 40821004,40706059)
文摘This paper presents two comparisons or tests for a Lagrangian model of zooplankton dispersion:numerical schemes and time steps.Firstly,we compared three numerical schemes using idealized circulations.Results show that the precisions of the advanced Adams-Bashfold-Moulton(ABM) method and the Runge-Kutta(RK) method were in the same order and both were much higher than that of the Euler method.Furthermore,the advanced ABM method is more efficient than the RK method in computational memory requirements and time consumption.We therefore chose the advanced ABM method as the Lagrangian particle-tracking algorithm.Secondly,we performed a sensitivity test for time steps,using outputs of the hydrodynamic model,Symphonie.Results show that the time step choices depend on the fluid response time that is related to the spatial resolution of velocity fields.The method introduced by Oliveira et al.in 2002 is suitable for choosing time steps of Lagrangian particle-tracking models,at least when only considering advection.
基金China Scholarship Council(CSC)The University of Queensland for a Ph D fellowship。
文摘The long-term stability of the roof is particularly important in designing underground rock structures.To estimate the durability of roof strata in underground excavation,a computation scheme of subcritical crack growth is proposed in this study.By adopting the proposed method,the potential collapse location of strata is derivable in accordance with a static model,the durability of roof strata can be estimated,a dynamic time step control strategy is achieved to balance the accuracy and speed of computing,and the initial crack size of rock can be estimated.In addition to the above,a mechanical model of underground excavation with non-uniformly distributed loads and partially yielded foundation is presented as the prototypical case.A set of case studies is carried out that showcase a power correlation between applied stress and roof durability.The allowable applied tensile stress for a 100-year life cycle is about 76%of the tensile strength.By using the proposed subcritical crack growth computation scheme,the roof stability in an underground excavation can be identified not only from the spatial view but also from the temporal perspective.
