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.     

Desingularization of Vortex Rings and Shallow Water Vortices by a Semilinear Elliptic Problem

Steady vortices for the three-dimensional Euler equation for inviscid incompressible flows and for the shallow water equation are constructed and shown to tend asymptotically to singular vortex filaments. The construction is based on a study of solutions to the semilinear elliptic problem $$ \left\{ \begin{aligned} -{\rm div} \left(\frac{\nabla u_{\varepsilon}}{b}\right) & = \frac{1}{\varepsilon^2} b f \left(u_{\varepsilon} - \log \tfrac{1}{\varepsilon} q \right) & & \text{ in } \; \Omega, \\u_\varepsilon & = 0 & & \text{ on } \; \partial \Omega, \end{aligned}\right.$$ for small values of ${\varepsilon > 0}$ .     

Determination of turbulent pressure loss of non-newtonian oil flow in rough pipes

Summary A detailed theoretical analysis is carried out of the turbulent flow of non-newtonian fluids in rough pipes in which approximate assumptions are of negligible importance. A new friction-factor equation (BNS-equation) is obtained which is valid for turbulent flow of power-law type non-newtonian fluids in the transition region between smooth and wholly rough wall turbulence.The accuracy of the deduced formula was checked by in situ measurements on a crude oil pipeline having a diameter of 305 mm and a length of 161 km. For this pipeline the Reynolds numbers ranged between 10⁴ and 10⁵. The computed and the measured friction factors showed good agreement. The mean of the absolute values of the relative error was 4.1%, while the standard deviation was 0.81%, these results appear to give better fitting than other similar equations published earlier.The BNS-equation is the first analytically deduced relationship for non-newtonian crude oils which includes the effects of both Reynolds number and relative roughness of the pipe.With 6 figures and 2 tables     

Miniaturized Compression Test at Very High Strain Rates by Direct Impact

A modified miniaturized version of the Direct Impact Compression Test (DICT) technique is described in this paper. The method permits determination of the rate-sensitive plastic properties of materials up to strain rate ∼10⁵ s⁻¹. Miniaturization of the experimental setup with specimen dimensions: diameter d _S = 2.0 mm and thickness l _S = 1.0 mm, Hopkinson bar diameter 5.2 mm, with application of a novel optical arrangement in measurement of specimen strain, makes possible compression tests at strain rates from ∼10³ s⁻¹ to ∼10⁵ s⁻¹. In order to estimate the rate sensitivity of a low-alloy construction steel, quasi-static, Split Hopkinson Pressure Bar (SHPB) and DICT tests have been performed at room temperature within the rate spectrum ranging from 5*10⁻⁴ s⁻¹ to 5*10⁴ s⁻¹. Adiabatic heating and friction effects are analyzed and the final true stress versus true strain curves at different strain rates are corrected to a constant temperature and zero friction. The results have been analyzed in the form of true stress versus the logarithm of strain rate and they show two regions of a constant rate sensitivity : relatively low up to the strain rate threshold ∼50 s⁻¹, and relatively high above the threshold, up to strain rate ∼4.5*10⁴ s⁻¹.     

Extension of k–ε models for a turbulent falling film

An extension applicable to various k–ε models is proposed to account for the damping of turbulence due to the surface tension. It consists of a sink in the equation for $ {{\overline{D} k} \mathord{\left/ {\vphantom {{\overline{D} k} {\overline{D} t}}} \right. \kern-0em} {\overline{D} t}} $ and a source in the equation for $ {{\overline{D} \varepsilon } \mathord{\left/ {\vphantom {{\overline{D} \varepsilon } {\overline{D} t}}} \right. \kern-0em} {\overline{D} t}} $ both derived from dimensional analysis. First numerical experiments are undertaken with a commercial CFD software. Reasons are given why the tuning of the model should be performed with thermal experiments. The proposed model is practical, since it only needs the programming of source terms. On account of it’s mathematical structure it needs a very small closure coefficient. The intention of this article is to stimulate the numerical turbulence research in a way that is applicable to the engineering practice. A comparison to experimental data is not done here.     

A-Priori Direct Numerical Simulation Modelling of Co-variance Transport in Turbulent Stratified Flames

Three-dimensional Direct Numerical Simulations of statistically planar turbulent stratified flames at global equivalence ratios <???>?=?0.7 and <???>?=?1.0 have been carried out to analyse the statistical behaviour of the transport of co-variance of the fuel mass fraction Y _F and mixture fraction ξ (i.e. $\widetilde{Y_F^{\prime\prime} \xi ^{\prime\prime}}={\overline {\rho Y_F^{\prime\prime} \xi^{\prime\prime}} } \Big/ {\overline \rho })$ for Reynolds Averaged Navier Stokes simulations where $\overline q $ , $\tilde{q} ={\overline {\rho q} } \big/ {\overline \rho }$ and $q^{\prime\prime}= q-\tilde{q}$ are Reynolds averaged, Favre mean and Favre fluctuation of a general quantity q with ρ being the gas density and the overbar suggesting a Reynolds averaging operation. It has been found that existing algebraic expressions may not capture the statistical behaviour of $\widetilde{Y_F^{\prime\prime} \xi^{\prime\prime}}$ with sufficient accuracy in low Damköhler number combustion and therefore, a transport equation for $\widetilde{Y_F^{\prime\prime} \xi^{\prime\prime}}$ may need to be solved. The statistical behaviours of $\widetilde{Y_F^{\prime\prime} \xi^{\prime\prime}}$ and the unclosed terms of its transport equation (i.e. the terms originating from turbulent transport T ₁ , reaction rate T ₄ and molecular dissipation $\left( {-D_2 } \right))$ have been analysed in detail. The contribution of T ₁ remains important for all cases considered here. The term T ₄ acts as a major contributor in <???>?=?1.0 cases, but plays a relatively less important role in <???>?=?0.7 cases, whereas the term $\left( {-D_2 } \right)$ acts mostly as a leading order sink. Through an a-priori DNS analysis, the performances of the models for T ₁ , T ₄ and $\left( {-D_2 } \right)$ have been addressed in detail. A model has been identified for the turbulent transport term T ₁ which satisfactorily predicts the corresponding term obtained from DNS data. The models for T ₄ , which were originally proposed for high Damköhler number flames, have been modified for low Damköhler combustion. Predictions of the modified models are found to be in good agreement with T ₄ obtained from DNS data. It has been found that existing algebraic models for $D_2 =2\overline {\rho D\nabla Y_F^{\prime\prime} \nabla \xi^{\prime\prime}} $ (where D is the mass diffusivity) are not sufficient for low Damköhler number combustion and therefore, a transport equation may need to be solved for the cross-scalar dissipation rate $\widetilde{\varepsilon }_{Y\xi } ={\overline {\rho D\nabla Y_F^{\prime\prime} \nabla \xi^{\prime\prime}} } \big/ {\overline \rho }$ for the closure of the $\widetilde{Y_F^{\prime\prime} \xi^{\prime\prime}}$ transport equation.     

Hadamard Variational Formula for the Green’s Function of the Boundary Value Problem on the Stokes Equations

For every ${\varepsilon > 0}$ , we consider the Green’s matrix ${G_{\varepsilon}(x, y)}$ of the Stokes equations describing the motion of incompressible fluids in a bounded domain ${\Omega_{\varepsilon} \subset \mathbb{R}^d}$ , which is a family of perturbation of domains from ${\Omega\equiv \Omega_0}$ with the smooth boundary ${\partial\Omega}$ . Assuming the volume preserving property, that is, ${\mbox{vol.}\Omega_{\varepsilon} = \mbox{vol.}\Omega}$ for all ${\varepsilon > 0}$ , we give an explicit representation formula for ${\delta G(x, y) \equiv \lim_{\varepsilon\to +0}\varepsilon^{-1}(G_{\varepsilon}(x, y) - G_0(x, y))}$ in terms of the boundary integral on ${\partial \Omega}$ of ${G_0(x, y)}$ . Our result may be regarded as a classical Hadamard variational formula for the Green’s functions of the elliptic boundary value problems.     

A comparison of time integration methods in an unsteady low‐Reynolds‐number flow

This paper describes three different time integration methods for unsteady incompressible Navier–Stokes equations. Explicit Euler and fractional‐step Adams–Bashford methods are compared with an implicit three‐level method based on a steady‐state SIMPLE method. The implicit solver employs a dual time stepping and an iteration within the time step. The spatial discretization is based on a co‐located finite‐volume technique. The influence of the convergence limits and the time‐step size on the accuracy of the predictions are studied. The efficiency of the different solvers is compared in a vortex‐shedding flow over a cylinder in the Reynolds number range of 100–1600. A high‐Reynolds‐number flow over a biconvex airfoil profile is also computed. The computations are performed in two dimensions. At the low‐Reynolds‐number range the explicit methods appear to be faster by a factor from 5 to 10. In the high‐Reynolds‐number case, the explicit Adams–Bashford method and the implicit method appear to be approximately equally fast while yielding similar results. Copyright © 2002 John Wiley & Sons, Ltd.     

Singular perturbations of two-point boundary problems for systems of ordinary differential equations

Asymptotic solutions of linear systems of ordinary differential equations are employed to discuss the relationship of the solution of a certain “complete” boundary problem.

$$\begin{gathered} \left\{ \begin{gathered} {\text{ }}\frac{{d{\text{ }}x_1 }}{{d{\text{ }}t}} = A_{11} (t,\varepsilon ){\text{ }}x_1 (t,\varepsilon ){\text{ }} + \cdots + A_{1p} (t,\varepsilon ){\text{ }}x_p (t,\varepsilon ) \hfill \\ \varepsilon ^{h_2 } \frac{{d{\text{ }}x_2 }}{{d{\text{ }}t}} = A_{21} (t,\varepsilon ){\text{ }}x_1 (t,\varepsilon ){\text{ }} + \cdots + A_{2p} (t,\varepsilon ){\text{ }}x_p (t,\varepsilon ) \hfill \\ {\text{ }} \vdots {\text{ }} \vdots {\text{ }} \vdots \hfill \\ \varepsilon ^{h_p } \frac{{d{\text{ }}x_2 }}{{d{\text{ }}t}} = A_{p1} (t,\varepsilon ){\text{ }}x_1 (t,\varepsilon ){\text{ }} + \cdots + A_{pp} (t,\varepsilon ){\text{ }}x_p (t,\varepsilon ) \hfill \\ \end{gathered} \right\} \hfill \\ {\text{ }}R(\varepsilon ){\text{ }}x(a,{\text{ }}\varepsilon ){\text{ }} + {\text{ }}S(\varepsilon ){\text{ }}x(b,{\text{ }}\varepsilon ) = c(\varepsilon ){\text{ }} \hfill \\ \end{gathered}$$     

Asymptotic analysis of linearly elastic shells. III. Justification of Koiter's shell equations

We consider as in Parts I and II a family of linearly elastic shells of thickness 2?, all having the same middle surfaceS=?(?)?R ³, whereω?R ² is a bounded and connected open set with a Lipschitz-continuous boundary, and? ∈ ?³ (?;R ³). The shells are clamped on a portion of their lateral face, whose middle line is?(γ ₀), whereγ ₀ is a portion of?ω withlength γ ₀>0. For all?>0, let $\zeta _i^\varepsilon$ denote the covariant components of the displacement $u_i^\varepsilon g^{i,\varepsilon }$ of the points of the shell, obtained by solving the three-dimensional problem; let $\zeta _i^\varepsilon$ denote the covariant components of the displacement $\zeta _i^\varepsilon$ a ⁱ of the points of the middle surfaceS, obtained by solving the two-dimensional model ofW.T. Koiter, which consists in finding $$\zeta ^\varepsilon = \left( {\zeta _i^\varepsilon } \right) \in V_K (\omega ) = \left\{ {\eta = (\eta _\iota ) \in {\rm H}^1 (\omega ) \times H^1 (\omega ) \times H^2 (\omega ); \eta _i = \partial _v \eta _3 = 0 on \gamma _0 } \right\}$$ such that $$\begin{gathered} \varepsilon \mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \gamma _{\sigma \tau } (\zeta ^\varepsilon )\gamma _{\alpha \beta } (\eta )\sqrt a dy + \frac{{\varepsilon ^3 }}{3} \mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \rho _{\sigma \tau } (\zeta ^\varepsilon )\rho _{\alpha \beta } (\eta )\sqrt a dy \hfill \\ = \mathop \smallint \limits_\omega p^{i,\varepsilon } \eta _i \sqrt a dy for all \eta = (\eta _i ) \in V_K (\omega ), \hfill \\ \end{gathered}$$ where $a^{\alpha \beta \sigma \tau }$ are the components of the two-dimensional elasticity tensor ofS, $\gamma _{\alpha \beta }$ (η) and $\rho _{\alpha \beta }$ (η) are the components of the linearized change of metric and change of curvature tensors ofS, and $p^{i,\varepsilon }$ are the components of the resultant of the applied forces. Under the same assumptions as in Part I, we show that the fields $\frac{1}{{2_\varepsilon }}\smallint _{ - \varepsilon }^\varepsilon u_i^\varepsilon g^{i,\varepsilon } dx_3^\varepsilon$ and $\zeta _i^\varepsilon$ a ⁱ , both defined on the surfaceS, have the same principal part as? → 0, inH ¹ (ω) for the tangential components, and inL ²(ω) for the normal component; under the same assumptions as in Part II, we show that the same fields again have the same principal part as? → 0, inH ¹ (ω) for all their components. For “membrane” and “flexural” shells, the two-dimensional model ofW.T. Koiter is therefore justified.     

On the efficiency of semi‐implicit and semi‐Lagrangian spectral methods for the calculation of incompressible flows

Classical semi‐implicit backward Euler/Adams–Bashforth time discretizations of the Navier–Stokes equations induce, for high‐Reynolds number flows, severe restrictions on the time step. Such restrictions can be relaxed by using semi‐Lagrangian schemes essentially based on splitting the full problem into an explicit transport step and an implicit diffusion step. In comparison with the standard characteristics method, the semi‐Lagrangian method has the advantage of being much less CPU time consuming where spectral methods are concerned. This paper is devoted to the comparison of the ‘semi‐implicit’ and ‘semi‐Lagrangian’ approaches, in terms of stability, accuracy and computational efficiency. Numerical results on the advection equation, Burger's equation and finally two‐ and three‐dimensional Navier–Stokes equations, using spectral elements or a collocation method, are provided. Copyright © 2001 John Wiley & Sons, Ltd.     

Existence and Concentration Result for the Kirchhoff Type Equations with General Nonlinearities

In this paper we study the existence and concentration behaviors of positive solutions to the Kirchhoff type equations $$- \varepsilon^2 M \left(\varepsilon^{2-N}\!\!\int_{\mathbf{R}^N}|\nabla u|^2\,\mathrm{d} x \right) \Delta u \!+\! V(x) u \!=\! f(u) \quad{\rm in}\ \mathbf{R}^N, \quad u \!\in\! H^1(\mathbf{R}^N), \ N \!\geqq\!1,$$ where M and V are continuous functions. Under suitable conditions on M and general conditions on f, we construct a family of positive solutions \({(u_\varepsilon)_{\varepsilon \in (0,\tilde{\varepsilon}]}}\) which concentrates at a local minimum of V after extracting a subsequence (ε _k ).     

A note on the soil-water conductivity of a fractal soil

We show that for a fractal soil the soil-water conductivity, K, is given by $$\frac{K}{{K_\varepsilon }} = (\Theta /\varepsilon )^{2D/3 + 2/(3 - D)}$$ where $K_\varepsilon$ is the saturated conductivity, θ the water content, ? its saturated value and D is the fractal dimension obtained from reinterpreting Millington and Quirk's equation for practical values of the porosity ?, as $$D = 2 + 3\frac{{\varepsilon ^{4/3} + (1 - \varepsilon )^{2/3} - 1}}{{2\varepsilon ^{4/3} \ln ,{\text{ }}\varepsilon ^{ - 1} + (1 - \varepsilon )^{2/3} \ln (1 - \varepsilon )^{ - 1} }}$$ .     

Parallel FEM LES with one-equation subgrid-scale model for incompressible flows

This article develops a parallel large-eddy simulation (LES) with a one-equation subgrid-scale (SGS) model based on the Galerkin finite element method and three-dimensional (3D) brick elements. The governing filtered Navier–Stokes equations were solved by a second-order accurate fractional-step method, which decomposed the implicit velocity–pressure coupling in incompressible flow and segregated the solution to the advection and diffusion terms. The transport equation for the SGS turbulent kinetic energy was solved to calculate the SGS processes. This FEM LES model was applied to study the turbulence of the benchmark open channel flow at a Reynolds number Re_τ = 180 (based on the friction velocity and channel height) using different model constants and grid resolutions. By comparing the turbulence statistics calculated by the current model with those obtained from direct numerical simulation (DNS) and experiments in literature, an optimum set of model constants for the current FEM LES model was established. The budgets of turbulent kinetic energy and vertical Reynolds stress were then analysed for the open channel flow. Finally, the flow structures were visualised to further reveal some important characteristics. It was demonstrated that the current model with the optimum model constants can predict well the organised structure near the wall and free surface, and can be further applied to other fundamental and engineering applications.     

On the logarithmic-law constants and the turbulent boundary layer at low Reynolds numbers

It is shown that a family of formally derived similarity solutions describe to leading order the outer region of a turbulent boundary layer for all Reynolds numbers for which the layer satisfies the logarithmic law-of-the-wall. The family includes Coles' [1] hypothesis. For consistency with this hypothesis and the logarithmic law-of-the-wall, it is further shown that the constants in the latter form the product κC=2+O(ε), suggesting the logarithmic law of the wall be written $${U \mathord{\left/ {\vphantom {U {U_\tau = \kappa ^{ - 1} }}} \right. \kern-\nulldelimiterspace} {U_\tau = \kappa ^{ - 1} }}\ln \left( {e^2 U_\tau {y \mathord{\left/ {\vphantom {y \nu }} \right. \kern-\nulldelimiterspace} \nu }} \right) + O\left( \in \right).$$ A range of data are reprocessed to determine the skin friction coefficientC _f using κC = 2 and these collapse well when plotted against momentum thickness Reynolds number, Re _θ . It is also shown that the form parameter, Π, in Coles hypothesis is not unique but is determined by history effects peculiar to the boundary layer. Expressions are derived forC _f (Re _θ ) and the shape factorH (Re _θ ); both agree closely with the data and are valid over all Reynolds numbers for which the logarithmic law of the wall is satisfied.     

Experimental and numerical study of turbulent flow and heat transfer inside hexagonal duct

Flow and heat transfer characteristics in transition and turbulent regions are studied experimentally and numerically in a horizontal smooth regular hexagonal duct under constant wall temperature boundary condition covering a range of Reynolds number from 2.3 × 10³ to 52 × 10³. Two types of k-omega (standard and shear stress transport (SST)) and three types of k-ε (standard, renormalization (RNG), and realizable) turbulence model are employed for transition and turbulent regions, respectively. Both average and fully developed Darcy friction factor and Nusselt number are presented as a function of Reynolds number. It is seen that k-omega SST and k-ε realizable turbulence models gave the best agreement with the experimental data in transition and turbulent regions, respectively. All the experimental results are correlated within an accuracy of ±13 % and ±7 % for Nusselt number and Darcy friction factor, respectively. Results obtained in this study are compared with circular duct results using hydraulic diameter.     

KdV Limit of the Euler–Poisson System

Consider the scaling ${\varepsilon^{1/2}(x-Vt) \to x, \varepsilon^{3/2}t \to t}$ in the Euler–Poisson system for ion-acoustic waves (1). We establish that as ${\varepsilon \to 0}$ , the solutions to such Euler–Poisson systems converge globally in time to the solutions of the Korteweg–de Vries equation.     

Singular Perturbation of Boundary Value Problem for Quasilinear Third-Order Ordinary Differential Equations Involving Two Small Parameters

ＩｎｔｒｏｄｕｃｔｉｏｎＡｂｏｕｔｓｉｎｇｕｌａｒｐｅｒｔｕｒｂａｔｉｏｎｏｆｂｏｕｎｄａｒｙｖａｌｕｅｐｒｏｂｌｅｍｆｏｒｓｅｃｏｎｄ＿ｏｒｄｅｒｏｒｄｉｎａｒｙｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｗｏｒｋｓｏｆａｌａｒｇｅｎｕｍｂｅｒｗｅｒｅｃｏｎｓｉｄｅｒｅｄｏｎｌｙｆｏｒｔｈｅｃａｓｅｏｆｉｎｖｏｌｖｉｎｇｏｎｅｓｍａｌｌｐａｒａｍｅｔｅｒ[1- 4].Ｏｎｌｙａｆｅｗｗｏｒｋｓｗｅｒｅｃｏｎｓｉｄｅｒｅｄｆｏｒｔｈｅｃａｓｅｏｆｉｎｖｏｌｖｉｎｇｔｗｏｓｍａｌｌｐａｒａｍｅｔｅｒｓ[5 - 8]…     

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.     

Homogenization of Non-linear Scalar Conservation Laws

We study the limit as ε → 0 of the entropy solutions of the equation . We prove that the sequence u ^ε two-scale converges toward a function u(t, x, y), and u is the unique solution of a limit evolution problem. The remarkable point is that the limit problem is not a scalar conservation law, but rather a kinetic equation in which the macroscopic and microscopic variables are mixed. We also prove a strong convergence result in .