Search NASA⌕ Search

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.

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

At least 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Agent-Based Modeling of Microbes in Space

Space is tough on organisms. Microorganisms traveling to space experience stress from environmental features such as ionizing radiation and lack of normal gravity, and much remains unknown about the mechanisms by which those environmental features affect microbial physiology. Microbes experience changes in gravity not directly but rather through changes in their fluid environment, and deep-space particle radiation causes cell damage that is complex but rare. Computational modeling at the single-cell level (agent-based modeling) can allow us to probe the spatially heterogeneous processes that characterize space stresses, to gain insight into the relationships of microbial cells with their environments and with each other. Here we present two software packages for simulating microbial population dynamics in space conditions: CAMDLES and AMMPER. Microbes growing in liquid culture medium in the microgravity of an orbital space station experience a quiescent, poorly-mixed fluid environment. CAMDLES (CFD-DEM Artificial Microgravity Developments for Living Ecosystem Simulation) simultaneously simulates biological, chemical, and mechanical processes to predict microbial ecological dynamics in microgravity, and in the rotating culture vessels used to create an artificial microgravity environment in the lab. Initial results demonstrate that the growth of a cross-feeding microbial consortium, dependent on the exchange of soluble metabolites, is sensitive to the initial spatial distribution of cells, and grows differently in real versus artificial microgravity. Microbial populations exposed to deep-space radiation experience spatially and temporally heterogeneous damage from the traversal of high-energy particles. AMMPER (Agent-Based Model for Microbial Populations Exposed to Radiation) pairs a 3d model of energy deposition along a radiation particle track with a microbial population growth and damage model to predict the effects of localized radiation damage on population-level responses. It includes a user-friendly graphical interface. AMMPER growth curves recapitulate experimental results, and allow comparison between direct effects (DNA damage) and indirect effects (reactive oxygen species generation, metabolic impairment) of radiation.

microbiology↗
DOE OSTI2021

Super-Resolution Photothermal Patterning in Conductive Polymers Enabled by Thermally Activated Solubility

Doping induced solubility control (DISC) patterning is a recently developed technique which uses the change in polymer solubility upon doping, along with an optical dedoping process, to achieve high resolution optical patterning. DISC patterning is remarkable because the process resolution exceeds the linear Abbe’s diffraction limit, however, no mechanism has been proposed to explain such high resolution. Here, we use diffraction to spatially modulate the light intensity and determine the dissolution rate, revealing a superlinear dependence on light intensity, which yields features that are sharper than the profile of the writing photomask. This rate law is independent of wavelength, indicating that patterning resolution is not dominated by an optical dedoping reaction, as was previously proposed. Instead we show here that the optical patterning mechanism is primarily controlled by the thermal profile generated by the laser. To quantify this effect, the thermal profile and dissolution rate are modeled using a finite-element model and compared against patterned line cross sections as a function of wavelength, laser intensity, and dwell time. Our model reveals that although the laser-generated thermal profile is broadened considerably beyond the lasers resolution, the highly temperature dependent dissolution rate results in selective dissolution near the peak of the thermal profile. Therefore, the key factor in achieving super-resolution patterning is a strongly temperature dependent dissolution rate|a common feature of many polymers. In addition to suggesting several routes to improved resolution, our model also demonstrates that doping is not required for optical patterning of conjugated polymers, as was previously believed. Instead, we demonstrate that superlinear resolution optical patterning is attainable in any conjugated polymer simply by tuning the solvent quality during patterning, thus extending the applicability of our method to a wide class of materials. Here, we demonstrate the generality of photothermal patterning by writing sub-400 nm features into undoped Pff-BT4T-2OD.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimizing process-based models to predict current and future soil organic carbon stocks at high-resolution

From hillslope to small catchment scales (< 50 km 2 ), soil carbon management and mitigation policies rely on estimates and projections of soil organic carbon (SOC) stocks. Here we apply a process-based modeling approach that parameterizes the MIcrobial-MIneral Carbon Stabilization (MIMICS) model with SOC measurements and remotely sensed environmental data from the Reynolds Creek Experimental Watershed in SW Idaho, USA. Calibrating model parameters reduced error between simulated and observed SOC stocks by 25%, relative to the initial parameter estimates and better captured local gradients in climate and productivity. The calibrated parameter ensemble was used to produce spatially continuous, high-resolution (10 m 2 ) estimates of stocks and associated uncertainties of litter, microbial biomass, particulate, and protected SOC pools across the complex landscape. Here, subsequent projections of SOC response to idealized environmental disturbances illustrate the spatial complexity of potential SOC vulnerabilities across the watershed. Parametric uncertainty generated physicochemically protected soil C stocks that varied by a mean factor of 4.4 × across individual locations in the watershed and a – 14.9 to + 20.4% range in potential SOC stock response to idealized disturbances, illustrating the need for additional measurements of soil carbon fractions and their turnover time to improve confidence in the MIMICS simulations of SOC dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗
NASA NTRS2007

Evaluation of Integrating the Invasive Species Forecasting System to Support National Park Service Decisions on Fire Management Activities and Invasive Plant Species Control

The USGS and NASA, in conjunction with Colorado State University, George Mason University and other partners, have developed the Invasive Species Forecasting System (ISFS), a flexible tool that capitalizes on NASA's remote sensing resource to produce dynamic habitat maps of invasive terrestrial plant species across the United States. In 2006 ISFS was adopted to generate predictive invasive habitat maps to benefit noxious plant and fire management teams in three major National Park systems: The Greater Yellowstone Area (Yellowstone / Grand Tetons National Parks), Sequoia and Kings Canyon National Park, and interior Alaskan (between Denali, Gates of The Arctic and Yukon-Charley). One of the objectives of this study is to explore how the ISFS enhances decision support apparatus in use by National Park management teams. The first step with each park system was to work closely with park managers to select top-priority invasive species. Specific species were chosen for each study area based on management priorities, availability of observational data, and their potential for invasion after fire disturbances. Once focal species were selected, sources of presence/absence data were collected from previous surveys for each species in and around the Parks. Using logistic regression to couple presence/absence points with environmental data layers, the first round of ISFS habitat suitability maps were generated for each National Park system and presented during park visits over the summer of 2006. This first engagement provided a demonstration of what the park service can expect from ISFS and initiated the ongoing dialog on how the parks can best utilized the system to enhance their decisions related to invasive species control. During the park visits it was discovered that separate "expert opinion" maps would provide a valuable baseline to compare against the ISFS model output. Opinion maps are a means of spatially representing qualitative knowledge into a quantitative two-dimensional map. Furthermore, our approach combines the qualitative expert opinion habitat maps -- with the quantitative ISFS habitat maps in a difference map that shows where the two maps agree and disagree. The objective of the difference map is to help focus future field sampling and improve model results. This paper presents a demonstration of the habitat, expert opinion, and difference map for Yellowstone National Park.

Ma, Peter↗
NASA NTRS2010

Super Resolution Image Enhancement for a Flash Lidar: Back Projection Method

In this paper a new image processing technique for flash LIDAR data is presented as a potential tool to enable safe and precise spacecraft landings in future robotic or crewed lunar and planetary missions. Flash LIDARs can generate, in real-time, range data that can be interpreted as a 3-dimensional (3-D) image and transformed into a corresponding digital elevation map (DEM). The NASA Autonomous Landing and Hazard Avoidance (ALHAT) project is capitalizing on this new technology by developing, testing and analyzing flash LIDARs to detect hazardous terrain features such as craters, rocks, and slopes during the descent phase of spacecraft landings. Using a flash LIDAR for this application looks very promising, however through theoretical and simulation analysis the ALHAT team has determined that a single frame, or mosaic, of flash LIDAR data may not be sufficient to build a landing site DEM with acceptable spatial resolution, precision, size, or for a mosaic, in time, to meet current system requirements. One way to overcome this potential limitation is by enhancing the flash LIDAR output images. We propose a new super-resolution algorithm applicable to flash LIDAR range data that will create a DEM with sufficient accuracy, precision and size to meet current ALHAT requirements. The performance of our super-resolution algorithm is analyzed by processing data generated during a series of simulation runs by a high fidelity model of a flash LIDAR imaging a high resolution synthetic lunar elevation map. The flash LIDAR model is attached to a simulated spacecraft by a gimbal that points the LIDAR to a target landing site. For each simulation run, a sequence of flash LIDAR frames is recorded and processed as the spacecraft descends toward the landing site. Each run has a different trajectory profile with varying LIDAR look angles of the terrain. We process the output LIDAR frames using our SR algorithm and the results show that the achieved level of accuracy and precision of the SR generated landing site DEM is more than adequate for detecting hazardous terrain features and identifying safe areas.

Bulyshev, Alexander↗
NASA NTRS2013

Direct Simulation of Extinction in a Slab of Spherical Particles

The exact multiple sphere superposition method is used to calculate the coherent and incoherent contributions to the ensemble-averaged electric field amplitude and Poynting vector in systems of randomly positioned nonabsorbing spherical particles. The target systems consist of cylindrical volumes, with radius several times larger than length, containing spheres with positional configurations generated by a Monte Carlo sampling method. Spatially dependent values for coherent electric field amplitude, coherent energy flux, and diffuse energy flux, are calculated by averaging of exact local field and flux values over multiple configurations and over spatially independent directions for fixed target geometry, sphere properties, and sphere volume fraction. Our results reveal exponential attenuation of the coherent field and the coherent energy flux inside the particulate layer and thereby further corroborate the general methodology of the microphysical radiative transfer theory. An effective medium model based on plane wave transmission and reflection by a plane layer is used to model the dependence of the coherent electric field on particle packing density. The effective attenuation coefficient of the random medium, computed from the direct simulations, is found to agree closely with effective medium theories and with measurements. In addition, the simulation results reveal the presence of a counter-propagating component to the coherent field, which arises due to the internal reflection of the main coherent field component by the target boundary. The characteristics of the diffuse flux are compared to, and found to be consistent with, a model based on the diffusion approximation of the radiative transfer theory.

spheres↗
NASA NTRSDate not supplied

Deep Learning System for Efficient Processing of Geostationary Satellite Imagery

Improved capabilities of Earth monitoring satellites are enabling a wide range of studies on the environmental effects of climate change, often leveraging the recent advancements in machine learning. At the same time, the new capabilities, including higher spatial resolution and temporal frequency, are expanding the amount of data generated at exponential rates. Further, a large majority of archived datasets generated by scientific processing is never used. This motivates the development of an efficient machine learning system for end-to-end processing of multi-level satellite datasets, from level 1 top of atmosphere observations to user friendly environmental variables of interest. Using current generation geostationary satellites GOES-16/17 (NOAA/NASA), and Himawari-8/9 (JAXA), we present an interchangeable set of machine models to perform spectral adjustment among sensors, physical model emulation, LEO-GEO emulation, and optical flow in a high performance computing environment. We use these tools on the NASA Earth eXchange (NEX) to generate consistent virtual observations across sensors, perform atmospheric correction and cloud detection, and estimate surface reflectance, surface temperature and atmospheric winds. This approach aims to improve the robustness of remotely sensed data processing by learning from diverse sets of observations while enabling near real-time and on-demand capabilities.

Thomas Vandal↗
DOE OSTI2017

3D Geophysical Inversion Modeling of Ground Magnetic Data at Baker Hot Springs, Washington State, USA

Ground-based magnetic survey data were geophysically modelled in 3D as part of a play fairway exploration program at Baker hot springs. Here, a multi-step inversion method was used to address the issue of remanent magnetization in rocks of the study area. The geophysical modelling effort generated an effective susceptibility model that displays the 3D distribution and orientation of magnetization in the subsurface. Importantly, the magnetic modelling effort found that a strong, dipolar magnetic anomaly observed in the study area can be explained by a remanently magnetized body that is ~1¼ km long x ~200 m wide and extends from ~25-400 m below the surface. Geologically, this body may represent 1) an intrusion of reversely magnetized mafic rocks or 2) a magnetite-bearing fault zone. This rock body coincides spatially with Baker hot springs; the margins of this body or the fault structures within it are likely responsible for providing the permeability that gives rise to the hot springs.

Baker hot springs↗
DOE OSTI2021

Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations". We used field observations of biomass and leaf traits (Salmon et al., 2019b) across a gradient of plant communities on the Seward Peninsula in northwest Alaska (Breen et al., 2020; Iversen et al., 2017) to replace the original ELM configuration with nine, arctic-specific PFTs. Original arctic PFTs in the model were (1) broadleaf deciduous boreal shrub and (2) C3 arctic grass. New PFTs were (1) lichen, (2) bryophyte, (3) evergreen dwarf shrub, (4) deciduous dwarf shrub, (5) deciduous low shrub, (6) deciduous low to tall shrub, (7) alder shrub, (8) forb, and (9) graminoid. New PFTs were parameterized and evaluated using site-level trait and biomass measurements (Salmon et al., 2019b), root C:N ratios and rooting depth patterns (Salmon et al., 2019c), and root biomass (Salmon et al., 2020). This archive contains forcing data and model output from ELM simulations conducted at the Kougarok Hillslope site (Kougarok Road Mile Marker 64). Simulations were conducted for the Kougarok Hillslope site using meteorological driving data from the Scenarios Network for Alaska and Arctic Planning (SNAP) downscaled climate projection dataset using NCAR-CCSM forcing (Bieniek et al., 2020 dataset; Walsh et al., 2018). The archive includes three versions of model simulations representing a progression from the original model configuration to an Arctic-specific configuration (see Output section below). Simulation 1 used grid cell data from a global E3SM configuration, including the fractional area of E3SM PFTs assigned to the grid cell containing the Kougarok Hillslope site in global simulations. This simulation did not distinguish between different plant communities on the Kougarok Hillslope, but instead used a single point simulation to represent the entire area. Simulation 2 used default E3SM PFT definitions combined with adjusted depth to bedrock for plant communities with shallow rocky layers to represent the role of abiotic soil factors in driving site differences (based on depth measurements of Iversen et al., 2019a), and adjusted the relative areas of the broadleaf deciduous boreal shrub and C3 arctic grass PFTs to reflect the observed spatial coverage of shrub and non-shrub PFTs across Kougarok Hillslope plant communities. Simulation 3 used the new arctic PFT definitions and parameterizations based on the vegetation types present at the site, along with the adjusted soil depths used in Simulation 2. Included are *.pdf, *.nc (NetCDF), *tar.gz, *.py, and *.bash files. We would like to thank Mary's Igloo Native Corporation for allowing us to perform this research on their land and giving us the opportunity to share this research with their community. Thanks to Shawn Serbin and Alistair Rogers for helpful discussions. Holly Vander Stel, Joanne Childs, and Stan Wullschleger helped with collection of field data. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Waste Characterization Methods Based on Gamma-ray Imaging Spectrometers - 19454

Characterization of LLW, ILW, U, NORM and TENORM is critical to the nuclear industry's bottom line. Accurate measurements enable the disposal of waste at the lowest possible cost. In the past the most accurate technology available involved using HPGe detectors and extensive system modeling to estimate the total source term given what the user knows about the geometry of the source and any shielding between the source and the detectors. Now there are a variety of imaging detectors that are available commercially, including CZT and HPGe detectors that achieve energy resolution ≤1% FWHM at 662 keV. These devices couple the traditional high resolution spectroscopic analysis with next generation gamma-ray imaging information to more accurately characterize the radiation source term. Uncertainties in the source location will result in significant error in the characterization of an object's radioactivity. It is helpful if the object can be spun, such as a barrel, so that any asymmetry in the source distribution is averaged out from the detector's point of view. Even in the case of a spinning barrel it is important to determine if the radiation is concentrated near the center of the barrel or distributed near the edge. The closer one measures to a barrel, which usually dictates the time required to acquire data, the more these spatial variations in the source distribution will affect the calculated activity of each isotope. This paper will focus on showing the use of CZT imaging spectrometers to characterize complicated radiation source terms. The radiation images are used to determine where the source term should be defined in a radiation transport model that calculates the expected number of photopeak counts for each isotope detected. This is used to calculate the activity of the source. The transport model is used to determine the attenuation of the radiation source term through distance and shielding material, but the detection efficiency is based solely on radiation measurement data. These measurements are taken with a robotic source positioning device that allows for detailed directional efficiency calibration. It is a fully automated system where a detector is loaded into the device and all of the efficiency curves are captured as a function of direction by creating new measurements for each location the robot positions the source. These measurements are tagged with the true source location and the efficiency curve is calculated for each measurement. This is used to create a response function matrix as a function of energy and incident direction which is used to convert the direction dependent flux on the detector head determined from the simulation model to the expected number of counts per unit activity of each isotope of interest. The results of these source term reconstructions in both the spatial and spectral domain will be shown. The impact of improper specification of the source term will be shown for a few standard waste characterization scenarios to demonstrate the importance of proper source location determination. Finally, the underlying calibration method and some sample efficiency curves as a function of direction will be shown to demonstrate the importance of calibrating out directional response functions in imaging devices (which often have a strong forward sensitivity bias). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2018

Three-dimensional kinetic modeling of streamer propagation in a nitrogen/helium gas mixture

A fully resolved kinetic model (Particle-in-Cell and Direct Simulation Monte Carlo for particle/photon collisions) of a near atmospheric pressure ionization wave is presented here. Fully resolving the required numerical spatial (sub-um) and temporal scales (tens of fs) for atmospheric pressure discharges in three dimensions is still a challenging task on modern super computers. To keep the overall problem tractable, the total number of elements are reduced by only simulating a 10 degree wedge rather than a full 360 degree geometry. The ionization wave is generated in a needle-plane configuration with a gap size of 250 um and a background of nitrogen and helium gas. A voltage of 1500 V is applied to the anode and an initial electron and ion density of 10 9 cm 3 is seeded in a region near the anode electrode tip and extending towards the cathode. As these initial electrons are swept away, photoionization and photoemission create new electrons and allow the ionization front to propagate towards the cathode. Results from the 90% N 2 , 10% He discharge indicate that photoionization has minimal impact on plasma formation processes and cathode photoemission is the dominant mechanism for new electrons. In the 90% He, 10% N 2 discharge case, however, photoionization likely has an impact as the observed locations of photoionization occur far enough away from the ionization front to allow for sufficient avalanche processes that contribute to the propagation of the ionization wave. Additionally, the electron energy distribution functions in the 90% He, 10% N 2 case indicate that there is less energy loss to the low lying molecular N 2 electronic states as well as the vibrational and rotational modes. This leads to higher electron energies and faster plasma development times of ~0.4 ns for the 90% He, 10% N 2 case, and ~1.5 ns for the 90% N 2 , 10% He case. In addition to analysis of the ionization wave results, the overall challenges associated with simulations near atmospheric pressure discharges in three dimensions are discussed, including the limitations of the 10 degree wedge that produces, at least qualitatively, minimal 3D effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗
NASA NTRS1998

Prototype Local Data Integration System and Central Florida Data Deficiency

This report describes the Applied Meteorology Unit's (AMU) task on the Local Data Integration System (LDIS) and central Florida data deficiency. The objectives of the task are to identify all existing meteorological data sources within 250 km of the Kennedy Space Center (KSC) and the Eastern Range at Cape Canaveral Air Station (CCAS), identify and configure an appropriate LDIS to integrate these data, and implement a working prototype to be used for limited case studies and data non-incorporation (DNI) experiments. The ultimate goal for running LDIS is to generate products that may enhance weather nowcasts and short-range (less than 6 h) forecasts issued in support of the 45th Weather Squadron (45 WS), Spaceflight Meteorology Group (SMG), and the Melbourne National Weather Service (NWS MLB) operational requirements. The LDIS has the potential to provide added value for nowcasts and short term forecasts for two reasons. First, it incorporates all data operationally available in east central Florida. Second, it is run at finer spatial and temporal resolutions than current national-scale operational models. In combination with a suitable visualization tool, LDIS may provide users with a more complete and comprehensive understanding of evolving fine-scale weather features than could be developed by individually examining the disparate data sets over the same area and time. The utility of LDIS depends largely on the reliability and availability of observational data. Therefore, it is important to document all existing meteorological data sources around central Florida that can be incorporated by it. Several factors contribute to the data density and coverage over east central Florida including the level in the atmosphere, distance from KSC/CCAS, time, and prevailing weather. The central Florida mesonet consists of existing surface meteorological and hydrological data available from the Tampa NWS and data servers at Miami and Jacksonville. However the utility of these data for operational use is limited, mainly because there are relatively few additional meteorological observations within 50 km of KSC/CCAS to supplement existing METAR and KSC/CCAS tower reports.

Manobianco, John↗
DOE OSTI2024

A three-point velocity estimation method for two-dimensional coarse-grained imaging data

Time delay and velocity estimation methods have been widely studied subjects in the context of signal processing, with applications in many different fields of physics. The velocity of waves or coherent fluctuation structures is commonly estimated as the distance between two measurement points divided by the time lag that maximizes the cross correlation function between the measured signals, but this is demonstrated to result in erroneous estimates for two spatial dimensions. We present an improved method to accurately estimate both components of the velocity vector, relying on three non-aligned measurement points. We introduce a stochastic process describing the fluctuations as a superposition of uncorrelated pulses moving in two dimensions. Using this model, we show that the three-point velocity estimation method, using time delays calculated through cross correlations, yields the exact velocity components when all pulses have the same velocity. The two- and three-point methods are tested on synthetic data generated from realizations of such processes for which the underlying velocity components are known. The results reveal the superiority of the three-point technique. Finally, we demonstrate the applicability of the velocity estimation on gas puff imaging data of strongly intermittent plasma fluctuations due to the radial motion of coherent, blob-like structures at the boundary of the Alcator C-Mod tokamak.

Materials Science↗
DOE OSTI2025

Field-level reconstruction from foreground-contaminated 21-cm maps

Current and upcoming 21-cm experiments will soon be able to map 21-cm spatial fluctuations in three dimensions for a wide range of redshifts. However, bright foreground contamination and the nature of radio interferometry create significant challenges, making it difficult to access rich cosmological information from the Fourier modes that lie within the “foreground wedge”. Here, in this work, we introduce two approaches aiming to reconstruct the full 21-cm density field, including the missing modes in the wedge: (a) a field-level inference under an effective field theory (EFT) framework; (b) a diffusion-based deep generative model trained on simulations. Under the EFT framework, we implement a fully differentiable forward model that maps the initial conditions of matter fluctuations to the observed, foreground-filtered 21-cm maps. This enables a gradient-based sampler to simultaneously sample the initial conditions and bias parameters, allowing a physically motivated mode reconstruction. Alternatively, we apply a variational diffusion model to perform 21-cm density reconstruction at the map level. Our model is trained on semi-numerical simulations over a wide range of astrophysical parameters. Our results from both approaches should provide improved cosmological constraints from the field level and also enable cross-correlation between experiments that have little or no overlapping modes.

cosmological perturbation theory↗