RESUMEN
By analyzing the statistically stationary stage of propagation of a Huygens front in homogeneous, isotropic, constant-density turbulence, a length scale l_{0} is introduced to characterize the smallest wrinkles on the front surface in the case of a low constant speed u_{0} of the front when compared to the Kolmogorov velocity u_{K}. The length scale is derived following a hypothesis of dynamical similarity that highlights a balance between (i) creation of a front area due to advection and (ii) destruction of the front area due to propagation. Consequently, the front speed is compared with the magnitude of the fluid velocity difference in two points separated by a distance smaller than the Kolmogorov length scale. Appropriateness of the smallest wrinkle scale is demonstrated by applying a fractal approach to evaluating the mean area of the instantaneous front surface. Since the scales of the smallest and larger wrinkles belong to different subranges (dissipation and inertial, respectively) of the Kolmogorov turbulence spectrum, the front is hypothesized to be a bifractal characterized by two different fractal dimensions in the two subranges. Both fractal dimensions are evaluated adapting the aforementioned hypothesis of dynamical similarity. Such a bifractal model yields a linear relation between the mean fluid consumption velocity, which is equal to the front speed u_{0} multiplied with a ratio of the mean area of the instantaneous front surface to the transverse projected area, and the rms turbulent velocity u^{'} even if a ratio of u_{0}/u^{'} tends to zero.
RESUMEN
The present work aims at modeling the entire convection flux ρuW¯ in the transport equation for a mean reaction rate ρW¯ in a turbulent flow, which (equation) was recently put forward by the present authors. In order to model the flux, several simple closure relations are developed by introducing flow velocity conditioned to reaction zone and interpolating this velocity between two limit expressions suggested for the leading and trailing edges of the mean flame brush. Subsequently, the proposed simple closure relations for ρuW¯ are assessed by processing two sets of data obtained in earlier 3D Direct Numerical Simulation (DNS) studies of adiabatic, statistically planar, turbulent, premixed, single-step-chemistry flames characterized by unity Lewis number. One dataset consists of three cases characterized by different density ratios and is associated with the flamelet regime of premixed turbulent combustion. Another dataset consists of four cases characterized by different low Damköhler and large Karlovitz numbers. Accordingly, this dataset is associated with the thin reaction zone regime of premixed turbulent combustion. Under conditions of the former DNS, difference in the entire, ρuW¯ , and mean, uρW¯ , convection fluxes is well pronounced, with the turbulent flux, ρu''W''¯ , showing countergradient behavior in a large part of the mean flame brush. Accordingly, the gradient diffusion closure of the turbulent flux is not valid under such conditions, but some proposed simple closure relations allow us to predict the entire flux ρuW¯ reasonably well. Under conditions of the latter DNS, the difference in the entire and mean convection fluxes is less pronounced, with the aforementioned simple closure relations still resulting in sufficiently good agreement with the DNS data.
RESUMEN
In order to experimentally study whether or not the density ratio σ substantially affects flame displacement speed at low and moderate turbulent intensities, two stoichiometric methane/oxygen/nitrogen mixtures characterized by the same laminar flame speed SL = 0.36 m/s, but substantially different σ were designed using (i) preheating from Tu = 298 to 423 K in order to increase SL , but to decrease σ, and (ii) dilution with nitrogen in order to further decrease σ and to reduce SL back to the initial value. As a result, the density ratio was reduced from 7.52 to 4.95. In both reference and preheated/diluted cases, direct images of statistically spherical laminar and turbulent flames that expanded after spark ignition in the center of a large 3D cruciform burner were recorded and processed in order to evaluate the mean flame radius R¯ft and flame displacement speed St=σ-1dR¯fdt with respect to unburned gas. The use of two counter-rotating fans and perforated plates for near-isotropic turbulence generation allowed us to vary the rms turbulent velocity u' by changing the fan frequency. In this study, u' was varied from 0.14 to 1.39 m/s. For each set of initial conditions (two different mixture compositions, two different temperatures Tu , and six different u') , five (respectively, three) statistically equivalent runs were performed in turbulent (respectively, laminar) environment. The obtained experimental data do not show any significant effect of the density ratio on St . Moreover, the flame displacement speeds measured at u'/SL = 0.4 are close to the laminar flame speeds in all investigated cases. These results imply, in particular, a minor effect of the density ratio on flame displacement speed in spark ignition engines and support simulations of the engine combustion using models that (i) do not allow for effects of the density ratio on St and (ii) have been validated against experimental data obtained under the room conditions, i.e. at higher σ.
RESUMEN
The theory of turbulent diffusion of chemically reacting gaseous admixtures developed previously [T. Elperin et al., Phys. Rev. E 90, 053001 (2014)PLEEE81539-375510.1103/PhysRevE.90.053001] is generalized for large yet finite Reynolds numbers and the dependence of turbulent diffusion coefficient on two parameters, the Reynolds number and Damköhler number (which characterizes a ratio of turbulent and reaction time scales), is obtained. Three-dimensional direct numerical simulations (DNSs) of a finite-thickness reaction wave for the first-order chemical reactions propagating in forced, homogeneous, isotropic, and incompressible turbulence are performed to validate the theoretically predicted effect of chemical reactions on turbulent diffusion. It is shown that the obtained DNS results are in good agreement with the developed theory.
RESUMEN
Within the framework of the Maxwell-Cattaneo relaxation model extended to reaction-diffusion systems with nonlinear advection, travelling wave (TW) solutions are analytically investigated by studying a normalized reaction-telegraph equation in the case of the reaction and advection terms described by quadratic functions. The problem involves two governing parameters: (i) a ratio φ^{2} of the relaxation time in the Maxwell-Cattaneo model to the characteristic time scale of the reaction term, and (ii) the normalized magnitude N of the advection term. By linearizing the equation at the leading edge of the TW, (i) necessary conditions for the existence of TW solutions that are smooth in the entire interval of -∞<ζ<∞ are obtained, (ii) the smooth TW speed is shown to be less than the maximal speed φ^{-1} of the propagation of a substance, (iii) the lowest TW speed as a function of φ and N is determined. If the necessary condition of N>φ-φ^{-1} does not hold, e.g., if the magnitude N of the nonlinear advection is insufficiently high in the case of φ^{2}>1, then, the studied equation admits piecewise smooth TW solutions with sharp leading fronts that propagate at the maximal speed φ^{-1}, with the substance concentration or its spatial derivative jumping at the front. An increase in N can make the solution smooth in the entire spatial domain. Moreover, an explicit TW solution to the considered equation is found provided that N>φ. Subsequently, by invoking a principle of the maximal decay rate of TW solution at its leading edge, relevant TW solutions are selected in a domain of (φ,N) that admits the smooth TWs. Application of this principle to the studied problem yields transition from pulled (propagation speed is controlled by the TW leading edge) to pushed (propagation speed is controlled by the entire TW structure) TW solutions at N=N_{cr}=sqrt[1+φ^{2}], with the pulled (pushed) TW being relevant at smaller (larger) N. An increase in the normalized relaxation time φ^{2} results in increasing N_{cr}, thus promoting the pulled TW solutions. The domains of (φ,N) that admit either the smooth or piecewise smooth TWs are not overlapped and, therefore, the selection problem does not arise for these two types of solutions. All the aforementioned results and, in particular, the maximal-decay-rate principle or appearance of the piecewise smooth TW solutions, are validated by numerically solving the initial boundary value problem for the reaction-telegraph equation with natural initial conditions localized to a bounded spatial region.
RESUMEN
The problem of traveling wave (TW) speed selection for solutions to a generalized Murray-Burgers-KPP-Fisher parabolic equation with a strictly positive cubic reaction term is considered theoretically and the initial boundary value problem is numerically solved in order to support obtained analytical results. Depending on the magnitude of a parameter inherent in the reaction term (i) the term is either a concave function or a function with the inflection point and (ii) transition from pulled to pushed TW solution occurs due to interplay of two nonlinear terms; the reaction term and the Burgers convection term. Explicit pushed TW solutions are derived. It is shown that physically observable TW solutions, i.e., solutions obtained by solving the initial boundary value problem with a sufficiently steep initial condition, can be determined by seeking the TW solution characterized by the maximum decay rate at its leading edge. In the Appendix, the developed approach is applied to a non-linear diffusion-reaction equation that is widely used to model premixed turbulent combustion.