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 469 records · Page 26

Field surveying data of low-cost networked flood sensors in southeast Texas

Floods are common natural disasters worldwide and pose substantial risks to life, property, food production, and natural resources. Effective measures for flood mitigation and warning are essential. Southeast Texas is still at significant risk of flooding, and Lamar University is assisting the region with asset management of a flood sensor network for flooding events. This network provides real-time water stage information. Lamar University developed a survey program to measure elevation and coordinates at each sensor site location to make this data more useful for flood monitoring and mapping. This paper overviews the measurement of the elevation and coordinates of 74 networked flood sensors and various flood stage thresholds at critical points that flood decision-makers can use for reference at each site. In the first phase of this program, these sensors were deployed throughout a 7-county region spanning nearly 6,000 square miles in Southeast Texas. The latitude and longitude of the sensors and their elevations were determined using survey-grade Global Navigation Satellite System (GNSS) technology. Various Continually Operating Reference Stations (CORS) were utilized for post-processing to achieve sub-inch resolution. The flood stage thresholds, water level sensors elevation, and the elevations and positions of other critical surrounding points are viewable to the public through two online repositories and a web-based sensor management dashboard. The data is used to aid with decisions related to road closures or modeling efforts by mitigation decision-makers, emergency managers, and the public, including the Texas Department of Transportation, Houston Transtar, the National Weather Service, and the Sabine River Authority of Texas (SRA).

54 ENVIRONMENTAL SCIENCES↗

Spectroscopic identification of Ca-bearing uranyl silicates formed in C–S–H systems

Portland cement-based grouts used for radioactive waste immobilisation contain a Ca- and Si-rich binder phase, known as calcium–silicate–hydrate (C–S–H). Depending on the blend of cement used, the Ca/Si ratio can vary considerably. A range of C–S–H minerals with Ca/Si ratios from 0.6 to 1.6 were synthesised and contacted with aqueous U(VI) at 0.5 mM and 10 mM concentrations. Solid-state 29 Si MAS-NMR spectroscopy was applied to probe the Si coordination environment in U(VI)-contacted C–S–H minerals and, in conjunction with U L III -edge X-ray absorption spectroscopy analysis, inferences of the fate of U(VI) in these systems were made. At moderate or high Ca/Si ratios, uranophane-type uranyl silicates or Ca-uranates dominated, while at the lowest Ca/Si ratios, the formation of a Ca-bearing uranyl silicate mineral, similar to haiweeite (Ca[(UO 2 ) 2 Si 5 O 12 (OH) 2 ]·3H 2 O) or Ca-bearing weeksite (Ca 2 (UO 2 ) 2 Si 6 O 15 ·10H 2 O) was identified. This study highlights the influence of Ca/Si ratio on uranyl sequestration, of interest in the development of post-closure safety models for U-bearing radioactive waste disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Growing neoclassical tearing modes seeded via transient-induced-multimode interactions

Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗

Validation Assessment of Turbulent Reacting Flow Model Using The Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

The chemistry and diffusion of aircraft exhausts in the lower stratosphere during the first few hours after fly-by

An analysis of the hydrogen-nitrogen-oxygen reaction systems in the lower stratosphere as they are initially perturbed by individual aircraft engine exhaust plumes was conducted in order to determine whether any significant chemical reactions occur, either among exhaust chemical species, or between these species and the environmental ozone, while the exhaust products are confined to intact plume segments at relatively high concentrations. The joint effects of diffusive mixing and chemical kinetics on the reactions were also studied, using the techniques of second-order closure diffusion/chemistry models. The focus of the study was on the larger problem of the potential depletion of ozone by supersonic transport aircraft exhaust materials emitted into the lower stratosphere.

Hilst, G. R.↗

Viscous effects in aircraft trailing vortices

The mechanism of merging of like-signed aircraft vortices leading to a rapid redistribution of trailed vorticity in a wake through both convective and turbulent processes was investigated. Research was done experimentally in a small wind tunnel and analytically through the use of a code which computes turbulent transport using a second-order closure turbulent model. Computations are reported which demonstrate the merging phenomenon, and comparisons are made with experimental results. The usefulness of point vortex computations in predicting merging was explored. Limited computations showed that jet exhaust does not appreciably alter the merging phenomenon. The effect of ambient atmospheric turbulence on the aging of an aircraft wake was investigated at a constant turbulent dissipation rate. It was shown that under stable atmospheric conditions, when atmospheric macroscales are less than or equal to the vortex spacing, misleading results may be obtained.

Bilanin, A. J.↗

Second order closure for variable density free shear layer

A full second order closure turbulence model for the prediction of free shear flows with variable density is developed. The importance of density fluctuations due to inhomogeneities in fluid composition is emphasized. Some simplicity is achieved by using density weighted averaging to account implicitly for density fluctuation effects. Transport equations for various second order correlations - such as Reynolds stress, turbulent mass fluxes or density velocity correlations - are derived, modelled and solved, together with the mean equations for momentum and species conservation. Solutions of these equations are compared with the two different sets of experimental data for low speed free mixing layers (helium into nitrogen and freon 12 into air).

Vandromme, D.↗

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The main objective is to extend the boundaries within which large eddy simulations (LES) and direct numerical simulations (DNS) can be applied in computational analyses of high speed reacting flows. In the efforts related to LES, we were concerned with developing reliable subgrid closures for modeling of the fluctuation correlations of scalar quantities in reacting turbulent flows. In the work on DNS, we focused our attention to further investigation of the effects of exothermicity in compressible turbulent flows. In our previous work, in the first year of this research, we have considered only 'simple' flows. Currently, we are in the process of extending our analyses for the purpose of modeling more practical flows of current interest at LaRC. A summary of our accomplishments during the third six months of the research is presented.

Givi, P.↗

Turbulence and stall in plane diffusers - Computational study

The effect of free-stream turbulence level on stall conditions in plane diffusers is investigated with a finite-volume numerical method which utilizes a k-epsilon turbulence model or a differential Reynolds stress model for closure. Computed results show qualitative agreement with measured data. The ability to eliminate stall and thereby increase pressure recovery in diffusers at large included angles by utilizing high free-stream turbulence is correctly predicted.

Demuren, A. O.↗

Simulation of Crack Propagation in Engine Rotating Components under Variable Amplitude Loading

The crack propagation life of tested specimens has been repeatedly shown to strongly depend on the loading history. Overloads and extended stress holds at temperature can either retard or accelerate the crack growth rate. Therefore, to accurately predict the crack propagation life of an actual component, it is essential to approximate the true loading history. In military rotorcraft engine applications, the loading profile (stress amplitudes, temperature, and number of excursions) can vary significantly depending on the type of mission flown. To accurately assess the durability of a fleet of engines, the crack propagation life distribution of a specific component should account for the variability in the missions performed (proportion of missions flown and sequence). In this report, analytical and experimental studies are described that calibrate/validate the crack propagation prediction capability ]or a disk alloy under variable amplitude loading. A crack closure based model was adopted to analytically predict the load interaction effects. Furthermore, a methodology has been developed to realistically simulate the actual mission mix loading on a fleet of engines over their lifetime. A sequence of missions is randomly selected and the number of repeats of each mission in the sequence is determined assuming a Poisson distributed random variable with a given mean occurrence rate. Multiple realizations of random mission histories are generated in this manner and are used to produce stress, temperature, and time points for fracture mechanics calculations. The result is a cumulative distribution of crack propagation lives for a given, life limiting, component location. This information can be used to determine a safe retirement life or inspection interval for the given location.

Bonacuse, P. J.↗

Residual Strength Prediction of Fuselage Structures with Multiple Site Damage

This paper summarizes recent results on simulating full-scale pressure tests of wide body, lap-jointed fuselage panels with multiple site damage (MSD). The crack tip opening angle (CTOA) fracture criterion and the FRANC3D/STAGS software program were used to analyze stable crack growth under conditions of general yielding. The link-up of multiple cracks and residual strength of damaged structures were predicted. Elastic-plastic finite element analysis based on the von Mises yield criterion and incremental flow theory with small strain assumption was used. A global-local modeling procedure was employed in the numerical analyses. Stress distributions from the numerical simulations are compared with strain gage measurements. Analysis results show that accurate representation of the load transfer through the rivets is crucial for the model to predict the stress distribution accurately. Predicted crack growth and residual strength are compared with test data. Observed and predicted results both indicate that the occurrence of small MSD cracks substantially reduces the residual strength. Modeling fatigue closure is essential to capture the fracture behavior during the early stable crack growth. Breakage of a tear strap can have a major influence on residual strength prediction.

Chen, Chuin-Shan↗

Simulation of Crack Propagation in Engine Rotating Components Under Variable Amplitude Loading

The crack propagation life of tested specimens has been repeatedly shown to strongly depend on the loading history. Overloads and extended stress holds at temperature can either retard or accelerate the crack growth rate. Therefore, to accurately predict the crack propagation life of an actual component, it is essential to approximate the true loading history. In military rotorcraft engine applications, the loading profile (stress amplitudes, temperature, and number of excursions) can vary significantly depending on the type of mission flown. To accurately assess the durability of a fleet of engines, the crack propagation life distribution of a specific component should account for the variability in the missions performed (proportion of missions flown and sequence). In this report, analytical and experimental studies are described that calibrate/validate the crack propagation prediction capability for a disk alloy under variable amplitude loading. A crack closure based model was adopted to analytically predict the load interaction effects. Furthermore, a methodology has been developed to realistically simulate the actual mission mix loading on a fleet of engines over their lifetime. A sequence of missions is randomly selected and the number of repeats of each mission in the sequence is determined assuming a Poisson distributed random variable with a given mean occurrence rate. Multiple realizations of random mission histories are generated in this manner and are used to produce stress, temperature, and time points for fracture mechanics calculations. The result is a cumulative distribution of crack propagation lives for a given, life limiting, component location. This information can be used to determine a safe retirement life or inspection interval for the given location.

P J Bonacuse↗

Further Studies on Oceanic Biogeochemistry and Carbon Cycling

This TM consists of two chapters. Chapter I describes the development of a coupled, one-dimensional biogeochemical model using turbulence closure mixed layer (TCMLM) dynamics. The model is applied to the Sargasso Sea at the BATS (Bermuda Atlantic Time Series) site and the results are compared with a previous model study in the same region described in NASNTP-2001-209991. The use of the TCMLM contributed to some improvements in the model simulation of chlorophyll, PAR, nitrate, phosphate, and oxygen, but most importantly, the current model achieved good agreement with the data with much more realistic background eddy diffusivity. However, off-line calculations of horizontal transport of biogeochemical properties revealed that one-dimensional dynamics can only provide a limited assessment of the nutrient and carbon balances at BATS. Future studies in the BATS region will require comprehensive three-dimensional field studies, combined with three-dimensional eddy resolving numerical experiments, to adequately quantify the impact of the local and remote forcing on ecosystem dynamics and carbon cycling. Chapter II addresses the sensitivity of global sea-air CO, flux estimates to wind speed, temperature, and salinity. Sensitivity analyses of sea-air CO, flux to wind speed climatologies, gas transfer algorithms, SSS and SST were conducted for the global oceans and regional domains. Large uncertainties in the global sea-air flux are identified, primarily due to the different gas transfer algorithms used. The sensitivity of the sea-air flux to SST and SSS is similar in magnitude to the effect of using different wind climatologies. Globally, the mean ocean uptake of CO, changes by 5 to 16%, depending upon the combination of SST and SSS used.

Signorini, S. R.↗

Formation Ages of the Apollo 16 Regolith Breccias: Implications for Accessing the Bombardment History of the Moon

Regolith breccias are lithified samples of the regolith that have been fused together by impact shock and thermal metamorphism. In lunar regolith samples, the ratio of trapped 40Ar/36Ar is a useful indicator of antiquity and can be used to model the closure age/lifithication event of the regolith (i.e. the apparent time when Ar became trapped [1]), thus providing an important insight into specific times when that regolith was interacting with the the dynamic inner solar system space environment [2-4].

Joy, K. H.↗