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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

A Priori Subgrid Scale Modeling for a Droplet Laden Temporal Mixing Layer

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using a direct numerical simulation (DNS) database. The DNS is for a Reynolds number (based on initial vorticity thickness) of 600, with droplet mass loading of 0.2. The gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. Since Large Eddy Simulation (LES) of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be given by the filtered variables plus a correction based on the filtered standard deviation, which can be computed from the sub-grid scale (SGS) standard deviation. This model predicts unfiltered variables at droplet locations better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: Smagorinsky, gradient and scale-similarity. When properly calibrated, the gradient and scale-similarity methods give results in excellent agreement with the DNS.

Okongo, Nora↗
NASA NTRS2019

Pilot Injector Redesign to Reduce N+3 Cycle Emissions for a Gas-Turbine Combustor

An overview is given of an effort for the use of CFD analysis to complement design and configuration definition of third generation Lean-Direct Injection combustion concepts (LDI-3) for NASA’s N+3 program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for a three-cup, nineteen-element flame tube array with redesigned pilot injectors to improve spray and emissions characteristics when compared to a previous LDI-3 design. All computations were performed with a consistent approach to mesh-generation, spray modeling, ignition and kinetics modeling for a ‘medium-power’ cycle condition. Computational predictions of the aerodynamics of a new pre-filming pilot injector were used to arrive at an optimized aerothermal design that meets effective area and fuel-air mixing criteria. The newly designed pilot injectors were shown to provide considerable improvements in aerodynamic stability, flame-tube pattern factor and NOx emissions, when compared to the original design.

nitrogen oxides injectors↗
NASA NTRSDate not supplied

Modeling of Particle Transport and Impact in Transient, Multiphase Oxidizer-Rich Environments Using Computational Fluid Dynamics and Lagrangian Particle Tracking

Particle impacts from foreign object debris (FOD) in an oxidizer-rich environment have been shown to be a credible and significant risk for ignition and failure of components and vehicles. The risk of particle ignition is dependent on many factors, including particle impact velocity/energy and fluid composition (i.e. fluid phase, species presence, and state). Prediction of particle transport and particle impacts in relevant environments, including transient, multiphase flows, is important for assessing the risk and avoiding repeated expensive system level tests. The capability to model both the fluid dynamics and particle transport in complex three-dimensional geometries and multiphase environments has been developed at Marshall Space Flight Center (MSFC) in the propulsion fluid dynamic branch (ER42) by utilizing computational fluid dynamics (CFD) simulations and Lagrangian particle tracking. Using CFD tools, including density-based Loci/CHEM with a cavitation model and pressure-based Loci/STREAM with a volume of fluids model, MSFC/ER42 has provided particle impact predictions to multiple projects for oxidizer compatibility assessments. Simulations of transient, multiphase environments such as liquid slugs flowing through tubing and cavitating flow in valves have been completed alongside one-way coupled Lagrangian particle tracking to predict particle impacts.

M.V. Fischels↗
NASA NTRS1995

The anomalous nuclear component in the three-dimensional heliosphere

The anomalous nuclear component is neither of solar nor galactic cosmic ray origin. Its presence in the heliosphere is an independent probe for both interplanetary electrodynamical investigations--especially solar modulation--and probably the most direct means for determining the elemental and isotopic composition of those neutral atoms in the local interstellar medium that have high first ionization potentials (e.g., He, N, O, Ne, Ar, etc.). This report is a brief account of the evolution of our knowledge of this component for readers not specializing in the subject. Included are the initial discoveries of the component, its composition, spectra, heliospheric radial and latitudinal intensity gradients, modulation over the approximately 22 year solar magnetic cycle, trapping in the magnetosphere and its use for estimating the location of a heliospheric termination shock. Recent measurements from the ULYSSES spacecraft have provided conclusive evidence that incoming neutral atoms, after photoionization by solar uv, are picked up by the solar wind, thus lending support for the concept that after their acceleration--probably at a termination shock--they return to the inner heliosphere as pseudo-cosmic rays. ULYSSES spacecraft investigations extending to approximately 56 deg south latitude reveal, for both the anomalous nuclear component and the galactic cosmic rays, that there is a surprisingly small latitudinal intensity gradient. Thus, for the current phase of the solar cycle, modulation is much more spherically symmetric in the inner solar system than had generally been believed. A further surprise is the continual presence of approximately 26 day recurrent modulation at hight latitudes, without corresponding magnetic field compressions. These results are changing our ideas and requiring modification of our models for solar modulation in three dimensions.

Simpson, J. A.↗
DOE OSTI2021

Multitude of topological phase transitions in bipartite dice and Lieb lattices with interacting electrons and Rashba coupling

We report the results of a Hartree-Fock study applied to interacting electrons moving in two different bipartite lattices: the dice and the Lieb lattices, at half-filling. Both lattices develop ferrimagnetic order in the phase diagram U-λ, where U is the Hubbard on-site repulsion and λ is the Rashba spin-orbit coupling strength. Our main result is the observation of an unexpected multitude of topological phases for both lattices. All these phases are ferrimagnetic, but they differ among themselves in their set of six Chern numbers (six numbers because the unit cells have three atoms). Here, the Chern numbers |C| observed in our study range from 0 to 3, showing that large Chern numbers can be obtained by the effect of electronic correlations, adding to the recently discussed methodologies to increase |C| based on extending the hopping range in tight-binding models, using sudden quenches, or photonic crystals, all without including electronic interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Big Data Synchrophasor Monitoring and Analytics for Resiliency Tracking (BDSMART)

This report contains key findings from a project titled Big Data Synchrophasor Monitoring and Analytics for Resiliency Tracking (BDSMART), which was carried out through a collaborative effort of a team of researchers from Texas A&M Engineering Experiment Station, Temple University, and Quanta Technology, LLC. The in-kind support came from OSIsoft (acquired by AVEVA), which provided their PI Historian software to demonstrate the use case of streaming PMU data. The first section of the report describes the project goals and objectives related to the development of Machine Learning (ML) models capable of detecting and classifying events by processing phasor measurements captured in the field by Phasor Measurement Units (PMUs). The data for this study was contributed by the utilities/ISOs from the Western and Eastern interconnects and ERCOT, further referred to as Interconnect B (IC B), Interconnect A (IC A), and Interconnect C (IC C), respectively. The approach that the BDSMART Research Team proposed and the key research tasks defined by the team are outlined in this section. The next section describes the technical approach. We first discuss the data constraints related to the PMU measurements and data interpretation constraints imposed by the data contributors. They provided neither the topological information of the grid nor PMU placement locations and captured recorded data at very few locations in the system with the reporting rate of either 30 or 60 fps. The recordings are mostly positive sequence voltage, frequency, and ROCOF, and in some limited cases, three-phase voltages and currents. We then reflect on the bad data issues that stem from poor recording practices and vague definitions of the PMU status bits to supposedly be used for bad data identification. Finally, the data discovery points to imprecise time stamps with incomplete event start/end time, as well as inconsistent and incomplete event labeling, which combined make the implementation of the data models using supervising learning quite challenging. Following the data discovery study, we hypothesize that because the IC B data has the most complete labels, we should focus our model development on that data and then test it on data from other interconnects. We also define the common metrics used to evaluate the results from the ML algorithm tests. We concluded this section by summarizing the common ML models we used and explaining how we implemented and tested them. The issues from this section are expanded in the Training Dataset Report from this project. The final section of this report deals with the accomplishments and conclusions. As the accomplishments, we formulate the problem we are solving and what is achieved by solving the problem. We then reflect on each of the analytics tools we developed and point out the performance of each tool when applied to solving the mentioned problems. We reference this work for further details to the papers we published on each tool. In the conclusions, we give recommendations on how to improve future PMU recording practices to facilitate the ML algorithm implementation and guidance for the future standardization work aimed at clarifying the ambiguities associated with the PMU status bits. We finally list future tasks that can bring about further improvements in the proposed algorithms. The issues from this section are expanded in the Training, and Test Dataset Report filed at the project completion date.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Active learning emulators for nuclear two-body scattering in momentum space

In this work we extend the active learning emulators for two-body scattering in coordinate space with error estimation, recently developed by Maldonado et al. [Phys. Rev. C 112, 024002], to coupled-channel scattering in momentum space. Our full-order model (FOM) solver is based on the Lippmann-Schwinger integral equation for the scattering t-matrix as opposed to the radial Schrödinger equation. We use (Petrov-)Galerkin projections and high-fidelity calculations at a few snapshots across the parameter space of the interaction to construct efficient reduced-order models (ROMs), trained by a greedy algorithm for locally optimal snapshot selection. Both the FOM solver and the corresponding ROMs are implemented efficiently in Python using Google's JAX library. We present results for emulating scattering phase shifts in coupled and uncoupled channels and cross sections, and assess the accuracy of the developed ROMs and their computational speedup factors. We also develop emulator error estimation for both the t-matrix and the total cross section. The software framework for reproducing and extending our results is publicly available. Together with our recent advances in developing active-learning emulators for three-body scattering, these emulator frameworks set the stage for full Bayesian calibrations of chiral nuclear interactions and optical models against scattering data with quantified emulator errors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

VWC-BERT: Scaling Vulnerability–Weakness–Exploit Mapping on Modern AI Accelerators

Defending cybersystems needs accurate mapping of software and hardware vulnerabilities to generalized descriptions of weaknesses, and weaknesses to exploits. These mappings enable cyber defenders to build plans for effective defense and assessment of potential risks to a cybersystem. With close to 170k vulnerabilities, manual mapping is not a feasible option. However, automated mapping is challenging due to limited training data, computational intractability, and limitations in computational natural language processing. Tools based on breakthroughs in Transformer-based language models have been demonstrated to classify vulnerabilities with high accuracy. We make three key contributions in this paper: (1) We present a new framework, \VWCBERT, that augments the Transformer-based hierarchical multi-class classification framework of Das et al. (\textsc{V2W-BERT}) with the ability to map weaknesses to exploits. (2) We implement \VWCBERT~ on modern AI accelerator platforms using two data parallel techniques for the pre-training phase and demonstrate nearly linear speedups across NVIDIA and Graphcore accelerator platforms. We observe nearly linear speedups for up to 16 V100 and 8 A100 GPUs, and about 3.4$\times$ speedup for A100 relative to V100 GPUs. We also observe excellent speedups on Graphcore, with $5.7\times$ speedup on 128 IPUs relative to 16 IPUs. Enabled by scaling, we also demonstrate higher accuracy using a larger language model, RoBERTa-Large. We show up to 87\% accuracy for strict and up to 98\% accuracy for relaxed classification. (3) We develop a novel parallel link manager for the link prediction phase and demonstrate up to 21$\times$ speedup with 16 V100 GPUs relative to one V100 GPU, and thus reducing the runtime from 2.5 hours to 10 minutes. We believe that generalizability and scalability of \VWCBERT~ will benefit both the theoretical development and practical deployment of novel cyberdefense solutions and vulnerability classification.

Das, Siddhartha Shankar↗
NASA NTRS2011

Multiangle Implementation of Atmospheric Correction (MAIAC):: Radiative Transfer Basis and Look-up Tables - 1

This paper describes a radiative transfer basis of the algorithm MAIAC which performs simultaneous retrievals of atmospheric aerosol and bidirectional surface reflectance from the Moderate Resolution Imaging Spectroradiometer (MODIS). The retrievals are based on an accurate semianalytical solution for the top-of-atmosphere reflectance expressed as an explicit function of three parameters of the Ross-Thick Li-Sparse model of surface bidirectional reflectance. This solution depends on certain functions of atmospheric properties and geometry which are precomputed in the look-up table (LUT). This paper further considers correction of the LUT functions for variations of surface pressure/height and of atmospheric water vapor, which is a common task in the operational remote sensing. It introduces a new analytical method for the water vapor correction of the multiple ]scattering path radiance. It also summarizes the few basic principles that provide a high efficiency and accuracy of the LUT ]based radiative transfer for the aerosol/surface retrievals and optimize the size of LUT. For example, the single-scattering path radiance is calculated analytically for a given surface pressure and atmospheric water vapor. The same is true for the direct surface-reflected radiance, which along with the single-scattering path radiance largely defines the angular dependence of measurements. For these calculations, the aerosol phase functions and kernels of the surface bidirectional reflectance model are precalculated at a high angular resolution. The other radiative transfer functions depend rather smoothly on angles because of multiple scattering and can be calculated at coarser angular resolution to reduce the LUT size. At the same time, this resolution should be high enough to use the nearest neighbor geometry angles to avoid costly three ]dimensional interpolation. The pressure correction is implemented via linear interpolation between two LUTs computed for the standard and reduced pressure levels. A linear mixture and a modified linear mixture methods are used to represent different aerosol types in the aerosol/surface retrievals from several base models of the fine and coarse aerosol fractions. In summary, the developed LUT algorithm allows fast high-accuracy simulations of the outgoing radiance with full variability of the atmospheric and surface bidirectional reflectance properties for the aerosol/surface remote sensing.

Lyapustin, Alexei↗
NASA NTRS1998

Construction of a Direct Water-Injected Two-Stroke Engine for Phased Direct Fuel Injection-High Pressure Charging Investigations

The development of a water injected Orbital Combustion Process (OCP) engine was conducted to assess the viability of using the powerplant for high altitude NASA aircraft and General Aviation (GA) applications. An OCP direct fuel injected, 1.2 liter, three cylinder, two-stroke engine has been enhanced to independently inject water directly into the combustion chamber. The engine currently demonstrates low brake specific fuel consumption capability and an excellent power to weight ratio. With direct water injection, significant improvements can be made to engine power, to knock limits/ignition advance timing, and to engine NO(x) emissions. The principal aim of the testing was to validate a cyclic model developed by the Systems Analysis Branch at NASA Ames Research Center. The work is a continuation of Ames' investigations into a Phased Direct Fuel Injection Engine with High Pressure Charging (PDFI-ITPC).

Somsel, James P.↗
NASA NTRS1979

Two-phase choked flow of cryogenic fluids in converging-diverging nozzles

Data are presented for the two phase choked flow of three cryogenic fluids - nitrogen, methane, and hydrogen - in four converging-diverging nozzles. The data cover a range of inlet stagnation conditions, all single phase, from well below to well above the thermodynamic critical conditions. In almost all cases the nozzle throat conditions were two phase. The results indicate that the choked flow rates were not very sensitive to nozzle geometry. However, the axial pressure profiles, especially the throat pressure and the point of vaporization, were very sensitive to both nozzle geometry and operating conditions. A modified Henry-Fauske model correlated all the choked flow rate data to within + or - 10 percent. Neither the equilibrium model nor the Henry-Fauske model predicted throat pressures well over the whole range of data. Above the thermodynamic critical temperature the homogeneous equilibrium model was preferred for both flow rate and pressure ratio. The data of the three fluids could be normalized by the principle of corresponding states.

Simoneau, R. J.↗
DOE OSTI2024

Iridium Surface Oxide Affects the Nafion Interface in Proton-Exchange-Membrane Water Electrolysis

Proton-exchange-membrane water electrolyzer (PEMWE) catalyst layers consist of aggregates of catalyst particles (typically iridium) and ionomer (typically Nafion). Prior work suggests that the oxide form of Ir affects the kinetics of the oxygen-evolution reaction. However, because most catalyst-benchmarking studies are conducted ex situ in liquid electrolytes, it remains unclear how the ionomer is influenced by the catalyst oxide and affects overall cell performance. Using a suite of experimental techniques, we conduct fundamental investigations into model ink (catalyst and ionomer dispersed in solution) and thin-film systems to inform cell-level overpotential analysis as a function of three forms of Ir (metallic Ir m , oxyhydroxide IrOOH, and oxide IrO 2 ). Furthermore, nafion on Ir m has a high degree of phase separation and higher swelling, likely improving the ionic conductivity. Additionally, Nafion binds most strongly to IrOOH, likely yielding reduced kinetic overpotentials. These findings highlight the intricacies of the ionomer/Ir interface and provide insight into all catalyst-layer systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Multi-decadal variation of ENSO forecast skill since the late 1800s

Diagnosing El Niño-Southern Oscillation (ENSO) predictability within operational forecast models is hindered by computational expense and the need for initialization with three-dimensional fields generated by global data assimilation. We instead examine multi-year ENSO predictability since the late 1800s using the model-analog technique, which has neither limitation. We first draw global coupled model states from pre-industrial control simulations, from the Coupled Model Intercomparison Project Phase 6, that are chosen to initially match observed monthly sea surface temperature and height anomalies in the Tropics. Their subsequent 36-month model evolution are the hindcasts, whose 20th century ENSO skill is comparable to twice-yearly hindcasts generated by a state-of-the-art European operational forecasting system. Despite the so-called spring predictability barrier, present throughout the record, there is substantial second-year ENSO skill, especially after 1960. Overall, ENSO exhibited notably high values of both amplitude and skill towards the end of the 19th century, and again in recent decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Structure of a seeded palladium nanoparticle and its dynamics during the hydride phase transformation

Palladium absorbs large volumetric quantities of hydrogen at room temperature and ambient pressure, making the palladium hydride system a promising candidate for hydrogen storage. Here, we use Bragg coherent diffraction imaging to map the strain associated with defects in three dimensions before and during the hydride phase transformation of an individual octahedral palladium nanoparticle, synthesized using a seed-mediated approach. The displacement distribution imaging unveils the location of the seed nanoparticle in the final nanocrystal. By comparing our experimental results with a finite-element model, we verify that the seed nanoparticle causes a characteristic displacement distribution of the larger nanocrystal. During the hydrogen exposure, the hydride phase is predominantly formed on one tip of the octahedra, where there is a high number of lower coordinated Pd atoms. Our experimental and theoretical results provide an unambiguous method for future structure optimization of seed-mediated nanoparticle growth and in the design of palladium-based hydrogen storage systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Machine Learning in Environmental Chemistry: Application to Surface Complexation Modeling

Environmental chemistry – or biogeochemistry – is the scientific discipline typically invoked when examining and quantifying groundwater or surface water contamination, and nutrient cycling in the environment. Over the last three decades, there have been significant advances in mechanistic model development to describe and predict these complex biogeochemical processes. In particular, surface complexation models (SCMs) have been developed to describe the rock/soil surface reactions of metals and radionuclides, and their partitioning between various mobile species in the aqueous phase or immobile species sorbed on solid surfaces. Often represented by a simplified linear isotherm constant – Kd – in reactive transport models, these reactions play a critical role in many environmental science applications; particularly in contamination risk assessments and nuclear waste disposal performance assessments. In the past several decades, efforts by various institutions across the world have focused on developing SCMs based on datasets from laboratory measurements, including the identification of key parameters such as equilibrium constants.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

Global sensitivity to velocity errors at the libration points

The global effects on non-linear stability are investigated at and around the triangular libration points. The model of the circular restricted problem of three bodies is used, considering the earth and the moon as the primaries. Areas of the initial conditions in the configuration space leading to stability with zero initial velocity around the equilibrium points are compared to stable areas in the phase space with variable initial velocity and zero initial deviations at the equilibrium points. The behavior of the system concerning errors due to initial conditions in the configuration space is studied for 6 years and the effects of initial velocity errors are considered for 7.73 years.

Szebehely, V.↗
NASA NTRS1991

Modeling the growth phase of a substorm using the Tsyganenko model and multi-spacecraft observations - CDAW-9

The CDAW-9 Event C focused upon the early part of 3 May 1986 when a large substorm onset occurred at 0111 UT. By modifying the Tsyganenko 1989 magnetic field model, a model is constructed in which the near-earth current systems are enhanced with time to describe the observed development of the tail magnetic field during the growth phase. The cross-tail current intensity and the thickness of the current sheet are determined by comparison with three spacecraft in the near-earth tail. The location of the auroral bulge as recorded by the Viking imager is mapped to the equatorial current sheet. The degree of chaotization of the thermal electrons is estimated, and the consequences to the tail stability towards on tearing are discussed. It is concluded that the mapping of the brightening region in the auroral oval corresponds to the regions in the tail where the current sheet may be unstable towards ion tearing.

Pulkkinen, T. I.↗