Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive Flow”

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.

1,691 records · Page 94

Influence of Microgravity on Bacterial Pathogen Virulence and Immune Cell Function—Relevance for Infectious Disease Risk During Spaceflight

Spaceflight has measurable impacts on astronaut immune profiles as well as the virulence patterns of bacterial pathogens. Data with respect to human immunity indicate diminished T and NK cell function, altered cytokine profiles, persistent inflammation, and latent herpesvirus reactivation. Furthermore, evaluation of International Space Station (ISS) crewmembers gives evidence of compromised immunity, including atypical allergy, infectious disease, and dermatitis. Data with respect to certain human bacterial pathogens suggest modified virulence that may be enhanced. It is therefore critical to examine this interaction of immune dysregulation and increased microbial virulence and whether it might synergistically increase the risk of infectious disease to crew members. The goal of this study is to use modeled microgravity to study the impact of medically significant ISS bacteria that may have altered virulence on the immune response of the host. This study consists of two primary aims to assess this relationship. First, immune cells will be collected from healthy test subjects and cultured in static or in modeled microgravity conditions together with either control pathogens or with microgravity conditioned pathogens that were grown in bioreactors. Second, immune cells will be collected from astronauts before, during, and after flight and co-cultured with the control or the microgravity conditioned bacteria. Three pathogens were selected for this investigation: Pseudomonas aeruginosa , Salmonella enterica serovar Enteritidis and Burkholderia cepacia . Previously, the optimal cell to bacteria ratios that produce the greatest immune cell responses have been derived for these three bacteria. Cellular activation, as determined by the induction of cell surface activation markers and cytokine profiles, will be measured. Interactions between cells and bacteria will be assessed using fluorescent and electron microscopy techniques. This study will provide critical information to help understand how microgravity alters microbial virulence and the associated infectious disease risk to crewmembers during spaceflight missions. Over the past year, the Immunology and Microbiology laboratories at NASA Johnson Space Center have collaborated to process the astronaut subject samples to complete the last year of the flight portion of this study. This included the completion of four astronaut subjects full mission sampling sets. Sample processing included innate and adaptive cell flow cytometry as well as analysis of cytokine concentrations in the supernatant. The ground control segment of the study will take place during FY25 which involves parallel infections run under static and clinostat conditions. The acquired data sets over this 3-year study are now being analyzed to provide a comprehensive set of results and conclusions that will contribute to a more effective risk assessment for astronauts during spaceflight regarding this host-pathogen context.

Immunology↗

Impacts of Bond Type and Grafting Density on the Thermal, Structural, and Transport Behaviors of Nanoparticle Organic Hybrid Materials–Based Electrolytes

Nanoscale Organic Hybrid Materials (NOHMs) consist of polymers tethered to a nanoparticle surface, and NOHMs formed with an ionic bond between the polymer and nanoparticle have been proposed for electrochemical applications. NOHMs exhibit negligible vapor pressure, chemical tunability, oxidative thermal stability, and high ionic conductivity making them attractive in reactive and separation systems. In this study, NOHMs are synthesized by tethering Jeffamine M2070 (HPE) to SiO 2 nanocores via ionic (NOHM-I-HPE) and covalent (NOHM-C-HPE) bonding to investigate the effect of the bond type on the thermal, structural and transport properties of the tethered HPE. In the neat state, NOHM-C-HPE displays the highest thermal stability in a nitrogen atmosphere, while NOHM-I-HPE is the most stable under oxidative conditions. Small-angle neutron scattering (SANS) reveals the presence of multiple types of HPE polymers in aqueous solutions of NOHM-I-HPE (i.e., tethered, interacting, and free), whereas only tethered HPE is observed in NOHM-C-HPE systems. Moreover, the SANS profiles identify clustering of NOHM-C-HPE in aqueous solutions, but not in the corresponding NOHM-I-HPE solutions, suggesting that the free HPE chains stabilize the dispersion of NOHM-I-HPE. Furthermore, the results of this study elucidate how the bond type and grafting density can be used to tune the properties of NOHMs.

25 ENERGY STORAGE↗

Evaluation of Iodine and Water Co-adsorption on Silver-Functionalized Aerogel

The co-adsorption of iodine and water on silver-functionalized silica aerogel was evaluated using research quantities of the aerogel prepared at Pacific Northwest National Laboratory. A macro thermogravimetric analyzer was used to measure weight changes in a thin bed of the aerogel as air streams containing moisture and/or iodine flowed through the bed. The total gas flow in all tests was adjusted to achieve a velocity of 10 m/min through the bed. Initial heating of the aerogel sample to 150°C in dry air (dew point < -60°C) resulted in weight loss on the order of 1 wt %. Because subsequent phases of the tests indicated that the aerogel did not sorb water, the weight loss was speculated to be caused by volatilization of material that may have been used in the manufacture of the aerogel. Follow-on tests are recommended to ascertain if the aerogel sorbs water at low temperatures, such as 25°C—this would help to determine if a manufacturing reagent is volatilized when raising the sample to iodine sorption temperatures (i.e., 150°C). Alternatively, the off-gas released upon heating could be condensed and analyzed. In tests using humidified air streams with dew points between -60°C and +10°C, the aerogel did not adsorb water at test temperatures of either 75°C or 150°C. Because water alone was not adsorbed, tests were performed to ascertain if iodine sorption varied between cases when nearly no water was present in the gas stream and when water was present. In one experiment, the gas stream consisted of dry air (dew point < -60°C) with an iodine concentration of 50 ppmv. In the other experiment, the gas stream consisted of air with a water dew point of 0°C and an iodine concentration of 50 ppmv. The ultimate weight gain was nearly the same in both experiments and was presumed to be due entirely to adsorption of iodine. However, a reduced rate of iodine uptake is indicated when water vapor is present. The reason for this is unclear. In both cases, purging the sample with dry, iodine-free air resulted in some weight loss (~2 wt %). It is recommended that the gas phase chemistry of the iodine–water system be reviewed to ascertain if iodine species that are less reactive with silver than molecular iodine are present; this could explain the reduced sorption rates when water is present. Further evaluation of the desorption is recommended, especially to determine how raising the temperature would affect desorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Concept of Operations for Management by Trajectory

This document describes Management by Trajectory (MBT), a concept for future air traffic management (ATM) in which every flight operates in accordance with a four-dimensional trajectory (4DT) that is negotiated between the airspace user and the Federal Aviation Administration (FAA) to respect the airspace user's goals while complying with National Airspace System (NAS) constraints. In the present-day NAS, the ATM system attempts to predict the trajectory for each flight based on the approved flight plan and scheduled or controlled departure time. However, once the aircraft starts to move, controllers tactically manage the aircraft to implement traffic management restrictions, separate otherwise conflicting aircraft, and address arising NAS constraints. Tactical controller actions are not directly communicated to the automation systems or other stakeholders. Furthermore, the initial trajectory prediction does not anticipate these disruptions or how they will impact the flight. Consequently, and compounded by gaps in required data and models, trajectory predictions are less accurate than possible, which affects Traffic Flow Management (TFM) performance. A cornerstone of the MBT concept is that all air vehicles have, at all times, an assigned 4DT from their current state to their destination. These assigned trajectories consist of trajectory constraints and descriptions. Pilots and air traffic controllers, with the aid of automation, operate the aircraft to comply with the assigned trajectory, unless first negotiating a revision. Equipped aircraft have substantial responsibility for complying with the assigned trajectory without controller intervention. To maximize the operational flexibility available to the airspace user, the assigned trajectory only imposes trajectory constraints as required to achieve the ATM goals of NAS constraint compliance and aircraft separation. Trajectory descriptions are added to the assigned trajectory to ensure sufficient predictability. To further improve trajectory prediction accuracy, airspace users supplement the assigned trajectory by broadcasting intent information and updating it as necessary. Air vehicle intent is a more detailed description of the airspace user's plan for how the flight will fly the assigned trajectory. Air vehicle intent can change freely, without negotiation, as long as it remains in compliance with the assigned trajectory. Aircraft assigned trajectories, air vehicle intent, and predicted trajectories are shared, creating a common view among stakeholders. A NAS Constraint Service gathers and publishes information about all known NAS constraints, enabling airspace users to be informed participants in trajectory negotiation. Trajectory constraints in the assigned trajectory are mapped to NAS constraints to facilitate identifying which aircraft are affected when NAS constraints change. To support efficient trajectory negotiation, all aircraft provide current information about air vehicle capabilities. Assigned trajectories are constructed to satisfy all known NAS constraints, improving trajectory stability and predictability. Uncertainty and disruptions are handled by modifying the assigned trajectory as far in advance as possible. By proactively negotiating changes to the assigned trajectory, rather than relying on controller-selected tactical actions such as vectors to resolve traffic conflicts or implement miles-in-trail restrictions, MBT keeps aircraft on closed trajectories that are fully known to all stakeholders. Since reactive air traffic control actions cannot be predicted in advance, the downstream trajectory cannot be accurately predicted until they happen. Reliable trajectory predictions allow the system to identify needed modifications to trajectories further in advance, where they can be negotiated and communicated as amendments (i.e., additional or altered trajectory constraints) to the assigned trajectory. Decision Support Tools (DSTs) aid controllers in rapidly defining and communicating closed trajectories to the aircraft and support all stakeholders in trajectory negotiation. Anticipated MBT benefit mechanisms include more accurate trajectory predictions, improved ATM performance and robustness to off-nominal conditions, increased flexibility and operational efficiency, reduced impediments to emerging classes of airspace users accessing NAS resources, reduced environmental impacts, and enhanced safety.

ConOps↗

Effects-Based Monitoring of Geomagnetically-Induced Current Using a Convolutional Neural Network

Geomagnetically-induced current (GIC) due to space weather can flow in the power grid causing undesirable effects such as transformer overheating, misoperation of protection devices, and potential blackouts. It is therefore important to monitor GIC in the power grid to improve online situational awareness and decision-making of system operators during a geomagnetic disturbance. To avoid the costly installation of GIC monitors at transformers’ neutrals, it is desirable to find correlations between GIC and already-monitored parameters. Hence, this work proposed the use of a convolutional neural network (CNN) to compute GIC amplitudes from learned patterns in the time-series data of transformer even harmonic currents. Using an electromagnetic transient program, GIC injection simulations were performed for a modeled Dominion Energy Virginia (DEV) substation with two 504 MVA, 500/230 kV transformers. Data collected from these offline simulations were used to train the CNN to provide online GIC monitoring. Testing the CNN performance involved using real GIC measurements from published literature and from a physical GIC monitor in the DEV area. Finally, the results showed that the proposed method was able to provide GIC readings with a root mean squared error of 1.56 A/phase (equivalent to an average accuracy of 94%) for these real GIC waveforms.

42 ENGINEERING↗

Modeling the impact of riparian hollows on river corridor nitrogen exports, Frontiers in Water: Dataset

This package contains all reactive transport simulation outputs used in "Modeling the impact of riparian hollows on river corridor nitrogen exports" (Rogers, D.B., et al., 2021). The simulation outputs reflect hydrological and biogeochemical behavior of a single, generalized riparian hollow located in the intra-meander region of a sub-alpine floodplain under transient hydrological conditions. The context for, purpose of, and conclusions drawn from this model are detailed in Rogers, D.B., et al., 2021.The .zip file contains simulation outputs for each of the five time periods analyzed (see metadata for specific dates), with each folder representing one continuous simulation period. Within each folder are the raw simulation outputs, which are comprised of hundreds of files. Files ending in ".vtk" are best visualized and analyzed using ParaView (https://www.paraview.org). All other files are best viewed or analyzed using R. For a detailed description of the meaning of each file, refer to the ".fls" file located in each folder.

54 ENVIRONMENTAL SCIENCES↗

Separating Propulsive Mass and Energy for Space Applications

Initially and traditionally, space access and in-space propulsion utilized combustion and expulsion of chemical “fuels” carried on board, with performance governed by the rocket equation. These chemicals produced energy, and after combustion, constituted the propulsive mass/momentum. The best such chemicals in terms of Isp that are deemed safe engineering/mission wise are hydrogen and oxygen, producing some 450 seconds of Isp. There are more reactive chemicals, such as fluorine, which have higher Isp, but also have serious safety issues. Propulsion in atmospheres can ingest atmospheric constituents (aka “airbreathing propulsion”) which provide additional propulsive mass when heated and a component of the combustion/energy generation process. This utilization of non-stored/carried propulsive mass provides partial separation of propulsive mass and energy and produces higher Isp. The other approach to separating propulsive mass and energy is to, either on board or added from offboard, supply additional energy, such as from nuclear or solar sources. This is the approach for fission nuclear thermal and electric propulsion, which can provide an Isp in excess of 800 seconds of Isp. Other processes or in addition to thermal expansion, such as electro-magnetics, can be employed to increase exit velocity and Isp. Separating propulsive mass and energy for space faring is commonly referred to as a means to circumvent the rocket equation and is capable of producing major benefits for the development of commercial deep space and space faring in general including affordable fast transits to mitigate the human health impacts of galactic cosmic radiation (GCR) and microG. The key to higher than chemical Isp, beyond H2-O2, and beyond the radiated energy and systems limitations of solar for space faring, is a light weight, high energy density source, either on board or via utilization of energy beaming to the vehicle. Traditional energy sources include chemical, heat, electrics, mechanical, photons, and nuclear. Propulsive mass can be carried on board, sourced beyond the surface of Earth via in-situ resource utilization (ISRU), and includes harvesting from atmospheres. High Isp via electromagnetic related propulsive mass acceleration requires sufficient ionization and conductivity. The major metrics for space propulsion are costs, safety, Isp, weight, and thrust level, the latter dependent upon mission requirements. High thrust for human missions is needed to reduce time exposed to radiation and microG and high thrust is required for space access. Costs of space access are reducing via reusability, printing manufacture, and robotization of manufacturing and operation. Chemical rockets provide high thrust from the expansion of the heated mass constituents at high mass flow. The other high thrust propulsion approach is magnetohydrodynamics (MHD), which, in addition to high thrust, has a high Isp of over 2,000 seconds of Isp. The VASIMIR engine offers some 5,000 seconds of Isp at high thrust [ref. 3]. Electric propulsion cycles are capable of Isp much higher than that but at low thrust levels using available energy sources. With cost as a major metric, reusable rockets and improved manufacturing and operations are rapidly greatly lowering the costs of space access, which could provide/supply in space fuel depots and affordable chemical propulsion for the desired human fast transits (e.g., some 200-day round trips to Mars). The other option for fast transits is VASIMIR, given a nuclear on-board energy source that has the requisite many megawatts of power and an alpha, kgs of weight/kW of energy produced, on the order of one (i.e., a light weight, high energy density energy source). The purpose of this report is to examine the options beyond traditional rocket engines, where propulsive mass and energy are combined, specifically the separation of propulsive mass and energy. This report considers the spectrum of advanced energetics, sources of propulsive mass, conductivity enhancement approaches, energy beaming possibilities, and candidate propulsion cycles. Suggestions are made for various combinatorial, system level, beyond traditional combustion rocket, space propulsion approaches for human deep space missions given the changing conditions of increased knowledge of deep space resources for ISRU, revolutionary energetics, and technology advancements writ large.

Dennis M. Bushnell↗

Spatiotemporal Analyses of Groundwater and Shoreline Cr(VI) Concentrations in the 100 Areas at Hanford

Cleanup efforts have been ongoing since the late 1990s to remediate contaminated waste sites and groundwater in the 100 Areas at the U.S. Department of Energy (DOE) Hanford Site. One of the primary contaminants of concern is hexavalent chromium (Cr(VI)), which was used as a corrosion inhibitor in cooling water for nuclear reactors that formerly operated along the shoreline of the Columbia River. Cleanup efforts have included 1) removal, treatment (as needed), and disposal of contaminated sediments; 2) in situ redox manipulation as a permeable reactive barrier; 3) pump-and-treat; 4) soil flushing; and 5) monitored natural attenuation. DOE’s annual groundwater monitoring reports document the significant reductions in Cr(VI) plume areas that have occurred over the past 10 years or more as a result of these cleanup efforts. The Record of Decision for the 100-HR-3 operable unit specified a cleanup level (CUL) for Cr(VI) in groundwater of 48 µg/L to protect human receptors, and a surface water CUL of 10 µg/L to protect aquatic organisms in the Columbia River. The Record of Decision did not specify point-of-compliance locations for the surface water CUL. Data for 2019 from the six groundwater operable units (OUs) in the 100 Areas indicate that the 48 μg/L groundwater CUL has been achieved in 100% of the wells in the 100-BC and 100-NR OUs, and in 89- 97% of the wells in the other OUs (100-KR, 100-HR-D, 100-HR-H, 100-FR). Data for 2019 indicate that 100% of the aquifer tubes monitored for Cr(VI) in the 100 Areas have concentrations below the 48 μg/L groundwater CUL. However, the 10 μg/L standard has not yet been consistently achieved for both inland groundwater monitoring wells and shoreline aquifer tubes. This report describes a series of data analyses performed to identify consistent relationships, if any, between inland well and shoreline Cr(VI) concentrations within the 100 Areas. To this end, select monitoring data for Cr(VI) measured in groundwater and aquifer tubes at the 100 Areas were analyzed for a 10-year period—2010 to 2019. Relationships between inland groundwater plumes and surface-water points of discharge in and along the Columbia River were examined through several analyses that included inland well and aquifer tube concentrations as a function of distance from the shoreline, evaluation of cumulative probability plots, trend analysis, correlation analysis, cluster analysis, and identification of plume trajectories for each of the 100 Areas. The analyses did not identify consistent relationships between inland groundwater Cr(VI) concentrations and shoreline concentrations within the 100 Areas due to several confounding factors influencing groundwater flow directions and Cr(VI) concentrations. The proximity of groundwater Cr(VI) plumes to the river, and the highly dynamic nature of the river, influence the transport behavior of the plumes and create challenges for quantifying attenuation of Cr(VI) between the inland monitoring wells and shoreline concentrations. Other factors contributing to temporal and spatial Cr(VI) concentrations, as supported by some of the data analyses, include the presence of vadose zone sources, variable sorption behavior, and complexities associated with Cr(VI) mass transfer between the upper and lower aquifers and their interactions with the river. Hence, monitoring to assess compliance with target CULs will need to be determined for each area individually since several factors influence Cr(VI) concentrations in the 100 Areas.

54 ENVIRONMENTAL SCIENCES↗

Interactions between phosphate and arsenic in iron/biochar-treated groundwater: Corrosion control insights from column experiments

An increasing number of studies have reported the coexistence of arsenic (As) and phosphorus at high concentrations in groundwater, which threatens human health and increases the complexity of groundwater remediation. However, limited work has been done regarding As interception in the presence of phosphate in flowing systems. In this study, a series of experiments were conducted to evaluate the interactions between phosphate and As during As removal by iron (Fe)-based biochar (FeBC). The addition of phosphate promoted As removal by FeBC in the batch and column experiments. X-ray absorption near edge structure (XANES) analysis provided evidence of simultaneous oxidation and reduction of trivalent arsenic in the FeBC column experiment, accompanied by corrosive Fe oxidation. However, the addition of phosphate enhanced As stabilization, attributed to the As-incorporated Fe-Ca-phosphates precipitates. The involvement of phosphate decelerated the Fe corrosion and the formation of secondary minerals in the column, mediating the risk of passivation and clogging. The As retained by Fe-Ca-phosphate precipitates was more readily oxidized, resulting in higher proportions of pentavalent arsenic. In conclusion, the results of this work identify the corrosion control and sustained-release roles of phosphate in FeBC application, informing the perspective of FeBC in As-contaminated groundwater remediation and providing new insights into the interactions between phosphate and As.

54 ENVIRONMENTAL SCIENCES↗

Dynamic simulation of CO 2 -injection-induced fault rupture with slip-rate dependent friction coefficient

It is known that poro-elastic stress and effective stress reduction associated with deep underground fluid injection can potentially trigger shear rupture along pre-existing faults. We modeled an idealized CO 2 injection scenario, to assess the effects on faults in the first phase of a generic CO 2 aquifer storage operation. We used coupled multiphase fluid flow and geomechanical numerical modeling to evaluate the stress and pressure perturbations induced by fluid injection and the response of a nearby normal fault. Slip-rate dependent friction and inertial effects have been taken into account during rupture. Contact elements have been used to take into account the frictional behavior of the rupture plane. We investigated different scenarios of injection rate to induce rupture on the fault, employing various fault rheologies. Published laboratory data on CO 2 -saturated intact and crushed rock samples, representative of a potential target aquifer, sealing formation and fault gouge, have been used to define a scenario where different fault rheologies apply at different depths. Nucleation of fault rupture takes place at the bottom of the reservoir, in agreement with analytical poro-elastic stress calculations, depending on injection-induced reservoir inflation and the tectonic stress scenario. For the stress state considered here, the first triggered rupture always produces the largest rupture length and slip magnitude, both of which correlate with the fault rheology. Velocity weakening produces larger ruptures and generates larger magnitude seismic events. Heterogeneous faults have been considered including velocity-weakening or velocity strengthening sections inside and below the aquifer, with the upper sections being velocity-neutral. Nucleation of rupture in a velocity-strengthening section results in a limited rupture extension, both in terms of maximum slip and rupture length. For a heterogeneous fault with nucleation in a velocity-weakening section, the rupture may propagate into the overlying velocity-neutral section, if the extent of velocity-weakening and associated friction drop are large enough.

58 GEOSCIENCES↗

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗

Enhancement of EM Pump Performance through Modeling and Testing

Liquid sodium has seen development as a nuclear reactor coolant in the United States for over 50 years. While sodium has many desirable characteristics as a heat transfer fluid, its reactivity poses numerous engineering challenges for even the most basic plant components like pumps. Despite the challenges, mechanical pumps have seen success in sodium service albeit at the expense of simplicity and flexibility. In fact, significant failures of these pumps have been recorded which can be attributed to their complexity. Notwithstanding the complexity, mechanical pumps are still widely proposed for the main circuit pump of sodium cooled reactors. However, an alternative to mechanical pumps are Electromagnetic Induction Pumps (EMIPs) which use the principle of the Lorentz Force. EMIPs avoid mechanically moving parts in contact with sodium which simplifies the design by avoiding mechanical wear and seals. Therefore, these pumps are generally safer and more reliable than their mechanical counterparts. That being said, EMIPs suffer from relatively low efficiencies that range from 5% to 45%. Some of the inefficiencies are inherent to the EMIPs pumps; namely, inefficiencies arise in the form of resistive heating losses in the fluid itself. However, other inefficiencies arise from Edge Effects which include Finite Length and Finite Width effects. Due to the significant safety advantages EMIPs have compared to mechanical pumps, there is a strong motivation to integrate these technologies into the design of Generation IV liquid-metal-cooled nuclear power plants. However, EMIPs have low efficiencies partly attributed to Edge Effects. Much of the literature focusing on Edge Effects is dated or simply does not exist in the public domain. More recently, numerical models of EMIPs have shown promising results. However, these models are often custom solvers or use complex software suits which limit their application. This report describes the work completed to understand Edge Effects in two EMIPs designs. These designs are the Annular Linear Induction Pump (ALIP) and Permanent Magnet Induction Pump (PMIP) also known as a Moving Magnet Pump (MMP). The methodology developed during this project aim at measuring key performance parameters (pressure-flowrate and efficiency-slip performance) using high temperature liquid-sodium. Baseline experimental measurements were then used to verify the accuracy of numerical models developed in the software suit ANSYS/FLUENT. On the experimental side, instrumentation was developed to characterize the pressure-flowrate and efficiency-slip performance of both the ALIP and PMIP in high-temperature liquid-sodium. Additionally, instrumentation was developed to measure potential output pressure pulsations and the magnetic field distributions. Moreover, a specialized ALIP was constructed to study the impacts magnetic field tapering, also known as coil grading, on the Finite Length Effect. Lastly, experimental performance curves of the PMIP were collected at 300 oC and 600 oC and compared to analytical theory which accounted for Finite Width Effects. These results showed that the Finite Width Effect and non-uniform magnetic field in a PMIP impact the performance. On the numerical side, models were developed using the software suit FLUENT/ANSYS. These methods were initially compared against a simple MHD flow in a Permanent Magnet Electromagnet Flowmeter (PMEMFM). The results were significant in that they showed agreeance between the numerical model and demonstrated that PMEMFM’s could be calibrated numerically. These modeling methods were then applied to a prototypic PMIP. Results showed that this simple model provided reasonably accurate results of the PMIP p-Q performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Fundamental Studies of Bifunctional Catalysts for Tandem Reactions

During the past few years of our funding from DOE Basic Energy Sciences, we have focused on the catalytic conversion of biomass-derived 5-hydroxymethyl furfural (HMF) as a platform molecule for sustainable synthesis of value-added chemicals. The combination of the diene and aldehyde functionalities in HMF enables catalytic production of acetalized HMF derivatives with diol or epoxy reactants to allow reversible synthesis of norcantharimide derivatives upon Diels-Alder reaction with maleimides. Whereas the electron-withdrawing nature of aldehyde group inhibits the reactivity of diene group in HMF to undergo Diels-Alder reactions with dienophiles, by converting the aldehyde group in HMF to an electron-donating hydroxyl group, HMF can append maleimide-based chemicals through Diels-Alder coupling. Acetalization of HMF not only produces an electron-rich diene for Diels-Alder reaction, but it also alters the reversibility of acetal formation and thereby allows for the controlled release of chemicals by hydrolysis. Therefore, we synthesized various norcantharimide derivatives from biomass-derived HMF by acetalization over an acid catalyst, Amberlyst-15, followed by Diels-Alder reaction with maleimide. The norcantharimides release the starting materials by retro Diels-Alder reaction that is triggered by acetal hydrolysis under acidic (≤pH 3) conditions. In subsequent work, we explored the synthesis and properties of functional polyurethanes and polyesters with tunable properties from biomass-derived (HMF)-Acetone-HMF (HAH) monomers. HAH can be selectively hydrogenated over Cu and Ru catalysts to produce partially-hydrogenated HAH (PHAH) and fully-hydrogenated HAH (FHAH). HAH contains functional groups, including hydroxyl, furan, enone, and ketone functionalities, that can be exploited to further tune the polymer properties. The π-electron conjugation between the enone and furan groups in HAH monomer stabilizes the furan group to prevent the uncontrolled degradation of furans during reactions, such as etherification. Therefore, HAH is a renewable monomer that can be used to synthesize polymers, having high molecular weight, symmetric functionalities, an inexpensive production price, and being derived from renewable resources without competition with food resources. We demonstrated the synthesis of functional polyurethanes and a renewable polyester from HAH-derived monomers, and 4,4’-methylenebis(phenyl isocyanate) (MDI). We showed that it is essential to be able to selective hydrogenate C=C double bonds in HAH to be able to achieve effective Diels Alder reactions of the furan ring. Therefore, we studied the hydrogenation at temperatures from 313 – 393 K of HAH over Pd, Ru, and Cu based catalysts. HAH was selectively hydrogenated to produce partially-hydrogenated monomers (PHAH) over Cu and Ru catalysts and to fully-hydrogenated HAH monomers (FHAH) over the Ru catalyst. Pd based catalysts yielded a mixture of partially and fully hydrogenated monomers. Reaction kinetics models were employed to quantify the kinetic behavior for hydrogenation over Ru, Cu, and Pd catalysts. A 5-step pathway exhibited over Pd and Ru catalysts consists of both series and parallel reaction steps, where HAH is both converted to fully hydrogenated products sequentially via series reactions of partially hydrogenated intermediates, as well as converted directly in parallel reactions to form the fully hydrogenated products. In contrast, a 3-step pathway over the Cu catalyst consists only of the consecutive reaction steps, where the final product was formed via series reactions of intermediate products. Additionally, reaction over the Cu catalyst did not hydrogenate the furan rings of the HAH molecule and yielded a different final product than those hydrogenation over Pd and Ru catalysts. Using the results from our detailed reaction kinetics studies, we developed universal conditions for the maximum production of the various hydrogenated products for both batch and plug flow systems. In subsequent work, we employed reaction schemes to optimize the yields of products from selective hydrogenations of HAH in isopropanol solvent and in the presence of liquid water. Reaction schemes consisting of 7, 9, and 11 steps were examined to describe the rates of formation of the observed products and reaction intermediates for hydrogenation of HAH over Ru and Pd catalysts, and a 3-step scheme was studied over Cu catalysts. Rate constants and activation energies were calculated using these reaction schemes, and we then apply these reaction schemes to explore the effects of water addition on the hydrogenation pathways. The effects of water addition to isopropanol (IPA) solvents on the hydrogenation of HAH were markedly different over Pd, Ru, and Cu catalysts. Over the Pd catalyst, the addition of water to IPA increased hydrogenation rates and promoted hydrogenation of furan rings. The addition of water to IPA yielded significant carbon losses over the Ru catalyst, and slowed hydrogenation steps over Cu, while significantly inhibiting hydrogenation of the ketone group. This behavior opened routes toward increased production rates of PHAH=O (a partially hydrogenated form of HAH containing a C=O bond), a product in which the diene groups of the furan rings were not hydrogenated. Addition of water also allowed increased feed concentrations of HAH that were previously not possible in pure IPA solvents.

09 BIOMASS FUELS↗

The impact of depth-dependent water content on steady state weathering and eroding systems

Models for weathering and regolith formation are normally built on the assumption of constant rates of water advection in the zone of water-saturated pores, and constant water content of those pores, but it is common that weathering occurs in the water-unsaturated zone where lateral flow occurs. Hence, water content in pores varies with depth. Here we model mineral weathering profiles while accounting for depth-dependent water content. Like previous models, a mineral equilibrates with water over a length ξ that depends on dissolution and advection rates, but a new lengthscale λ is introduced to describe the decrease of water content with depth. Steady states of the regolith thickness can be attained for any finite λ and non-zero velocity v E of erosion at the land surface. The type of mineral depletion profile developed over geological timescales depends on coupling between weathering and erosion: for slow erosion, a completely developed profile (CDP) is observed, in which the mineral-depleted zone at the top of the regolith has thickness of order λ ; as v E increases, there is a transition to an incompletely developed profile (IDP), in which partial mineral depletion at the land surface is constrained to a narrow range of velocities; when v E exceeds the advance rate of weathering under far from equilibrium conditions, the profile transitions to an unstable regime that exposes bedrock. In general, the reaction front (RF) thickness and the velocity where CPD transitions to IDP depend on the interplay of both water infiltration and chemical equilibration over the timescale of water residence in regolith. The RF thickness roughly equals a correlation length χ defined as half the harmonic average of ξ and λ . In cases of limited water infiltration, water-mineral equilibration is achieved within the length λ , so that the RF thickness is controlled by hydrological properties and is independent of dissolution and advection rates. In the opposite endmember, water infiltrates to large depths and the effects of physical and chemical parameters on RF thickness are the same as in previous geochemical models. The relaxation time for reaching a steady state is shown to be ~ χ / v E . We discuss the effects of physical and chemical parameters in CDPs and IDPs in those endmembers and show an application to a CDP in granitic regolith.

54 ENVIRONMENTAL SCIENCES↗

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↗

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↗