Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equilibrium programming”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Parametric study and speciation analysis of rare earth precipitation using oxalic acid in a chloride solution system

Oxalic acid precipitation is a common step in the purification of rare earth elements (REE) from a concentrated pregnant leach solution (PLS). However, the presence of contaminants such as Al, Fe, and Ca in given amounts decreases the REE precipitation efficiency and product purity while also increasing the amount of oxalic acid needed to maximize recovery. As such, a statistically designed test program was performed to identify the optimal conditions necessary for a relatively low REE content PLS containing elevated concentrations of contaminant ions. The performance objective was maximization of REE precipitation efficiency while minimizing the oxalic acid dosage. A central composite design was utilized to quantify performance impacts and identify the ultimate set of parameter values for oxalic acid dosage, Fe(III) contamination concentration, solution pH, and reaction temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicate that increasing the oxalic acid concentration from 0 g/L to 80 g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. Here, to further the understanding of the REE-oxalate system, a fundamental solution chemistry study was performed using the equilibrium constants of the reactions. The study resulted in the development of oxalate speciation diagrams and provided an analysis of the REE precipitation characteristics at various oxalate anion concentrations and Fe(III) contamination levels using MINTEQ software. The dominant Fe(III) species in the solution system were found to be Fe-(C 2 O 4 ) 3 3- , Fe-(C 2 O 4 ) 2- , and Fe-(C 2 O 4 ) + , which consume the majority of the oxalate anions. The simulated model was found to be in agreement with the experimental findings and helped to explain the adverse impact of increased iron concentrations on REE precipitation efficiency.

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

Thermal Scattering Law Data Development for Paraffin Wax

Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).

Nuclear Criticality Safety Program (NCSP)↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Fuel Fabrication Specification Impact Analysis for NBSR LEU Conversion

As part of a national initiative to enhance nuclear security and reduce proliferation risks, significant efforts have been undertaken by the National Nuclear Security Administration Material Management and Minimization Office of Reactor Conversion Program to convert U.S. high performance research reactors (USHPRRs) from the use of highly enriched uranium (HEU) to low-enriched uranium (LEU), including the National Bureau of Standards Reactor (NBSR). The current plan is to procure LEU fuel assemblies from commercial fabricators according to fuel specifications tailored for each USHPRR. The analysis conducted at Brookhaven National Laboratory was part of an effort to identify the sources of uncertainty in the fuel specifications that may impact the performance of the NBSR core after its conversion and, in particular, to assess the range of acceptable tolerance limits from the perspective of core safety and reactor performance. Using the stochastic neutronics code MCNP 6.2, the variations in important NBSR neutronics characteristics were analyzed as a function of the specification parameters independently and in combination. The important NBSR specification parameters analyzed were the fuel isotopic composition, the amount of impurity content in cladding, the fuel plate thickness, and the fuel element 235U mass loading. The range of variation of each specification parameter was based on the technical specification limit or available as-fabricated assay data and uncertainties. The NBSR neutronics characteristics selected for analysis were the reactor reactivity characteristics at equilibrium and the equilibrium fuel cycle length. Results show that with variations in the fabrication parameters of the as-fabricated U-10Mo fuel within the specification limitations, the excess reactivity of the NBSR LEU core remains well below the 15% Δk/k technical specification limit, and the shutdown margin is always significantly greater than the required 0.68% Δk/k. This ensures that the NBSR can be operated safely and reliably shut down for all analyzed cases within the specified fabrication limits after the LEU conversion. In the prototypic case, the fuel cycle length was 1.5 days longer than the targeted 38.5 days. In a credible worst-case scenario, where all low-reactivity parameters were combined, the fuel cycle length was reduced to 35.5 days, which is still considered manageable for reactor operations. Variations in cycle length are primarily driven by changes in 235U loading, with other parameters having secondary effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2003

Flight Planning for the International Space Station - Levitation Observation of Dendrite Evolution in Steel Ternary Alloy Rapid Solidification (LODESTARS)

During rapid solidification, a molten sample is cooled below its equilibrium solidification temperature to form a metastable liquid. Once nucleation is initiated, growth of the solid phase proceeds and can be seen as a sudden rise in temperature. The heat of fusion is rejected ahead of the growing dendrites into the undercooled liquid in a process known as recalescence. Fe-Cr-Ni alloys may form several equilibrium phases and the hypoeutectic alloys, with compositions near the commercially important 316 stainless steel alloy, are observed to solidify by way of a two-step process known as double recalescence. During double recalescence, the first temperature rise is associated with formation of the metastable ferritic solid phase with subsequent conversion to the stable austenitic phase during the second temperature rise. Selection of which phase grows into the undercooled melt during primary solidification may be accomplished by choice of the appropriate nucleation trigger material or by control of the processing parameters during rapid solidification. Due to the highly reactive nature of the molten sample material and in order to avoid contamination of the undercooled melt, a containerless electromagnetic levitation (EML) processing technique is used. In ground-based EML, the same forces that support the weight of the sample against gravity also drive convection in the liquid sample. However, in microgravity, the force required to position the sample is greatly reduced, so convection may be controlled over a wide range of internal flows. Space Shuttle experiments have shown that the double recalescence behavior of Fe-Cr-Ni alloys changes between ground and space EML experiments. This program is aimed at understanding how melt convection influences phase selection and the evolution of rapid solidification microstructures.

Flemings, Merton C.↗
NASA NTRS2000

X-38 Experimental Aerothermodynamics

The X-38 program seeks to demonstrate an autonomously returned orbital test flight vehicle to support the development of an operational Crew Return Vehicle for the International Space Station. The test flight, anticipated in 2002, is intended to demonstrate the entire mission profile of returning Space Station crew members safely back to earth in the event of medical or mechanical emergency. Integral to the formulation of the X-38 flight data book and the design of the thermal protection system, the aerothermodynamic environment is being defined through a synergistic combination of ground based testing and computational fluid dynamics. This report provides an overview of the hypersonic aerothermodynamic wind tunnel program conducted at the NASA Langley Research Center in support of the X-38 development. Global and discrete surface heat transfer force and moment, surface streamline patterns, and shock shapes were measured on scaled models of the proposed X-38 configuration in different test gases at Mach 6, 10 and 20. The test parametrics include angle of attack from 0 to 50 degs, unit Reynolds numbers from 0.3 x 10 (exp 6) to 16 x 10 (exp 6)/ ft, rudder deflections of 0, 2, and 5 deg. and body flap deflections from 0 to 30 deg. Results from hypersonic aerodynamic screening studies that were conducted as the configuration evolved to the present shape at, presented. Heavy gas simulation tests have indicated that the primary real gas effects on X-38 aerodynamics at trim conditions are expected to favorably influence flap effectiveness. Comparisons of the experimental heating and force and moment data to prediction and the current aerodynamic data book are highlighted. The effects of discrete roughness elements on boundary layer transition were investigated at Mach 6 and the development of a transition correlation for the X-38 vehicle is described. Extrapolation of ground based heating measurements to flight radiation equilibrium wall temperatures at Mach 6 and 10 were made and generally compared to within 50 deg F of flight prediction.

Horvath, Thomas J.↗
DOE OSTI2023

An Incremental Gradient Method for Optimization Problems With Variational Inequality Constraints

We consider minimizing a sum of agent-specific nondifferentiable merely convex functions over the solution set of a variational inequality (VI) problem in that each agent is associated with a local monotone mapping. This problem finds an application in computation of the best equilibrium in nonlinear complementarity problems arising in transportation networks. We develop an iteratively regularized incremental gradient method where at each iteration, agents communicate over a directed cycle graph to update their solution iterates using their local information about the objective and the mapping. The proposed method is single-timescale in the sense that it does not involve any excessive hard-to-project computation per iteration. We derive nonasymptotic agent-wise convergence rates for the suboptimality of the global objective function and infeasibility of the VI constraints measured by a suitably defined dual gap function. Finally, the proposed method appears to be the first fully iterative scheme equipped with iteration complexity that can address distributed optimization problems with VI constraints over cycle graphs.

convergence↗
DOE OSTI2023

Pros and Cons Analysis of HALEU Utilization in Example Fuel Cycles

The Systems Analysis and Integration campaign assessed the pros and cons of high-assay low-enriched uranium (HALEU) utilization in advanced reactors and associated fuel cycles. The assessment was done for three example fuel cycles at equilibrium states: once-through, limited recycle, and continuous recycle (CR) starting with HALEU. Front- and back-end fuel cycle parameters and the Levelized Cost of Fuel (LCF), which is the Levelized Cost of Electricity excluding reactor cost, of the three example fuel cycles were calculated using a single Analysis Example Reactor. The pros and cons of HALEU utilization were assessed by normalizing the fuel cycle parameters and LCF to a unit of electricity generation (GWe-year) and comparing them with a Basis of Comparison. In this study, a sodium-cooled reactor with sodium-bonded metallic fuel having a burnup of ~100 GWd/t was used as the Analysis Example Reactor because its technology readiness level is high, and the burnup and fuel enrichment are in the middle of those ranges of advanced reactor concepts that are under development. The current once-through Light Water Reactors (OT-LWRs) with <5% low-enriched uranium and 50 GWd/t burnup were used as the Basis of Comparison. In addition, a series of sensitivity analyses was conducted by varying burnup, enrichment, fuel forms, and reactor types to capture the design variations in two once-through Advanced Reactor Demonstration Program (ARDP) reactors, Natrium with sodium-free metallic fuel having a burnup of ~150 GW/t and Xe-100 with Tristructural-Isotropic (TRISO) pebble fuel having a burnup of ~168 GWd/t.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2012

A Peculiar Class of Debris Disks from Herschel/DUNES: A Steep Fall Off in the Far Infrared

Context. The existence of debris disks around old main sequence stars is usually explained by continuous replenishment of small dust grains through collisions from a reservoir of larger objects. Aims. We present photometric data of debris disks around HIP 103389 (HD199260), HIP 100350 (HN Peg, HD206860), and HIP 114948 (HD 219482), obtained in the context of our Herschel Open TIme Key Program DUNES (DUst around NEarby Stars). Methods. We used Herschel/PACS to detect the thermal emission of the three debris disks with a 30 sigma sensitivity of a few mJy at l00 micron and 160 micron. In addition, we obtained Herschel/PACS photometric data at 70 micron for HIP 103389. These observations are complemented by a large variety of optical to far-infrared photometric data. Two different approaches are applied to reduce the Herschel data to investigate the impact of data reduction on the photometry. We fit analytical models to the available spectral energy distribution (SED) data using the fitting method of simulated therma1 annealing as well as a classical grid search method. Results. The SEDs of the three disks potentially exhibit an unusually steep decrease at wavelengths >= 70 micron. We investigate the significance of the peculiar shape of these SEDs and the impact on models of the disks provided it is real. Using grain compositions that have been applied successfully for modeling of many other debris disks, our modeling reveals that such a steep decrease of the SEDs in the long wavelength regime is inconsistent with a power-law exponent of the grain size distribution -3.5 expected from a standard equilibrium collisional cascade. In contrast, a steep grain size distribution or, alternatively an upper grain size in the range of few tens of micrometers are implied. This suggests that a very distinct range of grain sizes would dominate the thermal. emission of such disks. However, we demonstrate that the understanding of the data of faint sources obtained with Herschel is still incomplete and that the significance of our results depends on the version of the data reduction pipeline used. Conclusions. A new mechanism to produce the dust in the presented debris disks, deviations from the conditions required for a standard equilibrium collisional cascade (grain size exponent of -3.5), and/or significantly different dust properties would be necessary to explain the potentially steep SED shape of the three debris disks presented.

Ertel, S.↗
NASA NTRS1989

Photo-sensitized ignition of hydrogen/oxygen mixtures for hypersonic flight vehicles

An approach toward improving supersonic ignition and combustion efficiency was investigated under a feasibility program with NASA-Ames. The idea is to selectively irradiate targetted species in the combustion zone. The ensuing photodissociation reactions produce highly reactive radicals which modify the gas phase kinetics. The overall objective of this research was to demonstrate photosensitized ignition of H2/O2 mixtures, and to gain a preliminary understanding of the ignition mechanism, especially as it pertains to the concept of supersonic combustion. To this end, a series of experiments were performed with a model sensitizer, NO2. At subatmospheric pressures and moderate temperatures (above 100 F), mixtures of H2/O2/NO2 are ignitable by a pulsed eximer laser at 308 nm, 115 to 200 mJ. Ignitability was found to be a function of NO2 pressure, temperature, laser energy per pulse, and laser lens focal length (fluence at focal volume). Mixtures could not be ignited in the absence of the sensitizer. While increasing temperature improved sensitized ignition in the lower temperature range, a retarding effect was observed at the highest temperatures. To more thoroughly understand these thermal effects, the UV-Vis and IR spectra of NO2 and H2/O2/NO2 mixtures were studied over the temperature range 70 to 900 F. Initially, increasing temperature shifts the N2O4 in a reversible reaction to 2 NO2 equilibrium toward NO2, favoring photochemical ignition. However, at still higher temperature, NO2 converts to non-absorbing NO + 1/2 O2, thus preventing photosensitized ignition.

Lavid, Moshe↗
NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗
NASA NTRS2003

MRI of the lung gas-space at very low-field using hyperpolarized noble gases

In hyperpolarized (HP) noble-gas magnetic resonance imaging, large nuclear spin polarizations, about 100,000 times that ordinarily obtainable at thermal equilibrium, are created in 3He and 129Xe. The enhanced signal that results can be employed in high-resolution MRI studies of void spaces such as in the lungs. In HP gas MRI the signal-to-noise ratio (SNR) depends only weakly on the static magnetic field (B(0)), making very low-field (VLF) MRI possible; indeed, it is possible to contemplate portable MRI using light-weight solenoids or permanent magnets. This article reports the first in vivo VLF MR images of the lungs in humans and in rats, obtained at a field of only 15 millitesla (150 Gauss).

NASA Discipline Life Sciences Technologies↗
NASA NTRS2010

Gas Flow Detection System

This system provides a portable means to detect gas flow through a thin-walled tube without breaking into the tubing system. The flow detection system was specifically designed to detect flow through two parallel branches of a manifold with only one inlet and outlet, and is a means for verifying a space shuttle program requirement that saves time and reduces the risk of flight hardware damage compared to the current means of requirement verification. The prototype Purge Vent and Drain Window Cavity Conditioning System (PVD WCCS) Flow Detection System consists of a heater and a temperature-sensing thermistor attached to a piece of Velcro to be attached to each branch of a WCCS manifold for the duration of the requirement verification test. The heaters and thermistors are connected to a shielded cable and then to an electronics enclosure, which contains the power supplies, relays, and circuit board to provide power, signal conditioning, and control. The electronics enclosure is then connected to a commercial data acquisition box to provide analog to digital conversion as well as digital control. This data acquisition box is then connected to a commercial laptop running a custom application created using National Instruments LabVIEW. The operation of the PVD WCCS Flow Detection System consists of first attaching a heater/thermistor assembly to each of the two branches of one manifold while there is no flow through the manifold. Next, the software application running on the laptop is used to turn on the heaters and to monitor the manifold branch temperatures. When the system has reached thermal equilibrium, the software application s graphical user interface (GUI) will indicate that the branch temperatures are stable. The operator can then physically open the flow control valve to initiate the test flow of gaseous nitrogen (GN2) through the manifold. Next, the software user interface will be monitored for stable temperature indications when the system is again at thermal equilibrium with the test flow of GN2. The temperature drop of each branch from its "no flow" stable temperature peak to its stable "with flow" temperature will allow the operator to determine whether a minimum level of flow exists. An alternative operation has the operator turning on the software only long enough to record the ambient temperature of the tubing before turning on the heaters and initiating GN2 flow. The stable temperature of the heated tubing with GN2 flow is then compared with the ambient tubing temperature to determine if flow is present in each branch. To help quantify the level of flow in the manifolds, each branch will be bench calibrated to establish its thermal properties using the flow detection system and different flow rates. These calibration values can then be incorporated into the software application to provide more detailed flow rate information.

Moss, Thomas↗
NASA NTRS2013

Implicit Coupling Approach for Simulation of Charring Carbon Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption

Ablation↗
NASA NTRS1986

Modeling the temporal and spatial variations of the vertical structure of Jupiter's atmosphere using observations of the 3-0 hydrogen quadrupole lines

An observational program was established in 1983 to monitor the spatial and temporal variations in the Jovian atmosphere over short and long time scales. The program involves tracking several different longitudes as they rotate around the planet from one limb to another. This tracking experiment was done at many different wavelengths including the 3-0 S(1) and S(0) hydrogen quadrupole lines as well as several broad band methane absorptions. The June 1983 hydrogen quadrupole data was reduced and equivalent widths were measured for approximately 25 east-west positions across the planet at 7 different latitudes for both wavelengths. The data for the South Tropical Zone (20 deg. S) was modeled extensively and the effects of the various model parameters on the value of the calculated equivalent widths of both lines was measured as a longitude rotated from the east (or morning) limb to the west (or evening) limb. The value of the equivalent width is also quite sensitive to the height of the NH3 cloud top and to the value used for the single scattering albedo. A combination of these parameters changing on a diurnal time scale could also explain these observations. This gradual increase from one limb to the other appears in the data for both the North and South Equatorial Belts as well as the equatorial region and the North Tropical Zone. Models that used only normal hydrogen and models that used only equilibrium hydrogen were studied.

Cunningham, C. C.↗
DOE OSTI2020

Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magnetohydrodynamic codes M3D-C 1 , JOREK, and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with a rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. We report the benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C 1 simulations and from the linear phase of axisymmetric nonlinear JOREK, NIMROD, and M3D-C 1 simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared, and an excellent agreement is found regarding the plasma location and plasma currents, as well as eddy and halo currents in the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗