Search NASA⌕ Search

SEARCH · Search NASA

Results for “BHR”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗
NASA NTRS1991

Fluid seals technology at BHR Group Ltd.

A mathematical model was developed for ESA, which couples the rigid body dynamics in four degrees of freedom of a flexibly mounted stator (with Coulomb damping) of a turbopump face seal with the fluid film behavior. The interfacial film is assumed gaseous, and the model incorporates laminar/turbulent flow, inertia and sonic choking. Seal faces may be wavy and coned, while the rotating ring can be eccentrically mounted, misaligned with the shaft and oscillate axially at arbitrary frequency. The model is currently being coded into software. An experimental program using a high speed cryogenic/hot gas test rig is being conducted to verify the output of the computer program.

Leefe, Simon↗
DOE OSTI2024

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hemerythrins enhance aerobic respiration in Methylomicrobium alcaliphilum 20ZR, a methane-consuming bacterium

ABSTRACT Numerous hemerythrins, di-iron proteins, have been identified in prokaryote genomes, but in most cases their function remains elusive. Bacterial hemerythrin homologs (bacteriohemerythrins, Bhrs) may contribute to various cellular functions, including oxygen sensing, metal binding and antibiotic resistance. It has been proposed that methanotrophic Bhrs support methane oxidation by supplying oxygen to a core enzyme, particulate methane monooxygenase. In this study, the consequences of the overexpression or deletion of the Bhr gene (bhr) in Methylomicrobiam alcaliphillum 20ZR were investigated. We found that the bhrknockout (20ZRΔbhr) displays growth kinetics and methane consumption rates similar to wild type. However, the 20ZRΔbhr accumulates elevated concentrations of acetate at aerobic conditions, indicating slowed respiration. The methanotrophic strain overproducing Bhr shows increased oxygen consumption and reduced carbon-conversion efficiency, while its methane consumption rates remain unchanged. These results suggest that the methanotrophic Bhr proteins specifically contribute to oxygen-dependent respiration, while they have minimal, if any, input of oxygen for the methane oxidation machinery.

Nariya, Snehal↗
DOE OSTI2020

Synthesis and flash sintering of (Hf 1-x Zr x )B 2 solid solution powders

(Hf 1-x Zr x )B 2 solid solution powders were synthesized by two methods. First, solution-based processing of HfCl 4 , ZrCl 4 , sucrose, and H 3 BO 3 was conducted followed by heat treatment in Argon to carry out the carbothermal reduction (CTR) reaction to form the diboride powders. Alternatively, in the so-called borohydride reduction (BHR) method, HfCl 4 , ZrCl 4 and NaBH 4 were mixed in an Argon glove box followed by heat treatment in Argon at 700-1500°C. The synthesized powders were characterized by XRD, SEM, TEM, EDS, and TGA, and the influence of different parameters such as starting composition, heat treatment temperature and time on products characteristics were revealed. Both CTR and BHR solid solution powders were then consolidated within ~5 min in a homemade flash sintering (FS) setup. Here, the composition, microstructure, hardness, and thermal-oxidation properties of flash sintered ceramics were characterized, and the implication of this study and directions for future research were discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A passive model for the evolution of subgrid-scale instabilities in turbulent flow regimes

In this work, a zero-dimensional model for the development of turbulent mixing layers driven by Rayleigh–Taylor (RT) or Richtmyer–Meshkov (RM) instabilities is presented, for application in Reynolds-Averaged Navier–Stokes (RANS) simulations. The model approximates the behavior of the turbulent mixing layer at a passively advected tracer particle that sits on the material interface of interest in the simulation, and is used to provide initial conditions for the RANS model once the mixing layer width has developed to a size that may be resolved on the computational mesh. The model relies on approximations to the Besnard–Harlow–Rauenzahn (BHR) RANS closure equations and is thus consistent the turbulent development of the mixing layer, unlike previous models of this type which have focused on the early-time laminar development of RT or RM instabilities. These instabilities will often grow through an initially laminar state before developing into turbulence, and criteria for switching from an existing laminar model to the proposed turbulent model are considered. The proposed model reproduces the mixing layer development of RT and RM instabilities seen in mixed width resolving BHR RANS simulations to a reasonable degree of accuracy, and relaxes the need for the RANS simulation to resolve the mixed width as soon as the mixing layer becomes turbulent, which may be impractical in cases where the mixing layer grows through a wide range of scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗
DOE OSTI2026

Development and assessment of models for turbulent Rayleigh-Taylor mixing using the macroscopic forcing method

Reynolds-Averaged Navier Stokes (RANS) simulations are a popular method for designing ICF experiments, and accurate mixing models are crucial for these simulations to give good predictions. To this end, the present work seeks to demonstrate the Macroscopic Forcing Method (MFM) as a tool for both improving existing RANS models as well as assessing RANS model forms. First, MFM analysis from Lavacot et al. (Phys. Rev. Fluids, 2025) is used to develop the k–L–F model, an extension of the k–L model of Dimonte and Tipton (Phys. Fluids, 2006) that incorporates nonlocality through addition of a turbulent species flux transport equation. MFM is then applied to the k–L–F model along with the k–L and BHR–4 models to assess their forms and compare the model-implied eddy diffusivity moments to those measured from high-fidelity simulations. Furthermore, the analysis reveals that models incorporating nonlocality (k–L–F and BHR–4) match the high-fidelity simulation data better than purely local models (k–L), both in terms of mean fields and eddy diffusivity moments. However, all of the considered RANS models struggle to match temporal moments at high Atwood numbers, highlighting the importance of temporal nonlocality in these regimes and the need for additional improvement even among models incorporating nonlocality.

general physics↗
DOE OSTI2021

A multispecies turbulence model for the mixing and de-mixing of miscible fluids

A de-mix capable model for turbulence in compressible, variable density flows is proposed. The model is based on the Besnard-Harlow-Rauenzahn (BHR) family of models (Besnard D, Harlow F, Rauenzahn R, et al. Turbulence transport equations for variable-density turbulence and their relationship to two-field models. NM (United States): Los Alamos National Laboratory; 1992 (Technical Report LA-12303-MS), but is extended to track the evolution of the turbulent fluxes and fluctuations in the material mass fractions for each species present. The new evolution equations are introduced without requiring additional closures or new empirically tuned coefficients relative to previous BHR models, and are shown to improve the model’s ability to reproduce the behaviour of simulations containing mixing layers that are constrained by a stabilising force. As a result, the model is tested in a range of canonical flows including Rayleigh-Taylor driven, shock driven, and shear driven turbulence, and is shown to produce reasonable agreement with simulations and experiments in these scenarios.

42 ENGINEERING↗
DOE OSTI2023

Mshock HED Council Request FY24 [Slides]

ICF and HEDP experiments experience complicated loadings including successive shocks, but modeling and analytic work has mostly engaged the simpler shock or shock-reshock case. This is because the co propagating case is difficult to achieve with conventional (non-HED) drivers. Successive shocks will be challenging to model in BHR or modal model interpretations: 1) Ex: Consider the case where long wavelength modes re-invert coherently, but short wavelength nonlinear modes are spun up turbulently. 2) Necessitates more advanced diagnostics and analysis than simple “mixing layers.” High resolution spectral information (concomitant with Ω-EP campaigns) and higher-order moment analyses will be necessary to stress and validate the reduced BHR-type models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

The tracer-BHR4 model for subgrid mixing layers

A model for the passive evolution of turbulent subgrid-scale mixing layers is developed based on the model of tracer-BHR model of Braun & Gore. The new model, tracer-BHR4, mainly differs from the previous iteration in that it aims to replicate the behavior of BHR4, instead of BHR3.1. The coefficients from BHR4 are used and a new model for mixing layer growth is developed to allow the tracer model to reproduce the de-mixing behavior seen in BHR4. Additionally, an improved timestep size controller is introduced. The new model is tested on problems containing Rayleigh-Taylor or Richtmyer-Meshkov mixing layers, and is able to reproduce the behavior of well resolved BHR4 simulations with reasonable accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Additive Manufacturing Interface Profiles for Jet Studies on the Vertical Shock Tube

Understanding the mixing between two fluids of disparate densities (i.e., variable density mixing) is important for the lab's mission. Experiments are required to validate our codes, such as RAGE, FLAG, and the BHR mix model. The role initial conditions (ICs) have on instability growth and the transition to turbulence is not completely understood and comparisons with experiments are needed. The Extreme Fluids team focuses on experiments at Los Alamos that utilize gases and low shock velocities (Mach <2 in air). These experiments have the advantage of high fidelity diagnostics capable of temporally and spatially resolved velocity and density measurements. One main difficulty is setting up a specific interface profile between the two gases. While attempts have been made to do so without the aid of a solid structure (i.e., a membrane), the results are limited in the range of possible shapes and reproducibility. Membranes provide the best way to provide a deterministic initial condition. However, fragments from the membrane can interfere with the flow and instability growth, as well as the diagnostics. In this project, we focused on exploring the materials and techniques available via modern additive manufacturing (i.e., 3D printing) to test if they could address some of the shortcomings of prior membrane designs. Our goal was to find a material that would have minimal impact on our experiment while also providing us with a means to produce complex initial conditions. This report details our work to explore the range of additive manufacturing capabilities and the jet-like structures we were able to generate using three-dimensional initial conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine-Learning for Rapid Optimization of Turbulence Models

1. Background and Motivation: What is turbulence; Why is turbulence important; How turbulence is modeled; BHR model; Motivation for a ML framework 2. DNS Database from OES-C4; 3. Machine learning framework: Formulation; Verification; Case scenarios; 4. Summary and Conclusions

42 ENGINEERING↗
DOE OSTI2020

Synthesis and Flash Sintering of (Hf1-xZrx)B2 Solid Solution Fine Powders

Fine powders of (Hf1-xZrx)B2 solid solution were synthesized by two methods. In the first one, solution-based processing of HfCl4, ZrCl4, sucrose and H3BO3 was carried out followed by heat treatment (e.g., at 1500 °C for 1 h) in Argon to achieve the carbothermal reduction (CTR) reaction to form the boride solid solution. In the second one, called boron hydride reduction (BHR) method, HfCl4, ZrCl4 and NaBH4 were directly mixed in a glove box followed by heat treatment in Argon at elevated temperatures from 700 to 1500 °C. In addition, the powders synthesized via both methods were flash sintered without sample preheating in a homemade setup. The synthesized powders as well as the flash sintered bulk ceramics were characterized using different techniques including XRD, SEM, EDS, TEM, TGA-DSC, and Vickers hardness test to reveal the inter-relationships between starting materials composition, processing conditions, and the resulting materials microstructure and physical/chemical properties.

Belisario, Jose↗
DOE OSTI2022

Spring 2022 update on the status of the Local Wavenumber Model (LWN) in xRAGE

An updated implementation of the Local Wavenumber Model (LWN) is discussed, primarily differing from recent versions by placing a greater focus on capturing a wide variety of turbulent flows including compressible flows. New models are introduced for spectral backscatter effects, the effect of bulk compression on the spectra, and incorporating the multispecies variables tracked in BHR4. Methods for reducing the compurational expense of tracking spectra for turbulent quantites are also investigated. Like recent versions of BHR, we test LWN in a number of canonical flows using a single set of coefficients, but additional coefficient tuning is likely to be required to improve the agreement with some of these flows.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗