SEARCH · Search NASA

Results for “Model initialization”

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 973 records · Page 54

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Uncertainty Evaluation of Computational Model Used to Support the Integrated Powerhead Demonstration Project

NASA and the U.S. Air Force are working on a joint project to develop a new hydrogen-fueled, full-flow, staged combustion rocket engine. The initial testing and modeling work for the Integrated Powerhead Demonstrator (IPD) project is being performed by NASA Marshall and Stennis Space Centers. A key factor in the testing of this engine is the ability to predict and measure the transient fluid flow during engine start and shutdown phases of operation. A model built by NASA Marshall in the ROCket Engine Transient Simulation (ROCETS) program is used to predict transient engine fluid flows. The model is initially calibrated to data from previous tests on the Stennis E1 test stand. The model is then used to predict the next run. Data from this run can then be used to recalibrate the model providing a tool to guide the test program in incremental steps to reduce the risk to the prototype engine. In this paper, they define this type of model as a calibrated model. This paper proposes a method to estimate the uncertainty of a model calibrated to a set of experimental test data. The method is similar to that used in the calibration of experiment instrumentation. For the IPD example used in this paper, the model uncertainty is determined for both LOX and LH flow rates using previous data. The successful use of this model is then demonstrated to predict another similar test run within the uncertainty bounds. The paper summarizes the uncertainty methodology when a model is continually recalibrated with new test data. The methodology is general and can be applied to other calibrated models.

Steele, W. G.↗

NASA NTRS2023

Harnessing the Digital Transformation for Development of Electrified Aircraft Propulsion Control Systems

Hybrid electric aircraft propulsion is an emerging technology that presents a variety of potential benefits along with technical integration challenges. Developing these new propulsion architectures with their complex control systems, and ultimately proving their benefit, is a multistep process. This process includes concept development and analysis, dynamic simulation, hardware-in-the-loop testing, full-scale testing, and so on. This effort is being revolutionized and indeed enabled by new digital tools that support increasing the technology readiness level throughout the maturation process. As part of this Digital Transformation, NASA has developed a suite of publicly available digital tools that facilitate the path from concept to implementation. This paper describes the NASA-developed tools and puts them in the context of control system development for hybrid electric aircraft propulsion. The three MATLAB®-based software packages are the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS), the Electrical Modeling and Thermal Analysis Toolbox (EMTAT), and the Thermal Systems Analysis Toolbox (TSAT). These tools are interactive, complementary, and compatible with each other. T-MATS is a modular thermodynamic modeling framework designed for creating custom component level models of jet engines. EMTAT is a modeling framework used to simulate a variety of power electronic devices, using both physics-based and power flow calculations. TSAT is a framework for modeling and analysis of dynamic heat transfer. These packages all consist of graphical, drag-and-drop, parameterizable building blocks representing various components of the system to be modeled, e.g., compressors, turbines, motors, energy storage devices, etc. They are designed to enable the user to model and simulate the end-to-end dynamic operation of a hybrid electric gas turbine engine powertrain at the timescale of the turbomachinery, capturing mechanical, electrical, and thermal interactions. This paper demonstrates through multiple examples how these tools have been used successfully in a variety of applications, including several of the early stages of hybrid electric gas turbine engine propulsion system development, from the initial system modeling to real-time interactive pilot-in-the-loop simulation to physical hardware-in-the-loop testing, each step bringing the technology closer to fruition.

hybrid electric propulsion↗

DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2017

Monte Carlo Computational Modeling of Atomic Oxygen Interactions

Computational modeling of the erosion of polymers caused by atomic oxygen in low Earth orbit (LEO) is useful for determining areas of concern for spacecraft environment durability. Successful modeling requires that the characteristics of the environment such as atomic oxygen energy distribution, flux, and angular distribution be properly represented in the model. Thus whether the atomic oxygen is arriving normal to or inclined to a surface and whether it arrives in a consistent direction or is sweeping across the surface such as in the case of polymeric solar array blankets is important to determine durability. When atomic oxygen impacts a polymer surface it can react removing a certain volume per incident atom (called the erosion yield), recombine, or be ejected as an active oxygen atom to potentially either react with other polymer atoms or exit into space. Scattered atoms can also have a lower energy as a result of partial or total thermal accommodation. Many solutions to polymer durability in LEO involve protective thin films of metal oxides such as SiO2 to prevent atomic oxygen erosion. Such protective films also have their own interaction characteristics. A Monte Carlo computational model has been developed which takes into account the various types of atomic oxygen arrival and how it reacts with a representative polymer (polyimide Kapton H) and how it reacts at defect sites in an oxide protective coating, such as SiO2 on that polymer. Although this model was initially intended to determine atomic oxygen erosion behavior at defect sites for the International Space Station solar arrays, it has been used to predict atomic oxygen erosion or oxidation behavior on many other spacecraft components including erosion of polymeric joints, durability of solar array blanket box covers, and scattering of atomic oxygen into telescopes and microwave cavities where oxidation of critical component surfaces can take place. The computational model is a two dimensional model which has the capability to tune the interactions of how the atomic oxygen reacts, scatters, or recombines on polymer or nonreactive surfaces. In addition to the specification of atomic oxygen arrival details, a total of 15 atomic oxygen interaction parameters have been identified as necessary to properly simulate observed interactions and resulting polymer erosion that have been observed in LEO. The tuning of the Monte Carlo model has been accomplished by adjusting interaction parameters so the erosion patterns produced by the model match those from several actual LEO space experiments. Surface texturing in LEO can also be predicted by the model. Such comparison of space tests with ground laboratory experiments have enabled confidence in ground laboratory lifetime prediction of protected polymers. Results of Monte Carlo tuning, examples of surface texturing and undercutting erosion prediction, and several examples of how the model can be used to predict other LEO and Mars orbital space results are presented.

Monte Carlo Modeling↗

NASA NTRS2017

Monte Carlo Computational Modeling of Atomic Oxygen Interactions

Atomic Oxygen↗

DOE OSTI2024

Dynamic Simulation Modeling and Control of a Desiccant Assisted Direct-expansion Air Handling Unit

Desirable built environments demand simultaneous regulation of thermal comfort and indoor air quality (IAQ) with energy-efficient operation of heating, ventilation and air conditioning (HV AC) systems, which involves controls of temperature, humidity and airborne contaminants simultaneously. This paper presents the efforts of dynamic modeling and initial development control strategy for a desiccant-assisted multi-functional air handling unit (AHU) coupled with a direct-expansion rooftop unit (RTU) system, which aims to achieve multiple functions for indoor environment conditioning with energy efficient control. The RTU-AHU system includes a desiccant wheel for dehumidification and a conceptual direct air capture (DAC) filtering device for CO2 regulation. A Modelica-based dynamic model is developed for this conceptual system, and a simple decentralized control strategy is designed, which combines a differential-enthalpy based AHU return-air ratio control, a demand-controlled ventilation, and supply-air temperature humidity control via the RTU and DW controls. The proposed control method is evaluated with the Modelica simulation model for a selected set of scenarios

Pan, Chao↗

DOE OSTI2018

Monte Carlo Study of the Effect of Initial States and Structural Defects on Nonequilibrium Critical Behavior of the Three-Dimensional Ising Model

The effect of various initial magnetizations m{sub 0} and structural defects the nonequilibrium critical behavior of the three-dimensional Ising model is numerically studied. Based on an analysis of the time dependence of the magnetization and the two-time dependence of the autocorrelation function and dynamic susceptibility, the significant effect of initial states on relaxation magnetizations and aging effects characterized by anomalous relaxation inhibition and correlation in the system with increasing waiting time was revealed. The fluctuation–dissipation theorem violation was studied, and the values of the limit fluctuation–dissipation ratio (FDR) are calculated. It is shown that two universality subclasses can be distinguished in the nonequilibrium critical behavior of the three-dimensional Ising model with random initial magnetization m{sub 0} These subclasses correspond to the system evolution from the high-temperature (m{sub 0} = 0) and low-temperature (m{sub 0} = 1) initial states with limit FDRs characteristic of these states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2004

Sensor Web Dynamic Measurement Techniques and Adaptive Observing Strategies

Sensor Web observing systems may have the potential to significantly improve our ability to monitor, understand, and predict the evolution of rapidly evolving, transient, or variable environmental features and events. This improvement will come about by integrating novel data collection techniques, new or improved instruments, emerging communications technologies and protocols, sensor mark-up languages, and interoperable planning and scheduling systems. In contrast to today's observing systems, "event-driven" sensor webs will synthesize real- or near-real time measurements and information from other platforms and then react by reconfiguring the platforms and instruments to invoke new measurement modes and adaptive observation strategies. Similarly, "model-driven" sensor webs will utilize environmental prediction models to initiate targeted sensor measurements or to use a new observing strategy. The sensor web concept contrasts with today's data collection techniques and observing system operations concepts where independent measurements are made by remote sensing and in situ platforms that do not share, and therefore cannot act upon, potentially useful complementary sensor measurement data and platform state information. This presentation describes NASA's view of event-driven and model-driven Sensor Webs and highlights several research and development activities at the Goddard Space Flight Center.

Talabac, Stephen J.↗

NASA NTRS2018

Acceleration Modes and Transitions in Pulsed Plasma Accelerators

Pulsed plasma accelerators typically operate by storing energy in a capacitor bank and then discharging this energy through a gas, ionizing and accelerating it through the Lorentz body force. Two plasma accelerator types employing this general scheme have typically been studied: the gas-fed pulsed plasma thruster and the quasi-steady magnetoplasmadynamic (MPD) accelerator. The gas-fed pulsed plasma accelerator is generally represented as a completely transient device discharging in approximately 1-10 microseconds. When the capacitor bank is discharged through the gas, a current sheet forms at the breech of the thruster and propagates forward under a j (current density) by B (magnetic field) body force, entraining propellant it encounters. This process is sometimes referred to as detonation-mode acceleration because the current sheet representation approximates that of a strong shock propagating through the gas. Acceleration of the initial current sheet ceases when either the current sheet reaches the end of the device and is ejected or when the current in the circuit reverses, striking a new current sheet at the breech and depriving the initial sheet of additional acceleration. In the quasi-steady MPD accelerator, the pulse is lengthened to approximately 1 millisecond or longer and maintained at an approximately constant level during discharge. The time over which the transient phenomena experienced during startup typically occur is short relative to the overall discharge time, which is now long enough for the plasma to assume a relatively steady-state configuration. The ionized gas flows through a stationary current channel in a manner that is sometimes referred to as the deflagration-mode of operation. The plasma experiences electromagnetic acceleration as it flows through the current channel towards the exit of the device. A device that had a short pulse length but appeared to operate in a plasma acceleration regime different from the gas-fed pulsed plasma accelerators was developed by Cheng, et al. The Coaxial High ENerGy (CHENG) thruster operated on the 10-microseconds timescales of pulsed plasma thrusters, but claimed high thrust density, high efficiency and low electrode erosion rates, which are more consistent with the deflagration mode of acceleration. Separate work on gas-fed pulsed plasma thrusters (PPTs) by Ziemer, et al. identified two separate regimes of performance. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency (ratio of jet kinetic energy to input electrical energy) as a function of mass bit. In the second regime at very low mass bits (termed Mode II), the efficiency increased with decreasing mass bit. Work by Poehlmann et al. and by Sitaraman and Raja sought to understand the performance of the CHENG thruster and the Mode I / Mode II performance in PPTs by modeling the acceleration using the Hugoniot Relation, with the detonation and deflagration modes representing two distinct sets of solutions to the relevant conservation laws. These works studied the proposal that, depending upon the values of the various controllable parameters, the accelerator would operate in either the detonation or deflagration mode. In the present work, we propose a variation on the explanation for the differences in performance between the various pulsed plasma accelerators. Instead of treating the accelerator as if it were only operating in one mode or the other during a pulse, we model the initial stage of the discharge in all cases as an accelerating current sheet (detonation mode). If the current sheet reaches the exit of the accelerator before the discharge is completed, the acceleration mode transitions to the deflagration mode type found in the quasi-steady MPD thrusters. This modeling method is used to demonstrate that standard gas-fed pulsed plasma accelerators, the CHENG thruster, and the quasi-steady MPD accelerator are variations of the same device, with the overall acceleration of the plasma depending upon the behavior of the plasma discharge during initial transient phase and the relative lengths of the detonation and deflagration modes of operation.

Polzin, Kurt A.↗

DOE OSTI2020

Weak Lensing And Cosmology In The Dark Energy Survey

In this PhD Thesis, we present studies which aim at understanding and testing the fundamental assumptions of the $\Lambda$CDM cosmological model. We initially show that weak gravitational lensing can be derived as a consequence of General Relativity, and that it can be applied to astronomical data from the Dark Energy Survey (DES) in order to constrain the amplitude of matter fluctuations in the universe. This constraint can be used as an end-to-end cosmological test when it is compared to inferences from the Cosmic Microwave Background. We further extend the cosmological model to account for non-standard parameters which can point to novel physics, and use DES data to assess their statistical significance. We finally show that the Cold Dark Matter (CDM) paradigm can also be fundamentally tested with observational data using morphological signatures of disk galaxies in cluster environments if CDM is self-interacting.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Domain Aware Deep-learning Algorithms Integrated with Scientific-computing Technologies (DADAIST)

This technical report summarized the contribution of the DADAIST project funded by the Data Model Convergence Initiative via the Laboratory Directed Research and Development (LDRD) investments at Pacific Northwest National Laboratory (PNNL). Specifically, we report the development of the NeuroMANCER (Neural Modules with Adaptive Nonlinear Constraints and Efficient Regularizations), a new open-source Scientific Machine Learning library for formulating and solving parametric constrained optimization problems, physics-informed system identification, and parametric optimal control problems. NeuroMANCER is using differentiable programming to combine modern data-driven models and optimization modeling language into a coherent algorithmic and software framework. NeuroMANCER is a Pytorch-based framework and adopts much of its philosophy focused on research and development, rapid prototyping, and streamlined deployment. Strong emphasis is given to extensibility, interoperability with the PyTorch ecosystem, and quick adaptability to custom domain problems. Neuromancer repository contains a comprehensive library of differentiable modules, including custom activation functions, matrix factorizations, deep learning architectures, neural differential equations, differential equation solvers, implicit layers such as iterative solvers, high-level API for symbolic expressions, API for modeling and control of dynamical systems, and extensive set of tutorial code examples in the form of python scripts and jupyter notebooks.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Enhancing segmentation fairness through curriculum learning and progressive loss: a centralized and federated perspective on radiograph analysis

Bias in medical image segmentation can lead to unequal performance across demographic subgroups, raising concerns about fairness and reliability in clinical AI systems. While deep learning models have achieved high segmentation accuracy, ensuring equitable performance across race and gender remains a significant challenge, particularly in privacy-sensitive healthcare environments. This study investigates fairness-aware medical image segmentation for hip and knee radiographs using deep learning models evaluated in both centralized and Federated Learning (FL) settings. We introduce Curriculum Learning (CL) strategies and Progressive Loss (PL) functions to regulate sample difficulty during training. In addition, we propose two novel fairness-oriented federated learning algorithms, Federated Intersection over Union (FedIoU) and Federated Intersection over Union with Outlier Analysis (FedIoUoutlier). Experiments are conducted using multiple segmentation backbones and simulated multi-site data partitions derived from the Osteoarthritis Initiative dataset. Model performance is evaluated using Intersection over Union (IoU), IoU standard deviation, Skewed Error Ratio (SER), and Min-Max Disparity across race and gender subgroups. Statistical significance was verified using paired t-tests to compare per-sample IoU performance against baseline configurations. Across both hip and knee segmentation tasks, curriculum learning and progressive loss strategies consistently improved segmentation accuracy and reduced demographic performance disparities in centralized training. In federated settings, fairness-aware aggregation further enhanced performance. Notably, FedIoUoutlier combined with balanced curriculum learning and tiered progressive loss achieved the highest mean IoU while yielding the lowest SER and Min-Max Disparity, indicating improved fairness without sacrificing accuracy. In several configurations, federated models matched or exceeded the performance of optimized centralized models, with statistically significant improvements in per-sample IoU over baseline configurations. The results demonstrate that structured training strategies and fairness-aware federated aggregation can jointly improve accuracy, stability, and demographic fairness in medical image segmentation. By integrating curriculum learning, progressive loss, and novel FL algorithms, this work provides a practical pathway toward equitable and privacy-preserving AI systems for medical imaging.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2015

IPRT Polarized Radiative Transfer Model Intercomparison Project-Phase A

The polarization state of electromagnetic radiation scattered by atmospheric particles such as aerosols, cloud droplets, or ice crystals contains much more information about the optical and microphysical properties than the total intensity alone. For this reason an increasing number of polarimetric observations are performed from space, from the ground and from aircraft. Polarized radiative transfer models are required to interpret and analyse these measurements and to develop retrieval algorithms exploiting polarimetric observations. In the last years a large number of new codes have been developed, mostly for specific applications. Benchmark results are available for specific cases, but not for more sophisticated scenarios including polarized surface reflection and multi-layer atmospheres. The International Polarized Radiative Transfer (IPRT) working group of the International Radiation Commission (IRC) has initiated a model intercomparison project in order to fill this gap. This paper presents the results of the first phase A of the IPRT project which includes ten test cases, from simple setups with only one layer and Rayleigh scattering to rather sophisticated setups with a cloud embedded in a standard atmosphere above an ocean surface. All scenarios in the first phase A of the intercomparison project are for a one-dimensional plane-parallel model geometry. The commonly established benchmark results are available at the IPRT website

radiative transfer↗

NASA NTRSDate not supplied

Quantum Hardware-Enabled Molecular Dynamics via Transfer Learning

The ability to perform ab initio molecular dynamics simulations using potential energy surfaces provided by quantum computers would open the door to virtually exact dynamics for a variety of chemical and biochemical systems, with impacts on catalysis and biophysics. Nonetheless, performing molecular dynamics on surfaces produced by quantum hardware has been hampered by the noisy energies typically produced by quantum computers and challenges associated with computing gradients and scaling to large systems interest. A recent set of advances in machine learning, known as transfer learning, provides a new path forward for molecular dynamics simulations on quantum hardware. Transfer learning offers a workaround, where one first trains models on larger, less accurate classical datasets and then refines them on smaller, more accurate quantum datasets. We explore this approach by training machine learning models to predict a molecule's potential energy based on its geometric structure using Behler-Parrinello neural networks. When successfully trained, the model enables energy gradient predictions necessary for dynamic simulations. To reduce the quantum resources needed, the model is initially trained with data derived from classical density functional theory and subsequently refined with a smaller dataset obtained from a variational quantum eigensolver optimization of the unitary coupled cluster ansatz. We show that this approach significantly reduces the size of the needed quantum training dataset while capturing the high accuracies needed within quantum chemistry simulations. The success of this two-step training method opens more opportunities to apply machine learning models on quantum data, a significant stride towards efficient quantum-classical hybrid computational models.

quantum computing↗

DOE OSTI2017

Modeling of the chemistry in oxidation flow reactors with high initial NO

Oxidation flow reactors (OFRs) are increasingly employed in atmospheric chemistry research because of their high efficiency of OH radical production from low-pressure Hg lamp emissions at both 185 and 254 nm (OFR185) or 254 nm only (OFR254). OFRs have been thought to be limited to studying low-NO chemistry (in which peroxy radicals (RO 2 ) react preferentially with HO 2 ) because NO is very rapidly oxidized by the high concentrations of O 3 , HO 2 , and OH in OFRs. However, many groups are performing experiments by aging combustion exhaust with high NO levels or adding NO in the hopes of simulating high-NO chemistry (in which RO 2 + NO dominates). This work systematically explores the chemistry in OFRs with high initial NO. Using box modeling, we investigate the interconversion of N-containing species and the uncertainties due to kinetic parameters. Simple initial injection of NO in OFR185 can result in more RO 2 reacted with NO than with HO 2 and minor non-tropospheric photolysis, but only under a very narrow set of conditions (high water mixing ratio, low UV intensity, low external OH reactivity (OHR ext ), and initial NO concentration (NO in ) of tens to hundreds of ppb) that account for a very small fraction of the input parameter space. These conditions are generally far away from experimental conditions of published OFR studies with high initial NO. In particular, studies of aerosol formation from vehicle emissions in OFRs often used OHR ext and NO in several orders of magnitude higher. Due to extremely high OHR ext and NO in , some studies may have resulted in substantial non-tropospheric photolysis, strong delay to RO 2 chemistry due to peroxynitrate formation, VOC reactions with NO 3 dominating over those with OH, and faster reactions of OH–aromatic adducts with NO 2 than those with O 2 , all of which are irrelevant to ambient VOC photooxidation chemistry. Some of the negative effects are the worst for alkene and aromatic precursors. To avoid undesired chemistry, vehicle emissions generally need to be diluted by a factor of >100 before being injected into an OFR. However, sufficiently diluted vehicle emissions generally do not lead to high-NO chemistry in OFRs but are rather dominated by the low-NO RO 2 + HO 2 pathway. To ensure high-NO conditions without substantial atmospherically irrelevant chemistry in a more controlled fashion, new techniques are needed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2010

A Combined Observational and Modeling Approach to Study Modern Dust Transport from the Patagonia Desert to East Antarctica

The understanding of present atmospheric transport processes from Southern Hemisphere (SH) landmasses to Antarctica can improve the interpretation of stratigraphic data in Antarctic ice cores. In addition, long range transport can deliver key nutrients normally not available to marine ecosystems in the Southern Ocean and may trigger or enhance primary productivity. However, there is a dearth of observational based studies of dust transport in the SH. This work aims to improve current understanding of dust transport in the SH by showing a characterization of two dust events originating in the Patagonia desert (south end of South America). The approach is based on a combined and complementary use of satellite retrievals (detectors MISR, MODIS, GLAS ,POLDER, OMI,), transport model simulation (HYSPLIT) and surface observations near the sources and aerosol measurements in Antarctica (Neumayer and Concordia sites). Satellite imagery and visibility observations confirm dust emission in a stretch of dry lakes along the coast of the Tierra del Fuego (TdF) island (approx.54deg S) and from the shores of the Colihue Huapi lake in Central Patagonia (approx.46deg S) in February 2005. Model simulations initialized by these observations reproduce the timing of an observed increase in dust concentration at the Concordia Station and some of the observed increases in atmospheric aerosol absorption (here used as a dust proxy) in the Neumayer station. The TdF sources were the largest contributors of dust at both sites. The transit times from TdF to the Neumayer and Concordia sites are 6-7 and 9-10 days respectively. Lidar observations and model outputs coincide in placing most of the dust cloud in the boundary layer and suggest significant de- position over the ocean immediately downwind. Boundary layer dust was detected as far as 1800 km from the source and approx.800 km north of the South Georgia Island over the central sub-Antarctic Atlantic Ocean. Although the analysis suggests the presence of dust at approx.1500 km SW of South Africa five days after, the limited capabilities of existing satellite platforms to differentiate between aerosol types do not permit a definitive conclusion. In addition, the model simulations show dust lifting to the free troposphere as it travels south but it could not be confirmed by the satellite observations due to cloudiness. This work demonstrates that complementary information from existing transport models, satellite and surface data can yield a consistent picture of the dust transport from the Patagonia desert to Antarctica. It also illustrates the limitation of using any of these approaches individually to characterize the transport of dust in a heavily cloudy area.

Gasso, S.↗

NASA NTRS1982

Finite-element analysis of initiation, stable crack growth and instability using a crack-tip-opening displacement criterion

An elastic-plastic (incremental and small strain) finite element analysis was used with a crack growth criterion to study crack initiation, stable crack growth, and instability under monotonic loading to failure of metallic materials. The crack growth criterion was a critical crack-tip-opening displacement (CTOD) at a specified distance from the crack tip, or equivalently, a critical crack-tip-opening angle (CTOA). Whenever the CTOD (or CTOA) equaled or exceeded a critical value, the crack was assumed to grow. Single values of critical CTOD were used in the analysis to model crack initiation, stable crack growth, and instability for 7075-T651 and 2024-T351 aluminum alloy compact specimens. Calculated and experimentally measured CTOD values at initiation agreed well for both aluminum alloys. These critical CTOD values were also used to predict failure loads on center-crack tension specimens and a specially-designed three-hole-crack tension specimen made of the two aluminum alloys and of 304 stainless steel. All specimens were 12.7 mm thick. Predicted failure loads for 7075-T651 aluminum alloy and 304 stainless steel specimens were generally within + or - 15 percent of experimental failure loads, whereas the predicted failure loads for 2024-T351 aluminum alloy specimens were generally within + or - 5 percent of the experimental loads.

Newman, J. C., Jr.↗

DOE OSTI2023

Joint velocity and density reconstruction of the Universe with nonlinear differentiable forward modeling

Reconstructing the initial conditions of the Universe from late-time observations has the potential to optimally extract cosmological information. Due to the high dimensionality of the parameter space, a differentiable forward model is needed for convergence, and recent advances have made it possible to perform reconstruction with nonlinear models based on galaxy (or halo) positions. In addition to positions, future surveys will provide measurements of galaxies' peculiar velocities through the kinematic Sunyaev-Zel'dovich effect (kSZ), type Ia supernovae, the fundamental plane relation, and the Tully-Fisher relation. Here we develop the formalism for including halo velocities, in addition to halo positions, to enhance the reconstruction of the initial conditions. We show that using velocity information can significantly improve the reconstruction accuracy compared to using only the halo density field. Here we study this improvement as a function of shot noise, velocity measurement noise, and angle to the line of sight. We also show how halo velocity data can be used to improve the reconstruction of the final nonlinear matter overdensity and velocity fields. We have built our pipeline into the differentiable Particle-Mesh FlowPM package, paving the way to perform field-level cosmological inference with joint velocity and density reconstruction. This is especially useful given the increased ability to measure peculiar velocities in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗