Advection upstream splitting method

Last updated

The Advection Upstream Splitting Method (AUSM) is a numerical method used to solve the advection equation in computational fluid dynamics. It is particularly useful for simulating compressible flows with shocks and discontinuities.

Contents

The AUSM is developed as a numerical inviscid flux function for solving a general system of conservation equations. It is based on the upwind concept and was motivated to provide an alternative approach to other upwind methods, such as the Godunov method, flux difference splitting methods by Roe, and Solomon and Osher, flux vector splitting methods by Van Leer, and Steger and Warming.

The AUSM first recognizes that the inviscid flux consist of two physically distinct parts, i.e., convective and pressure fluxes. The former is associated with the flow (advection) speed, while the latter with the acoustic speed; or respectively classified as the linear and nonlinear fields. Currently, the convective and pressure fluxes are formulated using the eigenvalues of the flux Jacobian matrices. The method was originally proposed by Liou and Steffen [1] for the typical compressible aerodynamic flows, and later substantially improved in [2] [3] to yield a more accurate and robust version. To extend its capabilities, it has been further developed in [4] [5] [6] for all speed-regimes and multiphase flow. Its variants have also been proposed. [7] [8]

Features

The Advection Upstream Splitting Method has many features. The main features are:

Since the method does not specifically require eigenvectors, it is especially attractive for the system whose eigen-structure is not known explicitly, as the case of two-fluid equations for multiphase flow.

Applications

The AUSM has been employed to solve a wide range of problems, low-Mach to hypersonic aerodynamics, large eddy simulation and aero-acoustics, [9] [10] direct numerical simulation, [11] multiphase flow, [12] galactic relativistic flow [13] etc.

See also

Related Research Articles

<span class="mw-page-title-main">Computational fluid dynamics</span> Analysis and solving of problems that involve fluid flows

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

The finite volume method (FVM) is a method for representing and evaluating partial differential equations in the form of algebraic equations. In the finite volume method, volume integrals in a partial differential equation that contain a divergence term are converted to surface integrals, using the divergence theorem. These terms are then evaluated as fluxes at the surfaces of each finite volume. Because the flux entering a given volume is identical to that leaving the adjacent volume, these methods are conservative. Another advantage of the finite volume method is that it is easily formulated to allow for unstructured meshes. The method is used in many computational fluid dynamics packages. "Finite volume" refers to the small volume surrounding each node point on a mesh.

<span class="mw-page-title-main">Level-set method</span>

Level-set methods (LSM) are a conceptual framework for using level sets as a tool for numerical analysis of surfaces and shapes. The advantage of the level-set model is that one can perform numerical computations involving curves and surfaces on a fixed Cartesian grid without having to parameterize these objects. Also, the level-set method makes it very easy to follow shapes that change topology, for example, when a shape splits in two, develops holes, or the reverse of these operations. All these make the level-set method a great tool for modeling time-varying objects, like inflation of an airbag, or a drop of oil floating in water.

<span class="mw-page-title-main">Smoothed-particle hydrodynamics</span> Method of hydrodynamics simulation

Smoothed-particle hydrodynamics (SPH) is a computational method used for simulating the mechanics of continuum media, such as solid mechanics and fluid flows. It was developed by Gingold and Monaghan and Lucy in 1977, initially for astrophysical problems. It has been used in many fields of research, including astrophysics, ballistics, volcanology, and oceanography. It is a meshfree Lagrangian method, and the resolution of the method can easily be adjusted with respect to variables such as density.

In numerical methods, total variation diminishing (TVD) is a property of certain discretization schemes used to solve hyperbolic partial differential equations. The most notable application of this method is in computational fluid dynamics. The concept of TVD was introduced by Ami Harten.

Flux limiters are used in high resolution schemes – numerical schemes used to solve problems in science and engineering, particularly fluid dynamics, described by partial differential equations (PDEs). They are used in high resolution schemes, such as the MUSCL scheme, to avoid the spurious oscillations (wiggles) that would otherwise occur with high order spatial discretization schemes due to shocks, discontinuities or sharp changes in the solution domain. Use of flux limiters, together with an appropriate high resolution scheme, make the solutions total variation diminishing (TVD).

<span class="mw-page-title-main">High-resolution scheme</span> Scheme used in the numerical solution of partial differential equations

High-resolution schemes are used in the numerical solution of partial differential equations where high accuracy is required in the presence of shocks or discontinuities. They have the following properties:

In numerical analysis and computational fluid dynamics, Godunov's scheme is a conservative numerical scheme, suggested by Sergei Godunov in 1959, for solving partial differential equations. One can think of this method as a conservative finite volume method which solves exact, or approximate Riemann problems at each inter-cell boundary. In its basic form, Godunov's method is first order accurate in both space and time, yet can be used as a base scheme for developing higher-order methods.

In the study of partial differential equations, the MUSCL scheme is a finite volume method that can provide highly accurate numerical solutions for a given system, even in cases where the solutions exhibit shocks, discontinuities, or large gradients. MUSCL stands for Monotonic Upstream-centered Scheme for Conservation Laws, and the term was introduced in a seminal paper by Bram van Leer. In this paper he constructed the first high-order, total variation diminishing (TVD) scheme where he obtained second order spatial accuracy.

<span class="mw-page-title-main">Bram van Leer</span> Dutch mathematician

Bram van Leer is Arthur B. Modine Emeritus Professor of aerospace engineering at the University of Michigan, in Ann Arbor. He specializes in Computational fluid dynamics (CFD), fluid dynamics, and numerical analysis. His most influential work lies in CFD, a field he helped modernize from 1970 onwards. An appraisal of his early work has been given by C. Hirsch (1979)

<span class="mw-page-title-main">Riemann solver</span> Numerical method used to solve a Riemann problem

A Riemann solver is a numerical method used to solve a Riemann problem. They are heavily used in computational fluid dynamics and computational magnetohydrodynamics.

Flux-corrected transport (FCT) is a conservative shock-capturing scheme for solving Euler equations and other hyperbolic equations which occur in gas dynamics, aerodynamics, and magnetohydrodynamics. It is especially useful for solving problems involving shock or contact discontinuities. An FCT algorithm consists of two stages, a transport stage and a flux-corrected anti-diffusion stage. The numerical errors introduced in the first stage are corrected in the anti-diffusion stage.

In computational fluid dynamics, shock-capturing methods are a class of techniques for computing inviscid flows with shock waves. The computation of flow containing shock waves is an extremely difficult task because such flows result in sharp, discontinuous changes in flow variables such as pressure, temperature, density, and velocity across the shock.

<span class="mw-page-title-main">Volume of fluid method</span> Free-surface modelling technique

In computational fluid dynamics, the volume of fluid (VOF) method is a free-surface modelling technique, i.e. a numerical technique for tracking and locating the free surface. It belongs to the class of Eulerian methods which are characterized by a mesh that is either stationary or is moving in a certain prescribed manner to accommodate the evolving shape of the interface. As such, VOF is an advection scheme—a numerical recipe that allows the programmer to track the shape and position of the interface, but it is not a standalone flow solving algorithm. The Navier–Stokes equations describing the motion of the flow have to be solved separately. The same applies for all other advection algorithms.

Gretar Tryggvason is Department Head of Mechanical Engineering and Charles A. Miller Jr. Distinguished Professor at Johns Hopkins University. He is known for developing the front tracking method to simulate multiphase flows and free surface flows. Tryggvason was the editor-in-chief of Journal of Computational Physics from 2002–2015.

The moving particle semi-implicit (MPS) method is a computational method for the simulation of incompressible free surface flows. It is a macroscopic, deterministic particle method developed by Koshizuka and Oka (1996).

Vorticity confinement (VC), a physics-based computational fluid dynamics model analogous to shock capturing methods, was invented by Dr. John Steinhoff, professor at the University of Tennessee Space Institute, in the late 1980s to solve vortex dominated flows. It was first formulated to capture concentrated vortices shed from the wings, and later became popular in a wide range of research areas. During the 1990s and 2000s, it became widely used in the field of engineering.

<span class="mw-page-title-main">Phillip Colella</span> American mathematician

Phillip Colella is an American applied mathematician and a member of the Applied Numerical Algorithms Group at the Lawrence Berkeley National Laboratory. He has also worked at Lawrence Livermore National Laboratory. He is known for his fundamental contributions in the development of mathematical methods and numerical tools used to solve partial differential equations, including high-resolution and adaptive mesh refinement schemes. Colella is a member of the US National Academy of Sciences.

In applied mathematics, the name finite pointset method is a general approach for the numerical solution of problems in continuum mechanics, such as the simulation of fluid flows. In this approach the medium is represented by a finite set of points, each endowed with the relevant local properties of the medium such as density, velocity, pressure, and temperature.

The finite point method (FPM) is a meshfree method for solving partial differential equations (PDEs) on scattered distributions of points. The FPM was proposed in the mid-nineties in, and with the purpose to facilitate the solution of problems involving complex geometries, free surfaces, moving boundaries and adaptive refinement. Since then, the FPM has evolved considerably, showing satisfactory accuracy and capabilities to deal with different fluid and solid mechanics problems.

References

  1. ↑ Liou, M.-S. and Steffen, C., “A New Flux Splitting Scheme,” J. Comput. Phys., Vol. 107, 23-39, 1993.
  2. ↑ Liou, M.-S., “A Sequel to AUSM: AUSM+” J. Comput. Phys., Vol. 129, 364-382, 1996.
  3. ↑ Wada, Y. and Liou, M.-S., “An Accurate and Robust Flux Splitting Scheme for Shock and Contact Discontinuities,” SIAM J. Scientific Computing, Vol. 18, 633-657, 1997.
  4. ↑ Liou, M.-S., “A Sequel to AUSM, Part II: AUSM+-up” J. Comput. Phys., Vol. 214, 137- 170, 2006.
  5. ↑ Edwards, J. R., Franklin, R., and Liou, M.-S., “Low-Diffusion Flux-Splitting Methods for Real Fluid Flows with Phase Transitions,” AIAA J., Vol. 38, 1624-1633, 2000.
  6. ↑ Chang, C.-H. and Liou, M.-S., “A New Approach to the Simulation of Compressible Multifluid Flows with AUSM+ Scheme,” AIAA Paper 2003-4107, 16th AIAA CFD Conference, Orlando, FL, June 23–26, 2003.
  7. ↑ Edwards, J. R. and Liou, M.-S., “Low-Diffusion Flux-Splitting Methods for Flows at All Speeds,” AIAA J., Vol. 36, 1610-1617, 1998.
  8. ↑ Kim, K. H., Kim, C., and Rho, O., “Methods for the Accurate Computations of Hypersonic Flows I. AUSMPW+ Scheme,” J. Comput. Phys., Vol. 174, 38-80, 2001.
  9. ↑ Mary, I. and Sagaut, P., “Large Eddy Simulation of Flow Around an Airfoil Near Stall,” AIAA J., Vol. 40, 1139-1145, 2002.
  10. ↑ Manoha, E., Redonnet, S., Terracol, M., and Guenanff, G., “Numerical Simulation of Aerodynamics Noise,” ECCOMAS 24–28 July 2004.
  11. ↑ Billet, G. and Louedin, O., “Adaptive Limiters for Improving the Accuracy of the MUSCL Approach for Unsteady Flows,” J. Comput. Phys., Vol. 170, 161-183, 2001.
  12. ↑ Center for Risk Studies and Safety Archived April 24, 2006, at the Wayback Machine , University of California (Santa Barbara)
  13. ↑ Wada, K. and Koda, J., “Instabilities of Spiral Shock – I. Onset of Wiggle Instability and its Mechanism,” Monthly Notices of the Royal Astronomical Society, Vol. 349, 270-280 (11), 2004.
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.