Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure model”

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 OSTI2021

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Early Signal of Emerging Nuclear Collectivity in Neutron-Rich 129Sb

Radioactive 129Sb, which can be treated as a proton plus semimagic 128Sn core within the particle-core coupling scheme, was studied by Coulomb excitation. Reduced electric quadrupole transition probabilities, B(E2), for the 2+Ⓧπg7/2 multiplet members and candidate πd5/2 state were measured. The results indicate that the total electric quadrupole strength of 129Sb is a factor of 1.39(11) larger than the 128Sn core, which is in stark contrast to the expectations of the empirically successful particle-core coupling scheme. Shell-model calculations performed with two different sets of nucleon-nucleon interactions suggest that this enhanced collectivity is due to constructive quadrupole coherence in the wave functions stemming from the proton-neutron residual interactions, where adding one nucleon to a core near a double-shell closure can have a pronounced effect. The enhanced electric quadrupole strength is an early signal of the emerging nuclear collectivity that becomes dominant away from the shell closure.

Gray, Tim↗
DOE OSTI2025

Charge Radii Measurements of Exotic Tin Isotopes in the Proximity of 𝑁 =50 and 𝑁 =82

We report nuclear charge radii for the isotopes 104–134 Sn , measured using two different collinear laser spectroscopy techniques at ISOLDE-CERN. These measurements clarify the archlike trend in charge radii along the isotopic chain and reveal an odd-even staggering that is more pronounced near the 𝑁 =50 and 𝑁 =82 shell closures. The observed local trends are well described by both nuclear density functional theory and valence space in-medium similarity renormalization group calculations. Both theories predict appreciable contributions from beyond-mean-field correlations to the charge radii of the neutron-deficient tin isotopes. The models, however, fall short of reproducing the magnitude of the known 𝐵⁡(𝐸⁢2) transition probabilities, highlighting the remaining challenges in achieving a unified description of both ground-state properties and collective phenomena.

laser spectroscopy↗
DOE OSTI2024

Predicting the Slowing of Stellar Differential Rotation by Instability-driven Turbulence

Abstract Differentially rotating stars and planets transport angular momentum (AM) internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich–Schubert–Fricke instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing AM gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or “parasitic saturation” models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline-instability-driven turbulence in certain regimes.

Astronomy & Astrophysics↗
DOE OSTI2025

Charge radii measurements of exotic tin isotopes in the proximity of $N=50$ and $N=82$

We report nuclear charge radii for the isotopes $^{104-134}$Sn, measured using two different collinear laser spectroscopy techniques at ISOLDE-CERN. These measurements clarify the archlike trend in charge radii along the isotopic chain and reveal an odd-even staggering that is more pronounced near the $N=50$ and $N=82$ shell closures. The observed local trends are well described by both nuclear density functional theory and valence space in-medium similarity renormalization group calculations. Both theories predict appreciable contributions from beyond-mean-field correlations to the charge radii of the neutron-deficient tin isotopes. The models fall short, however, of reproducing the magnitude of the known $B(E2)$ transition probabilities, highlighting the remaining challenges in achieving a unified description of both ground-state properties and collective phenomena.

FOS: Physical sciences↗
DOE OSTI2026

Northwest Colorado Local Energy Analysis: October 2025 Results

The U.S. Department of Energy's (DOE) Communities Local Energy Action Program (Communities LEAP) provided customized technical assistance to the Northwest Colorado region to assist in developing strategies that facilitate long-term economic benefits through locally driven energy planning and development. The region sought technical assistance to model alternatives energy projects in response to the expected retirement of coal-fired powers plant and associated coal mines by 2030. Northwest Colorado faces a potential economic gap due to the loss of tax revenue, labor income, and jobs within the region due to the impending closure of the Craig and Hayden Power Stations and related mining operations at Trapper and Colowyo mines.

08 HYDROGEN↗
NASA NTRS1977

A second-order closure analysis of turbulent diffusion flames

A complete second-order closure computer program for the investigation of compressible, turbulent, reacting shear layers was developed. The equations for the means and the second order correlations were derived from the time-averaged Navier-Stokes equations and contain third order and higher order correlations, which have to be modeled in terms of the lower-order correlations to close the system of equations. In addition to fluid mechanical turbulence models and parameters used in previous studies of a variety of incompressible and compressible shear flows, a number of additional scalar correlations were modeled for chemically reacting flows, and a typical eddy model developed for the joint probability density function for all the scalars. The program which is capable of handling multi-species, multistep chemical reactions, was used to calculate nonreacting and reacting flows in a hydrogen-air diffusion flame.

Varma, A. K.↗
NASA NTRS1982

The JT9D Jet Engine Diagnostics Program

The various engine deterioration phenomena that affect JT9D performance retention were studied, and approaches to improve performance retention of engines were identified. The program included surveys of historical data, monitoring of in service engines, ground and flight testing of instrumented engines, analysis, and analytical modeling. Performance deterioration is made up of both short and long term modes, both of which are flight cycle related phenomena. Short term deterioration occurs primarily during airplane acceptance testing prior to delivery to the airline. This effect is caused by flight load and power induced clearance closures and engine deflections with resulting rubbing of airfoils and seals. Long term deterioration is caused by erosion of airfoils and gas path seals during ground operation and take off and by cyclic induced thermal distortion of the high pressure turbine airfoils. Studies of possible remedial approaches have shown that performance retention within 1 to 2 percent of initial revenue service performance can be achieved with a proper program of hot section and cold section maintenance.

Olsson, W. J.↗
DOE OSTI2021

System Model Calculations for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

Two lined trenches (Trench 31 and 34) located within the 200 West Low-Level Burial Ground (LLBG) area of the U.S. Department of Energy’s (DOE) Hanford Site are designated for permanent disposal of low-level radioactive waste (LLW) and mixed low-level radioactive waste (MLLW). In accordance with DOE O 435.1, Radioactive Waste Management, radioactive waste shall be managed and disposed in a manner that is protective of worker, public health and safety, and the environment. DOE O 435.1 requires a site-specific radiological performance assessment (PA) that includes calculations of potential releases and subsequent doses to members of the public for a period of 1,000 years after closure of a low-level waste (LLW) disposal facility. The purpose of this environmental calculation file (ECF) is to document the various dose calculations performed to support DOE’s LLBG PA. The dose assessments are performed to evaluate the potential exposure of an all-pathway representative person to radionuclide contaminants of potential concern (COPCs) that may be released from the LLBG Trenches 31 and 34 to the point of calculation (POC) located at the outer edge of a 100 m buffer zone surrounding the trench boundary. This ECF uses inputs from other ECFs and their associated model package reports (MPRs) supporting the LLBG PA. The five major objectives of this ECF are to present the results of the following: (1) Deterministic dose assessments for both groundwater and atmospheric exposure scenarios; (2) Radon-222 (Rn-222) flux emanating at the surface from the disposed wastes; (3) Uncertainty analyses of the groundwater pathway annual dose due to uncertainty in input parameters; (4) Sensitivity analyses of the input parameters that could potentially impact the fate and transport of the contaminants and doses; (5) Doses arising from inadvertent intruder exposure scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1982

V/STOL Tandem Fan transition section model test

An approximately 0.25 scale model of the transition section of a tandem fan variable cycle engine nacelle was tested in the NASA Lewis Research Center 10-by-10 foot wind tunnel. Two 12-inch, tip-turbine driven fans were used to simulate a tandem fan engine. Three testing modes simulated a V/STOL tandem fan airplane. Parallel mode has two separate propulsion streams for maximum low speed performance. A front inlet, fan, and downward vectorable nozzle forms one stream. An auxilliary top inlet provides air to the aft fan - supplying the core engine and aft vectorable nozzle. Front nozzle and top inlet closure, and removal of a blocker door separating the two streams configures the tandem fan for series mode operations as a typical aircraft propulsion system. Transition mode operation is formed by intermediate settings of the front nozzle, blocker door, and top inlet. Emphasis was on the total pressure recovery and flow distortion at the aft fan face. A range of fan flow rates were tested at tunnel airspeeds from 0 to 240 knots, and angles-of-attack from -10 to 40 deg for all three modes. In addition to the model variables for the three modes, model variants of the top inlet were tested in the parallel mode only. These lip variables were: aft lip boundary layer bleed holes, and Three position turning vane. Also a bellmouth extension of the top inlet side lips was tested in parallel mode.

Simpkin, W. E.↗
NASA NTRS2006

Remote Sensing Reflectance and Inherent Optical Properties in the Mid-mesohaline Chesapeake Bay

We used an extensive set of bio-optical data and radiative transfer (RT) model simulations of radiation fields to investigate relationships between inherent optical properties and remotely sensed quantities in the optically complex, mid-mesohaline Chesapeake Bay waters. Field observations showed that the chlorophyll algorithms used by the MODIS (MODerate resolution Imaging Spectroradiometer) ocean color sensor (i.e. Chlor_a, chlor_MODIS, chlor_a_3 products) do not perform accurately in these Case 2 waters. This is because, when applied to waters with high concentrations of chlorophyll, all MODIS algorithms are based on empirical relationships between chlorophyll concentration and blue-green wavelength remote sensing reflectance (Rrs) ratios that do not account for the typically strong blue-wavelength absorption by non-covarying, dissolved and non-algal particulate components. Stronger correlation was observed between chlorophyll concentration and Rrs ratios in the red (i.e. Rrs(677)/Rrs(554)) where dissolved and non-algal particulate absorption become exponentially smaller. Regionally-specific algorithms that are based on the phytoplankton optical properties in the red wavelength region provide a better basis for satellite monitoring of phytoplankton blooms in these Case 2 waters. Good optical closure was obtained between independently measured Rrs spectra and the optical properties of backscattering, b(sub b), and absorption, a, over the wide range of in-water conditions observed in the Chesapeake Bay. Observed variability in the quantity f/Q (proportionality factor in the relationship between Rrs and the water inherent optical properties ratio b(sub b)/(a+b(sub b)) was consistent with RT model calculations for the specific measurement geometry and water bio-optical characteristics. Data and model results showed that f/Q values in these Case 2 coastal waters are not considerably different from those estimated in previous studies for Case 1 waters. Variation in surface backscattering significantly affected Rrs magnitude across the visible spectrum and was most strongly correlated (R(sup 2)=0.88) with observed variability in Rrs at 670 nm. Surface values of particulate backscattering were strongly correlated with non-algal particulate absorption, a(sub nap), in the blue wavelengths (R(sup 2)=0.83). These results, along with the measured values of backscattering fraction magnitude and non-algal particulate absorption spectral slope, suggest that suspended non-algal particles with high inorganic content are the major water constituents regulating b(sub b) variability in the mid-mesohaline Chesapeake Bay. Remote retrieval of surface b(sub b) and (a(sub nap), from Rrs(670) can be used in regionally-specific satellite algorithms to separate contribution by non-algal particles and dissolved organic matter to total light absorption in the blue, and monitor non-algal suspended particle concentration and distribution in these Case 2 waters.

Tzortziou, Maria↗
NASA NTRS1994

Radiation and scattering from cylindrically conformal printed antennas

Microstrip patch antennas offer considerable advantages in terms of weight, aerodynamic drag, cost, flexibility, and observables over more conventional protruding antennas. These flat patch antennas were first proposed over thirty years ago by Deschamps in the United States and Gutton and Baisinot in France. Such antennas have been analyzed and developed for planar as well as curved platforms. However, the methods used in these designs employ gross approximations, suffer from extreme computational burden, or require expensive physical experiments. The goal of this thesis is to develop accurate and efficient numerical modeling techniques which represent actual antenna structures mounted on curved surfaces with a high degree of fidelity. In this thesis, the finite element method is extended to cavity-backed conformal antenna arrays embedded in a circular, metallic, infinite cylinder. Both the boundary integral and absorbing boundary mesh closure conditions will be used for terminating the mesh. These two approaches will be contrasted and used to study the scattering and radiation behavior of several useful antenna configurations. An important feature of this study will be to examine the effect of curvature and cavity size on the scattering and radiation properties of wraparound conformal antenna arrays.

Kempel, Leo C.↗
NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗
NASA NTRS2018

Columbus IFHX Ammonia Leak Analysis

After the Columbus Moderate Temperature Loop (MTL) InterFace Heat eXchanger (IFHX) low temperature event of GMT 345-2013, NASA investigated relevant transient scenarios involving IFHX rupture after water freezing and subsequent thawing. NASA recommended development of a Fault Detection Isolation and Recovery (FDIR) plan that would, in the event of a heat exchanger freeze event, close the Water On/Off Valves (WOOVs) to isolate the heat exchanger and prevent ammonia from the external flow loops from spreading into the cabin. NASA performed a preliminary simplified analysis for the reference case of IFHX rupture, but for a deeper understanding TAS developed detailed SINDA-FLUINT models of the Columbus ITCS that were built and run through the SINAPS GUI. This allowed simulation of the ammonia leakage physics including the variation of environmental parameters, thus providing more accurate and specific input to the FDIR under development. The result was finalization of the IFHX WOOVs closure sequence and wait times to contain the ammonia propagation to Columbus and allow identification of the leaking IFHX. In addition, the analysis results provided reference pressure profiles to be used on console and by the Engineering as support for the telemetry data assessment in case of failure.This paper gives an overview on the issue and focuses on the analytical aspects of the multiphase fluid dynamics involved.

Tilloca, A.↗
NASA NTRS2020

Radiative Heating from Biomass Burning Aerosol and its Impact on Cloud Structure in the Southeast Atlantic

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. In many cases, the complex physical relationships between marine boundary cloud morphology and the environmental conditions in which the clouds exist are not well understood. Such uncertainties arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean. It is likely that the absorbing and heating properties of these aerosols influence the microphysical composition and macrophysical arrangement of marine stratocumulus and trade cumulus in this region; however, this has yet to be quantified. The deployment of the Atmospheric Radiation Measurement Mobile Facility #1 (AMF1) in support of LASIC (Layered Atlantic Smoke Interactions with Clouds) provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96 S, 14.35 W). Thermodynamic profiles will be analyzed through the unique combination of sounding data from radiosonde launches and microwave profiling radiometers, giving observations of additional quantities important for cloud development such as CAPE and CIN at a fine temporal resolution. The thermodynamic profiles will be presented in conjunction with detailed observations of the cloud structure over the site from a K-band cloud radar, micropulse lidar, and laser ceilometer. The observed thermodynamic and cloud profiles will be used as input forcing, alongside aerosols from the Modern Era Retrospective analysis for Research and Applications, version 2 (MERRA-2), for the Rapid Radiative Transfer Model (RRTM) to gain information regarding the radiative heating profiles. Idealized experiments using RRTM with and without aerosols will be used to quantify the impact of biomass burning carbonaceous aerosol plumes as they pass over the site. Due to documented discrepancies in the single scatter albedo (SSA) between models and observations, further sensitivity experiments will demonstrate the importance of the optical properties of biomass burning aerosol in accurately representing heating within the column. Finally, the heating rates will be put into context of the cloud structure over the site from the perspective of the mass flux closure from the University of Washington shallow convective scheme.

Collow, Allison↗
DOE OSTI2023

MRT 7365: Power flow physics and key physics phenomena

The Z accelerator at Sandia National Laboratories conducts z-pinch experiments at 26 MA in support of DOE missions in stockpile stewardship, dynamic materials, fusion, and other basic sciences. Increasing the current delivered to the z-pinch would extend our reach in each of these disciplines. To achieve increases in current and accelerator efficiency, a fraction of Z’s shots are set aside for research into transmission-line power flow. These shots, with supporting simulations and theory, are incorporated into this Advanced Diagnostics milestone report. The efficiency of Z is reduced as some portion of the total current is shunted across the transmission-line gaps prior to the load. This is referred to as “current loss”. Electrode plasmas have long been implicated in this process, so the bulk of dedicated power-flow experiments are designed to measure the plasma environment. The experimental analyses are enhanced by simulations conducted using realistic hardware and Z voltage pulses. In the same way that diagnostics are continually being improved for sensitivity and resolution, the modeling capability is continually being improved to provide faster and more realistic simulations. The specifics of the experimental hardware, diagnostics, simulations, and algorithm developments are provided in this report. The combined analysis of simulation and data confirms that electrode plasmas have the most detrimental impact on current delivery. Experiments over the last three years have tested the theoretical current-loss mechanisms of enhanced ion current, plasma gap closure, and Hall-related current. These mechanisms are not mutually exclusive and may be coincident in the final feed as well as in upstream transmission lines. The final-feed geometries tested here, however, observe lower-density plasmas without dominant ion currents which is consistent with a Hall-related current. The picture of plasma formation and transport formed from experiment and simulation is informing hardware designs being fielded on Z now and being proposed for the Next-Generation Pulsed Power (NGPP) facility. In this picture, the strong magnetic fields that heat the electrodes above particle emission thresholds also confine the charged particles near the surface. Some portion of the plasmas thus formed is transported into the transmission-line gap under the force of the electric field, with aid from plasma instabilities. The gap plasmas are then transported towards the load by a cross-field drift, where they accumulate and contribute to a likely Hall-related cross-gap current. The achievements in experimental execution, model validation, and physical analysis presented in this report set the stage for continued progress in power flow and load diagnostics on Z. The planned shot schedule for Z and Mykonos will provide data for extrapolation to higher current to ensure the predicted performance and efficiency of a NGPP facility.

43 PARTICLE ACCELERATORS↗
NASA NTRS1990

A locally interactive laminar separation bubble model

A laminar separation bubble model has been developed which approximates the viscous/inviscid interaction in the bubble region. An implicit transition criterion allows the accurate prediction of the bubble length and of the shear layer growth along the bubble. Recent insights into the bubble flow field have uncovered the need to modify existing turbulent closure correlations between the boundary-layer variables. As the experimental data base necessary for their development is lacking, possible modifications to correlations in use are proposed and preliminary results discussed.

Maughmer, Mark D.↗