Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

A TRMM Rainfall Estimation Method Applicable to Land Areas

Methods developed to estimate rain rate on a footprint scale over land with the satellite-borne multispectral dual-polarization Special Sensor Microwave Imager (SSM/1) radiometer have met with limited success. Variability of surface emissivity on land and beam filling are commonly cited as the weaknesses of these methods. On the contrary, we contend a more significant reason for this lack of success is that the information content of spectral and polarization measurements of the SSM/I is limited. because of significant redundancy. As a result, the complex nature and vertical distribution C, of frozen and melting ice particles of different densities, sizes, and shapes cannot resolved satisfactorily. Extinction in the microwave region due to these complex particles can mask the extinction due to rain drops. Because of these reasons, theoretical models that attempt to retrieve rain rate do not succeed on a footprint scale. To illustrate the weakness of these models, as an example we can consider the brightness temperature measurement made by the radiometer in the 85 GHz channel (T85). Models indicate that T85 should be inversely related to the rain rate, because of scattering. However, rain rate derived from 15-minute rain gauges on land indicate that this is not true in a majority of footprints. This is also supported by the ship-borne radar observations of rain in the Tropical Oceans and Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA-COARE) region over the ocean. Based on these observations. we infer that theoretical models that attempt to retrieve rain rate do not succeed on a footprint scale. We do not follow the above path of rain retrieval on a footprint scale. Instead, we depend on the limited ability of the microwave radiometer to detect the presence of rain. This capability is useful to determine the rain area in a mesoscale region. We find in a given rain event that this rain area is closely related to the mesoscale-average rain rate. Based on this observation, in this study we have developed a method to estimate the mesoscale-average rain rate over land utilizing microwave radiometer data. Because of the high degree of geographic and seasonal variability in the nature and intensity of rain, this method requires some tuning with 15-minute rain gauge data on land. After tuning the method, it can be applied to an independent set of rain events that are close in time and space. We find that the mesoscale rain rates retrieved over the period of a month on land with this method shows a correlation of about 0.85 with respect to the surface rain-gauge observations. This mesoscale-average rain rate estimation method can be useful to extend the spatial and temporal coverage of the rainfall data provided by the Precipitation Radar on board the Tropical Rainfall Measuring Mission (TRMM) satellite.

Prabhakara, C.↗
DOE OSTI2021

Scale-dependent spatial variabilities of hydrological exchange flows and transit time in a large regulated river

Hydrological exchange flows (HEF) across the river-aquifer interface and the associated residence time of river water in the aquifer have important implications for contaminant plume migration and biogeochemical processes in the river corridor. HEFs and residence time are influenced by both subsurface physical features and hydrologic forcing related to the transport process, which can exhibit complex spatial and temporal variations. In this study, we used a massively parallel subsurface flow model and a particle-tracking model to study the influences of different control factors on spatial variability of HEFs and residence time distributions (RTD) in the Hanford Reach of the Columbia River in Washington State. A total number of 100M particles were randomly injected in time and space and then tracked in a model domain that covers a 51-km 2 area (15.1M model cells). We used hourly river stages and groundwater levels to drive the model to provide dynamic velocity fields for the particle tracking in the simulation period that was longer than 2 years. The groundwater flow simulation and particle-tracking results provide the first comprehensive assessment of the spatial distribution of HEFs and residence time in large complex river corridors. Overall, our results show that the aquifer hydrogeological structure has the strongest correlation with the extent and magnitude of exchange flux. The residence time exhibits complex patterns that are impacted by all the river geomorphologic, hydrodynamic, and hydrogeologic factors and are strongly correlated with the downwelling ratio of exchange flux. The new insights gained through this study can be used to support the development of reduced-order models of HEFs and RTDs for large complex river systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Local-scale heterogeneity of soil thermal dynamics and controlling factors in a discontinuous permafrost region

In permafrost regions, the strong spatial and temporal variability in soil temperature cannot be explained by the weather forcing only. Understanding the local heterogeneity of soil thermal dynamics and their controls is essential to understand how permafrost systems respond to climate change and to develop process-based models or remote sensing products for predicting soil temperature. In this study, we analyzed soil temperature dynamics and their controls in a discontinuous permafrost region on the Seward Peninsula, Alaska. We acquired one-year temperature time series at multiple depths (at 5 or 10 cm intervals up to 85 cm depth) at 45 discrete locations across a 2.3 km 2 watershed. We observed a larger spatial variability in winter temperatures than that in summer temperatures at all depths, with the former controlling most of the spatial variability in mean annual temperatures. We also observed a strong correlation between mean annual ground temperature at a depth of 85 cm and mean annual or winter season ground surface temperature across the 45 locations. We demonstrate that soils classified as cold, intermediate, or warm using hierarchical clustering of full-year temperature data closely match their co-located vegetation (graminoid tundra, dwarf shrub tundra, and tall shrub tundra, respectively). We show that the spatial heterogeneity in soil temperature is primarily driven by spatial heterogeneity in snow cover, which induces variable winter insulation and soil thermal diffusivity. These effects further extend to the subsequent summer by causing variable latent heat exchanges. Finally, we discuss the challenges of predicting soil temperatures from snow depth and vegetation height alone by considering the complexity observed in the field data and reproduced in a model sensitivity analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Implementation of extrinsic cohesive zone model (ECZM) in 2D finite-discrete element method (FDEM) using node binding scheme

The combined finite-discrete element method (FDEM) has been widely used for rock fracturing simulations. Conventionally, FDEM is realized using the intrinsic cohesive zone model (ICZM); however, it has the drawback of artificial compliance and high computational expense. As a complement, the extrinsic cohesive zone model (ECZM) is seen to be realized in FDEM recently, whereas the node splitting scheme utilized is cumbersome. Here, within the framework of ICZM-based FDEM, we propose a node binding scheme to efficiently bind the pre-discretized finite elements and thus guarantee the continuum behavior of materials in the elastic stage. The yield surfaces, controlled by ECZM, are dynamically embedded by invoking the pre-inserted cohesive elements. The effectiveness and efficiency of the proposed approach are validated and tested by performing a suite of numerical experiments. Compared with ICZM-based FDEM, the proposed approach can correctly capture material deformation and reduce the computation cost. In contrast to the existing ECZM-based FDEM, the proposed approach can overcome the frequent and complex element topology updating. Finally, this work provides a novel perspective that fully inherits the advantages of both ICZM and ECZM, but circumvents their shortcomings, which guarantees a more efficient and effective simulation of brittle material evolution from continuum to discontinuum.

58 GEOSCIENCES↗
DOE OSTI2020

Erratum to: Soft- and reactive landing of ions onto surfaces: concepts and applications

Soft- and reactive landing of mass-selected ions is gaining attention as a promising approach for the precisely-controlled preparation of materials on surfaces that are not amenable to deposition using conventional methods. A broad range of ionization sources and mass filters are available that make ion soft-landing a versatile tool for surface modification using beams of hyperthermal (<100?eV) ions. The ability to select the mass-to-charge ratio of the ion, its kinetic energy and charge state, along with precise control of the size, shape, and position of the ion beam on the deposition target distinguishes ion soft landing from other surface modification techniques. Soft- and reactive landing have been used to prepare interfaces for practical applications as well as precisely-defined model surfaces for fundamental investigations in chemistry, physics, and materials science. For instance, soft- and reactive landing have been applied to study the surface chemistry of ions isolated in the gas-phase, prepare arrays of proteins for high-throughput biological screening, produce novel carbon-based and polymer materials, enrich the secondary structure of peptides and the chirality of organic molecules, immobilize electrochemically-active proteins and organometallics on electrodes, create thin films of complex molecules, and immobilize catalytically active organometallics as well as ligated metal clusters. In addition, soft landing has enabled investigation of the size-dependent behavior of bare metal clusters in the critical subnanometer size regime where chemical and physical properties do not scale predictably with size. The morphology, aggregation, and immobilization of larger bare metal nanoparticles, which are directly relevant to the design of catalysts as well as improved memory and electronic devices, have also been studied using ion soft landing. This review article begins in section 1 with a brief introduction to the existing applications of ion soft- and reactive landing. Section 2 provides an overview of the ionization sources and mass filters that have been used to date for soft landing of mass-selected ions. A discussion of the competing processes that occur during ion deposition as well as the types of ions and surfaces that have been investigated follows in section 3. Section 4 discusses the physical phenomena that occur during and after ion soft landing, including retention and reduction of ionic charge along with factors that impact the efficiency of ion deposition. The influence of soft landing on the secondary structure and biological activity of complex ions is addressed in section 5. Lastly, an overview of the structure and mobility as well as the catalytic, optical, magnetic, and redox properties of bare ionic clusters and nanoparticles deposited onto surfaces is presented in section 6. © 2015 Wiley Periodicals, Inc. Mass Spec Rev 35:439–479, 2016.

Johnson, Grant E.↗
DOE OSTI2024

Quasilinear theory and modelling of gyrokinetic turbulent transport in tokamaks

The theory, development, and validation of reduced quasilinear models of gyrokinetic turbulent transport in the closed flux surface core of tokamaks is reviewed. In combination with neoclassical collisional transport, these models are successful in accurately predicting core tokamak plasma temperature, density, rotation, and impurity profiles in a variety of confinement regimes. Refined experimental tests have been performed to validate the predictions of the quasilinear models, probing changes in the dominant gyrokinetic instabilities, as reflected in fluctuation measurements, cross-phases, and transport properties. These tests continue to produce a deeper understanding of the complex mix of instabilities at both electron and ion gyroradius scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Signal Decomposition for Intrusion Detection in Reliability Assessment in Cyber Resilience (Summary Report)

The complexity of assuring cyber resilience for physical process interactions in connected systems such as energy grids increases dramatically as the coupling between processes becomes more direct and responsive. An example of this growing complexity is provided by Integrated Energy Systems (IES), in which various processes such as nuclear heat generation and commodity production are being directly coupled for increased responsiveness to highly variable signals such as market pricing or electricity demand. As such, the potential attack surface of the coupled processes is larger than the two processes independently. Securing these complex systems requires two-fold monitoring: cybersecure monitoring for potential malicious incursion, and physics monitoring for system tampering. Physics monitoring includes analyzing the behavior of the signals within the system for anomalous behavior. This analysis has been shown to be insufficient if approached by only data-driven machine learning and artificial intelligence (MLAI) techniques or only low-level model comparison. Previous efforts at Purdue University suggested combining high-fidelity models with MLAI algorithms as a basis for a software tool for detecting anomalies in physical processes. This work built on that suggestion, developing an advanced library for signal decomposition and analysis using both MLAI and high-fidelity physics algorithms for greatly improved anomaly detection, especially false data injection. This software can be used as part of a secure imbedded intelligence (SEI) system designed under Consequence-driven Cyber-informed Engineering (CCE) for complex coupled systems. This library established a foundation for online and posteriori analysis of digital signals for the purpose of detecting potential malicious tampering in digital signals representing physical processes. Demonstrations carried out throughout the development highlight the effective use of characterization algorithms to detect signal perturbations, particularly triangle attack-style perturbations, in three wide-ranging applications: seismic monitoring, nuclear thermal hydraulics system simulation, and custom manufacturing.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Spectral identification/elimination of molecular species in spacecraft glow

Computer models of molecular electronic and vibrational emission intensities were developed. Known radiative emission rates (Einstein coefficients) permit the determination of relative excited state densities from spectral intensities. These codes were applied to the published spectra of glow above shuttle surface and to the Spacelab 1 results of Torr and Torr. The theoretical high-resolution spectra were convolved with the appropriate instrumental slit functions to allow accurate comparison with data. The published spacelab spectrum is complex but N2+ Meinel emission can be clearly identified in the ram spectrum. M2 First Positive emission does not correlate well with observed features, nor does the CN Red System. Spectral overlay comparisons are presented. The spectrum of glow above shuttle surfaces, in contrast to the ISO data, is not highly structured. Diatomic molecular emission was matched to the observed spectral shape. Source excitation mechanisms such as (oxygen atom)-(surface species) reaction product chemiluminescence, surface recombination, or resonance fluorescent re-emission will be discussed for each tentative assignment. These assignments are the necessary first analytical step toward mechanism identification. Different glow mechanisms will occur above surfaces under different orbital conditions.

Green, B. D.↗
NASA NTRS2015

Smoke Dispersion Modeling Over Complex Terrain Using High-Resolution Meteorological Data and Satellite Observations: The FireHub Platform

A total number of 20,212 fire hot spots were recorded by the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument over Greece during the period 2002e2013. The Fire Radiative Power (FRP) of these events ranged from 10 up to 6000 MW at 1 km resolution, and many of these fire episodes resulted in long-range transport of smoke over distances up to several hundred kilometers. Three different smoke episodes over Greece are analyzed here using real time hot-spot observations from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) satellite instrument as well as from MODIS hot-spots. Simulations of smoke dispersion are performed with the FLEXPART-WRF model and particulate matter emissions are calculated directly from the observed FRP. The modeled smoke plumes are compared with smoke stereo-heights from the Multiangle Imaging Spectroradiometer (MISR) instrument and the sensitivities to atmospheric and modeling parameters are examined. Driving the simulations with high resolution meteorology (4 4 km) and using geostationary satellite data to identify the hot spots allows the description of local scale features that govern smoke dispersion. The long-range transport of smoke is found to be favored over the complex coastline environment of Greece due to the abrupt changes between land and marine planetary boundary layers (PBL) and the decoupling of smoke layers from the surface.

fire smoke↗
NASA NTRSDate not supplied

Methodology for Validation of Unsteady Pressure-Sensitive Paint Measurements using Pressure Transducers

Unsteady Pressure-Sensitive Paint (uPSP) technology provides simultaneously high-frequency and high-spatial resolution surface pressure fluctuation measurements. It has demonstrated value in wind tunnel test facilities to evaluate aerospace vehicle unsteady aerodynamics. The calibration of measurements from digital pixel counts into corresponding physical units relies on both lab-based and in-situ calibration steps, and its validation is both challenging and critical to adoption of the technology. We present a rigorous methodology for direct comparison of uPSP-derived surface pressure measurements to reference pressure transducer data based on a model of the camera-based measurement error, a model of the flow field pressure signal, and their propagation through the uPSP data reduction algorithm. The methodology is demonstrated on data collected in a small-scale, subsonic wind tunnel test as well as data collected from a Space Launch System unsteady aerodynamics characterization test at the NASA Ames Research Center Unitary Plan Wind Tunnel Complex.

Pressure-Sensitive Paint↗
NASA NTRSDate not supplied

Validating uPSP with Transducers

Unsteady Pressure-Sensitive Paint (uPSP) technology provides simultaneously high-frequency and high-spatial resolution surface pressure fluctuation measurements. It has demonstrated value in wind tunnel test facilities to evaluate aerospace vehicle unsteady aerodynamics. The calibration of measurements from digital pixel counts into corresponding physical units relies on both lab-based and in-situ calibration steps, and its validation is both challenging and critical to adoption of the technology. We present a rigorous methodology for direct comparison of uPSP-derived surface pressure measurements to reference pressure transducer data based on a model of the camera-based measurement error, a model of the flow field pressure signal, and their propagation through the uPSP data reduction algorithm. The methodology is demonstrated on data collected in a small-scale, subsonic wind tunnel test as well as data collected from a Space Launch System unsteady aerodynamics characterization test at the NASA Ames Research Center Unitary Plan Wind Tunnel Complex.

Pressure-Sensitive Paint↗
NASA NTRSDate not supplied

Statistical Wind-Tunnel Experimentation Advancements for eVTOL Aircraft Aero-Propulsive Model Development

Compared to many conventional aircraft, electric vertical takeoff and landing (eVTOL) vehicle designs present greater aero-propulsive complexity and many interacting factors requiring development of novel testing and model development strategies. This paper describes several advances in statistically-rigorous experiment design methods for wind-tunnel testing of eVTOL aircraft enabling accurate aero-propulsive model development. The techniques are applied to a subscale tilt-rotor eVTOL aircraft configuration built at NASA Langley Research Center that has 24 independent control effectors. The general testing approach, leveraging design of experiments and response surface methodology, is outlined and followed by a detailed description of the experiment design and modeling approach employed for different phases of the wind-tunnel testing. A new approach for gravitational tare modeling is developed and validated, which reduces the required wind-tunnel test time by nearly 50%. A mostly automated trim envelope determination approach is formulated to rapidly determine the aircraft flight envelope to inform transition testing for aero-propulsive model development. Efficient powered-airframe characterization tests are conducted using a nested I-optimal design approach formulated in previous simulation-based testing. The collective wind-tunnel testing approach facilitates rapid characterization of eVTOL aircraft configurations and produces accurate aero-propulsive models that can be used in flight dynamics simulations.

design of experiments↗
DOE OSTI2020

Contaminant Migration Modeling to Support the In Situ Decommissioning of Hardened Nuclear Facilities - 20215

The decommissioning of hardened nuclear facilities provides a unique challenge in balancing current worker dose with exposure to future receptors and media. In situ decommissioning (ISD) is a cost-effective and safe option for the closure of such facilities but requires a detailed understanding of the threat that radioactive and hazardous constituents remaining in the facilities pose to groundwater and surrounding surface water. The Savannah River Site (SRS) extensively employs contaminant migration models within the decommissioning process to develop site-specific removal, grouting, and monitoring strategies in support of safe and effective final end states of nuclear facilities. This work discusses the use of contaminant migration modeling in the ISD process at SRS, including specific examples from the closure of P-Reactor and R-Reactor, Building 235-F, F-Canyon and F/H Laboratory Complex. Migration models developed during the closure of P-Reactor and R-Reactor were divided into four source areas: the reactor vessel, the process area, the disassembly basin, and the purification wing. Each source area was assigned a specific inventory and migration pathway and was then subject to varying hypothetical removal, capping, and grouting schemes. Modeling indicated that groundwater Maximum Concentration Limits (MCLs) might be exceeded if no action was taken for eleven and ten constituents of concern (COCs) at P-Reactor and R-Reactor, respectively. These exceedances could occur in as few as 200 to 500 years. Alternatively, the migration modeling demonstrated that the selected ISD actions reduced contaminant mobility which allowed for significant radioactive decay and resulted in fewer predicted exceedances of groundwater MCLs (five COCs for P-Reactor; eight COCs for R-Reactor). In response to the modeling results for P-Reactor and R-Reactor, effectiveness monitoring programs were developed to target contaminants, specific to each reactor, that may migrate to groundwater. Contaminant migration modeling also revealed that roof collapse was a large factor in the release of COCs to the environment, giving rise to roof improvements, and an inspection and vegetation control program to ensure roof stability over time. Contaminant migration modeling is also aiding in the closure planning for hardened facilities in F Area. Building 235-F housed the Actinide Billet Line, which produced Np-237 billets for irradiation in SRS reactors, and the Plutonium Fuel Form (PuFF) facility that produced Pu-238 heat sources for the space program. As a result of these missions, areas within Building 235-F contain considerable residual amounts of both Pu-238 and Np-237. Contaminant migration modeling of Building 235-F was originally performed in 2012 to identify the feasibility of ISD and the amount of radioactive material removal required to prevent the exceedance of groundwater MCLs. The original model indicated that a 60% reduction in the PuFF facility Pu-238 inventory could keep groundwater concentrations below standards, while Np-237 did not pose a threat to groundwater. However, updates to the model with an emphasis on source impact pathways revealed that, due to the orientation of the source areas relative to groundwater flow, no amount of reasonable removal of Pu-238 would keep groundwater concentrations below MCLs and that Np-237 could be a large contributor to localized MCL exceedances. With this insight, the refined 2019 model is being used to assist in the development of grouting plans specific to each facility source area, where bentonite may be utilized to slow the migration of Pu-238 and its daughter products from the PuFF facility and a reducing grout may decrease Np-237 transport by ensuring the nuclide remains in the less mobile +IV oxidation state. The beginning phases of contaminant migration modeling are underway for F-Canyon and associated facilities using lessons learned from the Reactors and Building 235-F. A contaminant migration pathway, similar to the pathway used for the reactor vessels, is being developed for the hot and warm canyons. The F-Canyon inventories are being spatially refined to identify specific and localized source zones, comparable to Building 235-F, that may impact groundwater. Contaminant migration modeling provides vital input within the ISD process, from facility investigation to post closure effectiveness monitoring, making it a valuable tool in the development of safe and effective end states for hardened nuclear facilities at SRS. (authors)

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

Simulation experiments for gamma-ray mapping of planetary surfaces: Scattering of high-energy neutrons

The concentration and distribution of certain elements in surface layers of planetary objects specify constraints on models of their origin and evolution. This information can be obtained by means of remote sensing gamma-ray spectroscopy, as planned for a number of future space missions, i.e., Mars, Moon, asteroids, and comets. To investigate the gamma-rays made by interactions of neutrons with matter, thin targets of different composition were placed between a neutron-source and a high-resolution germanium spectrometer. Gamma-rays in the range of 0.1 to 8 MeV were accumulated. In one set of experiments a 14-MeV neutron generator using the T(d,n) reaction as neutron-source was placed in a small room. Scattering in surrounding walls produced a spectrum of neutron energies from 14 MeV down to thermal. This complex neutron-source induced mainly neutron-capture lines and only a few scattering lines. As a result of the set-up, there was a considerable background of discrete lines from surrounding materials. A similar situation exists under planetary exploration conditions: gamma-rays are induced in the planetary surface as well as in the spacecraft. To investigate the contribution of neutrons with higher energies, an experiment for the measurement of prompt gamma radiation was set up at the end of a beam-line of an isochronous cyclotron.

Brueckner, J.↗
NASA NTRS1999

Development of Response Surface Models for Rapid Analysis and Multidisciplinary Optimization of Launch Vehicle Design Concepts

Multidisciplinary design optimization (MDO) is an important step in the design and evaluation of launch vehicles, since it has a significant impact on performance and lifecycle cost. The objective in MDO is to search the design space to determine the values of design parameters that optimize the performance characteristics subject to system constraints. Vehicle Analysis Branch (VAB) at NASA Langley Research Center has computerized analysis tools in many of the disciplines required for the design and analysis of launch vehicles. Vehicle performance characteristics can be determined by the use of these computerized analysis tools. The next step is to optimize the system performance characteristics subject to multidisciplinary constraints. However, most of the complex sizing and performance evaluation codes used for launch vehicle design are stand-alone tools, operated by disciplinary experts. They are, in general, difficult to integrate and use directly for MDO.

Unal, Resit↗
NASA NTRS2001

Experiments Developed to Study Microgravity Smoldering Combustion

The overall objective of the Microgravity Smoldering Combustion (MSC) research program is to understand and predict smoldering combustion under normal and microgravity (near-zero-gravity) conditions to help prevent and control smolder-originated fires, in both environments. Smoldering is defined as a nonflaming, self-sustaining, propagating, exothermic surface reaction. If a material is sufficiently permeable, smoldering is not confined to its outer surface, but can propagate as a reaction wave through the interior of the material. The MSC program will accomplish its goals by conducting smolder experiments on the ground and in a space-based laboratory, and developing theoretical models of the process. Space-based experiments are necessary because smoldering is a very slow process and, consequently, its study in a microgravity environment requires extended periods of time that can only be achieved in space. Smoldering can occur in a variety of processes ranging from the smolder of porous insulating materials to underground coal combustion. Many materials can sustain smoldering, including wood, cloth, foams, tobacco, other dry organic materials, and charcoal. The ignition, propagation, transition to flaming, and extinction of the smolder reaction are controlled by complex, thermochemical mechanisms that are not well understood. As with many forms of combustion, gravity affects the availability of the oxidizer and the transport of heat, and therefore, the rate of combustion. The smoldering combustion of porous materials has been studied both experimentally and theoretically, usually in the context of fire safety. Smoldering encompasses a number of fundamental processes, including heat and mass transfer in a porous media; endothermic pyrolysis of combustible material; ignition, propagation, and extinction of heterogeneous exothermic reactions at the solid-gas pore interface; and the onset of gas phase reactions (flaming) from existing surface reactions. Smoldering presents a serious fire risk because the combustion can propagate slowly in a material's interior and go undetected for long periods of time. It typically yields a substantially higher conversion of fuel to toxic compounds than does flaming (though more slowly), and may undergo a sudden transition to flaming.

Vergilii, Franklin↗
DOE OSTI2023

Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗