Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BUBBLE”

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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A reacting multiphase computational flow model for 2,3-butanediol synthesis in industrial-scale bioreactors

The microbe, Zymomonas mobilis, can efficiently convert sugars to 2,3 butanediol (BDO), an important intermediate for downstream chemical products, only in a well-controlled microaerated environment. However, controlling oxygen distribution in industrial-scale bioreactors is challenging, and at-scale BDO production is hard to optimize using traditional engineering methods. This study takes a step towards addressing this problem through a computational model for reacting multiphase flows in large-scale bioreactors. A phenomenological metabolic model was first developed and validated against experiments, and then it was coupled to a multiphase computational fluid dynamics (CFD) solver using a subcycling algorithm for simulating long conversion times (~30 hours). Large-scale (~500 m 3 ) bubble column simulations using this coupled approach demonstrated a 25% improvement in BDO yield with low sparging rates (0.002 m/s) and low height-to-diameter ratio (0.875) compared to a baseline bubble column (35 m height, 5 m diameter) with a superficial gas velocity of 0.06 m/s.

09 BIOMASS FUELS↗
DOE OSTI2021

Unified descriptor for enhanced critical heat flux during pool boiling of hemi-wicking surfaces

Boiling heat transfer is dictated by interfacial phenomena at the three-phase contact line where vapor bubbles form on the surface. Structured surfaces have shown significant enhancement in critical heat flux (CHF) during pool boiling by tailoring interfacial phenomena. This CHF enhancement has been primarily explained by two structural effects: roughness, which extends the contact line length at the bubble base, and wickability, the ability to imbibe liquid through surface structures by capillary pumping. Here, in this work, we show that CHF enhancement on structured surfaces cannot be described by roughness or wickability alone. This result was confirmed using systematically designed micropillar surfaces with controlled roughness and wickability. Further, we performed a scaling analysis and derived a unified descriptor, which represents the combined effects of thin film density and volumetric wicking rate. This unified descriptor shows a reasonable correlation with CHF values with our experiments and literature data. This work provides important insights in understanding the role of surface structures on CHF enhancement, thereby providing guidelines for the systematic design of surface structures for enhanced pool boiling heat transfer.

42 ENGINEERING↗
DOE OSTI2020

Deep learning interfacial momentum closures in coarse-mesh CFD two-phase flow simulation using validation data

Multiphase flow phenomena have been widely observed in the industrial applications while it remains a challenging yet unsolved problems. Three-dimensional computational fluid dynamics (CFD) approaches resolve the flow fields on a finer special and temporal scales which can complement the dedicated experimental study. However, closures have to be introduced to reflect the underlying physics in multiphase flow. Among them, the interfacial forces, including drag, lift, turbulent dispersion and wall lubrication forces, play in important role on the bubble’s distribution and migration in liquid-vapor two-phase flow. Development of those closures traditionally rely on the experimental data and analytical derivation with simplified assumptions which usually cannot deliver a universal solution across wide range of flow conditions. In this paper, a data-driven approach, named as Feature Similarity Measurement (FSM), is developed and applied to improve the simulation capability of two-phase flow with coarse-mesh CFD approach. Interfacial momentum transfer in adiabatic bubbly flow serves as the focus of the present study. Both a mature and a simplified set of interfacial closures are taken as the low fidelity data. Experimental data and fine mesh CFD simulations results are adopted as high-fidelity data. Qualitative and quantitative analysis are performed in this paper which reveals that FSM can substantially improve the prediction of coarse mesh CFD model regardless of the choice of interfacial closures and it provides scalability and consistency across discontinuous flow regimes. Furthermore, it demonstrates that data-driven method can aid the multiphase flow modeling by exploring the connections between local physical features and simulation errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Probing thermodynamics of radiogenic helium and defects in $δ$-plutonium alloys and interactions with adsorbed environmental gases

Differential scanning calorimetry coupled with simultaneous evolved gas analysis (DSC-EGA) on aged δ-Pu samples shows that most radiogenic helium remains trapped within the Pu matrix at temperatures very close to, or slightly above, the melting temperature. Our results indicate that helium release from 50-year-old δ-Pu occurs as a burst just below the melting temperature (>0.994 T m ), with subsequent pressure oscillations as temperature increases. Subordinate quantities of H 2 were also released along with helium. Here, the helium emission tails off and ceases above ~720 – 750°C. In a δ-Pu alloy aged 6 years, the initial helium burst occurs slightly above melting (~1.015 – 1.042 T m ), with a discrete, larger helium spike occurring between 670 and 686°C. The proximity of helium release to the liquidus transition presented challenges in the deconvolution of overlapping process enthalpies, the liquidus endotherm, and the exotherm resulting from bubble collapse, annealing and gas expulsion. Helium’s strong affinity for vacancy binding in a 2He-vac configuration is predicted by Density Functional Theory (DFT) modeling. The measured stored energy associated with the He release events in a 50-year-old δ alloy is on the order of ~ 10–11 J/g, which is significantly higher than stored energies measured in the sub-solidus regimes (~2 J/g) that are related to the solid-state annealing of processing- and radiation-induced defects. This implies that aged δ Pu alloys have a remarkable resilience to accommodate the lattice strain produced by the internal pressure of the helium bubbles and provides further insight into the thermodynamic behavior of aged δ Pu.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Effect of cold-cap-bottom shear stress on primary foam stability

Here, we analyze the heat transfer between the melt pool and the layer of reacting materials called cold cap, which is the determining factor for the glass production rate in Joule-heated air-bubbled nuclear waste glass melters. We consider the effects of bubbling rate, melt viscosity, melt thermal conductivity, conversion kinetics, and shear rate at the cold-cap bottom, and develop simple relationships for the cold-cap bottom temperature and glass production rate. We show that the melt convection enhances the heat transfer to the cold cap by suppressing the thermal boundary layer and primary foam, while the conversion kinetics opposes these effects. Additionally, we mention how dissolving silica particles affects the primary foam stability and discuss the limitations of the presented approach for vigorously foaming feeds.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The use of U3Si2/Al dispersion fuel for high power research reactors

In order to investigate the irradiation performance of U3Si2/Al dispersion fuel at relatively aggressive reactor conditions, compared to how this fuel is typically used in a research or test reactor, post irradiation examination was performed on two fuel plates that were irradiated at ~270 W/cm2 average surface heat flux and to a maximum local burnup greater than 8 × 1021 fissions/cm3 in the Advanced Test Reactor as part of the RERTR-8 experiment. As part of the non-destructive examinations that were performed on both fuel plates, detailed thickness measurements using a high-fidelity measurement bench were performed on one fuel plate and a segment of the second fuel plate. This paper provides the first reported results of such high-resolution thickness measurements for U3Si2/Al dispersion fuel plates, containing primarily the U3Si2 fuel phase, irradiated at such power levels to high burnup. For destructive examination using optical metallography, samples were generated from both irradiated plates. Localized regions of high swelling (characteristic of internal blistering) were observed on the surface of both fuel plates, suggesting that the performance limits had been reached for each plate. Destructive examinations showed that the fuel particles that are primarily U3Si2, which become amorphous during irradiation, in both plates contained relatively large fission gas bubbles with evidence of fission gas bubble interconnection.

Wachs, Daniel M.↗
DOE OSTI2020

Microstructure evolution in U-10Zr alloy irradiated by swift Xe ions at 700 °C

The post-irradiation examination (PIE) results of U-10Zr binary metallic fuel that was irradiated by 84 MeV Xe ions at 700 °C are reported here. The U-10Zr specimen was synthesized by homogenizing elemental U and Zr powders through arc-melting. The as-fabricated U-10Zr was dominated by the sub-micron acicular α+ δ phase. After being irradiated by 84 MeV Xe ions at 700 °C up to 2.19 ×10 17 ions/cm 2 fluence, the U-10Zr was speculated to evolve into micro-scale γ 1 + γ 2 phases as predicted by the U-Zr phase diagram. During post-irradiation cooling, the U-rich γ 1 phase was mainly transformed to the αphase, while the Zr-rich γ 2 phase was partially retained along with formation of hexagonal structure phase (δ or ω). Xe bubbles formed in the γ 1 and γ 2 phases have salient different morphologies, implying faster Xe diffusion in the γ 1 phase than the γ 2 phase. The quantitative bubble morphology data collected in this study can provide valuable references for the development and validation of fuel performance models for this fuel.

36 MATERIALS SCIENCE↗
DOE OSTI2020

α-U and ω-UZr 2 in neutron irradiated U-10Zr annular metallic fuel

Here, to develop metallic fuel with ultra-high burnup of 30%-40%, an annular U-10Zr fuel with 55% smear density was fabricated through a casting route and irradiated at the Advanced Test Reactor at Idaho National Laboratory. The annular fuel design also serves as a demonstration of the feasibility to replace sodium bond with a helium bond to benefit the geological disposal of irradiated fuel, cut the cost of fuel fabrication, and boost the overall metallic fuel economy. This paper reports the results from transmission electron microscopy (TEM) based post-irradiation examination of this fuel type irradiated to a burnup of 3.3% fissions per initial heave metal atoms for initial screening purpose. After irradiation, the initial U-10Zr separated into an a-U annular region and an UZr 2+x center region with nanoscale spinodal decomposed microstructure. Because of the provided large amount of coherent interface areas, the fission gas atoms and vacancies generated in UZr 2+x phase are possibly pinned at the interface areas, leading to 20 times smaller fission gas bubbles than those in the neighboring a-U. The large bubbles in a-U become connected and merged into large pores that provide fast release paths for fission gas which prevents further fuel swelling. The fuel center still has open space to accommodate further fuel swelling from solid fission products at higher burnup. Other neutron irradiation induced phase and microstructure change are also characterized and compared with traditional solid fuel designs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Radiation induced hardening of beryllium during low temperature He implantation

The effect of ion irradiation on evolution of microstructure and hardening of beryllium with different impurity levels was investigated using TEM and nanoindentation. High purity S-65 grade and less-pure S-200-F grade were implanted by helium ions at temperatures of 50°C and 200°C. 11 different energies were used, so as to create a quasi-homogeneous 3 µm irradiated layer with average radiation damage of 0.1 dpa and average He content of 2000 appm.Nanoindentation experiments demonstrated that before irradiation, the S-200-F and S-65 grades have an average hardness of 3.7$\pm$0.8 GPa and 3.4$\pm$0.8 GPa correspondently. After implantation the hardness of both grades increased by about 60% for the 200°C irradiation and 100% for the 50°C irradiation. The crystallographic analysis of indented grains demonstrated that in the as-received materials the hardness is about 2.5 times higher when the indentation direction is close to the [0001] c-axis of beryllium compared to indentation perpendicular to [0001]. Hardness anisotropy significantly decreased after irradiation: the “soft orientation” was most sensitive to irradiation-induced hardening, with hardness increasing by about 140% after irradiation at 50°C and 100% after irradiation at 200°C, compared to about 15 - 20% for the “hard” orientation at both irradiation temperatures. The higher purity grade had smaller increase of the “soft orientation” hardness: 2.5$\pm$0.3 GPa for the S-65 and 2.9$\pm$0.2 GPa for the S-200-F.At both temperatures in both grades, under TEM investigation the radiation damage appears as “black dots” which are likely to be small dislocation loops with the number density of ~ 10 22 m -3 . No bubbles were observed by TEM inside grains and at grain boundaries. Analysis of the possible hardening contribution demonstrated that the observed “black dots” could be responsible for up to half of the measured hardening, while the rest of the hardening should originate from helium bubbles with the size below the TEM resolution (at or below 1.5 nm).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational determination of a primary diffusion mode in γ U-10Mo under irradiation

Low enriched uranium (< 20 % 235 U)-molybdenum (U-Mo) monolithic fuel is the primary candidate for high performance research and test reactors, and is in the process of being qualified to replace highly enriched uranium (≥ 20 % 23 5U) fuel. As part of the qualification process, it is critical to understand and predict the behavior of fission gas bubbles under irradiation, which affects fuel swelling and fuel failure. Mechanistic fuel models are being developed that can both reproduce the existing experimental data for fuel swelling, and be further applied to irradiation conditions beyond the experimental scope. Diffusion of species under irradiation conditions is an important parameter in the mechanistic fuel models; however, no temperature-relevant experimental diffusion data exists. In the present work, radiation-enhanced diffusion coefficients of U, Mo, and Xe in λU-10wt.%Mo were calculated in the temperature range between 300 K and 1400 K via rate-theory models and molecular dynamics simulations with an embedded-atom method interatomic potential for the U-Mo-Xe system. Accordingly, total diffusion coefficients under relevant irradiation conditions are determined using previously obtained intrinsic thermal diffusion and radiation-driven diffusion coefficients, as well as the newly calculated radiation-enhanced diffusion coefficients presented herein. Radiation-enhanced diffusion of U and Mo was dominant in the intermediate temperature range, whereas radiation-enhanced diffusion of Xe did not significantly contribute to total diffusion of Xe at the relevant fission rate densities. Radiation-enhanced diffusion of Xe became faster than both intrinsic thermal diffusion and radiation-driven diffusion at a fission rate density of 5 x 10 22 fissions/m 3 /s, which is higher than the typical fission rate density range in research reactors. The temperature regime that radiation-enhanced diffusion of each element dominated was dependent on the fission rate density. Finally, the total diffusion coefficients of U, Mo, and Xe, updated in this work, will be utilized as parameters in the mechanistic fuel models to help predict the behavior of fission gas bubbles under irradiation more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Integrated simulation of U-10Mo monolithic fuel swelling behavior

Here, a separate computational branch has been implemented within the DART (Dispersion Analysis Research Tool) computational code to simulate the swelling behavior of U-10Mo monolithic fuel under the operating conditions of high-power research and test reactors (RTRs). The monolithic branch of the DART code implements a mechanistic rate-theory-based fission-gas-behavior model for the calculation of fission gas swelling, as well as a suite of thermal, physical, and mechanical models to take into account various processes occurring in RTR fuels during irradiation. In order to accurately simulate and eventually predict U-10Mo monolithic fuel irradiation behavior, the code uses materials properties calculated with lower length-scale computational methods, such as gas atom diffusivity and U-Mo surface energy from atomic simulations and grain-morphology-specific recrystallization kinetics (recrystallized fuel volume fractions vs. fission density) predicted using the phase-field method. The remainder of fission gas behavior parameters used in the model were calibrated with measured intergranular bubble size distributions. With this integrated simulation approach, the swelling behavior of U-10Mo monolithic fuel was simulated for various initial grain sizes at different operating conditions and compared with measured data. Furthermore, because limited experimental data exist for parameter calibration detailed sensitivity studies for the important parameters used in the fission gas behavior model were performed in order to examine their impact on both intergranular gas bubble morphology at low fission density, and on total porosity at high fission density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Microstructural evolution of compositionally complex solid-solution alloys under in-situ dual-beam irradiation

Here, this work attempts to link the microstructural evolution of single-phase compositionally complex alloy (CCA) compositions under dual-beam irradiation to their Mn-content via the stacking-fault energy (SFE) and vacancy migration energies. Two alloys, Cr 18 Fe 27 Mn 27 Ni 28 and Cr 15 Fe 35 Mn 15 Ni 35 , along with less compositionally complex pure Ni and Fe 56 Ni 44 binary were irradiated at 500 and 600 °C under dual-beam 1 MeV Kr 2+ and 16 keV He + heavy-ions up to 7 displacements per atom (dpa) with a He/dpa ratio of 0.75 %/dpa using in situ transmission electron microscopy (TEM). Due to the bubble-stabilizing effect of implanted He, bubbles were observed in all irradiations, and populations of faulted interstitial loops were characterized in Cr 18 Fe 27 Mn 27 Ni 28 and Cr 15 Fe 35 Mn 15 Ni 35 . A reduction in swelling was observed in the two CCAs compared to pure Ni and Fe 56 Ni 44 . Although swelling increased from 500 to 600 °C in Fe 56 Ni 44 , Cr 15 Fe 35 Mn 15 Ni 35 and Cr 18 Fe 27 Mn 27 Ni 28 both swelled slightly more at 500 °C. This was attributed to the difference in vacancy mobility, stronger pinning effect of vacancies on He, and the sink strength of faulted dislocation loops. Faulted interstitial loops nucleated with a higher number density and dislocation line density in Cr 15 Fe 35 Mn 15 Ni 35 at both temperatures, and at 600 °C in both materials. The differences in faulted loop population and temperature effect on swelling are correlated to the Mn-content and the measured SFE (20.2 ± 6.7 mJ/m 2 for Cr 18 Fe 27 Mn 27 Ni 28 and 9.2 ± 3.4 mJ/m 2 for Cr 15 Fe 35 Mn 15 Ni 35 ).

36 MATERIALS SCIENCE↗
DOE OSTI2023

An integrated statistical-thermodynamic model for fission gas release and swelling in nuclear fuels

Here, we propose a new model for burst fission gas release induced by microcracking in ceramic nuclear fuels such as uranium dioxide. The model stipulates that the densities of defects in the fuel material, such as microcracks and fission gas bubbles on grain boundaries, evolve in accordance with the second law of thermodynamics. Central to the model is the notion of an effective temperature, conjugate to the configurational entropy of the fuel material, and directly linked to the burnup. The model predicts that microcracking, driven by the internal stress state of the fuel material, reduces the bubble storage capacity of grain boundaries, and accounts for burst fission gas release during rapid temperature transients that simulate power transients, reactor startup, and loss-of-coolant accident conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Mesoscale Modeling of the Effects of Accelerated Burnup on UO2 Microstructural Evolution

Accelerating the nuclear fuel qualification process will rely on some combination of advanced modeling and simulation techniques with accelerated irradiation testing and separate effects experiments to enable the development of new fuel concepts in a shorter time frame. One of the key challenges to successfully leveraging accelerated irradiation tests will be understanding the artifacts that may be introduced with accelerated accumulation of dose and/or burnup. This work presents phase field (MARMOT) simulations of the evolution of representative 2D UO2 microstructures up to 40 MWd/kgU. Simulations were performed under both commercial light water reactor fuel conditions as well as those that would be expected for highly accelerated (~10x) burnup conditions similar to those used in the MiniFuel irradiations in Oak Ridge National Laboratory’s High Flux Isotope Reactor. The phase field model was coupled with a discrete nucleation algorithm to model re- structuring at high burnup. The effect of the different fission rates in both microstructures was investigated at two temperatures: 650?C and 800?C. The lower temperature simulations both showed an onset of restructuring at nearly 60 MWd/kgU. More extensive restructuring was obtained in the MiniFuel microstructure compared with that of the PWR fuel. At 800?C, no restructuring was obtained as a result of the thermally activated diffusion of Xe atoms and U vacancies to fission gas bubbles, which reduces the nucleation driving force. These results highlight the importance of using modeling and simulation tools to inform the environmental conditions during targeted accelerated irradiation tests to extract the most useful fuel performance data.

accelerated fuel qualification, Phase Field, Restr↗
DOE OSTI2022

Grain subdivision and structural modifications by high-energy heavy ions in UO 2 with different initial grain size

Irradiation with 84 MeV Xe ions at peak fluence above 7 × 10 17 ions/cm 2 has been carried out on two UO 2 samples with different initial grain sizes at 300 °C and 500 °C. Characterization of both grain structure and bubble morphology has been performed using transmission electron microscopy, image analysis and precession electron diffraction. The bubble size in the large-grained sample was smaller than that in the standard-grain sample, confirming fission gas precipitation resistance in the former. The ion irradiation in the high electronic energy deposition range has caused formation of a large density of dislocations which resulted in grain polygonization for the sample with standard grain size, which was observed in multiple lamellas extracted next to each other at peak dose. The sample with large grains showed inconsistent behavior. In one case, dislocation pile-up occurred, but no sub-grain formation was observed. An adjacent location showed sub-grain formation, even if to less extent compared to the standard grain-size sample. Furthermore, these findings suggest a limited polygonization resistance in large-grained UO 2 compared to standard-grain UO 2 . Recrystallization of both sample surfaces was present, leaving unclear if this is a surface effect or a consequence of the highest energy deposition rate in this area.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Simulation scaling studies of reactor core two-phase flow using direct numerical simulation

Tremendous growth in supercomputing power in recent years has resulted in the emergence of high-resolution flow analysis methods as an advanced research tool to evaluate single and two-phase flow behavior. In particular, unstructured mesh-based methods have been applied to analyze flows in complex reactor core geometries, including those of light water reactors (LWR). The finite-element based code, PHASTA, is utilized to perform large-scale simulations of two-phase bubbly flows in LWR geometries. Given the large computational cost of direct numerical simulation (DNS) coupled with interface tracking methods (ITM), typical domains encompass a portion of a single subchannel. In the presented research, the state-of-the-art analysis of turbulent two-phase flows in complex LWR subchannel geometries are demonstrated at both prototypical reactor parameters as well as scaled low pressure conditions. Three different cases are studied, a high-pressure simulation in prototypical reactor subchannel geometry, a low-pressure case in prototypical geometry and a final low-pressure case in a geometry scaled up to conserve the ratio between the bubble size and the domain pitch. Utilizing advanced statistical processing tools, these simulation conditions are compared to shed light on the relevancy of two-phase flow characteristics given the significant differences between LWR and low-pressure conditions. These findings can lead to the generation of useful guiding principles when researchers need to scale the two-phase flow behavior captured at low pressure and temperature conditions to those at reactor operating conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗