Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “driving behavior”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermal Stability of π-Conjugated n -Ethylene-Glycol-Terminated Quaterthiophene Oligomers: A Computational and Experimental Study

This work represents a joint computational and experimental study on a series of $\textit{n}$-ethylene glycol (PEO$\textit{n}$)-terminated quaterthiophene (4T) oligomers for 1 < $\textit{n}$ < 10 to elucidate their self-assembly behavior into a smectic-like lamellar phase. This study builds on an earlier study for $\textit{n}$ = 4 that showed that our model predictions were consistent with experimental data on the melting behavior and structure of the lamellar phase, with the latter consisting of crystal-like 4T domains and liquid-like PEO4 domains. The present study aims to understand how the length of the terminal PEO$\textit{n}$ chains modulates the disordering temperature of the lamellar phase and hence the relative stability of the ordered structure. In this work, a simplified bilayer model, where the 4T domains are not explicitly described, is put forward to efficiently estimate the disordering effect of the PEO domains with increasing n; this method is first validated by correctly predicting that layers of alkyl (PE)-capped 4T oligomers (for 1 < $\textit{n}$ < 10) stay ordered at room temperature. Both 4T-domain implicit and explicit model simulations reveal that the order-disorder temperature decreases with the length of the PEO capping chains, as the associated increase in conformational entropy drives a tendency toward disorder that overtakes the cohesive energy, keeping the ordered packing of the 4T domains.

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

On the dissolution of a borosilicate glass with the use of isotopic tracing – Insights into the mechanism for the long-term dissolution rate

An understanding of the mechanisms responsible for controlling the long-term corrosion rate of nuclear waste glass is paramount if reliable glass dissolution models are to be used to calculate the controlled release of radionuclides from nuclear waste glass under geologic disposal conditions. Additionally, understanding silicate glass dissolution is also important for natural glasses to understand the role from the dissolution of these glasses has on the composition of natural aquatic systems. Two general mechanisms appear to be responsible for the elemental release from glass to the surrounding biosphere - ion exchange (release principally of alkalis) and matrix dissolution (release of the structural components of the glass). Key unknowns related to these mechanisms are the impact of surface layers on the altered glass and how these layers participate in corrosion. To better understand the impact of these layers on glass corrosion, we made two glasses of the same elemental composition, but with selected elements enriched in specific isotopes, e.g. more 10 B than 11 B. We reacted these glasses for 1 y then the solutions from each were exchanged and the chemistry of the solutions and the solids were followed for another 2 y and in some cases 3 y. This allowed us to follow changes in the glasses with respect to these elements and how these elements were transported through the alteration layers with time and, for at least one element, into the glass. Results from these experiments demonstrate that the release behavior of different elements is strongly dependent on their structural role in the glass (e.g. network formers and modifiers), but, more importantly, the role of water transport and subsequent ion exchange has in the long-term dissolution of the glass studied here. In this article, we highlight the behavior of four elements: lithium, sodium, silicon, and boron. Lithium and sodium, network modifying elements, have similar chemistries. The behavior of Li could be tracked in more detail by following both 6 Li and 7 Li; Li penetrates through the gel layer in both directions without hinderance and into and released from the glass deeper than sodium. Silicon, a network forming element, reacts with the silica-rich alteration layer. Boron, a network forming element, does not accumulate in the gel or the pristine glass and has a very sharp elemental profile between the pristine glass and alteration layer. Boron, lithium, and sodium elemental profiles suggest that there is little if any transport control within the alteration layer. The correlation between the Li and Na profiles and that of a water species suggest that the limiting release is a steady state between the diffusion of a water species and matrix dissolution that results in a low release rate for alkali and this provides a steady driving force for the long-term dissolution of glass.

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

Signatures of polarized chiral spin disproportionation in rare earth nickelates

In rare earth nickelates (RENiO 3 ), electron-lattice coupling drives a concurrent metal-to-insulator and bond disproportionation phase transition whose microscopic origin has long been the subject of active debate. Of several proposed mechanisms, here we test the hypothesis that pairs of self-doped ligand holes spatially condense to provide local spin moments that are antiferromagnetically coupled to Ni spins. These singlet-like states provide a basis for long-range bond and spiral spin order. Using magnetic resonant X-ray scattering on NdNiO 3 thin films, we observe the chiral nature of the spin disproportionated state, with spin spirals propagating along the crystallographic (101) ortho direction. These spin spirals are found to preferentially couple to X-ray helicity, establishing the presence of a hitherto-unobserved macroscopic chirality. The presence of this chiral magnetic configuration suggests a potential multiferroic coupling between the noncollinear magnetic arrangement and improper ferroelectric behavior as observed in prior studies on NdNiO 3 (101) ortho films and RENiO 3 single crystals. Experimentally constrained theoretical double-cluster calculations confirm the presence of an energetically stable spin-disproportionated state with Zhang-Rice singlet-like combinations of Ni and ligand moments

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Integrating Intermediate Traits in Phylogenetic Genotype-to-Phenotype Studies

A major goal of research in evolution and genetics is linking genotype to phenotype. This work could be direct, such as determining the genetic basis of a phenotype by leveraging genetic variation or divergence in a developmental, physiological, or behavioral trait. The work could also involve studying the evolutionary phenomena (e.g., reproductive isolation, adaptation, sexual dimorphism, behavior) that reveal an indirect link between genotype and a trait of interest. When the phenotype diverges across evolutionarily distinct lineages, this genotype-to-phenotype problem can be addressed using phylogenetic genotype-to-phenotype (PhyloG2P) mapping, which uses genetic signatures and convergent phenotypes on a phylogeny to infer the genetic bases of traits. The PhyloG2P approach has proven powerful in revealing key genetic changes associated with diverse traits, including the mammalian transition to marine environments and transitions between major mechanisms of photosynthesis. However, there are several intermediate traits layered in between genotype and the phenotype of interest, including but not limited to transcriptional profiles, chromatin states, protein abundances, structures, modifications, metabolites, and physiological parameters. Each intermediate trait is interesting and informative in its own right, but synthesis across data types has great promise for providing a deep, integrated, and predictive understanding of how genotypes drive phenotypic differences and convergence. We argue that an expanded PhyloG2P framework (the PhyloG2P matrix) that explicitly considers intermediate traits, and imputes those that are prohibitive to obtain, will allow a better mechanistic understanding of any trait of interest. Furthermore, this approach provides a proxy for functional validation and mechanistic understanding in organisms where laboratory manipulation is impractical.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Characteristics of Ultralow Aspect Ratio Toroidal Plasmas and Non-Solenoidal Startup and Edge Stability Studies at Near-Unity Aspect Ratio (Final Report)

This grant supported a long-term research and development program to explore the physics and technology of spherical tokamak (ST) plasma behavior at near-unity aspect ratio ( A ). The PEGASUS Toroidal Experiment is a university-scale ( R ) ~ 0.40 m, α ~ .35 m, I p ≤ 0.3 MA, B T ≤ 0.15 T) experiment that accesses A → 1, to provide large values of the normalized plasma current, I N = I p /α B T , or equivalently the toroidal field utilization factor I p / I TF . The research activities focused on: 1) developing the technique of Local Helicity Injection (LHI) as a viable means of initiating and growing an ST plasma without the use of a central solenoid; 2) exploring the equilibrium and stability properties of spherical tokamak plasmas at extremely low aspect ratio; and 3) obtaining new insights in the properties of the enhanced confinement H-mode regime at A ~ 1. The initial (Phase I) facility employed a unique high-stress central solenoid Ohmic for current drive and heating with simple magnet current waveforms. A subsequent major upgrade (Phase II) offered more flexible control of the plasma evolution along with a powerful LHI system for plasma initiation, heating and sustainment. The LHI technique consists of injection of high-density current streams into vacuum toroidal and vertical magnetic fields. These streams become unstable and support magnetic reconnection to allow the system to relax into a tokamak-like magnetic configuration. This program focused on: 1) development of electron current sources that provide sufficiently high injection currents in the challenging edge plasma region; 2) validation of theoretical limits to achievable I p arising from magnetic helicity and energy conservation plus absolute limits arising from relaxation to a minimal energy states; 3) tests of a power balance model to describe the evolution of the plasma current I p (t); and 4) study of magnetic fluctuations and their role in the I p evolution. This culminated in the production of toroidal currents of I p ≥ 0.2 MA with only ~8 kA of injected current. Operation with Ohmic induction alone in Phase I created plasmas with soft limits on the achievable I p and field utilization ( I p / I TF ~ 1) due to the onset of large low-order (2/1 or 3/2) internal resonant tearing modes. With the expanded capabilities in Phase II operation, methods employed to mitigate those modes and achieve I p / I TF > 1 were: 1) I TF ramp-downs; 2) non-solenoidal startup via LHI; and 3) extremely low I TF Ohmic operations assisted by LHI current injector pre-ionization. Extensive use of LHI-driven startup and drive for non-solenoidal operation produced favorable current profiles to stably operate tokamak plasmas with unprecedentedly high I N = I p / aB T ~ 15 and I p / I TF > 2. With strong reconnection-driven anomalous ion heating, this in turn provided access to tokamak plasmas with B t up to 100%, a world record 2 times higher than that previously achieved in advanced or spherical tokamaks. This high B t regime had a large fraction (~50%) of the plasma volume in a region of an absolute minimum- B configuration. The MHD stability limit in this regime was the ideal external kink mode without a stabilizing wall. These results provided the first demonstration of access to tokamak plasmas with very high B t , high i N , high elongation κ, and low internal inductance ℓ i , and as such served as a proof-ofprinciple for access to this long-sought ST regime. The naturally low B T needed for stability at A ≤ 1.2 reduces the L-H mode power threshold, P LH , so that H-mode plasmas were achieved with Ohmic heating alone. The power threshold for the transition to the H-mode regime was found to be 10× the accepted values from the ITER international scaling, confirming initial trends indicated by experiments on MAST and NSTX, and suggesting that the physics behind this transition is not yet fully understood. The density dependence of the transition and the observed insensitivity to whether the plasma is bounded by a limiter or magnetic divertor are consistent with the recent FM 3 model. Multiple unstable toroidal magnetic modes were measured during ELM (Edge Localized Mode) crashes, with generally lower toroidal mode numbers than seen at higher A , which was attributed to the larger peeling mode instability drive at low A . First-ever edge current profile measurements through an ELM crash showed current-carrying filaments were generated, as seen in ELM simulations and in electromagnetic blob transport theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2009

System for Automated Calibration of Vector Modulators

Vector modulators are used to impose baseband modulation on RF signals, but non-ideal behavior limits the overall performance. The non-ideal behavior of the vector modulator is compensated using data collected with the use of an automated test system driven by a LabVIEW program that systematically applies thousands of control-signal values to the device under test and collects RF measurement data. The technology innovation automates several steps in the process. First, an automated test system, using computer controlled digital-to-analog converters (DACs) and a computer-controlled vector network analyzer (VNA) systematically can apply different I and Q signals (which represent the complex number by which the RF signal is multiplied) to the vector modulator under test (VMUT), while measuring the RF performance specifically, gain and phase. The automated test system uses the LabVIEW software to control the test equipment, collect the data, and write it to a file. The input to the Lab - VIEW program is either user-input for systematic variation, or is provided in a file containing specific test values that should be fed to the VMUT. The output file contains both the control signals and the measured data. The second step is to post-process the file to determine the correction functions as needed. The result of the entire process is a tabular representation, which allows translation of a desired I/Q value to the required analog control signals to produce a particular RF behavior. In some applications, corrected performance is needed only for a limited range. If the vector modulator is being used as a phase shifter, there is only a need to correct I and Q values that represent points on a circle, not the entire plane. This innovation has been used to calibrate 2-GHz MMIC (monolithic microwave integrated circuit) vector modulators in the High EIRP Cluster Array project (EIRP is high effective isotropic radiated power). These calibrations were then used to create correction tables to allow the commanding of the phase shift in each of four channels used as a phased array for beam steering of a Ka-band (32-GHz) signal. The system also was the basis of a breadboard electronic beam steering system. In this breadboard, the goal was not to make systematic measurements of the properties of a vector modulator, but to drive the breadboard with a series of test patterns varying in phase and amplitude. This is essentially the same calibration process, but with the difference that the data collection process is oriented toward collecting breadboard performance, rather than the measurement of output from a network analyzer.

Lux, James↗
NASA NTRSDate not supplied

Icing Physics Studies Using the 3D SIDRM Test Article: Aerodynamic and Supercooled Liquid Icing Analysis

In-flight icing is an important safety issue and is a factor that affects aircraft design and performance. Newer regulations are driving a need for improvements in airframe and engine icing simulation capability. Experimental data is required for development of icing physics models and simulation validation. To that end, this paper discusses analyzed data from icing tests conducted with a 3D test article in 2022 at the NASA Icing Research Tunnel. The test article, whose geometry is representative of an inter-compressor duct and strut region of a turbofan engine, has been designed and constructed to study the physics of supercooled water icing and ice crystal icing. The surfaces of the Simulated Inter-compressor Duct Research Model (SIDRM) can be heated to simulate the warm surfaces of the turbofan inter-compressor duct. The test article is instrumented with pressure taps, heaters, heat flux gauges, and thermocouples, while a 3D laser scanner, cameras, and a scale to measure ice mass were utilized to characterize the icing behavior. The aim of these tests was to generate ice accretions on the SIDRM test article under well-characterized supercooled liquid icing and ice crystal icing conditions. This paper discusses measurements related to aerodynamic tests and supercooled liquid icing tests that were conducted. Aerodynamic measurements were analyzed and compared to computational simulations and were found to be in good agreement for the range of airspeeds (50 to 230 knots) and angles of attack (0 to 4°) tested. Various parametric sweeps were conducted during the supercooled liquid icing portion of the test entry (cloud median volumetric diameter ranged from 15 to 90 µm, total air temperature from -3 to -17 °C, angle of attack from 0 to 4°, and accretion time from 5 to 20 min). These sweeps were performed to measure that parameter’s impact on ice accretion size, location, characteristics (such as glaze/rime ice and shedding behavior), and test article surface temperature. Analysis of the test data showed that clouds composed of larger drops, colder air temperatures, smaller angles of attack, and longer spray times were the primary parameters that resulted in accretions with greater ice mass. Test article angle of attack and cloud droplet size influenced the location of ice accretion as these two parameters directly impact collection efficiency. With respect to icing characteristics, total air temperature dictated icing type, and smaller cloud drop size along with warmer air temperatures resulted in greater amounts of ice shedding. Surface temperature increased during ice accretion from the release of latent (fusion) heat, where total air temperature and cloud drop size impacted the amount of surface temperature change. The icing measurements collected during the SIDRM tests will be used to develop and validate 3D computational engine icing tools, such as GlennICE, that predictively assesses the onset and growth of ice. One of the goals of the sponsoring NASA project is to develop simulation models and tools that can assist in the design and certification of engines for flight in icing conditions in a cost‑effective way.

supercooled liquid icing↗
DOE OSTI2020

Dynamic Structure and Phase Behavior of a Block Copolymer Electrolyte under dc Polarization

There is an important consideration when designing lithium battery electrolytes for advanced applications is how the electrolyte facilitates ion transport at fast charge and discharge rates. Large current densities are accompanied by large salt concentration gradients across the electrolyte. Nanostructured composite electrolytes have been proposed to enable the use of high energy density lithium metal anodes, but many questions about the interplay between the electrolyte morphology and the salt concentration gradient that forms under dc polarization remain unanswered. To address these questions, we use an in situ small-angle X-ray scattering technique to examine the nanostructure of a polystyrene- block -poly(ethylene oxide) copolymer electrolyte under dc polarization with spatial and temporal resolution. In the quiescent state, the electrolyte exhibits a lamellar morphology. The passage of ionic current in a lithium symmetric cell leads to the formation of concurrent phases: a disordered morphology near the negative electrode, lamellae in the center of the cell, and coexisting lamellae and gyroid near the positive electrode. The most surprising result of this study was obtained after the applied electric field was turned off: a current-induced gyroid phase grows in volume for 6 h in spite of the absence of an obvious driving force. We show that this reflects the formation of localized pockets of salt-dense electrolyte, termed concentration hotspots, under dc polarization. Our methods may be applied to understand the dynamic structure of composite electrolytes at appreciable current densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Driving forces: Slab subduction and mantle convection

Mantle convection is the mechanism ultimately responsible for most geological activity at Earth's surface. To zeroth order, the lithosphere is the cold outer thermal boundary layer of the convecting mantle. Subduction of cold dense lithosphere provides tha major source of negative buoyancy driving mantle convection and, hence, surface tectonics. There are, however, importnat differences between plate tectonics and the more familiar convecting systems observed in the laboratory. Most important, the temperature dependence of the effective viscosity of mantle rocks makes the thermal boundary layer mechanically strong, leading to nearly rigid plates. This strength stabilizes the cold boundary layer against small amplitude perturbations and allows it to store substantial gravitational potential energy. Paradoxically, through going faults at subduction zones make the lithosphere there locally weak, allowing rapid convergence, unlike what is observed in laboratory experiments using fluids with temperature dependent viscosities. This bimodal strength distribution of the lithosphere distinguishes plate tectonics from simple convection experiments. In addition, Earth has a buoyant, relatively weak layer (the crust) occupying the upper part of the thermal boundary layer. Phase changes lead to extra sources of heat and bouyancy. These phenomena lead to observed richness of behavior of the plate tectonic style of mantle convection.

Hager, Bradford H.↗
NASA NTRS2023

The Hazard Components of Representative Key Risks. The Physical Climate Perspective

The framework of Representative Key Risks (RKRs) has been adopted by the Intergovernmental Panel on Climate Change Working Group II (WGII) to categorize, assess and communicate a wide range of regional and sectoral key risks from climate change. These are risks expected to become severe due to the potentially detrimental convergence of changing climate conditions with the exposure and vulnerability of human and natural systems. Other papers in this special issue treat each of eight RKRs holistically by assessing their current status and future evolution as a result of this convergence. However, in these papers, such assessment cannot always be organized according to a systematic gradation of climatic changes. Often the big-picture evolution of risk has to be extrapolated from either qualitative effects of “low”, “medium” and “high” warming, or limited/focused analysis of the consequences of particular mitigation choices (e.g., benefits of limiting warming to 1.5 or 2C), together with consideration of the socio-economic context and possible adaptation choices. In this study we offer a representation – as systematic as possible given current literature and assessments – of the future evolution of the hazard components of RKRs. We identify the relevant hazards for each RKR, based upon the WGII authors’ assessment, and we report on their current state and expected future changes in magnitude, intensity and/or frequency, linking these changes to Global Warming Levels (GWLs) to the extent possible. We draw on the assessment of changes in climatic impact-drivers relevant to RKRs described in the 6th Assessment Report by Working Group I supplemented when needed by more recent literature. For some of these quantities - like regional trends in oceanic and atmospheric temperature and precipitation, some heat and precipitation extremes, permafrost thaw and Northern Hemisphere snow cover - a strong and quantitative relationship with increasing GWLs has been identified. For others - like frequency and intensity of tropical cyclones and extra-tropical storms, and fire weather - that link can only be described qualitatively. For some processes - like the behavior of ice sheets, or changes in circulation dynamics - large uncertainties about the effects of different GWLs remain, and for a few others - like ocean pH and air pollution - the composition of the scenario of anthropogenic emissions is most relevant, rather than the warming reached. In almost all cases, however, the basic message remains that every small increment in CO2 concentration in the atmosphere and associated warming will bring changes in climate phenomena that will contribute to increasing risk of impacts on human and natural systems, in the absence of compensating changes in these systems’ exposure and vulnerability, and in the absence of effective adaptation. Our picture of the evolution of RKR-relevant climatic impact-drivers complements and enriches the treatment of RKRs in the other papers in at least two ways: by filling in their often only cursory or limited representation of the physical climate aspects driving impacts, and by providing a fuller representation of their future potential evolution, an important component – if never the only one – of the future evolution of risk severity.

Representative Key Risks↗
NASA NTRS2016

Using Tabulated Experimental Data to Drive an Orthotropic Elasto-Plastic Three-Dimensional Model for Impact Analysis

An orthotropic elasto-plastic-damage three-dimensional model with tabulated input has been developed to analyze the impact response of composite materials. The theory has been implemented as MAT 213 into a tailored version of LS-DYNA being developed under a joint effort of the FAA and NASA and has the following features: (a) the theory addresses any composite architecture that can be experimentally characterized as an orthotropic material and includes rate and temperature sensitivities, (b) the formulation is applicable for solid as well as shell element implementations and utilizes input data in a tabulated form directly from processed experimental data, (c) deformation and damage mechanics are both accounted for within the material model, (d) failure criteria are established that are functions of strain and damage parameters, and mesh size dependence is included, and (e) the theory can be efficiently implemented into a commercial code for both sequential and parallel executions. The salient features of the theory as implemented in LS-DYNA are illustrated using a widely used composite - the T800S/3900-2B[P2352W-19] BMS8-276 Rev-H-Unitape fiber/resin unidirectional composite. First, the experimental tests to characterize the deformation, damage and failure parameters in the material behavior are discussed. Second, the MAT213 input model and implementation details are presented with particular attention given to procedures that have been incorporated to ensure that the yield surfaces in the rate and temperature dependent plasticity model are convex. Finally, the paper concludes with a validation test designed to test the stability, accuracy and efficiency of the implemented model.

Polymer Matrix Composites↗
DOE OSTI2026

Pseudo–Jahn–Teller Distortion in a One-Dimensional π-Conjugated Polymer

Structural distortions in low-dimensional π-conjugated systems profoundly influence their electronic properties, but the control of such behavior in laterally extended systems remains challenging. Here, in this study, we demonstrate that a one-dimensional conjugated polymer─ poly -(difluorenoheptalene-ethynylene) (PDFHE)─undergoes a pronounced out-of-plane backbone distortion, equivalent to a spontaneous symmetry breaking (SSB) of its mirror symmetry. We synthesized PDFHE on noble metal surfaces and characterized its structure and electronic states using low-temperature scanning tunneling microscopy (STM). Rather than adopting a planar, high-symmetry conformation, PDFHE relaxes into nonplanar isomers stabilized by a pseudo–Jahn–Teller (PJT) distortion having mirror-odd out-of-plane character. The distortion lowers the total energy and increases the band gap, providing a concise rationale for the observed symmetry breaking. Density functional theory calculations corroborate these findings, providing a microscopic explanation for the SSB. Our results show that even in mechanically robust extended π-systems, subtle electron–lattice coupling can spontaneously drive significant structural rearrangements, even in mechanically robust extended π-systems.

conformation↗
NASA NTRSDate not supplied

The Influence of Viability, Independence, and Self-Governance on Trust and Public Acceptance of Uncrewed Air Vehicle Operations

Trust is expected to be a critical construct that drives successful use of advanced air mobility (AAM) technologies. As yet, though, the role of trust in human-autonomy interaction is underexplored. Kaber (2018) argues that autonomy requires the highest level of three independent dimensions -- viability, independence, and self-governance. The present study examined whether trust varies across the three dimensions of autonomy under varying levels of risk. Participants in the high-risk group read a series of vignettes on a drone that delivers medical supplies over a city where the current study was conducted. Participants in the low-risk group read a series of vignettes on a drone that delivers fast food over a fictitious city. Each vignette described a drone that is either autonomous (i.e., possesses all dimensions) or automated (i.e., one of the dimensions is compromised). Results imply that the three dimensions of autonomy do not equally influence human-technology trust and behavior.

Advanced Air Mobility↗
NASA NTRSDate not supplied

Launching Liquid Slugs Through Tube Networks: Investigating Multiphase Liquid Slug Acceleration Behaviors in Spacecraft Propulsion Systems using 3D Multiphase VOF CFD Simulations and High-Speed Flow-Visualization Validation Experiments

Propellant management in spacecraft propulsion systems is complex. In many current and future crewed and robotic spacecraft propulsion systems, a number of unintended or unavoidable scenarios can occasionally lead to liquid “slugs” of propellant accumulating in unwanted locations and thereafter being rapidly accelerated through complex and arbitrary tubing networks otherwise filled with vapor or gas. To develop better understanding of the multiphase fluid dynamics associated with these “slug launch” scenarios, NASA Marshal Space Flight Center’s (MSFC) ER42 branch – the propulsion fluid dynamics research branch – conducted 3D multiphase transient Volume of Fluid (VOF) CFD simulations using an in-house tool called LOCI-Stream VOF on semi-canonical tube networks with straight sections, bends, and bend complexes. In addition, high-speed flow visualization validation experiments using propellants were conducted on similar tube network geometries at the NASA White Sands Test Facility (WSTF). These simulations and flow-visualization experiments found these scenarios lead to a variety of unique multiphase fluid dynamic effects involving the leading and trailing edges of the slug in both straight sections and curved bends as well as the progressive erosion, acceleration, and eventual “punch-through” of the driving gas through the liquid slug. Understanding these multiphase fluid dynamic effects has led to unique and often counter-intuitive observations which have been of critical relevance to important engineering concerns associated with these “slug launch” scenarios.

Multiphase Fluid Dynamics↗
DOE OSTI2024

FY24 Report on Water NSTF Testing: Lower Tank Inlet Piping Configuration

The following report serves as a summary of the accomplishments and testing results by the Natural convection Shutdown heat removal Test Facility (NSTF) experimental program over the past 12-month period. A major activity included reconfiguration of the chimney piping geometry which altered the discharge position into the tank from the original 50% tank height to a new lower position at a 10.9% tank height. This modification increases the volume of available coolant thus extending long-term operating capacity, however also results in a decreased liquid driving head which has the potential to reduce the natural circulation efficiency and decrease overall performance. Examination of the tradeoffs for these two configurations is an area of interest for RCCS designers and drove planned test activities this year. Ten matrix tests were performed in FY24, totaling 195 hours of heated operations and 12.9 MWh of electrical heating throughout the year, with eight classified as Accepted per NQA-1 and two classified as Trending. Testing prior to the facility reconfiguration examined the facility response to throttling at the tank inlet, demonstrating increased sensitivity of flow instabilities to throttling within the two-phase region when compared to throttling at the lower sensitivity single-phase inlet region. Testing after the chimney reconfiguration began with a baseline test at conditions of 80% inventory fill and prototypic thermal power input of 2.1 MW t . In addition to establishing a reference for nominal system behavior and performance, repeat testing with multiple subsequent runs demonstrated strong repeatability between tests performed at the same conditions in this new configuration. These tests also examined if influences would occur to system behavior with installation of new higher-resolution instrumentation within the upper chimney. In this critical region where boiling and flashing phenomena dominate, the test results provided confidence that the new instrumentation does not uncharacteristically influence the observed behaviors. An inventory parametric series was then initiated to examine system behavior with the new chimney configuration at six varying initial inventory levels, ranging from high 80% to low 20% fill. Two-phase oscillations with similar peak and mean flow rates were observed when comparing to the previous mid tank configuration. Generally, similar system response trends with inventory were also observed, such as two phase oscillations suppressing as inventories were lowered. However, in one absolute inventory comparison at the highest fill of 80%, flow oscillations saw a gradual growth over the 4 hours of two-phase operation in the new lower tank inlet configuration, opposite of the gradual dampening observed in the mid tank inlet configuration. This can be attributed to, in part, a greater hydrostatic head pressure above the two-phase discharge region where the boiling front is developed. Furthermore, the change resulted in greater sensitivity to liquid degassing phenomena during single-phase heating, causing loop instabilities to form which trigged moderate flow degradation during the period approaching saturation and boiling conditions. This behavior had been observed previously under some conditions but was a common occurrence in recent testing with the newer lower tank inlet configuration. Initial observations from these first data sets suggest an overall larger window of stability for the mid tank inlet configuration compared to the lower tank inlet. Lastly, the lowest inventory fill test was repeated over an extended testing window to examine depletion behavior. Natural circulation flow and effective heat removal performance were observed during most of the testing period; only after the tank became fully drained (0% fill) did flow stagnate and violent geysering events occur. This early observation confirms one relative advantage over the mid tank inlet configuration, which stagnated under comparable conditions at ~20% inventory remaining in the tank.

42 ENGINEERING↗
DOE OSTI2021

The Heating and Pulsations of V386 Serpentis after Its 2019 Dwarf Nova Outburst

Following the pulsation spectrum of a white dwarf through the heating and cooling involved in a dwarf nova outburst cycle provides a unique view of the changes to convective driving that take place on timescales of months versus millennia for nonaccreting white dwarfs. In 2019 January the dwarf nova V386 Ser (one of a small number containing an accreting, pulsating white dwarf) underwent a large-amplitude outburst. Hubble Space Telescope ultraviolet spectra were obtained 7 and 13 months after outburst along with optical ground-based photometry during this interval and high-speed photometry at 5.5 and 17 months after outburst. The resulting spectral and pulsational analysis shows a cooling of the white dwarf from 21,020 to 18,750 K (with a gravity $\mathrm{log}(g)=8.1$) between the two UV observations, along with the presence of strong pulsations evident in both UV and optical at a much shorter period after outburst than at quiescence. The pulsation periods consistently lengthened during the year following outburst, in agreement with pulsation theory. However, it remains to be seen if the behavior at longer times past outburst will mimic the unusual nonmonotonic cooling and long periods evident in the similar system GW Lib.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The hazard components of representative key risks. The physical climate perspective

The framework of Representative Key Risks (RKRs) has been adopted by the Intergovernmental Panel on Climate Change Working Group II in order to categorize, assess and communicate a wide range of regional and sectoral key risks, expected to become severe as a consequence of the potentially detrimental convergence of changing climate conditions, exposure and vulnerability of human and natural systems, and responses through (mal)adaptation measures. Other papers in this special issue treat each of eight RKRs by holistically assessing their current status and future evolution as a result of this convergence. However, their treatment cannot always organize such assessment according to a systematic gradation of climatic changes, often times having to extrapolate the big-picture evolution of risk from either qualitative effects of “low” medium” and “high” warming, or limited/focused analysis of the consequences of particular mitigation choices (e.g., benefits of limiting warming to 1.5 or 2C). In this study we offer as systematic a representation as possible – given current literature and assessments – of the future evolution of the hazard components of RKRs. We draw on the assessment of changes in climatic impact-drivers relevant to RKRs described in the 6th Assessment Report by Working Group 1 (WGI) supplemented when needed by more recent literature. We identify the relevant hazards for each RKR, based upon the WGII authors’ assessment, and we report on their current state and expected future changes in magnitude, intensity and/or frequency, linking these changes to Global Warming Levels to the extent possible. For some of these quantities - like regional trends in oceanic and atmospheric temperature and precipitation, some heat and precipitation extremes, permafrost thaw and Northern Hemisphere snow cover - a strong and quantitative relationship with global average temperature change has been assessed. For others - like frequency and intensity of tropical cyclones and extra-tropical storms, and fire weather - the link to increasing global temperatures can only be described qualitatively. For some processes - like the behavior of ice-sheets, or changes in circulation dynamics - large uncertainties about the effects of different levels of warming remain, and for a few others - like ocean pH and air pollution - the composition of the scenario of anthropogenic emissions is most relevant, rather than the warming reached. In almost all cases, however, the basic message remains that every small increment in CO2 concentration in the atmosphere and associated warming will bring changes in climate phenomena that will contribute to increasing risk of impacts on human and natural systems when these changes happen in conjunction with these systems’ conditions of exposure and vulnerability, and in lack of effective adaptation. Our picture of the evolution of RKR-relevant climatic impact-drivers complements and enriches the treatment of RKRs in the other papers in at least two ways: by filling in an often only cursory or limited representation of the physical climate aspects driving impacts, and by providing a fuller representation of their future possible evolution which represents an important component – if never the only one – of the potential future evolution of risk severity.

54 ENVIRONMENTAL SCIENCES↗