SEARCH · Search NASA

Results for “reference cell 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 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Simulation and Experimental Measurements of Inductively Coupled CF4 and CF4/Ar Plasmas

The recently developed code SEMS (semiconductor equipment modeling software)is applied to the simulation of CF4 and CF4/Ar inductively coupled plasmas (ICP). This work builds upon the earlier nitrogen, transformer coupled plasma (TCP) SEMS research by demonstrating its accuracy for more complex reactive mixtures, moving closer to the realization of a virtual plasma reactor. Attention is given to the etching of and/or formation of carbonaceous films on the quartz dielectric window and diagnostic aperatures. The simulations are validated through comparisons with experimental measurements using FTIR (Fourier Transform Infrared) and UV absorption spectroscopy for CFx and SiFx neutral radicals, QMS (quadrupole mass spectrometry) for the ions, and Langmuir probe measurements of electron number density and temperature in an ICP GEC reference cell.

Hash, D. B.↗

DOE OSTI2017

Washington Clean Energy Fund: Energy Storage System Performance Test Plans and Data Requirements

This report documents the test plans, including detailed duty cycles, used in evaluating the technical performance of five energy storage systems (ESSs) sponsored by the Washington State Clean Energy Fund (CEF). These ESSs were installed at three participating utilities’ (Avista, Snohomish Public Utility District or SnoPUD, and Puget Sound Energy or PSE) designated sites. Under each test plan, a number of baseline and use case tests are defined. This report documents the test plan for each utility, including duty cycles, data and performance reporting requirements. Baseline tests defined in this report mainly focus on determining an initial or reference performance for each ESS, which may be compared with the performance at any later time to assess the state of health of the ESS. Parameter values of a given ESS obtained from baseline tests are applied for different modeling purposes and use case tests. The stored energy capacity test is the first test conducted in the baseline test program, which generates data to calculate round trip efficiency (RTE). The response time and ramp rate tests provide the time required for an ESS to change from zero to full charging/discharging rate and hence the ramp rate, which is important in understanding ESS performance for applications calling for fast response. Internal resistance testing of ESS cells is also conducted in conjunction with the response time and ramp rate tests. Finally, the baseline tests include frequency regulation (FR) and peak shaving (PS) tests based on real-world duty cycles. Upon completion of the baseline tests, use case tests are conducted at each of the storage sites with duty cycles developed to meet the site-specific criteria and conditions. The use cases are grouped according to their intended target benefits within the electric infrastructure topology (e.g., transmission versus distribution). Although an ESS may be located on the low-voltage side of a substation that provides power to a distribution feeder, a use case that addresses bulk power services could still be provided and would be grouped under the transmission or bulk-power benefits. Use cases for ESS applications deep into the distribution circuit would be categorized under the distribution system cases.

25 ENERGY STORAGE↗

DOE OSTI2022

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

BSEC VPRM 10m Hourly Biogenic Fluxes in Baltimore (2021)

Model outputs from the Vegetation Photosynthesis and Respiration Model (VPRM: version from Horne et al. in prep). Model remote sensing inputs come from Sential 2-derived EVI and LSWI. Model meteorological inputs for two-meter air temperature and shortwave incoming come from the BSEC WRF 2021 Control Run (Foust, W. 2023). Plant functional Types (PFTs) are spatially classified using the Chesapeake Bay Program 2018 land use land cover product. The final biogenic flux (µmol CO2 m^-2 s^-1) outputs of NEE, RESP, and GEE are a weighted average based on the portion of PFTs within the cell. Individual PFT outputs are saved inside PFT directories (e.g., Crops, Grass, etc.) inside the specific month directory. Model outputs are denoted as a negative flux into the land system (i.e., photosynthesis) and a positive flux as a net release into the overlying atmosphere. Respiration (RESP) fluxes are positive and combine heterotrophic (only soil) and autotrophic sources. Gross ecosystem exchange (GEE) is a negative flux driven by only photosynthetic activity from vegetation, and the Net ecosystem exchange (NEE) is the sum of the two (i.e., NEE=RESP+GEE). Data Characteristics Spatial Resolution: 10m Temporal Resolution: Hourly File Format: VPRM_ _BSEC. .tif (Hour is in UTC) For more information on the model results, please email Jason Horne (jph6488@psu.edu). References: Foust, W. (2023). BSEC WRF 2021 Control Run Output (v0.1.0) [Data set]. MSD-LIVE Data Repository. https://data.msdlive.org/records/m0e6m-vvq17

Baltimore↗

NASA NTRS2000

Redistribution of annexin in gravistimulated pea plumules

We used immunocytochemistry to investigate the effects of gravistimulation on annexin localization in etiolated pea plumule shoots. In longitudinal sections, an asymmetric annexin immunostaining pattern was observed in a defined group of cells located just basipetal to apical meristems at the main shoot apex and at all of the axillary buds, an area classically referred to as the leaf gap. The pattern was observed using both protein-A-purified anti-annexin and affinity-purified anti-annexin antibodies for the immunostaining. A subset of the cells with the annexin staining also showed an unusually high level of periodic acid Schiff (PAS) staining in their cell walls. Prior to gravistimulation, the highest concentration of annexin was oriented toward the direction of gravity along the apical end of these immunostained cells. In contrast, both at 15 and 30 min after gravistimulation, the annexin immunostain became more evenly distributed all around the cell and more distinctly cell peripheral. The asymmetry along the lower wall of these cells was no longer evident. In accord with current models of annexin action, we interpret the results to indicate that annexin-mediated secretion in the leaf gap area is preferentially toward the apical meristem prior to gravistimulation, and that gravistimulation results in a redirection of this secretion. These data are to our knowledge the first to show a correlation between the vector of gravity and the distribution of annexins in the cells of flowering plants. c 2000 Editions scientifiques et medicales Elsevier SAS.

Non-NASA Center↗

DOE OSTI2022

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Parts & Pools: A Framework for Modular Design of Synthetic Gene Circuits

Published in 2008, Parts & Pools represents one of the first attempts to conceptualize the modular design of bacterial synthetic gene circuits with Standard Biological Parts (DNA segments) and Pools of molecules referred to as common signal carriers (e.g., RNA polymerases and ribosomes). The original framework for modeling bacterial components and designing prokaryotic circuits evolved over the last years and brought, first, to the development of an algorithm for the automatic design of Boolean gene circuits. This is a remarkable achievement since gene digital circuits have a broad range of applications that goes from biosensors for health and environment care to computational devices. More recently, Parts & Pools was enabled to give a proper formal description of eukaryotic biological circuit components. This was possible by employing a rule-based modeling approach, a technique that permits a faithful calculation of all the species and reactions involved in complex systems such as eukaryotic cells and compartments. In this way, Parts & Pools is currently suitable for the visual and modular design of synthetic gene circuits in yeast and mammalian cells too.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

CICE on a C-grid: new momentum, stress, and transport schemes for CICEv6.5

Abstract. This article presents the C-grid implementation of the CICE sea ice model, including the C-grid discretization of the momentum equation, the boundary conditions (BCs), and the modifications to the code required to use the incremental remapping transport scheme. To validate the new C-grid implementation, many numerical experiments were conducted and compared to the B-grid solutions. In idealized experiments, the standard advection method (incremental remapping with C-grid velocities interpolated to the cell corners) leads to a checkerboard pattern. A modal analysis demonstrates that this computational noise originates from the spatial averaging of C-grid velocities at corners. The checkerboard pattern can be eliminated by adjusting the departure regions to match the divergence obtained from the solution of the momentum equation. We refer to this novel approach as the edge flux adjustment (EFA) method. The C-grid discretization with edge flux adjustment allows for transport in channels that are one grid cell wide – a capability that is not possible with the B-grid discretization nor with the C-grid and standard remapping advection. Simulation results match the predicted values of a novel analytical solution for one-grid-cell-wide channels.

Lemieux, Jean-François (ORCID:0000000320845759)↗

DOE OSTI2025

Cooper-pair density modulation state in an iron-based superconductor

Superconducting (SC) states that break space-group symmetries of the underlying crystal can exhibit nontrivial spatial modulation of the order parameter. Previously, such states were intimately associated with the breaking of translational symmetry, resulting in the density-wave orders, with wavelengths spanning several unit cells. However, a related basic concept has long been overlooked: when only intra-unit-cell symmetries of the space group are broken, the SC states can show a distinct type of nontrivial modulation preserving long-range lattice translation. Here, in this study, we refer to this new concept as the pair density modulation (PDM) and report the first observation of a PDM state in exfoliated thin flakes of the iron-based superconductor FeTe 0.55 Se 0.45 . Using scanning tunnelling microscopy (STM), we discover robust SC gap modulation with the wavelength corresponding to the lattice periodicity and the amplitude exceeding 30% of the gap average. Notably, we find that the observed modulation originates from the large difference in SC gaps on the two nominally equivalent iron sublattices. The experimental findings, backed up by model calculations, suggest that, in contrast to the density-wave orders, the PDM state is driven by the interplay of sublattice symmetry breaking and a peculiar nematic distortion specific to the thin flakes. Our results establish new frontiers for exploring the intertwined orders in strong-correlated electronic systems and open a new chapter for iron-based superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

Quantitative Methods for Comparing Different Polyline Stream Network Models

Two techniques for exploring relative horizontal accuracy of complex linear spatial features are described and sample source code (pseudo code) is presented for this purpose. The first technique, relative sinuosity, is presented as a measure of the complexity or detail of a polyline network in comparison to a reference network. We term the second technique longitudinal root mean squared error (LRMSE) and present it as a means for quantitatively assessing the horizontal variance between two polyline data sets representing digitized (reference) and derived stream and river networks. Both relative sinuosity and LRMSE are shown to be suitable measures of horizontal stream network accuracy for assessing quality and variation in linear features. Both techniques have been used in two recent investigations involving extracting of hydrographic features from LiDAR elevation data. One confirmed that, with the greatly increased resolution of LiDAR data, smaller cell sizes yielded better stream network delineations, based on sinuosity and LRMSE, when using LiDAR-derived DEMs. The other demonstrated a new method of delineating stream channels directly from LiDAR point clouds, without the intermediate step of deriving a DEM, showing that the direct delineation from LiDAR point clouds yielded an excellent and much better match, as indicated by the LRMSE.

58 GEOSCIENCES↗

DOE OSTI2024

Bifacial Vertical Testbed and Ground Irradiance Data in Golden, Colorado

This data was collected for Tonita et al., “Vertical bifacial photovoltaic system model validation: study with field data, various orientations, and latitudes,” for validation of optical models for vertically-oriented photovoltaics under high albedo. Ground irradiance data for vertical PV arrays modeling in agrivoltaics is also provided. The dataset is provided for further use or study as open source. For any questions on the dataset, email silvana.ovaitt@nlr.gov.

14 SOLAR ENERGY↗

DOE OSTI2016

Confirming a pyrolitic lower mantle using self-consistent pressure scales and new constraints on CaSiO 3 perovskite

In this study, we have examined the lower mantle composition and mineralogy by modeling the density ( ρ ), bulk sound velocity ( V Φ ), and dln ρ /dln V Φ profiles of candidate lower mantle minerals using literature and new experimental equation of state (EoS) results. Additionally, for CaSiO 3 perovskite, complimentary synchrotron X-ray diffraction measurements in a laser-heated diamond anvil cell were conducted up to 156 GPa between 1200 K and 2600 K to provide more reliable constraints on the thermal EoS parameters. These new experimental results as well as literature P-V-T data sets are systematically analyzed using an internally self-consistent pressure scale. We have modeled ρ , V Φ , and dln ρ /dln V Φ profiles of the lower mantle with representative pyrolitic and chondritic compositional models in which the effect of Fe spin transition in ferropericlase is also taken into account. Our modeling results show that a pyrolitic lower mantle with an aggregate mineralogy of 75 vol % bridgmanite, 17 vol % ferropericlase, and 8 vol % CaSiO 3 perovskite produces ρ and V Φ profiles in better agreement with preliminary reference Earth model than a lower mantle with a chondritic composition. The modeled ρ , V Φ , and dln ρ /dln V Φ are mainly affected by the relative ratio of bridgmanite and ferropericlase but are not sensitive to the variation of the CaSiO 3 perovskite content. In addition, the spin crossover of Fe in ferropericlase can greatly raise the value of dln ρ /dln V Φ in the middle lower mantle, which is useful in detecting the presence of ferropericlase in the region. Based on these new mineral physical constraints and radial seismic structure, our study suggests the lower mantle is pyrolitic, which is chemically indistinguishable from the upper mantle.

58 GEOSCIENCES↗

DOE OSTI2025

Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection

Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection Description This dataset contains input and output data for the manuscript Mongird, K. et al. (under review) titled "Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection". Input data corresponds to gridded spatial siting attributes that are necessary to conduct a random forest machine learning analysis of siting feature importance. Output data includes SHAP feature analysis outputs, and classification report values. For data on power plant siting results referred to in the manuscript, please refer to the CERF: IM3 Projected Western US Power Plant Locations data download page. The downloadable data includes values for eight different future scenarios for the Western US. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 Technical Information The dataset includes two sets of data files: (1) CERF gridded siting parameters and (2) Feature analysis outputs and classification reports. All downloadable data is in csv file format. Files with x/y coordinate information use the Albers Equal Area Conic projection (ESRI:102003). 1. CERF Gridded Siting Parameters This directory provides a balanced sample of gridded CERF siting parameters data for eight different scenarios for the Western US through 2055, seven different technologies, and eight timesteps. This data serves as input to the feature analysis. It contains the following parameters. region_name - name of region (i.e., state) sited - binary value representing whether the grid cell received a siting of that technology type (1=True) rcp - binary value representing scenario resource concentration pathway (0 = RCP4.5, 1 = RCP8.5) ssp - binary value representing scenario shared socioeconomic pathway (0 = SSP3, 1 = SSP5) climate - binary value representing cooler (0) or hotter (1) GCM forcing tech_name - generation technology name sited_year - year that values correspond to transmission_cost - cost of transmission interconnection pipeline_cost - cost of natural gas pipeline interconnection interconnection_cost - total interconnection cost (sum of transmission cost and gas pipeline cost) lmp - associated locational marginal value ($/MWh) associated with the grid cell, timestep, scenario, and technology xcoord - x-coordinate of location ycoord - y-coordinate of location 2a. Feature Analysis Output The dataset includes the feature analysis shap output for locational marginal price and interconnection cost. It contains the following parameters. technology - generator technology name scenario - name of scenario feature - name of feature, either locational_marginal_price or interconnection_cost value - the mean of absolute value of SHAP values for given feature 2b. Feature Analysis Classification Report This download includes the classification report associated with each random forest model. The dataset contains the following parameters. technology - generation technology name scenario - name of scenario test - one of precision (the proportion of predicted positives that are actually correct), recall (the proportion of actual positives that were correctly identified), f1-score (the harmonic mean of precision and recall) 0.0 - value of test for classification of 0 (grid cell not chosen for siting) 1.0 - value of test for classification of 1 (grid cell chosen for siting) accuracy - accuracy of model (i.e., fraction of all predictions that were right) macro avg - Simple average of test values for all classes weighted avg - Weighted average of test values for all classes, weighted based on Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall [Pacific Northwest National Labor↗

DOE OSTI2014

Adsorption Properties of p -Methyl Red Monomeric-to-Pentameric Dye Aggregates on Anatase (101) Titania Surfaces: First-Principles Calculations of Dye/TiO 2 Photoanode Interfaces for Dye-Sensitized Solar Cells

The optical and electronic properties of dye aggregates of p-methyl red on a TiO2 anatase (101) surface were modeled as a function of aggregation order (monomer to pentameric dye) via first principles calculations. A progressive red-shifting and intensity increase toward the visible region in UV/vis absorption spectra is observed from monomeric-to-tetrameric dyes, with each molecule in a given aggregate binding to one of the four possible TiO2 (101) adsorption sites. The pentamer exhibits a blue-shifted peak wave- length in the UV/vis absorption spectra and less absorption intensity in the visible region in comparison; a corresponding manifestation of H-aggregation occurs since one of these five molecules cannot occupy an adsorption site. This finding is consistent with experiment. Calculated density of states (DOS) and partial DOS spectra reveal similar dye…TiO2 nanocomposite conduction band characteristics but different valence band features. Associated molecular orbital distributions reveal dye-to-TiO2 interfacial charge transfer in all five differing aggregate orders; meanwhile, the level of intramolecular charge transfer in the dye becomes progressively localized around its azo- and electron-donating groups, up to the tetrameric dye/TiO2 species. Dye adsorption energies and dye coverage levels are calculated and compared with experiment. Overall, the findings of this case study serve to aid the molecular design of azo dyes towards better performing DSSC devices wherein they are incorporated. In addition, they provide a helpful example reference for understanding the effects of dye aggregation on the adsorbate…TiO2 interfacial optical and electronic properties.

36 MATERIALS SCIENCE↗

DOE OSTI2013

Theory Aided Design of Active and Durable Nanoscale Cathode Catalysts (Final Report)

One of the greatest technological hurdles to deployment of fuel cells relates to the sluggish activity, low durability and the high cost of the catalysts that are currently employed. For automotive PEM fuel cells to become commercially viable, the Pt-specific power density would need to be reduced to less than 0.2gPt/kW (at cell voltages >0.65 V). This would require the Pt loadings to be less than 0.15 mg Pt /cm 2 MEA within the membrane electrode assembly. This could be achieved by enhancing the catalytic activity at the cathode, thus lowering its overpotential. Various different Pt-alloys have shown 2-4 times enhanced activities over Pt alone but still suffer some of the same durability issues as that of the pure Pt. There is a general loss of active Pt due to dissolution and sintering. While there have been a number of elegant fundamental experimental and theoretical studies on ideal single crystal Pt and Pt alloy surfaces which have helped to elucidate the factors that control the activity, there have been very few fundamental studies focused on the stability, reactivity and durability of well-defined Pt nanoparticles. We carried out ab initio density functional theory together with a novel double reference method that we developed to simulate constant potential electrochemical systems in order to model the electrocatalytic reduction of oxygen over model Pt alloy surfaces and nanoparticles. These simulations were used to probe the factors that control the electrocatalytic activity, guide the potential selection of new materials and test their stability under reaction conditions. Ab initio calculations were used to determine the reaction energies and activation barriers for a comprehensive array of different elementary adsorption, desorption, surface reaction and diffusion steps over Pt and Pt alloys involved in the electroreduction of oxygen as a function of surface coverage, the electrochemical potential and temperature. The calculations were used to simulate the potential dependent adsorption and surface reaction energies along with activation barriers in order to determine the kinetics for different surface structures and structural features (step edge and corner sites) to provide the necessary input for kinetic Monte Carlo simulations to follow the rates of reaction. More coarse-grained simulated annealing methods were used to help establish the lowest energy structures and morphologies for different Pt and Pt-alloy nanoparticles that form. The results from the DFT calculations were used to establish an ab initio-derived kinetic database that was used in both 2D and 3D kinetic Monte Carlo simulations that were used to follow the electrocatalytic performance over different particle sizes, shapes and compositions. The ab initio calculations together with the kinetic Monte Carlo simulations were used to complete the following objectives: 1) Determine the controlling elementary reaction pathways and intrinsic kinetics involved in ORR and their potential dependent behavior; 2) Establish the influence of the extrinsic reaction environment including surface structure, alloy composition and spatial arrangement and the humidity on ORR kinetics; 3) Elucidate the effects of particle size and morphology as well as the atomic structure and composition of nanoparticles of Pt, Pt 3 Co, Pt 3 Ni, Pt 3 Fe and other Pt-alloys on the electrocatalytic activity and reaction selectivity; 5) Construct surface Pourbaix phase diagrams for different Pt and Pt-alloys in order to map out the most stable surface phases for different material compositions and reaction conditions; 6) Elucidate the mechanisms that control dissolution and re-deposition of Pt.

08 HYDROGEN↗

DOE OSTI2025

Modeling neutral defects in III-V ternary alloys with a special quasirandom structure: Analysis of As- and III-site point defects in InGaAs

While first-principles density functional theory modeling has become a vital tool to investigate defect properties in semiconductors, the lack of crystalline periodicity in pseudobinary random composition alloys, such as In 1−𝑥 ⁢Ga 𝑥 ⁢As, complicates such analyses. We present a simulation strategy to systematically take into account the variability in the local defect environment in order to predict statistical properties of neutral intrinsic defects in In 1−𝑥⁢ Ga 𝑥 ⁢As. We use a comprehensive sampling from a modest-sized 64-atom special quasirandom structure (SQS) to define a statistically representative set of defects, and use a 512-atom hypercell, a 2 × 2 × 2 supercell of SQS supercells, to achieve cell-size convergence. We articulate an equivalent site principle and describe how it constrains atomic chemical reference energies in computation of defect formation energies in pseudobinary alloys. A simple protocol for estimating reference energies for the Ga and In atoms sharing the III site succeeds in obtaining the equivalence of defects at Ga-sites and In sites in the SQS supercell, (<30 meV differences in average formation energies). For III-site defects, such as the As antisite As III , the statistical variability in formation energies is modest, ≈ 0.1–0.2 eV. The variability in formation energy at As-site defects, such as the As vacancy 𝑣 As , can be much larger, >1 eV. The As antisite is shown to be a low-energy defect and the most likely to be present in as-grown materials, just as in GaAs. All other defects are higher-energy defects unlikely to be important in native material, but potentially important in radiation-damaged material. With a strong variability in defect energies, especially on the As-site, explicit consideration of statistical variability due to compositional randomness will be imperative for meaningful and quantitative comparisons to experiment.

Density functional theory↗

DOE OSTI2021

Sensitivity Analysis of MFiX-PIC Parameters Using Nodeworks, PSUADE, and DAKOTA

The study presented in this report was aimed to demonstrate UQ analysis performed not only with Nodeworks, but also two other well-established UQ software tools from the U.S. DOE’s National Laboratories (PSUADE from Lawrence Livermore National Laboratory and DAKOTA from Sandia National Laboratory). It is important to emphasize that the motivation of this study was not to determine the best UQ software, but to verify if the global sensitivity analyses from the end-to-end workflow in Nodeworks are consistent with the results of other two UQ software. The components of Nodeworks from Python’s ecosystem have been tested as standalone libraries. However, an assessment study for the complete workflow targeting a specific UQ analysis has not been performed for Nodeworks. Hence, this study is expected to serve as an equivalent of solution verification for Nodeworks using other established UQ tools as reference solution. For this purpose, three distinct flow configurations (i.e., settling bed, bubbling fluidized, and circulating fluidized bed) have been used as representative multiphase flow problems of interest. The results of the systematic simulation campaigns performed in an earlier study using the particle-in-cell (PIC) approach in the Multiphase Flow with Interphase eXchanges (MFIX) suite of solvers (i.e., MFiX-PIC) was utilized. The same set of tabulated results was provided as input to the different UQ software for global sensitivity analysis. Results for the three cases indicate that based on the Sobol’ Sensitivity Indices method the order of importance ranking determined by Nodeworks for the Sobol’ Total Sensitivity Indices is consistent with PSUADE and DAKOTA in each case for the five model parameters considered. The input files for Nodeworks for the three cases are also shared through NETL’s Gitlab repository for the reader interested in reproducibility and further analysis (See Section 1.2).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Subject-specific multi-scale modeling of the fate of inhaled aerosols

Determining the fate of inhaled aerosols in the respiratory system is essential in assessing the potential toxicity of inhaled airborne materials, responses to airborne pathogens, or in improving inhaled drug delivery. The availability of high-resolution clinical lung imaging and advances in the reconstruction of lung airways from CT images have led to the development of subject-specific in-silico 3D models of aerosol dosimetry, often referred to as computational fluid-particle-dynamics (CFPD) models. As CFPD models require extensive computing resources, they are typically confined to the upper and large airways. These models can be combined with lower-dimensional models to form multiscale models that predict the transport and deposition of inhaled aerosols in the entire respiratory tract. Understanding where aerosols deposit is only the first of potentially several key events necessary to predict an outcome, being a detrimental health effect or a therapeutic response. To that end, multiscale approaches that combine CFPD with physiologically-based pharmacokinetics (PBPK) models have been developed to evaluate the absorption, distribution, metabolism, and excretion (ADME) of toxic or medicinal chemicals in one or more compartments of the human body. CFPD models can also be combined with host cell dynamics (HCD) models to assess regional immune system responses. Here, this paper reviews the state of the art of these different multiscale approaches and discusses the potential role of personalized or subject-specific modeling in respiratory health.

60 APPLIED LIFE SCIENCES↗