Flame acceleration due to wall friction: Accuracy and intrinsic limitations of the formulations

2015 ◽  
Vol 29 (32) ◽  
pp. 1550205 ◽  
Author(s):  
Berk Demirgok ◽  
Hayri Sezer ◽  
V’yacheslav Akkerman

The analytical formulations on the premixed flame acceleration induced by wall friction in two-dimensional (2D) channels [Bychkov et al., Phys. Rev. E 72 (2005) 046307] and cylindrical tubes [Akkerman et al., Combust. Flame 145 (2006) 206] are revisited. Specifically, pipes with one end closed are considered, with a flame front propagating from the closed pipe end to the open one. The original studies provide the analytical formulas for the basic flame and fluid characteristics such as the flame acceleration rate, the flame shape and its propagation speed, as well as the flame-generated flow velocity profile. In the present work, the accuracy of these approaches is verified, computationally, and the intrinsic limitations and validity domains of the formulations are identified. Specifically, the error diagrams are presented to demonstrate how the accuracy of the formulations depends on the thermal expansion in the combustion process and the Reynolds number associated with the flame propagation. It is shown that the 2D theory is accurate enough for a wide range of parameters. In contrast, the zeroth-order approximation for the cylindrical configuration appeared to be quite inaccurate and had to be revisited. It is subsequently demonstrated that the first-order approximation for the cylindrical geometry is very accurate for realistically large thermal expansions and Reynolds numbers. Consequently, unlike the zeroth-order approach, the first-order formulation can constitute a backbone for the comprehensive theory of the flame acceleration and detonation initiation in cylindrical tubes. Cumulatively, the accuracy of the formulations deteriorates with the reduction of the Reynolds number and thermal expansion.

1979 ◽  
Vol 46 (3) ◽  
pp. 510-512 ◽  
Author(s):  
M. B. Stewart ◽  
F. A. Morrison

Low Reynolds number flow in and about a droplet is generated by an electric field. Because the creeping flow solution is a uniformly valid zeroth-order approximation, a regular perturbation in Reynolds number is used to account for the effects of convective acceleration. The flow field and resulting deformation are predicted.


Author(s):  
Orlando Ugarte ◽  
Sinan Demir ◽  
Berk Demirgok ◽  
V’yacheslav Akkerman ◽  
Vitaly Bychkov ◽  
...  

Flame dynamics in micro-pipes have been observed to be strongly affected by the wall boundary conditions. In this respect, two mechanisms of flame acceleration are related to the momentum transferred in these regions: 1) that associated with flame stretching produced by wall friction forces; and 2) when obstacles are placed at the walls, as a result of the delayed burning occurring between them, a jet-flow is formed, intensively promoting the flame spreading. Wall thermal conditions have usually been neglected, thus restricting the cases to adiabatic wall conditions. In contrast, in the present work, the effect of the boundary conditions on the flame propagation dynamics is investigated, computationally, with the effect of wall heat losses included in the consideration. In addition, the powerful flame acceleration attained in obstructed pipes is studied in relation to the obstacle size, which determines how different this mechanism is from the wall friction. A parametric study of two-dimensional (2D) channels and cylindrical tubes, of various radiuses, with one end open is performed. The walls are subjected to slip and non-slip, adiabatic and constant temperature conditions, with different fuel mixtures described by varying the thermal expansion coefficients. Results demonstrate that higher wall temperatures promote slower propagation as they reduce the thermal expansion rate, as a result of the post-cooling of the burn matter. In turn, smaller obstacle sizes generate weaker flame acceleration, although the mechanism is noticed to be stronger than the wall friction-driven, even for the smaller sizes considered.


Author(s):  
Neville de Mestre

AbstractA perturbation model is used to predict the distance jumped by a long-jumper for a range of tailwinds and headwinds. The zeroth-order approximation is based on gravity being the only force present, the effects of drag and lift only being included in the first-order corrections. The difference in predicted distances produced by the zeroth and first-order approximations is less than 2% for headwinds or tailwinds upto 4 ms−1. Most increases or decreases due to wind are caused by changes in the run-up speed, and consequently the take-off angle and speed.


2017 ◽  
Vol 20 (05) ◽  
pp. 1750034 ◽  
Author(s):  
COLIN TURFUS ◽  
ALEXANDER SHUBERT

We present a new model for pricing contingent convertible (CoCo) bonds which facilitates the calculation of equity, credit and interest rate risk sensitivities. We assume a lognormal equity process and a Hull–White (normal) short rate process for the conversion intensity with a downward jump in the equity price on conversion. We are able to derive an approximate solution in closed form based on the assumption that the conversion intensity volatility is asymptotically small. The simple first-order approximation is seen to be accurate for a wide range of market conditions, although particularly for longer maturities higher order terms in the asymptotic expansion may be needed.


2012 ◽  
Vol 140 (2) ◽  
pp. 637-649 ◽  
Author(s):  
Max Yaremchuk ◽  
Matthew Carrier

Many background error correlation (BEC) models in data assimilation are formulated in terms of a smoothing operator [Formula: see text], which simulates the action of the correlation matrix on a state vector normalized by respective BE variances. Under such formulation, [Formula: see text] has to have a unit diagonal and requires appropriate renormalization by rescaling. The exact computation of the rescaling factors (diagonal elements of [Formula: see text]) is a computationally expensive procedure, which needs an efficient numerical approximation. In this study approximate renormalization techniques based on the Monte Carlo (MC) and Hadamard matrix (HM) methods and on the analytic approximations derived under the assumption of the local homogeneity (LHA) of [Formula: see text] are compared using realistic BEC models designed for oceanographic applications. It is shown that although the accuracy of the MC and HM methods can be improved by additional smoothing, their computational cost remains significantly higher than the LHA method, which is shown to be effective even in the zeroth-order approximation. The next approximation improves the accuracy 1.5–2 times at a moderate increase of CPU time. A heuristic relationship for the smoothing scale in two and three dimensions is proposed for the first-order LHA approximation.


Author(s):  
Ozer Elbeyli ◽  
J. Q. Sun

Semi-discretization is an effective method for analysis and control design of time-invariant as well as periodic linear systems with time delay. This paper briefly describes various approximation schemes in conjunction to the semi-discretization method. Comparison measures making use of the stability bounds of control gains and the controlled response of the system are presented. The accuracy and efficiency of the method with the zeroth order, improved zeroth order and first order approximations are compared through numerical simulations. Another first order approximation of the system with multiple time delays under a non-delayed feedback leads to integro-differential equations where a simple approximation is utilized to generate discrete maps. It is found that the first order approximation provides substantial improvements in accuracy and efficiency for feedback control design of LTI and periodic systems.


AIAA Journal ◽  
1998 ◽  
Vol 36 ◽  
pp. 1721-1727
Author(s):  
Prasanth B. Nair ◽  
Andrew J. Keane ◽  
Robin S. Langley

2021 ◽  
Vol 76 (3) ◽  
pp. 265-283
Author(s):  
G. Nath

Abstract The approximate analytical solution for the propagation of gas ionizing cylindrical blast (shock) wave in a rotational axisymmetric non-ideal gas with azimuthal or axial magnetic field is investigated. The axial and azimuthal components of fluid velocity are taken into consideration and these flow variables, magnetic field in the ambient medium are assumed to be varying according to the power laws with distance from the axis of symmetry. The shock is supposed to be strong one for the ratio C 0 V s 2 ${\left(\frac{{C}_{0}}{{V}_{s}}\right)}^{2}$ to be a negligible small quantity, where C 0 is the sound velocity in undisturbed fluid and V S is the shock velocity. In the undisturbed medium the density is assumed to be constant to obtain the similarity solution. The flow variables in power series of C 0 V s 2 ${\left(\frac{{C}_{0}}{{V}_{s}}\right)}^{2}$ are expanded to obtain the approximate analytical solutions. The first order and second order approximations to the solutions are discussed with the help of power series expansion. For the first order approximation the analytical solutions are derived. In the flow-field region behind the blast wave the distribution of the flow variables in the case of first order approximation is shown in graphs. It is observed that in the flow field region the quantity J 0 increases with an increase in the value of gas non-idealness parameter or Alfven-Mach number or rotational parameter. Hence, the non-idealness of the gas and the presence of rotation or magnetic field have decaying effect on shock wave.


Author(s):  
Afshin Anssari-Benam ◽  
Andrea Bucchi ◽  
Giuseppe Saccomandi

AbstractThe application of a newly proposed generalised neo-Hookean strain energy function to the inflation of incompressible rubber-like spherical and cylindrical shells is demonstrated in this paper. The pressure ($P$ P ) – inflation ($\lambda $ λ or $v$ v ) relationships are derived and presented for four shells: thin- and thick-walled spherical balloons, and thin- and thick-walled cylindrical tubes. Characteristics of the inflation curves predicted by the model for the four considered shells are analysed and the critical values of the model parameters for exhibiting the limit-point instability are established. The application of the model to extant experimental datasets procured from studies across 19th to 21st century will be demonstrated, showing favourable agreement between the model and the experimental data. The capability of the model to capture the two characteristic instability phenomena in the inflation of rubber-like materials, namely the limit-point and inflation-jump instabilities, will be made evident from both the theoretical analysis and curve-fitting approaches presented in this study. A comparison with the predictions of the Gent model for the considered data is also demonstrated and is shown that our presented model provides improved fits. Given the simplicity of the model, its ability to fit a wide range of experimental data and capture both limit-point and inflation-jump instabilities, we propose the application of our model to the inflation of rubber-like materials.


Author(s):  
Amarjot Singh Bhullar ◽  
Gospel Ezekiel Stewart ◽  
Robert W. Zimmerman

Abstract Most analyses of fluid flow in porous media are conducted under the assumption that the permeability is constant. In some “stress-sensitive” rock formations, however, the variation of permeability with pore fluid pressure is sufficiently large that it needs to be accounted for in the analysis. Accounting for the variation of permeability with pore pressure renders the pressure diffusion equation nonlinear and not amenable to exact analytical solutions. In this paper, the regular perturbation approach is used to develop an approximate solution to the problem of flow to a linear constant-pressure boundary, in a formation whose permeability varies exponentially with pore pressure. The perturbation parameter αD is defined to be the natural logarithm of the ratio of the initial permeability to the permeability at the outflow boundary. The zeroth-order and first-order perturbation solutions are computed, from which the flux at the outflow boundary is found. An effective permeability is then determined such that, when inserted into the analytical solution for the mathematically linear problem, it yields a flux that is exact to at least first order in αD. When compared to numerical solutions of the problem, the result has 5% accuracy out to values of αD of about 2—a much larger range of accuracy than is usually achieved in similar problems. Finally, an explanation is given of why the change of variables proposed by Kikani and Pedrosa, which leads to highly accurate zeroth-order perturbation solutions in radial flow problems, does not yield an accurate result for one-dimensional flow. Article Highlights Approximate solution for flow to a constant-pressure boundary in a porous medium whose permeability varies exponentially with pressure. The predicted flowrate is accurate to within 5% for a wide range of permeability variations. If permeability at boundary is 30% less than initial permeability, flowrate will be 10% less than predicted by constant-permeability model.


Sign in / Sign up

Export Citation Format

Share Document