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.

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

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

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.

RETHINKING TO FINITE DIFFERENCE TIME-STEP INTEGRATIONS

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.

IMEX Large Time Step Finite Volume Methods for Low Froude Number Shallow Water Flows

We present new large time step methods for the shallow water flows in the lowFroude number limit.In order to take into accountmultiscale phenomena that typically appear in geophysical flows nonlinear fluxes are split into a linear part governing the gravitational waves and the nonlinear advection.We propose to approximate fast linear waves implicitly in time and in space bymeans of a genuinely multidimensional evolution operator.On the other hand,we approximate nonlinear advection part explicitly in time and in space bymeans of themethod of characteristics or some standard numerical flux function.Time integration is realized by the implicit-explicit(IMEX)method.We apply the IMEX Euler scheme,two step Runge Kutta Cranck Nicolson scheme,as well as the semi-implicit BDF scheme and prove their asymptotic preserving property in the low Froude number limit.Numerical experiments demonstrate stability,accuracy and robustness of these new large time step finite volume schemes with respect to small Froude number.