On the <Emphasis Type="Italic">H</Emphasis><Superscript><Emphasis Type="Italic">s</Emphasis></Superscript> Theory of Hydrostatic Euler Equations

In this paper we study the two-dimensional hydrostatic Euler equations in a periodic channel. We prove the local existence and uniqueness of H ^s solutions under the local Rayleigh condition. This extends Brenier’s (Nonlinearity 12(3):495–512, 1999) existence result by removing an artificial condition and proving uniqueness. In addition, we prove weak–strong uniqueness, mathematical justification of the formal derivation and stability of the hydrostatic Euler equations. These results are based on weighted H ^s a priori estimates, which come from a new type of nonlinear cancellation between velocity and vorticity.     

Adjoint and Compensated Compactness Methods for Hamilton–Jacobi PDE

We investigate the vanishing viscosity limit for Hamilton–Jacobi PDE with nonconvex Hamiltonians, and present a new method to augment the standard viscosity solution approach. The main idea is to introduce a solution σ ^ε of the adjoint of the formal linearization, and then to integrate by parts with respect to the density σ ^ε . This procedure leads to a natural phase space kinetic formulation and also to a new compensated compactness technique.     

Large Eddy Simulation of an Initially-Confined Triangular Oscillating Jet

This paper reports the first large eddy simulation (LES) of a self-excited oscillating triangular jet (OTJ) issuing from a fluidic nozzle that consists of a small triangular orifice inlet followed by a large circular chamber and an orifice outlet. The case simulated is identical to that measured experimentally by England et al. (Exp Fluids 48(1):69–80, 2010). The present prediction agrees well with the previous measurement. The simulation reveals that the central oscillating jet exhibits axis-switching in the cross-section and rotates by 60° approximately over a downstream distance of x = 0.5D (chamber diameter). Three strong longitudinal vortices occur associated with the three vertices of the inlet triangle. These vortices strongly interact with the central jet and also the surroundings, in the region at x/D ≤ 1, and appear to merge finally with the outer secondary swirling flow. These observations are consistent with the deduction from previous experiments.     

The Erpenbeck High Frequency Instability Theorem for Zeldovitch–von Neumann–D?ring Detonations

The rigorous study of spectral stability for strong detonations was begun by Erpenbeck (Phys. Fluids 5:604–614 1962). Working with the Zeldovitch–von Neumann–D?ring (ZND) model (more precisely, Erpenbeck worked with an extension of ZND to general chemistry and thermodynamics), which assumes a finite reaction rate but ignores effects such as viscosity corresponding to second order derivatives, he used a normal mode analysis to define a stability function V(t,e){V(\tau,\epsilon)} whose zeros in ${\mathfrak{R}\tau > 0}${\mathfrak{R}\tau > 0} correspond to multidimensional perturbations of a steady detonation profile that grow exponentially in time. Later in a remarkable paper (Erpenbeck in Phys. Fluids 9:1293–1306, 1966; Stability of detonations for disturbances of small transverse wavelength, 1965) he provided strong evidence, by a combination of formal and rigorous arguments, that for certain classes of steady ZND profiles, unstable zeros of V exist for perturbations of sufficiently large transverse wavenumber e{\epsilon} , even when the von Neumann shock, regarded as a gas dynamical shock, is uniformly stable in the sense defined (nearly 20 years later) by Majda. In spite of a great deal of later numerical work devoted to computing the zeros of V(t,e){V(\tau,\epsilon)} , the paper (Erpenbeck in Phys. Fluids 9:1293–1306, 1966) remains one of the few works we know of [another is Erpenbeck (Phys. Fluids 7:684–696, 1964), which considers perturbations for which the ratio of longitudinal over transverse components approaches ∞] that presents a detailed and convincing theoretical argument for detecting them. The analysis in Erpenbeck (Phys. Fluids 9:1293–1306, 1966) points the way toward, but does not constitute, a mathematical proof that such unstable zeros exist. In this paper we identify the mathematical issues left unresolved in Erpenbeck (Phys. Fluids 9:1293–1306, 1966) and provide proofs, together with certain simplifications and extensions, of the main conclusions about stability and instability of detonations contained in that paper. The main mathematical problem, and our principal focus here, is to determine the precise asymptotic behavior as e?￥{\epsilon\to\infty} of solutions to a linear system of ODEs in x, depending on e{\epsilon} and a complex frequency τ as parameters, with turning points x _* on the half-line [0,∞).     

Derivation of a Darcy’s Law for a Porous Medium Composed of Two Solid Phases Saturated by a Single- Phase Fluid: A Homogenization Approach

The objective of this article is to derive a macroscopic Darcy’s law for a fluid-saturated moving porous medium whose matrix is composed of two solid phases which are not in direct contact with each other (weakly coupled solid phases). An example of this composite medium is the case of a solid matrix, unfrozen water, and an ice matrix within the pore space. The macroscopic equations for this type of saturated porous material are obtained using two-space homogenization techniques from microscopic periodic structures. The pore size is assumed to be small compared to the macroscopic scale under consideration. At the microscopic scale the two weakly coupled solids are described by the linear elastic equations, and the fluid by the linearized Navier–Stokes equations with appropriate boundary conditions at the solid–fluid interfaces. The derived Darcy’s law contains three permeability tensors whose properties are analyzed. Also, a formal relation with a previous macroscopic fluid flow equation obtained using a phenomenological approach is given. Moreover, a constructive proof of the existence of the three permeability tensors allows for their explicit computation employing finite elements or analogous numerical procedures.     

Dynamics of a Non-Autonomous ODE System Occurring in Coagulation Theory

We consider a constant coefficient coagulation equation with Becker–D?ring type interactions and power law input of monomers J ₁(t) = α t ^ω, with α > 0 and . For this infinite dimensional system we prove solutions converge to similarity profiles as t and j converge to infinity in a similarity way, namely with either or constants, where is a function of t only. This work generalizes to the non-autonomous case a recent result of da Costa et al. (2004). Markov Processes Relat. Fields 12, 367–398. and provides a rigorous derivation of formal results obtained by Wattis J. Phys. A: Math. Gen. 37, 7823–7841. The main part of the approach is the analysis of a bidimensional non-autonomous system obtained through an appropriate change of variables; this is achieved by the use of differential inequalities and qualitative theory methods. The results about rate of convergence of solutions of the bidimensional system thus obtained are fed into an integral formula representation for the solutions of the infinite dimensional system which is then estimated by an adaptation of methods used by da Costa et al. (2004). Markov Processes Relat. Fields 12, 367–398.      

Weak in Space, Log in Time Improvement of the Lady?enskaja�CProdi�CSerrin Criteria

In this article we present a Ladyženskaja–Prodi–Serrin Criteria for regularity of solutions for the Navier–Stokes equation in three dimensions which incorporates weak L ^p norms in the space variables and log improvement in the time variable.     

Laminar-turbulent transition in the boundary layer on cones in a hypersonic flow at high Reynolds numbers per meter

The laminar-turbulent transition is experimentally studied in boundary-layer flows on cones with a rectangular axisymmetric step in the base part of the cone and without the step. The experiments are performed in an A-1 two-step piston-driven gas-dynamic facility with adiabatic compression of the working gas with Mach numbers at the nozzle exit M _∞ = 12–14 and pressures in the settling chamber P₀ = 60–600 MPa. These values of parameters allow obtaining Reynolds numbers per meter near the cone surface equal to Re _1e = (53–200) · 10⁶ m ⁻¹. The transition occurs at Reynolds numbers Re _tr = (2.3–5.7) · 10⁶. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 76–83, May–June, 2007.     

Insights into the Relationships Among Capillary Pressure,Saturation, Interfacial Area and Relative Permeability Using Pore-Network Modeling

To gain insight in relationships among capillary pressure, interfacial area, saturation, and relative permeability in two-phase flow in porous media, we have developed two types of pore-network models. The first one, called tube model, has only one element type, namely pore throats. The second one is a sphere-and-tube model with both pore bodies and pore throats. We have shown that the two models produce distinctly different curves for capillary pressure and relative permeability. In particular, we find that the tube model cannot reproduce hysteresis. We have investigated some basic issues such as effect of network size, network dimension, and different trapping assumptions in the two networks. We have also obtained curves of fluid–fluid interfacial area versus saturation. We show that the trend of relationship between interfacial area and saturation is largely influenced by trapping assumptions. Through simulating primary and scanning drainage and imbibition cycles, we have generated two surfaces fitted to capillary pressure, saturation, and interfacial area (P _c –S _w –a _nw ) points as well as to relative permeability, saturation, and interfacial area (k _r –S _w –a _nw ) points. The two fitted three-dimensional surfaces show very good correlation with the data points. We have fitted two different surfaces to P _c –S _w –a _nw points for drainage and imbibition separately. The two surfaces do not completely coincide. But, their mean absolute difference decreases with increasing overlap in the statistical distributions of pore bodies and pore throats. We have shown that interfacial area can be considered as an essential variable for diminishing or eliminating the hysteresis observed in capillary pressure–saturation (P _c –S _w ) and the relative permeability–saturation (k _r –S _w ) curves.     

Stability for the 3D Navier–Stokes Equations with Nonzero far Field Velocity on Exterior Domains

Concerning to the non-stationary Navier–Stokes flow with a nonzero constant velocity at infinity, just a few results have been obtained, while most of the results are for the flow with the zero velocity at infinity. The temporal stability of stationary solutions for the Navier–Stokes flow with a nonzero constant velocity at infinity has been studied by Enomoto and Shibata (J Math Fluid Mech 7:339–367, 2005), in L ^p spaces for p ≥ 3. In this article, we first extend their result to the case \frac32 < p{\frac{3}{2} < p} by modifying the method in Bae and Jin (J Math Fluid Mech 10:423–433, 2008) that was used to obtain weighted estimates for the Navier–Stokes flow with the zero velocity at infinity. Then, by using our generalized temporal estimates we obtain the weighted stability of stationary solutions for the Navier–Stokes flow with a nonzero velocity at infinity.     

Onset of Buoyancy-driven Instability in Porous Media Melted from Below

When a nonhomogeneous solid is melting from below, convection may be induced in a thermally–unstable melt layer. In this study, the onset of buoyancy-driven convection during time-dependent melting is investigated by using similarly transformed disturbance equations. The critical Darcy–Rayleigh numbers based on the melt-layer thickness, Ra _H,c, are found numerically for various conditions. For small superheats, the present predictions show that Ra _H,c is located between 27.1 and 4π ² and it approaches the well-known results of the original Horton–Rogers–Lapwood problem. However, for high superheats, it is dependent on the phase change rate λ and the relation of Ra _H,c λ = 25.89 is shown.     

Deduction of Saint-Venant"s Principle for the Asymptotic Residue of Elastic Rods

Let u(ε) be a rescaled 3-dimensional displacement field solution of the linear elastic model for a free prismatic rod Ω^ε having cross section with diameter of order ε, and let u ⁽⁰⁾ –Bernoulli–Navier displacement – and u ⁽²⁾ be the two first terms derived from the asymptotic method. We analyze the residue r(ε) = u(ε) − (u ⁽⁰⁾ + ε² u ⁽²⁾) and if the cross section is star-shaped, we prove such residue presents a Saint-Venant"s phenomenon near the ends of the rod. This revised version was published online in August 2006 with corrections to the Cover Date.     

Specific features of microparticle deformation upon impact on a rigid barrier

Based on experimental data and numerical modeling, it is shown that a lamina of melted metal of thickness of order0.01 d, in which the temperature is close to the melting point of the particle material, can be formed upon high-speed impact (v ₀≈500–1200 m/sec) of a fine metal particle (d=1–50 μm) on a rigid undeformable barrier near the contact surface. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 1, pp. 204–209, January–February, 2000.     

Multibody Formulation with Large Bending Finite Elements (LBFE)

The goal of this paper is to present a flexible multibody formulation for Euler-Bernoulli beams involving large displacements. This method is based on a discretisation of internal and kinetic energies. The beam is represented by its line of centroids and each section is oriented by a frame defined by three Euler angles. We apply a finite element formulation to describe the evolution of these angles along the neutral fibre and describe the internal energy. The kinetic energy is approximated as the one of two rigid bars tangent to the neutral fibre at the ends of the beam. We derive the equations of motion from a Lagrange formulation. These equations are solved using the Newmark method or/and the Newton-Raphson technique. We solve some very classic problems taken from the literature as the curved beam presented by Simo [Simo, J. C., ‘A three-dimensional finite-strain rod model. the three-dimensional dynamic problem. Part I’, Comput. Meths. Appl. Mech. Engrg. 49, 1985, 55–70; Simo, J. C. and Vu-Quoc, L., ‘A three-dimensional finite-strain rod model, Part II: Computationals aspects’, Comput. Meths. Appl. Mech. Engrg. 58, 1988, 79–116] and Lee [Lee, Kisu, ‘Analysis of large displacements and large rotations of three-dimensional beams by using small strains and unit vectors’, Commun. Numer. Meth. Engrg. 13, 1997, 987–997] or the rotational rod presented by Avello [Avello, A., Garcia de Jalon, J., and Bayo, E., ‘Dynamics of flexible multibody systems using cartesian co-ordinates and large displacement theory’, Int. J. Num. Methods in Engineering 32, 1991, 1543–1563] and Simo [Simo, J. C. and Vu-Quoc, L., ‘On the dynamics of flexible beams under large overall motions – the planar case. Part I’ Jour. of Appl. Mech. 53, 1986, 849–854; Simo, J. C. and Vu-Quoc, L., ‘On the dynamics of flexible beams under large overall motions – the planar case. Part II’, Jour. of Appl. Mech. 53, 1986, 855–863].     

Testing of Model Equations for the Mean Dissipation using Kolmogorov Flows

Direct Numerical Simulations (DNS) of Kolmogorov flows are performed at three different Reynolds numbers Re _λ between 110 and 190 by imposing a mean velocity profile in y-direction of the form U(y) = F sin(y) in a periodic box of volume (2π)³. After a few integral times the turbulent flow turns out to be statistically steady. Profiles of mean quantities are then obtained by averaging over planes at constant y. Based on these profiles two different model equations for the mean dissipation ε in the context of two-equation RANS (Reynolds Averaged Navier–Stokes) modelling of turbulence are compared to each other. The high Reynolds number version of the k-ε-model (Jones and Launder, Int J Heat Mass Transfer 15:301–314, 1972), to be called the standard model and a new model by Menter et al. (2006), to be called the Menter–Egorov model, are tested against the DNS results. Both models are solved numerically and it is found that the standard model does not provide a steady solution for the present case, while the Menter–Egorov model does. In addition a fairly good quantitative agreement of the model solution and the DNS data is found for the averaged profiles of the kinetic energy k and the dissipation ε. Furthermore, an analysis based on flow-inherent geometries, called dissipation elements (Wang and Peters, J Fluid Mech 608:113–138, 2008), is used to examine the Menter–Egorov ε model equation. An expression for the evolution of ε is derived by taking appropriate moments of the equation for the evolution of the probability density function (pdf) of the length of dissipation elements. A term-by-term comparison with the model equation allows a prediction of the constants, which with increasing Reynolds number approach the empirical values.     

Topological Invariants in a Model of a Time-Delayed Chua's Circuit

In the last 30 years, some authors have been studying several classes of boundary value problems (BVP) for partial differential equations (PDE) using the method of reduction to obtain a difference equation with continuous argument which behavior is determined by the iteration of a one-dimensional (1D) map (see, for example, Romanenko, E. Yu. and Sharkovsky, A. N., International Journal of Bifurcation and Chaos 9(7), 1999, 1285–1306; Sharkovsky, A. N., International Journal of Bifurcation and Chaos 5(5), 1995, 1419–1425; Sharkovsky, A. N., Analysis Mathematica Sil 13, 1999, 243–255; Sharkovsky, A. N., in “New Progress in Difference Equations”, Proceedings of the ICDEA'2001, Taylor and Francis, 2003, pp. 3–22; Sharkovsky, A. N., Deregel, Ph., and Chua, L. O., International Journal of Bifurcation and Chaos 5(5), 1995, 1283–1302; Sharkovsky, A. N., Maistrenko, Yu. L., and Romanenko, E. Yu., Difference Equations and Their Applications, Kluwer, Dordrecht, 1993.). In this paper we consider the time-delayed Chua's circuit introduced in (Sharkovsky, A. N., International Journal of Bifurcation and Chaos 4(5), 1994, 303–309; Sharkovsky, A. N., Maistrenko, Yu. L., Deregel, Ph., and Chua, L. O., Journal of Circuits, Systems and Computers 3(2), 1993, 645–668.) which behavior is determined by properties of one-dimensional map, see Sharkovsky, A. N., Deregel, Ph., and Chua, L. O., International Journal of Bifurcation and Chaos 5(5), 1995, 1283–1302; Maistrenko, Yu. L., Maistrenko, V. L., Vikul, S. I., and Chua, L. O., International Journal of Bifurcation and Chaos 5(3), 1995, 653–671; Sharkovsky, A. N., International Journal of Bifurcation and Chaos 4(5), 1994, 303–309; Sharkovsky, A. N., Maistrenko, Yu. L., Deregel, Ph., and Chua, L. O., Journal of Circuits, Systems and Computers 3(2), 1993, 645–668. To characterize the time-evolution of these circuits we can compute the topological entropy and to distinguish systems with equal topological entropy we introduce a second topological invariant.     

A Numerical Study Promoting Algebraic Models for the Lewis Number Effect in Atmospheric Turbulent Premixed Bunsen Flames

This numerical investigation carried out on turbulent lean premixed flames accounts for two algebraic – the Lindstedt–Vaos (LV) and the classic Bray–Moss–Libby (BML) – reaction rate models. Computed data from these two models is compared with the experimental data of Kobayashi et al. on 40 different methane, ethylene and propane Bunsen flames at 1 bar, where the mean flame cone angle is used for comparison. Both models gave reasonable qualitative trend for the whole set of data, in overall. In order to characterize quantitatively, firstly, corrections are made by tuning the model parameters fitting to the experimental methane–air (of Le = 1.0) flame data. In case of the LV model, results obtained by adjusting the pre-constant, i.e., reaction rate parameter, C_R, from its original value 2.6 to 4.0, has proven to be in good agreement with the experiments. Similarly, for the BML model, with the tuning of the exponent n, in the wrinkling length scale, L_y = C_l⋅ l_x(s_L/u′)ⁿ from value unity to 1.2, the outcome is in accordance with the measured data. The deviation between the measured and calculated data sharply rises from methane to propane, i.e., with increasing Lewis number. It is deduced from the trends that the effect of Lewis number (for ethylene–air mixtures of Le = 1.2 and propane–air mixtures of Le = 1.62) is missing in both the models. The Lewis number of the fuel–air mixture is related to the laminar flame instabilities. Second, in order to quantify for its influence, the Lewis number effect is induced into both the models. It is found that by setting global reaction rate inversely proportional to the Lewis number in both the cases leads to a much better numerical prediction to this set of experimental flame data. Thus, by imparting an important phenomenon (the Lewis number effect) into the reaction rates, the generality of the two models is enhanced. However, functionality of the two models differs in predicting flame brush thickness, giving scope for further analysis.     

The Regularity of Weak Solutions of the 3D Navier–Stokes Equations in {B^{-1}_{\infty,\infty}}

We show that if a Leray–Hopf solution u of the three-dimensional Navier–Stokes equation belongs to C((0,T]; B^-1_￥,￥){C((0,T]; B^{-1}_{\infty,\infty})} or its jumps in the B^-1_￥,￥{B^{-1}_{\infty,\infty}}-norm do not exceed a constant multiple of viscosity, then u is regular for (0, T]. Our method uses frequency local estimates on the nonlinear term, and yields an extension of the classical Ladyzhenskaya–Prodi–Serrin criterion.     

A Two-Scale Model for Coupled Electro-Chemo-Mechanical Phenomena and Onsager’s Reciprocity Relations in Expansive Clays: II Computational Validation

In Part I Moyne and Murad [Transport in Porous Media 62, (2006), 333–380] a two-scale model of coupled electro-chemo-mechanical phenomena in swelling porous media was derived by a formal asymptotic homogenization analysis. The microscopic portrait of the model consists of a two-phase system composed of an electrolyte solution and colloidal clay particles. The movement of the liquid at the microscale is ruled by the modified Stokes problem; the advection, diffusion and electro-migration of monovalent ions Na⁺ and Cl⁻ are governed by the Nernst–Planck equations and the local electric potential distribution is dictated by the Poisson problem. The microscopic governing equations in the fluid domain are coupled with the elasticity problem for the clay particles through boundary conditions on the solid–fluid interface. The up-scaling procedure led to a macroscopic model based on Onsager’s reciprocity relations coupled with a modified form of Terzaghi’s effective stress principle including an additional swelling stress component. A notable consequence of the two-scale framework are the new closure problems derived for the macroscopic electro-chemo-mechanical parameters. Such local representation bridge the gap between the macroscopic Thermodynamics of Irreversible Processes and microscopic Electro-Hydrodynamics by establishing a direct correlation between the magnitude of the effective properties and the electrical double layer potential, whose local distribution is governed by a microscale Poisson–Boltzmann equation. The purpose of this paper is to validate computationally the two-scale model and to introduce new concepts inherent to the problem considering a particular form of microstructure wherein the clay fabric is composed of parallel particles of face-to-face contact. By discretizing the local Poisson–Boltzmann equation and solving numerically the closure problems, the constitutive behavior of the diffusion coefficients of cations and anions, chemico-osmotic and electro-osmotic conductivities in Darcy’s law, Onsager’s parameters, swelling pressure, electro-chemical compressibility, surface tension, primary/secondary electroviscous effects and the reflection coefficient are computed for a range particle distances and sat concentrations.