Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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 145 records · Page 8

Massive Galaxies Impede Massive Outflows

A set of 66 3D hydrodynamical simulations explores how galactic stellar mass affects three-phase, starburst-driven outflows. Simulated velocities are compared to two basic analytic models: with and without a gravitational potential. For stellar mass <1010 M ⊙, simulated velocities match those of both analytical models and are unaffected by the potential; above they reduce significantly as expected from the analytic model with gravity. Gravity also affects total outflow mass and each of the three phases differently. Outflow masses in the hot, warm, and cold phases each scale with stellar mass as $\mathrm{log}{M}_{* }\,=$ −0.25, −0.97, and −1.70, respectively. Thus, the commonly used Chevalier & Clegg analytic model should be modified to include gravity when applied to higher-mass galaxies. In particular, using M82 as the canonical galaxy to interpret hydrodynamical simulations of starburst-driven outflows from higher-mass galaxies will underestimate the retarding effect of gravity. Using the analytic model of Johnson & Axford with realistic thermalization efficiency and mass loading, I find that only galaxy masses lesssim10(exp 11.5) M(⊙) can outflow.

Ryan Tanner↗

Hierarchical Coupling of Molecular Dynamics and Micromechanics to Predict the Elastic Properties of Three-Phase and Four-Phase Silicon Carbide Composites

The results obtained from previously conducted molecular dynamics analysis of silicon carbide (-SiC (6H, 4H, & 2H-SiC), -SiC (3C SiC)), silicon and boron nitride, were utilized as inputs in the MAC/GMC micromechanics software to model and evaluate the elastic properties of three-phase SiC/BN/SiC and four-phase SiC/BN/Si/SiC composites. This method of analysis eliminates the need for back-calculation of the apparent properties of the base constituents from the measured ceramic matrix composites properties. The multiscale models are validated against the available data in literature.

Aluko, Olanrewaju↗

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An international benchmark for wind plant wakes from the American WAKE ExperimeNt (AWAKEN)

This article introduces the first benchmark study within the International Energy Agency Wind Task 57 framework, focusing on wind plant wakes. Leveraging data from the American WAKE ExperimeNt (AWAKEN), the benchmark aims to assess the accuracy of simulation tools in modeling wind plant wakes and their impact on the downstream flow under diverse inflow conditions. The AWAKEN field campaign, conducted in Oklahoma from 2022 to 2024, provides unprecedented observations of wind plant-atmosphere interactions, thus offering a large dataset to validate numerical models of different complexity. The benchmark will include three phases—code calibration, blind comparison, and iteration—allowing participants to refine their numerical models based on the feedback from the benchmark team. This article describes the benchmark case study selected from observations providing details on atmospheric conditions, wake evidence, and wind turbine operation. The benchmark’s structure and timeline, along with the expected publication of results, are discussed as well. This collaborative effort aims to enhance the accuracy of wind plant wake simulations, thus contributing to the improvement of wind energy production estimates.

17 WIND ENERGY↗

Importance of Site Diversity and Connectivity in Electrochemical CO Reduction on Cu

Electrochemical CO 2 reduction on Cu is a promising approach to produce value-added chemicals using renewable feedstocks, yet various Cu preparations have led to differences in activity and selectivity toward single and multicarbon products. Here, we find, surprisingly, that the effective catalytic activity toward ethylene improves when there is a larger fraction of less active sites acting as reservoirs of *CO on the surface of Cu nanoparticle electrocatalysts. In an adaptation of chemical transient kinetics to electrocatalysis, we measure the dynamic response of a gas diffusion electrode (GDE) cell when the feed gas is abruptly switched between Ar (inert) and CO. When switching from Ar to CO, CO reduction (COR) begins promptly, but when switching from CO to Ar, COR can be maintained for several seconds (delay time) despite the absence of the CO reactant in the gas phase. A three-site microkinetic model captures the observed dynamic behavior and shows that Cu catalysts exhibiting delay times have a less active *CO reservoir that exhibits fast diffusion to active sites. The observed delay times and the estimated *CO reservoir sizes are affected by catalyst preparation, applied potential, and microenvironment (electrolyte cation identity, electrolyte pH, and CO partial pressure). Notably, we estimate that the *CO reservoir surface coverage can be as high as 88 ± 7% on oxide-derived Cu (OD-Cu) at high overpotentials (–1.52 V vs SHE) and this increases in reservoir coverage coincide with increased turnover frequencies to ethylene. We also estimate that *CO can travel substantial distances (up to 10s of nm) prior to desorption or reaction. It appears that active C–C coupling sites by themselves do not control selectivity to C 2+ products in electrochemical COR; the supply of CO to those sites is also a crucial factor. More generally, the overall activity of Cu electrocatalysts cannot be approximated from linear combinations of individual site activities. Future designs must consider the diversity of the catalyst network and account for intersite transportation pathways.

30 DIRECT ENERGY CONVERSION↗

Metallothionein expression in Eastern mosquitofish ( Gambusia holbrooki ) exposed to 137 Cs and lead: implications for using metallothionein as a metal exposure biomarker

Given the range of threats facing aquatic ecosystems, it is critical we have monitoring tools to quickly and efficiently identify stressors within ecosystems. Metallothionein’s (MTs) are a group of proteins widely used as biomarkers of metal exposure; however, their suitability has been questioned due to their influenceable nature and ability to be induced by various stressors. Through our study, we sought to (1) evaluate induction of MT in eastern mosquitofish (Gambusia holbrooki) exposed to a radionuclide ( 137 Cs) and metal (Pb), (2) assess the influence of in situ 137 Cs pre-exposure on Pb-induced MT expression, and (3) quantify changes in MT concentrations during Pb depuration. We conducted our experiment in three phases using mosquitofish as a model species: Phase I—pre-exposure to 137 Cs in the field, Phase II—exposure to Pb in the lab immediately following 137 Cs exposure, and Phase III—depuration in the laboratory. We quantified MT concentrations over time and assessed differences among treatments. We did not find in situ exposure of mosquitofish to 137 Cs to induce MT, nor did pre-exposure to 137 Cs appear to influence induction of MT following Pb exposure. Metallothionein was rapidly induced in mosquitofish following Pb exposure; however, MT induction was also observed in control fish. Additionally, there was no apparent relationship between MT and Pb body burdens. During depuration, we found MT concentrations to be highly variable among fish demonstrating no clear response. Furthermore, our study demonstrates the high variability exhibited by MT, the need for further validation of MT as a metal exposure biomarker, and the value in additional research on these proteins under varying exposure scenarios, particularly environmentally realistic concentrations.

54 ENVIRONMENTAL SCIENCES↗

Comparison of Spacecraft Contamination Models with Well-Defined Flight Experiment

The report presents analyzed surface areas on particular experiment trays from the Long Duration Exposure Facility (LDEF) for silicone-based molecular contamination. The trays for examination were part of the Ultra-Heavy Cosmic Ray Experiment (UHCRE). These particular trays were chosen because each tray was identical to the others in construction, and the materials on each tray were well known, documented, and characterized. In particular, a known specific source of silicone contamination was present on each tray. Only the exposure conditions varied from tray to tray. The results of post-flight analyses of surfaces of three trays were compared with the predictions of the three different spacecraft molecular contamination models. Phase one tasks included: 1) documenting the detailed geometry of the hardware; 2) determining essential properties of the anodized aluminum, Velcro(Tm), silverized Teflon(Tm), silicone gaskets, and DC6-1104(Tm) silicone adhesive materials used to make the trays, tray covers, and thermal control blankets; 3) selecting and removing areas from each tray; and 4) beginning surface analysis of the selected tray walls. Phase two tasks included: 1) completion of surface analysis measurements of the selected tray surface, 2) obtaining auger depth profiles at selected locations, and 3) running versions of the ISEM, MOFLUX, and PLIMP (Plume Impingement) contamination prediction models and making comparisons with experimental results.

Pippin, G. H.↗

Evolution of Dendritic Extended 3D Patterns during Directional Solidification: Microgravity Experiments in DECLIC-DSI onboard ISS and Phase-field Modeling

To characterize the dynamical formation of three-dimensional (3D) arrays of dendrites under diffusive growth conditions, in situ monitoring of a series of experiments on a transparent succinonitrile – 0.46 wt% camphor model alloy was carried out under low gravity in the DECLIC Directional Solidification Insert onboard the International Space Station. These experiments offer continuous interface observation and enable the construction of space-time evolution maps of dendrite location and primary spacing during directional solidification. We present the evolution of dendritic extended 3D patterns for a large range of parameters displaying different levels of side branching in both experiments and phase-field simulations. The combined and separate effects of macroscopic interface curvature and crystal mis-orientation are evidenced and detailed. In addition, the aligned dendritic array caused by the invasion of dendrites formed at the edge of the crucible border is investigated.

Kaihua Ji↗

Three-dimensional aspects of radiative transfer in remote sensing of precipitation: Application to the 1986 COHMEX storm

Several multifrequency techniques for passive microwave estimation of precipitation based on the absorption and scattering properties of hydrometers have been proposed in the literature. In the present study, plane-parallel limitations are overcome by using a model based on the discrete-ordinates method to solve the radiative transfer equation in three-dimensional rectangular domains. This effectively accounts for the complexity and variety of radiation problems encountered in the atmosphere. This investigation presents result for plane-parallel and three-dimensional radiative transfer for a precipitating system, discusses differences between these results, and suggests possible explanations for these differences. Microphysical properties were obtained from the Colorado State University Regional Atmospehric Modeling System and represent a hailstorm observed during the 1986 Cooperative Huntsville Meteorological Experiment. These properties are used as input to a three-dimensional radiative transfer model in order to simulate satellite observation of the storm. The model output consists of upwelling brightness temperatures at several of the frequencies on the Special Sensor Microwave/Imager. The radiative transfer model accounts for scattering and emission of atmospheric gases and hydrometers in liquid and ice phases. Brightness temperatures obtained from the three-dimensional model of this investigation indicate that horizontal inhomogeneities give rise to brightness temperature fields that can be quite different from fields obtained using plane-parallel radiative transfer theory. These differences are examined for various resolutions of the satellite sensor field of view. In adddition, the issue of boundary conditions for three-dimensional atmospheric radiative transfer is addressed.

Haferman, J. L.↗

One Earth Energy Static and Dynamic Reservoir Modeling

This report presents the static and dynamic reservoir modeling conducted for the CarbonSAFE Phase III Illinois Storage Corridor project to assess the feasibility of commercial-scale CO 2 storage in the Mt. Simon Sandstone at the One Earth Energy (OEE) site in McLean County, Illinois. Three-dimensional geocellular models of the Mt. Simon storage complex were developed in Petrel ® by integrating petrophysical log data, core analyses, and seismic surveys from the OEE #1 stratigraphic test well and two nearby wells, with multiple model versions created as new data became available. Dynamic reservoir simulations, performed using Landmark's Nexus software, progressed through three phases (preliminary, sensitivity, and UIC Class VI permit studies) evaluating injection scenarios across varying rates, well orientations, permeability models, and multi-well configurations. Results demonstrate that commercial-scale storage is feasible: three injection wells spaced approximately one mile apart can store a total of 90 million tonnes of CO 2 over 20 years, producing a combined plume with an equivalent radius of 3.2 miles and a maximum pressure-front-defined Area of Review of 178 mi 2 at the end of injection that diminishes to 34 mi 2 after 50 years of post-injection monitoring. Sensitivity analyses indicate that a 20% change in porosity or permeability yields approximately a 7% change in AoR radius, and that perforating the high-permeability arkosic zone minimizes the pressure front compared to injection in the upper Mt. Simon Sandstone.

09 BIOMASS FUELS↗

Large-Signal Stability of Phase-Balanced Equilibria in Single-Phase Grid-Forming Inverter Systems

This article explores the setup where large numbers of single-phase grid-forming inverters with droop control across distribution networks self-organize into a stable and balanced system with 120° phase offsets across aggregates in the absence of balanced three-phase generating resources or external communication. A suite of circuit- and system-theoretic notions are leveraged to derive a dynamical model for phase-angle differences across aggregates of inverters connected in the three phases. Focusing on this model, large-signal stability is established and the region of attraction of the phase-balanced equilibria is determined with the aid of a Lyapunov function. Experimental validation for a bench-top prototype network is included to support the analytical developments. Altogether, the effort supports the vision of facilitating balanced operation of distribution networks with grid-forming inverters during service disruptions at the bulk transmission network.

Droop control↗

Lattice model theory of the equation of state covering the gas, liquid, and solid phases

The three stable states of matter and the corresponding phase transitions were obtained with a single model. Patterned after Lennard-Jones and Devonshires's theory, a simple cubic lattice model containing two fcc sublattices (alpha and beta) is adopted. The interatomic potential is taken to be the Lennard-Jones (6-12) potential. Employing the cluster variation method, the Weiss and the pair approximations on the lattice gas failed to give the correct phase diagrams. Hybrid approximations were devised to describe the lattice term in the free energy. A lattice vibration term corresponding to a free volume correction is included semi-phenomenologically. The combinations of the lattice part and the free volume part yield the three states and the proper phase diagrams. To determine the coexistence regions, the equalities of the pressure and Gibbs free energy per molecule of the coexisting phases were utilized. The ordered branch of the free energy gives rise to the solid phase while the disordered branch yields the gas and liquid phases. It is observed that the triple point and the critical point quantities, the phase diagrams and the coexistence regions plotted are in good agreement with the experimental values and graphs for argon.

Bonavito, N. L.↗

Modeling framework for evaluating the impacts of hydrodynamic pressure on hydrologic exchange fluxes and residence time for a large-scale river section over a long-term period

Quantifying hydrologic exchange fluxes (HEF) at the river and subsurface interface and their residence times (RT) in subsurface are important for managing the water quality and ecosystem health in dynamic river corridor systems. In this study, a modeling framework is developed for coupling the three-dimensional (3D) multi-phase surface, subsurface flow transport, and numerical tracer model for RT in a large-scale river section over a long period. The framework is utilized to evaluate the impacts of hydrodynamic pressure on HEFs and RT for a 30 km section of the Columbia River in Washington State over a three-year period. Based on comparisons between model simulations with and without considering hydrodynamic pressure, we found that hydrodynamic pressure increases the net HEFs by 7% with river gaining water from the subsurface domain, and leads to slight reduction of RT.

54 ENVIRONMENTAL SCIENCES↗

Flame spread across liquid pools

For flame spread over liquid fuel pools, the existing literature suggests three gravitational influences: (1) liquid phase buoyant convection, delaying ignition and assisting flame spread; (2) hydrostatic pressure variation, due to variation in the liquid pool height caused by thermocapillary-induced convection; and (3) gas-phase buoyant convection in the opposite direction to the liquid phase motion. No current model accounts for all three influences. In fact, prior to this work, there was no ability to determine whether ignition delay times and flame spread rates would be greater or lesser in low gravity. Flame spread over liquid fuel pools is most commonly characterized by the relationship of the initial pool temperature to the fuel's idealized flash point temperature, with four or five separate characteristic regimes having been identified. In the uniform spread regime, control has been attributed to: (1) gas-phase conduction and radiation; (2) gas-phase conduction only; (3) gas-phase convection and liquid conduction, and most recently (4) liquid convection ahead of the flame. Suggestions were made that the liquid convection was owed to both vuoyancy and thermocapillarity. Of special interest to this work is the determination of whether, and under what conditions, pulsating spread can and will occur in microgravity in the absence of buoyant flows in both phases. The approach we have taken to resolving the importance of buoyancy for these flames is: (1) normal gravity experiments and advanced diagnostics; (2) microgravity experiments; and (3) numerical modelling at arbitrary gravitational level.

Ross, Howard↗

VRA Modeling, phase 1

The destruction of organic contaminants in waste water for closed systems, such as that of Space Station, is crucial due to the need for recycling the waste water. A co-current upflow bubble column using oxygen as the gas phase oxidant and packed with catalyst particles consisting of a noble metal on an alumina substrate is being developed for this process. The objective of this study is to develop a plug-flow model that will predict the performance of this three phase reactor system in destroying a multicomponent mixture of organic contaminants in water. Mass balances on a series of contaminants and oxygen in both the liquid and gas phases are used to develop this model. These mass balances incorporate the gas-to-liquid and liquid-to-particle mass transfer coefficients, the catalyst effectiveness factor, and intrinsic reaction rate. To validate this model, a bench scale reactor has been tested at Michigan Technological University at elevated pressures (50-83 psig,) and a temperature range of 200 to 290 F. Feeds consisting of five dilute solutions of ethanol (approx. 10 ppm), chlorobenzene (approx. 20 ppb), formaldehyde (approx. 100 ppb), dimethyl sulfoxide (DMSO approx. 300 ppb), and urea (approx. 20 ppm) in water were tested individually with an oxygen mass flow rate of 0.009 lb/h. The results from these individual tests were used to develop the kinetic parameter inputs necessary for the computer model. The computer simulated results are compared to the experimental data obtained for all 5 components run in a mixture on the differential test column for a range of reactor contact times.

Kindt, Louis M.↗

Illinois State Geological Survey (ISGS), Illinois Basin - Decatur Project (IBDP) Geological Models, July 7, 2021. Midwest Geological Sequestration Consortium (MGSC) Phase III Data Sets. DOE Cooperative Agreement No. DE-FC26-05NT42588.

Three geological models in Petrel (Mark of Schlumberger) and data output from shallow groundwater modeling using TOUGH Codes from Lawrence Berkeley National Laboratory, included under folders: Static_Geologic_Model, Dynamic_Reservoir_Model, Geomechanical_Model, and Groundwater_Model.

3D Geological Model↗

Predictions of the structure of turbulent, particle-laden, round jets

Models of gas and particle motion in turbulent, particle-laden, round jets were evaluated using existing measurements of flow structure. Three models were considered: (1) a locally homogeneous flow model, where velocities and turbulent mixing properties of both phases were assumed to be equal; (2) a deterministic separated flow model, where interphase slip was considered but effects of turbulent dispersion were ignored; and (3) a stochastic separated flow model where effects of interphase slip and turbulent dispersion were considered using random sampling techniques. In all three cases, mean and turbulent properties of the continuous phase were found with a well-calibrated k-epsilon model. The locally homogeneous flow and deterministic separated flow models over- and underestimated particle spread and flow development rates, respectively. The stochastic separated flow model, however, yielded satisfactory predictions of flow structure - except at high particle loadings. Uncertainties in initial conditions for the measurements and possible effects of turbulence modulation by the particles are proposed as the reason for these errors.

Shuen, J.-S.↗