SEARCH · Engineering Papers

Results for “Steady state process”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,351 records · Page 75

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Universal Behavior of a Cyclic Oxidation Model

A mathematical model has been generated to represent the iterative, discrete growth and spallation processes associated with cyclic oxidation. Parabolic growth kinetics (k(sub p)) over and a constant spall area (F(sub A)) were assumed, with spalling occurring interfacially at the thickest regions of the scale. Although most models require numerical techniques, the regularity and simplicity of this progression permitted an approximation by algebraic expressions. Normalization could now be performed to reflect all parametric effects, and a universal cyclic oxidation response was generated: W(sub u) = 1/2 {3J(sub u)(sup 1/2)+ J(sub u)(sup 3/2)} where W, is weight change normalized by the maximum and J(sub u) is the cycle number normalized by the number to reach maximum. Similarly, the total amount of metal consumed was represented by a single normalized curve. The factor [(S(sub c)-l)(raised dot)sqrt(F(sub A)k(sub p)DELTAt)] was identified as a general figure of merit, where S(sub c) is the mass ratio of oxide to oxygen and DELTAt is the cycle duration. A cyclic oxidation failure map was constructed, in normalized k(sub p)-F(sub A) space, as defined by the locus of points corresponding to a critical amount of metal consumption in a given time. All three constructions describe behavior for every value of growth rate, spall fraction, and cycle duration by means of single curves, but with two branches corresponding to the times before and after steady state is achieved.

Smialek, James L.↗

NASA NTRS2006

Real-Time Diagnosis of Faults Using a Bank of Kalman Filters

A new robust method of automated real-time diagnosis of faults in an aircraft engine or a similar complex system involves the use of a bank of Kalman filters. In order to be highly reliable, a diagnostic system must be designed to account for the numerous failure conditions that an aircraft engine may encounter in operation. The method achieves this objective though the utilization of multiple Kalman filters, each of which is uniquely designed based on a specific failure hypothesis. A fault-detection-and-isolation (FDI) system, developed based on this method, is able to isolate faults in sensors and actuators while detecting component faults (abrupt degradation in engine component performance). By affording a capability for real-time identification of minor faults before they grow into major ones, the method promises to enhance safety and reduce operating costs. The robustness of this method is further enhanced by incorporating information regarding the aging condition of an engine. In general, real-time fault diagnostic methods use the nominal performance of a "healthy" new engine as a reference condition in the diagnostic process. Such an approach does not account for gradual changes in performance associated with aging of an otherwise healthy engine. By incorporating information on gradual, aging-related changes, the new method makes it possible to retain at least some of the sensitivity and accuracy needed to detect incipient faults while preventing false alarms that could result from erroneous interpretation of symptoms of aging as symptoms of failures. The figure schematically depicts an FDI system according to the new method. The FDI system is integrated with an engine, from which it accepts two sets of input signals: sensor readings and actuator commands. Two main parts of the FDI system are a bank of Kalman filters and a subsystem that implements FDI decision rules. Each Kalman filter is designed to detect a specific sensor or actuator fault. When a sensor or actuator fault occurs, large estimation errors are generated by all filters except the one using the correct hypothesis. By monitoring the residual output of each filter, the specific fault that has occurred can be detected and isolated on the basis of the decision rules. A set of parameters that indicate the performance of the engine components is estimated by the "correct" Kalman filter for use in detecting component faults. To reduce the loss of diagnostic accuracy and sensitivity in the face of aging, the FDI system accepts information from a steady-state-condition-monitoring system. This information is used to update the Kalman filters and a data bank of trim values representative of the current aging condition.

Kobayashi, Takahisa↗

NASA NTRS1997

Candle Flames in Microgravity

This work is a study of a candle flame in a microgravity environment. The purpose of the work is to determine if a steady (or quasi-steady) flame can exist in a microgravity environment, study the characteristics of the steady flame, investigate the pre-extinction flame oscillations observed in a previous experiment in more detail, and finally, determine the nature of the interactions between two closely spaced candle flames. The candle flame is used as a model combustion system, in that in microgravity it is one of the only examples of a non-propagating, steady-state, pure diffusion flame. Others have used the candle to study a number of combustion phenomena including flame flicker, flame oscillations, electric field effects and enhanced and reduced gravitational effects in flames. The present work is a continuation of a small-scale Shuttle experiment on candle flames. That study showed that the candle flame lifetimes were on the order of 40 seconds, the flames were dim blue after a transient ignition period, and that just prior to extinction the flames oscillated spontaneously for about five seconds at a frequency of 1 Hz. The authors postulated that the gas phase in the immediate vicinity of the flame was quasi-steady. Further away from the flame, however, the assertion of a quasi-steady flame was less certain, thus the authors did not prove that a steady flame could exist. They also speculated that the short lifetime of the candle flame was due to the presence of the small, weakly perforated box that surrounded the candle. The Candle Flames in Microgravity (CFM) experiment, with revised hardware, was recently flown aboard the Mir orbiting station, and conducted inside the glovebox facility by Dr. Shannon Lucid. In addition to the purposes described above, the experiments were NASA's first ability to ascertain the merits of the Mir environment for combustion science studies. In this article, we present the results of that experiment. We are also in the process of developing a numerical model of the microgravity candle flame. The status and results of the modeling efforts to date are also presented.

Dietrich, Daniel L.↗

DOE OSTI2021

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ~2 km, as opposed to ~0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ~2 km (arguably 4 km) when the connectivity is low or moderate and ~1 km (arguably 2 km) when the connectivity is strong.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1995

Characteristics of ionospheric convection and field-aligned current in the dayside cusp region

The assimilative mapping of ionospheric electrodynamics (AMIE) technique has been used to estimate global distributions of high-latitude ionospheric convection and field-aligned current by combining data obtained nearly simultaneously both from ground and from space. Therefore, unlike the statistical patterns, the 'snapshot' distributions derived by AMIE allow us to examine in more detail the distinctions between field-aligned current systems associated with separate magnetospheric processes, especially in the dayside cusp region. By comparing the field-aligned current and ionospheric convection patterns with the corresponding spectrograms of precipitating particles, the following signatures have been identified: (1) For the three cases studied, which all had an IMF with negative y and z components, the cusp precipitation was encountered by the DMSP satellites in the postnoon sector in the northern hemisphere and in the prenoon sector in the southern hemisphere. The equatorward part of the cusp in both hemispheres is in the sunward flow region and marks the beginning of the flow rotation from sunward to antisunward. (2) The pair of field-aligned currents near local noon, i.e., the cusp/mantle currents, are coincident with the cusp or mantle particle precipitation. In distinction, the field-aligned currents on the dawnside and duskside, i.e., the normal region 1 currents, are usually associated with the plasma sheet particle precipitation. Thus the cusp/mantle currents are generated on open field lines and the region 1 currents mainly on closed field lines. (3) Topologically, the cusp/mantle currents appear as an expansion of the region 1 currents from the dawnside and duskside and they overlap near local noon. When B(sub y) is negative, in the northern hemisphere the downward field-aligned current is located poleward of the upward current; whereas in the southern hemisphere the upward current is located poleward of the downward current. (4) Under the assumption of quasi-steady state reconnection, the location of the separatrix in the ionosphere is estimated and the reconnection velocity is calculated to be between 400 and 550 m/s. The dayside separatrix lies equatorward of the dayside convection throat in the two cases examined.

Lu, G.↗

DOE OSTI2022

Cryptic Methane-Cycling by Methanogens During Multi-Year Incubation of Estuarine Sediment

As marine sediments are buried, microbial communities transition from sulfate-reduction to methane-production after sulfate is depleted. When this biogenic methane diffuses into the overlying sulfate-rich sediments, it forms a sulfate-methane transition zone (SMTZ) because sulfate reducers deplete hydrogen concentrations and make hydrogenotrophic methanogenesis exergonic in the reverse direction, a process called the anaerobic oxidation of methane (AOM). Microbial participation in these processes is often inferred from geochemistry, genes, and gene expression changes with sediment depth, using sedimentation rates to convert depth to time. Less is known about how natural sediments transition through these geochemical states transition in real-time. We examined 16S rRNA gene amplicon libraries and metatranscriptomes in microcosms of anoxic sediment from the White Oak River estuary, NC, with three destructively sampled replicates with methane added (586-day incubations) and three re-sampled un-amended replicates (895-day incubations). Sulfate dropped to a low value (∼0.3 mM) on similar days for both experiments (312 and 320 days, respectively), followed by a peak in hydrogen, intermittent increases in methane-cycling archaea starting on days 375 and 362 (mostly Methanolinea spp. and Methanosaeta spp., and Methanococcoides sp. ANME-3), and a methane peak 1 month later. However, methane δ 13 C values only show net methanogenesis 6 months after methane-cycling archaea increase and 4 months after the methane peak, when sulfate is consistently below 0.1 mM and hydrogen increases to a stable 0.61 ± 0.13 nM (days 553–586, n = 9). Sulfate-reducing bacteria (mostly Desulfatiglans spp. and Desulfosarcina sp. SEEP-SRB1) increase in relative abundance only during this period of net methane production, suggesting syntrophy with methanogens in the absence of sulfate. The transition from sulfate reduction to methane production in marine sediments occurs through a prolonged period of methane-cycling by methanogens at low sulfate concentrations, and steady growth of sulfate reducers along with methanogens after sulfate is depleted.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2000

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗

NASA NTRS1997

Eye Movement Latencies to Direction Change for Different Classes of Motion

In the analysis of visual motion, local features such as orientation are analyzed early in the cortical processing stream (V1), while integration across orientation and space is thought to occur in higher cortical areas such as MT, MST, etc. If all areas provide inputs to eye movement control centers, we would expect that local properties would drive eye movements with relatively short latencies, while global properties would require longer latencies. When such latencies are observed, they can provide information about when (and where?) various stimulus properties are analyzed. To this end, a stimulus was employed in which local and global properties determining perceived direction-of-motion could be manipulated independently: an elliptical Gabor patch with a drifting carrier, with variable orientation of the carrier grating and the contrast window. We have previously demonstrated that the directional percepts evoked by this stimulus vary between the "grating direction" (the normal to the grating's orientation) and the "window direction" (ARVO 91, 94), and that similar effects can be observed in reflexive eye movements (ARVO 95). Subjects viewed such a stimulus while attempting to maintain steady fixation on the center of the pattern, and the small reflexive eye movements ("stare OKN") were recorded. In the middle of the trial, the orientation of either the grating or the window was rotated smoothly by 30 degrees. Responses to the shift of both grating orientation and window orientation are seen in the average OKN slow phase velocity. Grating rotations produce a rapid OKN rotation to the grating direction (100 ms latency, 300 ms time constant), followed by a slower rebound to the steady state perceived direction midway between the grating and window directions. Window rotations, on the other hand, evoke a slower response (200 ms latency, 500 ms time constant). The results demonstrate multiple cortical inputs to eye movement control: a taste early input driven by orientation, and a slower input from higher areas sensitive to global stimulus properties.

Mulligan, Jeffrey B.↗

NASA NTRS2017

Biosignatures in the Context of Low Energy Flux

Many of the features that are thought of as biosignatures - including the mediation of chemical and physical processes with speed, specificity, and selectivity - result directly or indirectly from life's unique capability to mediate and direct energy flux. As such, it is important to consider the impact that differences in energy flux may have on the quantity and quality of evidence for life. Earth differs from every other body in our solar system in the magnitude of biologically-usable energy flux into a liquid water environment. On a global basis, the capture of light energy into photosynthesis and the flux of chemical energy represented in the products of that photosynthesis (organic material + O2) are about six and four orders of magnitude larger, respectively, than the flux of energy represented in geochemical sources. Our conception of what an inhabited world "looks like" and our intuition about how to search for life are based in this high-energy context. Energy fluxes on worlds beyond Earth may be better approximated by the million-fold smaller flux provided to Earth's biosphere by geochemical sources. As a result, the nature, abundance, and quality of evidence for life that could be expected on an inhabited extraterrestrial world within our solar system may differ profoundly from that found on Earth. Understanding this potential difference in quantitative terms provides important context for the formulation of life detection strategies. The influence of energy flux on biosignatures can be evaluated through reference to the two basic purposes into which life partitions energy flux: (1) Life expends energy to sustain existing biomass in a metabolic steady state (metabolically functional but non-growing). The formal representation of this relationship in the traditional microbiology literature equates biomass directly with energy flux. The direct implication is that worlds having lower energy flux will have correspondingly lower potential to support biomass. Life detection strategies that directly target extant organisms should therefore be prepared to encounter average biomass densities that may be many orders of magnitude smaller than those found in most of Earth's surface environments (2) Life expends energy to synthesize new biomass. An end-member case in which new biomass is created at the energy-limited rate and the corresponding cells are immediately destroyed (so that the energy partitioned to cell maintenance is minimized) establishes an upper bound on the rate at which biological material can enter a bulk global pool. For a specified bulk concentration [i] of any particular biological compound, i, or for biologically produced matter overall, this synthesis rate, R (sub i), defines a characteristic time scale tau (sub i) equals [i] divided by R (sub i). tau (sub i) can be thought of as (a) the minimum time required for biosynthesis to yield a specific bulk concentration (e.g., a detection threshold) of i, and (b) the average residence time of i within a bulk pool when [i] is held in steady state through a balance between biosynthesis and attrition by physical, chemical, or biological consumption. tau (sub i) becomes an important quantity in considering the potential utility of enantiomeric excess (as a product of homochiral biosynthesis) as a biosignature. Spontaneous racemization of amino acids acts to "erase" the signature of homochiral synthesis over time scales that may range from hundreds to hundreds of thousands of years, depending on temperature. For environments in which low energy flux translates to low rates of biosynthesis, including the synthesis of homochiral amino acids, amino acid residence times in pools having detectable concentrations may compare to or significantly exceed the time scale for racemization. This and similar consequences of long residence times should be considered in the formulation of life detection strategies based on detection of biologically-produced species. Fluxes of biologically-useful energy on potentially habitable worlds within our solar system are, at present, not well constrained. Improving such constraint has the potential to inform priorities in the formulation and targeting of search-for-life strategies, based on the implications of energy flux for the abundance and quality of biosignatures overall, and in specific categories.

ENERGY FLUX↗

DOE OSTI2026

Modeling and analysis of synthetic liquid fuel production from CO 2 and nuclear energy using methanol-to-diesel process

Electrofuels (e-fuels) are synthetic fuels produced from carbon dioxide (CO 2 ) and electricity for blending with or replacing petroleum fuels. Nuclear energy is an attractive energy feedstock for e-fuel production because of its low environmental footprint and its ability to provide steady heat and power essential for e-fuels production. We modeled and evaluated the cost and environmental footprint of e-fuels production in the distillate range for three nuclear power scales, 100, 500, and 1000 MWe, through methanol and olefins intermediates leveraging commercial or high technology readiness level (TRL) processes. Compared to the commonly studied e-fuels from Fischer Tropsch process that has a distillate yield of <70% with the rest being low value naphtha, the proposed process via methanol intermediate increases the product selectivity with distillate yield of 96% and only 4% naphtha. The modeled process has a carbon conversion ratio of 98%, and a process energy efficiency of 56% relative to the total equivalent nuclear electricity input. The e-fuel plant economics and GHG emissions were estimated by considering CO 2 collected from ethanol plants adjacent to nuclear power plants. The estimated minimum fuel selling prices (MFSP) of e-fuel is in the range of $5.7-$9.1/gal depending on e-fuel plant scale, electricity cost, and CO 2 transportation distance. The corresponding e-fuels life cycle GHG emissions is estimated in the range of 5-6 gCO 2 e/MJ of liquid fuel using the R&D Greenhouse gases, Regulated Emissions, and Energy use in Technologies (R&D GREET) model.

10 SYNTHETIC FUELS↗

NASA NTRS2014

Landform Evolution Modeling of Sedimentary Processes on Icy Worlds: the Cases of Hyperion and Helene

Sedimentary activity (erosion, transportation and deposition) is a major landscape-shaping process on many icy worlds of the outer solar system for which an understanding of this process is central to their geologic characterizations. Several recent studies, lead by the authors, have identified the effects of sedimentary activity on icy satellite landforms and recognized the role of mass wasting and volatile loss and redistribution in the evolution of these features. We apply state of the art, physics based, landform evolution modeling to icy satellite sedimentary landforms in order to fully quantitatively characterize how their morphologies change with time, and what their present appearances imply for the initial abundance and distribution of loose material, volatiles, and refractories in original pristine landforms. We will present results of our landform evolution modeling to the development of the landscapes of the Saturnian moons of Hyperion and Helene. Our current conclusions for Hyperion are that its unique appearance can be explained in part by the loss to space of ballistic ejecta during impact events. In order to create the smooth surfaces and the reticulate, honeycomb pattern of narrow divides between old craters, appreciable subsequent modification of crater morphology must occur through mass-wasting processes accompanied by sublimation, probably facilitated by the loss of CO2 as a component of the relief-supporting matrix of the bedrock. This mass wasting effectively destroys small craters, at least in part accounting for the paucity of sub-kilometer craters on Hyperion. Helene’s unusual morphology consists of broad depressions (modified large craters) and a generally smooth surface patterned with streaks and grooves. The streaks appear to be oriented down-gradient, as are the grooves. This pattern suggests intensive mass-wasting as a dominant process. Our initial modeling of this surface suggests a Bingham-like rheological behavior for the loose down-slope moving material. Interestingly, as a Bingham flow, the models indicate that aperiodic "intermittent" behavior to be present, suggesting that periods of quiet steady landform evolution are punctuated short periods of active surface readjustment.

Saturnian moons↗

DOE OSTI2025

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Models implemented in the methodological approach to design the initial STEP first wall contour

The official Spherical Tokamak for Energy Production mission aims to demonstrate the ability to generate net electricity from fusion with the STEP Prototype Power plant. One of the key technological and engineering challenges in fusion power plants is managing the loads on the first wall within acceptable limits. Therefore, the conceptual design development of the STEP Prototype Power plant needs to be based on load estimates derived using legitimate plasma physics assumptions through dynamic and flexible tools. The current design foresees the STEP main chamber first wall to withstand steady-state heat loads of up to ~1 MW/m 2 , excluding critical regions expected to receive higher heat loads such as the baffle regions approaching the divertors. These critical areas will require ad hoc assessments and will be designed with the presence of limiters. This article focuses on the models and methodology adopted for designing the 2-D poloidal contour of the STEP first wall, based on the anticipated charged particle and radiation heat loads during normal operation. Firstly, the models adopted for calculating the charged particle and radiation heat loads are introduced. The first model is validated through benchmarking against the particle tracing code SMARDDA, while the second model is verified by comparing it with data from the MAST-U experiment. Secondly, the model used to design the 2-D first wall contour according to the heat loads is explained. We acknowledge that this preliminary design stage assumes certain simplifications, notably an axisymmetric geometry, for computational efficiency and clarity in presentation. It is understood that subsequent design phases will address the complexities of real-world engineering, including non-axisymmetric effects, transient plasma scenarios, and the impact of disruptions on the first wall design. Finally, an automatic procedure based on these models is presented for defining the 2-D poloidal contour of the STEP first wall to minimize heat loads, taking into account the need to radiate most of the alpha-particle and auxiliary heating power. Here, by providing an overview of the models, methodology, and an automatic procedure, this paper contributes to the design process of the STEP first wall, addressing the engineering challenges associated with fusion power plant development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Isoenergetic two-photon excitation enhances solvent-to-solute excited-state proton transfer

Two-photon excitation (TPE) is an attractive means for controlling chemistry in both space and time. Since isoenergetic one- and twophoton excitations (OPE and TPE) in non-centrosymmetric molecules are allowed to reach the same excited state, it is usually assumed that they produce similar excited-state reactivity. We compare the solvent-to-solute excited-state proton transfer of the super photobase FR0-SB following isoenergetic OPE and TPE. We find up to 62% increased reactivity following TPE compared to OPE. From steady-state spectroscopy, we rule out the involvement of different excited states and find that OPE and TPE spectra are identical in non-polar solvents but not in polar ones. We propose that differences in the matrix elements that contribute to the two-photon absorption cross sections lead to the observed enhanced isoenergetic reactivity, consistent with the predictions of our high-level coupledcluster-based computational protocol. We find that polar solvent configurations favor greater dipole moment change between ground and excited states, which enters the probability for TPE as the absolute value squared. This, in turn, causes a difference in the Franck– Condon region reached via TPE compared to OPE. We conclude that a new method has been found for controlling chemical reactivity via the matrix elements that affect two-photon cross sections, which may be of great utility for spatial and temporal precision chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Microemulsions in the driven Widom-Rowlinson lattice gas

An investigation of the two-dimensional Widom-Rowlinson lattice gas under an applied drive uncovered a remarkable nonequilibrium steady state in which uniform stripes (reminiscent of an equilibrium lamellar phase) form perpendicular to the drive direction [R. Dickman and R. K. P. Zia, Phys. Rev. E 97, 062126 (2018)]. Here we study this model at low particle densities in two and three dimensions, where we find a disordered phase with a characteristic length scale (a “microemulsion”) along the drive direction. Further we develop a continuum theory of this disordered phase to derive a coarse-grained field-theoretic action for the nonequilibrium dynamics. The action has the form of two coupled driven diffusive systems with different characteristic velocities, generated by an interplay between the particle repulsion and the drive. We then show how fluctuation corrections in the field theory may generate the characteristic features of the microemulsion phase, including a peak in the static structure factor corresponding to the characteristic length scale. This work lays the foundation for understanding the stripe phenomenon more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2001

Fusion for Space Propulsion and Plasma Liner Driven MTF

The need for fusion propulsion for interplanetary flights is discussed. For a propulsion system, there are three important system attributes: (1) The absolute amount of energy available, (2) the propellant exhaust velocity, and (3) the jet power per unit mass of the propulsion system (specific power). For human exploration and development of the solar system, propellant exhaust velocity in excess of 100 km/s and specific power in excess of 10 kW/kg are required. Chemical combustion cannot meet the requirement in propellant exhaust velocity. Nuclear fission processes typically result in producing energy in the form of heat that needs to be manipulated at temperatures limited by materials to about 2,800 K. Using the energy to heat a low atomic weight propellant cannot overcome the problem. Alternatively the energy can be converted into electricity which is then used to accelerate particles to high exhaust velocity. The necessary power conversion and conditioning equipment, however, increases the mass of the propulsion system for the same jet power by more than two orders of magnitude over chemical system, thus greatly limits the thrust-to-weight ratio attainable. If fusion can be developed, fusion appears to have the best of all worlds in terms of propulsion - it can provide the absolute amount, the propellant exhaust velocity, and the high specific jet power. An intermediate step towards pure fusion propulsion is a bimodal system in which a fission reactor is used to provide some of the energy to drive a fusion propulsion unit. The technical issues related to fusion for space propulsion are discussed. There are similarities as well as differences at the system level between applying fusion to propulsion and to terrestrial electrical power generation. The differences potentially provide a wider window of opportunities for applying fusion to propulsion. For example, pulsed approaches to fusion may be attractive for the propulsion application. This is particularly so in the light of significant development of the enabling pulsed power component technologies that have occurred in the last two decades because of defense and other energy requirements. The extreme states of matter required to produce fusion reactions may be more readily realizable in the pulsed states with less system mass than in steady states. Significant saving in system mass may result in pulsed fusion systems using plasmas in the appropriate density regimes. Magnetized target fusion, which attempts to combine the favorable attributes of magnetic confinement and inertial compression-containment into one single integrated fusion scheme, appears to have benefits that are worth exploring for propulsion application.

Thio, Y.C. Francis↗

DOE OSTI2026

Single-particle chemical analysis reveals organic-rich detonation soot products

The detonation of high explosives produces a wide variety of particulate matter (PM) with distinct properties, not all of which are traditionally studied for chemical composition, formation processes, and forensic applications. We report particle-resolved measurements of Composition B detonation soot using soot particle aerosol mass spectrometry (SP-AMS), identifying carbonaceous species and metals not previously characterized on a single-particle basis. Results are combined with scanning electron microscopy with energy dispersive X-ray spectroscopy (SEM-EDS) to enhance source-dependent signatures. Black carbon, including graphitic carbon and detonation nanodiamonds, contributed 50.5-71.4% of PM 2.5 mass, while 22.5-43.4% was non-refractory organic carbon, a previously overlooked component that exhibited a complex and varying composition. Detonations were performed with and without PMMA confinement and under steady and overdriven conditions. Distinct particles enriched with polycyclic aromatic hydrocarbons (PAHs) were observed in experiments utilizing confinement, with quantities dependent on manufacturing method. SEM-EDS validated SP-AMS findings of metals internally mixed with carbonaceous species and extended the particle size range to 100 µm. This work makes detecting non-refractory organics using single-particle techniques more feasible for detonation forensics and understanding high-energy soot formation. While this analysis uses offline sample collection, SP-AMS could be deployed for in-situ measurements of detonation PM 2.5 transported in the atmosphere.

36 MATERIALS SCIENCE↗