SEARCH · Search NASA

Results for “Spatial Generative Model”

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 1,675 records · Page 93

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Modeling non-equilibrium discharge and validating transient plasma characteristics at above-atmospheric pressure

Non-equilibrium plasma generated from positive-pulsed nanosecond electrical discharges into desiccated air is simulated in this paper using a multi-dimensional, multi-physics plasma solver. A pin-to-pin electrode configuration is used with a fixed 5.2 mm gap spacing. Peak pulse voltages range between 10.2 and 22.5 kV. Care is taken to match the exact electrode profile from the experiments, and adjust the electron collision frequency so that breakdown limits closely match those from corresponding experimental results. The optimized numerical simulations predict qualitative streamer structure that is in close agreement with experimental observations. Quantitative measurements of atomic oxygen at the anode tip and qualitative estimates of streamer gas heating are closely matched by simulations. The model results are used to provide insight into the spatial and temporal development of the transient plasma. Furthermore, the work performed in this paper delivers a numerical tool that can be extremely useful to link the post-discharge plasma properties to low-temperature plasma ignition mechanisms that are of great interest for the automotive industry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Large-scale effects on the regulation of tropical sea surface temperature

The dominant terms in the surface energy budget of the tropical oceans are absorption of solar radiation and evaporative cooling. If it is assumed that relative humidity in the boundary layer remains constant, evaporative cooling will increase rapidly with sea surface temperature (SST) because of the strong temperature dependence of saturation water vapor pressure. The resulting stabilization of SST provided by evaporative cooling is sufficient to overcome positive feedback contributed by the decrease of surface net longwave cooling with increasing SST. Evaporative cooling is sensitive to small changes in boundary-layer relative humidity. Large and negative shortwave cloud forcing in the regions of highest SST are supported by the moisture convergence associated with largescale circulations. In the descending portions of these circulations the shortwave cloud forcing is suppressed. When the effect of these circulations is taken into account by spatial averaging, the area-averaged cloud forcing shows no sensitivity to area-averaged SST changes associated with the 1987 warming event in the tropical Pacific. While the shortwave cloud forcing is large and important in the convective regions, the importance of its role in regulating the average temperature of the tropics and in modulating temperature gradients within the tropics is less clear. A heuristic model of SST is used to illustrate the possible role of large-scale atmospheric circulations on SST in the tropics and the coupling between SST gradients and mean tropical SST. The intensity of large-scale circulations responds sensitivity to SST gradients and affects the mean tropical SST by supplying dry air to the planetary boundary layer. Large SST gradients generate vigorous circulations that increase evaporation and reduce the mean SST.

Hartmann, Dennis L.↗

NASA NTRS2011

4D Dynamic Required Navigation Performance Final Report

New advanced four dimensional trajectory (4DT) procedures under consideration for the Next Generation Air Transportation System (NextGen) require an aircraft to precisely navigate relative to a moving reference such as another aircraft. Examples are Self-Separation for enroute operations and Interval Management for in-trail and merging operations. The current construct of Required Navigation Performance (RNP), defined for fixed-reference-frame navigation, is not sufficiently specified to be applicable to defining performance levels of such air-to-air procedures. An extension of RNP to air-to-air navigation would enable these advanced procedures to be implemented with a specified level of performance. The objective of this research effort was to propose new 4D Dynamic RNP constructs that account for the dynamic spatial and temporal nature of Interval Management and Self-Separation, develop mathematical models of the Dynamic RNP constructs, "Required Self-Separation Performance" and "Required Interval Management Performance," and to analyze the performance characteristics of these air-to-air procedures using the newly developed models. This final report summarizes the activities led by Raytheon, in collaboration with GE Aviation and SAIC, and presents the results from this research effort to expand the RNP concept to a dynamic 4D frame of reference.

Finkelsztein, Daniel M.↗

NASA NTRS2017

Very High Resolution Tree Cover Mapping for Continental United States using Deep Convolutional Neural Networks

Uncertainties in input land cover estimates contribute to a significant bias in modeled above ground biomass (AGB) and carbon estimates from satellite-derived data. The resolution of most currently used passive remote sensing products is not sufficient to capture tree canopy cover of less than ca. 10-20 percent, limiting their utility to estimate canopy cover and AGB for trees outside of forest land. In our study, we created a first of its kind Continental United States (CONUS) tree cover map at a spatial resolution of 1-m for the 2010-2012 epoch using the USDA NAIP imagery to address the present uncertainties in AGB estimates. The process involves different tasks including data acquisition ingestion to pre-processing and running a state-of-art encoder-decoder based deep convolutional neural network (CNN) algorithm for automatically generating a tree non-tree map for almost a quarter million scenes. The entire processing chain including generation of the largest open source existing aerial satellite image training database was performed at the NEX supercomputing and storage facility. We believe the resulting forest cover product will substantially contribute to filling the gaps in ongoing carbon and ecological monitoring research and help quantifying the errors and uncertainties in derived products.

High Resolution↗

NASA NTRSDate not supplied

All-sky Retrieval of Atmospheric Temperature, Water Vapor, Clouds, Trace Gases, and Surface Properties from Operational Hyperspectral IR Sounders

Operational IR sounders such AIRS, CrIS, and IASI provide high quality hyperspectral measurements for weather and climate applications. We will describe a new all-sky Single Field-of-view Sounder Atmospheric Product (SiFSAP). The uniqueness of this product is that it uses all available channels from hyperspectral sounders and the optimal estimation retrieval is done at a single FOV spatial resolution. The SiFSAP includes atmospheric temperature, water vapor, clouds, trace gases, surface skin, and surface emissivity and will be produced operationally at NASA GES DISC. We will describe the core component of the SiFSAP algorithm, which is the Principal Component-based Radiative Transfer Model (PCRTM), and will show example applications of the SiFSAP product for various atmospheric weather and dynamics studies. We also describe a new Climate Fingerprinting Sounder Product (ClimFiSP), which is derived from spatiotemporally averaged level-1 hyperspectral radiances directly. The ClimFiSP algorithm uses consistent radiative kernels and a robust spectral fingerprinting method. It provides fast and accurate data fusion products from multiple satellite sensors. We have applied this method to both AIRS and CrIS (on SNPP and on NOAA 20) data and generated two decades climate data records for atmospheric temperature, water vapor, cloud, trace gases, and surface skin temperature. The ClimFiSP are being transitioned to NASA data centers for routine generations level-3 products.

pcrtm↗

NASA NTRS2021

Magnetic Reconnection in a Sheared Magnetic Flux Tube: Slippage Versus Tearing

The process of magnetic reconnection in a flux tube can occur in a time-stationary fashion as slippage reconnection or in a time-dependent manner based on, for example, the tearing instability. However, it is not well known under which conditions a system can sustain slippage reconnection. Likewise, it is unclear whether systems can exhibit slippage reconnection and time-dependent reconnection simultaneously. In order to investigate these questions, we employ a set of 3D magnetohydrodynamic simulation. Using these simulations, we model a twisted flux tube with a spatially localized resistive region in the center of the simulation box, and an applied driver on one end plane of the simulation box. As a result, the Poynting flux injected at the boundary propagates to the resistive region and dissipates there. If the driver is nearly circular, we find a predominance of slippage reconnection, which decouples the applied driver entirely from one half of the flux tube. Second, a tearing-like instability occurs with slippage reconnection in the resistive region when the resistivity is strong enough and the cross section of the flux tube, as defined by the shape of the applied velocity perturbation, is elliptical enough. In this case we find the surprising feature that magnetic field perturbation generated by the tearing-like instability radiates away from the resistive region to influence the entire flux tube shape.

Hidetaka Kuniyoshi↗

DOE OSTI2017

Effects of layered crust on the coseismic slip inversion and related CFF variations: Hints from the 2012 Emilia Romagna earthquake

The 2012 Emilia Romagna (Italy) seismic sequence has been extensively studied given the occurrence of two mainshocks, both temporally and spatially close to each other. The recent literature accounts for several fault models, obtained with different inversion methods and different datasets. Several authors investigated the possibility that the second event was triggered by the first mainshock with elusive results. Here, we consider all the available InSAR and GPS datasets and two planar fault geometries, which are based on both seismological and geological constraints. We account for a layered, elastic half-space hosting the dislocation and compare the slip distribution resulting from the inversion and the related changes in Coulomb Failure Function (CFF) obtained with both a homogeneous and layered half-space. Finally, we focus on the interaction between the two main events, discriminating the contributions of coseismic and early postseismic slip of the mainshock on the generation of the second event and discuss the spatio-temporal distribution of the seismic sequence. When accounting for both InSAR and GPS geodetic data we are able to reproduce a detailed coseismic slip distribution for the two mainshocks that is in accordance with the overall aftershock seismicity distribution. Additonally, we see that an elastic medium with depth dependent rigidity better accounts for the lack of the shallow seismicity, amplifying, with respect to the homogeneous case, the mechanical interaction of the two mainshocks.

58 GEOSCIENCES↗

DOE OSTI2023

Progenitors and explosion properties of supernova remnants hosting central compact objects: II. A global systematic study with a comparison to nucleosynthesis models

ABSTRACT Core-collapse explosions of massive stars leave behind neutron stars, with a known diversity that includes the ‘Central Compact Objects’ (CCOs). Typified by the neutron star discovered near the centre of the Cas A supernova remnant (SNR), CCOs have been observed to shine only in X-rays. To address their supernova progenitors, we perform a systematic study of SNRs that contain a CCO and display X-ray emission from their shock-heated ejecta. We make use of X-ray data primarily using the Chandra X-ray observatory, complemented with XMM–Newton. This study uses a systematic approach to the analysis of each SNR aimed at addressing the supernova progenitor as well as the explosion properties (energy and ambient density). After fitting for the ejecta abundances estimated from a spatially resolved spectroscopic study, we compare the data to six nucleosynthesis models making predictions on supernova ejecta yields in core-collapse explosions. We find that the explosion models commonly used by the astrophysics community do not match the ejecta yields for any of the SNRs, suggesting additional physics, for example multidimensional explosion models or updated progenitor structures, are required. Overall we find low-mass (≤25 solar masses) progenitors among the massive stars population and low-energy explosions (<1051 ergs). We discuss degeneracies in our model fitting, particularly how altering the explosion energy affects the estimate of the progenitor mass. Our systematic study highlights the need for improving on the theoretical models for nucleosynthesis predictions as well as for sensitive, high-resolution spectroscopy observations to be acquired with next generation X-ray missions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Detecting Important Drivers of Gridded Population Modeling With Machine Learning

High-resolution population datasets have been lever-aged across a broad swath of domains, such as climate change, public policy, humanitarian aid, and rescue operations, among others. Machine learning methods were adopted to generate high-resolution or gridded population estimates by using various geospatial input features such as buildings, roads, and nighttime lights. In this study, we evaluate the importance of population features using Random Forest models across three levels of analysis, utilizing permutation measures. Our research aims to address key questions to enhance our understanding of high-resolution population modeling, such as: Are certain features globally (10 countries collectively) more important than others? Do optimal features vary by country? Within each country, do feature importance differ across administrative units? What similarities exist in feature importance at the global, country, and administrative unit levels? To answer these questions, we leverage the Kneedle algorithm to automate the selection of optimum features. We find that there are patterns displayed by features across spatial boundaries, evidenced by the same feature being the most important indicator of population across 7 of the 10 countries modeled. Our findings indicate that while important features may vary across geographies, certain features consistently hold greater importance than others agnostic of geography.

Lebakula, Viswadeep [ORNL] (ORCID:0000000152935914↗

NASA NTRS2000

A Preview of AMSR: Airborne C-band Microwave Radiometer (ACMR) Observations from SGP99

Although L-band is generally considered ideal for passive microwave sensing of soil moisture, near-future satellite observing systems such as Advanced Mechanically Scanned Radiometer (AMSR) will provide C-band data for several years before any L-band data might become available. The Southern Great Plains'99 (SGP99) Experiment was designed to generate C-band observations suitable for testing and refinement of AMSR-era soil moisture retrieval algorithms. C-band data collected using the Airborne C-band Microwave Radiometer (ACMR), a new high-accuracy NASA/GSFC instrument, clearly demonstrated a strong response to a 9-day drydown event as well as to differences between the northern (cooler & wetter) and southern (warmer & dryer) areas covered by the P-3 flights. For example, the H-polarized brightness temperatures observed during the first three days of the drydown increased up to 50 K in the northern areas. These observations represent a preview of what we can expect from AMSR, albeit at 3-km spatial resolution vs. approximately 60 km for AMSR. Initial results of soil-vegetation microwave modeling will also be presented to estimate the relative contributions of soil physical temperature, canopy physical temperature, soil moisture, and canopy moisture. Significant radio-frequency interference (RFI) was evident during the experiment, and amelioration strategies will be discussed. The net effect of RFI (an upward bias in brightnesses) when averaged over an AMSR footprint is expected to be more subtle.

Kim, Edward↗

NASA NTRS2003

Double Layers in Expanding Plasmas and Their Relevance to the Auroral Plasma Processes

When a dense plasma consisting of a cold and a sufficiently warm electron population expands, a rarefaction shock forms [Bezzerides et al., 1978]. In the expansion of the polar wind in the magnetosphere, it has been previously shown that when a sufficiently warm electron population also exists, in addition to the usual cold ionospheric one, a discontinuity forms in the electrostatic potential distribution along the magnetic field lines [Barakat and Schunk, 1984]. Despite the lack of spatial resolution and the assumption of quasi-neutrality in the polar wind models, such discontinuities have been called double layers (DLs). Recently similar discontinuities have been invoked to partly explain the auroral acceleration of electrons and ions in the upward current region [Ergun et al., 2000]. By means of one-dimensional Vlasov simulations of expanding plasmas, for the first time we make here the connection between (1) the rarefaction shocks, (2) the discontinuities in the potential distributions, and (3) DLs. We show that when plasmas expand from opposite directions into a deep density cavity with a potential drop across it and when the plasma on the high-potential side contains hot and cold electron populations, the temporal evolution of the potential and the plasma distribution generates evolving multiple double layers with an ,extended density cavity between them. One of the DLs is the rarefaction-shock (RFS) and it forms by the reflections of the cold electrons coming from the high-potential side; it supports a part of the potential drop approximately determined by the hot electron temperature. The other DLs evolve from charge separations arising either from reflection of ions coming from the low-potential side or stemming from plasma instabilities; they support the rest of the potential drop. The instabilities forming these additional double layers involve electron-ion (e-i) Buneman or ion-ion (i-i) two-stream interactions. The electron-electron two-stream interactions on the high-potential side of the RFS generate electron-acoustic waves, which evolve into electron phase-space holes. The ion population originating from the low-potential side and trapped by the RFS is energized by the e-i and i-i instabilities and it eventually precipitates into the high-potential plasma along with an electron beam. Applications of these findings to the auroral plasma physics are discussed.

Singh, Nagendra↗

DOE OSTI2023

Stochastic Models, Indices & Optimization Algorithms for Pricing & Hedging Reliability Risks in Modern Power Grids: Data Plan - Princeton

We collected and cleaned the synthetic grid data produced by NREL for the Texas and New York synthetic grids. We developed a high dimensional joint stochastic model for load at the zone level, and solar and wind power productions at the asset level, capturing the spatial and temporal dependencies between all the variables, and demonstrated how such a model could be fitted to historical data. We designed and implemented a simulation engine which can produce Monte Carlo scenarios for the hourly day-ahead values of load, and solar and wind power productions at the spatial and temporal resolutions of the historical data used to fit the model. Finally we developed an open-source Python package which can, from an input grid model, efficiently use forecasts and large numbers of Monte Carlo scenarios to provide unit commitment and economic dispatch for each of these scenarios. The high dimensional stochastic model and the subsequent Monte Carlo simulation engine were implemented in the package PGscen and the corresponding UC and ED optimization programs in the package Vatic.

14 SOLAR ENERGY↗

DOE OSTI2017

Thermal fluctuations of ferroelectric nanodomains in a ferroelectric-dielectric PbTiO 3 / SrTiO 3 superlattice

Ferroelectric-dielectric superlattices consisting of alternating layers of ferroelectric PbTiO 3 and dielectric SrTiO 3 exhibit a disordered striped nanodomain pattern, with characteristic length scales of 6 nm for the domain periodicity and 30 nm for the in-plane coherence of the domain pattern. Spatial disorder in the domain pattern gives rise to coherent hard x-ray scattering patterns exhibiting intensity speckles. We show here using variable-temperature Bragg-geometry x-ray photon correlation spectroscopy that x-ray scattering patterns from the disordered domains exhibit a continuous temporal decorrelation due to spontaneous domain fluctuations. The temporal decorrelation can be described using a compressed exponential function, consistent with what has been observed in other systems with arrested dynamics. The fluctuation speeds up at higher temperatures and the thermal activation energy estimated from the Arrhenius model is 0.35±0.21 eV. As a result, the magnitude of the energy barrier implies that the complicated energy landscape of the domain structures is induced by pinning mechanisms and domain patterns fluctuate via the generation and annihilation of topological defects similar to soft materials such as block copolymers.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Energy Conversion and Electron Acceleration and Transport in 3D Simulations of Solar Flares

Recent observations and simulations indicate that solar flares undergo extremely complex 3D evolution, making 3D particle transport models essential for understanding electron acceleration and interpreting flare emissions. In this study, we investigate this problem by solving Parker’s transport equation with 3D MHD simulations of solar flares. By examining energy conversion in the 3D system, we evaluate the roles of different acceleration mechanisms, including reconnection current sheet (CS), termination shock (TS), and supra-arcade downflows (SADs). We find that large-amplitude turbulent fluctuations are generated and sustained in the 3D system. The model results demonstrate that a significant number of electrons are accelerated to hundreds of keV and even a few MeV, forming power-law energy spectra. These energetic particles are widely distributed, with concentrations at the TS and in the flare looptop region, consistent with results derived from recent hard X-ray (HXR) and microwave (MW) observations. By selectively turning particle acceleration on or off in specific regions, we find that the CS and SADs effectively accelerate electrons to several hundred keV, while the TS enables further acceleration to MeV. However, no single mechanism can independently account for the significant number of energetic electrons observed. Instead, the mechanisms work synergistically to produce a large population of accelerated electrons. Our model provides spatially and temporally resolved electron distributions in the whole flare region and at the flare footpoints, enabling synthetic HXR and MW emission modeling for comparison with observations. These results offer important insights into electron acceleration and transport in 3D solar flare regions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Data and scripts associated with “Allometric scaling of hyporheic respiration across basins in the Pacific Northwest USA"

This data package is associated with the publication “Allometric scaling of hyporheic respiration across basins in the Pacific Northwest USA” submitted to JGR-Biogeosciences (Regier et al. 2025).This study used reach-scale modeled estimates of hyporheic aerobic respiration made by the River Corridor Model (Fang et al. 2020) and watershed characteristics across the Willamette and Yakima River basins to explore potential allometric scaling (i.e., power-law relationships between size and function) of cumulative hyporheic respiration across catchment-to-basin scales. Scaling was explored quantitatively via the R2, slope, and y-intercept of relationships between cumulative hyporheic respiration and watershed area, divided into hyporheic exchange flux (HEF) quantiles. We also explored relationships between allometric scaling and other watershed characteristics through linear regression, spatial patterns, and mutual information analyses. Our results also suggest variability of hyporheic respiration allometry for middle exchange flux quantiles, and in relation to land-cover. Our findings provide initial evidence that allometric scaling may be useful for predicting hyporheic biogeochemical dynamics across watersheds from reach to basin scales. This data package is associated with the GitHub repository found at https://github.com/peterregier/rc_wrb_yrb_scaling. The data package is organized into several key directories. The “data” folder contains multiple CSV files, including landscape heterogeneity, scaling analysis, and watershed boundary data. The “figures” folder has all figure files in both PDF and PNG formats. Core analysis scripts and figure generation scripts are in the “scripts” directory, systematically numbered for sequential execution. The root directory includes essential project files; please see the file ending in “flmd.csv” for a list and description of all files contained in this data package and the file ending in “dd.csv” for data dictionaries used to describe tabular column headers.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

RCSB protein data Bank: Next‐generation advanced search for exploration of experimental structures and computed structure models

Abstract The Protein Data Bank (PDB), established in 1971, is the primary global, open‐access archive for experimentally determined 3D macromolecular structures (proteins, RNA, DNA). The research‐focused RCSB.org web‐portal provides access to these data alongside more than one million machine‐learning‐predicted structure models, greatly expanding the available structural landscape. Rapid growth of both experimental and computational structures has increased the need for powerful yet accessible search tools that serve a broad and diverse scientific community. Herein, we describe a redesigned RCSB Protein Data Bank RCSB.org Advanced Search capability that supports intuitive discovery of 3D structures through a unified interface. This interface integrates annotation‐, sequence‐, and 3D structure‐based searches, embeds an interactive 3D viewer, and incorporates curated biological knowledge, such as catalytic site definitions from Mechanism and Catalytic Site Atlas and ligand‐guided structural motifs, for constructing geometry‐driven queries. A new Chemical Search tool allows definition of chemical queries via an integrated drawing tool or standard identifiers, seamlessly combining them with annotation filters. By allowing query definition directly within spatial and chemical contexts, these search interfaces reduce the need for detailed knowledge of residue numbering, chain identifiers, or external cheminformatics software. This capability enables efficient exploration of structures, chemical diversity, and structure–function relationships across all life domains. The redesigned interfaces can be accessed directly at rcsb.org/search/advanced for Advanced Search and rcsb.org/search/chemical for Chemical Search.

Rose, Yana [Research Collaboratory for Structural ↗

NASA NTRS2004

Distributive On-line Processing, Visualization and Analysis System for Gridded Remote Sensing Data

The ability to use data stored in the current Earth Observing System (EOS) archives for studying regional or global phenomena is highly dependent on having a detailed understanding of the data's internal structure and physical implementation. Gaining this understanding and applying it to data reduction is a time- consuming task that must be undertaken before the core investigation can begin. This is an especially difficult challenge when science objectives require users to deal with large multi-sensor data sets that are usually of different formats, structures, and resolutions, for example, when preparing data for input into modeling systems. The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) has taken a major step towards meeting this challenge by developing an infrastructure with a Web interface that allows users to perform interactive analysis online without downloading any data, the GES-DISC Interactive Online Visualization and Analysis Infrastructure or "Giovanni." Giovanni provides interactive, online, analysis tools for data users to facilitate their research. There have been several instances of this interface created to serve TRMM users, Aerosol scientists, Ocean Color and Agriculture applications users. The first generation of these tools support gridded data only. The user selects geophysical parameters, area of interest, time period; and the system generates an output on screen in a matter of seconds. The currently available output options are: Area plot averaged or accumulated over any available data period for any rectangular area; Time plot time series averaged over any rectangular area; Time plots image view of any longitude-time and latitude-time cross sections; ASCII output for all plot types; Image animation for area plot. In the future, we will add correlation plots, GIS-compatible outputs, etc. This allow user to focus on data content (i.e. science parameters) and eliminate the need for expensive learning, development and processing tasks that are redundantly incurred by an archive's user community. The current implementation utilizes the GrADS-DODS Server (GDS), a stable, secure data server that provides subsetting and analysis services across the Internet for any GrADS-readable dataset. The subsetting capability allows users to retrieve a specified temporal and/or spatial subdomain from a large dataset, eliminating the need to download everything simply to access a small relevant portion of a dataset. The analysis capability allows users to retrieve the results of an operation applied to one or more datasets on the server. In our case, we use this approach to read pre-processed binary files and/or to read and extract the needed parts from HDF or HDF-EOS files. These subsets then serve as inputs into GrADS processing and analysis scripts. It can be used in a wide variety of Earth science applications: climate and weather events study and monitoring; modeling. It can be easily configured for new applications.

Leptoukh, G.↗

DOE OSTI2025

PIAFS: A 2D nonlinear hydrodynamics code to model gaseous optics

The survivability of final optics is expected to be a major challenge for all future inertial fusion energy concepts. Due to their higher damage threshold, gaseous optics have been identified as a promising solution to this problem. Gaseous optics can be created through the photoabsorption of spatially modulated UV light, which induces various chemical processes that heat the gas. This heating leads to a pressure perturbation, which in turn launches a density perturbation that can imprint a refractive index modulation such as a grating. In this article, we introduce a parallel C/C++ code to simulate gaseous optics. PIAFS2D is a high-order conservative finite-difference code to solve the compressible Navier–Stokes equations along with the photochemical heating sources on Cartesian grids. The simulations are validated by the linear theory derived in a previous paper [Michel et al., Phys. Rev. Appl. 22, 024014 (2024)]. For larger perturbations, the behavior of the system—particularly the evolution of the generated acoustic wave—demonstrates strong nonlinearity. PIAFS2D allows the study of nonlinear behaviors and can be used for the design of high-efficiency gaseous optics elements in realistic experimental conditions.

Oudin, A. [Lawrence Livermore National Laboratory ↗