Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity steps”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Organic Iodide Sorption from Dilute Gas Streams

Reprocessing used nuclear fuel releases volatile radionuclides, including 129 iodine (I), into the off-gas of a processing plant. Volatile radioiodine could be present in several forms, depending on the chemistry of the process used and the off-gas stream. Inorganic I 2 is expected to be the predominant I species in the dissolver off-gas (DOG), with minor organic iodides present. The bulk of the I is expected to volatilize into the DOG in parts-per-million-level (ppm) concentrations. In contrast, in the vessel off-gas (VOG), most of the volatile I is expected to be found as organic iodides, such as CH 3 I, C 4 H 9 I, and C 12 H 25 I. These species are expected to be present in parts-per-billion-level (ppb) concentrations but require abatement, even at low expected concentrations, to meet regulatory emissions limits in the United States. Historically, studies of I abatement by Ag-functionalized sorbents have focused on inorganic I in the DOG, but in the last few years, more research attention has been given to organic iodides, especially longer chain species, such as C 4 H 9 I, and C 12 H 25 I. This report has three main goals: (1) to present new data generated at Oak Ridge National Laboratory (ORNL) in FY21 on the sorption behavior of organic iodides on AgZ, (2) to summarize and synthesize organic iodide data produced by ORNL and Idaho National Laboratory (INL) over the last 4 years to answer questions on organic iodides behavior outlined in the 2018 joint test plan (Jubin et al. 2018), and (3) to propose a VOG abatement system design that can provide the capture efficiencies required to meet I emission limits. The 2021 ORNL experimental campaign tested the effects of organic iodide speciation and concentration in the off-gas, superficial velocity of the off-gas, and effects of aging on AgZ sorbent capacity. These studies found that the sorption rate of organic iodides by AgZ depends on the hydrocarbon chain length and the concentration in the off-gas. Higher molecular weight organic iodides adsorb to AgZ more slowly than I. At a concentration of 50 ppm concentration in the off-gas, CH 3 I loads 8% slower, C 4 H 9 I loads 20% slower, and C 12 H 25 I loads 40% slower than I. The lowest concentration loading rates calculated were in 5 ppm organic iodide gas streams in which AgZ gained on average 0.14 mg I/g sorbent/hour in the bench scale test system. Thus, longer sorbent beds might be needed to accommodate slower loading rates onto AgZ in lower concentration gas streams. Although sorption rate varies as a function of hydrocarbon chain length, the saturation concentration of the sorbent for these I-bearing species does not vary. Aging AgZ in a humid air stream for 9 months drops the overall sorbent capacity by ~35% for CH 3 I, ~50% for C 4 H 9 I, and ~40% for C 12 H 25 I. This results in a saturation capacity between 35 and 70 mg I/g sorbent for the aged AgZ. In conjunction with recent data produced by INL, these data are used to estimate the mass transfer zone (MTZ) and decontamination factor (DF) for sorbent beds of AgZ. Sorption tests performed with iodide gas concentrations of about 1 ppm and higher at a superficial gas velocity of 10 m/min, indicate that MTZ depths for these conditions tend to range between about 8-20 cm. Tests performed at lower concentrations between 50-90 ppb and at gas superficial velocities of 1, 10, and 20 m/min indicate that the MTZ depth increases with increasing superficial gas velocity. The 20 m/min test indicates that the MTZ for those conditions was at least 22 cm, and doubling the superficial gas velocity from 10 to 20 m/min could roughly double or triple the MTZ depth. Doubling and tripling the bounding MTZ depth of 20 cm for the body of MTZ estimates made at with 10 m/min superficial velocity would extend the MTZ for a superficial velocity of 20 m/min to 40-60 cm. This bounding limit applies to all of the organic iodides that have been tested. These results also indicate that the sorption rate-limiting step is not sensitive to the superficial gas velocity; otherwise the MTZ depth would not have increased approximately in proportion to the increase in the gas superficial velocity. This further suggests that the rate limiting step is not associated with mass transfer of the sorbate to the sorbent surface, or mass transfer of the reaction byproducts from the sorbent surface, but is associated with sorption or chemical reactions on the sorbent surface or in sorbent pores. Deep bed testing at INL has established DFs of >2,000 for I, CH 3 I, and C 4 H 9 I under a range of conditions (Soelberg et al. 2021, Bruffey et al. 2019). DF does not seem to be affected by the concentration of the organic iodide in the gas stream over the range of 1 to 50 ppm. Thus, if the MTZ is accommodated in sorbent bed design for the DOG and VOG, then regulatory DFs will be met. To meet the third objective outlined in this report, these experimental data were used to update an engineering evaluation of the VOG first completed in 2016. The updated VOG design can be found in an accompanying document (Welty et al., 2021; INL- LTD-21-64587). This report finds that the VOG will decrease in both size and complexity, relative to previous designs, and will still meet regulatory requirements for all iodine forms.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of k{sub eff} in the unfolding process of the LOTUS reactor design

An innovative small thermal reactor concept named LOTUS, which is a pool-type reactor aimed at heat supply. The passive concept applied in the LOTUS and it has a separable core. It means that its core can passively unfold into four parts along its central axis only with the buoyancy and gravity when some accidents happen, especially the wall rupture of the pool caused by earthquake, aircraft impact and so on. This idea can effectively enhance the ability of pool-type reactor for resisting some severe accidents. We have executed some pre-computations in order to check a non-supercritical system for the transient simulation. In this paper, the development of k{sub eff} have been calculated by both stepwise steady-state and transient modes using Serpent, as well as including some sensitivity analyses about the k{sub eff} to unfolding velocity, inserting velocity of control rod (CR) and time step. The result showed that k{sub eff} value decreased definitely as the reactor unfolded and the consistent change exposed between the stepwise steady state and transient result when the unfolding angle varied from 0 to 1 degree. However, the stepwise steady state k{sub eff} negatively declined after that angle, which demonstrated that for such a dynamic unfolding change the transient mode was more suitable, peculiarly for considering delayed neutrons and neutron moderation. For sensitivity calculation, the k{sub eff} decreased more rapidly as the unfolding velocity and inserting velocity of CR increased, and a smaller time step corresponded to a more stable change. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Notes on Interface Growth Kinetics 50 Years after Burton, Cabrera and Frank

This overview is devoted to some unresolved basic problems in crystal growth kinetics. The density wave approach to propagation of a spatially diffuse interface between a growing crystal and its simple (e.g., metallic) melt is discussed is Section 2. This approach allows for the calculation of kinetic coefficients and is an alternative to the localized interface concept in which each atom belongs to either a solid or a liquid. Sections 3 and 4 deals mainly with layer growth from solution. Mutual retardation of the growth steps via their bulk and d a c e diffusion fields is the major subject. The influence of solution flow on step bunching (Section 4) suggests the essential influence of bulk diffusion on the surface morphology. The flow within the solution boundary layer enhances step-step interaction, influences the step bunching process and the resulting step pattern morphology on the growing surface. Recent experiments on the rates at which strongly polygonized steps on protein and small molecule crystals propagate during growth from solution are analyzed in Section 5 . We have shown that the step segments may be "singular" and that "one-dimensional nucleation" may be the rate limiting stage for the segments that are shorter or comparable in length to the thermodynamically equilibrium interlink distance. In this case, the reciprocal dependence of the segment propagation rate on the segment length that follow from the Gibbs-Thompson law, should be replaced by an abrupt switch from zero to a finite constant velocity. Until recently, the Kossel crystal remained the only model used in crystal growth theory. In such Kossel Gibbs-Thomson law, should be replaced by an abrupt switch &om zero to a finite constant velocity. crystals, all kinks at the steps are identical and the kink rate is a linear function of the supersaturation. In the non-Kossel crystals, there may be several kink configurations characterized by different geometries and energies. These configurations should appear in a specific sequence when each new lattice unit cell is filled. As a result of such a cooperative interaction within the unit cell, a non-linear dependence of the kink rate on the vapor pressure or solution concentration in excess over the equilibrium value should be expected.

Chernov, A. A.↗

Inflow velocity perturbations due to fuselage effects in the presence of a fully interactive wake

The velocity field at locations above the plane of the rotor disk of a generic helicopter configuration has been measured experimentally and modeled using two computational models. The first computational model represents only the portion of the velocity field due to the fuselage, using only the potential due to freestream onset flow. The second model represents all elements of the problem in an incompressible field: the rotor, the rotor wake, and the non-lifting fuselage. The first model is used to assess the fidelity requirement for fuselage paneling. Comparisons are made between a much simplified body shape, an accurate body shape, a body shape with a crude hub/shaft model and wind-tunnel measurements. The second model is used to assess the relative magnitude of the velocity perturbations at the rotor disk due to the effects of the fuselage. A fully interactive time-stepping wake method is used to predict the velocities at the rotor disk. These velocities are computed with and without the fuselage in the model and compared with measured wind-tunnel inflow velocities. The difference between the velocities computed with and without the fuselage is then compared with the perturbation velocity field computed by the first method.

Berry, John D.↗

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The contribution of Humboldt Glacier, northern Greenland, to sea-level rise through 2100 constrained by recent observations of speedup and retreat

Abstract. Humboldt Glacier, northern Greenland, has retreated and accelerated through the 21st century, raising concerns that it could be a significant contributor to future sea-level rise. We use a data-constrained ensemble of three-dimensional higher-order ice sheet model simulations to estimate the likely range of sea-level rise from the continued retreat of Humboldt Glacier. We first solve for basal traction using observed ice thickness, bed topography, and ice surface velocity from the year 2007 in a PDE-constrained (partial differential equation) optimization. Next, we impose calving rates to match mean observed retreat rates from winter 2007–2008 to winter 2017–2018 in a transient calibration of the exponent in the power-law basal friction relationship. We find that power-law exponents in the range of 1/7–1/5 – rather than the commonly used 1/3–1 – are necessary to reproduce the observed speedup over this period. We then tune an iceberg calving parameterization based on the von Mises stress yield criterion in another transient-calibration step to approximate both observed ice velocities and terminus position in 2017–2018. Finally, we use the range of basal friction relationship exponents and calving parameter values to generate the ensemble of model simulations from 2007–2100 under three climate forcing scenarios from CMIP5 (two RCP8.5 forcings, Representative Concentration Pathway) and CMIP6 (one SSP5-8.5 forcing, Shared Socioeconomic Pathway). Our simulations predict 5.2–8.7 mm of sea-level rise from Humboldt Glacier, significantly higher than a previous estimate (∼ 3.5 mm) and equivalent to a substantial fraction of the 40–140 mm predicted by ISMIP6 from the whole Greenland Ice Sheet. Our larger future sea-level rise prediction results from the transient calibration of our basal friction law to match the observed speedup, which requires a semi-plastic bed rheology. In many simulations, our model predicts the growth of a sizable ice shelf in the middle of the 21st century. Thus, atmospheric warming could lead to more retreat than predicted here if increased surface melt promotes hydrofracture of the ice shelf. Our data-constrained simulations of Humboldt Glacier underscore the sensitivity of model predictions of Greenland outlet glacier response to warming to choices of basal shear stress and iceberg calving parameterizations. Further, transient calibration of these parameterizations, which has not typically been performed, is necessary to reproduce observed behavior. Current estimates of future sea-level rise from the Greenland Ice Sheet could, therefore, contain significant biases.

54 ENVIRONMENTAL SCIENCES↗

COBE nonspinning attitude propagation

The Cosmic Background Explorer (COBE) spacecraft will exhibit complex attitude motion consisting of a spin rate of approximately -0.8 revolution per minute (rpm) about the x-axis and simultaneous precession of the spin axis at a rate of one revolution per orbit (rpo) about the nearly perpendicular spacecraft-to-Sun vector. The effect of the combined spinning and precession is to make accurate attitude propagation difficult and the 1-degree (3 sigma) solution accuracy goal problematic. To improve this situation, an intermediate reference frame is introduced, and the angular velocity divided into two parts. The nonspinning part is that which would be observed if there were no rotation about the X-axis. The spinning part is simply the X-axis component of the angular velocity. The two are propagated independently and combined whenever the complete attitude is needed. This approach is better than the usual one-step method because each of the two angular velocities look nearly constant in their respective reference frames. Since the angular velocities are almost constant, the approximations made in discrete time propagation are more nearly true. To demonstrate the advantages of this nonspinning method, attitude is propagated as outlined above and is then compared with the results of the one-step method. Over the 100-minute COBE orbit, the one-step error grows to several degrees while the nonspinning error remains negligible.

Chu, D.↗

Structure of confined laminar spray diffusion flames: Numerical investigation

The structure of confined laminar spray diffusion flames is investigated numerically by solving the gas-phase conservation equations for mass species, continuity, momentum, and energy and the liquid-phase equations for droplet position, velocity, size, and temperature. A one-step global reaction scheme along with six equilibrium reactions are employed to model the flame chemistry. Monodisperse as well as polydisperse sprays are considered. The numerical results demonstrate that liquid spray flames substantially differ from gaseous flames in their structure, i.e., temperature, concentration, and velocity fields, shape, and dimensions under the same conditions. Spray flames are predicted to be taller and narrower than their counterpart gaseous ones and their shapes are almost cylindrical. This is in agreement with experimental observations. The numerical computations also show that the use of the equilibrium reactions with the one-step reaction scheme decreases the flame temperature compared to the one-step reaction scheme without the equilibrium reactions and more importantly increases the surface area of the flame zone due to a phenomenon termed 'equilibrium broadening.' The spray flames also possess a finite thickness with minimal overlap of the fuel and oxygen species. A case for which a fuel-mixture consisting of 20 to 80 percent gas-liquid by mass is introduced into the combustor is also investigated and compared with predictions using only gaseous or liquid fuel.

Mawid, M. A.↗

Applying Motion Constraints Based on Test Data

MSC ADAMS is a simulation software that is used to analyze multibody dynamics. Using user subroutines, it is possible to apply motion constraints to the rigid bodies so that they match the motion profile collected from test data. This presentation describes the process of taking test data and passing it to ADAMS using user subroutines, and uses the Morpheus free-flight 4 test as an example of motion data used for this purpose. Morpheus is the name of a prototype lander vehicle built by NASA that serves as a test bed for various experimental technologies (see backup slides for details) MSC.ADAMS"TM" is used to play back telemetry data (vehicle orientation and position) from each test as the inputs to a 6-DoF general motion constraint (details in backup slides) The MSC.ADAMS"TM" playback simulations allow engineers to examine and analyze flight trajectory as well as observe vehicle motion from any angle and at any playback speed. This facilitates the development of robust and stable control algorithms, increasing reliability and reducing development costs of this developmental engine The simulation also incorporates a 3D model of the artificial hazard field, allowing engineers to visualize and measure performance of the developmental autonomous landing and hazard avoidance technology ADAMS is a multibody dynamics solver. It uses forces, constraints, and mass properties to numerically integrate equations of motion. The ADAMS solver will ask the motion subroutine for position, velocity, and acceleration values at various time steps. Those values must be continuous over the whole time domain. Each degree of freedom in the telemetry data can be examined separately; however, linear interpolation of the telemetry data is invalid, since there will be discontinuities in velocity and acceleration.

Burlone, Michael↗

Analysis of neutrally stable atmospheric flow over a two-dimensional forward facing step

A mathematical model describing ground-wind induced flow fields around surface obstructions such as buildings, bridges or other man-made structures is presented. The Navier-Stokes equations with a two-equation turbulence model are used to analyze the flow over a two-dimensional forward facing step. A study of the assumed logarithmic approach velocity profile shows that an increase of surface roughness produces a rise to higher turbulence levels in the shear layer originating from the step corner. This rise leads to higher Reynolds stress and to faster reattachment of the separated flow.

Bitte, J.↗

Homologous mutations in human β, embryonic, and perinatal muscle myosins have divergent effects on molecular power generation

Mutations at a highly conserved homologous residue in three closely related muscle myosins cause three distinct diseases involving muscle defects: R671C in β-cardiac myosin causes hypertrophic cardiomyopathy, R672C and R672H in embryonic skeletal myosin cause Freeman–Sheldon syndrome, and R674Q in perinatal skeletal myosin causes trismus-pseudocamptodactyly syndrome. It is not known whether their effects at the molecular level are similar to one another or correlate with disease phenotype and severity. To this end, we investigated the effects of the homologous mutations on key factors of molecular power production using recombinantly expressed human β, embryonic, and perinatal myosin subfragment-1. We found large effects in the developmental myosins but minimal effects in β myosin, and magnitude of changes correlated partially with clinical severity. The mutations in the developmental myosins dramatically decreased the step size and load-sensitive actin-detachment rate of single molecules measured by optical tweezers, in addition to decreasing overall enzymatic (ATPase) cycle rate. In contrast, the only measured effect of R671C in β myosin was a larger step size. Our measurements of step size and bound times predicted velocities consistent with those measured in an in vitro motility assay. Finally, molecular dynamics simulations predicted that the arginine to cysteine mutation in embryonic, but not β, myosin may reduce pre-powerstroke lever arm priming and ADP pocket opening, providing a possible structural mechanism consistent with the experimental observations. This paper presents direct comparisons of homologous mutations in several different myosin isoforms, whose divergent functional effects are a testament to myosin’s highly allosteric nature.

59 BASIC BIOLOGICAL SCIENCES↗

Analysis of three-dimensional unsteady flow around oscillating wings

A method based on the Navier-Stokes equations was developed for determining analytically the three-dimensional unsteady flow patterns around oscillating wings. The Helmholz vorticity transport equations were discretized in three-dimensional finite element form from a variational formulation and integrated numerically. At each time step of the numerical integration the velocity field was calculated from the representation of the three-dimensional wing by a system of optimized distribution of vortices in space. During the numerical integration of the vorticity transport equations the time-dependent boundary conditions on the wing were specified as external constraint conditions. Examples of obtained results describing the three-dimensional unsteady flow around a wing were presented.

Bratanow, T.↗

Interaction of decaying trailing vortices in ground shear

The drift of trailing vortices in a cross wind near to the ground is modeled by an unsteady, two-dimensional, rotational flowfield with a concentration of large vorticity in spots having finite total strength. The problem is analyzed by a combination of the method of matched asymptotic analyses for the decay of the vortical spots and the Euler solution for the unsteady rotational flow. A primary feature of the numerical scheme is that the grid size and time step depend only the length and velocity scales of the background flow and is independent of the effective core size of a vortical spot which can be much smaller than the grid size. Numerical results are presented to demonstrated the strong interaction between the trajectjories of the vortical spots and the redistribution of vorticity in the background flowfield.

Liu, C. H.↗

Asymptotic structure and similarity solutions for three-dimensional turbulent boundary layers

The asymptotic structure of the three-dimensional turbulent boundary layer is investigated in the limit of large Reynolds numbers. A self-consistent, but relatively complex, two-layer structure exists and the simplest situation, corresponding to a plane of symmetry, is considered in this paper as a first step. The adjustment of the streamwise velocity to relative rest, through an outer defect layer and then an inner wall layer, is similar to that in two-dimensional flow. The adjustment of the cross-streamwise velocity is more complicated and it is shown that two terms in the expansion are required to obtain useful results, and in particular to obtain the velocity skew angle at the wall near the symmetry plane. The conditions under which self-similarity is achieved near a plane of symmetry are investigated. A set of ordinary differential equations is developed which describe the streamwise and cross-streamwise velocities near a plane of symmetry in a self-similar flow through two orders of magnitude. Calculated numerical solutions of these equations yield trends which are consistent with experimental observations.

Degani, A. T.↗

Exercise countermeasures for bed-rest deconditioning

The purpose for this 30-day bed rest study was to investigate the effects of short-term, high intensity isotonic and isokinetic exercise training on maintenance of working capacity (peak oxygen uptake), muscular strength and endurance, and on orthostatic tolerance, posture and gait. Other data were collected on muscle atrophy, bone mineralization and density, endocrine analyses concerning vasoactivity and fluid-electrolyte balance, muscle intermediary metabolism, and on performance and mood of the subjects. It was concluded that: The subjects maintained a relatively stable mood, high morale, and high esprit de corps throughout the study. Performance improved in nearly all tests in almost all the subjects. Isotonic training, as opposed to isokinetic exercise training, was associated more with decreasing levels of psychological tension, concentration, and motivation; and improvement in the quality of sleep. Working capacity (peak oxygen uptake) was maintained during bed rest with isotonic exercise training; it was not maintained with isokinetic or no exercise training. In general, there was no significant decrease in strength or endurance of arm or leg muscles during bed rest, in spite of some reduction in muscle size (atrophy) of some leg muscles. There was no effect of isotonic exercise training on orthostasis, since tilt-table tolerance was reduced similarly in all three groups following bed rest. Bed rest resulted in significant decreases of postural stability and self-selected step length, stride length, and walking velocity, which were not influenced by either exercise training regimen. Most pre-bed rest responses were restored by the fourth day of recovery.

Greenleaf, John↗

Turbulent Heat Transfer Characteristics in the Shear Layer of a Separated Flow

Experiments were performed to study the evolution of the heat transfer structure in a separated free shear layer region of an incompressible separated turbulent boundary layer flow behind a backward-facing step. While there is an abundance of velocity field measurements of separated flows, heat transfer measurements are rather scarce, thus limiting assessment of the heat transfer physics and its accurate modeling. The purpose of the paper is twofold: to improve an understanding of effects of flow separation on heat transfer characteristics, and to provide data for turbulence modeling and computation. The boundary layer upstream of the step was turbulent and fully developed. A constant temperature surface boundary condition was imposed upstream and downstream of the step for the heat transfer study. An internal mixing-layer like flow forms and grows from the step lip within the original boundary layer. The turbulent structure of the flow evolving downstream, however, does not switch immediately to that of a mixing layer over the entire shear layer thickness. Measurements of mean and fluctuating velocity and temperature fields indicate that the internal layer spreads gradually in the transverse direction while the outer part of the original boundary layer is effectively unperturbed. The results in this paper have not been previously reported.

Jovic, S.↗

WINDII, the wind imaging interferometer on the Upper Atmosphere Research Satellite

The WIND imaging interferometer (WINDII) was launched on the Upper Atmosphere Research Satellite (UARS) on September 12, 1991. This joint project, sponsored by the Canadian Space Agency and the French Centre National d'Etudes Spatiales, in collaboration with NASA, has the responsibility of measuring the global wind pattern at the top of the altitude range covered by UARS. WINDII measures wind, temperature, and emission rate over the altitude range 80 to 300 km by using the visible region airglow emission from these altitudes as a target and employing optical Doppler interferometry to measure the small wavelength shifts of the narrow atomic and molecular airglow emission lines induced by the bulk velocity of the atmosphere carrying the emitting species. The instrument used is an all-glass field-widened achromatically and thermally compensated phase-stepping Michelson interferometer, along with a bare CCD detector that images the airglow limb through the interferometer. A sequence of phase-stepped images is processed to derive the wind velocity for two orthogonal view directions, yielding the vector horizontal wind. The process of data analysis, including the inversion of apparent quantities to vertical profiles, is described.

Shepherd, G. G.↗

Combustion in a turbulent mixing layer formed at a rearward-facing step

A premixed propane/air flame was stabilized in a turbulent mixing layer formed at a rearward-facing step. The mean and rms averages of the turbulent velocity flowfield were determined by laser velocimetry for both reacting (phi = 0.57) and nonreacting flows (Re = 15,000-37,000 based on step height). The reacting-flow was visualized by high-speed schlieren photography. Large-scale structures dominate the reacting mixing layer. The growth of the large-scale structures was tied to the propagation of the flame. The linear growth rate of the reacting mixing layer defined by the mean velocity profiles was unchanged by combustion but the virtual origin moves downstream. The reacting mixing layer boundaries based on the mean velocity profiles were shifted toward the recirculation zone and reattachment lengths were shortened by 30 percent. The edge of the flame controlled by the large-scale structure development propagated faster into the incoming reactants than the boundary of the mixing layer given by the mean velocity flowfield. Thus, the region of high velocity gradient did not coincide with the region of high reaction and heat transfer.

Pitz, R. W.↗