Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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 217 records · Page 12

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗

Qualification of ANSI/HPS N13.1-2011 Mixing Criteria by Computational Fluid Dynamics Modeling for the 3430 Building Fan Addition and Increased Ventilation Capacity

Additional ventilation capacity has been designed for the 3430 Building filtered exhaust stack system. The updated system will increase the number of fans from two to three and include new ductwork with a larger diameter to integrate the new fan into the existing stack. Stack operations will involve running various fan combinations at any given time. The air monitoring system of the existing two-fan stack previously was found to comply with the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard. Full-scale, three-dimensional computational fluid dynamics (CFD) modeling was used to evaluate the modified three-fan system for compliance with the ANSI/HPS N13.1-2011 standard, which essentially is equivalent to the ANSI/HPS N13.1-1999 standard. The four mixing criteria evaluated are 1) flow angle, 2) velocity, 3) gas tracer, and 4) particle tracer. Benchmarking of the CFD modeling methodology showed good agreement with previous testing used to qualify the stack, and modeling of the existing two-fan system showed good agreement with test data collected from the 3430 Building stack. Modeling was performed to develop a suitable three-fan design. Initial modeling of the three-fan design and basic ductwork showed that flow angles and velocity uniformity were acceptable; however, the gas tracer and particle tracer mixing results were not acceptable. To meet ANSI/HPS N13.1-2011 criteria, an air blender was added to the stack design.

42 ENGINEERING↗

Qualification of ANSI/HPS N13.1-2011 Mixing Criteria by Computational Fluid Dynamics Modeling for the 3430 Building Fan Addition and Increased Ventilation Capacity

Additional ventilation capacity has been designed for the 3430 Building filtered exhaust stack system. The updated system will increase the number of fans from two to three and include new ductwork with a larger diameter to integrate the new fan into the existing stack. Stack operations will involve running various fan combinations at any given time. The air monitoring system of the existing two-fan stack previously was found to comply with the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard. Full-scale, three-dimensional computational fluid dynamics (CFD) modeling was used to evaluate the modified three-fan system for compliance with the ANSI/HPS N13.1-2011 standard, which essentially is equivalent to the ANSI/HPS N13.1-1999 standard. The four mixing criteria evaluated are 1) flow angle, 2) velocity, 3) gas tracer, and 4) particle tracer. Benchmarking of the CFD modeling methodology showed good agreement with previous testing used to qualify the stack, and modeling of the existing two-fan system showed good agreement with test data collected from the 3430 Building stack. Modeling was performed to develop a suitable three-fan design. Initial modeling of the three-fan design and basic ductwork showed that flow angles and velocity uniformity were acceptable; however, the gas tracer and particle tracer mixing results were not acceptable. To meet ANSI/HPS N13.1-2011 criteria, an air blender was added to the stack design. This revision models the individual maximum fan reduced flow capacity from 38,000 cfm to 31,200 cfm; no changes to the duct design are made.

42 ENGINEERING↗

An AI-based approach to structural damage identification by modal analysis

Flexible-structure damage is presently addressed by a combined model- and parameter-identification approach which employs the AI methodologies of classification, heuristic search, and object-oriented model knowledge representation. The conditions for model-space search convergence to the best model are discussed in terms of search-tree organization and initial model parameter error. In the illustrative example of a truss structure presented, the use of both model and parameter identification is shown to lead to smaller parameter corrections than would be required by parameter identification alone.

Glass, B. J.↗

Face Sheet/Core Disbond Growth in Honeycomb Sandwich Panels Subjected to Ground-Air-Ground Pressurization and In-Plane Loading

Typical damage modes in light honeycomb sandwich structures include face sheet/core disbonding and core fracture, both of which can pose a threat to the structural integrity of a component. These damage modes are of particular interest to aviation certification authorities since several in-service occurrences, such as rudder structural failure and other control surface malfunctions, have been attributed to face sheet/core disbonding. Extensive studies have shown that face sheet/core disbonding and core fracture can lead to damage propagation caused by internal pressure changes in the core. The increasing use of composite sandwich construction in aircraft applications makes it vitally important to understand the effect of ground-air-ground (GAG) cycles and conditions such as maneuver and gust loads on face sheet/core disbonding. The objective of the present study was to use a fracture mechanics based approach developed earlier to evaluate the loading at the disbond front caused by ground-air-ground pressurization and in-plane loading. A honeycomb sandwich panel containing a circular disbond at one face sheet/core interface was modeled with three-dimensional (3D) solid finite elements. The disbond was modeled as a discrete discontinuity and the strain energy release rate along the disbond front was computed using the Virtual Crack Closure Technique (VCCT). Special attention was paid to the pressure-deformation coupling which can decrease the pressure load within the disbonded sandwich section significantly when the structure is highly deformed. The commercial finite element analysis software, Abaqus/Standard, was used for the analyses. The recursive pressure-deformation coupling problem was solved by representing the entrapped air in the honeycomb cells as filled cavities in Abaqus/Standard. The results show that disbond size, face sheet thickness and core thickness are important parameters that determine crack tip loading at the disbond front. Further, the pressure-deformation coupling was found to have an important load decreasing effect [6]. In this paper, a detailed problem description is provided first. Second, the analysis methodology is presented. The fracture mechanics approach used is described and the specifics of the finite element model, including the fluid-filled cavities, are introduced. Third, the initial model verification and validation are discussed. Fourth, the findings from a closely related earlier study [6] are summarized. These findings provided the basis for the current investigation. Fifth, an aircraft ascent scenario from 0 to 12192 m (0 to 40000 ft) is considered and the resulting crack tip loading at the disbond front is determined. In-plane loading to simulate maneuvers and gust conditions are also considered. Sixth, the results are shown for a curved panel, which was used to simulate potential fuselage applications. Finally, a brief summary of observations is presented and recommendations for improvement are provided.

Chen, Zhi M.↗

Magnetic Shield Design Modeling and Validation for SWOT Spacecraft Ka-Band Extended Interaction Klystron

Two Extended Interaction Klystrons (EIKs) containing strong permanent magnets were modeled magnetically in a representative spacecraft geometry using commercial finite element modeling techniques and were validated against measurements made at varying distances. Initial modeling results for the 63 A-m2 (amperes per square meter) dipole moment magnets showed that magnetic shields would be necessary in order to meet magnetic field requirements for the Surface Water and Ocean Topography (SWOT) spacecraft, which contains components that are susceptible to external DC (Direct Current) magnetic fields. JPL (Jet Propulsion Laboratory) and the EIK vendor proposed cold rolled steel and mu-metal as potential shield materials along with proposed thicknesses of 0.5 millimeters and 1.5 millimeters. Magnetic shields made from each of these materials were designed and modeled in software, taking high-field saturation into account. Prototype magnetic shields with these parameters were then built, measured with an existing EIK, and compared against modeling results. For single-axis field measurements along the dipole axis, modeling results were within 7 gauss of the measured values at 10 centimeters from the magnet, and converged to less than 1.5 gauss at distances greater than 14 centimeters from the magnet. Three-axis field measurements at locations of interest showed that model correlation improved to within 4 gauss at 11 centimeters and 2 gauss for distances ranging between 15 centimeters and 36 centimeters.

Gonzales, Edward↗

Magnetic Shield Design Modeling and Validation for SWOT Spacecraft Ka-Band Extended Interaction Klystron

Two Extended Interaction Klystrons (EIKs) containing strong permanent magnets were modeled magnetically in a representative spacecraft geometry using commercial finite element modeling techniques and were validated against measurements made at varying distances. Initial modeling results for the 63 A-m2 dipole moment magnets showed that magnetic shields would be necessary in order to meet magnetic field requirements for the Surface Water and Ocean Topography (SWOT) spacecraft, which contains components that are susceptible to external DC magnetic fields. JPL and the EIK vendor proposed cold rolled steel and mu-metal as potential shield materials along with proposed thicknesses of 0.5 mm and 1.5 mm. Magnetic shields made from each of these materials were designed and modeled in software, taking highfield saturation into account. Prototype magnetic shields with these parameters were then built, measured with an existing EIK, and compared against modeling results. For single-axis field measurements along the dipole axis, modeling results were within 7 gauss of the measured values at 10 cm from the magnet, and converged to less than 1.5 gauss at distances greater than 14 cm from the magnet. Three-axis field measurements at locations of interest showed that model correlation improved to within 4 gauss at 11 cm and 2 gauss for distances ranging between 15 cm and 36 cm.

Gonzales, Edward↗

Summary of 2020 edge-plasma modeling for FNSFand plans for 2021

This year, following the study of liquid-wall options, our edge-plasma modeling has returned to analyzing issues associated with a solid tungsten-wall, Fusion Nuclear Science Facility (FNSF)-size device, primarily using the 2D UEDGE plasma/neutral transport code. The focus is to predict plasma/neutral heat and particle fluxes to various wall components and determine the concentration of injected and wall-sputtered impurities that reach the core plasma region. Our models have been advanced in a number of areas compared to the initial modeling done in the 2017 timeframe: this year we include each individual charge-state of impurity ions instead of assuming a fixed concentration of impurities relative to the deuterium-tritium (DT) fuel plasma, including the interaction between inner and outer divertor legs. Here edge fueling of the impurity and pumping in the private-flux region is incorporated. Increasingly detailed atomic-physics rate files for impurity radiation have been compared, as well as the role of molecules. The impact of spatial variation of the turbulence-driven radial plasma transport is considered to contribute to the understanding of the heat-flux-limit effects of choosing a double-null or single-null divertor configuration. Detailed plasma/neutral wall fluxes, densities, and temperatures have been provided to the ERO group for their more detailed studies of wall sputtering and near-surface impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-resolution heat-transfer-coefficient maps applicable to compound-curve surfaces using liquid crystals in a transient wind tunnel

Tests were performed in a transient heat transfer tunnel in which the model under test was preheated prior to allowing room temperature air to be suddenly drawn over the model. The resulting movement of isothermal contours on the model is revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record is obtained of a temperature and time data pair for all points on the model during a single test. Experiments on a duct model are reported in which the model was preheated using a hot air stream. A manner in which initial model temperature nonuniformities could be taken into account was investigated. The duct model was also tested with a steady-state measurement technique and results were compared with the transient measurements, but recognizing that differences existed between the upstream thermal boundary conditions. The steady-state and transient measurements were shown to be consistent with predicted values. The main advantage of this transient heat transfer technique using liquid crystals is that since the test model need not be actively heated, high-resolution measurements on surfaces with complex shapes may be obtained.

Jones, T. V.↗

High-resolution heat-transfer-coefficient maps applicable to compound-curve surfaces using liquid crystals in a transient wind tunnel

Tests were performed in a transient heat transfer tunnel in which the model under test was preheated prior to allowing room temperature air to be suddenly drawn over the model. The resulting movement of isothermal contours on the model is revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record is obtained of a temperature and time data pair for all points on the model during a single test. Experiments on a duct model are reported in which the model was preheated using a hot air stream. A manner in which initial model temperature nonuniformities could be taken into account was investigated. The duct model was also tested with a steady-state measurement technique and results were compared with the transient measurements, but recognizing that differences existed between the upstream thermal boundary conditions. The steady-state and transient measurements were shown to be consistent with predicted values. The main advantage of this transient heat transfer technique using liquid crystals is that since the test model need not be actively heated, high-resolution measurements on surfaces with complex shapes may be obtained.

Jones, Terry V.↗

Monsoonal MCS Initiation, Rainfall, and Diurnal Gravity Waves over the Bay of Bengal: Observation and a Linear Model

Abstract Previous observational studies have indicated that mesoscale convective systems (MCSs) contribute the majority of precipitation over the Bay of Bengal (BoB) during the summer monsoon season, yet their initiation and propagation remain incompletely understood. To fill this knowledge gap, we conducted a comprehensive study using a combination of 20-yr satellite observations, MCS tracking, reanalysis data, and a theoretical linear model. Satellite observations reveal clear diurnal propagation signals of MCS initiation frequency and rainfall from the west coast of the BoB toward the central BoB, with the MCS rainfall propagating slightly slower than the MCS initiation frequency. Global reanalysis data indicate a strong association between the offshore-propagating MCS initiation frequency/rainfall and diurnal low-level wind perturbations, implying the potential role of gravity waves. To verify the hypothesis, we developed a 2D linear model that can be driven by realistic meteorological fields from reanalysis. The linear model realistically reproduces the characteristics of offshore-propagating diurnal wind perturbations. The wind perturbations, as well as the offshore propagation signals of MCS initiation frequency and rainfall, are associated with diurnal gravity waves emitted from the coastal regions, which in turn are caused by the diurnal land–sea thermal contrast. The ambient wind speed and vertical wind shear play crucial roles in modulating the timing, propagation, and amplitude of diurnal gravity waves. Using the linear model and satellite observations, we further show that the stronger monsoonal flows lead to faster offshore propagation of diurnal gravity waves, which subsequently control the offshore propagation signals of MCS initiation and rainfall. Significance Statement Rainfall over the Bay of Bengal (BoB) is primarily contributed by large and organized rainfall systems in the summer monsoon season. During this season, these systems are commonly observed over the east coast of India around midnight, the western BoB in the morning, and the central BoB in the afternoon. This eastward rainfall propagation is confirmed by observations, reanalysis data, and a theoretical model to have a strong association with atmospheric diurnal gravity waves. These waves are caused by land–sea thermal contrast and can trigger rainfall systems over the offshore regions. We also found that the diurnal gravity waves, as well as their triggered rainfall systems, can be greatly modulated by the large-scale monsoonal flows. All the above findings improve our understanding of diurnal rainfall cycle over the BoB and may contribute to the future improvement of rainfall forecast over the region.

Meteorology & Atmospheric Sciences↗

Implications of rootless geothermal models: Missing processes, parameter compensation, and imposter convection

Numerical models of geothermal systems commonly capture only the top of a reservoir. Deeper areas of the reservoir are simplified to a boundary at the base of the model domain. Commonly, the basal boundary is given either a heat source or a source of mass and enthalpy. Here we developed and present simple numerical experiments which demonstrate that these approaches do not produce the correct model behavior in comparison to a model that captures the entire convecting domain with a heat flux only. We describe a variety of incorrect types of model behavior that arise directly from the choice of boundary condition, independent of the specific parameterization of the model. Heat sources are sensitive to the thickness of the domain and parameters take unphysical values to compensate for the reduced height. The combined mass/heat boundary can produce temperatures that show similarities to circulating geothermal systems, but with incorrect fluid flow and a strong boundary layer focused at the base of the simulated clay cap. These errors likely cause parameters to adjust their values to compensate for the incorrect physics. We highlight these issues and show an example from a developed reservoir model. Initial calibrations to natural state temperature were unsuitable for history matching. The 3D model required parameter adjustments to achieve a more realistic production model. Appropriate mitigation measures should be considered to reduce parameter compensation and improve decision-support models.

15 GEOTHERMAL ENERGY↗

Data for Robust Paths to Net Greenhouse Gas Mitigation and Negative Emissions via Advanced Biofuels

This zip file contains a UNIX-format DayCent model executable, input files, automation code, and associated directory structure necessary to re-produce the DayCent analysis underlying the manuscript. The main script “autodaycent.py” (written for Python 2.7) opens an interactive command line routine that facilitates: Calibrating the DayCent pine growth model; Initializing DayCent for a set of case studies sites; Executing an ensemble of model runs representing case study site reforestation, grassland restoration, or conversion to switchgrass cultivation; and Results analysis & generation of manuscript Fig. 3. Note that the interactive analysis code requires that all input files to be contained in the directory structure as uploaded, without modification. Executable versions of the DayCent model compatible with other operating systems are available upon request.

Feedstock Production↗

Ozone Sensitivity to Varying Greenhouse Gases and Ozone-Depleting Substances in CCMI-1 Simulations

Ozone fields simulated for the first phase of the Chemistry-Climate Model Initiative (CCMI-1) will be used as forcing data in the 6th Coupled Model Intercomparison Project. Here we assess, using reference and sensitivity simulations produced for CCMI-1, the suitability of CCMI-1 model results for this process, investigating the degree of consistency amongst models regarding their responses to variations in individual forcings. We consider the influences of methane, nitrous oxide, a combination of chlorinated or brominated ozone-depleting substances, and a combination of carbon dioxide and other greenhouse gases. We find varying degrees of consistency in the models' responses in ozone to these individual forcings, including some considerable disagreement. In particular, the response of total-column ozone to these forcings is less consistent across the multi-model ensemble than profile comparisons. We analyse how stratospheric age of air, a commonly used diagnostic of stratospheric transport, responds to the forcings. For this diagnostic we find some salient differences in model behaviour, which may explain some of the findings for ozone. The findings imply that the ozone fields derived from CCMI-1 are subject to considerable uncertainties regarding the impacts of these anthropogenic forcings. We offer some thoughts on how to best approach the problem of generating a consensus ozone database from a multi-model ensemble such as CCMI-1.

CCMI↗

GMI-IPS: Processing & Visualization Software Used in ATom DC-8 Aircraft Studies

NASA's Atmospheric Tomography Mission (ATom) deployed in each of the four seasons during 2016-2018, the DC-8 aircraft in order to establish global-scale datasets intended to improve the representation of chemically reactive gases in global atmospheric chemistry models (ACMs). The Global Modeling Initiative (GMI) executed simulations for each ATom flight using the GMI Chemistry Transport Model (GMI-CTM) to provide species concentrations of chemical gases along the DC-8 flight transects. To solve the problem of translating the GMI-CTM simulation data to the unique spatial resolutions of each ATom flight, the GMI ICARTT Processing Software (GMI-IPS) was developed.The GMI-IPS is written in Python and provides data processing, flight extraction, and visualization support for aircraft research projects using ICARTT format, which is a standard format for airborne instrument data. Additionally, the GMI-IPS interpolates global gridded model data from Hierarchical Data Format (HDF) to ICARTT compatible flight transects. Software classes for instruments and collections provided by the ATom DC-8 aircraft such as MER10, MMS, etc. are derived from a common base class. Other functionality provided by the GMI-IPS are: deriving missing flight entries along a transect, reading ICARTT entries from file, and providing Python data structures for storing flight and model information, and more.The GMI-IPS is GIT source controlled, has approximately 30,000 lines of code, and supports parallelization across data collections. It delivered GMI-CTM data for more than forty distinct DC-8 aircraft flights that took place under ATom. The output ICARTT files adhere to format standard V1.1, and pass the scan utility provided by NASA LaRC Airborne Science Data for Atmospheric Composition. This presentation will include a software and methods overview, and results from ATom, including assessments using the GMI-CTM showing how well observations from ATom flight transects represent a broader region.

Damon, M. R.↗

Comparative Analysis of Static and Dynamic Probabilistic Risk Assessment

This study examines three different methodologies for producing loss-of-mission (LOM) and loss-of-crew (LOC) risks estimates for probabilistic risk assessments (PRA) of crewed spacecraft. The three bottom-up, component-based PRA approaches examined are a traditional static fault tree, a dynamic Monte Carlo simulation, and a fault tree hybrid that incorporates some dynamic elements. These approaches were used to model the reaction control system thruster pod of a generic crewed spacecraft and mission, and a comparative analysis of the methods is presented. The methodologies are assessed in terms of the process of modeling a system, the actionable information produced for the design team, and the overall fidelity of the quantitative risk evaluation generated. The system modeling process is compared in terms of the effort required to generate the initial model, update the model in response to design changes, and support mass-versus-risk trade studies. The results are compared by examining the top-level LOM/LOC estimates and the relative risk driver rankings at the failure mode level. The fidelity of each modeling methodology is discussed in terms of its capability to handle real-world system dynamics such as cold-sparing, changes in mission operations due to loss of redundancy, and common cause failure modes. The paper also discusses the applicability of each methodology to different phases of system development and shows that a single methodology may not be suitable for all of the many purposes of a spacecraft PRA. The fault tree hybrid approach is shown to be best suited to the needs of early assessments during conceptual design phases. As the design begins to mature, the level of detail represented in the risk model must go beyond redundancy and nominal mission operations to include dynamic, time- and state-dependent system responses as well as diverse system capabilities. This is best accomplished using the dynamic simulation approach, since these phenomena are not easily captured by static methods. Ultimately, once the design has been finalized and the goal of the PRA is to provide design validation and requirement verification, more traditional, static fault tree approaches may become as appropriate as the simulation method.

Mattenberger, Christopher J.↗

Improving Air Quality Forecasts with AURA Observations

Past studies have identified model initial and boundary conditions as sources of reducible errors in air-quality simulations. In particular, improving the initial condition improves the accuracy of short-term forecasts as it allows for the impact of local emissions to be realized by the model and improving boundary conditions improves long range transport through the model domain, especially in recirculating anticyclones. During the August 2006 period, we use AURA/OMI ozone measurements along with MODIS and CALIPSO aerosol observations to improve the initial and boundary conditions of ozone and Particulate Matter. Assessment of the model by comparison of the control run and satellite assimilation run to the IONS06 network of ozonesonde observations, which comprise the densest ozone sounding campaign ever conducted in North America, to AURA/TES ozone profile measurements, and to the EPA ground network of ozone and PM measurements will show significant improvement in the CMAQ calculations that use AURA initial and boundary conditions. Further analyses of lightning occurrences from ground and satellite observations and AURA/OMI NO2 column abundances will identify the lightning NOx signal evident in OMI measurements and suggest pathways for incorporating the lightning and NO2 data into the CMAQ simulations.

Newchurch, M. J.↗

Visualization of Hemispherical Post‐Detonation Fireball Internal Structures

Hemispherical charges are initiated above a transparent plate in a half‐plane configuration, allowing optical access to the internal regions of the luminous fireball. High‐speed visualization enables characterization of the internal luminous structure, showing clearly the bright shell and dark core regions. Spectroscopic and pyrometric diagnostics are applied to these flows, generating quantitative information on the spatial distribution of temperature inside the fireball. Both cased and uncased charges can be examined using this methodology, and several different high explosives are tested. This approach can be useful in validation of detailed explosive fireball models. Initial comparisons with such models are presented.

detonation↗