1,853 research outputs found
Effects of discrete energy and helicity conservation in numerical simulations of helical turbulence
Helicity is the scalar product between velocity and vorticity and, just like
energy, its integral is an in-viscid invariant of the three-dimensional
incompressible Navier-Stokes equations. However, space-and time-discretization
methods typically corrupt this property, leading to violation of the inviscid
conservation principles. This work investigates the discrete helicity
conservation properties of spectral and finite-differencing methods, in
relation to the form employed for the convective term. Effects due to
Runge-Kutta time-advancement schemes are also taken into consideration in the
analysis. The theoretical results are proved against inviscid numerical
simulations, while a scale-dependent analysis of energy, helicity and their
non-linear transfers is performed to further characterize the discretization
errors of the different forms in forced helical turbulence simulations
Effects of discrete energy and helicity conservation in numerical simulations of helical turbulence
Helicity is the scalar product between velocity and vorticity and, just like
energy, its integral is an in-viscid invariant of the three-dimensional
incompressible Navier-Stokes equations. However, space-and time-discretization
methods typically corrupt this property, leading to violation of the inviscid
conservation principles. This work investigates the discrete helicity
conservation properties of spectral and finite-differencing methods, in
relation to the form employed for the convective term. Effects due to
Runge-Kutta time-advancement schemes are also taken into consideration in the
analysis. The theoretical results are proved against inviscid numerical
simulations, while a scale-dependent analysis of energy, helicity and their
non-linear transfers is performed to further characterize the discretization
errors of the different forms in forced helical turbulence simulations
Conservative finite-volume framework and pressure-based algorithm for flows of incompressible, ideal-gas and real-gas fluids at all speeds
A conservative finite-volume framework, based on a collocated variable
arrangement, for the simulation of flows at all speeds, applicable to
incompressible, ideal-gas and real-gas fluids is proposed in conjunction with a
fully-coupled pressure-based algorithm. The applied conservative discretisation
and implementation of the governing conservation laws as well as the definition
of the fluxes using a momentum-weighted interpolation are identical for
incompressible and compressible fluids, and are suitable for complex geometries
represented by unstructured meshes. Incompressible fluids are described by
predefined constant fluid properties, while the properties of compressible
fluids are described by the Noble-Abel-stiffened-gas model, with the
definitions of density and specific static enthalpy of both incompressible and
compressible fluids combined in a unified thermodynamic closure model. The
discretised governing conservation laws are solved in a single linear system of
equations for pressure, velocity and temperature. Together, the conservative
finite-volume discretisation, the unified thermodynamic closure model and the
pressure-based algorithm yield a conceptually simple, but versatile, numerical
framework. The proposed numerical framework is validated thoroughly using a
broad variety of test-cases, with Mach numbers ranging from 0 to 239, including
viscous flows of incompressible fluids as well as the propagation of acoustic
waves and transiently evolving supersonic flows with shock waves in ideal-gas
and real-gas fluids. These results demonstrate the accuracy, robustness and the
convergence, as well as the conservation of mass and energy, of the numerical
framework for flows of incompressible and compressible fluids at all speeds, on
structured and unstructured meshes
Effects of discrete energy and helicity conservation in numerical simulations of helical turbulence
Helicity is the scalar product between velocity and vorticity and, just like energy, its integral is an inviscid invariant of the three-dimensional incompressible Navier-Stokes equations. However, space- and time-discretization methods typically corrupt this prop- erty, leading to violation of the inviscid conservation principles. This work investigates the discrete helicity conservation properties of spectral and finite-differencing methods, in relation to the form employed for the convective term. Effects due to Runge-Kutta time- advancement schemes are also taken into consideration in the analysis. The theoretical results are proved against inviscid numerical simulations, while a scale-dependent analy- sis of energy, helicity and their non-linear transfers is performed to further characterize the discretization errors of the different forms in forced helical turbulence simulations.Postprint (published version
A symmetry-preserving second-order time-accurate PISO-based method
A new conservative symmetry-preserving second-order time-accurate PISO-based pressure-velocity coupling for solving the incompressible Navier-Stokes equations on unstructured collocated grids is presented in this paper. This new method for implicit time stepping is an extension of the conservative symmetry-preserving incremental-pressure projection method for explicit time stepping and unstructured collocated meshes of Trias et al. [35]. In order to assess and compare both methods, we have implemented them within one unified solver in the open source code OpenFOAM where we use a Butcher array to prescribe the Runge-Kutta method. Thus, by changing the entries of the Butcher array, explicit and diagonally implicit Runge-Kutta schemes can be combined into one solver. We assess the energy conservation properties of the implemented discretisation methods and the temporal consistency of the selected Runge-Kutta schemes using Taylor-Green vortex and lid-driven cavity flow test cases. Finally, we use a more complex turbulent channel flow test case in order to further assess the performance of the presented new conservative symmetry-preserving incremental-pressure PISO-based method. Although both implemented methods are based on a symmetry-preserving discretisation, we show they still produce a small amount of numerical dissipation when the total pressure is directly solved from a Poisson equation. When an incremental-pressure approach is used, where a pressure correction is solved from a Poisson equation, both methods are effectively fully-conservative. For high-fidelity simulations of incompressible turbulent flows, it is highly desirable to use fully-conservative methods. For such simulations, the presented numerical methods are therefore expected to have large added value, since they pave the way for the execution of truly energy-conservative high-fidelity simulations in complex geometries. Furthermore, both methods are implemented in OpenFOAM, which is widely used within the CFD community, so that a large part of this community can benefit from the developed and implemented numerical methods
Center for Modeling of Turbulence and Transition (CMOTT). Research briefs: 1990
Brief progress reports of the Center for Modeling of Turbulence and Transition (CMOTT) research staff from May 1990 to May 1991 are given. The objectives of the CMOTT are to develop, validate, and implement the models for turbulence and boundary layer transition in the practical engineering flows. The flows of interest are three dimensional, incompressible, and compressible flows with chemistry. The schemes being studied include the two-equation and algebraic Reynolds stress models, the full Reynolds stress (or second moment closure) models, the probability density function models, the Renormalization Group Theory (RNG) and Interaction Approximation (DIA), the Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS)
- …