Fluid penetration effects in porous media contact

The contact boundary conditions at the interface between two fluid-saturated porous bodies are derived. The general derivation is performed within the well-founded framework of the Theory of Porous Media (TPM) based on the constituent balance relations of mass, momentum, and energy accounting for finite discontinuities at the contact surface. Particular attention is drawn to the effects associated with the interstitial fluid flux across the interface. The derived contact conditions include two kinematic continuity conditions for the solid velocity and the fluid seepage velocity as well as two jump conditions for the effective solid stress and the pore-fluid pressure. As an application, the common case of biphasic porous media contact proceeding from materially incompressible constituents and inviscid fluid properties is discussed in detail.      

Diffuse interface model for high speed cavitating underwater systems

High speed underwater systems involve many modelling and simulation difficulties related to shocks, expansion waves and evaporation fronts. Modern propulsion systems like underwater missiles also involve extra difficulties related to non-condensable high speed gas flows. Such flows involve many continuous and discontinuous waves or fronts and the difficulty is to model and compute correctly jump conditions across them, particularly in unsteady regime and in multi-dimensions. To this end a new theory has been built that considers the various transformation fronts as ‘diffuse interfaces’. Inside these diffuse interfaces relaxation effects are solved in order to reproduce the correct jump conditions. For example, an interface separating a compressible non-condensable gas and compressible water is solved as a multiphase mixture where stiff mechanical relaxation effects are solved in order to match the jump conditions of equal pressure and equal normal velocities. When an interface separates a metastable liquid and its vapor, the situation becomes more complex as jump conditions involve pressure, velocity, temperature and entropy jumps. However, the same type of multiphase mixture can be considered in the diffuse interface and stiff velocity, pressure, temperature and Gibbs free energy relaxation are used to reproduce the dynamics of such fronts and corresponding jump conditions. A general model, based on multiphase flow theory is thus built. It involves mixture energy and mixture momentum equations together with mass and volume fraction equations for each phase or constituent. For example, in high velocity flows around underwater missiles, three phases (or constituents) have to be considered: liquid, vapor and propulsion gas products. It results in a flow model with 8 partial differential equations. The model is strictly hyperbolic and involves waves speeds that vary under the degree of metastability. When none of the phase is metastable, the non-monotonic sound speed is recovered. When phase transition occurs, the sound speed decreases and phase transition fronts become expansion waves of the equilibrium system. The model is built on the basis of asymptotic analysis of a hyperbolic total non-equilibrium multiphase flow model, in the limit of stiff mechanical relaxation. Closure relations regarding heat and mass transfer are built under the examination of entropy production. The mixture equation of state (EOS) is based on energy conservation and mechanical equilibrium of the mixture. Pure phases EOS are used in the mixture EOS instead of cubic one in order to prevent loss of hyperbolicity in the spinodal zone of the phase diagram. The corresponding model is able to deal with metastable states without using Van der Waals representation.     

A theory of mixtures

This paper is concerned with a dynamical theory of mixtures, composed of n reactive constituents in relative motion to each other. The theory is developed in terms of the constituent ingredients using a balance of energy and an entropy production inequality for each constituent of the mixture, together with invariance requirements under superposed rigid body motions of the whole mixture. The balance of energy and the entropy production inequality for each of the constituents, which include contributions arising from interactions, combine to yield a single energy equation and a single entropy production inequality in terms of the ingredients of the mixture as a whole; the relations between the thermodynamical variables of the mixture and those of its constituents depend, in general, on the past history of the temperature and the kinematic variables. Full thermodynamical restrictions are deduced, and the theory is applied to the special case of a mixture of two ideal fluids.     

Balance relations for classical mixtures containing a moving non-material surface with application to phase transitions

This work is concerned with an extension of classical mixture theory to the case in which the mixture contains an evolving non-material surface on which the constituents may interact, as well as be created and/or annihilated. The formulation of constituent and mixture jump balance relations on/across such a non-material surface proceed by analogy with the standard volume or bulk constituent and mixture balance relations. On this basis, we derive various forms of the constituent mass, momentum, energy and entropy balances assuming (1), that the constituent in question is present on both sides of the moving, non-material surface, and (2), that it is created or annihilated on this surface, as would be the case in a phase transition. In particular, we apply the latter model to the transition between cold and temperate ice found in polythermal ice masses, obtaining in the process the conditions under which melting or freezing takes place at this boundary. On a more general level, one of the most interesting aspects of this formulation is that it gives rise to certain combinations of the limits of constituent and mixture volume fields on the moving mixture interface which can be interpreted as the corresponding surface form of these fields, leading to the possibility of exploiting the surface entropy inequality to obtain restrictions on surface constitutive relations.     

Mechanics and thermodynamics of surface growth viewed as moving discontinuities

Surface growth is presently described as the motion of a moving interface of vanishing thickness, physically representing the generating cells, separating a zone not yet affected by growth from a domain in which growth has occurred. The jump conditions of density, velocity, momentum, energy, and entropy over the moving front are expressed from the general balance laws of open systems in both physical and material format. The writing of the jump of the internal entropy production in material format allows the identification of a driving force for surface growth, thermodynamically conjugated to the material velocity of the moving front.     

Two-layer debris mixture flows on arbitrary terrain with mass exchange at the base and the interface

Granular-fluid gravity driven flow down arbitrary topographic terrain is modelled as a two-layer system of a solid-fluid mixture layer overlain by a slurry, consisting of a particle-laden fluid. The lower layer is dynamically treated as a two-phase flow with two constituent mass and momentum balance laws. By contrast, the slurry is described by mass and momentum balances for the mixture as a whole and a diffusive mass balance for the suspended particle phase. At the base, the mixture interacts with the stagnant base by solid-fluid deposition or erosion. At the mixture-slurry interface, solid and fluid mass exchanges are equally taken into account, but the free surface is treated as material and tractionless. The dynamical equations are formulated in three-dimensional form as general balance laws of mass and momentum in each layer. Intrinsic expressions of the jump conditions of mass and momentum are given for the basal and interface surfaces. The field equations are put into dimensionless form and presented relative to topography adjusted coordinates. These equations are further simplified and approximated by a depth-averaging procedure using an order parameter ${\varepsilon = H/L}$ , where H and L are typical thickness and length scales of the gravity current. Detailed proposals are worked out for the parameterizations of the solid and fluid mass flows across the basal surface and layer interface.     

Hydraulic jumps in two-layer flow with a free surface

This paper presents a closure relation which describes hydraulic jumps in two-layer flows with a free surface over a flat bottom. This relation is derived from the momentum equations for each layer, which, subject to the condition of conservation of the total momentum and mass of each layer, become conservative in a sense. It is shown that use of this relation provides a reduction in the total energy at the jump.     

Generalization of the Stefan model to allow for both velocity and temperature jumps

We obtain an expression for the energy dissipation due to an evolving nonmaterial interface across which the mass density, velocity, stress, energy density, heat flux, entropy density, and temperature may be discontinuous. This expression is a sum of three terms: the product of the interfacial mass flux with the interfacial energy release; the scalar product of the interfacial velocity slip with the interfacial friction; and, the product of the interfacial temperature jump, scaled by the interfacial temperature average, with the interfacial heating. When the surface in question is a phase interface, we propose, on the basis of the interfacial dissipation inequality, supplemental relations that determine the interfacial energy release, the interfacial friction, and the interfacial heating constitutively as functions of the interfacial mass flux, the interfacial velocity slip, and the scaled interfacial temperature jump. As a step toward an understanding of the role that such interfacial relations may serve in theories for phase transitions, we investigate a problem involving the solidification of a pure substance in the absence of flow. Received February 17, 1999     

The concept of material forces in phase transition problems within the level-set framework

The dynamics of a phase transition front in solids using the level set method is examined in this paper. Introducing an implicit representation of singular surfaces, a regularized version of the sharp interface model arises. The interface transforms into a thin transition layer of nonzero thickness where all quantities take inhomogeneous expressions within the body. It is proved that the existence of an inhomogeneous energy of the material predicts inhomogeneity forces that drive the singularity. The driving force is a material force entering the canonical momentum equation (pseudo-momentum) in a natural way. The evolution problem requires a kinetic relation that determines the velocity of the phase transition as a function of the driving force. Here, the kinetic relation is produced by invoking relations that can be considered as the regularized versions of the Rankine–Hugoniot jump conditions. The effectiveness of the method is illustrated in a shape memory alloy bar.     

Jump conditions and entropy sources in two-phase systems. Local instant formulation

The continuum mechanics of two-phase systems involving surface tension and surface properties is discussed.The integral forms of the balance laws are given for the following quantities: mass, linear momentum, angular momentum, total energy and entropy. Starting from these integral balance laws, the jump conditions and the entropy source at the interface are derived.     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CONTINUUM THEORIES (Ⅵ)—CONSERVATION LAWS OF MASS AND INERTIA

IntroductionThisworkisadirectcontinuationandasupplementofRefs .[1～8] .InRefs.[1～8]thecoupledbalancelawsandequationsofmomentum ,angularmomentumandenergyaswellasthenewHamiltonprinciple,principleofvirtualpowerandNoethertheoremhavebeenpresented .However,thecoupledconservationlawsofmassandinertiahavenotbeenreestablishedyet.Thepurposeofthispaperistoreestablishtheconservationlawsandequationsofmassandinertiaandtocombinethemwiththecoupledbalancelawsandequationsofmomentum ,angularmomentum ,energyand…     

Direct and indirect methods of calculating entropy generation rates in turbulent convective heat transfer problems

Computational fluid dynamics (CFD) solutions of turbulent convective heat transfer problems based on the mass, momentum and energy conservation principle provide all information to calculate the entropy production rate in such a transfer process. It can be determined in the post processing phase of a CFD calculation. Two methods are discussed in detail which can provide the information about the entropy production with different degrees of accuracy.     

Two-degree-of-freedom vortex-induced vibrations of a circular cylinder at Re=3900

The vortex-induced vibrations of an elastically mounted circular cylinder are investigated on the basis of direct numerical simulations. The body is free to move in the in-line and cross-flow directions. The natural frequencies of the oscillator are the same in both directions. The Reynolds number, based on the free stream velocity and cylinder diameter, is set to 3900 and kept constant in all simulations. The behavior of the coupled flow-structure system is analyzed over a wide range of the reduced velocity (inverse of the natural frequency) encompassing the lock-in range, i.e. where body motion and flow unsteadiness are synchronized. The statistics of the structural responses and forces are in agreement with prior experimental results. Large-amplitude vibrations develop in both directions. The in-line and cross-flow oscillations are close to harmonic; they exhibit a frequency ratio of 2 and a variable phase difference across the lock-in range. Distinct trends are noted in the force-displacement phasing mechanisms in the two directions: a phase difference jump associated with a sign change of the effective added mass and a vibration frequency crossing the natural frequency is observed in the cross-flow direction, while no phase difference jump occurs in the in-line direction. Higher harmonic components arise in the force spectra; their contributions become predominant when the cylinder oscillates close to the natural frequency. The force higher harmonics are found to impact the transfer of energy between the flow and the moving body, in particular, by causing the emergence of new harmonics in the energy transfer spectrum.     

Discontinuities in helium II

The jump conditions at surface of discontinuity are derived for the two-fluid model of helium II from postulated balance laws for the total energy, the linear momentum of the superfluid and from an entropy production inequality. These conditions are used to discuss a contact surface and the propagation of a weak shock.     

Galilean relativistic fluid mechanics

Single-component nonrelativistic dissipative fluids are treated independently of reference frames and flow-frames. First the basic fields and their balances are derived, then the related thermodynamic relations and the entropy production are calculated and the linear constitutive relations are given. The usual basic fields of mass, momentum, energy and their current densities, the heat flux, pressure tensor and diffusion flux are the time- and spacelike components of the third-order mass–momentum–energy density-flux four-tensor. The corresponding Galilean transformation rules of the physical quantities are derived. It is proved that the non-equilibrium thermodynamic frame theory, including the thermostatic Gibbs relation and extensivity condition and also the entropy production, is independent of the reference frame and also the flow-frame of the fluid. The continuity-Fourier–Navier–Stokes equations are obtained almost in the traditional form if the flow of the fluid is fixed to the temperature. This choice of the flow-frame is the thermo-flow. A simple consequence of the theory is that the relation between the total, kinetic and internal energies is a Galilean transformation rule.     

The intersection marker method for 3D interface tracking of deformable surfaces in finite volumes

Currently, the majority of computational fluid dynamics (CFD) codes use the finite volume method to spatially discretise the computational domain, sometimes as an array of cubic control volumes. The Finite volume method works well with single‐phase flow simulations, but two‐phase flow simulations are more challenging because of the need to track the surface interface traversing and deforming within the 3D grid. Surface area and volume fraction details of each interface cell must be accurately accounted for, in order to calculate for the momentum exchange and rates of heat and mass transfer across the interface. To attain a higher accuracy in two‐phase flow CFD calculations, the intersection marker (ISM) method is developed. The ISM method is a hybrid Lagrangian–Eulerian front‐tracking algorithm that can model an arbitrary 3D surface within an array of cubic control volumes. The ISM method has a cell‐by‐cell remeshing capability that is volume conservative and is suitable for the tracking of complex interface deformation in transient two‐phase CFD simulations. Copyright © 2015 John Wiley & Sons, Ltd.     

Experimental study of the flow regimes resulting from the impact of an intermittent gasoline spray

The present paper reports a complete set of measurements made with a two-component phase Doppler anemometer of the two-phase flow generated at the impact of a transient gasoline spray onto a flat surface. The spray is generated by a pintle injector and the fuel used was gasoline. The measurements of droplet size–velocity were processed to provide time fluxes of number, mass, normal momentum, and energy of the poly-dispersion of droplets ejected at impact, and analyzed based on predictive tools available in the literature. The results show that splash is the dominant mechanism by which secondary droplets are ejected from the surface, either in the stagnation region or in the core region of the spray. In the stagnation region, a large fraction of each incident droplet adheres to the surface and the axial incident momentum contributes with a larger parcel than tangential momentum. As a result, the normal velocity of ejected droplets is much smaller than that of the original incident droplets, while tangential velocity is enhanced. The region near the stagnation point is immediately flooded upon impact of the leading front of the spray, forming a liquid film that is forced to move radially outwards as droplets continue to impinge during the steady period. Spray/wall interaction in the core region thus occurs in the presence of a moving thin liquid film, which enhances transfer of tangential momentum. As a result, film spreading and dynamics as a result of impingement forces are crucial to accurate model spray/wall interaction. The outer region of the spray is dominated by the vortical structure induced by shear forces, which entrains small responsive secondary droplets to re-impinge. Furthermore, prediction of the outcome of spray impact requires a precise knowledge of the two-phase flow in the presence of the target.