Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure 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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Finite-Time Ensemble Method for Mixed Layer Model Comparison

Here, this work evaluates the fidelity of various upper-ocean turbulence parameterizations subject to realistic monsoon forcing and presents a finite-time ensemble vector (EV) method to better manage the design and numerical principles of these parameterizations. The EV method emphasizes the dynamics of a turbulence closure multimodel ensemble and is applied to evaluate 10 different ocean surface boundary layer (OSBL) parameterizations within a single-column (SC) model against two boundary layer large-eddy simulations (LES). Both LES include realistic surface forcing, but one includes wind-driven shear turbulence only, while the other includes additional Stokes forcing through the wave-average equations that generate Langmuir turbulence. The finite-time EV framework focuses on what constitutes the local behavior of the mixed layer dynamical system and isolates the forcing and ocean state conditions where turbulence parameterizations most disagree. Identifying disagreement provides the potential to evaluate SC models comparatively against the LES. Observations collected during the 2018 monsoon onset in the Bay of Bengal provide a case study to evaluate models under realistic and variable forcing conditions. The case study results highlight two regimes where models disagree 1) during wind-driven deepening of the mixed layer and 2) under strong diurnal forcing.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Convex set and linear mixing model

A major goal of optical remote sensing is to determine surface compositions of the earth and other planetary objects. For assessment of composition, single pixels in multi-spectral images usually record a mixture of the signals from various materials within the corresponding surface area. In this report, we introduce a closed and bounded convex set as a mathematical model for linear mixing. This model has a clear geometric implication because the closed and bounded convex set is a natural generalization of a triangle in n-space. The endmembers are extreme points of the convex set. Every point in the convex closure of the endmembers is a linear mixture of those endmembers, which is exactly how linear mixing is defined. With this model, some general criteria for selecting endmembers could be described. This model can lead to a better understanding of linear mixing models.

Xu, P.↗
DOE OSTI2021

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗
NASA NTRS1982

Turbulent transport modelling of separating and reattaching shear flows

The improvement of capabilities for computer simulation of turbulent recirculating flows was investigated. Attention has been limited to two dimensional flows and principally to statistically stationary motion. Improvement of turbulence modeling explored the treatment of the near wall sublayer and of the exterior fully turbulent region, working within the framework of turbulence closures requiring the solution of transport equations for the turbulence energy and its dissipation rate. The work on the numerical procedure, based on the Gosman-Pun program TEACH, addressed the problems of incorporating the turbulence model as well as the extension to time dependent flows, the incorporation of a third order approximation of convective transport, and the treatment of non-orthogonal boundaries.

Launder, B. E.↗
NASA NTRS1995

Lead-induced atmospheric circulations

Leads provide a significant source of heat and moisture to the Arctic winter atmosphere, and plumes from wide leads have been observed to penetrate the Arctic inversion. We have developed a two-dimensional, high-resolution, deep anelastic numerical model to investigate the atmospheric convection from leads with widths ranging from 100 m to 10 km. A second-order turbulence closure scheme is used to parameterize the atmospheric turbulence in the horizontally inhomogeneous system. This study describes how the lead-induced circulations can enhance the vertical transport of heat into the atmospheric boundary layer. This model is compared with large-eddy simulation results and with lidar observations of a lead-induced ice crystal plume. The model is used to study the effect of varying lead widths and ambient atmospheric conditions on the resultant convection from leads, and some preliminary results are described.

Alam, Afshan↗
NASA NTRS2007

Panel-Stiffener Debonding and Analysis Using a Shell/3D Modeling Technique

A shear loaded, stringer reinforced composite panel is analyzed to evaluate the fidelity of computational fracture mechanics analyses of complex structures. Shear loading causes the panel to buckle. The resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. The panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot, web and noodle as well as the panel skin near the delamination front were modeled with a local 3D solid model. Across the width of the stringer foot, the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. The objective was to study the effect of the fidelity of the local 3D finite element model on the computed mixed-mode strain energy release rates and the failure index.

Krueger, Ronald↗
DOE OSTI2023

Assessing turbulence and mixing parameterizations in the gray-zone of multiscale simulations over mountainous terrain during the METEX21 field experiment

Multiscale numerical weather prediction models transition from mesoscale, where turbulence is fully parameterized, to microscale, where the majority of highly energetic scales of turbulence are resolved. The turbulence gray-zone is situated between these two regimes and multiscale models must downscale through these resolutions. Here, we compare three multiscale simulations which vary by the parameterization used for turbulence and mixing within the gray-zone. The three parameterizations analyzed are the Mellor-Yamada Nakanishi and Niino (MYNN) Level 2.5 planetary boundary layer scheme, the TKE-1.5 large eddy simulation (LES) closure scheme, and a recently developed three-dimensional planetary boundary layer scheme based on the Mellor-Yamada model. The simulation domain includes complex (i.e., mountainous) terrain in Nevada that was instrumented with meteorological towers, profiling and scanning lidars, a tethered balloon, and a surface flux tower. Simulations are compared to each other and to observations, with assessment of model skill at predicting wind speed, wind direction and TKE, and qualitative evaluations of transport and dispersion of smoke from controlled releases. This analysis demonstrates that microscale predictions of transport and dispersion can be significantly influenced by the choice of turbulence and mixing parameterization in the terra incognita , particularly over regions of complex terrain and with strong local forcing. This influence may not be apparent in the analysis of model skill, and motivates future field campaigns involving controlled tracer releases and corresponding modeling studies of the turbulence gray-zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Catalytic Autoxidation for Depolymerization of Multilayer Plastic Films

Recycling multilayer plastic films is challenged by a diversity of polymers, prompting development of new recycling methods. For the depolymerization of mixed polymers like those in multilayer films, metal-catalyzed autoxidation offers a versatile chemical recycling method to deconstruct multiple polymers to useful oxygenates. Here, we demonstrate that catalytic autoxidation is effective for depolymerizing multilayer films across diverse chemistries. We investigated conditions for a model polyethylene substrate using a Co, Mn, and Br cocatalyst system, achieving full carbon closure with oxygenated small molecules contributing up to 48 mol% carbon. Subsequently, we characterized product distributions for several common polymers used in multilayer films using high-resolution mass spectrometry (HRMS) and developed analytical methods to quantify the resulting complex product streams. Optimized conditions for polyethylene were applied to 11 multilayer plastic films containing 10 different polymers, including films with nonpolymeric potential disrupters like aluminum foil and titanium dioxide, showing that catalytic autoxidation is effective across a broad range of polymer types and is resistant to disrupters and additives. The generation of CO 2 in these reactions overall suggests that both reaction engineering and modifications to the reaction conditions will be required to achieve higher yields of soluble oxygenated products.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nonadiabatic Dynamics of Photoexcited cis-Stilbene Using Ab Initio Multiple Spawning

The photochemistry of cis-stilbene proceeds through two pathways: cis-trans isomerization and ring closure to 4a,4b-dihydrophenanthrene (DHP). Despite serving for many decades as a model system for photoisomerization, the photodynamics of cis-stilbene is still not fully understood. In this work, we use ab initio multiple spawning on a SA-2-CASSCF(2,2) potential energy surface to simulate the nonadiabatic dynamics of isolated cis-stilbene. We find the cyclization (to DHP and cis-stilbene) and isomerization (to trans- and cis-stilbene) reaction coordinates to be orthogonal; branching between the two pathways is determined on the S1 excited state within 150 fs of photoexcitation. Furthermore ,trajectory basis functions (TBFs) undergoing cyclization decay rapidly to the ground state in 250 fs, while TBFs moving along the isomerization coordinate remain on the excited state longer, with the majority decaying between 300 and 500 fs. We observe three avoided crossing regions in the dynamics: two along the isomerization coordinate (displaying pyramidalization and migration of an ethylenic hydrogen or phenyl group), and one DHP-like conical intersection along the cyclization coordinate. The isomeric form of the vibrationally hot photoproducts (as determined by measurement 2 ps after photoexcitation) is determined within less than 50 fs of decay to the ground state mediated by passage through a conical intersection. Excess vibrational energy of ground state cis- and trans-stilbene is channelled into phenyl torsions (with mostly opposing directionality). Finally, our simulations are validated by direct comparison to experiment for the absorption spectrum, branching ratio of the three photoproducts (44:52:4 cis-stilbene:trans-stilbene:DHP), and excited state lifetime (520 ± 40 fs).

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

Feasibility Study for Disposal Control Rod Assemblies Using UNF-ST&DARDS As-Loaded Zion Dual Purpose Cask Models

This report documents an initial evaluation to support future use of disposal control rod assemblies (DCRAs) for post-closure criticality control in dual purpose canisters (DPCs). The work described herein is an extension of previous efforts performed by Walker using inputs generated by the Used Nuclear Fuel-Storage, Transportation & Disposal Analysis Resource and Data System (UNFST&DARDS) for the Zion site DPCs with as-loaded isotopic compositions. The results of this analysis demonstrate that there are multiple pathways to using DCRA for post-closure criticality control. Various configurations of DCRA material, diameter, number of rods per DCRA, and number and location of DCRAs within a DPC were shown to be effective in varying degrees for the set of DPCs analyzed by Walker. Because of the variations in DPC as-loaded isotopic compositions considered in an array of DCRA parametric sweeps, it can be concluded that a DPC-specific methodology is feasible (i.e., a one-size-fits-all approach may not be needed). Instead, the utility program created for this work can be expanded to develop capabilities to provide DPC-specific DCRA arrangements to limit cost and weight and to allow for operational considerations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Nuclide Inventory Benchmark for BWR Spent Nuclear Fuel: Challenges in Evaluation of Modeling Data Assumptions and Uncertainties

This work discusses challenges and approaches to uncertainty analyses associated with the development of a nuclide inventory benchmark for fuel irradiated in a boiling water reactor. The benchmark under consideration is being developed based on experimental data from the SFCOMPO international database. The focus herein is on how to address missing data in fuel design and operating conditions that are important for adequately simulating the time-dependent changes in fuel during irradiation in the reactor. The effects of modeling assumptions and uncertainties in modeling parameters on the calculated nuclide inventory were analyzed and quantified through computational models developed using capabilities in the SCALE code system. Particular attention was given to the impact of the power history and water coolant density on the calculated nuclide inventory, as well as to the effect of geometry modeling considerations not usually addressed in a nuclide inventory benchmark. These considerations include gap closure, channel bow, and channel corner radius, which do not usually apply to regular reactor operation but are relevant for assessing impacts of potential anomalous operating scenarios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

On the Existence of the Plateau Emission in High-energy Gamma-Ray Burst Light Curves Observed by Fermi-LAT

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) shows long-lasting high-energy emission in many gamma-ray bursts (GRBs), similar to X-ray afterglows observed by the Neil Gehrels Swift Observatory (Swift). Some LAT light curves (LCs) show a late-time flattening reminiscent of X-ray plateaus. We explore the presence of plateaus in LAT temporally extended emission analyzing GRBs from the second Fermi-LAT GRB Catalog from 2008 to 2016 May with known redshifts, and check whether they follow closure relations corresponding to four distinct astrophysical environments predicted by the external forward shock model. We find that three LCs can be fit by the same phenomenological model used to fit X-ray plateaus and show tentative evidence for the existence of plateaus in their high-energy extended emission. The most favorable scenario is a slow-cooling regime, whereas the preferred density profile for each GRBs varies from a constant-density interstellar medium to an r –2 wind environment. We also compare the end time of the plateaus in γ-rays and X-rays using a statistical comparison with 222 Swift GRBs with plateaus and known redshifts from 2005 January to 2019 August. Within this comparison, the case of GRB 090510 shows an indication of chromaticity at the end time of the plateau. Finally, we update the 3D fundamental plane relation among the rest-frame end time of the plateau, its correspondent luminosity, and the peak prompt luminosity for 222 GRBs observed by Swift. We find that these three LAT GRBs follow this relation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Aerosol-induced closure of marine cloud cells: enhanced effects in the presence of precipitation

Abstract. The Weather Research Forecasting (WRF) version 4.3 model is configured within a Lagrangian framework to quantify the impact of aerosols on evolving cloud fields. Kilometer-scale simulations utilizing meteorological boundary conditions are based on 10 case study days offering diverse meteorology during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA). Measurements from aircraft, the ground-based Atmosphere Radiation Measurement (ARM) site at Graciosa Island in the Azores, and A-Train and geostationary satellites are utilized for validation, demonstrating good agreement with the WRF-simulated cloud and aerosol properties. Higher aerosol concentration leads to suppressed drizzle and increased cloud water content in all case study days. These changes lead to larger radiative cooling rates at cloud top, enhanced vertical velocity variance, and increased vertical and horizontal wind speed near the base of the lower-tropospheric inversion. As a result, marine cloud cell area expands, narrowing the gap between shallow clouds and increasing cloud optical thickness, liquid water content, and the top-of-atmosphere outgoing shortwave flux. While similar aerosol effects are observed in lightly to non-raining clouds, they tend to be smaller by comparison. These simulations show a relationship between cloud cell area expansion and the radiative adjustments caused by liquid water path and cloud fraction changes. The adjustments are positive and scale as 74 % and 51 %, respectively, relative to the Twomey effect. While higher-resolution large-eddy simulations may provide improved representation of cloud-top mixing processes, these results emphasize the importance of addressing mesoscale cloud-state transitions in the quantification of aerosol radiative forcing that cannot be attained from traditional climate models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1976

Second-order closure analysis of turbulent reacting flows

A reacting shear layer program for a complete second-order closure of turbulent reacting flows is developed. The program uses a typical eddy box model and many of the fluid mechanical turbulence models developed in earlier studies. The computer program actually solves a total of 30 equations, where the numerical integration of the equations is performed by a forward-time-centered space quasi-implicit upwind finite-difference scheme. Results for nonreacting and simple isothermal reacting flows are presented along with preliminary results for propane-air flames. It is shown that the program predictions for simple nonreacting flows are in good agreement with available data. The calculations for isothermal reacting flows demonstrate the importance of including a mixedness correlation in finite rate chemistry calculations, where significant errors in the flowfield properties can result from neglect of the interaction between turbulence and the chemistry. Results for propane-air diffusion flames further support this conclusion.

Varma, A. K.↗
NASA NTRS1979

The diffuse X-ray background spectrum from 3 to 50keV

The spectrum of the extragalactic diffuse X-ray background was measured with the GSFC cosmic X-ray experiment on HEAO-1 for regions of the sky away from known point sources and more than 20 deg from the galactic plane. A total exposure of 80 sq m-sec-sr is available at present. Free-free emission from an optically thin plasma of 40 plus or minus 5 keV provides an excellent description of the observed spectrum from 3 to 50 keV. This spectral shape is confirmed by measurements from 5 separate layers of three independent detectors. With an estimated absolute precision of about 10 percent, the intensity of the emission at 10 keV is 3.2 keV/keV-sq cm-sec-sr, a value consistent with the average of previously reported spectra. No other spectral features, such as iron line emission, are evident. This spectrum is not typical of known extragalactic objects. A uniform hot intergalactic medium of approximately 36 percent of the closure density of the universe would produce such a flux, although non-uniform models indicating less total matter are probably more realistic.

Marshall, F. E.↗
NASA NTRS1981

The potential role of aerobic biological waste treatment in regenerative life support systems

The purpose of the paper is to make a preliminary assessment of the feasibility of using aerobic biological waste treatment in closed systems. Issues that are addressed in this paper are: (1) how high a degree of material balance is possible, (2) how much might such a system weigh, and (3) how would system closure and weight be affected if animals were included in the system. A computer model has been developed to calculate for different scenarios the compositions and amounts of the streams entering or leaving the waste treatment system and to estimate the launch weight of such a system. A bench scale apparatus has been built to mimic the proposed waste treatment system; the experiments are used to verify model predictions and to improve model parameter estimations.

Shuler, M. L.↗
NASA NTRS1990

Corrosion fatigue crack propagation in metals

This review assesses fracture mechanics data and mechanistic models for corrosion fatigue crack propagation in structural alloys exposed to ambient temperature gases and electrolytes. Extensive stress intensity-crack growth rate data exist for ferrous, aluminum and nickel based alloys in a variety of environments. Interactive variables (viz., stress intensity range, mean stress, alloy composition and microstructure, loading frequency, temperature, gas pressure and electrode potential) strongly affect crack growth kinetics and complicate fatigue control. Mechanistic models to predict crack growth rates were formulated by coupling crack tip mechanics with occluded crack chemistry, and from both the hydrogen embrittlement and anodic dissolution/film rupture perspectives. Research is required to better define: (1) environmental effects near threshold and on crack closure; (2) damage tolerant life prediction codes and the validity of similitude; (3) the behavior of microcrack; (4) probes and improved models of crack tip damage; and (5) the cracking performance of advanced alloys and composites.

Gangloff, Richard P.↗
NASA NTRS1992

Effect of afterbody geometry on aerodynamic characteristics of isolated nonaxisymmetric afterbodies at transonic Mach numbers

A parametric study was conducted in the Langley 16-Foot Transonic Tunnel on an isolated nonaxisymmetic fuselage model that simulates a twin-engine fighter. The effects of aft-end closure distribution (top/bottom) nozzle-flap boattail angle versus nozzle-sidewall boattail angle) and afterbody and nozzle corner treatment (sharp or radius) were investigated. Four different closure distributions with three different corner radii were tested. Tests were conducted over a range of Mach numbers from 0.40 to 1.25 and over a range of angles of attack from -3 to 9 degrees. Solid plume simulators were used to simulate the jet exhaust. For a given closure distribution in the range of Mach numbers tested, the sharp-corner nozzles generally had the highest drag, and the 2-in. corner-radius nozzles generally had the lowest drag. The effect of closure distribution on afterbody drag was highly dependent on configuration and flight condition.

Bangert, Linda S.↗