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,405 records · Page 78

Initial Simulations of Empty Room Collapse and Reconsolidation at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is a geologic repository for defense-related nuclear waste. If left undisturbed, the virtually impermeable rock salt surrounding the repository will isolate the nuclear waste from the biosphere. If humans accidentally intrude into the repository in the future, then the likelihood of a radionuclide release to the biosphere will depend significantly on the porosity and permeability of the repository itself. Room ceilings and walls at the WIPP tend to collapse over time, causing rubble piles to form on floors of empty rooms. The surrounding rock formation will gradually compact these rubble piles until they eventually become solid salt, but the length of time for a rubble pile to reach a certain porosity and permeability is unknown. This report details the first efforts to build models to predict the porosity and permeability evolution of an empty room as it closes. Conventional geomechanical numerical methods would struggle to model empty room collapse and rubble pile consolidation, so three different meshless methods, the Immersed Isogeometric Analysis Meshfree, Reproducing Kernel Particle Method (RKPM), and the Conformal Reproducing Kernel method, were assessed. First, the meshless methods and the finite element method each simulated gradual room closure, without ceiling or wall collapse. All three methods produced equivalent room closure predictions with comparable computational speed. Second, the Immersed Isogeometric Analysis Meshfree method and RKPM simulated two-dimensional empty room collapse and rubble pile consolidation. Both methods successfully simulated large viscoplastic deformations, fracture, and rubble pile rearrangement to produce qualitatively realistic results. In addition to geomechanical simulations, the flow channels in damaged salt and crushed salt were measured using micro-computed tomography, and input into a computational fluid dynamics simulation to predict the salt's permeability. Although room for improvement exists, the current simulation approaches appear promising.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Fatigue Life Analysis of Tapered Hybrid Composite Flexbeams

Nonlinear-tapered flexbeam laminates from a full-size composite helicopter rotor hub flexbeam were tested under combined constant axial tension and cyclic bending loads. The two different graphite/glass hybrid configurations tested under cyclic loading failed by delamination in the tapered region. A 2-D finite element model was developed which closely approximated the flexbeam geometry, boundary conditions, and loading. The analysis results from two geometrically nonlinear finite element codes, ANSYS and ABAQUS, are presented and compared. Strain energy release rates (G) obtained from the above codes using the virtual crack closure technique (VCCT) at a resin crack location in the flexbeams are presented for both hybrid material types. These results compare well with each other and suggest that the initial delamination growth from the resin crack toward the thick region of the flexbeam is strongly mode II. The peak calculated G values were used with material characterization data to calculate fatigue life curves and compared with test data. A curve relating maximum surface strain to number of loading cycles at delamination onset compared reasonably well with the test results.

Murri, Gretchen B.↗

Mapping and Monitoring Delmarva Fox Squirrel Habitat Using an Airborne LiDAR Profiler

Twenty five hundred thirty nine kilometers of airborne laser profiling and videography data were acquired over the state of Delaware during the summer of 2000. The laser ranging measurements and video from approximately one-half of that data set (1304 km) were analyzed to identify and locate forested sites that might potentially support populations of Delmarva fox squirrel (DFS, Sciurus niger cinereus). The DFS is an endangered species previously endemic to tall, dense, mature forests with open understories on the Eastern Shore of the Chesapeake Bay. The airborne LiDAR employed in this study can measure forest canopy height and canopy closure, but cannot measure or infer understory canopy conditions. Hence the LiDAR must be viewed as a tool to map potential, not actual, habitat. Fifty-three potentially suitable DFS sites were identified in the 1304 km of flight transect data. Each of the 53 sites met the following criteria according to the LiDAR and video record: (1 ) at least 120m of contiguous forest; (2) an average canopy height greater than 20m; (3) an average canopy closure of >80%; and (4) no roofs, impervious surface (e.g., asphalt, concrete), and/or open water anywhere along the 120m length of the laser segment. Thirty-two of the 53 sites were visited on the ground and measurements taken for a DFS habitat suitability model. Seventy eight percent of the sites (25 of 32) were judged by the model to be suited to supporting a DFS population. Twenty-eight of the 32 sites visited in the field were in forest cover types (hardwood, mixed wood, conifer, wetlands) according to a land cover GIS map. Of these, 23 (82%) were suited to support DFS. The remaining 4 sites were located in nonforest cover types - agricultural or residential areas. Two of the four, or 50% were suited to the DFS. All of the LiDAR flight data, 2539 km, were analyzed to

Nelson, Ross↗

Adaptative Site Management for a 115 Acre Chlorinated Solvent Plume with Two Separate Source Areas at Kennedy Space Center, Florida

Background/Objectives. During Resource Conservation and Recovery Act (RCRA) Facility Investigation (RFI) activities, Geosyntec delineated a chlorinated volatile organic compound (CVOC) plume at the National Aeronautics and Space Administration’s (NASA’s) Vehicle Assembly Building (VAB) area located at KSC, Florida. The RFI activities identified an approximate 115-acre dissolved plume (primarily vinyl chloride) and a trichloroethene (TCE) source area in an active aerospace complex that is surrounded by sensitive wetland/waterbodies. Due to the size of the impacted area, the Corrective Measure Design included a multi-component strategy: (i) address the source area via bioremediation; (ii) protect sensitive wetlands from impacted groundwater discharge via biosparging; and (iii) Long Term Monitoring (LTM) of the remaining dissolved plume. After the Corrective Measures implementation (CMI), NASA and Geosyntec worked with Florida Department of Environmental Protection (FDEP) to implement an adaptive site management for the complex, 115-acre site outside of the traditional RCRA process. The adaptive site management approach relied on performing supplemental assessments and implementing Interim Measures (IMs) to further assess and implement remedies over time while working within site and budget constraints, with an overall goal of achieving enough mass reduction to transition the entire site to LTM and eventually achieve site closure. Approach/Activities. After the biosparge barrier was operational and bioremediation within the source area (referred to as Hot Spot 1) achieved the Corrective Action Objective (CAO), supplemental assessment of the area between Hot Spot 1 and the biosparge barrier was performed. The conceptual site model was updated using the supplemental assessment results and an air sparge system IM was designed to treat an approximate 1.2 acre area (referred to as Hot Spot 2). After installation of the air sparge system, supplemental assessment within the remainder of the 115-acre dissolved plume was performed and a second TCE source area was identified. The TCE source area and associated areas with elevated CVOC concentrations (referred to as Hot Spot 3) were delineated and a bioremediation IM was implemented. Also, the downgradient impacts from Hot Spot 3 were adjacent to a sensitive waterbody, and negotiations with the FDEP allowed the area to be monitored using LTM. Results/Lessons Learned. The performance of supplemental assessment activities and implementation of remedial alternatives as IMs allowed NASA to successfully address groundwater impacts over time, while working within the FDEP regulatory framework. The implementation of the CMI and multiple IMs has achieved the following goals: (i) the biosparge barrier has mitigated the potential discharge of impacted groundwater to an adjacent wetland; (ii) enhanced bioremediation within Hot Spot 1 achieved the CAO within 2 years and transitioned the area into LTM; (iii) operation of an air sparge system within Hot Spot 2 removed TCE as a constituent of concern and contributed to a reduction (approximately 43%) in the impacted groundwater area outside the air sparge treatment area (plume collapse); and (iv) bioremediation within Hot Spot 3 removed approximately 80% of the CVOC mass and contributed to a reduction (approximately 47%) in the impacted groundwater area outside the bioremediation IM treatment area. Overall, the adaptive approach is protecting the sensitive water bodies surrounding the complex site and reducing the area of impacted groundwater, which is moving the entire site towards LTM.

Rebecca C Daprato↗

Data for Propagation Method and Planting Density Influence Canopy Developmental Transition and Biomass Productivity in Miscanthus × giganteus

Understanding how establishment practices influence the mechanisms underlying Miscanthus × giganteus (miscanthus) productivity and canopy development is critical for optimizing management. Data was collected during the juvenile (2011–2013) and mature (2024) phases of a long-term field experiment established in Urbana, Illinois, to evaluate the effects of propagation method (plug propagation [PP] and rhizome propagation [RP]), planting density (1.0, 0.75, and 0.25 plants m⁻²), and nitrogen application (0 and 67 kg N ha⁻¹) on end-of-season biomass yield, tiller mass, tiller density, and tiller height. Linear regression models identified the dominant predictors of yield across stand ages and management regimes. Planting density, nitrogen (N) application, and propagation method significantly influenced early yield and canopy development. During the juvenile phase, biomass yield was driven by tiller density due to canopy expansion; in the mature phase, yield became driven by tiller mass. The PP plots produced higher tiller density than the RP plots, resulting in faster canopy closure and higher juvenile-phase yields. Rhizome-propagated (RP) plots produced lower tiller density, but individual tillers were 3.3–6.4 g tiller−1 heavier than PP tillers. After the canopy reached equilibrium, the PP and RP yields were similar because greater RP tiller mass compensated for its lower tiller density. Higher planting density resulted in greater yield and tiller density during the second year (2012), but this effect was absent from the third year (2013) onward. In the juvenile phase, N fertilization enhanced yield by 1.6–3.4 Mg ha−1. Initiating fertilization in 2013 on unfertilized plots produced biomass similar to that in fertilized plots, suggesting yield recovery in the mature phase. These findings revealed that establishment strategies, including propagation method and planting density, influence juvenile miscanthus canopy development and productivity, transitioning from tiller-density- to mass-dominated yields, but not mature phase productivity.

Miscanthus↗

Diurnal rainfall response to the physiological and radiative effects of CO 2 in tropical forests in the Energy exascale earth system model v1

Understanding how the connection between rainfall and tropical forests will respond to increasing CO$_2$ concentrations is a key element in understanding how the tropical water cycle will respond to increasing CO$_2$. The plant physiological and radiative impacts of CO$_2$ on rainfall patterns over tropical forest regions are examined in the Energy Exascale Earth System Model version 1.1 (E3SMv1.1-BGC) biogeochemistry experiments. Composite analysis reveals a dampening of the diurnal cycle of rainfall over the Amazon, Congo, and Maritime Continent in response to rising CO$_2$ levels, regardless of the sign of total rainfall change. A full factorial model experiment confirms that the CO$_2$ radiative and CO$_2$ plant physiological effects can individually or jointly reduce the magnitude of the rainfall diurnal cycle, though the physical pathway giving rise to the reduction differs between the two effects. For the physiological response, stomatal closure reduces evapotranspiration, which dries the boundary layer and raises the lifting condensation level. These effects combine to reduce deep convective rainfall during its peak occurrence in the late daytime to early nighttime period. For the radiative response, a relative reduction in daytime CAPE (consistent with a reduction in the diurnal temperature range) leads to less frequent triggering of deep convection and a reduction of rainfall diurnal amplitude. These diurnal rainfall changes are structurally similar across seasons, and show little sensitivity to representation of nutrient coupling for the land biogeochemistry. In agreement with previous findings, the physiological response has only minor impact on extreme rainfall relative to the radiative response.

54 ENVIRONMENTAL SCIENCES↗

Proton Distribution Radii of 16–24 O: Signatures of New Shell Closures and Neutron Skin

The root mean square radii of the proton density distribution in 16–24 O derived from measurements of charge changing cross sections with a carbon target at ~900A MeV together with the matter radii portray thick neutron skin for 22–24 O despite 22,24 O being doubly magic. Imprints of the shell closures at N=14 and 16 are reflected in local minima of their proton radii that provide evidence for the tensor interaction causing them. The radii agree with ab initio calculations employing the chiral NNLO sat interaction, though skin thickness predictions are challenged. Furthermore, shell model predictions agree well with the data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Annual Summary Report (FY 2025) Performance Assessment for the Integrated Disposal Facility

The purpose of this Annual Summary Report (ASR) for fiscal year (FY) 2025 is to evaluate the continued adequacy of the Integrated Disposal Facility (IDF) Performance Assessment (PA) and Disposal Authorization Statement (DAS). This report consolidates relevant monitoring data, modeling analyses, and regulatory reviews to demonstrate a reasonable expectation that the PA objectives and performance measures will be met, as required under DOE O 435.1, Radioactive Waste Management. The ASR follows the guidance in DOE-STD-5002-2017, Disposal Authorization Statement and Tank Closure Documentation, which provides a framework for maintaining the validity of the DAS through periodic assessment of facility performance and compliance with waste disposal requirements.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Two-moment Neutrino Flavor Transformation with Applications to the Fast Flavor Instability in Neutron Star Mergers

Abstract Multi-messenger astrophysics has produced a wealth of data with much more to come in the future. This enormous data set will reveal new insights into the physics of core-collapse supernovae, neutron star mergers, and many other objects where it is actually possible, if not probable, that new physics is in operation. To tease out different possibilities, we will need to analyze signals from photons, neutrinos, gravitational waves, and chemical elements. This task is made all the more difficult when it is necessary to evolve the neutrino component of the radiation field and associated quantum-mechanical property of flavor in order to model the astrophysical system of interest—a numerical challenge that has not been addressed to this day. In this work, we take a step in this direction by adopting the technique of angular-integrated moments with a truncated tower of dynamical equations and a closure, convolving the flavor-transformation with spatial transport to evolve the neutrino radiation quantum field. We show that moments capture the dynamical features of fast flavor instabilities in a variety of systems, although our technique is by no means a universal blueprint for solving fast flavor transformation. To evaluate the effectiveness of our moment results, we compare to a more precise particle-in-cell method. Based on our results, we propose areas for improvement and application to complementary techniques in the future.

Astronomy & Astrophysics↗

Robust and High Throughput Tools for Characterization of Wet Cellulosic Biogas Feedstocks and Co-Products: Cooperative Research and Development Final Report, CRADA Number CRD-21-21472

This project will develop standardized analytical procedures to support understanding of the chemical composition of wet manure feedstocks, and the digestate formed following anaerobic digestion of organic waste to produce renewable natural gas (RNG). The developed methods will be validated to improve component mass closure, eliminate ambiguities around key constituents, and document the precision and accuracy of the methods. Following validation of the analytical laboratory methods for characterization of wet organic waste feedstocks and digestate, NLR will use the same samples to build and validate near-infrared (NIR) models for high throughput analysis, potentially deployable at commercial facilities, enabling increased production of low carbon intensity RNG from wet organic wastes.

09 BIOMASS FUELS↗

Discrete crack growth analysis methodology for through cracks in pressurized fuselage structures

A methodology for simulating the growth of long through cracks in the skin of pressurized aircraft fuselage structures is described. Crack trajectories are allowed to be arbitrary and are computed as part of the simulation. The interaction between the mechanical loads acting on the superstructure and the local structural response near the crack tips is accounted for by employing a hierarchical modeling strategy. The structural response for each cracked configuration is obtained using a geometrically nonlinear shell finite element analysis procedure. Four stress intensity factors, two for membrane behavior and two for bending using Kirchhoff plate theory, are computed using an extension of the modified crack closure integral method. Crack trajectories are determined by applying the maximum tangential stress criterion. Crack growth results in localized mesh deletion, and the deletion regions are remeshed automatically using a newly developed all-quadrilateral meshing algorithm. The effectiveness of the methodology and its applicability to performing practical analyses of realistic structures is demonstrated by simulating curvilinear crack growth in a fuselage panel that is representative of a typical narrow-body aircraft. The predicted crack trajectory and fatigue life compare well with measurements of these same quantities from a full-scale pressurized panel test.

Potyondy, David O.↗

Stimulation of isoprene emissions and electron transport rates are a key mechanism of thermal tolerance in the tropical species Vismia guianensis

Tropical forests absorb large amounts of atmospheric CO 2 through photosynthesis, but high surface temperatures suppress this absorption while promoting isoprene emissions. While mechanistic isoprene emission models predict a tight coupling to photosynthetic electron transport (ETR) as a function of temperature, direct field observations of these phenomenon are lacking in the tropics and are necessary to assess the impact of a warming climate on global isoprene emissions. In this paper, we demonstrate that in the early successional species Vismia guianensis in the central Amazon, ETR rates increased with temperature in concert with isoprene emissions, even as stomatal conductance ($g _s$) and net photosynthetic carbon fixation ($P_n$) declined. We observed the highest temperatures of continually increasing isoprene emissions yet reported (50°C). While $P_n$ showed an optimum value of 32.6 ± 0.4°C, isoprene emissions, ETR, and the oxidation state of PSII reaction centers (q L ) increased with leaf temperature with strong linear correlations for ETR (ρ = 0.98) and q L (ρ= 0.99) with leaf isoprene emissions. In contrast, other photoprotective mechanisms, such as non-photochemical quenching (NPQ), were not activated at elevated temperatures. Inhibition of isoprenoid biosynthesis repressed $P_n$ at high temperatures through a mechanism that was independent of stomatal closure. While extreme warming will decrease $g_s$ and $P_n$ in tropical species, our observations support a thermal tolerance mechanism where the maintenance of high photosynthetic capacity under extreme warming is assisted by the simultaneous stimulation of ETR and metabolic pathways that consume the direct products of ETR including photorespiration and the biosynthesis of thermoprotective isoprenoids. Our results confirm that models which link isoprene emissions to the rate of ETR hold true in tropical species and provide necessary "ground-truthing" for simulations of the large predicted increases in tropical isoprene emissions with climate warming.

54 ENVIRONMENTAL SCIENCES↗

The Functor system: a new on-the-fly take on Material Properties based on C++ functions

In the context of solving multiphysics problems, the discretization of the partial differential equations (PDE) at hand often takes the spotlight. However, for most engineering users and even application developers, the discretization of the equations has already been performed. Instead, they are tasked with implementing specific closure relations and material properties. MOOSE has long enabled this using the Materials system. This system relied on the pre-computation of all properties before they are used in the PDE or in postprocessing. In this talk we will introduce the Functor system, which was deployed in MOOSE in 2021, then present a few applications of functors in flow modeling simulations by the NEAMS program. Functors first offer great flexibility in their evaluation. Rather than storing various arrays for material properties, they are evaluated on the fly at the location and state, e.g. current or old value, requested. Unlike regular material properties, several operations such as the time derivative, the divergence and the curl can be requested from a functor. Similar to material properties, functors can be made to depend on arbitrary combinations of variables, functions, postprocessors and other properties. However, unlike material properties, any of these can be substituted for a functor material property. Thanks to this, objects no longer need to be duplicated based on the types of their parameters.

97 - MATHEMATICS AND COMPUTING↗

The accuracy of multi-group models for nonlocal electron transport in magnetized plasmas

In the extreme conditions of inertial confinement fusion experiments, heat flow plays a vital role, but local diffusive models frequently break down and overestimate the heat flow. The situation becomes more complicated again in the significant magnetic fields generated during laser–plasma interactions or in magnetized fusion schemes. Accurate non-local and magnetized heat flow computations can be carried out using Vlasov–Fokker–Planck (VFP) simulations, but these are computationally expensive. There is, therefore, significant interest in using faster multi-group models to accurately calculate the non-local heat flow in magnetized plasmas. We benchmark two such multi-group models for calculating the heat flow, M1 and hybrid-AWBS-BGK, against diffusive models and full VFP simulations, before applying the models to realistic example test cases, both magnetized and unmagnetized. We find that the multi-group models generally perform very well for moderate non-localities up to kλmfp∼0.01, but the computational cost increases dramatically. hybrid-AWBS-BGK performs more effectively than M1 at high non-localities, up to kλmfp∼1, due to its adaptive solver and robust P1 closure, but tends to fail in very strong magnetic fields. Both codes are much faster than VFP simulations but are still slow in steep temperature gradients.

Arran, C. (ORCID:0000000286448118)↗

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗

Basin Scale Estimates of Evapotranspiration Using GRACE and other Observations

Evapotranspiration is integral to studies of the Earth system, yet it is difficult to measure on regional scales. One estimation technique is a terrestrial water budget, i.e., total precipitation minus the sum of evapotranspiration and net runoff equals the change in water storage. Gravity Recovery and Climate Experiment (GRACE) satellite gravity observations are now enabling closure of this equation by providing the terrestrial water storage change. Equations are presented here for estimating evapotranspiration using observation based information, taking into account the unique nature of GRACE observations. GRACE water storage changes are first substantiated by comparing with results from a land surface model and a combined atmospheric-terrestrial water budget approach. Evapotranspiration is then estimated for 14 time periods over the Mississippi River basin and compared with output from three modeling systems. The GRACE estimates generally lay in the middle of the models and may provide skill in evaluating modeled evapotranspiration.

Rodell, M.↗

DEVELOPMENT AND APPLICATION OF BENCHMARK EXAMPLES FOR MIXED-MODE I/II FATIGUE GROWTH PREDICTIONS

The development of benchmark examples for assessing growth prediction capabilities is presented and demonstrated for a commercial code. The examples are based on finite element models of Mixed-Mode Bending (MMB) specimens for three different mode ratios. The examples are independent of the analysis software used and allow the assessment of the automated delamination propagation, onset and growth prediction capabilities in commercial finite element codes based on the virtual crack closure technique (VCCT). First, based on existing results for delamination propagation under static loading, new benchmark examples for delamination growth under fatigue loading were created. Second, starting from an initially straight front, the delamination was allowed to grow under fatigue loading. The number of cycles to delamination onset and the number of cycles during delamination growth for each growth increment were obtained from the automated analysis and compared to the benchmark examples. Good agreement between the results obtained from the growth analysis and the benchmark results could be achieved by selecting the appropriate input parameters. The benchmarking procedure proved valuable by highlighting the issues associated with choosing the input parameters of the particular implementation. Overall the results from automatic propagation are encouraging, however, selecting the appropriate input parameters was not straightforward and often required an iterative procedure.

Composites↗