Search results

This paper aims to provide a simple and accurate higher order predictor‐corrector integration which can be used in dynamic analysis and to compare it with the previous works.

The predictor‐corrector integration is defined by combining the higher order explicit and implicit integrations in which displacement and velocity are assumed to be functions of accelerations of several previous time steps. By studying the accuracy and stability conditions, the weighted factors and acceptable time step are determined.

Simplicity and vector operations plus accuracy and stability are the main specifications of the new predictor‐corrector method. This procedure can be used in linear and nonlinear dynamic analysis.

In the proposed integration, time step is assumed to be constant.

The numerical integration is the heart of a dynamic analysis. The result's accuracy is strongly influenced by the accuracy and stability of the numerical integration.

This paper presents simple and accurate predictor‐corrector integration based on accelerations of several previous time steps. This may be used as a routine in any dynamic analysis software to enhance accuracy and reduce computational time.

In the present work, a novel dual time stepping approach is applied to a quasi-implicit (QI) fractional step method and its performance is assessed against the classical versions of the QI procedure for the solution of incompressible Navier-Stokes equations. The paper aims to discuss these issues.

In the proposed method, a local time stepping algorithm is utilised to accelerate the solution to steady state, while the transient solution is recovered through the use of a dual time step. It is demonstrated that, unlike the classical fractional step method, the temporal convergence rate of the proposed method depends solely upon the choice of the time discretisation.

While additional stabilisation is the prerequisite for obtaining higher order accuracy in the standard QI methods, the proposed dual time stepping approach completely eliminates this requirement. In addition, the dual time stepping approach proposed achieves the correct formal accuracy in time for both velocity and pressure. It is also demonstrated that a time accuracy beyond second order for both pressure and velocity is possible. In summary, the proposed dual time approach to QI methods simplifies the algorithm, accelerates solution and achieves a higher order time accuracy.

The dual time stepping removed first order pressure error.

Gives introductory remarks about chapter 1 of this group of 31 papers, from ISEF 1999 Proceedings, in the methodologies for field analysis, in the electromagnetic community. Observes that computer package implementation theory contributes to clarification. Discusses the areas covered by some of the papers ‐ such as artificial intelligence using fuzzy logic. Includes applications such as permanent magnets and looks at eddy current problems. States the finite element method is currently the most popular method used for field computation. Closes by pointing out the amalgam of topics.

The purpose of this paper is to solve unsteady compressible Navier–Stokes equations without the commonly used dual-time loop. The authors would like to use an adaptive time-stepping (ATS)-based local error control instead of CFL-based time-stepping technique. Also, an all-speed flow algorithm is implemented with simple low dissipation AUSM convective scheme, which can be computed without preconditioning which in general destroys the time accuracy.

In transient flow computations, the time-step is generally determined from the CFL condition. In this paper, the authors demonstrate the usefulness of ATS based on local time-stepping previously used extensively in ordinary differential equations (ODE) integration. This method is implemented in an implicit framework to ensure the numerical domain of dependence always contains the physical domain of dependence.

In this paper, the authors limit their focus to capture the unsteady physics for three cases: Sod’s shock-tube problem, Stokes’ second problem and a circular cylinder. The use of ATS with local truncation error control enables the solver to use the maximum allowable time-step, for the prescribed tolerance of error. The algorithm is also capable of converging very rapidly to the steady state (if there is any) after the initial transient phase. The authors present here only the first-order time-stepping scheme. An algorithmic comparison is made between the proposed adaptive time-stepping method and the commonly used dual time-stepping approach that indicates the former will be more efficient.

The original method of ATS based on local error control is used extensively in ODE integration, whereas, this method is not so popular in the computational fluid dynamics (CFD) community. In this paper, the authors investigate its use in the unsteady CFD computations. The authors hope that it would provide CFD researchers with an algorithm based on an adaptive time-stepping approach for unsteady calculations.

The purpose of this paper is to investigate the effect of numerical boundary condition implementation on local error and convergence in L2-norm of a finite volume discretization of the transient heat conduction equation subject to several boundary conditions, and for cases with volumetric heat generation, using both fully implicit and Crank-Nicolson time discretizations. The goal is to determine which combination of numerical boundary condition implementation and time discretization produces the most accurate solutions with the least computational effort.

The paper studies several benchmark cases including constant temperature, convective heating, constant heat flux, time-varying heat flux, and volumetric heating, and compares the convergence rates and local to analytical or semi-analytical solutions.

The Crank-Nicolson method coupled with second-order expression for the boundary derivatives produces the most accurate solutions on the coarsest meshes with the least computation times. The Crank-Nicolson method allows up to 16X larger time step for similar accuracy, with nearly negligible additional computational effort compared with the implicit method.

The findings can be used by researchers writing similar codes for quantitative guidance concerning the effect of various numerical boundary condition approximations for a large class of boundary condition types for two common time discretization methods.

The paper provides a comprehensive study of accuracy and convergence of the finite volume discretization for a wide range of benchmark cases and common time discretization methods.

This paper aims to propose an efficient and accurate method to analyse the transient heat conduction in a periodic structure with moving heat sources.

The moving heat source is modelled as a localised Gaussian distribution in space. Based on the spatial distribution, the physical feature of transient heat conduction and the periodic property of structure, a special feature of temperature responses caused by the moving heat source is illustrated. Then, combined with the superposition principle of linear system, within a small time-step, computation of results corresponding to the whole structure excited by the Gaussian heat source is transformed into that of some small-scale structures. Lastly, the precise integration method (PIM) is used to solve the temperature responses of each small-scale structure efficiently and accurately.

Within a reasonable time-step, the heat source applied on a unit cell can only cause the temperature responses of a limited number of adjacent unit cells. According to the above feature and the periodic property of a structure, the contributions caused by the moving heat source for the most of time-steps are repeatable, and the temperature responses of the entire periodic structure can be obtained by some small-scale structures.

A novel numerical method is proposed for analysing moving heat source problems, and the numerical examples demonstrate that the proposed method is much more efficient than the traditional methods, even for larger-scale problems and multiple moving heat source problems.

The purpose of this paper is to investigate the influence of time steps, integration methods, and saturation modeling on the accuracy of the synchronous machine model. This model is compared with the PSCAD built‐in synchronous machine model in order to compare the accuracy of one of the most used synchronous machine models in a commercially available software versus a well‐documented and widely accepted state‐space synchronous machine model.

In the paper, a synchronous condenser with the saturation phenomenon is modeled using state‐space equations in the rotating dq‐reference frame and is implemented both in Matlab/Simulink and PSCAD. Integration methods of up to the fifth order are implemented for increased accuracy. The saturation modeling includes modeling of the saturation in both d‐ and q‐axis. A steady‐state and dynamic performance comparison towards the built‐in PSCAD synchronous machine model is performed. The saturation modeling does not include the saturation of the leakage fluxes.

When the forward Euler method is used, in order to obtain less than 5 percent error, the time step should not exceed 5 μs. The third‐order Runge‐Kutta method is the preferred choice and it provides desired accuracy when the time step is equal or smaller than 1,000 μs. The built‐in PSCAD model satisfies the error criteria for time steps smaller than 300 μs. A small discrepancy of 2 percent is found during the steady‐state test.

The paper presents the performance of the higher order integration methods in an EMTP‐type software environment where the trapezoidal integration method is most often used. It provides a good guide for building an owner‐defined model. A comparison of a dynamic performance between the publicly documented state‐space and a synchronous machine models commonly used for power system transient studies is presented.

This paper aims to present a modified precise integration time domain (PITD) method for the numerical solution of 2D scalar wave equation.

The split step (SS) scheme is applied to factorize the conventional PITD calculation into two sub-steps procedures and then field components can be updated along one spatial direction only in each sub-step. The perfectly matched layer (PML) absorber is extended to this method for modeling open region problems by using the stretched coordinate approach.

It is shown that this method requires less computation time and storage space in comparison with the conventional PITD method, yet maintains the numerical stability despite using large time steps.

The WE-PITD method requires the divergence free region, which may be a limit on its usage. Hence, there is a challenge of using this technique in the 3D problems.

Based on the SS scheme, the PITD method is used to solve the scale wave equation rather than Maxwell's equations, leading to a significant reduction in the computation time and memory usage.

The purpose of this paper is to design a simple and accurate a-posteriori Lagrangian-based error estimator is developed for the class of backward differentiation formula (BDF) algorithms with variable time step size, and the adaptive time-stepping in BDF algorithms is demonstrated for efficient time-dependent simulations in fluid flow and heat transfer.

The Lagrange interpolation polynomial is used to predict the time derivative, and then the accurate primary result is obtained by the Gauss integral, which is applied to evaluate the local error. Not only the generalized formula of the proposed error estimator is presented but also the specific expression for the widely applied BDF1/2/3 is illustrated. Two essential executable MATLAB functions to implement the proposed error estimator are appended for practical applications. Then, the adaptive time-stepping is demonstrated based on the newly proposed error estimator for BDF algorithms.

The validation tests show that the newly proposed error estimator is accurate such that the effectivity index is always close to unity for both linear and nonlinear problems, and it avoids under/overestimation of the exact local error. The applications for fluid dynamics and coupled fluid flow and heat transfer problems depict the advantage of adaptive time-stepping based on the proposed error estimator for time-dependent simulations.

In contrast to existing error estimators for BDF algorithms, the present work is more accurate for the local error estimation, and it can be readily extended to practical applications in engineering with a few changes to existing codes, contributing to efficient time-dependent simulations in fluid flow and heat transfer.

The purpose of this paper is to describe a novel universal error estimator and the adaptive time-stepping process in the generalized single-step single-solve (GS4-1) computational framework, applied for the fluid dynamics with illustrations to incompressible Navier–Stokes equations.

The proposed error estimator is universal and versatile that it works for the entire subsets of the GS4-1 framework, encompassing the nondissipative Crank–Nicolson method, the most dissipative backward differential formula and anything in between. It is new and novel that the cumbersome design work of error estimation for specific time integration algorithms can be avoided. Regarding the numerical implementation, the local error estimation has a compact representation that it is determined by the time derivative variables at four successive time levels and only involves vector operations, which is simple for numerical implementation. Additionally, the adaptive time-stepping is further illustrated by the proposed error estimator and is used to solve the benchmark problems of lid-driven cavity and flow past a cylinder.

The proposed computational procedure is capable of eliminating the nonphysical oscillations in GS4-1(1,1)/Crank–Nicolson method; being CPU-efficient in both dissipative and nondissipative schemes with better solution accuracy; and detecting the complex physics and hence selecting a suitable time step according to the user-defined error threshold.

To the best of the authors’ knowledge, for the first time, this study applies the general purpose GS4-1 family of time integration algorithms for transient simulations of incompressible Navier–Stokes equations in fluid dynamics with constant and adaptive time steps via a novel and universal error estimator. The proposed computational framework is simple for numerical implementation and the time step selection based on the proposed error estimation is efficient, benefiting to the computational expense for transient simulations.