Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPressure-tight and non-stiff volume penalization for compressible flows
99
0
0.0
(
0
)
Authors
Julius Reiss
Ask ChatGPT about the research
No Arabic abstract
Embedding geometries in structured grids allows a simple treatment of complex objects in fluid flows. Various methods are available. The commonly used Brinkman-volume-penalization models geometries as porous media, where in the limit of vanishing porosity a solid object is approximated. In the simplest form, the velocity equations are augmented by a term penalizing the fluid velocity, the body velocity. yielding good results in many applications. However, it induces numerical stiffness, especially if high pressure gradients need to be balanced. Here, we focus on the effect of the reduced effective volume (commonly called porosity) of the porous medium. An approach is derived, which allows to reduce the flux through objects to practically zero with little increase of numerical stiffness. Also, non-slip boundary conditions and adiabatic boundary conditions are easily constructed. The porosity terms allow to keep the skew symmetry of the underlying numerical scheme, by which the numerical stability is improved. Furthermore, a very good conservation of mass and energy in the non-penalized domain can be achieved. The scheme is tested for acoustic scenarios, near incompressible and strongly compressible flows.
rate research
Read More
A two-phase, low-Mach-number flow solver is proposed for variable-density liquid and gas with phase change. The interface is captured using a split Volume-of-Fluid method, which solves the advection of the reference phase, generalized for the case where the liquid velocity is not divergence-free and both phases exchange mass. A sharp interface is identified by using PLIC. Mass conservation is achieved in the limit of incompressible liquid, but not with the liquid compressibility and mass exchange. This is a relevant modeling choice for two-phase mixtures at near-critical and supercritical pressure conditions for the liquid but away from the mixture critical temperature. Under this thermodynamic environment, the dissolution of lighter gas species into the liquid phase is enhanced and vaporization or condensation can occur simultaneously at different interface locations. The numerical challenge of solving two-phase, supercritical-pressure flows is greater than simpler two-phase solvers because: a) local phase equilibrium is imposed at each interface cell to determine temperature, composition, or surface tension coefficient; b) a real-fluid thermodynamic model is used to obtain fluid properties; and c) necessary phase-wise values for certain variables are obtained via extrapolation techniques. To alleviate the increased numerical cost, the pressure Poisson equation (PPE) used to solve the low-Mach-number flow is split into a constant-coefficient implicit part and a variable-coefficient explicit part. Thus, a Fast Fourier Transform method can be used for the PPE. Various verification tests are performed to show the accuracy and viability of the present approach. The growth of surface instabilities in a binary system composed of liquid n-decane and gaseous oxygen at supercritical pressures for n-decane is analyzed. Other features of supercritical liquid injection are also shown.
We study numerically joint mixing of salt and colloids by a chaotic velocity field $mathbf{V}$, and how salt inhomogeneities accelerate or delay colloid mixing by inducing a velocity drift $mathbf{V}_{rm dp}$ between colloids and fluid particles as proposed in recent experiments cite{Deseigne2013}. We demonstrate that because the drift velocity is no longer divergence free, small variations to the total velocity field drastically affect the evolution of colloid variance $sigma^2=langle C^2 rangle - langle C rangle^2$. A consequence is that mixing strongly depends on the mutual coherence between colloid and salt concentration fields, the short time evolution of scalar variance being governed by a new variance production term $P=- langle C^2 abla cdot mathbf{V}_{rm dp} rangle/2$ when scalar gradients are not developed yet so that dissipation is weak. Depending on initial conditions, mixing is then delayed or enhanced, and it is possible to find examples for which the two regimes (fast mixing followed by slow mixing) are observed consecutively when the variance source term reverses its sign. This is indeed the case for localized patches modeled as gaussian concentration profiles.
We present a series of three-dimensional discrete Boltzmann (DB) models for compressible flows in and out of equilibrium. The key formulating technique is the construction of discrete equilibrium distribution function through inversely solving the kinetic moment relations that it satisfies. The crucial physical requirement is that all the used kinetic moment relations must be consistent with the non-equilibrium statistical mechanics. The necessity of such a kinetic model is that, with increasing the complexity of flows, the dynamical characterization of non-equilibrium state and the understanding of the constitutive relations need higher order kinetic moments and their evolution. The DB models at the Euler and Navier-Stokes levels proposed by this scheme are validated by several well-known benchmarks, ranging from one-dimension to three-dimension. Particularly, when the local Mach number, temperature ratio, and pressure ratio are as large as $10^2$, $10^4$, and $10^5$, respectively, the simulation results are still in excellent agreement with the Riemann solutions. How to model deeper thermodynamic non-equilibrium flows by DB is indicated. Via the DB method, it convenient to simulate nonequilibrium flows without knowing exact form of the hydrodynamic equations.
Wing flexibility plays an essential role in the aerodynamic performance of insects due to the considerable deformation of their wings during flight under the impact of inertial and aerodynamic forces. These forces come from the complex wing kinematics of insects. In this study, both wing structural dynamics and flapping wing motion are taken into account to investigate the effect of wing deformation on the aerodynamic efficiency of a bumblebee in tethered flight. A fluid-structure interaction solver, coupling a mass-spring model for the flexible wing with a pseudo-spectral code solving the incompressible Navier-Stokes equations, is implemented for this purpose. We first consider a tethered bumblebee flying in laminar flow with flexible wings. Compared to the rigid model, flexible wings generate smaller aerodynamic forces but require much less power. Finally, the bumblebee model is put into a turbulent flow to investigate its influence on the force production of flexible wings.
Scale-space energy density function, $E(mathbf{x}, mathbf{r})$, is defined as the derivative of the two-point velocity correlation. The function E describes the turbulent kinetic energy density of scale r at a location x and can be considered as the generalization of spectral energy density function concept to inhomogeneous flows. We derive the transport equation for the scale-space energy density function in compressible flows to develop a better understanding of scale-to-scale energy transfer and the degree of non-locality of the energy interactions. Specifically, the effects of variable-density and dilatation on turbulence energy dynamics are identified. It is expected that these findings will yield deeper insight into compressibility effects leading to improved models at all levels of closure for mass flux, density-variance, pressure-dilatation, pressure-strain correlation and dilatational dissipation processes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
