Search NASA⌕ Search

SEARCH · Search NASA

Results for “Step height”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
NASA NTRS1998

Streamer Evaporation

Evaporation is the consequence of slow plasma heating near the tops of streamers where the plasma is only weakly contained by the magnetic field. The form it takes is the slow opening of field lines at the top of the streamer and transient formation of new solar wind. It was discovered in polytropic model calculations, where due to the absence of other energy loss mechanisms in magnetostatic streamers, its ultimate endpoint is the complete evaporation of the streamer. This takes, for plausible heating rates, weeks to months in these models. Of course streamers do not behave this way, for more than one reason. One is that there are losses due to thermal conduction to the base of the streamer and radiation from the transition region. Another is that streamer heating must have a characteristic time constant and depend on the ambient physical conditions. We use our global Magnetohydrodynamics (MHD) model with thermal conduction to examine a few examples of the effect of changing the heating scale height and of making ad hoc choices for how the heating depends on ambient conditions. At the same time, we apply and extend the analytic model of streamers, which showed that streamers will be unable to contain plasma for temperatures near the cusp greater than about 2xl0(exp 6) K. Slow solar wind is observed to come from streamers through transient releases. A scenario for this that is consistent with the above physical process is that heating increases the near-cusp temperature until field lines there are forced open. The subsequent evacuation of the flux tubes by the newly forming slow wind decreases the temperature and heating until the flux tubes are able to reclose. Then, over a longer time scale, heating begins to again refill the flux tubes with plasma and increase the temperature until the cycle repeats itself. The calculations we report here are first steps towards quantitative evaluation of this scenario.

Suess, Steven T.↗
DOE OSTI2026

Quantum utility in simulating the real-time dynamics of the Fermi–Hubbard model using superconducting quantum computers

The Fermi–Hubbard model is a fundamental model in condensed matter physics that describes strongly correlated electrons. On the other hand, quantum computers are emerging as powerful tools for exploring the complex dynamics of these quantum many-body systems. In this work, we demonstrate the quantum simulation of the one-dimensional Fermi–Hubbard model using IBM's superconducting quantum computers, employing over 100 qubits. We introduce a first-order Trotterization scheme and extend it to an optimized second-order Trotterization for the time evolution in the Fermi–Hubbard model, specifically tailored for the limited qubit connectivity of quantum architectures, such as IBM's platforms. Notably, both Trotterization approaches are scalable and maintain a constant circuit depth at each Trotter step, regardless of the qubit count, enabling us to precisely investigate the relaxation dynamics in the Fermi–Hubbard model by measuring the expectation value of the Néel observable (staggered magnetization) for time-evolved quantum states. Lastly, our successful measurement of expectation values in such large-scale quantum many-body systems, especially at longer time scales with larger entanglement, highlights the quantum utility of superconducting quantum platforms over conventional classical approximation methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

I’ve Looked at Clouds from Both Sides, Now: Viewpoints from Surface and Spaceborne Lidar Systems

Clouds play a critical role in the Earth’s climate system because they are inextricably linked to the hydrological cycle and radiation budget. Information about cloud height, thickness, occurrence, and amount are critical inputs for a host of numerical applications involving climate research. Therefore, it is important to have highly accurate and quantitative data records of cloud properties that span several years and geographic regions. Verification of even the most basic modeling processes demands long term and continuous observations of global cloud occurrence, if there is to be any confidence in their fidelity. For a number of reasons, however, it is impossible to gauge the complexities of clouds from a single source. Fundamentally, an array of remote sensing methods is needed in order to provide a complete picture. Within the Micropulse Lidar Network (MPLNET), we have developed a new algorithm to improve the quality of our cloud products. The largest impact of the changes to the cloud detection algorithm is evident with high clouds (those with cloud base > 5 km). Furthermore, polarized measurements enable us to determine the cloud thermodynamic phase. Naturally, the next step is to investigate how these improvements compare with observations from spaceborne lidars (e.g. CALIOP) which have a better (unobstructed) view of high clouds. Recent advances in remote sensing have revealed that cirrus clouds are the most common cloud genus observed in the atmosphere. Furthermore, cirrus skew highly towards relatively low cloud optical depths, as observed from both surface and spaceborne viewpoints. The radiative impacts of these findings are quite significant, considering the cumulative effect cirrus exhibit when compared to low clouds.

Jasper R Lewis↗
NASA NTRS2015

NASA's Meteoroid Environments Office's Response to Three Significant Bolide Events Over North America

Being the only U.S. Government entity charged with monitoring the meteor environment, the Meteoroid Environment Office has deployed a network of all sky and wide field meteor cameras, along with the appropriate software tools to quickly analyze data from these systems. However, the coverage of this network is still quite limited, forcing the incorporation of data from other cameras posted to the internet in analyzing many of the fireballs reported by the public and media. A procedure has been developed that determines the analysis process for a given fireball event based on the types and amount of data available. The differences between these analysis process will be explained and outlined by looking at three bolide events, all of which were large enough to produce meteorites. The first example is an ideal event - a bright meteor that occurred over NASA's All Sky Camera Network on August 2, 2014. With clear video of the event from various angles, a high-accuracy trajectory, beginning and end heights, orbit and approximate brightness/size of the event are able to be found very quickly using custom software. The bolide had the potential to have dropped meteorites, so dark flight analysis and modeling was performed, allowing potential fall locations to be mapped as a function of meteorite mass. The second case study was a bright bolide that occurred November 3, 2014 over West Virginia. This was just north of the NASA southeastern all-sky network, and just south of the Ohio-Pennsylvania network. This case study showcases the MEO's ability to use social media and various internet sources to locate videos of the event from obscure sources (including the Washington Monument) for anything that will permit a determination of a basic trajectory and fireball light curve The third case study will highlight the ability to use doppler weather radar in helping locate meteorites, which enable a definitive classification of the impactor. The input data and analysis steps differ for each case study, but the goals remain the same - a trajectory, orbit, and mass estimate for the bolide within hours of the event, and, for events with a high probability of producing meteorites, a location of the strewn field within a day.

Blaauw, Rhiannon C.↗
NASA NTRS1997

Retrieval of Soil Moisture and Roughness from the Polarimetric Radar Response

The main objective of this investigation was the characterization of soil moisture using imaging radars. In order to accomplish this task, a number of intermediate steps had to be undertaken. In this proposal, the theoretical, numerical, and experimental aspects of electromagnetic scattering from natural surfaces was considered with emphasis on remote sensing of soil moisture. In the general case, the microwave backscatter from natural surfaces is mainly influenced by three major factors: (1) the roughness statistics of the soil surface, (2) soil moisture content, and (3) soil surface cover. First the scattering problem from bare-soil surfaces was considered and a hybrid model that relates the radar backscattering coefficient to soil moisture and surface roughness was developed. This model is based on extensive experimental measurements of the radar polarimetric backscatter response of bare soil surfaces at microwave frequencies over a wide range of moisture conditions and roughness scales in conjunction with existing theoretical surface scattering models in limiting cases (small perturbation, physical optics, and geometrical optics models). Also a simple inversion algorithm capable of providing accurate estimates of soil moisture content and surface rms height from single-frequency multi-polarization radar observations was developed. The accuracy of the model and its inversion algorithm is demonstrated using independent data sets. Next the hybrid model for bare-soil surfaces is made fully polarimetric by incorporating the parameters of the co- and cross-polarized phase difference into the model. Experimental data in conjunction with numerical simulations are used to relate the soil moisture content and surface roughness to the phase difference statistics. For this purpose, a novel numerical scattering simulation for inhomogeneous dielectric random surfaces was developed. Finally the scattering problem of short vegetation cover above a rough soil surface was considered. A general scattering model for grass-blades of arbitrary cross section was developed and incorporated in a first order random media model. The vegetation model and the bare-soil model are combined and the accuracy of the combined model is evaluated against experimental observations from a wheat field over the entire growing season. A complete set of ground-truth data and polarimetric backscatter data were collected. Also an inversion algorithm for estimating soil moisture and surface roughness from multi-polarized multi-frequency observations of vegetation-covered ground is developed.

Sarabandi, Kamal↗
DOE OSTI2021

Modelling a compressible packed bed flow-through washing and deacetylation reactor for corn stover pretreatment

The deacetylation process is a key step in the lignocellulosic biomass pretreatment process that currently is done in batch reactors. A flow-through diffuser type reactor is a promising approach to make the process more efficient and economically favorable. A washing step before deacetylation can further reduce chemical usage and increase the process efficiency. However, the highly compacted biomass particles can cause significant pressure-drop for fluid flow through the bed, which in the worst case can cause self-enhancing pressure-drop leading to reactor blocking. Detailed understanding of the fundamentals of flow through a compressible packed bed is necessary to scale up the washing and deacetylation process. This work closely examines the classical Kozeny Carman equation and proposes modifications to accurately predict the pressure-drop at different compaction states. Here, a compressible packed bed model is further developed based on the interaction between solid mechanics and fluid dynamics. Experimental measurements of water flow through a packed bed of corn stover show good agreement with the developed model. Corn stover washing process scale up simulations show that smaller reactor diameter, higher friction lining material, lower process throughput, and higher process temperature is favorable in stabilizing the bed from extensive compression when the bed height is fixed.

09 BIOMASS FUELS↗
DOE OSTI2020

An Update to the Letter of Intent for MATHUSLA: Search for Long-Lived Particles at the HL-LHC

We report on recent progress in the design of the proposed MATHUSLA Long Lived Particle (LLP) detector for the HL-LHC, updating the information in the original Letter of Intent (LoI), see CDS:LHCC-I-031, arXiv:1811.00927. A suitable site has been identified at LHC Point 5 that is closer to the CMS Interaction Point (IP) than assumed in the LoI. The decay volume has been increased from 20 m to 25 m in height. Engineering studies have been made in order to locate much of the decay volume below ground, bringing the detector even closer to the IP. With these changes, a 100 m x 100 m detector has the same physics reach for large c$\tau$ as the 200 m x 200 m detector described in the LoI and other studies. The performance for small c$\tau$ is improved because of the proximity to the IP. Detector technology has also evolved while retaining the strip-like sensor geometry in Resistive Plate Chambers (RPC) described in the LoI. The present design uses extruded scintillator bars read out using wavelength shifting fibers and silicon photomultipliers (SiPM). Operations will be simpler and more robust with much lower operating voltages and without the use of greenhouse gases. Manufacturing is straightforward and should result in cost savings. Understanding of backgrounds has also significantly advanced, thanks to new simulation studies and measurements taken at the MATHUSLA test stand operating above ATLAS in 2018. We discuss next steps for the MATHUSLA collaboration, and identify areas where new members can make particularly important contributions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1995

The interstellar disk-halo connection in the spiral galaxy NGC 3079

We discuss the morphology and excitation of ionized gas in the nearby Sc galaxy NGC 3079. The almost edge-on orientation is ideal for studying the vertical structure of the gaseous disk, and especially the diffuse ionized medium (DIM) found between the bright H II regions. We used the Hawaii Imaging Fabry-Perot Interferometer (HIFI) to map 150,000 H-alpha + (N II) lambda lambda 6548, 6583 emission-line profiles across the entire disk, with resolution 70 km/s at subarcsecond steps, down to a flux level of approximately 10(exp -17) ergs/s/sq cm (EM approximately equal to 4 cm(exp -6) pc). The DIM contributes approximately 30% of the total disk H-alpha emission within a radius of 10 kpc. The DIM has broader emission lines and larger (N II) H-alpha flux ratios than the adjacent H II regions. Within a radius of 5 kpc, we find that the X-shaped filaments reported in previous studies emerge from the inner (R approximately equal to 1.5 kpc) disk, and rise more than 4 kpc above the disk plane. The morphology, kinematics, and excitation of the filaments suggest that they form a biconic interface between the undisturbed disk gas, and gas entrained in the wide-angle outflow. The DIM beyond 5 kpc radius is more vertically extended than the thick ionized disk detected in our Galaxy and in a few nearby edge-on systems. After correcting for dust, the vertical profile of this DIM has an exponential scale height of about 1.1 kpc, similar to that of the H I disk. The (N II) lambda 6538/H-alpha flux ratio of the DIM increases monotonically with vertical height, reaching unity for absolute value of z greater than or approximately equal to 2.5 kpc. The flux required to keep the DIM ionized at R = 8 kpc is similar to that near the solar circle of our Galaxy. Highly dilute radiation from O stars in the galactic plane probably maintains the DIM. The total mass of the DIM is of order 10(exp 8) - 10(exp 9) solar mass, representing less than 1% of the total dynamical mass of NGC 3079. Mechanical energy from intense star formation in the disk probably lifts the DIM above the disk. The several bubbles and filaments within 1 kpc of the disk plane is direct evidence for gas flow between the disk and halo.

Veilleux, Sylvain↗
NASA NTRS1996

A Sundial-Atlas Precursor to the TIMED Mission: A Quick-Response Global Investigation into Coupled Lower Thermospheric, Ionospheric, and Mesospheric Physics

The SUNDIAL-ATLAS effort was a global-scale investigation which responded to the science priorities of the ITM Panel, the Integrated SPD Strategy Implementation Plan as a whole, and the need for potential cost-saving design criteria for the TIMED mission. The investigation focused on coupling processes in the ionospheric-thermospheric system, taking advantage of the timelines of the ATLAS-1 mission (March 1992), and the global-scale ground-based measurement and modeling activities of the SUNDIAL program. The collaborative SUNDIAL-ATLAS activity was the first opportunity for global measurements of the chemistry, kinetics, and electrodynamics which couple the E-, Fl-, and F2-regions into a single interactive system. As such, the program represented an important first step in studying global issues; and accordingly, was an important proof of concept experiment relevant to the strategic mission plans for the ITM community and the upcoming intermediate class satellite program called TIMED. To meet its projected goals, TIMED must perform a number of critical measurements and execute a number of correlations that were to be tried and tested for the first time in the SUNDIAL-ATLAS investigation. This was designed to include global correlations of thermospheric and ionospheric composition during quiet and disturbed conditions and the co-registration of global-scale ground-based measurements with along-track satellite diagnostics. The SUNDIAL component of the current investigation addressed this need by acquiring, reducing, and analyzing a multi-sensor database that complemented and extended that which was generated in the ATLAS mission (Atmospheric Laboratory for Applications and Science). The SUNDIAL data defined the state and condition of the global-scale ionosphere in the altitude range from 100 km to the F2-peak. These data specified the peak heights and densities of the E-, Fl-, and F2-regions, along with the global distributions of intermediate, descending, and sequential layers which play a critical role in the dynamo region of the lower ionospheric-thermospheric domain. The data were collected by the SUNDIAL network of more than 50 ground-based stations utilizing ionosondes, radars, photometers, Fabry-Perot interferometers, and total electron content measurements. The data were acquired during a three-week period centered on the eight-day ATLAS-1 mission, which provided image and photometric sensing of the altitude distributions of the major and minor ions and neutrals in the ITM system. This report focuses on the scientific contributions of the SUNDIAL component of the overall investigation. Specific findings are described in seven papers (attached) published in the Journal of Geophysical Research.

Szuszczewicz, E. P.↗
NASA NTRS2011

Understanding Geospace on a Grand Scale: The Global Ionoshphere/Thermosphere Constellation

We present the concept of a constellation of polar orbiting satellites equally spaced in longitude (local time) to systematically sample both the neutral and ionzed gas components of the Earth in circular orbits near 350 km, including their density, temperature, and velocities. The instrumentation would include techniques to measure the height of the ionospheric "F-peak" and its variations along the orbit. The number of satellites (l2? 24? 48?) and their configuration would be determined from modeling analysis and expected geophysical phenomena, including their drivers and characteristic time scales. Together with imaging data from separate satellites, the array of satellites with in situ probes would be expected to provide a new picture of (1) high latitude electrodynamics and atmospheric processes and associated coupling with magnetospheric mass and momentum input, (2) the response of the global ionosphere and thermosphere to magnetic storms, and (3) global neutral wind circulation patterns, neutral density structure, tides, planetary waves, and gravity waves. The comprehensive measurements gathered by the IT-Constellation envisioned here would provide a major leap forward in each of these areas, addressing global physical processes and providing fundamental, new knowledge of Geospace. In particular, by its very nature, the constellation addresses "system science", revealing how the ionosphere-thermosphere connects globally to the magnetosphere above and the troposphere below. We present this concept as the next logical step in observing the "whole" space environment using in situ probes in conjunction with imagers. We invite modelers to not only comment on this concept but also to become actively engaged in helping to define it.

Pfaff, Robert↗
NASA NTRS2012

The Surface Water and Ocean Topography Satellite Mission - An Assessment of Swath Altimetry Measurements of River Hydrodynamics

The Surface Water and Ocean Topography (SWOT) satellite mission, scheduled for launch in 2020 with development commencing in 2015, will provide a step-change improvement in the measurement of terrestrial surface water storage and dynamics. In particular, it will provide the first, routine two-dimensional measurements of water surface elevations, which will allow for the estimation of river and floodplain flows via the water surface slope. In this paper, we characterize the measurements which may be obtained from SWOT and illustrate how they may be used to derive estimates of river discharge. In particular, we show (i) the spatia-temporal sampling scheme of SWOT, (ii) the errors which maybe expected in swath altimetry measurements of the terrestrial surface water, and (iii) the impacts such errors may have on estimates of water surface slope and river discharge, We illustrate this through a "virtual mission" study for a approximately 300 km reach of the central Amazon river, using a hydraulic model to provide water surface elevations according to the SWOT spatia-temporal sampling scheme (orbit with 78 degree inclination, 22 day repeat and 140 km swath width) to which errors were added based on a two-dimension height error spectrum derived from the SWOT design requirements. Water surface elevation measurements for the Amazon mainstem as may be observed by SWOT were thereby obtained. Using these measurements, estimates of river slope and discharge were derived and compared to those which may be obtained without error, and those obtained directly from the hydraulic model. It was found that discharge can be reproduced highly accurately from the water height, without knowledge of the detailed channel bathymetry using a modified Manning's equation, if friction, depth, width and slope are known. Increasing reach length was found to be an effective method to reduce systematic height error in SWOT measurements.

Wilson, Matthew D.↗
DOE OSTI2020

GADRAS Batch Inject Tool User Guide

Gamma Detector Response and Analysis Software (GADRAS) is used by the radiation detection and emergency response community to perform modeling and spectral analysis for gamma detector systems. Built into GADRAS is the ability to define a detector, geometry, background characteristics and source composition to generate synthetic spectra for drills and exercises (injects). Consequence Management is currently in development of a sample result data simulator tool in which a deposition model is probed for source conditions at moments in time and locations in space. These values are used to generate realistic sample results for use in drills and exercises. In addition to sample results, there is a need to simulate the actual spectra that would be observed in the field by downlooking HPGe instruments given a deposition activity. This way, the FRMAC Gamma Spectroscopist can practice their process of generating quantified results from spectra on realistic data as well. Recognizing the decades of work done in GADRAS to accurately generate synthetic spectra, this team decided to build a link between the new simulator and GADRAS to generate these spectra quickly and easily. The simulator tool will generate a file that specifies the name of the spectra, its location, date/time of measurement, duration of measurement, height off the ground, and the deposition activity and age for every radionuclide in the simulation. Then, a new tool within the Inject Tab of GADRAS was developed to read in this file given a detector selection and generate In-Situ spectra for each row in the file in any file format the user chooses. This way, simulation cell staff can take these files and then upload them to the appropriate data system (RAMS or RadResponder) for use during drills and exercises. An advanced feature of this tool allows for generating any spectra given an appropriate model and mapping of source to model layer in the batch inject tool. This way, spectra from field sample counts, mobile laboratories, or even fixed laboratories can be generated in bulk given an estimate of the radioactivity concentration or total radioactivity in an import file. This expands the capabilities of this tool a great deal and will make it a more useful tool for CM and others to help estimate detector response for nearly any situation. This user guide will explain the steps needed to perform a batch inject file generation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

High precision magnetic field measurement and mapping of the LEReC 180 degree bending magnet using very low field NMR with hall combined probe (140 - 350 gauss)

The Relativistic Heavy Ion Collider (RHIC) at BNL are using the Low Energy RHIC Electron Cooling (LEReC) to conduct experiments that search for the quantum chromodynamic (QCD) critical point. The first ever electron cooling based on the RF acceleration of electron beams was experimentally demonstrated on April 5, 2019 at LEReC at BNL. The first critical step in obtaining successful 3D non-magnetized cooling of the Au ion bunches in the RHIC cooling section was matching the electron beam energy with a relative error less than 5*10 -4 to the ion beam energy. Part of the LEReC beamline is a dipole magnet that bends the electron beam 180 degree. One of the most outstanding measurement challenges is that the dipole field is so low (≈200 G). Most of the existing NMR probes can only measure fields >400 G. Lower signal-to-noise ratio at low fields is requires the use of larger sample volumes. Working with CAYLAR, a NMR probe has been redesigned and optimized for these low field measurements with high resolution. We report the methods, challenges, and results for extensive magnetic field mappings of the 180 dipole magnet. A combination of NMR and Hall sensors has been successfully implemented to measure uniform field regimes inside the magnet center area and non-uniform field regimes at the magnet ends. Detailed measurement and mapping have been performed at five radii and five heights along the beam trajectory. Meanwhile a finite element magnetic modeling simulation of the magnet using Opera software has been performed. The calculated and measured data are compared, and the calculated data is good reference for measured data over long length mapping from magnet edge to center. The measured magnetic measurement data is directly useful for beam instrumentation, diagnostics and operation.

Magnetic Field Measurement↗
DOE OSTI2020

Long-term effects of thinning and woody control on longleaf pine plantation development, understory abundance, and tree damage from an ice storm

Longleaf pine (Pinus palustris) plantations have been established across the southeastern United States as the first step in restoring the declining savanna ecosystem, yet their long-term stability and development to open canopied woodlands may depend on early density and vegetation management treatments. Commercial thinning (CT) was applied in 2012 (year 18) to a long-term study of longleaf pine community restoration at the Savannah River Site near Aiken, SC to increase abundance of herb species. Although precommercial thinning (PCT) and woody vegetation control increased herb species cover during the first five years of the study, measurements from years 17 to 22 indicated a continuing decline in their cover since year 9 despite application of CT in year 18. Herb species richness, however, increased from 21 to 43 species per 100 m2 area during the last six years of the study, likely a result of disturbances associated with CT and an ice storm in year 20. Application of PCT in 1994, when the plantations were 8–11 years old, reduced their susceptibility to stem bending and breakage from the ice storm, which occurred 20 years later in 2014. Measured 5 months before the storm, pine height: dbh ratio (HD) averaged 78 and 90 for trees growing with and without PCT, respectively, and stem breakage (% of trees) from the ice storm increased linearly (r2 = 0.53) with HD. Pine mortality from the combined effects of glaze ice and wind differed with (28%) versus without (45%) PCT, indicating that long-term stability of longleaf pine plantations depends on early density management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics↗
NASA NTRS2012

Atmospheric Dust Modeling from Meso to Global Scales with the Online NMMB/BSC-Dust Model Part 2: Experimental Campaigns in Northern Africa

The new NMMB/BSC-Dust model is intended to provide short to medium-range weather and dust forecasts from regional to global scales. It is an online model in which the dust aerosol dynamics and physics are solved at each model time step. The companion paper (Perez et al., 2011) develops the dust model parameterizations and provides daily to annual evaluations of the model for its global and regional configurations. Modeled aerosol optical depth (AOD) was evaluated against AERONET Sun photometers over Northern Africa, Middle East and Europe with correlations around 0.6-0.7 on average without dust data assimilation. In this paper we analyze in detail the behavior of the model using data from the Saharan Mineral dUst experiment (SAMUM-1) in 2006 and the Bodele Dust Experiment (BoDEx) in 2005. AOD from satellites and Sun photometers, vertically resolved extinction coefficients from lidars and particle size distributions at the ground and in the troposphere are used, complemented by wind profile data and surface meteorological measurements. All simulations were performed at the regional scale for the Northern African domain at the expected operational horizontal resolution of 25 km. Model results for SAMUM-1 generally show good agreement with satellite data over the most active Saharan dust sources. The model reproduces the AOD from Sun photometers close to sources and after long-range transport, and the dust size spectra at different height levels. At this resolution, the model is not able to reproduce a large haboob that occurred during the campaign. Some deficiencies are found concerning the vertical dust distribution related to the representation of the mixing height in the atmospheric part of the model. For the BoDEx episode, we found the diurnal temperature cycle to be strongly dependant on the soil moisture, which is underestimated in the NCEP analysis used for model initialization. The low level jet (LLJ) and the dust AOD over the Bodélé are well reproduced. The remaining negative AOD bias (due to underestimated surface wind speeds) can be substantially reduced by decreasing the threshold friction velocity in the model.

photometers↗