Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “material flow analysis”

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

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

At least 1,225 records · Page 68

Buoyancy Effects on Concurrent Flame Spread Over Thick PMMA

The flammability of combustible materials in a spacecraft is important for fire safety applications because the conditions in spacecraft environments differ from those on earth. Experimental testing in space is difficult and expensive. However, reducing buoyancy by decreasing ambient pressure is a possible approach to simulate on-earth the burning behavior inside spacecraft environments. The objective of this work is to determine that possibility by studying the effect of pressure on concurrent flame spread, and by comparison with microgravity data, observe up to what point low-pressure can be used to replicate flame spread characteristics observed in microgravity. Specifically, this work studies the effect of pressure and microgravity on upward/concurrent flame spread over 10 mm thick polymethyl methacrylate (PMMA) slabs. Experiments in normal gravity were conducted over pressures ranging between 100 and 40 kPa and a forced flow velocity of 200 mm/s. Microgravity experiments were conducted during NASA’s Spacecraft Fire Experiment (Saffire II), on board the Cygnus spacecraft at 100 kPa with an air flow velocity of 200 mm/s. Results show that reductions of pressure slow down the flame spread over the PMMA surface approaching that in microgravity. The data is correlated in terms of a non-dimensional mixed convection analysis that describes the convective heat transferred from the flame to the solid, and the primary mechanism controlling the spread of the flame. The extrapolation of the correlation to low pressures predicts well the flame spread rate obtained in microgravity in the Saffire II experiments. Similar results were obtained by the authors with similar experiments with a thin composite cotton/fiberglass fabric (published elsewhere). Both results suggest that reduced pressure can be used to approximately replicate flame behavior of untested gravity conditions for the burning of thick and thin solids. This work could provide guidance for potential ground-based testing for fire safety design in spacecraft and space habitats.

Environmental conditions↗

Comparative analysis of characteristics of a supersonic cw chemical HF laser on molecular fluorine and nitrogen trifluoride

We have experimentally studied the characteristics of a supersonic cw chemical HF laser with a flat nozzle array based on a nozzle – nozzle reagent mixing scheme operating with the use of F{sub 2} – D{sub 2} – He and NF{sub 3} – D{sub 2} – He fuel compositions in an atomic fluorine generator. The active medium flow field pattern, its gas-dynamic characteristics, and the laser radiation spectrum are considered. The comparison of these characteristics has made it possible to establish a strong impact on the flow field (lasing zone length) of the chemical composition of the combustion products from the atomic fluorine generator. The gas-dynamic parameters of the active medium formed by using a NF{sub 3} – D{sub 2} – He fuel have turned out to be less acceptable in terms of the gas flow pressure recovery than in the case of F{sub 2} – D{sub 2} – He fuel, while the radiation spectrum was more preferable from the viewpoint of passing through the atmosphere. (paper)

36 MATERIALS SCIENCE↗

Effects of shielding gas and spot reduction on humping suppression in high-speed laser welding

High-speed laser welding is a crucial process to enhance the production rate of fuel cell fabrication. However, exceeding a critical welding speed leads to defect formation such as humps and root concavities, which degrade weld quality and mechanical performance. This study employed shielding gas and spot reduction approaches to suppress these defects. The relationships between individual defect types and their influence on weld strength were systematically examined. Moreover, in situ high-speed synchrotron X-ray imaging was employed to comprehend the mechanisms of humping mitigation by shielding gas through analysis of instantaneous keyhole and melt pool dynamics. The observations reveal that shielding gas mitigates humping by shortening the molten pool through enhanced surface heat convection and by altering surface tension to prevent melt flow collisions. The effects of shielding gas and spot reduction were further incorporated into the dimensionless humping index (π h ) that describes the humping formation tendency. These findings provide mechanistic insights and quantitative guidance for defect control strategies in high-speed laser welding.

08 HYDROGEN↗

Weather Research and Forecasting Model Sensitivity Comparisons for Warm Season Convective Initiation

Mesoscale weather conditions can significantly affect the space launch and landing operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). During the summer months, land-sea interactions that occur across KSC and CCAFS lead to the formation of a sea breeze, which can then spawn deep convection. These convective processes often last 60 minutes or less and pose a significant challenge to the forecasters at the National Weather Service (NWS) Spaceflight Meteorology Group (SMG). The main challenge is that a "GO" forecast for thunderstorms and precipitation at the Shuttle Landing Facility is required at the 90 minute deorbit decision for End Of Mission (EOM) and at the 30 minute Return To Launch Site (RTLS) decision. Convective initiation, timing, and mode also present a forecast challenge for the NWS in Melbourne, FL (MLB). The NWS MLB issues such tactical forecast information as Terminal Aerodrome Forecasts (TAF5), Spot Forecasts for fire weather and hazardous materials incident support, and severe/hazardous weather Watches, Warnings, and Advisories. Lastly, these forecasting challenges can also affect the 45th Weather Squadron (45 WS), which provides comprehensive weather forecasts for shuttle launch, as well as ground operations, at KSC and CCAFS. The need for accurate mesoscale model forecasts to aid in their decision making is crucial. This study specifically addresses the skill of different model configurations in forecasting warm season convective initiation. Numerous factors influence the development of convection over the Florida peninsula. These factors include sea breezes, river and lake breezes, the prevailing low-level flow, and convergent flow due to convex coastlines that enhance the sea breeze. The interaction of these processes produces the warm season convective patterns seen over the Florida peninsula. However, warm season convection remains one of the most poorly forecast meteorological parameters. To determine which configuration options are best to address this specific forecast concern, the Weather Research and Forecasting (WRF) model, which has two dynamical cores - the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM) was employed. In addition to the two dynamical cores, there are also two options for a "hot-start" initialization of the WRF model - the Local Analysis and Prediction System (LAPS; McGinley 1995) and the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS; Brewster 1996). Both LAPS and ADAS are 3- dimensional weather analysis systems that integrate multiple meteorological data sources into one consistent analysis over the user's domain of interest. This allows mesoscale models to benefit from the addition of highresolution data sources. Having a series of initialization options and WRF cores, as well as many options within each core, provides SMG and MLB with considerable flexibility as well as challenges. It is the goal of this study to assess the different configurations available and to determine which configuration will best predict warm season convective initiation.

Watson, Leela R.↗

In-Lab Rapid Analytical Detection of Lunar Volatiles By Universal Gas Analyzer With Comparison to GC-MS System

Introduction: The curation of permanently shadowed regions (PSRs) [1] on the lunar surface is centered around studies based upon the observed volatiles from the LCROSS mission [2]. The rapid detection of important volatile gases and vapors present in planetary bodies and Astromaterials by a standalone analytical device is an area of intense research interest in our group and Planetary Exploration & Astromaterials Research Laboratory (PEARL) facility and this work is relevant to the future preparation of viable lunar simulants for testing curation efforts down the road. The groundbreaking results obtained from the LCROSS Mission [2] open the requirements for the direct detection of volatiles present in regolith materials collected from the lunar surface. The mass spectrometry of volatile chemicals is a general technique that utilizes a set of instruments that creates charged ions from a gaseous chemical species and measures the intensities vs. mass-to-charge ratio (m/z) [3]. In this context, we discuss in-lab experimental results and procedures for rapid qualitative analysis of main LCROSS volatiles (water, H2S, NH3, CO2, and CH3OH) by a Universal Gas Analyzer (UGA) instrument. Additionally, the instrument performance was evaluated by measuring the isotopic abundance ratio of atmospheric Ar-40 to Ar-36 present in room air since, argon is a relevant gas in planetary studies as it can provide an insight and reference point to isotope studies [4]. Additional, cross comparisons were attempted and made between the two instruments to develop a robust analytical technique by comparing mass spectral data for H2S headspace samples with a Trace-1310/ ISQ 7000 (ThermoFisher Scientific.) GC-MS system. Background: The benchtop UGA System is equipped with an SRS UGA 300 quadrupole mass spectrometer designed and built by Stanford Research Systems [5]. This system can be configured for several types of gaseous chemical analysis. The inlet line continuously samples gases at low flow rates (several milliliters per minute) through a capillary limiting the intake pressure making the instrument ideal for online analysis of select gases and/or room atmosphere. Moreover, in our current UGA system, a change in composition at the inlet can be detected in about 200 milliseconds and a complete spectrum is acquired (for a range of 1-100 amu) in under 45 sec with masses measured at rates up to 25 msec per point [5]. This system provides a quick upstream analytical data that we can then compare to results obtained by our GC-MS system. Sample Preparation: Small volume (2-4 mL) of analyte sample was taken in a 10 mL glass vial and sealed with a crimped cap and purged with pure Ar or N2 gas to displace air from the top. The headspace sample was scanned by the UGA instrument at analog, histogram, and pressure vs. time modes. The isotopic abundance ratio for 40Ar-to-36Ar was estimated by measuring partial pressure vs. time scans and setting the mass at 40 and 36 respectively. Results and Discussions: In this work, we have investigated the applicability of the UGA system by qualitative analysis of a series of LCROSS volatiles measured individually. Fig. 1 demonstrates a set of vertically offset spectra for the partial pressures measured as a function of mass-to-charge (m/z) ratios. The average acquisition time for each spectrum was less than a minute suggesting that the UGA system is ideal for quick analysis of geochemical volatiles. For the cross-comparison, we analyzed an H2S headspace sample by a Trace-1310/ISQ-7000 system and compared mass spectral data with previously measured UGA histogram scan data (Fig. 2). In both cases, major peak positions are the same, however, the intensities of fragment ions ([1H132S]+ and [32S]+) are higher for UGA suggesting that the fragment ionization process is stronger in UGA compared to that of GC-MS. To investigate how the integrated area under each chromatogram varies with the headspace sample volume, a set of five H2S headspace samples with increasing volumes was analyzed by the GC-MS system (Fig. 3, inset). A small volume (e.g., 200 to 1000 µL) of H2S/H2O vapor was withdrawn from a 20 mL stock sample vial containing ~5 mL of 0.4% H2S in water by a gas-tight syringe and added to another 20 mL vial filled with argon and analyzed by the GC-MS system. Finally, the UGA detector sensitivity was evaluated by calculating the atmospheric 40Ar-to-36Ar isotopic abundance ratio in room air by running a partial pressure vs. time scan with setting the atomic mass at 40, and 36. Fig. 4(a) shows a ~25 min duration “P vs. time” scan for 40Ar (plot for 36Ar is not shown). The partial pressure values (~100 points) were corrected by subtracting the corresponding background pressure value for 37Ar and utilized to calculate 40Ar-to-36Ar isotopic abundance ratios as shown by Fig. 4b. The average isotopic abundance ratio is ~306 with a 2*STDEV ~13. This abundance ratio is significantly close to the previously reported value of 298.56 [6] and the ratio obtained by our GC-MS system (303 for a UHP grade Ar sample). Conclusions: Our study strongly evidenced that the benchtop UGA system is a valuable analytical tool for the detection of major LCROSS volatiles. The rapid scanning capability, the inexpensiveness of the whole system, and impressive detection sensitivity prove its worthiness as an essential device for advanced geochemical applications. Moreover, cross comparisons with the GC-MS provide important bridges into advanced curatorial efforts into the future. References: [1] Bickel, V.T., et al. (2021) Nat Commun 12, 5607. [2] Colaprete, A., et al. (2010) Science, 330, 463-468. [3] Glavin, D. P. et al. (2012) 2012 IEEE Aerospace Conference, 1-11. [4] Willett, C. D., et al. (2022) Geochimica et Cosmochimica Acta 329, 119-134. [5] Operation Manual and Programming Reference. (2018) Universal gas Analyzers, Stanford Research Systems. [6] Lee, J. Y., et al. (2006) Geochimica et Cosmochimica Acta 70, 4507–4512. Notes: (4 figures are attached with text as shown by the attached file)

Curation↗

A Preliminary Assessment of Step Effects on the BOLT Geometry in Mach 6 Flow

An infrared thermography measurement technique has been installed in the Langley Aerothermodynamic Laboratory. The system is calibrated for use with commercially available thermoplastics and ceramic materials, and to the surface temperature ranges when exposed to the facility hypersonic flow. Herein an uncertainty quantification is performed, revealing 5% surface temperature and 20% heat transfer uncertainty bands to two-sigma confidence. As a first application of the test measurement technique, an investigation was conducted to explore the correlation between the height of forward facing (FF) and rearward facing (RF) two-dimensional steps, free stream Reynolds number, and the onset of boundary layer transition from laminar to turbulent in Mach 6 flow on a 50% scale Boundary Layer Transition (BOLT) flight test vehicle geometry. Computational predictions were generated to identify potential x/L locations of transition onset for various flow conditions with a smooth outer mold line geometry and are compared to experimental results. The BOLT Step Test was conducted in the NASA Langley Aerothermodynamic Laboratory 20-inch Mach 6 wind tunnel at Reynolds numbers between 1.4 and 11.3 million conditions for various step heights. Data were collected using a developmental infrared thermography test technique and reduced to heat transfer for comparison with CFD solutions. This is a preliminary assessment of the developmental infrared thermography test technique used to acquire and reduce the experimental data. insight into the ongoing comparative analysis of the effects of FF and RF steps on boundary layer transition for the BOLT geometry in Mach 6 flow.

BOLT↗

A Preliminary Assessment of Step Effects on the BOLT Geometry in Mach 6 Flow

An infrared thermography measurement technique has been installed in the Langley Aerothermodynamic Laboratory. The system is calibrated for use with commercially available thermoplastics and ceramic materials, and to the surface temperature ranges when exposed to the facility hypersonic flow. Herein an uncertainty quantification is performed, revealing 5% surface temperature and 20% heat transfer uncertainty bands to two-sigma confidence. As a first application of the test measurement technique, an investigation was conducted to explore the correlation between the height of forward facing (FF) and rearward facing (RF) two-dimensional steps, free stream Reynolds number, and the onset of boundary layer transition from laminar to turbulent in Mach 6 flow on a 50% scale Boundary Layer Transition (BOLT) flight test vehicle geometry. Computational predictions were generated to identify potential x/L locations of transition onset for various flow conditions with a smooth outer mold line geometry and are compared to experimental results. The BOLT Step Test was conducted in the NASA Langley Aerothermodynamic Laboratory 20-inch Mach 6 wind tunnel at Reynolds numbers between 1.4 and 11.3 million conditions for various step heights. Data were collected using a developmental infrared thermography test technique and reduced to heat transfer for comparison with CFD solutions. This is a preliminary assessment of the developmental infrared thermography test technique used to acquire and reduce the experimental data. insight into the ongoing comparative analysis of the effects of FF and RF steps on boundary layer transition for the BOLT geometry in Mach 6 flow.

BOLT↗

Computational Analysis of Static and Dynamic Behaviour of Magnetic Suspensions and Magnetic Bearings

Static modelling of magnetic bearings is often carried out using magnetic circuit theory. This theory cannot easily include nonlinear effects such as magnetic saturation or the fringing of flux in air-gaps. Modern computational tools are able to accurately model complex magnetic bearing geometries, provided some care is exercised. In magnetic suspension applications, the magnetic fields are highly three-dimensional and require computational tools for the solution of most problems of interest. The dynamics of a magnetic bearing or magnetic suspension system can be strongly affected by eddy currents. Eddy currents are present whenever a time-varying magnetic flux penetrates a conducting medium. The direction of flow of the eddy current is such as to reduce the rate-of-change of flux. Analytic solutions for eddy currents are available for some simplified geometries, but complex geometries must be solved by computation. It is only in recent years that such computations have been considered truly practical. At NASA Langley Research Center, state-of-the-art finite-element computer codes, 'OPERA', 'TOSCA' and 'ELEKTRA' have recently been installed and applied to the magnetostatic and eddy current problems. This paper reviews results of theoretical analyses which suggest general forms of mathematical models for eddy currents, together with computational results. A simplified circuit-based eddy current model proposed appears to predict the observed trends in the case of large eddy current circuits in conducting non-magnetic material. A much more difficult case is seen to be that of eddy currents in magnetic material, or in non-magnetic material at higher frequencies, due to the lower skin depths. Even here, the dissipative behavior has been shown to yield at least somewhat to linear modelling. Magnetostatic and eddy current computations have been carried out relating to the Annular Suspension and Pointing System, a prototype for a space payload pointing and vibration isolation system, where the magnetic actuator geometry resembles a conventional magnetic bearing. Magnetostatic computations provide estimates of flux density within airgaps and the iron core material, fringing at the pole faces and the net force generated. Eddy current computations provide coil inductance, power dissipation and the phase lag in the magnetic field, all as functions of excitation frequency. Here, the dynamics of the magnetic bearings, notably the rise time of forces with changing currents, are found to be very strongly affected by eddy currents, even at quite low frequencies. Results are also compared to experimental measurements of the performance of a large-gap magnetic suspension system, the Large Angle Magnetic Suspension Test Fixture (LAMSTF). Eddy current effects are again shown to significantly affect the dynamics of the system. Some consideration is given to the ease and accuracy of computation, specifically relating to OPERA/TOSCA/ELEKTRA.

Britcher, Colin P.↗

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance, Facility Pressure Effects, and Stability Characterization Tests of NASA's Hall Effect Rocket with Magnetic Shielding Thruster

NASAs Hall Effect Rocket with Magnetic Shielding (HERMeS) 12.5 kW Technology Demonstration Unit-1 (TDU-1) Hall thruster has been the subject of extensive technology maturation in preparation for flight system development. Part of the technology maturation effort included experimental evaluation of the TDU-1 thruster with conducting and dielectric front pole cover materials in two different electrical configurations. A graphite front pole cover thruster configuration with the thruster body electrically tied to cathode and an alumina front pole cover thruster configuration with the thruster body floating were evaluated. Both configurations were also evaluated at different facility background pressure conditions to evaluate background pressure effects on thruster operation. Performance characterization tests found that higher thruster performance was attained with the graphite front pole cover configuration with the thruster electrically tied to cathode. A total thrust efficiency of 68 and a total specific impulse of 2,820 s was demonstrated at a discharge voltage of 600 V and a discharge power of 12.5 kW. Thruster stability regimes were characterized with respect to the thruster discharge current oscillations and with maps of the current-voltage-magnetic field (IVB). Analysis of TDU-1 discharge current waveforms found that lower normalized discharge current peak-to-peak and root mean square magnitudes were attained when the thruster was electrically floated with alumina front pole covers. Background pressure effects characterization tests indicated that the thruster performance and stability was mostly invariant to changes in the facility background pressure for vacuum chamber pressure below 110-5 Torr-Xe (for thruster flow rate above 8 mgs). Power spectral density analysis of the discharge current waveform showed that increasing the vacuum chamber background pressure resulted in a higher discharge current dominant frequency. Finally the IVB maps of the TDU-1 thruster taken at elevated magnetic fields indicated that the discharge current became more oscillatory with increased facility background pressure at lower thruster mass flow rates, where thruster operation at higher flow rates resulted in less change to the thrusters IVB characteristics.

Hall Thrusters↗

Advanced Multi-phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

A Navier-Stokes code, finite difference Navier-Stokes (FDNS), is used to analyze the complicated internal flowfield of the SRM (solid rocket motor) to explore the impacts due to the effects of chemical reaction, particle dynamics, and slag accumulation on the solid rocket motor (SRM). The particulate multi-phase flowfield with chemical reaction, particle evaporation, combustion, breakup, and agglomeration models are included in present study to obtain a better understanding of the SRM design. Finite rate chemistry model is applied to simulate the chemical reaction effects. Hermsen correlation model is used for the combustion simulation. The evaporation model introduced by Spalding is utilized to include the heat transfer from the particulate phase to the gase phase due to the evaporation of the particles. A correlation of the minimum particle size for breakup expressed in terms of the Al/Al2O3 surface tension and shear force was employed to simulate the breakup of particles. It is assumed that the breakup occurs when the Weber number exceeds 6. A simple L agglomeration model is used to investigate the particle agglomeration. However, due to the large computer memory requirements for the agglomeration model, only 2D cases are tested with the agglomeration model. The VOF (Volume of Fluid) method is employed to simulate the slag buildup in the aft-end cavity of the redesigned solid rocket motor (RSRM). Monte Carlo method is employed to calculate the turbulent dispersion effect of the particles. The flowfield analysis obtained using the FDNS code in the present research with finite rate chemical reaction, particle evaporation, combustion, breakup, agglomeration, and VOG models will provide a design guide for the potential improvement of the SRM including the use of materials and the shape of nozzle geometry such that a better performance of the SRM can be achieved. The simulation of the slag buildup in the aft-end cavity can assist the designer to improve the design of the RSRM geometry.

Liaw, Paul↗

Towards the Ultimate Strength of Medium‐Entropy Alloys Through Pulsed Lasers

The tensile strength of metals at extreme strain rates is a key predictor of their performance in ballistic and structural impact applications. An important experimental method to reach these extreme strain rates is the use of high-amplitude, short-duration pulsed lasers. The Jupiter Laser Facility at the Lawrence Livermore National Laboratory enabled probing for the first time the mechanical response of several promising High Entropy Alloys at times on the order of nanoseconds (strain rates of ∼10 7 and ∼10 9 s −1 ). The measured strength is in the range of 6 to 10 GPa, ten times the quasistatic value. The mechanisms of plastic deformation and failure were identified and quantified through analysis and molecular dynamics simulation. The reflected wave amplitudes, obtained by VISAR, were used to determine the tensile (spall) stress. The high tensile strength obtained is due to two factors: the strain-rate dependence of plastic flow and the kinetics of void nucleation, growth, and coalescence. The experimental results are compared with an analytical prediction considering both grain-interior and grain-boundary void initiation. Molecular dynamics simulations, conducted at strain rates of 10 8 and 10 9 s −1 , rationalize the experimental results. In conclusion, they provide valuable information about the process of failure evolution, and reveal that grain boundary separation plays a pivotal role in spalling.

Materials science↗

THERMAL/STRUCTURAL ANALYSIS OF THE AXION QUANTUM METROLOGY CAVITY AND ITS COMPONENTS

This research was centered around maximizing the capability to cool dielectric material within a containment unit, or Photonic Band Gap (PBG) cavity, designed for detecting axion dark matter and identifying the unit’s thermal properties. There are multiple types of PBG cavities, but the latest version that axion researchers wish to use has been theorized to contain possible issues related to its thermal properties. Thermal conductivity is an issue with the dielectric material because it is made from alumina which is highly insulative. This is important since the research is being done in a cryogenic environment and the thermal noise affects the quantum bit used for detecting the axion to photon conversion process. Therefore, any improvements to this unit should be justified and implemented but are not entirely limited to thermal contact related aspects of PBG cavities. A prospect of using sapphire in place of alumina also exists, but this is a more expensive, less tested, and more elusive material to justify the creation of a full dielectric structure out of sapphire. Thermal aspects of the cavity were analyzed using finite element method (FEM) and an experiment designed to test different thermal joint materials. FEM was used to check the contraction of the cavity during cooling, the contact quality between the dielectric material and cavity wall, theheat flow rate through the assembly, and helps visualize the cavity’s reaction to different design changes. A simple comparison between thermal conductivity curves justified the usage of sapphire over alumina from a thermal property point of view. Thus, the method of experimentation had an additional dielectric material, sapphire, to test with and compare to alumina. The thermal test identified which material is best to use for a thermal joint but also simultaneously found the conductance of the joint and dielectric material as well as an estimate for what the temperature is inside the larger scale cavity since a thermometer cannot be placed inside the PBG when it cools. Initially, the cavity simulation was tested for structural deformation properties and stress distribution since it is made of copper which contracts heavily in comparison to most other materials. The reduction in volume of the cavity was less than anticipated which gave more room for a possible usage of a modified thermal joint for additional contact area and higher conductivity.

43 PARTICLE ACCELERATORS↗

Adsorptive denitrogenation of model aviation fuel using mesoporous silica in a packed bed adsorption system

This study aims to understand the effects of system process parameters such as flow rate, adsorbent particle size, and use of recycled adsorbent on denitrogenation performance of a model fuel using mesoporous silica gel. The goal is to reduce the nitrogen content of the model fuel from 1500 parts per million to single-digit ppm to meet ASTM specifications for drop-in fuels. This work was done with the intent of applying adsorptive denitrogenation to sustainable aviation fuel (SAF) product fractions produced via hydrothermal liquefaction (HTL). The adsorption performance of the silica is evaluated via packed column breakthrough data, with data generated from collecting from the column outlet and quantifying nitrogen content via gas chromatography. Select experiments use a significantly larger (2.5x column diameter and length) column to demonstrate linear scalability of the process. Thermogravimetric analysis data is collected to evaluate the effects of thermal calcination as a sorbent regeneration method. Effects of a more complex feed are also investigated using a known reference fuel with additional added nitrogen containing compounds. The results presented in this work successfully demonstrate up to 99.8 % removal of NCCs from a model fuel fraction at an original NCC concentration of approximately 1500 ppm to single-digit parts per million after treatment. We also examine calcination of sorbent materials to remove the adsorbed species to enable sorbent reuse and minimize waste generation and show that the calcined material can be reused up to 5 cycles with reduced adsorption capacity. Overall, this work indicates that adsorptive denitrogenation using silica gel is a viable solution to enable the integration of HTL-derived aviation fuels into existing fuel infrastructure.

Adsorption techniques↗

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗

Active Volcanoes of Kamchatka as Suitable Terrestrial Analogs for Future Landing Missions and for Change Detection Analysis on Venus

The next decade will see the return to Venus due to a number of missions which have been recently selected and proposed for launch. The selected missions NASA DAVINCI and VERITAS, Roscosmos Venera-D, and ESA EnVision, along with the proposed missions ISRO Shukrayaan- 1, and CNSA VOICE, will open a new era for the exploration of the Earth’s hellish twin planet. These missions will shed new light on the science questions such as: (a) whether recent volcanic activity is locally constrained, or active volcanism occurs on a global scale; (b) the rate of the present-day volcanic activity; (c) the styles of volcanism on Venus, including possible occurrences of pyroclastic/explosive volcanism; and (d) how does the extreme environment on Venus effect volcanic activity. In preparation for these future missions to Venus, the Analogs for VENus’ GEologically Recent Surfaces (AVENGERS) initiative has the scope to select and analyze a number of active terrestrial volcanoes as suitable analogs for the identification and analysis of active volcanism on Venus. The Roscosmos Venera-D is the only mission, among the currently selected or planned mission, equipped with a lander to analyze the elemental and mineralogical composition of the surface. Therefore, it is crucial to find applicable volcanic areas on Earth to perform operational tests such as surface drilling and in-situ elemental composition analysis. To this regard, we propose the analysis of active volcanoes of the Kamchatka Peninsula, which is in a geodynamic setting (subduction) contrasting to that of the volcanic rises on Venus (hot spot). Thus, Kamchatka may represent a helpful geodynamic endmember to put the expected landing site of the Venera-D mission within the context of the entire variety of volcanic landforms and surface chemistry. The Kamchatka Peninsula is located on the eastern margin of the Eurasian plate, in proximity to the Kuril-Kamchatka subduction zone. The frequent eruptions characterizing the volcanoes of the Kamchatka Peninsula make them a suitable terrestrial analog for the possible detection of ongoing eruptions on Venus by future missions, which can be achieved by comparing a sequence of radar images of the same volcano and its surroundings taken at different times. Moreover, it is possible to collect samples of both unweathered and weathered volcanic deposits, to analyze their near infrared spectra in the laboratory for direct comparison with the spectra to be provided by the future missions to Venus. The chemical analyses of Venusian surface materials from the previous Soviet Venera and Vega missions are consistent with that of tholeiitic basalts. Since the volcanism of the Kamchatka Peninsula is characterized by both basaltic lava flows eruptions and pyroclastic activity, performing operational tests in the Kamchatka Peninsula also offers the unique opportunity to provide us the tools to potentially interpret the diverse output given by both effusive and explosive volcanic products.

Piero D'Incecco↗

Conformal maps of viscous electron flow in the Gurzhi crossover

Herein we investigate the impact of geometric constriction on the viscous flow of electron liquid through quantum point contacts. We provide analysis on the electric potential distribution given the setup of a slit configuration and use the method of conformal mapping to obtain analytical results. The potential profile can be tested and contrasted experimentally with the scanning tunneling potentiometry technique. We discuss intricate physics that underlies the Gurzhi effect, i.e., the enhancement of conductivity in the viscous flow, and compare results for different boundary conditions. In addition, we calculate the temperature dependence of the momentum relaxation time as a result of impurity assisted quasiballistic interference effects and discuss various correlational corrections that lead to the violation of Matthiessen's rule in the hydrodynamic regime. We caution that spatially inhomogeneous profiles of current in the Gurzhi crossover between Ohmic and Stokes flows might also appear in the nonhydrodynamic limit where nonlocality plays an important role. This conclusion is corroborated by calculation of dispersive conductivity in the weakly impure limit.

2-dimensional systems↗

Burst Pressure Solutions of Thin and Thick-Walled Cylindrical Vessels

Pressure vessels (PVs) are widely used in the energy industry. Accurate burst pressure is critical to structural design and safe operation for both thin and thick-walled PVs. The traditional strength theories utilized a single-parameter material property, such as the yield stress or the ultimate tensile stress (UTS) to develop failure models for determining the yield or ultimate pressure carrying capacity in the PV design. The UTS-based Barlow formula is a typical burst pressure model developed from the Tresca strength theory that provides the basis for developing regulation rules and failure models for different industry design codes, such as ASME BPVC, ASME B31.3, and ASME B31G, among others. In order to reduce the conservatism of the Tresca strength model, ASME BPVC recently adapted failure models developed from the von Mises strength theory for the PV design and analysis. It has been commonly accepted that the burst pressure of pipelines depends on the UTS and strain hardening exponent, n, of the pipeline steel. An average shear stress yield theory was thus developed, and the Zhu-Leis solution of burst pressure was obtained as a function of UTS and n for thin-walled line pipes. Experiments showed that the Zhu-Leis solution provides an accurate, reliable prediction of burst pressure for defect-free thin-walled pipes. In order to extend the Zhu–Leis solution to thick-walled cylindrical PVs, this paper defined three new flow stresses, modified the traditional strength theories, and obtained three new burst pressure solutions that are valid for both thin and thick-walled cylindrical vessels. The proposed flow stresses are able to describe the tensile strength and the plastic flow response of PVs for a strain hardening steel. The associated strength theories were then developed in terms of the Tresca, von Mises, and Zhu-Leis yield criteria. From these new strength theories, three burst pressure solutions were obtained for thick-walled cylinders, where the von Mises solution is an upper bound prediction, the Tresca solution is a lower bound prediction, and the Zhu-Leis solution is an intermediate prediction of burst pressure for thick-walled cylinders. Lastly, the proposed burst pressure solutions were evaluated and validated by two large datasets of full-scale burst tests for thick-walled tubes and for thin-walled pipes.

42 ENGINEERING↗