No abstract

SUMMARY A semi‐implicit finite element scheme and a Newton‐like solver are developed for the stationary compressible Euler equations. Since the Galerkin discretization of the inviscid fluxes is potentially oscillatory and unstable, the troublesome antidiffusive part is constrained within the framework of algebraic flux correction. A generalization of total variation diminishing (TVD) schemes is employed to blend the original Galerkin scheme with its nonoscillatory low‐order counterpart. Unlike standard TVD limiters, the proposed limiting strategy is fully multidimensional and readily applicable to unstructured meshes. However, the nonlinearity and nondifferentiability of the limiter function makes efficient computation of stationary solutions a highly challenging task, especially in situations when the Mach number is large in some subdomains and small in other subdomains. In this paper, a semi‐implicit scheme is derived via a time‐lagged linearization of the Jacobian operator, and a Newton‐like method is obtained in the limit of infinite CFL numbers. Special emphasis is laid on the numerical treatment of weakly imposed characteristic boundary conditions. A boundary Riemann solver is used to avoid unphysical boundary states. It is shown that the proposed approach offers unconditional stability, as well as higher accuracy and better convergence behavior than algorithms in which the boundary conditions are implemented in a strong sense. The overall spatial accuracy of the constrained scheme and the benefits of the new boundary treatment are illustrated by grid convergence studies for 2D benchmark problems. Copyright © 2011 John Wiley & Sons, Ltd.

a b s t r a c tThe derivation of macroscopic models for particle-laden gas flows is reviewed. Semiimplicit and Newton-like finite element methods are developed for the stationary two-fluid model governing compressible particle-laden gas flows. The Galerkin discretization of the inviscid fluxes is potentially oscillatory and unstable. To suppress numerical oscillations, the spatial discretization is performed by a high-resolution finite element scheme based on algebraic flux correction. A multidimensional limiter of TVD type is employed. An important goal is the efficient computation of stationary solutions in a wide range of Mach numbers. This is a challenging task due to oscillatory correction factors associated with TVD-type flux limiters and the additional strong nonlinearity caused by interfacial coupling terms. A semi-implicit scheme is derived by a time-lagged linearization of the nonlinear residual, and a Newton-like method is obtained in the limit of infinite CFL numbers. The original Jacobian is replaced by a low-order approximation. Special emphasis is laid on the numerical treatment of weakly imposed boundary conditions. It is shown that the proposed approach offers unconditional stability and faster convergence rates for increasing CFL numbers. The strongly coupled solver is compared to operator splitting techniques, which are shown to be less robust.

Digital rock physics combines microtomographic imaging with advanced numerical simulations of effective material properties. It is used to complement laboratory investigations with the aim to gain a deeper understanding of relevant physical processes related to transport and effective mechanical properties. We apply digital rock physics to reticulite, a natural mineral with a strong analogy to synthetic open-cell foams. We consider reticulite an end-member for high-porosity materials with a high stiffness and brittleness. For this specific material, hydro-mechanical experiments are very difficult to perform. Reticulite is a pyroclastic rock formed during intense Hawaiian fountaining events. the honeycombed network of bubbles is supported by glassy threads and forms a structure with a porosity of more than 80%. Comparing experimental with numerical results and theoretical estimates, we demonstrate the high potential of in situ characterization with respect to the investigation of effective material properties. We show that a digital rock physics workflow, so far applied to conventional rocks, yields reasonable results for high-porosity rocks and can be adopted for fabricated foam-like materials with similar properties. Numerically determined porosities, effective elastic properties, thermal conductivities and permeabilities of reticulite show a fair agreement to experimental results that required exeptionally high experimental efforts.

<p>Microtomographic imaging techniques and advanced numerical simulations are combined by digital rock physics (DRP) to obtain effective physical material properties. The numerical results are typically used to complement laboratory investigations with the aim to gain a deeper understanding of physical processes related to transport (e.g. permeability and thermal conductivity) and effective elastic properties (e.g. bulk and shear modulus). The present study focuses on DRP and laboratory techniques applied to a rock called reticulite, which is considered as an end-member material with respect to porosity, stiffness and brittleness of the skeleton. Classical laboratory investigations on effective properties, such as ultrasonic transmission measurements and uniaxial deformation experiments, are very difficult to perform on this class of high-porosity and brittle materials.</p><p>Reticulite is a pyroclastic rock formed during intense Hawaiian fountaining events. The open honeycombed network has a porosity of more than 80 % and consists of bubbles that are supported by glassy threads. The natural mineral has a strong analogy to fabricated open-cell foams. By comparing experimental with numerical results and theoretical estimates we demonstrate the potential of digital material methodology with respect to the investigation of porosity, effective elastic properties, thermal conductivity and permeability</p><p>We show that the digital rock physics workflow, previously applied to conventional rock types, yields reasonable results for a high-porosity rock and can be adopted for fabricated foam-like materials. Numerically determined effective properties of reticulite are in good agreement with the experimentally determined results. Depending on the fields of application, numerical methods as well as theoretical estimates can become reasonable alternatives to laboratory methods for high porous foam-like materials.</p>