Investigation of Particle Deposition Characteristics in Vicinity of Laidback Fan-shaped Film Cooling Holes

To design efficient film cooling systems and mitigate particulate deposition, it is very important to know the influences of the design parameters of film cooling holes on particulate deposition. However, most previous research focused on round film cooling holes. Particle deposition characteristics near shaped film cooling holes need to be studied further. In the present study, numerical computations were carried out to simulate the particle deposition behavior on gas turbine disk samples with laidback fan-shaped film cooling holes by using CFD-DPM (Computational fluid dynamics-discrete particle method). The critical velocities for particle sticking and detachment were determined by EI-Batsh model. Compared with round holes, shaped holes mitigate the particle-wall collision for small particles (d_p≤2μm, ρ_p=990kg/m ³), but promote particle-wall collision for large particles (d_p≥4μm, ρ_p=990kg/m ³). Adding the lateral and forward expansion angle can both cause the decrease of particle deposition efficiency, however, the effect of lateral expansion angle on particle deposition is more active.     

Effect of Moisture Content on High Strain Rate Compressibility and Particle Breakage in Loose Sand

Soil-filled gabion structures are widely used to protect against the effects of blast and fragmentation. It is known that moisture content significantly affects the capability of such structures, but the behaviour of partially-saturated soils is not well characterised at the strain rates and stresses experienced in these events. In particular, little data is available for loose soils, whose compaction behaviour can have a substantial impact on structural stability and ballistic performance. This paper describes the use of split Hopkinson pressure bar experiments to characterise the pre- and post-saturation compressibility of a loose quartz sand at moisture contents of up to 15.0%. In contrast to dense soils, increases in moisture content between 0.0% and 7.5% led to a decrease in the stiffness of the sand. Above 7.5% moisture content, specimens reached full saturation during the experiment: the additional water had no further effect on the pre-saturation stiffness, but post-saturation behaviour was dominated by the stiffness of the pore water. Full saturation occurred at lower dry densities as moisture content increased, leading to a decrease in particle breakage.     

Exploring the relationship between critical state and particle shape for granular materials

The relationship between critical state and particle shape corresponds to the most fundamental aspect of the mechanics of granular materials. This paper presents an investigation into this relationship through macro-scale and micro-scale laboratory experiments in conjunction with interpretation and analysis in the framework of critical state soil mechanics. Spherical glass beads and crushed angular glass beads of different percentages were mixed with a uniform quartz sand (Fujian sand) to create a sequence of mixtures with varying particle shape. On the micro-scale, particle shape was accurately measured using a laser scanning technique, and was characterized by aspect ratio, sphericity and convexity; a new shape index, taken as the average of the three shape measures and referred to as overall regularity, was proposed to provide a collective characterization of particle shape. On the macro-scale, both undrained and drained triaxial tests were carried out to provide evidence that varying particle shape can alter the overall response as well as the critical states in both stress space and volumetric compression space. The mixtures of Fujian sand and spherical glass beads were found to be markedly more susceptible to liquefaction than the mixtures of Fujian sand and crushed angular glass beads. The change in liquefaction susceptibility was shown to be consistent with the change in the position of the critical state locus (CSL) in the compression space, manifested by a decrease in the intercept and gradient of the CSL due to the presence of spherical glass beads. Quantitative relationships have been established between each of the critical state parameters and each of the shape parameters, thereby providing a way to construct macro-scale constitutive models with intrinsic micro-scale properties built in.     

Application of the γ-Reθ Transition Model to Simulations of the Flow Past a Circular Cylinder

Based on the finite volume method, the flow past a two-dimensional circular cylinder at a critical Reynolds number (Re = 8.5 × 10⁵) was simulated using the Navier-Stokes equations and the γ-Re_θ transition model coupled with the SST k ? ω turbulence model (hereinafter abbreviated as γ-Re_θ model). Considering the effect of free-stream turbulence intensity decay, the SST k ? ω turbulence model was modified according to the ambient source term method proposed by Spalart and Rumsey, and then the modified SST k ? ω turbulence model is coupled with the γ-Re_θ transition model (hereinafter abbreviated as γ-Re_θ-SR model). The flow past a circular cylinder at different inlet turbulence intensities were simulated by the γ-Re_θ-SR model. At last, the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers were each simulated by the γ-Re_θ-SR model, and the three flow states were analyzed. It was found that compared with the SST k ? ω turbulence model, the γ-Re_θ model could simulate the transition of laminar to turbulent, resulting in better consistency with experimental result. Compared with the γ-Re_θ model, for relatively high inlet turbulence intensities, the γ-Re_θ-SR model could better simulate the flow past a circular cylinder; however the improvement almost diminished for relatively low inlet turbulence intensities The γ-Re_θ-SR model could well simulate the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers.     

Revisiting History Effects in Adverse-Pressure-Gradient Turbulent Boundary Layers

The goal of this study is to present a first step towards establishing criteria aimed at assessing whether a particular adverse-pressure-gradient (APG) turbulent boundary layer (TBL) can be considered well-behaved, i.e., whether it is independent of the inflow conditions and is exempt of numerical or experimental artifacts. To this end, we analyzed several high-quality datasets, including in-house numerical databases of APG TBLs developing over flat-plates and the suction side of a wing section, and five studies available in the literature. Due to the impact of the flow history on the particular state of the boundary layer, we developed three criteria of convergence to well-behaved conditions, to be used depending on the particular case under study. (i) In the first criterion, we develop empirical correlations defining the R e _?? -evolution of the skin-friction coefficient and the shape factor in APG TBLs with constant values of the Clauser pressure-gradient parameter β = 1 and 2 (note that β = δ ^?/τ _w dP _e /dx, where δ ^? is the displacement thickness, τ _w the wall-shear stress and dP _e /dx the streamwise pressure gradient). (ii) In the second one, we propose a predictive method to obtain the skin-friction curve corresponding to an APG TBL subjected to any streamwise evolution of β, based only on data from zero-pressure-gradient TBLs. (iii) The third method relies on the diagnostic-plot concept modified with the shape factor, which scales APG TBLs subjected to a wide range of pressure-gradient conditions. These three criteria allow to ensure the correct flow development of a particular TBL, and thus to separate history and pressure-gradient effects in the analysis.     

Wall slip and multi-tier yielding in capillary suspensions

Impact of wall slip on the yield stress measurement is examined for capillary suspensions consisting of cocoa powder as the dispersed phase, vegetable oil as the continuous primary fluid, and water as the secondary fluid using smooth and serrated parallel plates. Using dynamic oscillatory measurements, we investigated the yielding behavior of this ternary solid-fluid-fluid system with varying particle volume fraction, ?, from 0.45 to 0.65 and varying water volume fraction, ?_w, from 0.02 to 0.08. Yield stress is defined as the maximum in the elastic stress (G^′γ), which is obtained by plotting the product of elastic modulus (G^′) and strain amplitude (γ) as a function of applied strain amplitude. With serrated plates, which offer minimal slippage, capillary suspensions with ? ≥?0.45 and a fixed ?_w =?0.06 showed a two-step yielding behavior as indicated by two peaks in the plots of elastic stress as a function of strain amplitude. On the other hand with smooth plates, the capillary suspensions showed strong evidence of wall slip as evident by the presence of three distinct peaks and lowered first yield stresses for all ? and ?_w. These results can be interpreted based on the fact that a particle-depleted layer, which is known to be responsible for slip, is present in the vicinity of the smooth surfaces. The slip layer presents itself as an additional “pseudo-microstructure” (characteristic length scale) besides the two microstructures, aqueous bridges and solid particle agglomerates, that may occur in the system. With serrated plates, both the yield stresses (σ₁σ₂) and storage moduli plateau at lower strain (before the first yield point) and at higher strain (before the second yield point) (G\(^{\prime }_{p1}\), G\(^{\prime }_{p2}\)) were found to increase with ? (at a fixed ?_w =?0.06) following power-law dependences. Similarly with increasing ?_w (0.02 – 0.08) at a fixed ? =?0.62, the system behaved as a solid-like material in a jammed state with particles strongly held together as manifested by rapidly increasing σ₁ and σ₂. The usage of smooth surfaces primarily affected σ₁ which was reflected by an approximately 70–90% decrement in the measured σ₁ for all values of ?. By contrast, σ₂ and G\(^{\prime }_{p2}\) were found to be unaffected as shown by close agreement of values obtained using serrated geometry due to vanishing slip layers at higher strain amplitudes.     

Simulation of Receptivity and Induced Transition From Discrete Roughness Elements

Simulations have been carried out to predict the receptivity and growth of crossflow vortices created by Discrete Roughness Elements (DREs) The final transition to turbulence has also been examined, including the effect of DRE spacing and freestream turbulence. Measurements by Hunt and Saric (2011) of perturbation mode shape at various locations were used to validate the code in particular for the receptivity region. The WALE sub-grid stress (SGS) model was adopted for application to transitional flows, since it allows the SGS viscosity to vanish in laminar regions and in the innermost region of the boundary layer when transition begins. Simulations were carried out for two spanwise wavelengths: λ= 12mm (critical) and λ= 6mm (control) and for roughness heights (k) from 12 μm to 42 μm. The base flow considered was an ASU (67)-0315 aerofoil with 45 ⁰ sweep at -2.9 ⁰ incidence and with onset flow at a chord-based Reynolds number Re _c= 2.4x10 ⁶. For λ= 12mm results showed, in accord with the experimental data, that the disturbance amplitude growth rate was linear for k = 12 μm and 24 μm, but the growth rate was decreased for k = 36 μm Receptivity to λ= 6mm roughness showed equally good agreement with experiments, indicating that this mode disappeared after a short distance to be replaced by a critical wavelength mode. Analysis of the development of modal disturbance amplitudes with downstream distance showed regions of linear, non-linear, saturation, and secondary instability behaviour. Examination of breakdown to turbulence revealed two possible routes: the first was 2D-like transition (probably Tollmien-Schlichting waves even in the presence of crossflow vortices) when transition occurred beyond the pressure minimum; the second was a classical crossflow vortex secondary instability, leading to the formation of a turbulent wedge.     

Regularity and Singularities of Optimal Convex Shapes in the Plane

We focus here on the analysis of the regularity or singularity of solutions Ω ₀ to shape optimization problems among convex planar sets, namely:

$J(\Omega_{0})={\rm min} \{J(\Omega), \Omega \quad {\rm convex},\Omega \in \mathcal{S}_{\rm ad}\},$

where \({\mathcal{S}_{\rm ad}}\) is a set of 2-dimensional admissible shapes and \({J:\mathcal{S}_{\rm ad}\rightarrow\mathbb{R}}\) is a shape functional. Our main goal is to obtain qualitative properties of these optimal shapes by using first and second order optimality conditions, including the infinite dimensional Lagrange multiplier due to the convexity constraint. We prove two types of results:

1. i)
   under a suitable convexity property of the functional J, we prove that Ω ₀ is a W ^2,p -set, \({p\in[1, \infty]}\). This result applies, for instance, with p = ∞ when the shape functional can be written as J(Ω) = R(Ω) + P(Ω), where R(Ω) = F(|Ω|, E _f (Ω), λ₁(Ω)) involves the area |Ω|, the Dirichlet energy E _f (Ω) or the first eigenvalue of the Laplace–Dirichlet operator λ₁(Ω), and P(Ω) is the perimeter of Ω;
    

1. ii)
   under a suitable concavity assumption on the functional J, we prove that Ω ₀ is a polygon. This result applies, for instance, when the functional is now written as J(Ω) = R(Ω) ? P(Ω), with the same notations as above.
    

     

The Effect of Grain Size Distribution on Nonlinear Flow Behavior in Sandy Porous Media

The current study provides new experimental data on nonlinear flow behavior in various uniformly graded granular materials (20 samples) ranging from medium sands (\(d_{50 }>0.39\) mm) to gravel (\(d_{50}=6.3\) mm). Generally, theoretical equations relate the Forchheimer parameters a and b to the porosity, as well as the characteristic pore length, which is assumed to be the median grain size \((d_{50})\) of the porous medium. However, numerical and experimental studies show that flow resistance in porous media is largely determined by the geometry of the pore structure. In this study, the effect of the grain size distribution was analyzed using subangular-subrounded sands and approximately equal compaction grades. We have used a reference dataset of 11 uniformly graded filter sands. Mixtures of filter sands were used to obtain a slightly more well-graded composite sand (increased \(C_{u}\) values by a factor of 1.19 up to 2.32) with respect to its associated reference sand at equal median grain size \((d_{50})\) and porosity. For all composite sands, the observed flow resistance was higher than in the corresponding reference sand at equal \(d_{50}\), resulting in increased a coefficients by factors up to 1.68, as well as increased b coefficients by factors up to 1.44. A modified Ergun relationship with Ergun constants of 139.1 for A and 2.2 for B, as well as the use of \(d_{m}-\sigma \) as characteristic pore length predicted the coefficients a and b accurately.     

Slow gas motions and the negative drag of a strongly heated spherical particle

In the slow flows of a strongly and nonuniformly heated gas, in the continuum regime (Kn → 0) thermal stresses may be present. The theory of slow nonisothermal continuum gas flows with account for thermal stresses was developed in 1969–1974. The action of the thermal stresses on the gas results in certain paradoxical effects, including the reversal of the direction of the force exerted on a spherical particle in Stokes flow. The propulsion force effect is manifested at large but finite temperature differences between the particle and the gas. This study is devoted to the thermal-stress effect on the drag of a strongly heated spherical particle traveling slowly in a gas for small Knudsen numbers (M ～ Kn → 0), small but finite Reynolds numbers (Re ≤ 1), a linear temperature dependence of the transport coefficients µ ∝ T, and large but finite temperature differences ((T _w ? T _∞ )/T _M8 ～ 1). Two different systems of equations are solved numerically: the simplified Navier-Stokes equations and the modified Navier-Stokes equations with account for the thermal stresses.     

Dual yielding in capillary suspensions

Rheological measurements were performed to examine the yielding behavior of capillary suspensions prepared by mixing cocoa powder as dispersed phase, vegetable oil as the continuous primary fluid, and water as the secondary fluid. Here, we investigated the yielding behavior of solid-fluid-fluid systems with varying particle volume fraction, ?, spanning the regime from a low volume fraction (? = 0.25) to a highly filled regime (? = 0.65) using dynamic oscillatory measurements. While for ? ≤ 0.4 with a fixed water volume fraction (? _w ) of 0.06 as the secondary fluid, capillary suspensions exhibited a single yield point due to rupturing of aqueous capillary bridges between the particles, while capillary suspensions with ? ≥ 0.45 showed a two-step yielding behavior. On plotting elastic stress (G^′ γ) as a function of applied strain (γ), two distinct peaks, indicating two yield stresses, were observed. Both the yield stresses and storage modulus at low strains were found to increase with ? following a power law dependence. With increasing ? _w (0 – 0.08) at a fixed ? = 0.65, the system shifted to a frustrated, jammed state with particles strongly held together shown by rapidly increasing first and second yield stresses. In particular, the first yield stress was found to increase with ? _w following a power law dependence, while the second yield stress was found to increase exponentially with ? _w . Transient steady shear tests were also performed. The single stress overshoot for ? ≤ 0.4 with ? _w = 0.06 reflected one-step yielding behavior. In contrast, for high ? (≥ 0.45) values with ? _w = 0.06, two stress overshoots were observed in agreement with the two-step yielding behavior shown in the dynamic oscillatory measurements. Experiments on the effect of resting time on microstructure recovery demonstrated that aggregates could reform after resting under quiescent conditions.     

A simple criterion for finite time stability with application to impacted buckling of elastic columns

The existing theories of finite-time stability depend on a prescribed bound on initial disturbances and a prescribed threshold for allowable responses. It remains a challenge to identify the critical value of loading parameter for finite time instability observed in experiments without the need of specifying any prescribed threshold for allowable responses. Based on an energy balance analysis of a simple dynamic system, this paper proposes a general criterion for finite time stability which indicates that finite time stability of a linear dynamic system with constant coefficients during a given time interval [0, t _f ] is guaranteed provided the product of its maximum growth rate (determined by the maximum eigen-root p₁ >0) and the duration t _f does not exceed 2, i.e., p₁t _f <2. The proposed criterion (p₁t _f =2) is applied to several problems of impacted buckling of elastic columns: (i) an elastic column impacted by a striking mass, (ii) longitudinal impact of an elastic column on a rigid wall, and (iii) an elastic column compressed at a constant speed (“Hoff problem”), in which the time-varying axial force is replaced approximately by its average value over the time duration. Comparison of critical parameters predicted by the proposed criterion with available experimental and simulation data shows that the proposed criterion is in robust reasonable agreement with the known data, which suggests that the proposed simple criterion (p₁t _f =2) can be used to estimate critical parameters for finite time stability of dynamic systems governed by linear equations with constant coefficients.     

Large amplitude oscillatory shear rheology of three different shear-thickening particle dispersions

We present a large amplitude oscillatory shear rheology (LAOS) investigation of three different shear-thickening particle dispersions - fumed silica in polyethylene oxide (FLOC), fumed silica in polypropylene glycol (HydroC), and cornstarch in water (JAM). These systems shear-thicken by three different mechanisms - shear-induced formation of particle clusters flocculated by polymer bridging, hydrocluster formation, and jamming. The viscoelastic non-linearities of the three fluids were studied as a function of strain and strain-rate space through the use of Lissajous-Bowditch curves and local nonlinear viscoelastic moduli of an oscillatory shear cycle. The nonlinear behaviors of the three fluids were compared and contrasted to understand the nonlinear shear-thickening mechanism of each. Both HydroC and JAM dispersions were found to exhibit strong strain stiffening of the elastic moduli and strain thickening of the loss moduli behavior associated with possible hydrocluster formation and particle jamming. However, the FLOC dispersion, in contrast, showed strong strain softening and strain thinning behavior at large strain amplitudes associated with yielding of the microstructure. The expected thickening of the loss modulus of FLOC in LAOS with increasing strain was not observed even though viscosity of FLOC was found to shear-thicken in steady-shear measurements. This disagreement is likely due to very large strain amplitudes required for shear-thickening to occur by shear-induced polymer bridging mechanism. The hypothesis was confirmed through stress growth experiments. Conversely, the HydroC and JAM dispersions required relatively small applied strains for shear-thickening to occur by hydrocluster and jamming mechanism. The comparison of local intra-cycle nonlinearity through Lissajous-Bowditch plots and nonlinear viscoelastic parameters indicated that the elastic nonlinearities of all three systems are primarily driven by a strong dependence on the magnitude of the applied strain-rates within an oscillatory cycle rather than the amplitude of the applied strain. A close inspection of the LAOS data reveals strong differences in the viscoelastic nonlinearities of these three different shear-thickening dispersions which can be used to create a nonlinear rheological fingerprint for each and offers valuable new insights into the nonlinear dynamics associated with each of the shear-thickening mechanisms.     

Effects of the Reynolds number on a scale-similarity model of Lagrangian velocity correlations in isotropic turbulent flows

A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZHAO, X. Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence. Physical Review E, 80, 066313(2009)). The model can be expressed as a two-point Eulerian space correlation and the dispersion velocity V. The dispersion velocity denotes the rate at which one moving particle departs from another fixed particle. This paper numerically validates the robustness of the scale-similarity model at high Taylor micro-scale Reynolds numbers up to 373, which are much higher than the original values(R_λ = 66, 102). The effect of the Reynolds number on the dispersion velocity in the scale-similarity model is carefully investigated. The results show that the scale-similarity model is more accurate at higher Reynolds numbers because the two-point Lagrangian velocity correlations with different initial spatial separations collapse into a universal form compared with a combination of the initial separation and the temporal separation via the dispersion velocity.Moreover, the dispersion velocity V normalized by the Kolmogorov velocity V_η≡η/τ_η in which η and τ_η are the Kolmogorov space and time scales, respectively, scales with the Reynolds number R_λ as V/V_η∝ R_λ~(1.39) obtained from the numerical data.     

Fokker-Planck operator for free-molecular,nonspherical, thermally nonequilibrium Brownian particles

The Fokker-Planck collision operator entering into the kinetic equation for the velocity distribution function of convex nonspherical rigid homogeneous Brownian particles in a traveling inhomogeneous monatomic gas is derived. The regime of flow past the particles is free-molecular, i.e., the characteristic particle dimensions are much smaller than the average free path of the gas molecules, the interaction between the particles and their effect on the gas phase can be neglected, and a specular-diffuse law of interaction between the molecules and the particle surface holds. The particle temperatures T _p are the same and differ from the local gas temperature T. Thermal nonequilibrium (T _p ≠ T) leads to violation of the well-known relations between the diffusion coefficients in the spaces of the translational and angular velocities and the coefficients of forces and moments exerted on a particle. The coefficients in the unknown operator are calculated for the particles in the form of bodies of revolution with longitudinal symmetry. Data characterizing the effect of the nonsphericity of the particles, i.e., spheroids and sphere-cylinders, on the degree of violation are given.     

Operator method for solving the plane elasticity problem for a strip with periodic cuts

In the present paper, we use the conformal mapping z/c = ζ?2a sin ζ (a, c?const, ζ = u + iv) of the strip {|v| ≤ v ₀, |u| < ∞} onto the domain D, which is a strip with symmetric periodic cuts. For the domain D, in the orthogonal system of isometric coordinates u, v, we solve the plane elasticity problem. We seek the biharmonic function in the form F = C ψ ₀ + S ψ*₀ + x(C ψ ₁ ? S ψ ₂) + y(C ψ ₂ + S ψ ₁), where C(v) and S(v) are the operator functions described in [1] and ψ ₀(u), …, ψ ₂(u) are the desired functions. The boundary conditions for the function F posed for v = ±v ₀ are equivalent to two operator equations for ψ ₁(u) and ψ ₂(u) and to two ordinary differential equations of first order for ψ ₀(u) and ψ*₀(u) [2]. By finding the functions ψ _j (u) in the form of trigonometric series with indeterminate coefficients and by solving the operator equations, we obtain infinite systems of linear equations for the unknown coefficients. We present an efficient method for solving these systems, which is based on studying stable recursive relations. In the present paper, we give an example of analysis of a specific strip (a = 1/4, v ₀ = 1) loaded on the boundary v = v ₀ by a normal load of intensity p. We find the particular solutions corresponding to the extension of the strip by the longitudinal force X ^∞ and to the transverse and pure bending of the strip due to the transverse force Y ^∞ and the constant moment M ^∞, respectively. We also present the graphs of normal and tangential stresses in the transverse cross-section x = 0 and study the stress concentration effect near the cut bottom.     

Experimental Study on Three-Dimensional Deformation Field of Portevin–Le Chatelier Effect Using Digital Image Correlation

In view of its high precision and high efficiency, three-dimensional digital image correlation (3D-DIC) is widely used to accurately measure full-field deformation. A spatiotemporal experimental study using 3D-DIC to explore the Portevin–Le Chatelier (PLC) deformational behavior, provides a new insight into the whole 3D deformation field, including the out-of-plane displacement, and in particular the relationship between the serrations and the strain field in the deformation bands corresponding to individual serrations. Specimens 1, 2 and 3 mm thick of 5456 Al-based alloy were tested in uniaxial tension at room temperature at strain rates from 1.8 × 10^?4 to 9.1 × 10^?3s^?1. The spatial and temporal characteristics of the strain localization were quantitatively analyzed. The out-of-plane displacement increment field (w) of the localized bands was observed by 3D-DIC, and found to be related to the specimen thickness and the in-plane strain increment. The largest displacement increments were respectively 15, 10 and 5 μm for 3, 2 and 1 mm specimens at maximum strain increment of about 12000 με. The elastic shrinkage outside the deformation bands was found to be an essential characteristic of the PLC effect. The width of the PLC band (w_band) increased with increasing thickness; the angle of the PLC band (??_band) was not affected by either specimen thickness or serration amplitude. Temporally, the serrations in the plots both of in-plane strain and out-of-plane displacement vs. time coincided throughout the entire loading procedure. When PLC banding occurred, the serration amplitude within the bands was found to be proportional to the maximum strain increment in the direction of the applied tensile force (??_max).     

Effects of Fuel Lewis Number on Localised Forced Ignition of Globally Stoichiometric Stratified Mixtures: a Numerical Investigation

The influences of fuel Lewis number Le_F on localised forced ignition of globally stoichiometric stratified mixtures have been analysed using three-dimensional compressible Direct Numerical Simulations (DNS) for cases with Le_F ranging from 0.8 to 1.2. The globally stoichiometric stratified mixtures with different values of root-mean-square (rms) equivalence ratio fluctuation (i.e. ?^′= 0.2, 0.4 and 0.6) and the Taylor micro-scale l_? of equivalence ratio ? variation (i.e. l_?/l_f= 2.1, 5.5 and 8.3 with l_f being the Zel’dovich flame thickness of the stoichiometric laminar premixed flame) have been considered for different initial rms values of turbulent velocity u^′. A pseudo-spectral method is used to initialise the equivalence ratio variation following a presumed bi-modal distribution for prescribed values of ?^′ and l_?/l_f for global mean equivalence ratio 〈?〉=1.0. The localised ignition is accounted for by a source term in the energy transport equation that deposits energy for a stipulated time interval. It has been observed that the maximum values of temperature and the fuel reaction rate magnitude increase with decreasing Le_F during the period of external energy deposition. The initial values of Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f have been found to have significant effects on the extent of burning of the stratified mixtures following localised ignition. For a given value of u^′/S_b(?=1), the extent of burning decreases with increasing Le_F. An increase in u^′ leads to a monotonic reduction in the burned gas mass for all values of Le_F in all stratified mixture cases but an opposite trend is observed for the Le_F=0.8 homogeneous mixture. It has been found that an increase in ?^′ has adverse effects on the burned gas mass, whereas the effects of l_?/l_f on the extent of burning are non-monotonic and dependent on ?^′ and Le_F. Detailed physical explanations have been provided for the observed Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f dependences.     

Linearized Stability for Semiflows Generated by a Class of Neutral Equations,with Applications to State-Dependent Delays

We prove a principle of linearized stability for semiflows generated by neutral functional differential equations of the form x′(t) = g(? x _t , x _t ). The state space is a closed subset in a manifold of C ²-functions. Applications include equations with state-dependent delay, as for example x′(t) = a x′(t + d(x(t))) + f (x(t + r(x(t)))) with \({a\in\mathbb{R}, d:\mathbb{R}\to(-h,0), f:\mathbb{R}\to\mathbb{R}, r:\mathbb{R}\to[-h,0]}\).     

Fracture in Microscale SU-8 Polymer Thin Films

In this work, the fracture of spin coated SU-8 epoxy thin films was investigated under mode I loading using in situ optical experiments on specimens with double edge notched tensile geometry. A method was developed to fabricate 3 μm thick SU-8 films with tapered Chevron type notches using a combination of electron beam and ultra-violet lithography techniques. Subsequently, through speckle patterning under tensile loading, the local deformation fields around the crack tip were extracted using digital image correlation. Since the notches were tapered through the thickness, a crack nucleated from them and grew stably until it spanned the entire thickness before propagating unstably leading to catastrophic failure. As SU-8 underwent brittle fracture with no evidence of a large process zone, the critical energy release rate, J _{I C} was computed from deformation fields, and was found to be 106.6 ± 12.03 J /m ². As the film thickness was small compared to lateral dimensions, assuming plane stress conditions, the critical stress intensity factor was calculated as 0.57 ± 0.03 MPa\(\sqrt {m}\). Furthermore, to assess the validity of the experimental method, a finite element simulation on the exact specimen geometry was conducted with experimentally evaluated far field displacement boundary conditions. The strain fields and J-integral value obtained from the simulation were in good agreement with the experimental results, implying the validity of the in situ experimental method proposed given the challenges of small scale specimens. Furthermore, using fractography and optical imaging it was confirmed that the unstable crack propagation started once the crack front reached full thickness, thereby providing sharp crack at the time of failure, which is necessary for brittle materials for valid fracture toughness experiment. It is expected that the proposed methods of specimen preparation and fracture experiments on microscale polymer thin films can be used on other materials.