Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material flow analysis”

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

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

At least 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Design and demonstration of a direct air capture system with moisture-driven CO 2 delivery into aqueous medium

Two direct air capture (DAC) systems were designed and demonstrated to passively capture CO 2 from ambient air and use moisture to release the CO2 into an alkaline medium. A bench-scale system delivering ∼1 g CO 2 d –1 was demonstrated in a laminar flow hood, and a small pilot-scale system that could deliver ∼100 g CO 2 d –1 was operated outdoors in a 4.2 m 2 raceway pond. Novel elongated mesh-tube packets containing anion-exchange resin (AER) beads were found to reduce drying and CO 2 loading time 4.3-fold compared to larger mesh bags. Technoeconomic analysis (TEA) estimates the cost of capturing CO 2 into an alkaline solution, suitable for cultivating photosynthetic microorganisms, to be $\$229$ per tonne for a practical scenario based on current results and $\$72$ per tonne for an aspirational scenario considering improvements to sorbent capacity, hydrophobicity, and sorbent lifetime. TEA further estimates an additional $\$110$ per tonne to extract CO 2 from solution, purify it, and compress it to 15 MPa, suitable for sequestration. Furthermore, moisture-driven processes have the potential to use up to 87% less energy than thermal and/or vacuum swing DAC by using energy from water evaporation.

09 BIOMASS FUELS↗
DOE OSTI2025

B‐Staging and Crosslinking of Polycarbosilane at Room Temperature: Cure Mechanism and Properties

The viscosity of polycarbosilane (PCS) polymers is not advantageous for free forming. Various fillers and heat are used to obtain a formable paste. Due to the low yield stress, structures and preforms tend to slump or resin will flow out, especially during the curing of the polymer. B-staging of PCS allows for a more stable structure from room temperature until the full cure of the allyl groups occurs with heat. Additionally, hydrosilation is an effective means of crosslinking at room temperature and controlling viscosity. A network structure was formed in SMP-10 using silane and a vinyl-based crosslinker to bridge each polymer chain. Pt addition catalysts were added to enhance the increase in viscosity to make a lightly crosslinked gel to aid in thickening and forming and to improve the ceramic yield at 1000°C. In conclusion, this enables several options for controlling rheology and improving the properties of preceramic polymers to avoid slumping during curing.

inorganic polymers↗
DOE OSTI2023

Development of an All-Aqueous Thermally Regenerative Redox Flow Battery to Support Fossil Fuel Assets

Low-temperature thermal energy, a largely untapped resource, holds significant promise for large-scale electrical power generation globally. Various stationary sources, including industrial entities and thermal power plants, emit considerable low-temperature heat that currently remains unutilized. This energy is often overlooked because its low temperature renders it unsuitable for efficient power generation using conventional methods. However, current research is exploring diverse technologies capable of converting low-temperature heat into grid-scale power, aiming to enhance grid efficiency, further decarbonization initiatives, and facilitate a shift toward more decentralized power systems. One such innovative technology is the thermally regenerative battery (TRB), noted for its high power and energy densities compared to similar technologies, positioning it as a potential game-changer in power generation. TRBs integrate two scalable and well-established unit operations: a redox flow battery and a distillation column. This integration suggests that once an effective TRB chemistry is established, the pathway to commercialization could be expedited. The copper-based thermally regenerative ammonia battery (Cu aq -TRAB) stands out as the first TRB that circumvents electrodeposition/dissolution reactions, stabilizing Cu(I) and Cu(II) within the electrolyte and maintaining stability of all electroactive species in an aqueous phase. This stabilization has led to improvements in coulombic efficiency, open circuit potential, and copper solubility, thereby enhancing power density, energy density, and overall energy efficiency. Preliminary tests were conducted to determine the effects of various electrolyte species on the performance metrics of the battery, both theoretically and experimentally. These tests revealed that the solubility of copper in the Cu aq -TRAB electrolyte was constrained by the Cu(I)-NH 3 complex. Adjusting the background electrolyte to 5 M NH4Br and the ligand concentration to 4 M NH 3 enabled the copper concentration to reach a maximum of 0.6 M. This modification led to an estimated theoretical maximum energy density of 9.5 Wh L -1 for the Cu aq -TRAB. Additionally, full cell testing indicated a tradeoff between peak power and energy density with varying copper and ammonia concentrations. Increasing the applied current density during discharge linearly raised the average power output, with a minimal reduction in energy density due to a balance between higher ohmic overpotential and reduced time for undesirable ammonia crossover. Furthermore, a comprehensive numerical sensitivity analysis of the complete Cu aq -TRAB system was performed. This analysis aimed to assess how the battery and the distillation column responded to changes in system input parameters, providing insights into optimal research directions for enhancing system performance. The analysis revealed that at room temperature, battery power was significantly more sensitive to ohmic losses than to mass transfer, with reaction rates having minimal impact. This trend continued even at higher temperatures. Also, the thermal energy required for ammonia separation was studied, showing that increased temperatures generally reduced energy requirements, except in low-pressure scenarios above 65 °C. An investigation into membrane performance in the Cu aq -TRAB was undertaken, given the significant impact of ammonia transport control and ohmic losses on system performance. Various membranes were evaluated to identify key performance metrics. Among the tested membranes, Selemion CMVN exhibited the highest performance, with a peak power density of 84 mW cm -2 and average values of 26 ± 6.8 mW cm -2 for power density and 2.9 Wh L -1 for energy density at an applied current density of 50 mA cm -2 . An economic assessment indicated a levelized cost of storage at $410 per MWh under optimal conditions, highlighting the commercial potential of the Cu aq -TRAB when utilizing cost-effective, readily available materials.

25 ENERGY STORAGE↗
NASA NTRS1998

An Assessment of Molten Metal Detachment Hazards During Electron Beam Welding in Space

The safety issue has been raised with regards to potential molten metal detachments from the weld pool and cold filler wire during electron beam welding in space. This investigation was undertaken to evaluate if molten metal could detach and come in contact with astronauts and burn through the fabric of the astronauts' Extravehicular Mobility Unit (EMU) during electron beam welding in space. Molten metal detachments from either the weld/cut substrate or weld wire could present harm to a astronaut if the detachment was to burn through the fabric of the EMU. Theoretical models were developed to predict the possibility and size of the molten metal detachment hazards during the electron beam welding exercises at Low Earth Orbit (LEO). The primary molten metal detachment concerns were those cases of molten metal separation from the metal surface due to metal cutting, weld pool splashing, entrainment and release of molten metal due to filler wire snap-out from the weld puddle, and molten metal accumulation and release from the end of the weld wire. Some possible ways of obtaining molten metal drop detachments would include an impulse force, or bump, to the weld sample, cut surface, or filler wire. Theoretical models were developed for these detachment concerns from principles of impact and kinetic energies, surface tension, drop geometry, surface energies, and particle dynamics. The surface tension represents the force opposing the liquid metal drop from detaching whereas the weight of the liquid metal droplet represents a force that is tending to detach the molten metal drop. Theoretical calculations have indicated that only a small amount of energy is required to detach a liquid metal drop; however, much of the energy of an impact is absorbed in the sample or weld plate before it reaches the metal drop on the cut edge or surface. The tendency for detachment is directly proportional to the weld pool radius and metal density and inversely proportional to the surface tension of the liquid metal. For a detachment the initial kinetic energy of the weld pool with respect to the plate has to exceed the energy to form the extra surface required for the detachment of the pool. The difficulty is in transferring the energy from the point of impact through the plate and sample to the cut edge. It is likely that not all of the kinetic energy is available for detaching the pool; some may be sequestered in weld pool oscillations. The coefficient of restitution for the collision will be lower than one if irreversible deformation, for example plastic flow deformation, takes place during the collision. Thus determining the amount of energy from an impact that actually reaches the molten metal droplet is critical. Various molten metal detachment scenarios were tested experimentally in an enclosed vacuum chamber using the Ukrainian Universal Hand Tool, an electron beam welder designed for space welding. The experimental testing was performed in a 4 ft. X 4 ft. vacuum chamber at Marshall Space Flight Center, evacuated to vacuum levels of at least 50 microTorr, and also some welding garment material was utilized to observe the effect of the molten metal detachments on the material. A "carillon" apparatus consisting of four pendulum hammer strikers, each weighing approximately 3.65 lbs, raised to predetermined specific heights was used to apply an impact force to the weld sample/plate during electron beam welding and cutting exercises. The strikers were released by switching on an electric motor to rotate a pin holding wires retaining the strikers at desired heights. The specimens were suspended so as to be free to respond to the blows with a sudden velocity increment. The specimens were mounted on a hinged plate for minimizing effective mass with the option to fasten it down so as to raise its effective mass closer to that anticipated for an actual space welding scenario. Measurements were made of the impact energy and the horizontal fling distances of the detached metal drops. It was not particularly easy to generate the detachments for this experiment. This document presents the details of the theoretical modeling effort and a summary of the experimental effort to measure molten metal drop detachments from terrestrial electron beam welding in the enclosed vacuum chamber. The results of the experimental effort have shown that molten metal detachments can occur from the sample/weld plate only if a sufficiently large impact force is applied to the weld plate. A "weld pool detachment parameter" was determined to indicate whether detachment would occur. Detachment can be either full or partial (dripping), Partial detachment means that the weld pool detached from one side of the liquid-solid boundary so as to leave a hole at the puddle site but remained attached over part of the liquid-solid boundary and dripped down the plate with no fully detached material detected. Full detachment, however, does not necessarily mean that the whole pool fully detached; in some cases only a smaller portion of the pool detached, the remainder dripping down the plate. The weld pool detachment parameter according to theory and according to the empirical data allows a determination of whether full detachments might occur. Theoretical calculations indicated titanium alloy would be the most difficult from which to detach molten metal droplets followed by stainless steel and then by aluminum. The experimental results were for the most part consistent with the theoretical analysis and predictions. The above theory is applicable to other situations as desired for assessing the potential for molten metal detachments.

Fragomeni, James M.↗
DOE OSTI2023

The Influence of Carbon Fiber Composite Specimen Design Parameters on Artificial Lightning Strike Current Dissipation and Material Thermal Damage

Previous artificial lightning strike direct effect research has examined a broad range of specimen design parameters. No works have studied how such specimen design parameters and electrical boundary conditions impact the dissipation of electric current flow through individual plies. This article assesses the influence of carbon fiber composite specimen design parameters (design parameters = specimen size, shape, and stacking sequence) and electrical boundary conditions on the dissipation of current and the spread of damage resulting from Joule heating. Thermal-electric finite element (FE) modelling is used and laboratory scale (<1 m long) and aircraft scale (>1 m long) models are generated in which laminated ply current dissipation is predicted, considering a fixed artificial lightning current waveform. The simulation results establish a positive correlation between the current exiting the specimen from a given ply and the amount of thermal damage in that ply. The results also establish that the distance to ground, from the strike location to the zero potential boundary conditions (ground), is the controlling factor which dictates the electric current dissipation in each ply. Significantly, this distance to ground is dependent on each of the specimen shape, dimensions, stacking sequence, and location of ground boundary conditions. Therefore, it is not possible to decouple current dissipation and damage from specimen design and boundary condition setup. However, it is possible to define a specimen size for a given specimen shape, stacking sequence, and waveform which limit the influence of specimen dimensions on the resulting current distribution and damage. For a rectangular specimen design which appears in literature multiple times, as 100 × 150 mm and with a stacking sequence of [45/0/-45/90] 4s , a specimen design of greater than 300 × 200 mm is required to limit the influence of specimen dimensions on current distribution and damage.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Global Stratigraphy of Venus: Analysis of a Random Sample of Thirty-Six Test Areas

The age relations between 36 impact craters with dark paraboloids and other geologic units and structures at these localities have been studied through photogeologic analysis of Magellan SAR images of the surface of Venus. Geologic settings in all 36 sites, about 1000 x 1000 km each, could be characterized using only 10 different terrain units and six types of structures. These units and structures form a major stratigraphic and geologic sequence (from oldest to youngest): (1) tessera terrain; (2) densely fractured terrains associated with coronae and in the form of remnants among plains; (3) fractured and ridged plains and ridge belts; (4) plains with wrinkle ridges; (5) ridges associated with coronae annulae and ridges of arachnoid annulae which are contemporary with wrinkle ridges of the ridged plains; (6) smooth and lobate plains; (7) fractures of coronae annulae, and fractures not related to coronae annulae, which disrupt ridged and smooth plains; (8) rift-associated fractures; and (9) craters with associated dark paraboloids, which represent the youngest 1O% of the Venus impact crater population (Campbell et al.), and are on top of all volcanic and tectonic units except the youngest episodes of rift-associated fracturing and volcanism; surficial streaks and patches are approximately contemporary with dark-paraboloid craters. Mapping of such units and structures in 36 randomly distributed large regions (each approximately 10(exp 6) sq km) shows evidence for a distinctive regional and global stratigraphic and geologic sequence. On the basis of this sequence we have developed a model that illustrates several major themes in the history of Venus. Most of the history of Venus (that of its first 80% or so) is not preserved in the surface geomorphological record. The major deformation associated with tessera formation in the period sometime between 0.5-1.0 b.y. ago (Ivanov and Basilevsky) is the earliest event detected. In the terminal stages of tessera fon-nation, extensive parallel linear graben swarms representing a change in the style of deformation from shortening to extension were formed on the tessera and on some volcanic plains that were emplaced just after, and perhaps also during the latter stages of the major compressional phase of tessera emplacement. Our stratigraphic analyses suggest that following tessera formation, extensive volcanic flooding resurfaced at least 85% of the planet in the form of the presently-ridged and fractured plains. Several lines of evidence favor a high flux in the post-tessera period but we have no independent evidence for the absolute duration of ridged plains emplacement. During this time, the net state of stress in the lithosphere apparently changed from extensional to compressional, first in the form of extensive ridge belt development, followed by the formation of extensive wrinkle ridges on the flow units. Subsequently, there occurred local emplacement of smooth and lobate plains units which are presently essentially undefortned. The major events in the latest 10% of the presently preserved history of Venus (less than 50 m.y. ago) are continued rifting and some associated volcanism, and the redistribution of eolian material largely derived from impact crater deposits.

Basilevsky, Alexander T.↗
NASA NTRS2014

Performance of a Supersonic Over-Wing Inlet with Application to a Low-Sonic-Boom Aircraft

Development of commercial supersonic aircraft has been hindered by many related factors including fuel-efficiency, economics, and sonic-boom signatures that have prevented over-land flight. Materials, propulsion, and flight control technologies have developed to the point where, if over-land flight were made possible, a commercial supersonic transport could be economically viable. Computational fluid dynamics, and modern optimization techniques enable designers to reduce the boom signature of candidate aircraft configurations to acceptable levels. However, propulsion systems must be carefully integrated with these low-boom configurations in order that the signatures remain acceptable. One technique to minimize the downward propagation of waves is to mount the propulsion systems above the wing, such that the wing provides shielding from shock waves generated by the inlet and nacelle. This topmounted approach introduces a number of issues with inlet design and performance especially with the highly-swept wing configurations common to low-boom designs. A 1.79%-scale aircraft model was built and tested at the NASA Glenn Research Center's 8-by 6-Foot Supersonic Wind Tunnel (8x6 SWT) to validate the configuration's sonic boom signature. In order to evaluate performance of the top-mounted inlets, the starboard flow-through nacelle on the aerodynamic model was replaced by a 2.3%-scale operational inlet model. This integrated configuration was tested at the 8x6 SWT from Mach 0.25 to 1.8 over a wide range of angles-of-attack and yaw. The inlet was also tested in an isolated configuration over a smaller range of angles-of-attack and yaw. A number of boundary-layer bleed configurations were investigated and found to provide a substantial positive impact on pressure recovery and distortion. Installed inlet performance in terms of mass capture, pressure recovery, and distortion over the Mach number range at the design angle-of-attack of 4-degrees is presented herein and compared to that at 0- degrees, as well as the isolated inlet configuration to highlight installation effects. Performance of the installed inlet fell below that of the isolated inlet at Mach numbers of 1.4 and greater. The installed inlet demonstrated adequate operability over the expected range of angles-of-attack and yaw, but did exhibit definite angle-ofattack and yaw limits at supersonic conditions. At each supersonic flight Mach number, performance parameters near zero yaw angle were relatively insensitive to yaw, but in general the yaw angle yielding best performance was non-zero and varied with angle-of-attack. Performance of the installed inlet is also presented as functions of angle-of-attack and yaw to highlight these effects. Distortion at the aerodynamic interface plane ranged between 10 and 25% at the inlet critical points over the range of flight Mach numbers tested and did not decrease significantly for the isolated inlet. Although these distortion levels would be considered high for operation with a turbine engine, the over-wing installation is likely not as significant a contributor as the low test Reynolds number. This is demonstrated by comparing CFD analysis of the isolated inlet at test scale with that at intermediate and full scales.

Performance↗
NASA NTRSDate not supplied

Extinction of the Stagnation Point Diffusion Flame: Effect of Conductive Heat Loss into Solid Interior

The flammability boundary of a solid is U-shaped with the ordinate as ambient oxygen percentage (or total pressure) and the abscissa represents flow velocity, gravity level, stretch rate, etc. The left side of the flammability boundary is the quenching branch, and the right side is the blow-off branch. Quenching is due to a weak flame whose heat release from combustion cannot overcome the excessive heat losses from, for example, radiation. Blow-off is due to a short gas residence time that is smaller than the chemical reaction time (traditionally referred to as the small Damköhler number effect). For thick solids, the heat conduction into the solid interior has a long transient after ignition before the solid temperature reaches the steady state. During this slow transient, there is excessive (i.e. more than that at the steady state) gas phase heat conduction into the solid interior. From the viewpoint of the quasi-static gaseous flame, this is an additional heat loss that will affect the gas flame's flammability. In this work, a previously developed one-dimensional axisymmetric stagnation-point diffusion flame model is extended to account for the excessive heat loss into the solid beyond the steady burning state. A non-dimensional excessive conductive heat loss parameter Ψ is defined in the paper to account for the transient solid effect. For each Ψ, the gas phase is considered reaching a quasi-steady state. Extensive numerical computations have been performed to determine the extinction boundary and the characteristics of near-limit flames. The U-shaped flammability boundary is now a series of boundaries with Ψ as a parameter. Flame structures under the same environment condition but with different Ψ are compared. For a given oxygen ambient, the extinction boundary Ψ vs stretch rate is an inverted U-shape. The boundary consists of a blowoff branch and a radiative quenching branch. Extensive computed data including flame and pyrolysis temperatures, burning rate, flame standoff distance, species concentration, and reaction rate are presented along the extinction boundaries. Heat balance analysis is performed for both the gas phase and the gas-solid interface. Along the extinction boundary, the relative importance of conduction, convection, radiation, and chemical reaction are analyzed in detail.

Stagnation point diffusion flame↗
DOE OSTI2020

CO 2 to Bioplastics: Beneficial Re-use of Carbon Emissions from Coal-fired Power Plants using Microalgae

This project sought to address the technical and economic barriers to carbon dioxide (CO 2 ) capture and utilization using microalgae. Specifically, a dual photobioreactor (PBR)/open raceway pond (ORP) cultivation system was evaluated with respect to capital and operational costs, productivity, and culture health, and compared to an ORP-only cultivation system. In the dual (or “hybrid”) system, a cyclic flow PBR was used to provide inoculum for two open raceway ponds after each harvest, with the hypothesis that this PBR + pond system should result in improved performance compared to traditional raceway ponds. Two other raceway ponds were operated conventionally as a control, the experiments being performed at Duke Energy’s East Bend Station in northern KY using coal-derived flue gas as the CO 2 source. Although algae productivity achieved was mediocre due to the fact that the experiments had to be performed late in the growing season when climatic conditions were not optimal for high growth, the hybrid cultivation system showed a higher level of algae productivity (statistically significant) compared to the traditional ponds. In order to realize the maximum value of the algal biomass produced, fractionation of the biomass was examined. Bioplastic compounding and material characterization was subsequently performed using three feed stocks: whole biomass, lipid-extracted biomass and the proteinaceous residue from full fractionation. In general, the whole and lipid-extracted biomass gave similar results in terms of the properties of the resulting bioplastics, while the proteinaceous residue showed promise for the production of a polybutylene adipate terephthalate (PBAT) blend. Finally, sustainability assessment of a putative biorefinery system was conducted, including techno-economic and life cycle impact assessment. Nine different production scenarios were considered, comprising combinations of the three different biomass growth architectures (ORP, PBR and dual PBR-ORP systems) coupled with three different algae biomass processing pathways (drying only, lipid extraction and fractionation). Results show that the minimum selling price of the bioplastic feed stock (BPFS) is within the realm of economic competition with prices as low as $970 USD tonne -1 . Additionally, life cycle impact assessment results indicate drastic improvements in performance of the produced BPFS, with reductions in greenhouse gas emissions ranging between 67 and 116% compared to a petroleum based plastic feedstock.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1998

Selection and Assessment of Predictions of the Mars Pathfinder Landing Site

The successful landing of the Mars Pathfinder spacecraft on Mars allows the review of the process of selecting the landing site and assessing predictions made for the site based on Viking and Earth-based data. Selection of the landing site for Mars Pathfinder was a two-phase process. The first phase took place from October 1993 to June 1994 and involved: initial identification of engineering constraints, definition of environmental conditions at the site for spacecraft design, and evaluation of the scientific potential of different landing sites. This phase culminated with the first "Mars Pathfinder Landing Site Workshop", held at the Lunar and Planetary Institute in Houston, Texas on April 18-19, 1994, in which suggested approaches and landing sites were solicited from the entire scientific community. A preliminary site was selected by the project for design purposes in June 1994. The second phase took place from July 1994 to March 1996 and involved: developing criteria for evaluating site safety using images and remote sensing data, testing of the spacecraft and landing subsystems (with design improvements) to establish quantitative engineering constraints on landing site characteristics, evaluating all potential landing sites on Mars, and certification of the site by the project. This phase included a second open workshop, "Mars Pathfinder Landing Site Workshop II: Characteristics of the Ares Vallis Region and Field Trips in the Channeled Scabland, Washington" held in Spokane and Moses Lake September 24-30, 1995 and formal acceptance of the site by NASA Headquarters. Engineering constraints on Pathfinder landing sites were developed from the initial design of the spacecraft and the entry, descent and landing scenario. The site must be within 5 degrees of the subsolar latitude at the time of landing (15N for maximum solar power and flexible communications with Earth. It also must be below 0 km elevation to enable enough time for the parachute to bring the lander to the proper terminal velocity for landing. The entire landing ellipse, which is 70 km by 200 km due to navigational, ephemeris and atmospheric uncertainties, must be free of steep slopes, scarps and obvious hazards in Viking orbiter images, have acceptable radar reflectivity, moderate rock abundances and have little or no dust. Scientific considerations of the Mars Pathfinder payload and mission indicate that analyses of "grab bag" samples at the mouths of outflow channels can offer a first order assessment of a variety of rock types on Mars. Highland sites offer the advantage of in situ analysis of ancient rocks on Mars that record crustal differentiation and the nature of the early environment. Dark gray sites offer the potential of analyzing unweathered and unoxidized materials. Following a general assessment of the safety of different sites, a preliminary selection of a "grab bag" site was made. This site, Ares Vallis, is near the mouth of an outflow channel that may contain ancient Noachian terrain, Hesperian ridged plains, and reworked channel materials. All potential landing sites on Mars that met basic safety criteria were analyzed in detail. Sites (100 by 200 km target ellipses) were considered safe if they were below 0 km elevation, were free of obvious hazards (high relief surface features) in high-resolution (< 50 m/pixel) Viking orbiter images and had acceptable reflectivity and roughness at radar wavelengths, high thermal inertia, moderate rock abundance, low red to violet ratio, and low albedo. Only 4 sites on Mars met all the above criteria, which included 1995 opposition 3.5 cm delay-Doppler radar data. Complete data were evaluated for 7 sites and the Viking landing sites for comparison for all the above criteria as well as crater abundance, hill and mesa abundance, slopes over meter to kilometer scales, low altitude winds (from global circulation models and slopes), the size-frequency distribution of large rocks, as well as rover trafficability and science potential. Discussion of potential hazards at Ares Vallis using a variety of data sets (including radar) at a second open workshop, indicated this site cannot be shown to be any more hazardous than the Viking landing sites. Field trips to the Channeled Scabland and the Ephrata Fan, analogs for Ares Vallis and the landing site, respectively, provided valuable insight into possible geologic processes and potential surface characteristics. Three sites met all the data requirements and safety criteria for landing Pathfinder. Ares Vallis was selected by the project because it appeared acceptably safe (although it appeared to have greater rock abundances than other sites, its elevation was likely the best known) and offered the prospect of analyzing a variety of rock types expected to be deposited by catastrophic floods, which would enable addressing first-order scientific questions such as differentiation of the crust, the development of weathering products, and the nature of the early martian environment and its subsequent evolution. The selection was reviewed by an external board at a number of meetings and accepted, and the site was approved by NASA Headquarters. Data gathered by the Pathfinder lander' and rover provides the opportunity to test the predictions made for the site in the selection process based on remote observations from Earth, orbit, and the surface. The discussion below is taken from Golombek et al. to which the reader is referred for a more complete discussion and a complete list of references, which are omitted here for brevity. Many characteristics of the landing site are consistent with its being shaped and deposited by the Ares and Tiu catastrophic floods. The rocky surface is consistent a depositional plain comprising semi-rounded pebbles, cobbles and tabular boulders (some of which appear imbricated and/or inclined in the direction of flow) that appear similar to depositional plains in terrestrial catastrophic floods. The Twin Peaks appear to be streamlined hills in lander images, which is consistent with interpretations of larger hills in Viking orbiter images of the region that suggest the lander is on the flank of a broad, gentle ridge trending northeast from Twin Peaks. This ridge, which is the rise to the north of the lander, is aligned in the downstream direction from the Ares and Tiu Valles floods, and may be a debris tail deposited in the wake of the Twin Peaks. Channels visible throughout the scene may be a result of late stage drainage. As predicted by delay-Doppler radar measurements and tracking results, the average elevation of the center of the site was about the same as Viking Lander I relative to the 6.1 mbar geoid. The Doppler tracking and two-way ranging estimate for the elevation of the spacecraft is only 45 in lower than the Viking I Lander and within 100 in of that expected, which is within the uncertainties of the measurements. After landing, surface pressures and winds (5-10 m/s) were found to be similar to expectations based on Viking data, although temperatures were about 10 K warmer. The temperature profile below 50 km was also roughly 20 K warmer. As a result, predicted densities were 5% higher near the surface and up to 40% lower at 50 km but within the entry, descent and landing design margins. The populations of craters and small hills and the slopes of the hills measured in high-resolution (38 m/pixel) Viking orbiter images and the radar derived slopes of the landing site are all consistent with observations of these properties in the lander images. A rocky surface was expected from Viking Infra-Red Thermal Mapper (IRTM) observations and comparisons with the Viking landing sites. The observed cumulative fraction of area covered by rocks with diameters greater than 3 cm and heights greater than 0.5 in (potentially hazardous to landing) at Ares is similar to that predicted by IRTM observations and models of Viking lander and Earth analog rock size-frequency distributions. The IRTM prediction postulated an effective thermal inertia of 30 (10(exp -3) cgs units - cal/cubic cm/s(exp 0.5)/K) for the rock population, but we obtain a slightly different effective thermal inertia for the actual rock population. The validity of interpretations of radar echoes prior to landing are supported by a simple radar echo model, an estimate of the reflectivity of the soil from its bulk density, and the fraction of area covered by rocks. In the calculations, the soil produces the quasi-specular echo and the rocks produce the diffuse echo. The derived quasispecular cross section is comparable to the cross-sections and reflectivities reported for 3.5-cm wavelength observations. The model yields a diffuse echo that is modestly larger than the polarized diffuse echo reported for 3.5-cm wavelength observations. At 12.5-cm wavelength, similar rock populations at Ares and the Viking I site were expected because the diffuse echoes are comparable, but the large normal reflectivities suggests that bulk densities of the soils at depth are greater than those at the surface. We also obtain a fine-component inertia near 8.4 which agrees with the fine-component inertia of 8.7 (in 10(exp -3) cgs units) estimated from thermal observations from orbit by the IRTM; for this estimate, we used a bulk thermal inertia of 10.4 for the landing site, an effective thermal inertia near 40 (10(exp -3) cgs units) for the rock population, and a graphical representation of Kieffer's model. Color and albedo data for Ares suggested surfaces of materials at Ares Vallis would be relatively dust free or unweathered prior to landing compared with the materials at the Viking landing sites. This suggestion is supported by the abundance of relatively dark-gray rocks at Ares and their relative rarity at the Viking landing sites, where rocks are commonly coated with bright red dust. Finally, the 40 km long Ephrata Fan of the Channeled Scabland in Washington state, which was deposited where c

Golombek, M. P.↗
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↗
DOE OSTI2026

HDG-1 Experiment Irradiation Monitoring Data Qualification Final Report

SUMMARY The U.S. Department of Energy (DOE) Advanced Reactor Technologies (ART) Graphite Research and Development (GRD) Program is conducting a series of six experiments to quantify the effects of irradiation on nuclear-grade graphite. This report documents the qualification of irradiation monitoring data for the fifth experiment, High Dose Graphite-1 (HDG-1). Qualified monitoring data are required by the ART program to support the design and licensing of the first high-temperature reactor (HTR) nuclear plant. Data are classified as Qualified if they meet the usage requirements described in the experiment planning and quality assurance (QA) documents, Failed if they do not meet those requirements and provide no usable information, or Trend if they do not fully meet all requirements but still provide useful information subject to an assessment of how any deficiencies may affect a particular use of the data. HDG-1 irradiation began with Advanced Test Reactor (ATR) Cycle 168B on August 24, 2020, and concluded after Cycle 173C on January 27, 2025. The HDG-1 capsule was removed from the reactor core twice—during core internal change (CIC) Cycle 170A and powered axial locator mechanism (PALM) Cycle 172A—to prevent overheating of the graphite specimens during high-power PALM cycles. The capsule was therefore irradiated during a total of seven normal ATR cycles: 168B, 169A, 171A, 171B, 173A, 173B, and 173C. Irradiation monitoring data evaluated in this report include thermocouple (TC) temperature, gas flow rate, gas moisture, gas pressure, specimen load, and graphite stack displacement. Temperature. A total of 14,508,065 TC temperature records were captured. Of these, 13,901,785 (95.8%) are Qualified and 606,280 (4.2%) are Failed. The principal source of failed temperature data was the instrument failure of TC-9 (Zone 2) on June 24, 2024, and TC-10 (Zone 1) on July 5, 2024, near the end of Cycle 173A, which resulted in 595,554 Failed readings. An additional 379 missing values and 10,347 slightly negative values from TC-13 during ATR outages are also Failed. Neither TC-9 nor TC-10 was used as a temperature-control TC, and their failures did not compromise capsule condition monitoring. Correlation analysis of all 13 TCs found no evidence of virtual junction formation. Control chart analysis revealed clear downward drift of approximately 80°C for TC-6 (Zone 3) relative to other stable TCs, and possible downward drift of approximately 60°C for TC-13 relative to the Zone 5 control TC (TC-1), though TC-13 remained consistent with the Zone 2 control TC (TC-12). Gas flow. A total of 20,088,090 gas flow rate records were captured. Of these, 19,941,463 (99.3%) are Qualified and 146,627 (0.7%) are Failed due to missing values. All argon, helium, and total gas flow data were within expected ranges throughout the irradiation. Gas moisture. A total of 1,116,005 outlet gas moisture values were captured. Of these, 1,101,421 (98.7%) are Qualified and 14,584 (1.3%) are Failed, comprising 14,556 out-of-range values and 28 missing values. The out-of-range moisture values exceeded 22,000 ppmv for approximately 1 week at the beginning of Cycle 173A, when accumulated moisture evaporated after the capsule was retrieved from water storage during PALM Cycle 172A and reinserted into the east flux trap. Moisture levels returned to below 25 ppmv for the remaining three cycles, and the transient high-moisture event did not affect the integrity of specimen irradiation. Gas pressure. A total of 7,812,035 gas pressure values were captured. Of these, 6,642,048 (85.0%) are Qualified and 1,169,987 (15.0%) outlet pressure values are Failed, comprising 718,537 zero outlet pressure values due to sensor failure from Cycle 168B through Cycle 171B, 54,550 missing values, and 396,900 too-low outlet pressure values, ranging from 1.1 to 1.6 psia after sensor replacement during Cycle 173A. Load. A total of 6,696,030 load values were captured. Of these, 6,694,580 (99.98%) are Qualified and 1,450 (0.02%) are Failed due to missing values. Applied loads to the six specimen stacks were stable throughout the irradiation. Stack displacement. A total of 6,696,030 displacement values were captured. Of these, 5,713,297 (85.32%) are Qualified and 3,781 (0.06%) are Failed due to missing values. Stack displacement increased consistently throughout the irradiation, reaching approximately 3.08 in. for Channels 5 and 6 by the end of irradiation. 978,952 (14.62%) substantially elevated displacements observed for Channel 6 beginning in Cycle 171A and for Channel 5 beginning in Cycle 173A are assigned Trend status. Raising pressure. A total of 1,115,999 raising pressure values were captured. Of these, 1,115,430 (99.95%) are Qualified and 569 (0.05%) are Failed due to missing values. Ram pressure. A total of 6,696,030 ram pressure values were captured. Of these, 6,692,249 (99.95%) are Qualified and 3,484 (0.05%) are Failed due to missing values. Stack raising was perf

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Coal to Carbon Fiber (C2CF) Continuous Processing for High Value Composites (Final Report)

Coal tar is a condensed and recovered by-product of the coking of metallurgical coal for steel production. The heaviest fraction of distilled coal tar is an isotropic pitch largely used as a binder in the manufacturing of carbonaceous electrodes for primary aluminum smelting and in electric arc furnaces. Coal tar pitch offers high carbon yield upon carbonization. In this project, a process to convert the domestically sourced isotropic coal tar pitch, containing very low particulates (QI = 0.32 wt.%), to form flow-domain mesophase pitch amenable for melt spinning into precursor (green) fibers for carbon fiber, was developed and optimized. The final reproducible processing developed is reviewed in this report, along with several characterizations of the mesophase pitch. A final definition of the characteristics of a ‘spinnable’ mesophase pitch is presented. With this mesophase pitch, reproducible and stable multifilament melt spinning was developed, producing green fiber tows with filament diameters of approximately 20 m. The multifilament melt spinning was the most challenging aspect of the project and required the most effort. The best practices learned from this project for melt spinning are reviewed in this report. Once spun, the green fibers were oxidatively stabilized, carbonized and graphitized under inert gas atmosphere to form the final carbon fibers. Given the relatively high softening point of the mesophase pitch, no issues of interfilament fusion were observed during batch oxidation, and subsequent batch carbonization and graphitization went smoothly in all cases. An ~ 80 wt.% conversion of the green fiber mass to final carbon fiber was achieved. After graphitization, the carbon fibers showed high tensile moduli (most were ~ 600 GPa, or 87 Msi) consistent with commercially-available high-performance pitch-based carbon fiber. However, tensile strength and strain to failure were comparatively low. Further work to reduce defects in and on the fiber surfaces would increase these properties. SEM imaging of the graphitic fiber textures is presented herein. Rudimentary composites were fabricated from the carbon fibers and characterized showing similar modulus to baseline composites fabricated with commercial carbon fiber. Finally, a basic economic analysis was done showing the potential to increase the value of the isotropic coal tar pitch by up to 13.6 to 136 times based on a carbon fiber value of $\$$5/lb to $\$$50/lb, respectively. Moreover, the site case study suggests that the coal tar from the single integrated steel mill could supply production of up to 16 kt/yr of carbon fiber. Finally, a technological gap analysis was done which shed light on remaining technical challenges. These challenges included: recovery and utilization of condensates from the mesophase pitch processing, further advancing and increasing stability of the multifilament melt spinning processing, optimization of the oxidation processing, defect reduction for increased carbon fiber strength, and the development of a weaving process towards carbon fiber fabrics. To maximize the coal value chain, the primary objectives of this project were to (a) develop and scale efficient processing technology for producing melt-spinnable mesophase pitch from isotropic coal tar pitch, (b) clarify and simplify tedious continuous fiber processing technologies (particularly multifilament melt spinning of mesophase pitch) towards the efficient production of high performance carbon fiber, and (c) demonstrate and characterize representative composite parts derived from the final carbon fiber. Immense progress was made on all 3 objectives and is detailed in this report.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Mechanism for Acetone and Crotonaldehyde Production during Steam Reforming of Ethanol over La 0.7 Sr 0.3 MnO 3–x Perovskite: Evidence for a Shared C4 Aldol Addition Intermediate

Here a mechanistic study was conducted on the catalytic conversion of ethanol over La 0.7 Sr 0.3 MnO 3–x perovskite catalysts in the presence and absence of water. The study sought insights into the path of C–C coupling toward acetone and crotonaldehyde and also into clarifying whether the lack of previous reports of C–C coupling over La 0.7 Sr 0.3 MnO 3–x (100) could be due to a “pressure gap”. Several types of experiments were performed at 400–800 K: flow experiments with a torr range reactant gas flown over La 0.7 Sr 0.3 MnO 3–x powders; ultra-high vacuum experiments with continuous gas exposures to a La 0.7 Sr 0.3 MnO 3–x (100) single-crystal sample; and torr range continuous gas exposures to a La 0.7 Sr 0.3 MnO 3–x (100) single-crystal sample. When ethanol and water were flown over La 0.7 Sr 0.3 MnO 3–x powders at 400–800 K, the products detected were ethene, acetaldehyde, acetone, crotonaldehyde, CO, CO 2 , and H 2 . Acetone was catalytically produced over both the La 0.7 Sr 0.3 MnO 3–x powder and the La 0.7 Sr 0.3 MnO 3–x (100) single-crystal sample at temperatures of 700–800 K when reaction conditions were on the order of 1 Torr of reactant gas and with an excess of water relative to ethanol (1 ethanol/9 water). Isotopic labeling with deuterium was used to gain insights into the C–C coupling reaction mechanism and paths in species with three and four carbons (C 3 and C 4 species). Additionally, steady-state isotopic transient kinetic analysis (SSITKA) experiments + simulations using carbon labeling of the ethanol feed were performed. Three mechanistic paths were considered for the C–C coupling step: the first two paths, A and B, involve coupling between two intermediates which are both in oxygen vacancies; and the third path, C, involves coupling between one intermediate in an oxygen vacancy and one intermediate outside of an oxygen vacancy. The results suggest that the dominant path to the C 3 product, acetone, depends on the conditions. The less active path (attributed to path A or B) occurs at 600–700 K and involves coupling between two irreversibly bound species. The more active path (attributed to path C) requires an excess of water, becomes dominant at 600–800 K, and involves coupling between one irreversibly bound species and one reversibly bound species. Based on these various observations from experiments and simulations, an elementary step is proposed for acetone formation involving a previously unreported C 4 transition state that is formed after aldol addition. Density functional theory calculations were performed based on this hypothesis, and it confirmed that this specific and previously unreported aldol addition path to acetone does exist and that this path consistent with the experimental data. In this path, C–C formation occurs to create a C 4 intermediate that is bound to an oxygen vacancy, then a hydrogen transfer with C–C bond breaking occurs that results in the production of the acetone molecule. The proposed mechanism is also consistent with the experimental observation that acetone formation has a greater than first-order dependence on the water vapor pressure.

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

NNSS Soils Monitoring: Plutonium Valley (CAU 366) FY2019

Desert Research Institute (DRI) conducted a field assessment of the potential for radionuclide-contaminated soil to be transported from the Plutonium Valley Dispersion Sites Contamination Area (CA) because of both wind and storm water runoff. This activity supported U.S. Department of Energy (DOE) Environmental Management Nevada Program (EM NV) efforts to establish post-closure monitoring plans for the Plutonium Valley Dispersion Sites Corrective Action Unit (CAU) 366. The DRI task was intended to identify the likely mechanism(s) of transport and determine the meteorological conditions that might cause the movement of radionuclide-contaminated soils. The emphasis of the work was on collecting sediment transported by channelized storm runoff and measuring airborne dust concentrations and associated wind conditions. These data will facilitate an appropriate post-closure monitoring program. In 2011, DRI installed meteorological monitoring stations north and south of the Plutonium Valley CA—known as Pu Valley North and Pu Valley South, respectively—as well as a fluvial sediment sampling station within the CA. Since installation, temperature, wind speed, wind direction, relative humidity, precipitation, solar radiation, barometric pressure, soil temperature, volumetric soil moisture content, and airborne particulate concentrations have been collected at both meteorological stations. The maximum, minimum, and average or total (as appropriate) for each of these parameters were recorded for each 10-minute interval. The sediment sampling station included an automatically activated sampling pump with collection bottles for suspended sediment, which was activated when sufficient flow was present in the channel, and passive traps for bedload material that was transported down the channel during runoff events. This report presents data collected from these stations during fiscal year (FY) 2019. During the FY2019 (October 1, 2018, through September 30, 2019) reporting period, the warmest month was July and the coldest month was February. Monthly total solar radiation was highest in June or July (depending on the station) and lowest in December. Monthly mean wind speeds were highest during February. At the Pu Valley North station, winds were commonly northerly, northwesterly, and southerly. At the Pu Valley South station, the wind direction was more variable and winds frequently came from most directions, excluding westerly and northwesterly winds. Winds above 15.0 miles per hour (mph) (24.1 kilometers per hour [km/hr]) were frequently southerly at both stations, and less-frequent strong winds were northerly (at the Pu Valley South station) or northwesterly (at the Pu Valley North station). Monthly average relative humidity ranged from 18 percent in August to 62 percent in January and February. Monthly total precipitation ranged from zero at both stations in September to 2.23 inches (in) (56.6 millimeters [mm]) during February at the Pu Valley South station and 2.09 in (53.1 mm) during March at the Pu Valley North station. From October 1, 2018, through September 30, 2019, the total precipitation was 8.23 in (209 mm) and 8.15 in (207 mm) at the Pu Valley South and Pu Valley North stations, respectively. The largest daily precipitation totals of 1.46 in (37.1 mm) and 1.14 in (29.0 mm) occurred at the Pu Valley South and Pu Valley North stations, respectively, on March 6, 2019. However, the ISCO autosampler was not activated during any storm event in FY2019 because the channel water depth measurements were insufficient for suspended sediment sample collection. Additionally, because no significant flow through the channel was indicated, no bedload samples were collected. Light breezes of 0.0 mph to 5.0 mph (0.0 km/hr to 8.0 km/hr) occurred most frequently (approximately 55 percent to 57 percent of the time). The frequency of occurrence diminished exponentially as the wind speed increased such that winds in excess of 25.0 mph (40.2 km/hr) occurred less than 0.01 percent of the time. Winds in excess of 15.0 mph (24.1 km/hr) were most commonly southerly, which is likely controlled by the topography of the valley because topographic highs define the east and west sides of the valley and converge toward the north.The concentrations of PM2.5 (particulate matter with an aerodynamic diameter ≤2.5 micrometers [μm]) and PM10 (particulate matter with an aerodynamic diameter ≤10 μm) in the air generally increase as wind speed increases. However, the particle profiler at the Pu Valley South station regularly malfunctioned, which resulted in poor data quality and a limited period of record within FY2019. Therefore, discussion of the annual dust observations is mostly limited to the Pu Valley North station. At the Pu Valley North station, significant increases in windblown dust concentrations were observed when wind speeds exceeded 15.0 mph (24.1 km/hr), especially during southerly winds. When all wind speeds were considered, high dust concentrations at the Pu Valley North station were frequently associated with northerly, northwesterly, or southerly winds. The majority of PM10 transport at the Pu Valley North station occurred when winds were blowing northerly and northwesterly because of the high frequency of winds up to 15.0 mph (24.1 km/hr). The ratio of PM10 to PM2.5 is a qualitative indicator of the proximity of dust sources to the observation point and the values of the PM10 to PM2.5 ratio tend to be higher nearer to the source. The PM10 to PM2.5 ratio increased from 1.2 to 2.8 with higher wind speeds at the Pu Valley North station. This increase suggests that some locally sourced PM10 was present in the air when local winds exceeded 15.0 mph (24.1 km/hr). Local dust suspension is interpreted during major dust events when peaks in the PM10 to PM2.5 ratio occur simultaneously with the maximum wind speed and the peak PM10 concentration. Analysis of the five major dust events in FY2019 revealed that northerly winds were dominant in three events and southerly winds were dominant in two events. Local suspension of dust within Plutonium Valley is suspected to have occurred during two of the five major events (April 9, 2019, and September 28 to 29, 2019), and a nonlocal source of dust is presumed for the three other major dust events (October 14, 2018; May 16, 2019; and September 2, 2019). The three suspected regional dust transport events exhibited winds that typically did not exceed the 15.0 mph (24.1 km/hr) wind speed threshold for local dust transport and featured PM10 to PM2.5 ratios that indicate large-scale regional dust transport. Similar event characteristics (e.g., wind speed and PM10 concentration) were observed at air monitoring stations approximately 70 mi (113 km) to the northwest within the Tonopah Test Range (TTR). Therefore, these three regional events may not reflect transport from the Plutonium Valley CA. Routine monitoring of meteorological and hydrologic parameters at Plutonium Valley has been discontinued at the request of EM NV, effective at the conclusion of FY2019. Environmental data acquisition has been completed, and both Pu Valley North and Pu Valley South meteorological stations within Plutonium Valley as well as the ISCO autosampler station within the CA were removed on October 14, 2019.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Progress Report on Performance of A709 and G91 Steels in Sodium

Specimens in six different processing and heat treatment conditions of A709 H58776 were exposed to sodium at 550, 600, and 650°C, respectively for various exposure times. G91 base metal was tested in sodium at 550, 600, and 650°C, and G91 weldment at 550 and 600°C, respectively. Parallel thermal aging experiments were conducted on G91 and A709 to obtain thermal aging data for comparison with sodium exposure results to separate the thermal and sodium effects. The corrosion data obtained on A709 steel continue to show low corrosion rates, which indicate good compatibility of A709 with sodium when oxygen content is controlled. G91 has also shown acceptable corrosion rates over the temperature range of 550-650°C in sodium. Thermal aging or sodium exposures reduced the tensile strength, uniform elongation, and total elongation for A709 ESR, AOD specimens. These specimens consistently showed higher yield stress and tensile strength and lower uniform and total elongations after sodium exposure than after thermal aging under comparable testing conditions. It implies that there was an additional effect of sodium exposure that gave rise to an increase in tensile strength of A709. Dynamic strain aginginduced flow serrations were completely removed by either thermal aging or sodium exposure at 650°C, which imply that carbon or nitrogen in solution was removed from the solution and formed precipitates during aging or sodium exposures. The tensile data ruled out the possibility of decarburization under the sodium exposure conditions. The HOMO specimens behaved somewhat differently from the ESR and AOD specimens in sodium. While thermal aging and sodium exposures at 550 and 600°C decreased the yield stress and the ultimate tensile strength of G91, there was virtually no additional effect resulting from sodium exposures at these two temperatures. In contrast, sodium exposures at 650°C had a drastic effect on the yield stress and the ultimate tensile strength of G91. It is suggested that G91 experienced decarburization in the 650°C sodium environments. Carbon concentrations in sodium in the SMT-1 and SMT-2 loops were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8-1.2 ppm in the SMT-1 loop and 0.3-0.7 ppm in the SMT-2 loop. The carbon activity in sodium in the SMT-2 loop was estimated to be 0.03-0.08 at 600°C, and 0.08-0.2 at 550°C. The carbon activity in sodium at 650°C in the SMT-1 loop was 0.04-0.07. Equilibrium simulation of the carburization – decarburization processes was conducted for A709 and G91 steels exposed in sodium environments at temperatures of 550-700°C. The carbon activity-concentration relationship for G91 was re-evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC carbides dominated the process at low carbon activities. Formation of NbC increases the decarburization resistance of G91 steel. It remains to be understood whether the decarburization resistance provided by MC carbides in G91 can be maintained during long-term operations of SFRs. The carbon concentration-activity relationship for A709 was calculated based on the equilibrium of fcc-austenite and M 23 C 6 carbide phase. The calculations showed that A709 decarburizes at 650°C in the SMT-1 loop environments, which is different from the experimental findings. The carbon concentration-activity relationship in A709 was further evaluated by considering four phases, including fcc-austenite, M 23 C 6 , NbC, and TiC carbide phases. TiC is the most stable carbides among the three carbide phases. The carburization-decarburization process in A709 is dictated by M 23 C 6 carbides at high carbon activities and by TiC at low carbon activities. Formation of TiC increases the decarburization resistance of A709 steel. Because of the complexity of the precipitation process in A709 during sodium exposures, detailed characterization of precipitates in sodium-exposed specimens and improved thermodynamic models are needed to understand the carburization-decarburization behavior of A709. Kinetic analysis of carbon transfer will be investigated for G91 and A709 in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗