Study of high-order energy-stable discretization techniques for compressible flows

This work presents the development of an explicit/implicit particle finite element method (PFEM) for the 2D modeling of metal cutting processes. The purpose is to study the efficiency of implicit and explicit time integration schemes in terms of precision, accuracy and computing time. The formulation for implicit and explicit time marching schemes is developed, and a detailed study on the explicit solution steps is presented. The PFEM remeshing procedures for insertion and removal of particles have been improved to model the multiple scales of time and/or space of the solution. The detection and treatment of the rigid tool contact are presented for both, implicit and explicit schemes. The performance of explicit/implicit integration is studied with a set of different two-dimensional orthogonal cutting tests of AISI 4340 steel at cutting speeds ranging from 1 m/s up to 30 m/s. It was shown that if the correct selection of the time integration scheme is made, the computing time can decrease up to 40 times. It allows us to affirm that the computing time of the PFEM simulations can be excessive due to the used time marching scheme independently of the meshing process. As a practical result, a set of recommendations to select the time integration schemes for a given cutting speed are given. This is intended to minimize one of the negative constraints pointed out by the industry when using metal cutting simulators.

This study investigates the full penetration simulation of piles from the ground surface, focusing on frictional contact modeling without mesh distortion. To overcome issues related to mesh distortion and improve solution convergence, the Arbitrary Lagrangian–Eulerian (ALE) adaptive mesh technique was implemented within both explicit and implicit time integration schemes. The numerical model was validated against field experiments conducted at Bothkennar, Scotland, using the Imperial College instrumented displacement pile (ICP) in soft clay, where the soil behavior was effectively represented using the modified Cam-Clay model and the Mohr–Coulomb model. The primary objectives of this study are to evaluate the ALE method performance in handling mesh distortion; analyze the effects of soil–pile interface friction, pile dimensions, and various dilation angles on pile resistance; and compare the effectiveness of explicit and implicit time integration schemes in terms of stability, computational efficiency, and solution accuracy. The ALE method effectively modeled pile penetration in Bothkennar clay, validating the numerical model against field experiments. Comparative analysis revealed the explicit time integration method as more robust and computationally efficient, particularly for complex soil–pile interactions with higher friction coefficients.

We present implicit time integration schemes suitable for modeling free surface Stokes flow dynamics with marker in cell (MIC) based spatial discretization. Our target is for example thermal convection surrounded by deformable surface boundaries to simulate the long term planetary formation process. The numerical system becomes stiff when the dynamical balancing time scale for the increasing/decreasing load by surface deformation is very short compared with the time scale associated with thermal convection. Any explicit time integration scheme will require very small time steps; otherwise, serious numerical oscillation (spurious solutions) will occur. The implicit time integration scheme possesses a wider stability region than the explicit method; therefore, it is suitable for stiff problems. To investigate an efficient solution method for the stiff Stokes flow system, we apply first (backward Euler (BE)) and second order (trapezoidal method (TR) and trapezoidal rule—backward difference formula (TR-BDF2)) accurate implicit methods for the MIC solution scheme. The introduction of implicit time integration schemes results in nonlinear systems of equations. We utilize a Jacobian free Newton Krylov (JFNK) based Newton framework to solve the resulting nonlinear equations. In this work we also investigate two efficient implicit solution strategies to reduce the computational cost when solving stiff nonlinear systems. The two methods differ in how the advective term in the material transport evolution equation is treated. We refer to the method that employs Lagrangian update as “fully implicit” (Imp), whilst the method that employs Eulerian update is referred to as “semi-implicit” (SImp). Using a finite difference (FD) method, we have performed a series of numerical experiments which clarify the accuracy of solutions and trade-off between the computational cost associated with the nonlinear solver and time step size. In comparison with the general explicit Euler method, the second order accurate Imp methods reduce total computational cost successfully through the utilization of a large time step without sacrificing accuracy and stability. Moreover, the proposed SImp method is effective in reducing the computational cost associated with evaluating the nonlinear residual while obtaining a solution similar to the Imp method.

This paper is concerned with the accurate, conservative, and stable imposition of boundary conditions and inter-element coupling for multi-dimensional summation-by-parts (SBP) finite-difference operators. More precisely, the focus is on diagonal-norm SBP operators that are not based on tensor products and are applicable to unstructured grids composed of arbitrary elements. We show how penalty terms—simultaneous approximation terms (SATs)—can be adapted to discretizations based on multi-dimensional SBP operators to enforce boundary and interface conditions. A general SAT framework is presented that leads to conservative and stable discretizations of the variable-coefficient advection equation. This framework includes the case where there are no nodes on the boundary of the SBP element at which to apply penalties directly. This is an important generalization, because elements analogous to Legendre–Gauss collocation, i.e. without boundary nodes, typically have higher accuracy for the same number of degrees of freedom. Symmetric and upwind examples of the general SAT framework are created using a decomposition of the symmetric part of an SBP operator; these particular SATs enable the pointwise imposition of boundary and inter-element conditions. We illustrate the proposed SATs using triangular-element SBP operators with and without nodes that lie on the boundary. The accuracy, conservation, and stability properties of the resulting SBP–SAT discretizations are verified using linear advection problems with spatially varying divergence-free velocity fields.

Improved algorithms for the solution of the 3-D time dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time integration scheme using a multistage Runge-Kutta procedure or an implicit time integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the ONERA M6 wing to demonstrate applications of the new Euler solvers. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

A continuum damage model for concrete is developed with a focus on fatigue under compressive stresses. This includes the possibility to model stress redistributions and capture size effects. In contrast to cycle based approaches, where damage is accumulated based on the number of full stress cycles, a strain based approach is developed that can capture cyclic degradation under variable loading cycles including different amplitudes and loading frequencies. The model is designed to represent failure under static loading as a particular case of fatigue failure after a single loading cycle. As a consequence, most of the material parameters can be deduced from statictests. Only a limit set of additional constitutive parameters is required to accurately describe the evolution under fatigue loading. Another advantage of the proposed model is the possibility to directly incorporate other multi-physics effects such as creep and shrinkage or thermal loading on the constitutive level. A multiscale approach in time is presented to enable structural computations of fatigue failure with a reduced computational effort. The damage rate within the short time scale corresponding to a single cycle is computed based on a Fourier based approach. This evolution equation is then solved on the long time scale using different implicit and explicit time integration schemes. Their performance and some limitations for specific loading regimes is discussed.

Model reduction for multi-scale transport problems using model-form preserving least-squares projections with variable transformation

In this paper high order time integration schemes are investigated for the integration of uid-structure interaction (FSI) model problem with deforming domains. The integration algorithm consists of a partitioned scheme in which a combination of implicit and explicit (IMEX) high order Runge-Kutta time integration schemes is used: the uid and structure systems are decoupled at their interface and both systems are integrated using the implicit system under a given boundary condition at the interface. The inuence of this boundary condition (the coupling term) is integrated using the explicit Runge-Kutta scheme. The uid dynamics are written in the arbitrary Lagrange Eulerian formulation to cope with the deforming uid mesh. An expression is proposed for the mesh face velocities of the deforming mesh such that the discrete geometric conservation law is satised. For the model problem numerical investigations are made into the eciency of these high order time integration schemes compared to lower order methods. The IMEX scheme retains the order of the implicit and explicit schemes and obtains the same accuracy as a second order staggered scheme for less computational work.

An adaptive h‐version finite element procedure is presented for a 2D elastodynamic transient analysis using direct integration. The procedure modifies the spatial discretization, i.e. the mesh, in time so that the spatial discretization error can be controlled in each time instant. We analyse an a priori error estimate for the total energy which shows that an a posteriori error estimate developed by Zienkiewicz and Zhu for linear elastostatics can reasonably be extended to dynamic transient analysis. For the time integration, both explicit and implicit time integration schemes can be coupled with the procedure. However, the study in this paper is confined to the use of the Newmark scheme. Numerical examples are used to study the performance of the error estimate and to illustrate the adaptive procedure. Obtained numerical results indicate that: (i) the extension of the Zienkiewicz–Zhu estimate to dynamic transient analysis is useful; (ii) the presented adaptive procedure can, in an efficient way, provide a finite element solution with a user‐controlled accuracy.

Stability of semidiscrete formulations for parabolic problems at small time steps

In this thesis, the numerical solution of time-dependent partial differential equations (PDE) is studied. In particular high-order finite difference methods on Summation-by-parts (SBP) form are analysed and applied to model problems as well as the PDEs governing aerodynamics. The SBP property together with an implementation of boundary conditions called SAT (Simultaneous Approximation Term), yields stability by energy estimates. The first derivative SBP operators were originally derived for Cartesian grids. Since aerodynamic computations are the ultimate goal, the scheme must also be stable on curvilinear grids. We prove that stability on curvilinear grids is only achieved for a subclass of the SBP operators. Furthermore, aerodynamics often requires addition of artificial dissipation and we derive an SBP version. With the SBP-SAT technique it is possible to split the computational domain into a multi-block structure which simplifies grid generation and more complex geometries can be resolved. To resolve extremely complex geometries an unstructured discretisation method must be used. Hence, we have studied a finite volume approximation of the Laplacian. It can be shown to be on SBP form and a new boundary treatment is derived. Based on the Laplacian scheme, we also derive an SBP artificial dissipation for finite volume schemes. We derive a new set of boundary conditions that leads to an energy estimate for the linearised three-dimensional Navier-Stokes equations. The new boundary conditions will be used to construct a stable SBP-SAT discretisation. To obtain an energy estimate for the discrete equation, it is necessary to discretise all the second derivatives by using the first derivative approximation twice. According to previous theory that would imply a degradation of formal accuracy but we present a proof that this is not the case.

A phase field-immersed boundary-lattice Boltzmann coupling method for fluid-structure interaction analysis

Solving dynamic problems for fluid saturated porous media at large deformation regime is an interesting but complex issue. An implicit time integration scheme is herein developed within the framework of the u–w (solid displacement–relative fluid displacement) formulation for the Biot’s equations. In particular, liquid water saturated porous media is considered and the linearization of the linear momentum equations taking into account all the inertia terms for both solid and fluid phases is for the first time presented. The spatial discretization is carried out through a meshfree method, in which the shape functions are based on the principle of local maximum entropy LME. The current methodology is firstly validated with the dynamic consolidation of a soil column and the plastic shear band formulation of a square domain loaded by a rigid footing. The feasibility of this new numerical approach for solving large deformation dynamic problems is finally demonstrated through the application to an embankment problem subjected to an earthquake.

This study describes the application of two main improvements in highly viscous fluid simulations using the smoothed particle hydrodynamics (SPH) method: an implicit time integration scheme to overcome the problem of impractically small time step restriction and the introduction of air ghost particles to fix problems regarding the free surface treatment. This study adopts the incompressible SPH as a basis for the implementation of these improvements, which guarantees a stable and accurate pressure distribution. We verified the proposed implicit time integration scheme with simulations of pipe flow and the free surface treatment with a simple hydrostatic problem. As a result, the free surface of the hydrostatic problem became very smooth and stable. In addition, we conducted a variety of dam-break simulations to validate this proposed SPH method, as well as to analyze the density and divergence error. Finally, we demonstrate the potential of this method with the highly viscous vertical jet flow over a horizontal plate test, which features a complex viscous coiling behavior.

The lid-driven cavity test case is widely used to validate incompressible Navier–Stokes flow solvers. However, the rigorous treatment of discontinuous corner boundary conditions remains a challenge for high-order methods. This is the main goal of the paper. We write the incompressible Navier–Stokes equations in skew-symmetric form and we impose the boundary conditions weakly which leads to boundedness without any special treatment in the corners. The continuous procedure is mimicked in the discrete setting using high-order finite difference methods in Summation-By-Parts (SBP) form complemented with weak boundary conditions using the Simultaneous Approximation Term (SAT) technique. Stability is then formally proven using the SBP-SAT framework. Numerical tests commence with a method-of-manufactured solution and the scheme is shown to be high order accurate. The lid-driven cavity test case is finally studied using a 4th order accurate scheme and the results are compared to benchmark solutions. Accurate solutions are achieved that are devoid of spurious oscillations near the top corners and the velocities remain bounded, demonstrating the unique versatility of the weak SAT boundary treatment.