Computational uncertainty and optimal grid size and time step of the Lax-Friedrichs scheme for the 1D advection equation

本文对于应用Lax-Friedrichs格式数值求解一维平流方程,研究数值求解过程中产生的截断误差与舍入误差,以及两种误差逐层向高时间层传播的累积,得到新的数值解总误差上界的理论近似公式,以及最优格距和最优时间步长的理论公式.通过数值算例验证了所得公式的可靠性。然后,发现了两种不同机器精度下最优时间步长之比满足的一个仅与机器精度有关的普适关系.最后,理论验证了在网格比固定的情况下,此问题满足数值计算的不确定性原理,以及在机器有限精度下最优时间步长的必然存在。

Effect of Time Step Size and Turbulence Model on the Open Water Hydrodynamic Performance Prediction of Contra-Rotating Propellers

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.

Modified precise time step integration method of structural dynamic analysis

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.

Multiple time step molecular dynamics simulation for interaction between dislocations and grain boundaries

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.

Derivation and application of time step model in solidification process simulation

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 Stable Explicit Time Stepping Analysis with a New Enrichment Scheme by XFEM

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.

ENTROPY CONSISTENCY OF LARGE TIME STEP SCHEMES FOR ISENTROPIC EQUATIONS OF GAS DYNAMICS

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.

LARGE TIME STEP GENERALIZATION OF RANDOM CHOICE FINITE DIFFERENCE SCHEME FOR HYPERBOLIC CONSERVATION LAWS

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.

UTC TIME STEP on the 1st of July 1997

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:

Splitting Based Scheme for Three-dimensional MHD with Dual Time Stepping

A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avoid the occurrence of negative pressure in the reconstructed profiles and its updated value,a positivity preserving method is provided.Furthermore,the MHD equations are solved at each physical time step by advancing in pseudo time.The use of dual time stepping is beneficial in the computation since the use of dual time stepping allows the physical time step not to be limited by the corresponding values in the smallest cell and to be selected based on the numerical accuracy criterion.This newly established hybrid scheme combined with positivity preserving method and dual time technique has demonstrated the accurateness and robustness through numerical experiments of benchmark problems such as the 2D Orszag-Tang vortex problem and the3 D shock-cloud interaction problem.

Influence of Initial Value of Transient Time Step on Numerical Simulation of Blow Molding Balloon

The correlation between the initial time value of transient iterative parameters and the blowing pressure in the numerical simulation process of blowing balloon is investigated by POLYFLOW. The results show that: 1) As the blow molding pressure increases, the boundary value of the iterative time step decreases rapidly at first and then slowly. At the end of the first step of iterative calculation for each boundary value, the balloon parison is in the mold core cavity. 2) If the initial time value of transient iterative parameters is smaller than the boundary value of the iterative time step, the balloon parison is still in the mold core cavity at the end of the first iteration. However, if the iterative calculation continues, the calculation process may be interrupted when the time step is smaller than the initial time value of the transient iterative parameters, which makes the blow molding simulation of balloon unable to continue. 3) It is suggested that the initial time value of transient iterative parameters is one order of magnitude smaller than the boundary value of the iterative time step to complete smoothly the simulation of blow molding balloon.

Time Step Issue in Unit Hydrograph for Improving Runoff Prediction in Small Catchments

Unit hydrograph is a very practical tool in runoff prediction which has been used since decades ago and to date it remains useful. Unit hydrograph method is applied in Way Kuala Garuntang, an ungauged catchment in Lampung Province, Indonesia. To derive an observed unit hydrograph it requires rainfall and water level data with fine time scale which are obtained from automatic gauges. Observed unit hydrograph has an advantage that it is possible to derive it for various time steps including those with time step less than an hour. In order to get a more accurate unit hydrograph, it is necessary to derive a unit hydrograph with small time step for a small catchment such as those used in this study. The study area includes Way Kuala Garuntang and its tributaries, i.e. Way Simpur, Way Awi with areas are 60.52 km2, 3.691 km2, and 9.846 km2 respectively. The results of this study highlight the importance of time step selection on unit hydrograph, which are shown to have a significant impact on the resulting unit hydrograph’s variables such as peak discharge and time to peak.

Improvements for matrix-type implicit temporal scheme regarding entropy fix and local time step

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.

A Time Series Intrusion Detection Method Based on SSAE,TCN and Bi-LSTM

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.

On time-step in structural seismic response analysis under ground displacement/acceleration

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.

Theoretical development for DSMC local time stepping technique

The direct simulation Monte Carlo(DSMC) method is the most mature and wildly used approach for nonequilibrium gas flow simulation.The phenomenological nature of this method brings flexibility to the computation algorithms.In this study,the theoretical foundations to decouple the molecular motion and collision within a time step are discussed in detail,which can be treated as criterions for the DSMC algorithms.Based on the theoretical developments,an improved local time stepping scheme is proposed,which specifies the movement time attribute and the collision time attribute for each representative particle.A free flow about a sphere body is considered as an example,which is compared with the calculations using the published local time stepping technique.The results show that the improved local time scheme is valid and is promising in realizing flow structures with strong variations.

A Proper Time Integration with Split Stepping for the Explicit Free-Surface Modeling

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.

A DECOUPLING METHOD WITH DIFFERENT SUBDOMAIN TIME STEPS FOR THE NON-STATIONARY NAVIER-STOKES/DARCY MODEL

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.

Picard-Newton Iterative Method with Time Step Control for Multimaterial Non-Equilibrium Radiation Diffusion Problem

For a new nonlinear iterative method named as Picard-Newton(P-N)iterative method for the solution of the time-dependent reaction-diffusion systems,which arise in non-equilibrium radiation diffusion applications,two time step control methods are investigated and a study of temporal accuracy of a first order time integration is presented.The non-equilibrium radiation diffusion problems with flux limiter are considered,which appends pesky complexity and nonlinearity to the diffusion coef-ficient.Numerical results are presented to demonstrate that compared with Picard method,for a desired accuracy,significant increase in solution efficiency can be obtained by Picard-Newton method with the suitable time step size selection.

Optimization of the imaginary time step evolution for the Dirac equation

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.