Temporal Homogenization of Linear ODEs,with Applications to Parametric Super-Resonance and Energy Harvest

We consider the temporal homogenization of linear ODEs of the form \({\dot{x}=Ax+\epsilon P(t)x+f(t)}\), where P(t) is periodic and \({\epsilon}\) is small. Using a 2-scale expansion approach, we obtain the long-time approximation \({x(t)\approx {\rm exp}(At) \left( \Omega(t)+\int_0^t {\rm exp}(-A \tau) f(\tau) {\rm d}\tau \right)}\), where \({\Omega}\) solves the cell problem \({\dot{\Omega}=\epsilon B \Omega + \epsilon F(t)}\) with an effective matrix B and an explicitly-known F(t). We provide necessary and sufficient conditions for the accuracy of the approximation (over a \({{\mathcal{O}}(\epsilon^{-1})}\) time-scale), and show how B can be computed (at a cost independent of \({\epsilon}\)). As a direct application, we investigate the possibility of using RLC circuits to harvest the energy contained in small scale oscillations of ambient electromagnetic fields (such as Schumann resonances). Although a RLC circuit parametrically coupled to the field may achieve such energy extraction via parametric resonance, its resistance R needs to be smaller than a threshold \({\kappa}\) proportional to the fluctuations of the field, thereby limiting practical applications. We show that if n RLC circuits are appropriately coupled via mutual capacitances or inductances, then energy extraction can be achieved when the resistance of each circuit is smaller than \({n\kappa}\). Hence, if the resistance of each circuit has a non-zero fixed value, energy extraction can be made possible through the coupling of a sufficiently large number n of circuits (\({n\approx 1000}\) for the first mode of Schumann resonances and contemporary values of capacitances, inductances and resistances). The theory is also applied to the control of the oscillation amplitude of a (damped) oscillator.     

Topological Pressure of Generic Points for Amenable Group Actions

Let (X, G) be a G-action topological dynamical system (t.d.s. for short), where G is a countably infinite discrete amenable group. In this paper, we study the topological pressure of the sets of generic points. We show that when the system satisfies the almost specification property, for any G-invariant measure \(\mu \) and any continuous map \(\varphi \),

$$\begin{aligned} P\left( X_{\mu },\varphi ,\{F_n\}\right) = h_{\mu }(X)+\int \varphi d\mu , \end{aligned}$$

where \(\{F_n\}\) is a Følner sequence, \(X_{\mu }\) is the set of generic points of \(\mu \) with respect to (w.r.t. for short) \(\{F_n\}\), \(P(X_{\mu },\varphi ,\{F_n\})\) is the topological pressure of \(X_{\mu }\) for \(\varphi \) w.r.t. \(\{F_n\}\) and \(h_{\mu }(X)\) is the measure-theoretic entropy.

     

Steady Nearly Incompressible Vector Fields in Two-Dimension: Chain Rule and Renormalization

Given bounded vector field \({b : {\mathbb{R}^{d}} \to {\mathbb{R}^{d}}}\), scalar field \({u : {\mathbb{R}^{d}} \to {\mathbb{R}}}\), and a smooth function \({\beta : {\mathbb{R}} \to {\mathbb{R}}}\), we study the characterization of the distribution \({{\rm div}(\beta(u)b)}\) in terms of div b and div(ub). In the case of BV vector fields b (and under some further assumptions), such characterization was obtained by L. Ambrosio, C. De Lellis and J. Malý, up to an error term which is a measure concentrated on the so-called tangential set of b. We answer some questions posed in their paper concerning the properties of this term. In particular, we construct a nearly incompressible BV vector field b and a bounded function u for which this term is nonzero. For steady nearly incompressible vector fields b (and under some further assumptions), in the case when d = 2, we provide complete characterization of div(\({\beta(u)b}\)) in terms of div b and div(ub). Our approach relies on the structure of level sets of Lipschitz functions on \({{\mathbb{R}^{2}}}\) obtained by G. Alberti, S. Bianchini and G. Crippa. Extending our technique, we obtain new sufficient conditions when any bounded weak solution u of \({\partial_t u + b \cdot \nabla u=0}\) is renormalized, that is when it also solves \({\partial_t \beta(u) + b \cdot \nabla \beta(u)=0}\) for any smooth function \({\beta \colon{\mathbb{R}} \to {\mathbb{R}}}\). As a consequence, we obtain new a uniqueness result for this equation.     

Local and Global Behaviour of Solutions to Nonlinear Equations with Natural Growth Terms

In this paper we study the Dirichlet problem

$\left\{\begin{array}{lll}-\Delta_p{u} = \sigma |u|^{p-2}u + \omega \quad {\rm in}\;\Omega,\\ u = 0 \qquad\quad\qquad\quad\;\qquad{\rm on}\;\partial\Omega,\end{array}\right.$

, where σ and ω are nonnegative Borel measures, and \({\Delta_p{u} = \nabla \cdot (\nabla{u} \, |\nabla{u}|^{p-2})}\) is the p-Laplacian. Here \({\Omega \subseteq \mathbf{R}^n}\) is either a bounded domain, or the entire space. Our main estimates concern optimal pointwise bounds of solutions in terms of two local Wolff’s potentials, under minimal regularity assumed on σ and ω. In addition, analogous results for equations modeled by the k-Hessian in place of the p-Laplacian will be discussed.

     

Influence of hydrodynamic slip on convective transport in flow past a circular cylinder

The presence of a finite tangential velocity on a hydrodynamically slipping surface is known to reduce vorticity production in bluff body flows substantially while at the same time enhancing its convection downstream and into the wake. Here, we investigate the effect of hydrodynamic slippage on the convective heat transfer (scalar transport) from a heated isothermal circular cylinder placed in a uniform cross-flow of an incompressible fluid through analytical and simulation techniques. At low Reynolds (\({\textit{Re}}\ll 1\)) and high Péclet (\({\textit{Pe}}\gg 1\)) numbers, our theoretical analysis based on Oseen and thermal boundary layer equations allows for an explicit determination of the dependence of the thermal transport on the non-dimensional slip length \(l_s\). In this case, the surface-averaged Nusselt number, Nu transitions gradually between the asymptotic limits of \(Nu \sim {\textit{Pe}}^{1/3}\) and \(Nu \sim {\textit{Pe}}^{1/2}\) for no-slip (\(l_s \rightarrow 0\)) and shear-free (\(l_s \rightarrow \infty \)) boundaries, respectively. Boundary layer analysis also shows that the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) holds for a shear-free cylinder surface in the asymptotic limit of \({\textit{Re}}\gg 1\) so that the corresponding heat transfer rate becomes independent of the fluid viscosity. At finite \({\textit{Re}}\), results from our two-dimensional simulations confirm the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) for a shear-free boundary over the range \(0.1 \le {\textit{Re}}\le 10^3\) and \(0.1\le {\textit{Pr}}\le 10\). A gradual transition from the lower asymptotic limit corresponding to a no-slip surface, to the upper limit for a shear-free boundary, with \(l_s\), is observed in both the maximum slip velocity and the Nu. The local time-averaged Nusselt number \(Nu_{\theta }\) for a shear-free surface exceeds the one for a no-slip surface all along the cylinder boundary except over the downstream portion where unsteady separation and flow reversal lead to an appreciable rise in the local heat transfer rates, especially at high \({\textit{Re}}\) and Pr. At a Reynolds number of \(10^3\), the formation of secondary recirculating eddy pairs results in appearance of additional local maxima in \(Nu_{\theta }\) at locations that are in close proximity to the mean secondary stagnation points. As a consequence, Nu exhibits a non-monotonic variation with \(l_s\) increasing initially from its lowermost value for a no-slip surface and then decreasing before rising gradually toward the upper asymptotic limit for a shear-free cylinder. A non-monotonic dependence of the spanwise-averaged Nu on \(l_s\) is observed in three dimensions as well with the three-dimensional wake instabilities that appear at sufficiently low \(l_s\), strongly influencing the convective thermal transport from the cylinder. The analogy between heat transfer and single-component mass transfer implies that our results can directly be applied to determine the dependency of convective mass transfer of a single solute on hydrodynamic slip length in similar configurations through straightforward replacement of Nu and \({\textit{Pr}}\) with Sherwood and Schmidt numbers, respectively.     

Exponential Decay of the Vorticity in the Steady-State Flow of a Viscous Liquid Past a Rotating Body

We construct a Sobolev homeomorphism in dimension \({n \geqq 4,\,f \in W^{1,1}((0, 1)^n,\mathbb{R}^n)}\) such that \({J_f = {\rm det} Df > 0}\) on a set of positive measure and J _f  < 0 on a set of positive measure. It follows that there are no diffeomorphisms (or piecewise affine homeomorphisms) f _k such that \({f_k\to f}\) in \({W^{1,1}_{\rm loc}}\).     

Asymptotics of the Solutions of the Random Schrödinger Equation

We consider solutions of the Schrödinger equation with a weak time-dependent random potential. It is shown that when the two-point correlation function of the potential is rapidly decaying, then the Fourier transform \({\hat\zeta_\epsilon(t,\xi)}\) of the appropriately scaled solution converges point-wise in ξ to a stochastic complex Gaussian limit. On the other hand, when the two-point correlation function decays slowly, we show that the limit of \({\hat\zeta_\epsilon(t,\xi)}\) has the form \({\hat\zeta_0(\xi){\rm exp}(iB_\kappa(t,\xi))}\) where B _κ (t, ξ) is a fractional Brownian motion.     

Physical Clogging of Uniformly Graded Porous Media Under Constant Flow rates

This study investigated the physical clogging of uniformly graded porous media under constant flow rates using natural porous media and suspensions. The porous media selected for this experimental study was a fine-to-medium sandy soil fractioned into thirteen uniformly graded beds: seven unisize beds and six uniform beds. The physical clogging of the beds was studied using two types of silt suspensions as along with two suspension concentrations and three water discharges. It was found that the permeability reduction due to physical clogging \([(K_\mathrm{i} - K_\mathrm{t})/K_\mathrm{i}]\) increased with decreasing \({D}_{15}/{d}_{85}\) ratios until a critical value of \({D}_{15}/{d}_{85}\), after which a surface mat of suspension was formed on the porous media. It was also found that the value of reduced permeability at any time (at any number of pore volumes of injected suspension-laden water), \(K_\mathrm{t}\), is directly proportional to square of \({D}_{15}\) and inversely proportional to \({C}_{\mathrm{u}}\) of the porous media and \({d}_{85}\) of suspensions. The effects of suspension type and flow rates on physical clogging seemed to depend on the size of the pores in the porous media.     

The Viscosity Method for the Homogenization of Soft Inclusions

In this paper, we consider periodic soft inclusions T _ε with periodicity ε, where the solution, u _ε , satisfies semi-linear elliptic equations of non-divergence in \({\Omega_{\epsilon}=\Omega\setminus \overline{T}_\epsilon}\) with Neumann data on \({\partial T^{\mathfrak a}}\). The difficulty lies in the non-divergence structure of the operator where the standard energy method, which is based on the divergence theorem, cannot be applied. The main object is to develop a viscosity method to find the homogenized equation satisfied by the limit of u _ε , referred to as u, as ε approaches to zero. We introduce the concept of a compatibility condition between the equation and the Neumann condition on the boundary for the existence of uniformly bounded periodic first correctors. The concept of a second corrector is then developed to show that the limit, u, is the viscosity solution of a homogenized equation.     

Mixing and Bimolecular Reaction Kinetics in a Plane Poisseulle Flow

Mixing and a nonlinear bimolecular chemical reaction (reactant A + reactant B → product; reaction rate r?=?κc ₁ c ₂) in laminar shear flow are investigated. It is found that asymptotically the dominant balance between the rates of production and dissipation of the mean-squared concentration fluctuations \((\sigma_{c_1 }^2 ,\sigma_{c_2 }^2)\) and cross-covariance of concentration fluctuations \((\overline {c_1 c_2 })\) occurs under nonreactive and reactive conditions. Longitudinal dispersion of the cross-sectional averages (C ₁, C ₂), and variances and the cross-covariance of reactant concentrations can be asymptotically quantified by the classic Taylor dispersion coefficient (D) even under reactive conditions. The characteristic time-scale (τ) over which molecular diffusion dissipates concentration variance and the cross-covariance of reactant concentrations is also shown to be the same under nonreactive and reactive conditions. A variational estimate of τ is shown to be close to the values inferred from detailed numerical simulation. The production-dissipation balance implies that the cross-sectional averaged reaction rate follows \(\overline r =\kappa_{eff} C_1 C_2 \) and \(\kappa _{eff} \approx \kappa \left[ {1+2D\tau \left( {{\partial \ln C_1 } \mathord{\left/ {\vphantom {{\partial \ln C_1 } {\partial x}}} \right. \kern-\nulldelimiterspace} {\partial x}} \right)\left( {{\partial \ln C_2 } \mathord{\left/ {\vphantom {{\partial \ln C_2 } {\partial x}}} \right. \kern-\nulldelimiterspace} {\partial x}} \right)} \right]\). The effective reaction rate parameter (κ _eff ) is higher than that of well-mixed batch test reaction rate constant (κ) for initially overlapping species and κ _eff is smaller than κ for initially non-overlapping species.     

Impermeability Through a Perforated Domain for the Incompressible two dimensional Euler Equations

We study the asymptotic behavior of the motion of an ideal incompressible fluid in a perforated domain. The porous medium is composed of inclusions of size \({\varepsilon}\) separated by distances \({d_{\varepsilon}}\) and the fluid fills the exterior. If the inclusions are distributed on the unit square, the asymptotic behavior depends on the limit of \({\frac{d_{\varepsilon}}\varepsilon}\) when \({\varepsilon}\) goes to zero. If \({\frac{d_{\varepsilon}}\varepsilon \to \infty}\), then the limit motion is not perturbed by the porous medium, namely, we recover the Euler solution in the whole space. If, on the contrary, \({\frac{d_{\varepsilon}}\varepsilon \to 0}\), then the fluid cannot penetrate the porous region, namely, the limit velocity verifies the Euler equations in the exterior of an impermeable square. If the inclusions are distributed on the unit segment then the behavior depends on the geometry of the inclusion: it is determined by the limit of \({\frac{d_{\varepsilon}}{\varepsilon^{2+\frac1\gamma}}}\) where \({\gamma \in (0,\infty]}\) is related to the geometry of the lateral boundaries of the obstacles. If \({\frac{d_{\varepsilon}}{\varepsilon^{2+\frac1\gamma}} \to \infty}\), then the presence of holes is not felt at the limit, whereas an impermeable wall appears if this limit is zero. Therefore, for a distribution in one direction, the critical distance depends on the shape of the inclusions; in particular, it is equal to \({\varepsilon^{3}}\) for balls.     

The Mixed Convection Boundary Layer on a Vertical Melting Front in a Non-Darcian Porous Medium

The effect of surface melting on the dual solutions that can arise in the problem of the mixed convection boundary-layer flow past a vertical surface embedded in a non-Darcian porous medium is considered. The problem is described by M, melting parameter, \(\lambda \), mixed convection parameter, and \(\gamma \), the flow inertia coefficient, numerical results being obtained in terms of these three parameters. It is seen that the melting phenomenon reduces the heat transfer rate and enhances the boundary-layer separation at the solid–liquid interface. Asymptotic solutions for the forced convection, \(\lambda =0\), and free convection, large \(\lambda \), limits are derived.     

Focusing Quantum Many-body Dynamics: The Rigorous Derivation of the 1D Focusing Cubic Nonlinear Schrödinger Equation

We consider the dynamics of N bosons in 1D. We assume that the pair interaction is attractive and given by \({N^{\beta-1}V(N^{\beta}.) where }\) where \({\int V \leqslant 0}\). We develop new techniques in treating the N-body Hamiltonian so that we overcome the difficulties generated by the attractive interaction and establish new energy estimates. We also prove the optimal 1D collapsing estimate which reduces the regularity requirement in the uniqueness argument by half a derivative. We derive rigorously the 1D focusing cubic NLS with a quadratic trap as the \({N \rightarrow \infty}\) limit of the N-body dynamic and hence justify the mean-field limit and prove the propagation of chaos for the focusing quantum many-body system.     

Glimpses of Flow Development and Degradation During Type B Drag Reduction by Aqueous Solutions of Polyacrylamide B1120

Flow development and degradation during Type B turbulent drag reduction by 0.10 to 10 wppm solutions of a partially-hydrolysed polyacrylamide B1120 of MW \(=\) 18x10⁶ was studied in a smooth pipe of ID \(=\) 4.60 mm and L/D \(=\) 210 at Reynolds numbers from 10000 to 80000 and wall shear stresses Tw from 8 to 600 Pa. B1120 solutions exhibited facets of a Type B ladder, including segments roughly parallel to, but displaced upward from, the P-K line; those that attained asymptotic maximum drag reduction at low Re√ f but departed downwards into the polymeric regime at a higher retro-onset Re√ f; and segments at MDR for all Re√ f. Axial flow enhancement profiles of S\(^{\prime }\) vs L/D reflected a superposition of flow development and polymer degradation effects, the former increasing and the latter diminishing S\(^{\prime }\) with increasing distance downstream. Solutions that induced normalized flow enhancements S\(^{\prime }\)/S\(^{\prime }_{\mathrm {m}} <\) 0.4 developed akin to solvent, with L_e,p/D \(=\) L_e,n/D \(<\) 42.3, while those at maximum drag reduction showed entrance lengths L_e,m/D \(\sim \) 117, roughly 3 times the solvent L_e,n/D. Degradation kinetics were inferred by first detecting a falloff point (Re√f^∧, S\(^{{\prime }\wedge }\)), of maximum observed flow enhancement, for each polymer solution. A plot of S\(^{{\prime }\wedge }\)vs C revealed S\(^{{\prime }\wedge }\)linear in C at low C, with lower bound [S\(^{\prime }\)] \(=\) 5.0 wppm^??1, and S\(^{{\prime }\wedge }\) independent of C at high C, with upper bound S\(^{\prime }_{\mathrm {m}} =\) 15.9. The ratio S\(^{\prime }\)/S\(^{{\prime }\wedge }\) in any pipe section was interpreted to be the undegraded fraction of original polymer therein. Semi-log plots of (S\(^{\prime }\)/S\(^{{\prime }\wedge }\)) at a section vs transit time from pipe entrance thereto revealed first order kinetics, from which apparent degradation rate constants k_deg s^??1 and entrance severities ?ln(S\(^{\prime }\)/S\(^{{\prime }\wedge }\))₀ were extracted. At constant C, k_deg increased linearly with increasing wall shear stress Tw, and at constant Tw, k_deg was independent of C, providing a B1120 degradation modulus (k_deg/Tw) \(=\) (0.012 \(\pm \) 0.001) (Pa s)^??1 for 8 \(<\) Tw Pa \(<\) 600, 0.30 \(<\) C wppm \(<\) 10. Entrance severities were negligible below a threshold Twe \(\sim \) 30 Pa and increased linearly with increasing Tw for Tw \(>\) Twe. The foregoing methods were applied to Type A drag reduction by 0.10 to 10 wppm solutions of a polyethyleneoxide PEO P309, MW \(=\) 11x10⁶, in a smooth pipe of ID \(=\) 7.77 mm and L/D \(=\) 220 at Re from 4000 to 115000. P309 solutions that induced S\(^{\prime }\)/S\(^{\prime }_{\mathrm {m}} <\) 0.4 developed akin to solvent, with L_e,p/D \(=\) L_e,n/D \(<\) 23, while those at MDR had entrance lengths L_e,m/D \(\sim \) 93, roughly 4 times the solvent L_e,n/D. P309 solutions described a Type A fan distorted by polymer degradation. A typical trajectory departed the P-K line at an onset point Re√ f* followed by ascending and descending polymeric regime segments separated by a falloff point Re√f^∧, of maximum flow enhancement; for all P309 solutions, onset Re√ f* = 550 \(\pm \) 100 and falloff Re√f^∧ = 2550 \(\pm \) 250, the interval between them delineating Type A drag reduction unaffected by degradation. A plot of falloff S\(^{{\prime }\wedge }\) vs C for PEO P309 solutions bore a striking resemblance to the analogous S\(^{{\prime }\wedge }\) vs C plot for solutions of PAMH B1120, indicating that the initial Type A drag reduction by P309 after onset at Re√ f* had evolved to Type B drag reduction by falloff at Re√f^∧. Presuming that Type B behaviour persisted past falloff permitted inference of P309 degradation kinetics; k_deg was found to increase linearly with increasing Tw at constant C and was independent of C at constant Tw, providing a P309 degradation modulus (k_deg/Tw) \(=\) (0.011 \(\pm \) 0.002) (Pa s)^??1 for 4 \(<\) Tw Pa \(<\) 400, 0.10 \(<\) C wppm < 5.0. Comparisons between the present degradation kinetics and previous literature showed (k_deg/Tw) data from laboratory pipes of D \(\sim \) 0.01 m to lie on a simple extension of (k_deg/Tw) data from pipelines of D \(\sim \) 0.1 m and 1.0 m, along a power-law relation (k_deg/Tw) \(=\) 10^??5.4.D^??1.6. Intrinsic slips derived from PAMH B1120 and PEO P309-at-falloff experiments were compared with previous examples from Type B drag reduction by polymers with vinylic and glycosidic backbones, showing: (i) For a given polymer, [S\(^{\prime }\)] was independent of Re√ f and pipe ID, implying insensitivity to both micro- and macro-scales of turbulence; and (ii) [S\(^{\prime }\)] increased linearly with increasing polymer chain contour length Lc, the proportionality constant \(\beta =\) 0.053 \(\pm \) 0.036 enabling estimation of flow enhancement S\(^{\prime } =\) C.Lc.β for all Type B drag reduction by polymers.     

From the Highly Compressible Navier–Stokes Equations to Fast Diffusion and Porous Media Equations,Existence of Global Weak Solution for the Quasi-Solutions

We consider the compressible Navier–Stokes equations for viscous and barotropic fluids with density dependent viscosity. The aim is to investigate mathematical properties of solutions of the Navier–Stokes equations using solutions of the pressureless Navier–Stokes equations, that we call quasi solutions. This regime corresponds to the limit of highly compressible flows. In this paper we are interested in proving the announced result in Haspot (Proceedings of the 14th international conference on hyperbolic problems held in Padova, pp 667–674, 2014) concerning the existence of global weak solution for the quasi-solutions, we also observe that for some choice of initial data (irrotationnal) the quasi solutions verify the porous media, the heat equation or the fast diffusion equations in function of the structure of the viscosity coefficients. In particular it implies that it exists classical quasi-solutions in the sense that they are \({C^{\infty}}\) on \({(0,T)\times \mathbb{R}^{N}}\) for any \({T > 0}\). Finally we show the convergence of the global weak solution of compressible Navier–Stokes equations to the quasi solutions in the case of a vanishing pressure limit process. In particular for highly compressible equations the speed of propagation of the density is quasi finite when the viscosity corresponds to \({\mu(\rho)=\rho^{\alpha}}\) with \({\alpha > 1}\). Furthermore the density is not far from converging asymptotically in time to the Barrenblatt solution of mass the initial density \({\rho_{0}}\).     

Statistical Analysis of Scalar Dissipation Rate Transport in Turbulent Partially Premixed Flames: A Direct Numerical Simulation Study

Statistically planar turbulent premixed and partially premixed flames for different initial turbulence intensities are simulated for global equivalence ratios ??>?=?0.7 and ??>?=?1.0 using three-dimensional Direct Numerical Simulations (DNS) with simplified chemistry. For the simulations of partially premixed flames, a random distribution of equivalence ratio following a bimodal distribution of equivalence ratio is introduced in the unburned reactants ahead of the flame. The simulation parameters in all of the cases were chosen such that the combustion situation belongs to the thin reaction zones regime. The DNS data has been used to analyse the behaviour of the dissipation rate transports of both active and passive scalars (i.e. the fuel mass fraction Y _F and the mixture fraction ξ) in the context of Reynolds Averaged Navier–Stokes (RANS) simulations. The behaviours of the unclosed terms of the Favre averaged scalar dissipation rates of fuel mass fraction and mixture fraction (i.e. \(\widetilde {\varepsilon }_Y =\overline {\rho D\nabla Y_F^{\prime \prime } \cdot \nabla Y_F^{\prime \prime } } /\overline{\rho }\) and \(\widetilde {\varepsilon }_\xi =\overline {\rho D\nabla \xi ^{\prime \prime }\cdot \nabla \xi ^{\prime \prime }} /\overline {\rho })\) transport equations have been analysed in detail. In the case of the \(\widetilde {\varepsilon }_Y \) transport, it has been observed that the turbulent transport term of scalar dissipation rate remains small throughout the flame brush whereas the terms due to density variation, scalar–turbulence interaction, reaction rate and molecular dissipation remain the leading order contributors. The term arising due to density variation remains positive throughout the flame brush and the combined contribution of the reaction and molecular dissipation to the \(\widetilde {\varepsilon }_Y \) transport remains negative throughout the flame brush in all cases. However, the behaviour of scalar–turbulence interaction term of the \(\widetilde {\varepsilon }_Y \) transport equation is significantly affected by the relative strengths of turbulent straining and the straining due to chemical heat release. In the case of the \(\widetilde {\varepsilon }_\xi \) transport, the turbulent transport term remains small throughout the flame brush and the density variation term is found to be negligible in all cases, whilst the reaction rate term is exactly zero. The scalar–turbulence interaction term and molecular dissipation term remain the leading order contributors to the \(\widetilde {\varepsilon }_\xi \) transport throughout the flame brush in all cases that have been analysed in the present study. Performances of existing models for the unclosed terms of the transport equations of \(\widetilde {\varepsilon }_Y \) and \(\widetilde {\varepsilon }_\xi \) are assessed with respect to the corresponding quantities obtained from DNS data. Based on this exercise either suitable models have been identified or new models have been proposed for the accurate closure of the unclosed terms of both \(\widetilde {\varepsilon }_Y \) and \(\widetilde {\varepsilon }_\xi \) transport equations in the context of Reynolds Averaged Navier–Stokes (RANS) simulations.     

Approximation of Flat <Emphasis Type="Italic">W</Emphasis><Superscript>2,2</Superscript> Isometric Immersions by Smooth Ones

Let \({S\subset\mathbb{R}^2}\) be a bounded Lipschitz domain and denote by \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\) the set of mappings \({u\in W^{2,2}(S;\mathbb{R}^3)}\) which satisfy \({(\nabla u)^T(\nabla u) = Id}\) almost everywhere. Under an additional regularity condition on the boundary \({\partial S}\) (which is satisfied if \({\partial S}\) is piecewise continuously differentiable), we prove that the strong W ^2,2 closure of \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)\cap C^{\infty}(\overline{S};\mathbb{R}^3)}\) agrees with \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\).