SEARCH · Search NASA

Results for “formulation study”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Tuning the 3D Printability and Thermomechanical Properties of Radiation Shields

Additive manufacturing, with its rapid advances in materials science, allows for researchers and companies to have the ability to create novel formulations and final parts that would have been difficult or near impossible to fabricate with traditional manufacturing methods. One such 3D printing technology, direct ink writing, is especially advantageous in fields requiring customizable parts with high amounts of functional fillers. Nuclear technology is a prime example of a field that necessitates new material design with regard to unique parts that also provide radiation shielding. Indeed, much effort has been focused on developing new rigid radiation shielding components, but DIW remains a less explored technology with a lot of potential for nuclear applications. In this study, DIW formulations that can behave as radiation shields were developed and were printed with varying amounts of porosity to tune the thermomechanical performance.

36 MATERIALS SCIENCE↗

NASA NTRS2005

Evaluation of a Shuttle Derived Vehicle (SDV) for Cargo Transportation

In this new era of space exploration, a host of launch vehicles are being examined for possible use in transporting cargo and crew to low Earth orbit and beyond. Launch vehicles derived from the Space Shuttle Program (SSP), known as Shuttle Derived Vehicles (SDVs), are prime candidates for heavy-lift duty because of their potential to minimize non-recurring costs and because the Shuttle can leverage off proven high-performance flight systems with established ground and flight support. To determine the merits of SDVs, a detailed evaluation was performed. This evaluation included a trade study and risk assessment of options based on performance, safety reliability, cost, operations, and evolution. The purpose of this paper is to explain the approach, processes, and tools used to evaluate launch vehicles for heavy lift cargo transportation. The process included defining the trade space, characterizing the concepts, analyzing the systems, and scoring the options. The process also included a review by subject experts from NASA and industry to compare past and recent study data and assess the risks. A set of technical performance measures (TPMs) was generated based on the study requirements and constraints. Tools such as INTROS and POST were used to calculate performance, FIRST was used for prediction of reliability, and other software packages, both commercial and NASA-owned, were applied to study the trade space. By following a clear process and using the right tools a thorough assessment was performed. An SDV can be classified as either a side-mount vehicle (SMV) or an in-line vehicle OLV). An SMV is a Space Shuttle where the Orbiter is replaced by a cargo carrier. An ILV is comprised of a modified Shuttle External Tank (ET) with engines mounted to the bottom and cargo mounted atop. For both families of vehicles, Solid Rocket Boosters (SRBs) are attached to the ET. The first derivate of Shuttle is defined as the vehicle with minimum changes necessary to transform the Space Shuttle into an SDV. Deltas from the first derivate were also formulated to study more SDV options. Examples of deltas include replacing the SRBs with larger and/or more SRBs, adding an upper stage, increasing the size of the ET, changing the engines, and modifying the elements. Challenges for SDV range from tailoring infrastructure to meeting the exploration schedule. Although SDV is based on the Space Shuttle, it still includes development risk for designing and building a Cargo Carrier. There are also performance challenges in that Shuttle is not optimized for cargo-only missions, but it is a robust system built on reusability. Balancing the strengths and weaknesses of the Shuttle to meet Lunar and Mars mission objectives provides the framework for an informative trade study. SDV was carefully analyzed and the results of the study provide invaluable data for use in the new exploration initiative.

Roman, Jose M.↗

DOE OSTI2022

Electrification of industrial boilers in the USA: potentials, challenges, and policy implications

Decarbonization of the industrial heat demand through electrification could contribute significantly to climate change mitigation efforts. In the US industry, thermal processes accounted for 75% of the total final energy demand in 2018, of which 17% is consumed by conventional industrial boilers (excluding cogeneration) for steam generation. Electric boilers have a small share in the US industrial steam generation due to several techno-economic reasons. Here, this study employs a bottom-up approach to investigate the sector-level and state-level techno-enviro-economic potentials of deploying electric boilers in the US industry in different timeframes. The results show that the technical potential energy and CO 2 savings by electrifying industrial boilers are 595 PJ and 200 MtCO 2 per year in 2050, respectively; however, these incur additional costs in each sector. Although there may be individual cost-effective opportunities for electrifying boilers in specific industrial sites, the overall costs are high in all industrial sectors and states due to the large disparity between electricity and combustion fuel prices. To overcome the highlighted techno-economic barriers, a comprehensive action plan for different stakeholders is also formulated. This study provides novel insights that should inform policymakers’ and executives’ decisions about the electrification of the current and future US industrial boiler systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS2000

Modeling of Nonacoustic Combustion Instability in Simulations of Hybrid Motor Tests

A transient model of a hybrid motor was formulated to study the cause and elimination of nonacoustic combustion instability. The transient model was used to simulate four key tests out of a series of seventeen hybrid motor tests conducted by Thiokol, Rocketdyne, and Martin Marietta at NASA Marshall Space Flight Center (MSFC). These tests were performed under the Hybrid Propulsion Technology for Launch Vehicle Boosters (HPTLVB) program. The first test resulted in stable combustion. The second test resulted in large-amplitude, 6.5-Hz chamber pressure oscillations that gradually damped away by the end of the test. The third test resulted in large-amplitude, 7.5-Hz chamber pressure oscillations that were sustained throughout the test. The seventh test resulted in elimination of combustion instability with the installation of an orifice immediately upstream of the injector. Formulation and implementation of the model are the scope of this presentation. The current model is an independent continuation of modeling presented previously by joint Thiokol-Rocketdyne collaborators Boardman, Hawkins, Wassom. and Claflin. The previous model simulated an unstable independent research and development (IR&D) hybrid motor test performed by Thiokol. There was very good agreement between the model and test data. Like the previous model, the current model was developed using Matrix-x simulation software. However, tests performed at MSFC under the HPTLVB program were actually simulated. ln the current model, the hybrid motor, consisting of the liquid oxygen (lox) injector, the multiport solid fuel grain, and nozzle, was simulated. The lox feedsystem, consisting of the tank, venturi. valve, and feed lines, was also simulated in the model. All components of the hybrid motor and lox feedsystem are treated by a lumped-parameter approach. Agreement between the results of the transient model and actual test data was very good. This agreement between simulated and actual test data indicated that the combustion instability in the hybrid motor was due to two causes: 1. a lox feed system of insufficient stiffness, and 2. a lox injector with an impedance (it pressure drop that was too low to provide damping against the feed system oscillations. Also, it was discovered that testing with a new grain of solid fuel sustained the combustion instability. However, testing with a used grain of solid fuel caused the combustion instability to gradually decay.

Rocker, M.↗

NASA NTRS1999

Simulation of Non-Acoustic Combustion Instability in a Hybrid Rocket Motor

A transient model of a hybrid motor was formulated to study the cause and elimination of non-acoustic combustion instability. The transient model was used to simulate four key tests out of a series of seventeen hybrid motor tests conducted by Thiokol, Rocketdyne and Martin Marietta at NASA/Marshall Space Flight Center (NASA/MSFC). These tests were performed under the Hybrid Propulsion Technology for Launch Vehicle Boosters (HPTLVB) program. The first test resulted in stable combustion. The second test resulted in large-amplitude, 6.5 Hz chamber pressure oscillations that gradually damped away by the end of the test. The third test resulted in large-amplitude, 7.5 Hz chamber pressure oscillations that were sustained throughout the test. The seventh test resulted in the elimination of combustion instability with the installation of an orifice immediately upstream of the injector. The formulation and implementation of the model are the scope of this presentation. The current model is an independent continuation of modeling presented previously by joint Thiokol-Rocketdyne collaborators Boardman, Hawkins, Wassom, and Claflin. The previous model simulated an unstable IR&D hybrid motor test performed by Thiokol. There was very good agreement between the model and the test data. Like the previous model, the current model was developed using Matrix-x simulation software. However, the tests performed at NASA/MSFC under the HPTLVB program were actually simulated. In the current model, the hybrid motor consisting of the liquid oxygen (LOX) injector, the multi-port solid fuel grain and the nozzle was simulated. Also, simulated in the model was the LOX feed system consisting of the tank, venturi, valve and feed lines. All components of the hybrid motor and LOX feed system are treated by a lumped-parameter approach. Agreement between the results of the transient model and the actual test data was very good. This agreement between simulated and actual test data indicated that the combustion instability in the hybrid motor was due to two causes. The first cause was a LOX feed system of insufficient stiffness. The second cause was a LOX injector with an impedance or pressure drop that was too low to provide damping against the feed system oscillations. Also, it was discovered that testing with a new grain of solid fuel sustained the combustion instability. However, testing with a used grain of solid fuel caused the combustion instability to gradually decay.

Rocker, Marvin↗

NASA NTRS1999

Simulation of Non-Acoustic Combustion Instability in a Hybrid Rocket Motor

A transient model of a hybrid motor was formulated to study the cause and elimination of non-acoustic combustion instability. The transient model was used to simulate four key tests out of a series of seventeen hybrid motor tests conducted by Thiokol, Rocketdyne and Martin Marietta at NASA/Marshall Space Flight Center (NASAIMSFC). These tests were performed under the Hybrid Propulsion Technology for Launch Vehicle Boosters (HPTLVB) program. The first test resulted in stable combustion. The second test resulted in large-amplitude, 6.5 Hz chamber pressure oscillations that gradually damped away by the end of the test. The third test resulted in large-amplitude, 7.5 Hz chamber pressure oscillations that were sustained throughout the test. The seventh test resulted in the elimination of combustion instability with the installation of an orifice immediately upstream of the injector. The formulation and implementation of the model are the scope of this presentation. The current model is an independent continuation of modeling presented previously by joint Thiokol-Rocketdyne collaborators Boardman, Hawkins, Wassom, and Claflin. The previous model simulated an unstable IR&D hybrid motor test performed by Thiokol. There was very good agreement between the model and the test data. Like the previous model, the current model was developed using Matrix-x simulation software. However, the tests performed at NASA/MSFC under the HPTLVB program were actually simulated. In the current model, the hybrid motor consisting of the liquid oxygen (LOX) injector, the multi-port solid fuel grain and the nozzle was simulated. Also, simulated in the model was the LOX feed system consisting of the tank, venturi, valve and feed lines. All components of the hybrid motor and LOX feed system are treated by a lumped-parameter approach. Agreement between the results of the transient model and the actual test data was very good. This agreement between simulated and actual test data indicated that the combustion instability in the hybrid motor was due to two causes. The first cause was a LOX feed system of insufficient stiffness. The second cause was a LOX injector with an impedance or pressure drop that was too low to provide damping against the feed system oscillations. Also, it was discovered that testing with a new grain of solid fuel sustained the combustion instability. However, testing with a used grain of solid fuel caused the combustion instability to gradually decay.

Rocker, Marvin↗

DOE OSTI2021

Energy flexibility of residential buildings: A systematic review of characterization and quantification methods and applications

With building electric demand becoming increasingly dynamic, and a growing percentage of intermittent renewable power generation from solar photovoltaics and wind turbines, the power grid is facing increasing challenge to manage the real time balance between the supply and demand. With advancements in smart sensing and metering, smart appliances, electric vehicles, and energy storage technologies, demand side management of residential buildings can help the grid to improve stability by optimizing flexible loads. This paper reviews recent studies on residential building demand side management, with a focus on characterization and quantification of energy flexibility covering various types of flexible loads, metrics, methods, and applications. The reviewed studies showed four levels of applications: building level (45%), district or community level (29%), system level (19%), and building sector level (7%). Shifting loads is the dominant flexibility type in 60% of applications, followed by shedding (19%), generation (16%), and modulating (6%). Depending on the technology and application scope, flexible operations have a wide range of performance, with peak power reductions of 1%~65%, energy savings up to 60%, operational cost reduction of 1%~48%, and greenhouse gas emission reductions of up to29%. More than half (51%) of the studies employed control strategies to achieve flexibility; among those 72% used optimal controls, while 28% used rule-based controls. About 58% of the studies used mathematical formulation to quantify energy flexibility. Most studies were based on simulation, while less than 15% of the studies had measurements from experiments or field tests. The review reveals research opportunities to address significant gaps in the existing literature: (1) establishing a common definition and performance metrics for energy flexibility of buildings that are technology and application agnostic, (2) developing an ontology to standardize representation of flexibility resources for interoperability, (3) integrating occupant impacts into the quantification and optimization of energy flexibility, and (4) developing requirements and credits of energy flexibility in building energy codes and standards. Findings from the review can inform future research and development of energy flexible buildings which are essential to a reliable and resilient power grid.

Li, H↗

DOE OSTI2025

Networked Microgrid Topology Reconfiguration to Promote Fairness in Proactive Load Shedding

Increasing occurrences of natural disasters and grid emergency events consistently challenge the safe and reliable operations of power systems. During such emergency situations, system operators may proactively shed load to mitigate risks. However, uncoordinated implementation of load shedding may disrupt electricity supply and even lead to cascading failures. Meanwhile, it is crucial to address potential biases affecting different customers when executing load shedding. This paper addresses the dynamic topology reconfiguration problem for networked microgrids with distributed energy resources under emergency conditions. Specifically, we propose a novel rolling-horizon optimization model that integrates fairness-aware constraints into the networked microgrid topology reconfiguration. Unlike existing approaches that focus solely on efficiency or apply fairness considerations in static settings, our method explicitly incorporates temporal fairness constraints to restrict repeated or excessive load curtailment for load blocks. Moreover, the fairness-aware constraints are specifically developed for the context of dynamic networked microgrid topology reconfiguration, and are designed to be convex or amenable to linear reformulations, which offers a more tractable alternative to traditional models with non-convex formulations. Numerical studies on a modified IEEE 13-bus system and a larger-sized SMART-DS networked microgrid system demonstrate the performance of the proposed algorithm towards more fairness-aware networked microgrid topology reconfiguration decision-making.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Model development for flash boiling spray and validations with isooctane, hexane, ethanol and their binary mixtures

The objective of this study is to develop a comprehensive model for spray atomization under the consideration of flash boiling, which includes nucleation, bubble growth, and droplet breakup. Although flash boiling has been previously studied, there is still no well-established analytical and numerical model of their mechanisms, or general guidelines on the implementation of this atomization strategy. In this study, detailed formulations are elaborated for the unified spray atomization model under both non-flash and flash boiling conditions. Next, the model was implemented into CFD software and 3-D simulations were performed for model validations. Two fuel blends are tested: ethanol-isooctane blends and hexane-isooctane blends. Simulation results were compared with the experimental measurements conducted using two methodologies, DBI (Diffuse Backlight Imaging) system for spray images captured by a high-speed camera, and a PDA (Phase Doppler Anemometry) system for droplet size measurement during the spray injection process. It was found that the current flash boiling model can accurately capture the spray characteristics under different environment conditions, such as the spray plume interactions and liquid penetration. The averaged Sauter mean diameter (SMD) calculations match the measurement sufficiently well. The spray structures for both single- and multi-component fuels during the injection process are well represented by current model.

42 ENGINEERING↗

DOE OSTI2023

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

Here, extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov subspace approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. The current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semi-empirical quantum-mechanical treatment. Applied to molecular dynamics, simulation of 840 carbon atoms, one integration time step executes in 4 s on a single Nvidia RTX A6000 GPU.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Temperature Evolution of the Activation Barriers Leads to Meyer−Neldel Rules for Structural Relaxation and Transport in Polymers

Understanding activation barriers controlling structural relaxation in glass-forming liquids, molecular transport, and ionic conductivity in amorphous polymers is a grand challenge of fundamental scientific and materials engineering interest across disciplines. Over decades, intriguing but puzzling empirical correlations between the elementary time scale of activated barrier crossing and the apparent Arrhenius activation energy, the so-called Meyer−Neldel (MN) rules, have been discovered in diverse liquids and glasses. Here, in this study, we formulate and successfully apply a new experimental analysis and an explicitly dynamical theoretical framework which provides an understanding of the origin, validity, and failure of such correlations, that bridge and unify the three fields of structural relaxation, molecular transport, and ionic conductivity in liquids and quenched glasses. Distinct quasi-universal laws are predicted in equilibrated liquids and nonequilibrium glasses, consistent with experiments. Our analysis reveals that even if the relaxation appears Arrhenius over a limited temperature range, the physical activation barrier is generally temperature-dependent in polymeric systems even below glass transition temperature. In addition, we show that the approximate validity of classical MN rules hinges on a linear temperature dependence of this barrier and the temperature range probed in experiments. Our findings are relevant for controlling the activation barrier in functional soft polymeric materials relevant to molecular separations, barrier coatings, and charge transport, and also provide new constraints on the theoretical understanding of the mechanism underlying slow activated dynamics in glass-forming condensed matter.

Meyer-Neldel rules↗

DOE OSTI2021

Learning nonlinear operators via DeepONet based on the universal approximation theorem of operators

It is widely known that neural networks (NNs) are universal approximators of continuous functions. However, a less known but powerful result is that a NN with a single hidden layer can accurately approximate any nonlinear continuous operator. This universal approximation theorem of operators is suggestive of the structure and potential of deep neural networks (DNNs) in learning continuous operators or complex systems from streams of scattered data. Here, in this work, we thus extend this theorem to DNNs. We design a new network with small generalization error, the deep operator network (DeepONet), which consists of a DNN for encoding the discrete input function space (branch net) and another DNN for encoding the domain of the output functions (trunk net). We demonstrate that DeepONet can learn various explicit operators, such as integrals and fractional Laplacians, as well as implicit operators that represent deterministic and stochastic differential equations. We study different formulations of the input function space and its effect on the generalization error for 16 different diverse applications.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Stochastic Microgrid Scheduling With Chance‐Constrained Resilience Consideration

Traditionally, it is assumed that microgrids transition seamlessly from grid‐connected operation to islanded mode in the event of sudden main grid outages. In reality, the islanding process, especially unintentional islanding, is rarely seamless. Instead, it is subject to voltage and frequency fluctuations caused by the instantaneous disconnection of the point of common coupling (PCC) switch, variations in loads and renewable generation output and even the protection tripping of distributed energy resources (DERs). To mitigate these fluctuations and facilitate a smooth islanding process, we propose a stochastic microgrid scheduling model that incorporates chance‐constrained resilience measures. Specifically, the resilience measure is defined as the probability of successful islanding (PSI), that is, the probability that a microgrid can mitigate the generation‐demand imbalance caused by the disconnection of the PCC switch, variations in load and renewable generation and DER tripping. This measure is modelled using chance constraints. Unlike existing reliability and resilience indices, which typically neglect the possibility of microgrid/DER failure under extreme events and assume their survival while primarily focussing on reducing impact duration or magnitude, the proposed PSI‐based framework explicitly addresses microgrid and DER survival during the islanding transition. The formulated nonlinear chance constraints are approximated using a multiinterval approach and equivalently represented as a mixed‐integer linear programming (MILP) formulation. Case study results validate the proposed method, showing that the PSI estimation error is reduced to less than 8%, compared to approximately 28% with existing methods. Various sensitivity analyses on the DER tripping rate and PSI settings were performed to validate the robustness of the proposed method. In particular, the necessity of accounting for DER tripping in the PSI calculation was demonstrated.

chance constrained optimization↗

DOE OSTI2022

Deblurring galaxy images with Tikhonov regularization on magnitude domain

Abstract We propose a regularization-based deblurring method that works efficiently for galaxy images. The spatial resolution of a ground-based telescope is generally limited by seeing conditions and is much worse than space-based telescopes. This circumstance has generated considerable research interest in the restoration of spatial resolution. Since image deblurring is a typical inverse problem and often ill-posed, solutions tend to be unstable. To obtain a stable solution, much research has adopted regularization-based methods for image deblurring, but the regularization term is not necessarily appropriate for galaxy images. Although galaxies have an exponential or Sérsic profile, the conventional regularization assumes the image profiles to behave linearly in space. The significant deviation between the assumption and real situations leads to blurring of the images and smoothing out the detailed structures. Clearly, regularization on logarithmic domain, i.e., magnitude domain, should provide a more appropriate assumption, which we explore in this study. We formulate a problem of deblurring galaxy images by an objective function with a Tikhonov regularization term on a magnitude domain. We introduce an iterative algorithm minimizing the objective function with a primal–dual splitting method. We investigate the feasibility of the proposed method using simulation and observation images. In the simulation, we blur galaxy images with a realistic point spread function and add both Gaussian and Poisson noise. For the evaluation with the observed images, we use galaxy images taken by the Subaru HSC-SSP. Both of these evaluations show that our method successfully recovers the spatial resolution of the deblurred images and significantly outperforms the conventional methods. The code is publicly available from the GitHub 〈https://github.com/kzmurata-astro/PSFdeconv_amag〉.

Murata, Kazumi↗

DOE OSTI2022

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Symmetric Kondo Lattice States in Doped Strained Twisted Bilayer Graphene

We use the topological heavy fermion (THF) model and its Kondo lattice (KL) formulation to study the possibility of a symmetric Kondo (SK) state in twisted bilayer graphene. Via a large-N approximation, we find a SK state in the KL model at fillings v = 0,±1,±2 where a KL model can be constructed. In the SK state, all symmetries are preserved and the local moments are Kondo screened by the conduction electrons. At the mean-field level of the THF model at v = 0,±1,±2,±3 we also find a similar symmetric state that is adiabatically connected to the symmetric Kondo state. We study the stability of the symmetric state by comparing its energy with the ordered (symmetry-breaking) states found in [H. Hu et al., Phys. Rev. Lett. 131, 026502 (2023)., Z.-D. Song and B. A. Bernevig, Phys. Rev. Lett. 129, 047601 (2022).] and find the ordered states to have lower energy at v = 0,±1,±2. However, moving away from integer fillings by doping the light bands, our mean-field calculations find the energy difference between the ordered state and the symmetric state to be reduced, which suggests the loss of ordering and a tendency toward Kondo screening. In order to include many-body effects beyond the mean-field approximation, we also performed dynamical mean-field theory calculations on the THF model in the nonordered phase. The spin susceptibility follows a Curie behavior at v = 0,±1,±2 down to ~2 K where the onset of screening of the local moment becomes visible. This hints to very low Kondo temperatures at these fillings, in agreement with the outcome of our mean-field calculations. At noninteger filling v = ±0.5,±0.8,±1.2 dynamical mean-field theory shows deviations from a 1/T susceptibility at much higher temperatures, suggesting a more effective screening of local moments with doping. Lastly, we study the effect of a C 3⁢z -rotational-symmetry-breaking strain via mean-field approaches and find that a symmetric phase (that only breaks C 3⁢z symmetry) can be stabilized at sufficiently large strain at v = 0,±1,±2. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at noninteger fillings or by applying strain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

On the Role of Airborne Ice Nucleating Particles in Primary and Secondary Ice Formation Processes in Convective Midlatitude Clouds

Abstract Formation of ice in mixed-phase clouds may proceed via primary nucleation and secondary ice production (SIP). Primary nucleation involves the presence of ice nucleating particles (INPs), while SIP follows primary nucleation, resulting in ice crystal number concentrations N ice exceeding the number of INPs. The ice formation pathways in young congestus clouds are not well understood nor are the environmental conditions favorable for the onset of SIP. Coincident airborne measurements of INPs and N ice collected during the Secondary Production of Ice in Cumulus Experiment (SPICULE) campaign over the U.S. central Great Plains are reported for young congestus clouds. Number concentrations of INPs and N ice at temperatures between −10° and −20°C were used to categorize whether SIP was active in each analyzed cloud pass. Observational data suggested that fragmentation of freezing droplets (FFDs) was occurring and the most likely SIP mechanism for ice enhancement during the early stages of the sampled clouds, where higher cloud-base temperatures and stronger updrafts produced conditions favorable for the onset of FFD. Simple model simulations of a congestus cloud were used to compare predictions of N ice from newly implemented SIP mechanisms (including FFD) alongside those from the existing Hallett–Mossop (HM) mechanism. Although predicted N ice , when all SIP mechanisms were active, was in good agreement with observations, the dominant SIP mechanism in the model was predicted to be HM. The rate of heterogeneous freezing in raindrops was likely inhibiting realistic FFD rates in the simulated cloud. These discrepancies between the observations and simulations underscore the need for extended laboratory studies to formulate improved model representations of SIP in convective clouds.

Patnaude, Ryan J. [Department of Atmospheric Scien↗

DOE OSTI2024

Advancing attenuation estimation through integration of the Hessian in multiparameter viscoacoustic full-waveform inversion

Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.

Geochemistry & Geophysics↗