Search NASA⌕ Search

SEARCH · Search NASA

Results for “Domain Reduction”

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.

396 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Institutional Memory Preservation at NASA Glenn Research Center

In this era of downsizing and deficit reduction, the preservation of institutional memory is a widespread concern for U.S. companies and governmental agencies. The National Aeronautical and Space Administration faces the pending retirement of many of the agency's long-term, senior engineers. NASA has a marvelous long-term history of success, but the agency faces a recurring problem caused by the loss of these engineers' unique knowledge and perspectives on NASA's role in aeronautics and space exploration. The current work describes a knowledge elicitation effort aimed at demonstrating the feasibility of preserving the more personal, heuristic knowledge accumulated over the years by NASA engineers, as contrasted with the "textbook" knowledge of launch vehicles. Work on this project was performed at NASA Glenn Research Center and elsewhere, and focused on launch vehicle systems integration. The initial effort was directed toward an historic view of the Centaur upper stage which is powered by two RL-10 engines. Various experts were consulted, employing a variety of knowledge elicitation techniques, regarding the Centaur and RL-10. Their knowledge is represented in searchable Web-based multimedia presentations. This paper discusses the various approaches to knowledge elicitation and knowledge representation employed, and assesses successes and challenges in trying to perform large-scale knowledge preservation of institutional memory. It is anticipated that strategies for knowledge elicitation and representation that have been developed in this grant will be utilized to elicit knowledge in a variety of domains including the complex heuristics that underly use of simulation software packages such as that being explored in the Expert System Architecture for Rocket Engine Numerical Simulators.

Coffey, J.↗
NASA NTRS1994

Dual-Use Aspects of System Health Management

System Health Management functionality is an essential part of any space launch system. Health management functionality is an integral part of mission reliability, since it is needed to verify the reliability before the mission starts. Health Management is also a key factor in life cycle cost reduction and in increasing system availability. The degree of coverage needed by the system and the degree of coverage made available at a reasonable cost are critical parameters of a successful design. These problems are not unique to the launch vehicle world. In particular, the Intelligent Vehicle Highway System, commercial aircraft systems, train systems, and many types of industrial production facilities require various degrees of system health management. In all of these applications, too, the designers must balance the benefits and costs of health management in order to optimize costs. The importance of an integrated system is emphasized. That is, we present the case for considering health management as an integral part of system design, rather than functionality to be added on at the end of the design process. The importance of maintaining the system viewpoint is discussed in making hardware and software tradeoffs and in arriving at design decisions. We describe an approach to determine the parameters to be monitored in any system health management application. This approach is based on Design of Experiments (DOE), prototyping, failure modes and effects analyses, cost modeling and discrete event simulation. The various computer-based tools that facilitate the approach are discussed. The approach described originally was used to develop a fault tolerant avionics architecture for launch vehicles that incorporated health management as an integral part of the system. Finally, we discuss generalizing the technique to apply it to other domains. Several illustrations are presented.

Owens, P. R.↗
NASA NTRS2002

Aeroelastic Calculations of Quiet High- Speed Fan Performed

An advanced high-speed fan was recently designed under a cooperative effort between the NASA Glenn Research Center and Honeywell Engines & Systems. The principal design goals were to improve performance and to reduce fan noise at takeoff. Scale models of the Quiet High-Speed Fan were tested for operability, performance, and acoustics. During testing, the fan showed significantly improved noise characteristics, but a self-excited aeroelastic vibration known as flutter was encountered in the operating range. Flutter calculations were carried out for the Quiet High-Speed Fan using a three-dimensional, unsteady aerodynamic, Reynolds-averaged Navier-Stokes turbomachinery code named "TURBO." The TURBO code can accurately model the viscous flow effects that can play an important role in various aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), and flutter in the presence of shock and boundary-layer interaction. Initially, calculations were performed with no blade vibrations. These calculations were at a constant rotational speed and a varying mass flow rate. The mass flow rate was varied by changing the backpressure at the exit boundary of the computational domain. These initial steady calculations were followed by aeroelastic calculations in which the blades were prescribed to vibrate harmonically in a natural mode, at a natural frequency, and with a fixed interblade phase angle between adjacent blades. The AE-prep preprocessor was used to interpolate the in-vacuum mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh and to smoothly propagate the grid deformations from the blade surface to the interior points of the grid. The aeroelastic calculations provided the unsteady aerodynamic forces on the blade surface due to blade vibrations. These forces were vector multiplied with the structural dynamic mode shape to calculate the work done on the blade during one vibration period, then this result was converted to an aerodynamic damping. Flutter occurs when the aerodynamic damping becomes negative, if structural damping is ignored. The results of these aeroelastic calculations are summarized in this plot of aerodynamic damping versus mass flow rate at a constant rotational speed. As the backpressure is increased, the mass flow rate through the fan decreases and the fan operating point moves towards the stall line. The aeroelastic calculations showed that the aerodynamic damping decreases as the stall line is approached, as observed during testing. In addition, the aeroelastic calculations with the TURBO code correctly predicted the aeroelastic parameters: the most unstable vibration mode and interblade phase angle, as observed during testing. The Quiet High-Speed Fan demonstrated significant noise reductions during testing, but flutter imposed limits on its operating range. The accurate calculation of the aeroelastic characteristics using the TURBO code is a significant step toward eliminating flutter from the operating range and toward realizing the benefits of reduced fan noise. The aeroelastic calculations described here were performed under a grant by University of Toledo researchers in collaboration with Glenn's researchers.

Bakhle, Milind A.↗
NASA NTRS2022

Cloud Adjustments From Large-Scale Smoke–Circulation Interactions Strongly Modulate the Southeastern Atlantic Stratocumulus-to-Cumulus Transition

Smoke from southern Africa blankets the southeastern Atlantic Ocean from June to October, producing strong and competing aerosol radiative effects. Smoke effects on the transition between overcast stratocumulus and scattered cumulus clouds are investigated along a Lagrangian (air-mass-following) trajectory in regional climate and large eddy simulation models. Results are compared with observations from three recent field campaigns that took place in August 2017: ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES), CLouds and Aerosol Radiative Impacts and Forcing: Year 2017 (CLARIFY), and Layered Atlantic Smoke Interactions with Clouds (LASIC). The case study is set up around the joint ORACLES–CLARIFY flight that took place near Ascension Island on 18 August 2017. Smoke sampled upstream on an ORACLES flight on 15 August 2017 likely entrained into the marine boundary layer later sampled during the joint flight. The case is first simulated with the WRF-CAM5 regional climate model in three distinct setups: (1) FireOn, in which smoke emissions and any resulting smoke–cloud–radiation interactions are included; (2) FireOff, in which no smoke emissions are included; (3) RadOff, in which smoke emissions and their microphysical effects are included but aerosol does not interact directly with radiation. Over the course of the Lagrangian trajectory, differences in free tropospheric thermodynamic properties between FireOn and FireOff are nearly identical to those between FireOn and RadOff, showing that aerosol–radiation interactions are primarily responsible for the free tropospheric effects. These effects are non-intuitive: in addition to the expected heating within the core of the smoke plume, there is also a “banding” effect of cooler temperature (∼1–2 K) and greatly enhanced moisture (>2 g kg−1) at the plume top. This banding effect is caused by a vertical displacement of the former continental boundary layer in the free troposphere in the FireOn simulation resulting from anomalous diabatic heating due to smoke absorption of sunlight that manifests primarily as a few hundred meters per day reduction in large-scale subsidence over the ocean. A large eddy simulation (LES) is then forced with free tropospheric fields taken from the outputs for the WRF-CAM5 FireOn and FireOff runs. Cases are run by selectively perturbing one variable (e.g., aerosol number concentration, temperature, moisture, vertical velocity) at a time to better understand the contributions from different indirect (microphysical), “large-scale” semi-direct (above-cloud thermodynamic and subsidence changes), and “local” semi-direct (below-cloud smoke absorption) effects. Despite a more than 5-fold increase in cloud droplet number concentration when including smoke aerosol concentrations, minimal differences in cloud fraction evolution are simulated by the LES when comparing the base case with a perturbed aerosol case with identical thermodynamic and dynamic forcings. A factor of 2 decrease in background free tropospheric aerosol concentrations from the FireOff simulation shifts the cloud evolution from a classical entrainment-driven “deepening–warming” transition to trade cumulus to a precipitation-driven “drizzle-depletion” transition to open cells, however. The thermodynamic and dynamic changes caused by the WRF-simulated large-scale adjustments to smoke diabatic heating strongly influence cloud evolution in terms of both the rate of deepening (especially for changes in the inversion temperature jump and in subsidence) and in cloud fraction on the final day of the simulation (especially for the moisture “banding” effect). Such large-scale semi-direct effects would not have been possible to simulate using a small-domain LES model alone.

southeast Atlantic Ocean↗
DOE OSTI2024

Application of Modified Meshgraphnets for Subsurface Prediction during CO2 Sequestration

In the face of the increasingly dire consequences of anthropogenic climate change, capturing and storing carbon dioxide is paramount. However, several impediments exist to the safe and effective subsurface storage of CO2, such as cost of transport, identification of suitable sites for subsurface storage, and assessment of long-term risk from storage in subsurface aquifers. Accurate subsurface modeling is necessary to ensure that CO2 storage is both safe and effective. Still, such modeling has traditionally required either substantial time and computational power (numerical simulation) or a substantial amount of pre-existing data for training (machine learning models). Additionally, these models lack flexibility in dealing with both changes in discretization of the input data and generalizability beyond the data on which they are trained. In order to address these issues, this research applies graph neural networks (GNNs) to predict subsurface saturation and pressure during CO₂ injection in a model of the Illinois Basin-Decatur Project (IBDP). GNNs provide a flexible, intuitive method for representing and manipulating complex unstructured data, which is often found in many practical domain problems such as fluid flow and subsurface characterization. These unstructured grids are easily represented in GNNs by representing spatially-localized features such as permeability, porosity, saturation, and pressure as nodes in a graph and relationships between these properties as edges connecting these nodes. This research applies a specific GNN model called MeshGraphNets (MGN) to model the change in CO2 saturation and pressure over a 50-month time period (36 months of injection, 14 months post-injection). The MGN model leverages a message passing process that allows the network to learn both the spatial and temporal dynamics of this system simultaneously. Additionally, training on a limited dataset (64 realizations, 20 time points each) resulted in a high degree of accuracy in saturation prediction both within the same timeframe as the training (20 months, 0.039 average RMSE) and when projecting out to the end of injection (36 months, 0.053 average RMSE). Temporal predictions such as those generated by MGNs and other similar models are prone to accumulated error over time; in order to address this, a multi-step rollout (MSR) training process was applied to calculate training loss. This method mimics the forward prediction during inference by “rolling out” multiple time points in a single training step using the previous prediction as input to the MGN model. By calculating the loss several time steps forward from the current prediction, the model is forced to find a more stable state over time. Application of MSR to the MGN model resulted in an average 15% reduction in inference error over time during forward prediction. This study showcases the immense potential of GNNs as a game-changing methodology for predicting pressure and saturation evolution in CCS projects, ultimately paving the way for more sustainable and effective carbon storage solutions. Presentation prepared for the 2024 AiChE Annual Meeting, October 27 to November 1 2024, San Diego, CA.

Holcomb, Paul↗
NASA NTRS2000

Optimization of Car Body under Constraints of Noise, Vibration, and Harshness (NVH), and Crash

To be competitive on the today's market, cars have to be as light as possible while meeting the Noise, Vibration, and Harshness (NVH) requirements and conforming to Government-man dated crash survival regulations. The latter are difficult to meet because they involve very compute-intensive, nonlinear analysis, e.g., the code RADIOSS capable of simulation of the dynamics, and the geometrical and material nonlinearities of a thin-walled car structure in crash, would require over 12 days of elapsed time for a single design of a 390K elastic degrees of freedom model, if executed on a single processor of the state-of-the-art SGI Origin2000 computer. Of course, in optimization that crash analysis would have to be invoked many times. Needless to say, that has rendered such optimization intractable until now. The car finite element model is shown. The advent of computers that comprise large numbers of concurrently operating processors has created a new environment wherein the above optimization, and other engineering problems heretofore regarded as intractable may be solved. The procedure, shown, is a piecewise approximation based method and involves using a sensitivity based Taylor series approximation model for NVH and a polynomial response surface model for Crash. In that method the NVH constraints are evaluated using a finite element code (MSC/NASTRAN) that yields the constraint values and their derivatives with respect to design variables. The crash constraints are evaluated using the explicit code RADIOSS on the Origin 2000 operating on 256 processors simultaneously to generate data for a polynomial response surface in the design variable domain. The NVH constraints and their derivatives combined with the response surface for the crash constraints form an approximation to the system analysis (surrogate analysis) that enables a cycle of multidisciplinary optimization within move limits. In the inner loop, the NVH sensitivities are recomputed to update the NVH approximation model while keeping the Crash response surface constant. In every outer loop, the Crash response surface approximation is updated, including a gradual increase in the order of the response surface and the response surface extension in the direction of the search. In this optimization task, the NVH discipline has 30 design variables while the crash discipline has 20 design variables. A subset of these design variables (10) are common to both the NVH and crash disciplines. In order to construct a linear response surface for the Crash discipline constraints, a minimum of 21 design points would have to be analyzed using the RADIOSS code. On a single processor in Origin 2000 that amount of computing would require over 9 months! In this work, these runs were carried out concurrently on the Origin 2000 using multiple processors, ranging from 8 to 16, for each crash (RADIOSS) analysis. Another figure shows the wall time required for a single RADIOSS analysis using varying number of processors, as well as provides a comparison of 2 different common data placement procedures within the allotted memories for each analysis. The initial design is an infeasible design with NVH discipline Static Torsion constraint violations of over 10%. The final optimized design is a feasible design with a weight reduction of 15 kg compared to the initial design. This work demonstrates how advanced methodology for optimization combined with the technology of concurrent processing enables applications that until now were out of reach because of very long time-to-solution.

Kodiyalam, Srinivas↗
NASA NTRS2009

Relation Between Type II Bursts and CMEs Inferred from STEREO Observations

The inner coronagraph (COR1) of the Solar Terrestrial Relations Observatory (STEREO) mission has made it possible to observe CMEs in the spatial domain overlapping with that of the metric type II radio bursts. The type II bursts were associated with generally weak flares (mostly B and C class soft X-ray flares), but the CMEs were quite energetic. Using CME data for a set of type II bursts during the declining phase of solar cycle 23, we determine the CME height when the type II bursts start, thus giving an estimate of the heliocentric distance at which CME-driven shocks form. This distance has been determined to be approx. 1.5Rs (solar radii), which coincides with the distance at which the Alfven speed profile has a minimum value.We also use type II radio observations from STEREO/WAVES and Wind/WAVES observations to show that CMEs with moderate speed drive either weak shocks or no shock at all when they attain a height where the Alfven speed peaks (approx. 3Rs - 4Rs). Thus the shocks seem to be most efficient in accelerating electrons in the heliocentric distance range of 1.5Rs to 4Rs. By combining the radial variation of the CME speed in the inner corona (CME speed increase) and interplanetary medium (speed decrease) we were able to correctly account for the deviations from the universal drift-rate spectrum of type II bursts, thus confirming the close physical connection between type II bursts and CMEs. The average height (approx 1.5Rs) of STEREO CMEs at the time of type II bursts is smaller than that (2.2Rs) obtained for SOHO (Solar and Heliospheric Observatory) CMEs. We suggest that this may indicate, at least partly, the density reduction in the corona between the maximum and declining phases, so a given plasma level occurs closer to the Sun in the latter phase. In two cases, there was a diffuse shock-like feature ahead of the main body of the CME, indicating a standoff distance of 1Rs - 2Rs by the time the CME left the LASCO field of view.

Solar Energetic Particles↗
NASA NTRS2010

Relation Between Type II Bursts and CMEs Inferred from STEREO Observations

The inner coronagraph (COR1) of the Solar Terrestrial Relations Observatory (STEREO) mission has made it possible to observe coronal mass ejections (CMEs) a in the spatial domain overlapping with that of the metric type II radio bursts. The type II bursts were associated with generally weak flares (mostly B and C class soft X-ray flares), but the CMEs were quite energetic. Using CME data for a set of type II bursts during the declining phase of solar cycle 23, we determine the CME height when the type II bursts start, thus giving an estimate of the heliocentric distance at which CME-driven shocks form. This distance has been determined to be approximately 1.5Rs (solar radii), which coincides with the distance at which the Alfv?n speed profile has a minimum value. We also use type II radio observations from STEREO/WAVES and Wind/WAVES observations to show that CMEs with moderate speed drive either weak shocks or no shock at all when they attain a height where the Alfv?n speed peaks (?3Rs ? 4Rs). Thus the shocks seem to be most efficient in accelerating electrons in the heliocentric distance range of 1.5Rs to 4Rs. By combining the radial variation of the CME speed in the inner corona (CME speed increase) and interplanetary medium (speed decrease) we were able to correctly account for the deviations from the universal drift-rate spectrum of type II bursts, thus confirming the close physical connection between type II bursts and CMEs. The average height (approximately 1.5 Rs) of STEREO CMEs at the time of type II bursts is smaller than that (2.2 Rs) obtained for SOHO (Solar and Heliospheric Observatory) CMEs. We suggest that this may indicate, at least partly, the density reduction in the corona between the maximum and declining phases, so a given plasma level occurs closer to the Sun in the latter phase. In two cases, there was a diffuse shock-like feature ahead of the main body of the CME, indicating a standoff distance of 1Rs - 2Rs by the time the CME left the LASCO field of view.

Gopalswamy, N.↗
NASA NTRS2002

Skeletal Phenotype of Transgenic Mice Expressing the Beta1 Integrin Cytoplasmic Tail In Osteoblasts

To define the physiologic role of beta1 integrin in bone formation and mechanical loading, transgenic mice were generated by expressing the cytoplasmic tall and transmembrane domain of Beta1 integrin under the control of the osteocalcin promoter. In cultured cells, this truncated fragment of Beta1 can act as a dominant negative. Previously, the matrix of calvariae was shown to be abnormal in transgenic (TG) compared to wildtype (WT) mice. In this study, we analyzed appendicular bone in TG and WT, male and female mice at 14, 35, 63, 90 and 365 days old (n=8-12/gp). To assess beta1 integrin function in mechanical loading, a pilot study using hindlimb unloading by tail suspension was performed. 35d old TG and WT females were hindlimb unloaded for 4 wks (n=3-5). Body mass, bone mineral content, histomorphometric (distal femur) and biomechanical parameters were analyzed. Statistical significance (P less than.05) was defined by ANOVA using the Tukey-Kramer post-hoc test. We confirmed transgene expression by immunoprecipitating then immunoblotting bone lysates using an antibody against the beta1 tail. Body masses of TG mice at 63, 90 and 365d old were greater (16-25%) than WT. Some TG female mice at 365d appeared obese; mean abdominal fat mass was 415% greater in TG than WT mice. Tibiae were longer (5-7%) in TG than WT mice at 63 and 90d. Tibial mineral mass of 35d males was 7% lower in TG than WT mice, but at 63d was 21% higher. The % osteoblast surface in 35d TG mice was 20% higher than WT, and at 63d was 17% lower, while % osteoclast surface did not differ. In 365d mice, cancellous bone volume (125%) and endocortical mineral apposition rate (40%) were greater in TG than WT males but not females. In WT mice, hindlimb unloading caused a reduction in mineral mass of tibiae (-20%) and lumbar vertebrae (-22%) relative to normally loaded controls. Surprisingly, hindlimb unloading also caused a relative reduction (-13%) in humerus mass. The effects of hindlimb unloading on tibia and humerus mass were less obvious in TG than in WT mice. Since hindlimb unloading caused skeletal changes in both loaded and unloaded bones, systemic changes may contribute to bone responses observed using this animal model. In conclusion, transgene expression resulted in marked metabolic changes during growth and in the aged female. Our results demonstrate that expression of the Beta1 integrin cytoplasmic tail in vivo causes gender- and age-specific changes in select morphometric parameters, bone length, and bone mass.

Globus, R. K.↗
NASA NTRS2011

Performance Tests of Snow-Related Variables Over the Tibetan Plateau and Himalayas Using a New Version of NASA GEOS-5 Land Surface Model that Includes the Snow Darkening Effect

The snow darkening effect, i.e. the reduction of snow albedo, is caused by absorption of solar radiation by absorbing aerosols (dust, black carbon, and organic carbon) deposited on the snow surface. This process is probably important over Himalayan and Tibetan glaciers due to the transport of highly polluted Atmospheric Brown Cloud (ABC) from the Indo-Gangetic Plain (IGP). This effect has been incorporated into the NASA Goddard Earth Observing System model, version 5 (GEOS-5) atmospheric transport model. The Catchment land surface model (LSM) used in GEOS-5 considers 3 snow layers. Code was developed to track the mass concentration of aerosols in the three layers, taking into account such processes as the flushing of the compounds as liquid water percolates through the snowpack. In GEOS-5, aerosol emissions, transports, and depositions are well simulated in the Goddard Chemistry Aerosol Radiation and Transport (GO CART) module; we recently made the connection between GOCART and the GEOS-5 system fitted with the revised LSM. Preliminary simulations were performed with this new system in "replay" mode (i.e., with atmospheric dynamics guided by reanalysis) at 2x2.5 degree horizontal resolution, covering the period 1 November 2005 - 31 December 2009; we consider the final three years of simulation here. The three simulations used the following variants of the LSM: (1) the original Catchment LSM with a fixed fresh snowfall density of 150 kg m-3 ; (2) the LSM fitted with the new snow albedo code, used here without aerosol deposition but with changes in density formulation and melting water effect on snow specific surface area, (3) the LSM fitted with the new snow albedo code as same as (2) but with fixed aerosol deposition rates (computed from GOCART values averaged over the Tibetan Plateau domain [Ion.: 60-120E; lat.: 20-50N] during March-May 2008) applied to all grid points at every time step. For (2) and (3), the same setting on the fresh snowfall density as in (1) was used.

Yasunari, Tppei J.↗
NASA NTRS2017

Wavenumber-frequency Spectra of Pressure Fluctuations on a Generic Space Vehicle Measured via Unsteady Pressure-Sensitive Paint

Time histories of pressure fluctuations on a generic, hammerhead space vehicle model were measured using unsteady Pressure-Sensitive Paint (uPSP). The test was conducted in the 11-foot transonic wind tunnel of NASA Ames Research Center over a Mach number range of 0.6 M 1.2, and angles of attack of -4 4. The model was coated with a porous binder and PtTFPP-based porous polymer paint. An elaborate system of four high-speed cameras, and forty LED lamps was used for image acquisition. Various steps for image registration, reduction of shot noise, photogrammetry procedure to map images from the four cameras on a grid for the model, and finally a calibration procedure to convert the measured fluctuations in light intensity to fluctuating pressure, are discussed in the paper. The calibration process using a set of unsteady pressure sensors mounted on the model, was found to overcome some of the inherent problems of the fast response paint, such as rapid photo-degradation, non-linearity in pressure response, and significant temperature sensitivity. Comparison of spectra of pressure fluctuations between UPSP and pressure sensors demonstrated the ability of the paint to faithfully follow fluctuations up to 10 kHz, the maximum attempted. It was also found that the camera bit-depth and the illumination level limited the lowest measurable levels of pressure fluctuations to around 140dB. The large data set exposed various critical transonic flow physics not seen before, such as a coupling of the shock motion on the Payload Fairing (PF) with the separated flow region on the upper stage of the launch vehicle, and upstream convection of pressure fluctuation on PF at certain Mach numbers. The data also confirmed the expectation of a general lowering of the coefficient of pressure fluctuation with Mach number. The availability of the data set on a dense, regularly-spaced, surface grid allowed for the calculation of wavenumber-frequency (k-) spectra via straightforward applications of Fourier transform. The k- spectra were compared for the separated flow regions on the Second Stage, and the shock-boundary layer interactions on PF. The former showed self-similarity with Mach number while the latter was distinctly different, and confirmed the upstream propagation of pressure fluctuations. The k- spectra were dominated by the convected fluctuations; the acoustic domain was not discernable. These data, valuable for the vibro-acoustics analysis of aerospace vehicles, are believed to be the first obtained for the transonic flight regime, and pave the path for application on production models of aerospace vehicles.

buffet↗
NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗
NASA NTRS2019

Modelling above-ground biomass stock over Norway using national forest inventory data with ArcticDEM and Sentinel-2 data

Boreal forests constitute a large portion of the global forest area, yet they are undersampled through field surveys, and only a few remotely sensed data sources provide structural information wall-to-wall throughout the boreal domain. ArcticDEM is a collection of high-resolution (2 m) space-borne stereogrammetric digital surface models (DSM) covering the entire land area north of 60° of latitude. The free-availability of ArcticDEM data offers new possibilities for aboveground biomass mapping (AGB) across boreal forests, and thus it is necessary to evaluate the potential for these data to map AGB over alternative open-data sources (i.e., Sentinel-2). This study was performed over the entire land area of Norway north of 60° of latitude, and the Norwegian national forest inventory (NFI) was used as a source of field data composed of accurately geolocated field plots (n=7710) systematically distributed across the study area. Separate random forest models were fitted using NFI data, and corresponding remotely sensed data consisting of either: i) a canopy height model (ArcticCHM) obtained by subtracting a high-quality digital terrain model (DTM) from the ArcticDEM DSM height values, ii) Sentinel-2 (S2), or iii) a combination of the two (ArcticCHM+S2). Furthermore, we assessed the effect of the forest- and terrain-specific factors on the models’ predictive accuracy. The best model (,i.e., ArcticCHM+S2) explained nearly 60% of the variance of the training set, which translated in the largest accuracy in terms of root mean square error (RMSE=41.4 t/ha). This result highlights the synergy between 3D and multispectral data in AGB modelling. Furthermore, this study showed that despite the importance of ArcticCHM variables, the S2 model performed slightly better than ArcticCHM model. This finding highlights some of the limitations of ArcticDEM, which, despite the unprecedented spatial resolution, is highly heterogeneous due to the blending of multiple acquisitions across different years and seasons. We found that both forest- and terrain-specific characteristics affected the uncertainty of the ArcticCHM+S2 model and concluded that the combined use of ArcticCHM and Sentinel-2 represents a viable solution for AGB mapping across boreal forests. The synergy between the two data sources allowed for a reduction of the saturation effects typical of multispectral data while ensuring the spatial consistency in the output predictions due to the removal of artifacts and data voids present in ArcticCHM data. While the main contribution of this study is to provide the first evidence of the best-case-scenario (i.e., availability of accurate terrain models) that ArcticDEM data can provide for large-scale AGB modelling, it remains critically important for other studies to investigate how ArcticDEM may be used in areas where no DTMs are available as is the case for large portions of the boreal zone.

space-borne imagery↗
NASA NTRS2007

Primary Objective Grating Astronomical Telescope

It has been 370 years since a seventeenth century French mathematician, Mersenne, presciently sketched out an astronomical telescope based on dual parabolic reflectors. Since that time the concept of the primary objective has been virtually unchanged. Now a new class of astronomical telescope with a primary objective grating (POG) has been studied as an alternative. The POG competes with mirrors, in part, because diffraction gratings provide the very chromatic dispersion that mirrors defeat. The resulting telescope deals effectively with long-standing restrictions on multiple object spectroscopy (MOS). Other potential benefits include unprecedented apertures and collection areas. The new design also favors space deployment as a gossamer membrane. The inventor, Tom Ditto, first discovered that higher-order diffraction images contain hidden depth cues, for which he was granted a seminal range finding patent in 1987. Subsequently, he invented and patented 3D localizers, profilometers and microscopes using POGs. The POG telescope was placed in the public domain to expedite research. The function of a telescopes primary objective is to collect flux and to deliver images. Both functions dictate that size matters, and bigger is better. For that reason, there has been a steady push over the past century to ramp up the size of the primary mirror. However, for every doubling of mirror diameter, the elapsed time between initial effort and first light has also doubled. Meanwhile, costs escalated beyond the mirror alone, because larger instruments required larger enclosures and better pointing mechanisms. One key catalog of observation, spectrographic data, is far more difficult to amass than two-dimensional imagery. While the number of observable objects has increased with mirror size, the capacity to take spectra has not increased proportionately. In the best of circumstances, spectrograms are available for one per cent of the all objects surveyed. Spectroscopy was a historical afterthought introduced in the nineteenth century shortly after the invention of the diffraction grating and over a century after Newtons 1670 telescope. Spectroscopy is generally accomplished using a diffraction grating as the disperser in the secondary. The light being delivered to the spectrograph is first captured by a primary mirror which provides no chromatic magnification by itself. Sizeable spectrographs could not be deployed while diffraction gratings were rare commodities scribed using mechanical ruling engines that produced one grating line at a time. Today diffraction gratings are commonplace. Their recent availability is a product of both the invention of holography and the mass replication of surface microstructures. Holography permits all lines in a grating to be made simultaneously in a single photographic exposure. Holograms can then be reproduced by embossing processes. The improvement in replication is analogous to how Gutenberg changed the availability of books. The masters may be expensive, but the copies are not. Computer science is another technology that emerged in the second half of the twentieth century without which our proposed spectrographic instrument could not function due to the complexity of image processing required in data reduction. The employment of very large diffraction gratings as primary objectives for astronomical telescopes requires a novel

spectroscopy↗
DOE OSTI2025

Coal to Carbon Fiber (C2CF) Continuous Processing for High Value Composites (Final Report)

Coal tar is a condensed and recovered by-product of the coking of metallurgical coal for steel production. The heaviest fraction of distilled coal tar is an isotropic pitch largely used as a binder in the manufacturing of carbonaceous electrodes for primary aluminum smelting and in electric arc furnaces. Coal tar pitch offers high carbon yield upon carbonization. In this project, a process to convert the domestically sourced isotropic coal tar pitch, containing very low particulates (QI = 0.32 wt.%), to form flow-domain mesophase pitch amenable for melt spinning into precursor (green) fibers for carbon fiber, was developed and optimized. The final reproducible processing developed is reviewed in this report, along with several characterizations of the mesophase pitch. A final definition of the characteristics of a ‘spinnable’ mesophase pitch is presented. With this mesophase pitch, reproducible and stable multifilament melt spinning was developed, producing green fiber tows with filament diameters of approximately 20 m. The multifilament melt spinning was the most challenging aspect of the project and required the most effort. The best practices learned from this project for melt spinning are reviewed in this report. Once spun, the green fibers were oxidatively stabilized, carbonized and graphitized under inert gas atmosphere to form the final carbon fibers. Given the relatively high softening point of the mesophase pitch, no issues of interfilament fusion were observed during batch oxidation, and subsequent batch carbonization and graphitization went smoothly in all cases. An ~ 80 wt.% conversion of the green fiber mass to final carbon fiber was achieved. After graphitization, the carbon fibers showed high tensile moduli (most were ~ 600 GPa, or 87 Msi) consistent with commercially-available high-performance pitch-based carbon fiber. However, tensile strength and strain to failure were comparatively low. Further work to reduce defects in and on the fiber surfaces would increase these properties. SEM imaging of the graphitic fiber textures is presented herein. Rudimentary composites were fabricated from the carbon fibers and characterized showing similar modulus to baseline composites fabricated with commercial carbon fiber. Finally, a basic economic analysis was done showing the potential to increase the value of the isotropic coal tar pitch by up to 13.6 to 136 times based on a carbon fiber value of $\$$5/lb to $\$$50/lb, respectively. Moreover, the site case study suggests that the coal tar from the single integrated steel mill could supply production of up to 16 kt/yr of carbon fiber. Finally, a technological gap analysis was done which shed light on remaining technical challenges. These challenges included: recovery and utilization of condensates from the mesophase pitch processing, further advancing and increasing stability of the multifilament melt spinning processing, optimization of the oxidation processing, defect reduction for increased carbon fiber strength, and the development of a weaving process towards carbon fiber fabrics. To maximize the coal value chain, the primary objectives of this project were to (a) develop and scale efficient processing technology for producing melt-spinnable mesophase pitch from isotropic coal tar pitch, (b) clarify and simplify tedious continuous fiber processing technologies (particularly multifilament melt spinning of mesophase pitch) towards the efficient production of high performance carbon fiber, and (c) demonstrate and characterize representative composite parts derived from the final carbon fiber. Immense progress was made on all 3 objectives and is detailed in this report.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Microchannel-based Membrane-less Extraction of Li from Unconventional Lithium Sources & the Separation of REE

This final report provides an overview of the Project's entire duration, covering July 1, 2021 to December 31, 2023. It primarily focuses on the achievements, technological developments, and unique challenges the team faced while working on separating and extracting Lithium from produced waters. The project's primary aim was to create an integrated, high-throughput, membrane-less, and modular microfluidic platform that could extract Lithium from unconventional sources. We have successfully met all goals and milestones envisioned in the SOPO document. The most critical primary milestones, including the Go-No-Go milestone (refer to the Gantt chart in the Appendices), were successfully accomplished. We demonstrated phase separation (>90%) and extraction (>85%) performance in the MPSE using synthetic, and representative produced water composition feed at 50 ml/min total flow through MPSE 36. We have also performed a parametric study of the MPSE operations, beyond the scope of SOPO, exploring operating conditions of current and broader interest. The extended investigation of operational parameters is concurrent with our efforts to seek further development of the MPSE technology beyond the scope of the Project. Along these lines of development, we have made efforts to be responsive to DOE calls for technological developments of other types of resources (beyond PW) for the recovery of Critical Materials and higher TRL development (beyond TRL 4). During the work on this Project, we developed and implemented three innovative technical approaches that emerged from our efforts to successfully meet the Project milestones. The innovative & original technical approaches developed and implemented in this Project are now the contributions to process engineering that could be clearly credited to the Project. First, Convergent Design Approach is a comprehensive feedforward & feedback loop of four design phases: i) design for functionality, ii) design for manufacturing, iii) design for sustainability, and iv) design for market. Next was Process Intensification. A major aim of this Project was to create an innovative phase separation & extraction microscale-based technology for Li separation – thus the words microchannel-based in the Project title. A microscale-based technology is intrinsically in the center of the Process Intensification domain as defined by its unique principles. Therefore, Process Intensification was implicitly envisioned in the Project’s SOPO. Lastly, Time Scale Analysis is a novel tool for discovering the needs and directions of Process Intensification implementations in any process technology. This Project is fully credited for developing and implementing the three novel technical approaches mentioned above. These are general contributions to process engineering that emerged from this Project. Beyond the original SOPO scope, the OSU-U.Pitt research group utilized a Convergent Design methodology, integrating first-principles mathematical modeling with experimental validation on the Minimum Development Vehicle. By creating these Digital Twins, the team rapidly assessed manufacturing iterations to support TEA analysis. This framework further enabled the development of advanced Surface Modification Techniques, where hydrophobic and oleophobic coating strategies were optimized via Digital Twin tools and validated through rigorous 100-hour longevity testing. TEA Analysis: The closing efforts of this Project were focused on the TEA analysis. TEA analysis had two primary functions: i) enabling critical assessments of design variations withing 10 the Concurrent Design Approach, thus enabling evolution of the MPSE design to reach faster- better-cheaper alternatives; and ii) to create a bridge between the accomplishments of this Project and future projects of higher TRL, beyond TRL 6 level. It is important to note that the TEA model created in the Project stirred the technological solutions for the recovery of critical materials toward a vision of a very profitable modular plant that has unique zero-waste water discharge signature. More importantly, thanks to our experimental performance data and conservative assumptions, the TEA model predicts minimal technological and investment risks. Low cost of a modular unit of a nominal capacity of [1000 tons of Li 2 CO 3 /year] positions the MPSE based technology within the reach of community investors, thus offering a paradigm shift in the development of critical technologies. The project successfully navigated two primary challenges: solvent selection and manufacturing adaptation. Restricted by the SOPO to existing literature for lithium recovery, the team identified a critical need for a "material excellence program" to develop next-generation solvents, eventually concluding with a preliminary investigation into promising Ionic Liquids (ILs). Simultaneously, COVID-19 supply chain disruptions forced a pivot from traditional manufacturing to advanced additive methods at ATAMI-OSU. By transitioning from stainless steel to 3D-printed polymer substrates, the team achieved a transformative three-order-of- magnitude reduction in manufacturing costs and compressed prototyping timelines from several months to just two days. The MPSE technology offers significant energy, environmental, and economic advantages by overcoming the traditional bottlenecks of phase-separation hardware and contactor size. Unlike conventional mixer-settlers or membrane-based systems, MPSE operates without moving parts or fouling-prone membranes, achieving robust performance even with challenging, viscous, or particulate-heavy feeds. Key performance metrics include an energy intensity reduction of 5–50x (3–40 kJ/m 3 ) compared to incumbent technologies and a dramatic reduction of processing time to under 60 seconds, which drastically reduces the physical plant footprint. These technical efficiencies translate into superior economic outcomes; for a 100 t/year Li 2 CO 3 facility, implementing MPSE is projected to nearly halve contactor CAPEX (from $\$$6.08M to $\$$3.01M) and significantly increase the project's Net Present Value (NPV), derisking new investment and enabling distributed critical-mineral processing configurations. The commercialization of MPSE technology is being spearheaded by Vigsur Dynamics Inc., which has adopted a structured, parallel approach to technical and business development since its formation in January 2026. Following extensive customer discovery and engagement with the Oregon State University accelerator, Vigsur Dynamics is working to establish a business model that transitions from pilot demonstrations to modular hardware sales, ultimately aiming for a "build-own-operate" service strategy. Current technical milestones—including 100 hours of continuous operation, superior energy efficiency, and successful 6-unit modular scale-up— provide a foundation for this transition. Backed by ongoing IP licensing and a growing network of industrial and venture advisors, the company is actively de-risking the platform to replace conventional mixer-settler systems in the critical minerals market.

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

Improving the Quasi‐Biennial Oscillation via a Surrogate‐Accelerated Multi‐Objective Optimization

Accurate simulation of the quasi-biennial oscillation (QBO) is challenging due to uncertainties in representing convectively generated gravity waves. We develop an end-to-end uncertainty quantification workflow that calibrates these gravity wave processes in E3SM for a realistic QBO. Central to our approach is a domain knowledge-informed, compressed representation of high-dimensional spatio-temporal wind fields. By employing a parsimonious statistical model that learns the fundamental frequency from complex observations, we extract interpretable and physically meaningful quantities capturing key attributes. Building on this, we train a probabilistic surrogate model that approximates the fundamental characteristics of the QBO as functions of critical physics parameters governing gravity wave generation. Leveraging the Karhunen–Loève decomposition, our surrogate efficiently represents these characteristics as a set of orthogonal features, capturing cross-correlations among multiple physics quantities evaluated at different pressure levels and enabling rapid surrogate-based inference at a fraction of the computational cost of full-scale simulations. Finally, we analyze the inverse problem using a multi-objective approach. Our study reveals a tension between amplitude and period that constrains the QBO representation, precluding a single optimal solution. To navigate this, we quantify the bi-criteria trade-off and generate a set of Pareto optimal parameter values that balance the conflicting objectives. This integrated workflow improves the fidelity of QBO simulations and offers a versatile template for uncertainty quantification in complex geophysical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Non-intrusive reduced-order modeling for dynamical systems with spatially localized features

This work presents a non-intrusive reduced-order modeling framework for dynamical systems with spatially localized features characterized by slow singular value decay. The proposed approach builds upon two existing methodologies for reduced and full-order non-intrusive modeling, namely Operator Inference (OpInf) and sparse Full-Order Model (sFOM) inference. We decompose the domain into two complementary subdomains that exhibit fast and slow singular value decay. The dynamics of the subdomain exhibiting slow singular value decay are learned with sFOM while the dynamics with intrinsically low dimensionality on the complementary subdomain are learned with OpInf. The resulting, coupled OpInf-sFOM formulation leverages the computational efficiency of OpInf and the high resolution of sFOM, and thus enables fast non-intrusive predictions for conditions beyond those sampled in the training data set. A novel regularization technique with a closed-form solution based on the Gershgorin disk theorem is introduced to promote stable sFOM and OpInf models. We also provide a data-driven indicator for subdomain selection and ensure solution smoothness over the interface via a post-processing interpolation step. We evaluate the efficiency of the approach in terms of offline and online speedup through a quantitative, parametric computational cost analysis. We demonstrate the coupled OpInf-sFOM formulation for two test cases: a one-dimensional Burgers’ model for which accurate predictions beyond the span of the training snapshots are presented, and a two-dimensional parametric model for the Pine Island Glacier ice thickness dynamics, for which the OpInf-sFOM model achieves an average prediction error on the order of 1% with an online speedup factor of approximately 8$\times$ compared to the numerical simulation.

42 ENGINEERING↗