Higher order method for the numerical solution of the compressible Euler equations
This paper deals with the numerical study of a nonlinear, strongly anisotropic heat equation. The use of standard schemes in this situation leads to poor results, due to the high anisotropy. An Asymptotic-Preserving method is introduced in this paper, which is second-order accurate in both, temporal and spacial variables. The discretization in time is done using an L-stable Runge−Kutta scheme. The convergence of the method is shown to be independent of the anisotropy parameter , and this for fixed...
In this paper we present a methodology for constructing accurate and efficient hybrid central-upwind (HCU) type schemes for the numerical resolution of a two-fluid model commonly used by the nuclear and petroleum industry. Particularly, we propose a method which does not make use of any information about the eigenstructure of the jacobian matrix of the model. The two-fluid model possesses a highly nonlinear pressure law. From the mass conservation equations we develop an evolution equation which...
In this paper we present a methodology for constructing accurate and efficient hybrid central-upwind (HCU) type schemes for the numerical resolution of a two-fluid model commonly used by the nuclear and petroleum industry. Particularly, we propose a method which does not make use of any information about the eigenstructure of the Jacobian matrix of the model. The two-fluid model possesses a highly nonlinear pressure law. From the mass conservation equations we develop an evolution equation which...
We present a hybrid OpenMP/MPI parallelization of the finite element method that is suitable to make use of modern high performance computers. These are usually built from a large bulk of multi-core systems connected by a fast network. Our parallelization method is based firstly on domain decomposition to divide the large problem into small chunks. Each of them is then solved on a multi-core system using parallel assembling, solution and error estimation. To make domain decomposition for both, the...
In this work we describe an efficient model for the simulation of a two-phase flow made of a gas and a granular solid. The starting point is the two-velocity two-pressure model of Baer and Nunziato [Int. J. Multiph. Flow16 (1986) 861–889]. The model is supplemented by a relaxation source term in order to take into account the pressure equilibrium between the two phases and the granular stress in the solid phase. We show that the relaxation process can be made thermodynamically coherent with an...
In this work three branches of Immersed Boundary Methods (IBM) are described and validated for incompressible aerodynamics and fluid-structure interactions. These three approaches are: Cut Cell method, Vortex-Penalization method and Forcing method. The first two techniques are validated for external bluff-body flow around a circular obstacle. The last one is used to predict the deformations of an elastic membrane immersed in a fluid. The paper confirms the ability of this family of numerical schemes...
Initial problems for nonlinear hyperbolic functional differential systems are considered. Classical solutions are approximated by solutions of suitable quasilinear systems of difference functional equations. The numerical methods used are difference schemes which are implicit with respect to the time variable. Theorems on convergence of difference schemes and error estimates of approximate solutions are presented. The proof of the stability is based on a comparison technique with nonlinear estimates...
Classical solutions of quasilinear functional differential equations are approximated with solutions of implicit difference schemes. Proofs of convergence of the difference methods are based on a comparison technique. Nonlinear estimates of the Perron type with respect to the functional variable for given functions are used. Numerical examples are given.
Solutions of initial boundary value problems for parabolic functional differential equations are approximated by solutions of implicit difference schemes. The existence and uniqueness of approximate solutions is proved. The proof of the stability is based on a comparison technique with nonlinear estimates of the Perron type for given operators. It is shown that the new methods are considerably better than the explicit difference schemes. Numerical examples are presented.
We analyze a two-stage implicit-explicit Runge–Kutta scheme for time discretization of advection-diffusion equations. Space discretization uses continuous, piecewise affine finite elements with interelement gradient jump penalty; discontinuous Galerkin methods can be considered as well. The advective and stabilization operators are treated explicitly, whereas the diffusion operator is treated implicitly. Our analysis hinges on L2-energy estimates on discrete functions in physical space. Our main...
We analyze a two-stage implicit-explicit Runge–Kutta scheme for time discretization of advection-diffusion equations. Space discretization uses continuous, piecewise affine finite elements with interelement gradient jump penalty; discontinuous Galerkin methods can be considered as well. The advective and stabilization operators are treated explicitly, whereas the diffusion operator is treated implicitly. Our analysis hinges on L2-energy estimates on discrete functions in physical space. Our main...