13,433 research outputs found
A Fully Equivalent Global Pressure Formulation for Three-Phase Compressible Flow
We introduce a new global pressure formulation for immiscible three-phase
compressible flows in porous media which is fully equivalent to the original
equations, unlike the one introduced in \cite{CJ86}. In this formulation, the
total volumetric flow of the three fluids and the global pressure follow a
classical Darcy law, which simplifies the resolution of the pressure equation.
However, this global pressure formulation exists only for Total Differential
(TD) three-phase data, which depend only on two functions of saturations and
global pressure: the global capillary pressure and the global mobility. Hence
we introduce a class of interpolation which constructs such TD-three-phase data
from any set of three two-phase data (for each pair of fluids) which satisfy a
TD-compatibility condition
A Fully Equivalent Global Pressure Formulation for Three-Phase Compressible Flow
We introduce a new global pressure formulation for immiscible three-phase compressible flows in porous media which is fully equivalent to the original equations, unlike the one introduced in \cite{CJ86}. In this formulation, the total volumetric flow of the three fluids and the global pressure follow a classical Darcy law, which simplifies the resolution of the pressure equation. However, this global pressure formulation exists only for Total Differential (TD) three-phase data, which depend only on two functions of saturations and global pressure: the global capillary pressure and the global mobility. Hence we introduce a class of interpolation which constructs such TD-three-phase data from any set of three two-phase data (for each pair of fluids) which satisfy a TD-compatibility condition
A Numerical Study of Methods for Moist Atmospheric Flows: Compressible Equations
We investigate two common numerical techniques for integrating reversible
moist processes in atmospheric flows in the context of solving the fully
compressible Euler equations. The first is a one-step, coupled technique based
on using appropriate invariant variables such that terms resulting from phase
change are eliminated in the governing equations. In the second approach, which
is a two-step scheme, separate transport equations for liquid water and vapor
water are used, and no conversion between water vapor and liquid water is
allowed in the first step, while in the second step a saturation adjustment
procedure is performed that correctly allocates the water into its two phases
based on the Clausius-Clapeyron formula. The numerical techniques we describe
are first validated by comparing to a well-established benchmark problem.
Particular attention is then paid to the effect of changing the time scale at
which the moist variables are adjusted to the saturation requirements in two
different variations of the two-step scheme. This study is motivated by the
fact that when acoustic modes are integrated separately in time (neglecting
phase change related phenomena), or when sound-proof equations are integrated,
the time scale for imposing saturation adjustment is typically much larger than
the numerical one related to the acoustics
Lattice Boltzmann Methods for thermal flows: continuum limit and applications to compressible Rayleigh-Taylor systems
We compute the continuum thermo-hydrodynamical limit of a new formulation of
lattice kinetic equations for thermal compressible flows, recently proposed in
[Sbragaglia et al., J. Fluid Mech. 628 299 (2009)]. We show that the
hydrodynamical manifold is given by the correct compressible Fourier-
Navier-Stokes equations for a perfect fluid. We validate the numerical
algorithm by means of exact results for transition to convection in
Rayleigh-B\'enard compressible systems and against direct comparison with
finite-difference schemes. The method is stable and reliable up to temperature
jumps between top and bottom walls of the order of 50% the averaged bulk
temperature. We use this method to study Rayleigh-Taylor instability for
compressible stratified flows and we determine the growth of the mixing layer
at changing Atwood numbers up to At ~ 0.4. We highlight the role played by the
adiabatic gradient in stopping the mixing layer growth in presence of high
stratification and we quantify the asymmetric growth rate for spikes and
bubbles for two dimensional Rayleigh- Taylor systems with resolution up to Lx
\times Lz = 1664 \times 4400 and with Rayleigh numbers up to Ra ~ 2 \times
10^10.Comment: 26 pages, 13 figure
A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting
An explicit moving boundary method for the numerical solution of
time-dependent hyperbolic conservation laws on grids produced by the
intersection of complex geometries with a regular Cartesian grid is presented.
As it employs directional operator splitting, implementation of the scheme is
rather straightforward. Extending the method for static walls from Klein et
al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme
calculates fluxes needed for a conservative update of the near-wall cut-cells
as linear combinations of standard fluxes from a one-dimensional extended
stencil. Here the standard fluxes are those obtained without regard to the
small sub-cell problem, and the linear combination weights involve detailed
information regarding the cut-cell geometry. This linear combination of
standard fluxes stabilizes the updates such that the time-step yielding
marginal stability for arbitrarily small cut-cells is of the same order as that
for regular cells. Moreover, it renders the approach compatible with a wide
range of existing numerical flux-approximation methods. The scheme is extended
here to time dependent rigid boundaries by reformulating the linear combination
weights of the stabilizing flux stencil to account for the time dependence of
cut-cell volume and interface area fractions. The two-dimensional tests
discussed include advection in a channel oriented at an oblique angle to the
Cartesian computational mesh, cylinders with circular and triangular
cross-section passing through a stationary shock wave, a piston moving through
an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil
profile.Comment: 30 pages, 27 figures, 3 table
Recommended from our members
Numerical investigation of high-speed droplet impact using a multiscale two-fluid approach
A single droplet impact onto solid surfaces remains a fundamental and challenging topic in both experimental and numerical studies with significant importance in a plethora of industrial applications, ranging from printing technologies to fuel injection in internal combustion engines. Under high-speed impact conditions, additional complexities arise as a result of the prompt droplet splashing and the subsequent violent fragmentation; thus, different flow regimes and a vast spectrum of sizes for the produced secondary flow structures coexist in the flow field. The present work introduces a numerical methodology to capture the multiscale processes involved with respect to local topological characteristics. The proposed methodology concerns a compressible Σ-Υ two-fluid model with dynamic interface sharpening based on an advanced flow topology detection algorithm. The model has been developed in OpenFOAM® and provides the flexibility of dealing with the multiscale character of droplet splashing, by switching between a sharp and a diffuse interface within the Eulerian-Eulerian framework in segregated and dispersed flow regions, respectively. An additional transport equation for the interface surface area density (Σ) introduces important information for the sub-grid scale phenomena, which is exploited in the dispersed flow regions to provide an insight into the extended cloud of secondary droplets after impact on the target. A high-speed water droplet impact case has been examined and evaluated against new experimental data; these refer to a millimetre size droplet impacting a solid dry smooth surface at velocity as high as 150m/s, which corresponds to a Weber number of ~7.6×10^5. At the investigated impact conditions compressibility effects dominate the early stages of droplet splashing. A strong shock wave forms and propagates inside the droplet, where transonic Mach numbers occur; local Mach numbers up to 2.5 are observed for the expelled surrounding gas outside the droplet. The proposed numerical approach is found to capture relatively accurately the phenomena and provide significant information regarding the produced flow structure dimensions, which is not available from the experiment
An algebraic approach to modeling distributed multiphysics problems: The case of a DRI reactor
© 2015, IFAC (International Federation of Automatic Control) Hosting by Elsevier Ltd. All rights reserved.This paper deals with the problem of modelling a chemical reactor for the Direct Reduction of Iron ore (DRI). Such a process is being increasingly promoted as a more viable alternative to the classic Blast Furnace for the production of iron from raw minerals. Due to the inherent complexity of the process and the reactor itself, its effective monitoring and control requires advanced mathematical models containing distributed-parameter components. While classical approaches such as Finite Element or Finite Differences are still reasonable options, for accuracy and computational efficiency reasons, an algebraic approach is proposed. A full multi-physical, albeit one-dimensional model is addressed and its accuracy is analysed
- …