Toward a reduction of mesh imprinting

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/108641/1/fld3941.pd

Numerical simulation of spray coalescence in an eulerian framework : direct quadrature method of moments and multi-fluid method

The scope of the present study is Eulerian modeling and simulation of polydisperse liquid sprays undergoing droplet coalescence and evaporation. The fundamental mathematical description is the Williams spray equation governing the joint number density function f(v, u; x, t) of droplet volume and velocity. Eulerian multi-fluid models have already been rigorously derived from this equation in Laurent et al. (2004). The first key feature of the paper is the application of direct quadrature method of moments (DQMOM) introduced by Marchisio and Fox (2005) to the Williams spray equation. Both the multi-fluid method and DQMOM yield systems of Eulerian conservation equations with complicated interaction terms representing coalescence. In order to validate and compare these approaches, the chosen configuration is a self-similar 2D axisymmetrical decelerating nozzle with sprays having various size distributions, ranging from smooth ones up to Dirac delta functions. The second key feature of the paper is a thorough comparison of the two approaches for various test-cases to a reference solution obtained through a classical stochastic Lagrangian solver. Both Eulerian models prove to describe adequately spray coalescence and yield a very interesting alternative to the Lagrangian solver

Implementation of the LANS-alpha turbulence model in a primitive equation ocean model

This paper presents the first numerical implementation and tests of the Lagrangian-averaged Navier-Stokes-alpha (LANS-alpha) turbulence model in a primitive equation ocean model. The ocean model in which we work is the Los Alamos Parallel Ocean Program (POP); we refer to POP and our implementation of LANS-alpha as POP-alpha. Two versions of POP-alpha are presented: the full POP-alpha algorithm is derived from the LANS-alpha primitive equations, but requires a nested iteration that makes it too slow for practical simulations; a reduced POP-alpha algorithm is proposed, which lacks the nested iteration and is two to three times faster than the full algorithm. The reduced algorithm does not follow from a formal derivation of the LANS-alpha model equations. Despite this, simulations of the reduced algorithm are nearly identical to the full algorithm, as judged by globally averaged temperature and kinetic energy, and snapshots of temperature and velocity fields. Both POP-alpha algorithms can run stably with longer timesteps than standard POP. Comparison of implementations of full and reduced POP-alpha algorithms are made within an idealized test problem that captures some aspects of the Antarctic Circumpolar Current, a problem in which baroclinic instability is prominent. Both POP-alpha algorithms produce statistics that resemble higher-resolution simulations of standard POP. A linear stability analysis shows that both the full and reduced POP-alpha algorithms benefit from the way the LANS-alpha equations take into account the effects of the small scales on the large. Both algorithms (1) are stable; (2) make the Rossby Radius effectively larger; and (3) slow down Rossby and gravity waves.Comment: Submitted to J. Computational Physics March 21, 200

Multi-Dimensional, Compressible Viscous Flow on a Moving Voronoi Mesh

Numerous formulations of finite volume schemes for the Euler and Navier-Stokes equations exist, but in the majority of cases they have been developed for structured and stationary meshes. In many applications, more flexible mesh geometries that can dynamically adjust to the problem at hand and move with the flow in a (quasi) Lagrangian fashion would, however, be highly desirable, as this can allow a significant reduction of advection errors and an accurate realization of curved and moving boundary conditions. Here we describe a novel formulation of viscous continuum hydrodynamics that solves the equations of motion on a Voronoi mesh created by a set of mesh-generating points. The points can move in an arbitrary manner, but the most natural motion is that given by the fluid velocity itself, such that the mesh dynamically adjusts to the flow. Owing to the mathematical properties of the Voronoi tessellation, pathological mesh-twisting effects are avoided. Our implementation considers the full Navier-Stokes equations and has been realized in the AREPO code both in 2D and 3D. We propose a new approach to compute accurate viscous fluxes for a dynamic Voronoi mesh, and use this to formulate a finite volume solver of the Navier-Stokes equations. Through a number of test problems, including circular Couette flow and flow past a cylindrical obstacle, we show that our new scheme combines good accuracy with geometric flexibility, and hence promises to be competitive with other highly refined Eulerian methods. This will in particular allow astrophysical applications of the AREPO code where physical viscosity is important, such as in the hot plasma in galaxy clusters, or for viscous accretion disk models.Comment: 26 pages, 21 figures. Submitted to MNRA

Numerical wave propagation on the hexagonal C-grid

Copyright © 2008 Elsevier. NOTICE: this is the author’s version of a work that was accepted for publication in Journal of Computational Physics. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Journal of Computational Physics, Vol. 227, Issue 11 (2008), DOI: 10.1016/j.jcp.2008.02.010Inertio-gravity mode and Rossby mode dispersion properties are examined for discretizations of the linearized rotating shallow water equations on a regular hexagonal C-grid in planar geometry. It is shown that spurious non-zero Rossby mode frequencies found by previous authors in the f-plane case can be avoided by an appropriate discretization of the Coriolis terms. Three generalizations of this discretization that conserve energy even for non-constant Coriolis parameter are presented. A quasigeostrophic ββ-plane analysis is carried out to investigate the Rossby mode dispersion properties of these three schemes. The Rossby mode dispersion relation is found to have two branches. The primary branch modes are good approximations, in terms of both structure and frequency, to corresponding modes of the continuous governing equations, and offer some improvements over a quadrilateral C-grid scheme. The secondary branch modes have vorticity structures approximating those of small-scale modes of the continuous governing equations, suggesting that the hexagonal C-grid might have an advantage in terms of resolving extra Rossby modes; however, the frequencies of the secondary branch Rossby modes are much smaller than those of the corresponding continuous modes, so this potential advantage is not fully realized

Efficient Algorithm on a Non-staggered Mesh for Simulating Rayleigh-Benard Convection in a Box

An efficient semi-implicit second-order-accurate finite-difference method is described for studying incompressible Rayleigh-Benard convection in a box, with sidewalls that are periodic, thermally insulated, or thermally conducting. Operator-splitting and a projection method reduce the algorithm at each time step to the solution of four Helmholtz equations and one Poisson equation, and these are are solved by fast direct methods. The method is numerically stable even though all field values are placed on a single non-staggered mesh commensurate with the boundaries. The efficiency and accuracy of the method are characterized for several representative convection problems.Comment: REVTeX, 30 pages, 5 figure

Development of high-order realizable finite-volume schemes for quadrature-based moment method

Kinetic equations containing terms for spatial transport, gravity, fluid drag and particle-particle collisions can be used to model dilute gas-particle flows. However, the enormity of independent variables makes direct numerical simulation of these equations almost impossible for practical problems. A viable alternative is to reformulate the problem in terms of moments of velocity distribution. Recently, a quadrature-based moment method was derived by Fox for approximating solutions to kinetic equation for arbitrary Knudsen number. Fox also described 1st- and 2nd-order finite-volume schemes for solving the equations. The success of the new method is based on a moment-inversion algorithm that is used to calculate non-negative weights and abscissas from moments. The moment-inversion algorithm does not work if the moments are non-realizable, meaning they do not correspond to a distribution function. Not all the finite-volume schemes lead to realizable moments. Desjardins et al. showed that realizability is guaranteed with the 1 st-order finite-volume scheme, but at the expense of excess numerical diffusion. In the present work, the nonrealizability of the standard 2 nd-order finite-volume scheme is demonstrated and a generalized idea for the development of high-order realizable finite-volume schemes for quadrature-based moment methods is presented. This marks a significant improvement in the accuracy of solutions using the quadrature-based moment method as the use of 1st-order scheme to guarantee realizability is no longer a limitation

Some considerations for high‐order ‘incremental remap’‐based transport schemes: edges, reconstructions, and area integration

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/96654/1/fld3703.pd

Ignition and combustion development for high speed direct injection diesel engines under low temperature cold start conditions

Diesel engine cold start is an important issue for current technology at low (below 0 °C) temperatures and for future applications. The aim of this work is to develop a description of how, when and where does fuel spray ignition occur in a glow-plug assisted engine under simulated low temperature cold start conditions. In-cylinder pressure analysis is combined with high speed visualization in an optical engine. A pilot plus main injection strategy is used. Visualization results show that pilot ignition occurs in the vicinity of the glow plug, and strongly influences main combustion initiation. Main combustion starts from the pilot flame, and propagates to the rest of the combustion chamber showing a strong visible reaction zone. After end of main injection, the rapid leaning of the mixture suppresses the strong radiation, and OH radiation is observed to progress to the rest of the combustion chamber. The combustion process shows a strong scattering, which has been quantified by combustion parameters. At higher rail pressures scattering increases, which eventually inhibits combustion initiation. However, if ignition occurs at higher rail pressures, cycle performance is better.Authors thank the Spanish Ministry of Innovation and Science for the financial support through the project OPTICOMB (reference code: TRA2007-67961-C03-C01). Authors also thank Daniel Lerida Sanchez de las Heras for his outstanding work in the facility set-up and adaptation and for his support during the tests.Pastor Soriano, JV.; García Oliver, JM.; Pastor Enguídanos, JM.; Ramírez Hernández, JG. (2011). Ignition and combustion development for high speed direct injection diesel engines under low temperature cold start conditions. Fuel. 90(4):1556-1566. https://doi.org/10.1016/j.fuel.2011.01.008S1556156690