198 research outputs found

    Toward a Simple, Accurate Lagrangian Hydrocode.

    Full text link
    Lagrangian hydrocodes play an important role in the computation of transient, compressible, multi-material flows. This research was aimed at developing a simply constructed cell-centered Lagrangian method for the Euler equations that respects multidimensional physics while achieving second-order accuracy. Algorithms that can account for the multidimensional physics associated with acoustic wave propagation and vorticity transport are needed in order to increase accuracy and prevent mesh imprinting. Many of the building blocks of traditional finite volume schemes, such as Riemann solvers and spatial gradient limiters, have their foundations in one-dimensional ideas and so were not used here. Instead, multidimensional point estimates of the fluxes were computed with a Lax-Wendroff type procedure and then nonlinearly modified using a temporal flux limiting mechanism. The linear acoustic equations were used as a simplified test environment for the Lagrangian Euler system. Here Lax-Wendroff methods that exactly preserve vorticity were investigated and found to resist mesh imprinting. However, the dispersion properties of the schemes were poor and so third-order accurate vorticity preserving methods were developed to remedy the problem. The third-order methods guided the construction of a temporal limiting mechanism, which was then used in a vorticity preserving flux-corrected transport scheme. While the acoustic work was interesting in its own right, it also proved to be a useful stepping stone to Lagrangian hydrodynamics. The acoustics algorithms were extended to produce the Simple Lagrangian Method (SLaM). Standard test problems have shown that a first-order accurate version of the method is able to resist mesh imprinting and spurious vorticity despite its minimalistic structure. SLaM is capable of second-order accuracy with a simple parameter change and some preliminary work was done to extend the temporal flux limiting ideas from acoustics to the Lagrangian case. The limited SLaM method converges at second-order for smooth data and is able to capture shocks without producing large unphysical oscillations.PhDAerospace EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/113577/1/tblung_1.pd

    Stability of a Kirchhoff–Roe scheme for two-dimensional linearized Euler systems

    Get PDF
    International audienceBy applying Helmholtz decomposition, the unknowns of a linearized Euler system can be recast as solutions of uncoupled linear wave equations. Accordingly, the Kirchhoff expression of the exact solutions is recast as a time-marching, Lax–Wendroff type, numerical scheme for which consistency with one-dimensional upwinding is checked. This discretization, involving spherical means, is set up on a 2D uniform Cartesian grid, so that the resulting numerical fluxes can be shown to be conservative. Moreover, semi-discrete stability in the Hs norms and vorticity dissipation are established, along with practical second-order accuracy. Finally, some relations with former “shape functions” and “symmetric potential schemes” are highlighted

    Assessing Implicit Large Eddy Simulation for Two-Dimensional Flow

    Get PDF
    Implicit large eddy simulation (ILES) has been shown, in the literature, to have some success for three-dimensional flow (e.g. see [Grinstein, F. F., Margolin, L. G. and Rider, W. Implicit Large Eddy Simulation. {\it Cambridge}. 2007]), but it has not previously been examined for two-dimensional flow. This thesis investigates whether ILES can be applied successfully to two-dimensional flow. Modified equation analysis is used to demonstrate the similarities between the truncation errors of certain numerical schemes and the subgrid terms of the barotropic vorticity equation (BVE). This presents a theoretical motivation for the numerical testing. Burgers equation is first used as a model problem to develop the ideas and methodology. Numerical schemes that are known to model Burgers equation well (shock capturing schemes) are shown to be implicitly capturing the subgrid terms of the one-dimensional inviscid Burgers equation through their truncation errors. Numerical tests are performed on three equation sets (BVE, Euler equations and the quasi-geostrophic potential vorticity equation) to assess the application of ILES to two-dimensional flow. The results for each of these equation sets show that the schemes considered for ILES are able to capture some of the subgrid terms through their truncation errors. In terms of accuracy, the ILES schemes are comparable (or outperform) schemes with simple explicit subgrid models when comparing vorticity solutions with a high resolution reference vorticity solution. The results suggest that conservation of vorticity is important to the successful application of ILES to two-dimensional flow, whereas conservation of momentum is not. The schemes considered for ILES are able to successfully model the downscale enstrophy transfer, but none of the schemes considered for ILES (or the schemes with simple subgrid models) can model the correct upscale energy transfer from the subgrid to the resolved scales. Energy backscatter models are considered and are used with the ILES schemes. It is shown that it is possible to create an energy conserving and enstrophy dissipating scheme, composed of an ILES scheme and a backscatter model, that improves the accuracy of the vorticity solution (when compared with the corresponding ILES scheme without backscatter)

    ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement

    Full text link
    We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space--time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.Comment: With updated bibliography informatio

    An Entropy Stable Nodal Discontinuous Galerkin Method for the Two Dimensional Shallow Water Equations on Unstructured Curvilinear Meshes with Discontinuous Bathymetry

    Full text link
    We design an arbitrary high-order accurate nodal discontinuous Galerkin spectral element approximation for the nonlinear two dimensional shallow water equations with non-constant, possibly discontinuous, bathymetry on unstructured, possibly curved, quadrilateral meshes. The scheme is derived from an equivalent flux differencing formulation of the split form of the equations. We prove that this discretisation exactly preserves the local mass and momentum. Furthermore, combined with a special numerical interface flux function, the method exactly preserves the mathematical entropy, which is the total energy for the shallow water equations. By adding a specific form of interface dissipation to the baseline entropy conserving scheme we create a provably entropy stable scheme. That is, the numerical scheme discretely satisfies the second law of thermodynamics. Finally, with a particular discretisation of the bathymetry source term we prove that the numerical approximation is well-balanced. We provide numerical examples that verify the theoretical findings and furthermore provide an application of the scheme for a partial break of a curved dam test problem

    On the Acoustic Component of Active Flux Schemes for Nonlinear Hyperbolic Conservation Laws

    Full text link
    Current numerical methods used in production-level CFD codes are found to be lacking in many respects; they are only second-order accurate, rely on inherently one-dimensional solvers, and are ill-equipped to handle more complex fluid flow problems such as turbulence, aeroacoustics and vortical flows just to name a few. Recently, a new class of third-order methods known simply as Active Flux (AF) has been introduced to address some of these issues. The AF method is best understood as a finite-volume method with additional degrees of freedom (DOF) at the interface to independently evolve interface fluxes. It is a fully discrete, maximally stable method that uses continuous data representation, and because the interface fluxes are computed independently from the cell-average values, true multidimensional solvers can be used. This dissertation focuses on the development of the AF method aimed at solving conservation laws describing acoustic processes. The method is demonstrated for linear and nonlinear acoustic equations in two-dimensions as well as for the full Euler equations where we employ operator splitting between the advective and acoustic processes. Given its continuous representation, the AF method economically achieves third-order accuracy using only three DOF in two dimensions, which is comparable to the discontinuous Galerkin method using linear reconstruction (DG1). A direct comparison between the two methods for acoustic problems finds that the AF method is capable of matching the accuracy of DG1 with a mesh spacing about three times greater and uses time steps about 2.5 times longer. The AF solutions also display superior circular symmetry with significantly less scatter than DG1, which we attribute to the method being able to employ truly multidimensional solvers. In addition, we find that on the same grid and to achieve the same level of error, the computation time for the AF method is more than one magnitude less than DG1 and approximately 3 to 5 times less than DG with quadratic reconstruction (DG2). Finally, various boundary conditions are introduced and developed for the AF scheme including far-field and curved wall boundaries.PHDAerospace EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/140800/1/dufan_1.pd

    On adomian based numerical schemes for euler and navier-stokes equations, and application to aeroacoustic propagation

    Get PDF
    140 p.En esta tesis se ha desarrollado un nuevo método de integración en tiempo de tipo derivadas sucesivas (multiderivative), llamado ABS y basado en el algoritmo de Adomian. Su motivación radica en la reducción del coste de simulación para problemas en aeroacústica, muy costosos por su naturaleza transitoria y requisitos de alta precisión. El método ha sido satisfactoriamente empleado en ambas partes de un sistema híbrido, donde se distinguen la parte aerodinámica y la acústica.En la parte aerodinámica las ecuaciones de Navier-Stokes incompresibles son resueltas con unmétodo de proyección clásico. Sin embargo, la fase de predicción de velocidad ha sido modificadapara incluir el método ABS en combinación con dos métodos: una discretización espacial MAC devolúmenes finitos, y también con un método de alto orden basado en ADER. El método se ha validado respecto a los problemas (en 2D) de vórtices de Taylor-Green, y el desarrollo de vórticesde Karman en un cilindro cuadrado. La parte acústica resuelve la propagación de ondas descritaspor las ecuaciones linearizadas de Euler, empleando una discretización de Galerkin discontinua. El método se ha validado respecto a la ecuación de Burgers.El método ABS es sencillo de programar con una formulación recursiva. Los resultados demuestran que su sencillez junto con sus altas capacidades de adaptación lo convierten en un método fácilmente extensible a órdenes altos, a la vez que reduce el coste comparado con otros métodos clásicos
    corecore