Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mass density ratio”

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 415 records · Page 23

Wind variability in the Large Magellanic Cloud B(e) supergiant HD 34664

We present optical linear spectropolarimetry of the B(e) supergiant HD 34664. The polarization and position angle display significant spectral features which correspond to the H-alpha and H-beta emission lines. We use the line polarizations to separate the interstellar foreground polarization from the intrinsic polarization of HD 34664. The intrinsic polarization is consistent with electron scattering in a circumstellar disk seen at high inclination and provides further and independent evidence for the two-component wind model of the B(e) supergiants. We compare the polarization and spectrum of HD 34664 to observations published in the literature. In our data the continuum polarization is by a factor of two larger and the H-alpha line-to-continuum ratio is by a factor of at least four increased over previous observations. Since both the polarization and the line intensity are measures of the electron density in the stellar wind, we suggest that the observed variations can be attributed to a rise in the stellar mass-loss rate.

Schulte-Ladbeck, R. E.↗

ASCA Observations of "Type 2" LINERs Evidence for a Stellar Source of Ionization

We present ASCA observations of LINERs without broad H.alpha emission in their optical spectra. The sample of "type 2" LINERs consists of NGC 404, 4111, 4192, 4457, and 4569. We have detected X-ray emission from all the objects except for NGC 404; among the detected objects are two so-called transition objects (NGC 4192 and NGC 4569), which have been postulated to be composite nuclei having both an H II region and a LINER component. The images of NGC 4111 and NGC 4569 in the soft (0.5-2 keV) and hard (2-7 keV) X-ray bands are extended on scales of several kpc. The X-ray spectra of NGC 4111, NGC 4457 and NGC 4569 are well fitted by a two-component model that consists of soft thermal emission with kT approximately 0.65 keV and a hard component represented by a power law (photon index approximately 2) or by thermal bremsstrahlung emission (kT approximately several keV). The extended hard X-rays probably come from discrete sources, while the soft emission most likely originates from hot gas produced by active star formation in the host galaxy. We have found no clear evidence for the presence of active galactic nuclei (AGNs) in the sample. Using black hole masses estimated from host galaxy bulge luminosities, we obtain an upper limit on the implied Eddington ratios less than 5 x 10(exp -5). If an AGN component is the primary ionization source of the optical emission lines, then it must be heavily obscured with a column density significantly larger than 10(exp 23)/sq cm, since the observed X-ray luminosity is insufficient to drive the luminosities of the optical emission lines. Alternatively, the optical emission could be ionized by a population of exceptionally hot stars. This interpretation is consistent with the small [O I] lambda6300/H.alpha ratios observed in these sources, the ultraviolet spectral characteristics in the cases where such information exists, and the X-ray results reported here. We also analyze the X-ray properties of NGC 4117, a low-luminosity Seyfert 2 galaxy serendipitously observed in the field of NGC 4111.

Terashima, Yuichi↗

Connecting GRBs from binary neutron star mergers to nuclear properties of neutron stars

The fate of the binary neutron star (NS) merger remnants hinges sensitively upon the NS equation of state and the threshold mass, M ls , that separates a long-lived from a short-lived NS remnant. The nature of the electromagnetic counterparts is also influenced by the remnant type, particularly in determining whether a gamma-ray burst from a compact binary merger (cbGRB) is of short or long duration. We propose a novel approach to probe Mls by linking it to the estimated observed ratio of long to short cbGRBs. We find that current observations broadly favor a relatively high value for this transition, M ls ≃ 1.3M TOV , for which M TOV ≲ 2.6M ⊙ , consistent with numerical simulations, as also shown here. Furthermore, our results disfavor nuclear physics scenarios that would lead to catastrophic pressure loss at a few times nuclear density and temperatures of tens of MeV, leading to a rapid gravitational collapse of binaries with total mass M ≲ 1.3M TOV . Future individual gravitational wave events with on-axis cbGRBs can further bound Mls.

79 ASTRONOMY AND ASTROPHYSICS↗

The Precipitation Imaging Package: Assessment of Microphysical and Bulk Characteristics of Snow

Remote-sensing observations are needed to estimate the regional and global impacts of snow. However, to retrieve accurate estimates of snow mass and rate, these observations require augmentation through additional information and assumptions about hydrometeor properties. The Precipitation Imaging Package (PIP) provides information about precipitation characteristics and can be utilized to improve estimates of snowfall rate and accumulation. Here, the goal is to demonstrate the quality and utility of two higher-order PIP-derived products: liquid water equivalent snow rate and an approximation of volume-weighted density called equivalent density. Accuracy of the PIP snow rate and equivalent density is obtained through intercomparison with established retrieval methods and through evaluation with colocated ground-based observations. The results confirm the ability of the PIP-derived products to quantify properties of snow rate and equivalent density, and demonstrate that the PIP produces physically realistic snow characteristics. When compared to the National Weather Service (NWS) snow field measurements of six-hourly accumulation, the PIP-derived accumulations were biased only +2.48% higher. Additionally, this work illustrates fundamentally different microphysical and bulk features of low and high snow-to-liquid ratio events, through assessment of observed particle size distributions, retrieved mass coefficients, and bulk properties. Importantly, this research establishes the role that PIP observations and higher-order products can serve for constraining microphysical assumptions in ground-based and spaceborne remotely sensed snowfall retrievals.

Claire Pettersen↗

Investigation of the Solution Chemistry of Hybrid Organic–Inorganic Indium Halides for New Material Discovery

Recently, metal halide perovskites (MHPs) have emerged as a new class of materials for optical and electronic applications such as solar cells and ionizing radiation detectors. Although the solution-processability of MHPs is among their greatest advantages, the solution chemistries of most metal halide systems and their relationship with the observed structural and chemical diversity are poorly understood. In this work, we study the solution chemistry of a model indium halide system, methylammonium (MA)–In–Br, using a combination of the UV–vis spectroscopy, electrospray ionization mass spectrometry (ESI-MS) measurements, small-angle X-ray scattering (SAXS), and density functional theory (DFT) calculations. Our results show that indium could form either octahedral [InBr 6 3– ] or tetrahedral [InBr 4 – ] anions in solution or a combination of both, depending on the loading ratios of MABr and InBr 3 reactants. Understanding the solution chemistry of this system and recognizing the optical fingerprints of these polyanions allow for targeted crystallization of two novel compounds: MAInBr 4 featuring tetrahedral [InBr 4 – ] anions and MA 2 InBr 5 containing both octahedral [InBr 6 3– ] and tetrahedral [InBr 4 – ] anions. Further increase of the MABr content leads to the formation of previously reported MA 4 InBr 7 , containing only octahedral [InBr 6 3– ] anions separated by Br – anions. Furthermore, our results suggest that understanding the solution chemistry of multinary metal halide systems could be a valuable tool for discovering functional materials for practical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Low Cost Glass-Ceramic Matrix Composite Heat Exchanger

As part of ARPA-E’s High Intensity Thermal Exchange through Materials and Manufacturing Processes (HITEMMP) program, this project sought to develop novel heat exchanger (HX) capabilities to enable efficient and power dense power generation cycles. This class of HX comes under the category of ceramic/composite materials with the higher temperature goal in the program of ≥1100 °C inlet temperature operation. The enabling capability of this effort is the use of glass-ceramic matrix composite (GCMC) material which provides the high temperature durability of a ceramic, the flaw tolerance of a composite, a significantly faster and lower cost manufacturing process than conventional matrix CMCs and very low porosity levels < 0.5%. For thin-walled HX structures and the need to minimize leakage, the low porosity differentiator is particularly important. RTRC has prior experience with this material system and in the current project advanced the component design and manufacturing methods into new territory to produce features required for effective heat exchange under high pressures. In this approach, silicon carbide fiber is fabricated into a fiber preform using various textile processes. Graphite tooling is used both during the build-up of the fiber preform (interior tooling) and after the fiber preform has been completed (exterior tooling). This tooling assembly is heated to high temperature in an environment that has been evacuated and backfilled with inert gas. A reservoir of specialty glass is present and once the desired temperature has been reached to achieve the desired glass viscosity, an actuator distributes the glass throughout the fiber preform using passageways which are part of the tooling design in a process known as glass transfer molding. After the tooling has been removed, the composite is heat treated to convert the amorphous glass to a crystalline ceramic, providing improved properties. The project was divided into three phases focusing on the following: 1) 10 kW HX design and coupon-level tube sheet fabrication, 2) 10 kW HX fabrication, 3) 50 kW HX fabrication. During Budget Period 1 (BP1), additional risks were encountered and the need for additional funds was agreed upon by ARPA-E program leadership. Due to a variety of factors, the contract modification required nominally 18 months to execute at which time the HITEMMP program was effectively concluding. Because of this and the time that would be required to perform BP2 tasks, it was decided to conclude the project at the end of BP1. During the design of the 10 kW HX, manufacturing constraints were learned and incorporated, leading to a revised configuration for the fiber preform and HX. Heat exchange and pressure drop predictions also played a role in modifying the original design concept to be a higher aspect ratio shell-and-tube HX, simplifying the manufacturing process and improving the heat exchanger performance. Good gravimetric and volumetric thermal power densities of 11.2 kW/kg and 10,200 kW/m3 for the entire HX were projected that involved thermo-structural Finite Element Analysis to determine the structural mass needed for the high operation pressures of 250 bar cold inlet and 80 bar hot inlet. Fiber preforms using textile processes were produced for multiple headered tube sheets. Additional challenges were encountered during the glass transfer molding step for which solutions were identified, but programmatics did not allow them to be implemented in BP1. While complete HX test articles were not fabricated, the benefits of this GCMC material for a variety of high temperature applications remain.

30 DIRECT ENERGY CONVERSION↗

A new determination of the halo luminosity density of the galaxy

The distribution of tangential velocities for arbitrary halo velocity dispersions and rotation rates are explicitly computed, assuming a particular form of the halo velocity distribution function near the sun. It is found that for a specific model of the halo velocity dispersion tensor that Schmidt's (1975) estimate of the mass carried by halo stars of 0.25-0.70 solar masses should be revised by a factor of between 1.0 and 2.0, depending on the rotation rate of the halo. Although obtaining a mass density in all stars from the stellar density in that mass range is extremely uncertain, it is found, in agreement with earlier work, that even under rather extreme assumptions about the remnants the observed stars and expected remnants fail to provide the theoretically expected mass. Comparison of the luminosity to the theoretical mass yields a mass to light ratio consistent with that of rich clusters of galaxies.

Richstone, D. O.↗

Type II bursts, shock waves, and coronal transients - The event of 1980 June 29, 0233 UT

The metric Type-II solar burst event of June 29, 1980, is characterized on the basis of spatially resolved radioheliograph observations obtained at Culgoora, Australia, and visible-light observations obtained with the coronograph/polarimeter of the SMM satellite. The data are presented in images, diagrams, and graphs and discussed in detail. The Type-II emission is found to arise in the dense moving material behind the transient loops, which have sky-plane width 0.5 solar radius and line-of-sight depth 0.1-0.4 solar radius. A faint arc observed moving ahead of the transient loops at about 900 km/sec and not associated with the Type-II burst is attributed to a shock front, and the compression ratio and Alfven Mach number of the enhanced-density region are estimated as n2/n1 = 1.3-3 and M(A) = 1.2-3. The ambient material at 3 solar radii is determined to have Alfven speed 250-625 km/sec and magnetic-field strength 50-120 mG. The total mass of the event is calculated as 700 Tg; the total magnetic energy of the loops is (1.5-15) x 10 to the 29th ergs.

Gary, D. E.↗

In defense of Ertel's potential vorticity and its general applicability as a meteorological tracer

A potential vorticity theorem and its two summary statements published by Haynes and McIntyre are challenged conceptually by equations, discussions and examples. The apparent simplification proposed by the authors to convert from a mass to volume integral, i.e., by cancelling density against the specific volume in the potential vorticity, changes the physical significance of the integrand. It no longer is the potential vorticity. The resulting mean for either a bulk Eulerian or Lagrangian system is then not analogous to a mixing ratio and therefore not independent of the broad spectrum of internal waves, the independence that makes Ertel's potential vorticity so valuable either as a stratospheric tracer or as a predictive or diagnostic, large scale, meteorological variable.

Danielsen, Edwin F.↗

Parameterizations of Chromospheric Condensations in dG and dMe Model Flare Atmospheres

The origin of the near-ultraviolet and optical continuum radiation in flares is critical for understanding particle acceleration and impulsive heating in stellar atmospheres. Radiative-hydrodynamic (RHD) simulations in 1D have shown that high energy deposition rates from electron beams produce two flaring layers at T approximately 10 (exp 4) K that develop in the chromosphere: a cooling condensation (downflowing compression) and heated non-moving (stationary) flare layers just below the condensation. These atmospheres reproduce several observed phenomena in flare spectra, such as the red-wing asymmetry of the emission lines in solar flares and a small Balmer jump ratio in M dwarf flares. The high beam flux simulations are computationally expensive in 1D, and the (human) timescales for completing NLTE models with adaptive grids in 3D will likely be unwieldy for some time to come. We have developed a prescription for predicting the approximate evolved states, continuum optical depth, and emergent continuum flux spectra of RHD model flare atmospheres. These approximate prescriptions are based on an important atmospheric parameter: the column mass (m(sub ref)) at which hydrogen becomes nearly completely ionized at the depths that are approximately in steady state with the electron beam heating. Using this new modeling approach, we find that high energy flux density (>F11) electron beams are needed to reproduce the brightest observed continuum intensity in IRIS data of the 2014 March 29 X1 solar flare, and that variation in m(sub ref) from 0.001 to 0.02 g cm (exp -2) reproduces most of the observed range of the optical continuum flux ratios at the peak of M dwarf flares.

Kowalski, Adam F.↗

Predicting the rheology of limestone calcined clay cements (LC 3 ): Linking composition and hydration kinetics to yield stress through Machine Learning

The physicochemical characteristics of calcined clay influence yield stress of limestone calcined clay cements (LC 3 ), but the independent influences the clay's physical and chemical characteristics as well as the effect of other variables on LC 3 rheology are less well-understood. Further, a relationship between LC 3 hydration kinetics and yield stress – important for informing mixture design – has not yet been established. Here, rheological properties were determined in pastes with varying water-to-solid ratio (w/s), constituent mass ratios (PC:metakaolin:limestone), limestone particle size and gypsum content. From these data, an ML model developed allowed the independent examination of the different mechanisms by which metakaolin fraction influences yield stress of LC 3 , identifying four predictors – packing index, Al 2 O 3 /SO 3 , total particle density and metakaolin fraction relative to limestone (MK/LS) – most significant for predicting LC 3 yield stress. A methodology based on kernel smoothing also identified hydration kinetics parameters best correlated with yield stress.

36 MATERIALS SCIENCE↗

Fuel cells for single-aisle regional aircraft: System configuration, performance and cost

A hydrogen fuel cell propelled electric aircraft can compete with incumbent turbofan technologies for single-aisle regional aircraft by coupling design of stack, air handling, thermal management, propulsion, and airframe to optimize performance. The stack operates at 95°C to facilitate heat rejection during take-off and below 75°C during cruise to extend lifetime and is oversized to satisfy power requirements at end of life. A multi-stage turbocompressor with a compression ratio >10 is selected to reach high stack power density at 11,300-m cruise altitude. The propulsion system is configured to accommodate air handling within the core duct, an inclined heat exchanger in the outer duct to limit the nacelle size, and variable area nozzles to independently control mass flows through the core and bypass ducts. The airframe is modified for maximum lift coefficient and longer balanced field length for dramatically reduced thrust during take-off, and the fuselage is stretched by 20% to store liquid hydrogen (LH 2 ). Modularization of power systems promotes safety in one engine inoperative scenarios and allows reaching specific power metrics for stack, balance-of-plant and fuel cell system (FCS), necessary for acceptable take-off weight. In conclusion, cost parity requires increase in FCS lifetime, LH 2 cost reduction, and improved FCS specific power.

Catalyst durability↗

Shock wave relations in lunar ash flow

A detailed analysis of steady normal shock waves in a layer of lunar ash by the theory of two phase flow of a mixture of a gas and small solid particles is presented. New terms of pressure gradient and virtual mass forces in the particle momentum equation and the particle kinetic energy and work done in the particle energy equation are included in the fundamental equations and their influences are investigated. The flow variables as functions of the free stream Mach number, initial particle volume fraction and the density ratio G are presented. The thickness of the relaxation zone is found to increase with decreasing initial particle volume fraction and almost independent of G for any given values of free stream Mach number. A new empirical relation of the drag coefficient for the spherical particles in the mixture is proposed to fit experimental data.

Pai, S. I.↗

Generation and comparison of globally isotropic space-filling truss structures

The purposes of this paper are to present a rationale for obtaining space-filling truss structures that behave like a globally isotropic continuum and to use continuum modeling to investigate their relative structural efficiencies (e.g., modulus-to-density, strength-to-density, and part-count-to-volume ratios). The trusses considered herein are generated by replication of a characteristic truss cell uniformly through space. The characteristic cells are categorized by one of a set of possible geometric symmetry groups derived using the techniques of crystallography. The implied elastic symmetry associated with each geometric symmetry group is identified to simplify the task of determining stiffness tailoring rules for guaranteeing global isotropy. Four truss geometries are analyzed to determine stiffness tailoring necessary for isotropy. All geometries exhibit equivalent isotropic Poisson's ratios of 1/4 and equivalent modulus-to-density ratios of 1/6 times the modulus-to-density ratio of the material used in their members. The truss configuration that has the lowest percent difference in member lengths is shown to have the lowest component part-count-to-volume ratios of all geometries considered when compared on a basis of equal stiffness, equal strength, and equal mass.

Lake, Mark S.↗

Chemical and Physical Properties of Bulk Aerosols within Four Sectors Observed during TRACE-P

Chemical and physical aerosol data collected on the DC-8 during TRACE-P were grouped into four sectors based on back trajectories. The four sectors represent long-range transport from the west (WSW), regional circulation over the western Pacific and Southeast Asia (SE Asia), polluted transport from Northern Asia with substantial sea salt at low altitudes (NNW) and a substantial amount of dust (Channel). WSW has generally low mixing ratios at both middle and high altitudes, with the bulk of the aerosol mass due to non-sea-salt water-soluble inorganic species. Low altitude SE Asia also has low mean mixing ratios in general, with the majority of the aerosol mass comprised of non-sea-salts, however, soot is also relatively important m this region. "w had the highest mean sea salt mixing ratios, with the aerosol mass at low altitudes (a km) evenly divided between sea salts, mm-sea-salts, and dust. The highest mean mixing ratios of water-soluble ions and soot were observed at the lowest altitudes (a km) in the Channel sector. The bulk of the aerosol mass exported from Asia emanates h m Channel at both low and midaltitudes, due to the prevalence of dust compared to other sectors. Number densities show enhanced fine particles for Channel and NNW, while their volume distributions are enhanced due to sea salt and dust Low-altitude Channel exhibits the highest condensation nuclei ((34) number densities along with enhanced scattering coefficients, compared to the other sectors. At midaltitudes (2-7 km), low mean CN number densities coupled with a high proportion of nonvolatile particles (265%) observed in polluted sectors (Channel and NNW) are attributed to wet scavenging which removes hygroscopic CN particles. Low single scatter albedo m SE Asia reflects enhanced soot

Jordan, C. E.↗

Can the Streaming Instability Form the First Planetesimals in Globally Turbulent Protoplanetary Disks?

Understanding the formation of the first planetesimals remains key to deciphering the history of planet formation within our own solar system and beyond. Evidence from the vast meteorite record [e.g., 1] as well as observations [e.g., 2] strongly suggest that the first planetesimals, and perhaps giant planet core accretion occurred well within the first million years of disk evolution. Moreover, chemical and lithological mixing as well as observations of line-broadening in protoplanetary disks (PPDs) [e.g., 3-6] suggest that the solar system nebula in this epoch was weakly-to-moderately turbulent in the regions where particle growth is of the greatest interest [7-8]. It is a well-known though that global hydrodynamic turbulence complicates particle growth due to a slew of barriers that can slow or even stall particle or aggregate growth at pebble sizes (with corresponding small particle Stokes numbers St) that can lead to loss to the central star via radial drift before planetesimals can ever form [9,10], requiring that some other mechanism come into play that collects growth-frustrated pebbles into gravitationally bound multi-km bodies – objects that are “born big” [11]. The current leading candidate for such a “leap-frog” mechanism is the so-called Streaming Instability (SI), a gas-drag mediated momentum exchange resonance in which the relative velocity between the particle component and a rotating gaseous fluid can lead to high densities in the particle field [12], which has been invoked in a number of recent PPD models that include a turbulent intensity α [e.g., 13-14] as the defacto mechanism for planetesimal formation if conditions for the SI (depending on particle St and the solids-to-gas mass ratio) are satisfied. However, these previous works use conditions established from occurrence studies for the onset of SI in laminar disks [15-17] in which the only source of turbulence is that self-generated by the settling particle layer, and not externally driven global turbulence. Recent analytical theories of the SI subject to global turbulence predict much more stringent conditions for the effectiveness of the SI than the laminar case [18-19]. Thus, whether the efficient operation of the SI can be attained in realistic models of the solar nebula have yet to be established. In this work, we ask whether the conditions under which the SI can produce gravitationally bound particle overdensities can actually be met in the first million years of evolution in globally turbulent PPDs.

P R Estrada↗

Five Years of Mid-Infrared Evolution of the Remnant of SN 1987A: The Encounter Between the Blast Wave and the Dusty Equatorial Ring

We have used the Spitzer satellite to monitor the laid-IR evolution of SN 1987A over a 5 year period spanning the epochs between days 6000 and 8000 since the explosion. The supernova (SN) has evolved into a supernova remnant (SNR) and its radiative output, is dominated by the interaction of the SN blast wave with the pre-existing equatorial ring (ER). The mid-IR spectrum is dominated by emission from approximately 180 K silicate dust, collisionally-heated by the hot X-ray emitting gas with a temperature and density of 5 x 10(exp 6) K and approximately 3 x 10(exp 4) per cubic centimeter, respectively. The mass of the radiating dust is approximately 1.2 x 10(exp -6) solar mass on day 7554, and scales linearly with IR flux. Comparison of the IR data with the soft X-ray flux derived from Chandra observations shows that the IR-to-X-ray flux ratio, IRX, is roughly constant with a value of 2.5. Gas-grain collisions therefore dominate the cooling of the shocked gas. The constancy of IRX is most consistent with the scenario that very little grain processing or gas cooling have occurred throughout this epoch. The shape of the dust spectrum remained unchanged during the observations while the total flux increased by a factor of approximately 5 with a time dependence of t(sup '0.87 plus or minus 0.20), t' being the time since the first encounter between the blast wave and the ER. These observations are consistent with the transitioning of the blast wave from free expansion to a Sedov phase as it propagates into the main body of the ER, as also suggested by X-ray observations. The constant spectral shape of they IR, emission provides strong constraints on the density and temperature of the shocked gas in which the interaction takes place. The IR spectra also suggest the presence of a secondary population of very small, hot (T greater than or equal to 350 K), featureless dust. If these grains spatially coexists with the silicates, then they must have shorter lifetimes. The data show slightly different rates of increase of their respective fluxes, lending some support to this hypothesis. However, the origin of this emission component and the exact nature of its relation to the silicate emission is still a major unsolved puzzle.

Dwek, Eli↗

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux $a_y$ is calculated. In addition to mostly negative values of $a_y$ , typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of $a_y$ due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ( $\sim$ 60 %) is observed.

Mechanics↗