Search results

A consistent implementation of the general computational framework of unified second-order time accurate integrators via the well-known GSSSS framework in conjunction with the traditional Finite Difference Method is presented to improve the numerical simulations of reactive two-phase flows.

In the present paper, the phase interaction evaluation in the present implementation of the reactive two-phase flows has been derived and implemented to preserve the consistency of the correct time level evaluation during the time integration process for solving the two phase flow dynamics with reactions.

Numerical examples, including the classical Sod shock tube problem and a reactive two-phase flow problem, are exploited to validate the proposed time integration framework and families of algorithms consistently to second order in time accuracy; this is in contrast to the traditional practices which only seem to obtain first-order time accuracy because of the inconsistent time level implementation with respect to the interaction of two phases. The comparisons with the traditional implementation and the advantages of the proposed implementation are given in terms of the improved numerical accuracy in time. The proposed approaches provide a correct numerical simulation implementation to the reactive two-phase flows and can obtain better numerical stability and computational features.

The new algorithmic framework and the consistent time level evaluation extended with the GS4 family encompasses a multitude of past and new schemes and offers a general purpose and unified implementation for fluid dynamics.

This paper aims to investigate post-buckling responses of shell-like structures using an implicit conservative-decaying time integration dynamic scheme.

In this work, the authors have proposed the use of a four-node quadrilateral flat shell finite element with drilling rotational degree of freedom within the framework of an updated Lagrangian formulation mutually with an implicit conservative-dissipative time integration dynamic scheme.

Several numerical simulations were considered to evaluate the accuracy, robustness, stability and the capacity of the considered time integration scheme to dissipate numerical noise in the presence of high frequencies. The obtained results illustrate a very satisfying performance of the implicit conservative-dissipative direct time integration scheme conjointly with the quadrilateral flat shell finite element with drilling rotation.

The authors have investigated the potential of the implicit dynamic scheme to deal with unstable branches after limit points in the non-linear post-buckling response of shell structures with no need for structural damping. The capability of the studied algorithm to study buckling and post-buckling behaviour of thin shell structures is illustrated through several numerical examples.

A transient magneto-quasistatic vector potential formulation involving nonlinear material is spatially discretized using the finite element method of first and second polynomial order. By applying a generalized Schur complement the resulting system of differential algebraic equations is reformulated into a system of ordinary differential equations (ODE). The ODE system is integrated in time by using explicit time integration schemes. The purpose of this paper is to investigate explicit time integration for eddy current problems with respect to the performance of the first-order explicit Euler scheme and the Runge-Kutta-Chebyshev (RKC) method of higher order.

The ODE system is integrated in time using the explicit Euler scheme, which is conditionally stable by a maximum time step size. To overcome this limit, an explicit multistage RKC time integration method of higher order is used to enlarge the maximum stable time step size. Both time integration methods are compared regarding the overall computational effort.

The numerical simulations show that a finer spatial discretization forces smaller time step sizes. In comparison to the explicit Euler time integration scheme, the multistage RKC method provides larger stable time step sizes to diminish the overall computation time.

The explicit time integration of the Schur complement vector potential formulation of eddy current problems is accelerated by a multistage RKC method.

To propose novel predictor‐corrector time‐integration algorithms for pseudo‐dynamic testing.

The novel predictor‐corrector time‐integration algorithms are based on both the implicit and the explicit version of the generalized‐α method. In the non‐linear unforced case second‐order accuracy, stability in energy, energy decay in the high‐frequency range as well as asymptotic annihilation are distinctive properties of the generalized‐α scheme; while in the non‐linear forced case they are the limited error near the resonance in terms of frequency location and intensity of the resonant peak. The implicit generalized‐α algorithm has been implemented in a predictor‐one corrector form giving rise to the implicit IPC‐ρ_∞ method, able to avoid iterative corrections which are expensive from an experimental standpoint and load oscillations of numerical origin. Moreover, the scheme embodies a secant stiffness formula able to approximate closely the actual stiffness of a structure. Also an explicit algorithm has been implemented, the EPC‐ρ_b method, endowed with user‐controlled dissipation properties. The resulting schemes have been tested experimentally both on a two‐ and on a six‐degrees‐of‐freedom system, exploiting substructuring techniques.

The analytical findings and the tests have indicated that the proposed numerical strategies enhance the performance of the pseudo‐dynamic test (PDT) method even in an environment characterized by considerable experimental errors. Moreover, the schemes have been tested numerically on strongly non‐linear multiple‐degrees‐of‐freedom systems reproduced with the Bouc‐Wen hysteretic model, showing that the proposed algorithms reap the benefits of the parent generalized‐α methods.

Further developments envisaged for this study are the application of the IPC‐ρ_∞ method and of EPC‐ρ_b scheme to partitioned procedures for high‐speed pseudo‐dynamic testing with substructuring.

The implicit IPC‐ρ_∞ and the explicit EPC‐ρ_b methods allow a user to have defined dissipation which reduces the effects of experimental error in the PDT without needing onerous iterations.

The paper proposes novel time‐integration algorithms for pseudo‐dynamic testing. Thanks to a predictor‐corrector form of the generalized‐α method, the proposed schemes maintain a high computational efficiency and accuracy.

The purpose of this paper is to analyse algorithms for fluid‐structure interaction (FSI) from a purely algorithmic point of view.

First of all a 1D model problem is selected, for which both the fluid and structural behavior are represented through a minimum number of parameters. Different coupling algorithm and time integration schemes are then applied to the simplified model problem and their properties are discussed depending on the values assumed by the parameters. Both exact and approximate time integration schemes are considered in the same framework so to allow an assessment of the different sources of error.

The properties of staggered coupling schemes are confirmed. An insight on the convergence behavior of iterative coupling schemes is provided. A technique to improve such convergence is then discussed.

All the results are proved for a given family of time integration schemes. The technique proposed can be applied to other families of time integration techniques, but some of the analytical results need to be reworked under this assumption.

The problems that are commonly encountered in FSI can be justified by simple arguments. It can also be shown that the limit at which trivial iterative schemes experience convergence difficulties is very close to that at which staggered schemes become unstable.

All the results shown are based on simple mathematics. The problems are presented so to be independent of the particular choice for the solution of the fluid flow.

The distinct element method as proposed by Cundall and Strack uses the computationally efficient, explicit, central difference time integration scheme. A limitation of this scheme is that it is only conditionally stable, so small time steps must be used. Some researchers have proposed using an implicit time integration scheme to avoid the stability issues arising from the explicit time integrator typically used in these simulations. However, these schemes are computationally expensive and can require a significant number of iterations to form the stiffness matrix that is compatible with the contact state at the end of each time step. In this paper, a new, simple approach for calculating the critical time increment in explicit discrete element simulations is proposed. Using this approach, it is shown that the critical time increment is a function of the current contact conditions. Considering both two‐ and three‐dimensional scenarios, the proposed refined estimates of the critical time step indicate that the earlier recommendations contained in the literature can be unconservative, in that they often overestimate the actual critical time step. A three‐dimensional simulation of a problem with a known analytical solution illustrates the potential for erroneous results to be obtained from discrete element simulations, if the time‐increment exceeds the critical time step for stable analysis.

Although the existence of a close relationship between the areas of digital signal processing and time integration methodology is known, a systematic application of the concepts and methods of the first area to the second is missing. Such an approach is followed in this paper, arising from the fact that any time integration formula can be viewed as a digital filter of the applied excitation force, approximating as close as possible to the behaviour of a ‘prototype analogue filter’, which is in fact the semi discrete equations of motion of the system. This approach provides a universal framework for handling and analysing all various aspects of time integration formulae, such as analysis in the frequency domain, algebraic operations, accuracy and stability, aliasing, spurious oscillations generation, introduction of digital filters within the time integration formula, initial conditions handling and overshooting. Additionally it is shown that digital signal processing methods, such as pre‐ or post‐processing, time delays, etc. can be in certain cases a quite effective complement of the time integration scheme.

In the original effective heat capacity method the latent heat effect is approximated by a large effective heat capacity over a small temperature range. This method is simple in concept and easy to implement. However, it is so sensitive to the choice of the phase change temperature interval and the chosen time integration scheme that non‐convergence always occurs when it is implemented with implicit time integration schemes. Presents the cause and cure of non‐convergence in effective heat capacity methods and discusses sample results.

The geometric conservation law (GCL) is an important concept for moving grid techniques because it directly regulates the treatments of the fluid flow and grid movement. With the grid movement at every time instant, the Jacobian, associated with the volume of each element in curvilinear co‐ordinates, needs to be updated in a conservative manner. In this study, alternative GCL schemes for evaluating the Jacobian have been investigated in the context of a pressure‐based Navier‐Stokes solver, utilizing moving grid and the first‐order implicit time stepping procedure as well as the PISO scheme. GCL‐based on first and second‐order, implicit as well as time‐averaged, time integration schemes were considered. Accuracy and conservative properties were tested on steady‐state, laminar flow inside a 2D channel and time dependent, turbulent flow around a 3D elastic wing; both treated with moving grid techniques. It seems that the formal order of accuracy is not a decisive indicator. Instead, the speed of grid movement and the interplay between the flow solver and the GCL treatments make a more noticeable impact.

A numerical integration scheme using the Monte Carlo method is discussed to evaluate the singular integral in boundary elements. A numerical demonstration is given for some potential problems. Results evaluated by the Monte Carlo method are compared with the analytical ones for accuracy and computation time. Examination shows the validity and capability of the approach.