Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional Data Analysis”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Analysis of phase stability and chemical segregation in the Mo-V alloys using a generalized embedded atom method potential

A new interatomic potential for the Mo-V system is introduced to facilitate the study of phase stability and mechanical properties at lower temperatures. This potential is based on a generalization of the embedded atom method and includes contributions from embedding energy, explicit two- and three-body interactions and nonlocal many-body interaction terms. The parameters of the potential are optimized by using data from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, such as elastic constants, equation of states, phonon dispersion curves, point defect properties and melting temperatures for different compositions. Even though our potential is trained on a small dataset, its accuracy is comparable to available machine learning potentials for Mo and V. Furthermore, our results show that an ordered B2 phase is stable at low temperatures in alloys containing 50% V, but the solid solution phase is stable above 800 K. However, such long-range ordering is not observed in V-rich or Mo-rich alloys. In addition, our results show that V segregates to dislocation cores and grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Jensen–Shannon divergence based novel loss functions for Bayesian neural networks

Bayesian neural networks (BNNs) are state-of-the-art machine learning methods that can naturally regularize and systematically quantify uncertainties using their stochastic parameters. Kullback–Leibler (KL) divergence-based variational inference used in BNNs suffer from unstable optimization and challenges in approximating light-tailed posteriors due to the unbounded nature of the KL divergence. To resolve these issues, we formulate a novel loss function for BNNs based on a new modification to the generalized Jensen–Shannon (JS) divergence, which is bounded. In addition, we propose a Geometric JS divergence-based loss, which is computationally efficient since it can be evaluated analytically. We found that the JS divergence-based variational inference is intractable, and hence employed a constrained optimization framework to formulate these losses. Our theoretical analysis and empirical experiments on multiple regression and classification data sets suggest that the proposed losses perform better than the KL divergence-based loss, especially when the data sets are noisy or biased. Specifically, there are approximately 5% and 8% improvements in accuracy for a noise-added CIFAR-10 dataset and a regression dataset, respectively. There is about 13% reduction in false negative predictions of a biased histopathology dataset. Additionally, we quantify and compare the uncertainty metrics for the regression and classification tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Inelastic Neutron Scattering Study of Magnetic Exchange Pathways in MnS

We report an investigation of the magnetic structure and magnetic exchange pathways in MnS via neutron scattering methods, aided by density functional theory (DFT) modeling. The material has been confirmed to undergo antiferromagnetic (AFM) ordering at 152 K, with the magnetic structure representing AFM stacking of ferromagnetic (FM) (111) planes of Mn magnetic moments. Correspondingly, the magnetic structure is described by a propagation vector k = (1/2, 1/2, 1/2), with the volume of the magnetic unit cell being 8 times larger than the volume of the nuclear unit cell. Analysis of inelastic neutron scattering (INS) data collected on a powder sample of MnS revealed that the next-nearest-neighbor magnetic exchange constant (J 2 ) exceeds the nearest-neighbor exchange constant (J 1 ) by more than 3 times, while in the case of MnO, which exhibits the same nuclear and magnetic structures as MnS, the J 2 /J 1 ratio was reported to be below 1.5. Although for MnO the signs of both J 1 and J 2 indicated AFM exchange interactions, machine-learning INS data analysis in combination with DFT calculations suggests that the INS data collected on MnS are best described with J 1 < 0 and J 2 > 0, corresponding to FM and AFM exchange couplings, respectively. To achieve a satisfactory fit to the experimentally observed data, the Hamiltonian used to model the INS spectra also included the next-next-nearest-neighbor magnetic exchange constant (J 3 ). The best-fit model has been obtained with the values of the exchange constants J 1 = -0.27 meV, J 2 = 1.05 meV, and J 3 = -0.19 meV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The Exoproteome and Surfaceome of Toxigenic Corynebacterium diphtheriae 1737 and Its Response to Iron Restriction and Growth on Human Hemoglobin

Toxin-producing Corynebacterium diphtheriae strains are the etiological agents of the severe upper respiratory disease, diphtheria. A global phylogenetic analysis revealed that biotype gravis is particularly lethal as it produces diphtheria toxin and a range of other virulence factors, particularly when it encounters low levels of iron at sites of infection. Here, to gain insight into how it colonizes its host, we have identified iron-dependent changes in the exoproteome and surfaceome of C. diphtheriae strain 1737 using a combination of whole-cell fractionation, intact cell surface proteolysis, and quantitative proteomics. In total, we identified 1414 of the predicted 2265 proteins (62%) encoded by its reference genome. For each protein, we quantified its degree of secretion and surface exposure, revealing that exoproteases and hydrolases predominate in the exoproteome, while the surfaceome is enriched with adhesins, particularly DIP2093. Our analysis provides insight into how components in the heme-acquisition system are positioned, showing pronounced surface exposure of the strain-specific ChtA/ChtC paralogues and high secretion of the species-conserved heme-binding HtaA protein, suggesting it functions as a hemophore. Profiling the response of the exoproteome and surfaceome after microbial exposure to human hemoglobin and iron limitation reveals potential virulence factors that may be expressed at sites of infection. Data are available via ProteomeXchange with identifier PXD051674.

cell envelope↗
DOE OSTI2023

Tripogon loliiformis tolerates rapid desiccation after metabolic and transcriptional priming during initial drying

Abstract Crop plants and undomesticated resilient species employ different strategies to regulate their energy resources and growth. Most crop species are sensitive to stress and prioritise rapid growth to maximise yield or biomass production. In contrast, resilient plants grow slowly, are small, and allocate their resources for survival in challenging environments. One small group of plants, termed resurrection plants, survive desiccation of their vegetative tissue and regain full metabolic activity upon watering. However, the precise molecular mechanisms underlying this extreme tolerance remain unknown. In this study, we employed a transcriptomics and metabolomics approach, to investigate the mechanisms of desiccation tolerance in Tripogon loliiformis , a modified desiccation-tolerant plant, that survives gradual but not rapid drying. We show that T. loliiformis can survive rapid desiccation if it is gradually dried to 60% relative water content (RWC). Furthermore, the gene expression data showed that T. loliiformis is genetically predisposed for desiccation in the hydrated state, as evidenced by the accumulation of MYB, NAC, bZIP, WRKY transcription factors along with the phytohormones, abscisic acid, salicylic acid, amino acids (e.g., proline) and TCA cycle sugars during initial drying. Through network analysis of co-expressed genes, we observed differential responses to desiccation between T. loliiformis shoots and roots. Dehydrating shoots displayed global transcriptional changes across broad functional categories, although no enrichment was observed during drying. In contrast, dehydrating roots showed distinct network changes with the most significant differences occurring at 40% RWC. The cumulative effects of the early stress responses may indicate the minimum requirements of desiccation tolerance and enable T. loliiformis to survive rapid drying. These findings potentially hold promise for identifying biotechnological solutions aimed at developing drought-tolerant crops without growth and yield penalties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Characterization of contaminants in the Lyman-alpha forest auto-correlation with DESI

Baryon Acoustic Oscillations can be measured with sub-percent precision above redshift two with the Lyman-α (Lyα) forest auto-correlation and its cross-correlation with quasar positions. This is one of the key goals of the Dark Energy Spectroscopic Instrument (DESI) which started its main survey in May 2021. We present in this paper a study of the contaminants to the Lyα forest which are mainly caused by correlated signals introduced by the spectroscopic data processing pipeline as well as astrophysical contaminants due to foreground absorption in the intergalactic medium. Notably, an excess signal caused by the sky background subtraction noise is present in the Lyα auto-correlation in the first line-of-sight separation bin. We use synthetic data to isolate this contribution, we also characterize the effect of spectro-photometric calibration noise, and propose a simple model to account for both effects in the analysis of the Lyα forest. We then measure the auto-correlation of the quasar flux transmission fraction of low redshift quasars, where there is no Lyα forest absorption but only its contaminants. We demonstrate that we can interpret the data with a two-component model: data processing noise and triply ionized Silicon and Carbon auto-correlations. This result can be used to improve the modeling of the Lyα auto-correlation function measured with DESI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Magnetism of Ir 5 + -based double perovskites: Unraveling its nature and the influence of structure

A combination of x-ray-diffraction, x-ray absorption, x-ray magnetic circular dichroism, macroscopic magnetization, and muon-spin relaxation measurements is used to investigate the interplay between structure (both crystallography and microstructure) and magnetism in Sr 2 YIrO 6 as a function of both chemical (Ca-doping) and physical (hydrostatic) pressure. X-ray absorption spectroscopy clearly shows that the physical pressure is more effective in modifying the structure. On the other hand, the dichroic measurements evidence a constant magnetic signal with physical pressure and strong differences with Ca-doping. Muon-spin relaxation reveals the presence of magnetic order, even when this is hidden in the magnetization data. From the combined analysis, the magnetic results are explained in terms of the presence of Ir 6+ or Ir 4+ ions in magnetic clusters, most likely located at and triggered by antisite disorder. The measurements under high physical pressure indicate that the magnetic state is independent of the crystallographic details such as Ir-O distances and Ir-O-Y angles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

High-Quality Revision of the Israeli Seismic Bulletin

Seismic bulletins, with trustworthy phase picks, origin times, and source locations are key for regional seismic studies, such as travel-time (TT) tomography, attenuation tomography, and anisotropy studies. To lay the groundwork for such studies in Israel, we revised the seismic bulletin of Israel and the surrounding area and obtained a trustworthy TT data set. From the earthquake and explosion bulletins of the Geophysical Institute of Israel, we compiled a starting data set of about 123,000 earthquakes and explosions that occurred during the past 40 yr. After screening out the poorly recorded events, we were left with a data set of ~38,000 well-recorded events. We then revised the remaining data set in two consecutive steps. In the first, we reviewed and updated station metadata, including changes in station metadata parameters over time. In the second step, we jointly relocated a list of selected seismic events, using the Bayesian hierarchical location software package (BayesLoc) of Myers et al. (2007) that performs joint relocation of multiple events. We observed striking dissimilarities between the spatial distributions of the newly relocated catalog and the initial locations. Although the depth distribution of the starting catalog is trimodal with peaks at 0, 5, and 10 km, the distribution in this study is unimodal, with a broad peak between 7.5 and 12.5 km. By differencing the observed arrival times and the origin times obtained through relocation with BayesLoc, we obtained a revised TT database that consists of 261,336 Pg, 132,876 Pn, 114,816 Sg, and 60,394 Sn arrivals, from a set of 30,458 jointly relocated seismic sources. In this work, we compared prerevision and postrevision TTs as a function of epicentral distance and concluded that the revised data set contains far fewer outliers and inconsistencies than the original data set. The revised TT data set may be used for seismic studies, such as TT tomography, attenuation tomography, and anisotropy studies.

58 GEOSCIENCES↗
DOE OSTI2025

FRAM v.7.1

Fixed-energy Response-function Analysis with Multiple efficiency (FRAM) is a software code designed primarily for plutonium and uranium isotopic analysis. It is widely used in both the domestic and international safeguards community. FRAM can quickly and accurately determine the isotopic compositions of plutonium, uranium, and mixed oxides (MOX) when measuring with a high-purity germanium (HPGe), cadmium zinc telluride (CZT), or lanthanum bromide (LaBr3) detector. The capabilities of FRAM have been enhanced to analyze the data of the pixelated CZT detector (made by H3D) and to measure the mass of plutonium, uranium, and MOX. Both the isotopic composition and mass of the item can be quickly determined with one measurement using a gamma detector with FRAM v.7.1.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Robust Dark Energy Constraints with the Dark Energy Spectroscopic Survey (Final Technical Report)

This project developed and applied advanced theoretical, computational, and data-analysis methodologies to extract robust and precise cosmological constraints from the Dark Energy Spectroscopic Instrument (DESI). The work focused on maximizing the scientific return of DESI through optimized survey strategy, novel higher-order clustering statistics, improved modeling of small-scale structure, and rigorous mitigation of observational systematics. Over the award period, the project made substantial contributions to DESI science planning, produced new methods for bispectrum and three-point correlation function analyses, advanced constraints on primordial non-Gaussianity, and delivered widely used software tools. The project also played a major role in training graduate students and a postdoctoral researcher who contributed directly to DESI key projects. The results have significantly enhanced the cosmological reach of DESI and provide a strong foundation for future surveys such as DESI-II and Stage-V experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Discovery top-down proteomics in symbiotic soybean root nodules

Proteomic methods have been widely used to study proteins in complex biological samples to understand biological molecular mechanisms. Most well-established methods (known as bottom-up proteomics, BUP) employ an enzymatic digestion step to cleave intact proteins into smaller peptides for liquid chromatography (LC) mass spectrometry (MS) detection. In contrast, top-down proteomics (TDP) directly characterizes intact proteins including all possible post-translational modifications (PTMs), thus offering unique insights into proteoform biology where combinations of individual PTMs may play important roles. We performed TDP on soybean root nodules infected by the symbiotic Bradyrhizobium japonicum in both the wildtype bacterium and a nifH- mutant, which lacks the ability to fix nitrogen in the soybean root nodule. TDP captured 1648 proteoforms derived from 313 bacterial genes and 178 soybean genes. Leghemoglobin, the most abundant protein in the sample, existed in many truncated proteoforms. Interestingly, these truncated proteoforms were considerably more abundant in the wildtype relative to the nifH- mutant, implicating protease activity as an important factor in nitrogen fixation. Proteoforms with various PTMs and combinations thereof were identified using an unrestricted open modification search. This included less common PTMs such as myristoylation, palmitoylation, cyanylation, and sulfation. In parallel, we collected high resolution MS imaging (MSI) data of intact proteins and biopolymers (<20 kDa due to current technical limitations) from sections of the soybean root nodules using matrix-assisted laser desorption/ionization (MALDI) coupled to high resolution Orbitrap. Several detected proteoforms exhibited unique spatial distributions inside the infection zone and cortex, suggesting functional compartmentalization in these regions. A subset of peaks from the MALDI-MSI were assigned to proteoforms detected in TDP LCMS data based on matching accurate masses. Many of the proteins detected in both LCMS and MALDI-MSI are currently uncharacterized in UniProt: the PTM and spatial information presented here will be valuable in understanding their biological functions. Taken together, our study demonstrates how untargeted TDP approach can provide unique insights into plant proteoform biology. On-going technology developments are expected to further improve TDP coverage for more comprehensive high-throughput analysis of proteoforms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Thermal scattering law for structure-dependent-Doppler broadening in FLASSH

Traditional Doppler broadening methods assume that target nuclei velocities follow a Maxwell-Boltzmann distribution which inherently assumes that the atoms are a free gas or that their velocities are independent of each other. This approximation is reasonable at high target temperatures and high neutron energies. However, nuclei are not independent: the lattice temperature and chemical binding will define the atomic motion. The thermal scattering law (i.e., TSL or S(α,β)) is a material property which describes the chemical binding and temperature response in terms of available momentum and energy states. In the thermal energy range, this TSL will define the thermal scattering cross sections. In the resonance region, the TSL offers a first-principles description of the probability distribution function for the velocity of the target nuclei. Using the TSL in Doppler broadening captures the structure of the material consistently from the thermal range into the resonance cross sections. In this work, both traditional free-gas and precise TSL Doppler broadening methods for resonance cross section evaluation have been implemented within the Full Law Analysis Scattering System Hub (FLASSH). This framework provides a generalized formulation for evaluating low-lying resonance data and streamlines nuclear data processing methods. TSL data generated using ab initio lattice dynamics (AILD) can be used to broaden cross section libraries in ENDF format for high-fidelity input into reactor physics calculations. These capabilities are demonstrated for the lowest absorption resonances of {sup 238}U in UO{sub 2}, UC and UN. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Data for A Fluorescence-Based Transient Expression Assay for the Analysis of Upstream Open Reading Frames in Plants

Scripts for the manuscript "A fluorescence-based transient expression assay for the analysis of upstream open reading frames in plant" by Haas et al. Upstream open reading frames (uORFs) are regulatory elements present in the 5′ leaders of mRNA that can significantly impact downstream gene expression in eukaryotes. In crop engineering, editing of uORFs can provide an avenue to upregulate expression of native genes without the need to add persistent transgenic copies. Even with genome- wide methods to identify translated uORFs such as ribosome profiling, their functional characterization depends on validation through reporter gene assays and mutagenesis studies. Current screening methods for plants use luciferases or protoplasts to measure differential gene expression between wild- type and mutated transcript leaders, which requires tissue processing and/or substrate addition. Here, we present a time- and cost- efficient alternative to investigate transcript leaders by co- expression of two fluorescent proteins in Nicotiana benthamiana leaf tissue and test our assay on genes involved in photoprotection, editing of which could provide a pathway to increase CO2 assimilation during sun–shade transitions.

Gene Editing↗
DOE OSTI2021

COVID-19 dynamics across the US: A deep learning study of human mobility and social behavior

This paper presents a deep learning framework for epidemiology system identification from noisy and sparse observations with quantified uncertainty. The proposed approach employs an ensemble of deep neural networks to infer the time-dependent reproduction number of an infectious disease by formulating a tensor-based multi-step loss function that allows us to efficiently calibrate the model on multiple observed trajectories. The method is applied to a mobility and social behavior-based SEIR model of COVID-19 spread. The model is trained on Google and Unacast mobility data spanning a period of 66 days, and is able to yield accurate future forecasts of COVID-19 spread in 203 US counties within a time-window of 15 days. Interestingly, a sensitivity analysis that assesses the importance of different mobility and social behavior parameters reveals that attendance of close places, including workplaces, residential, and retail and recreational locations, has the largest impact on the effective reproduction number. Furthermore, the model enables us to rapidly probe and quantify the effects of government interventions, such as lock-down and re-opening strategies. Taken together, the proposed framework provides a robust workflow for data-driven epidemiology model discovery under uncertainty and produces probabilistic forecasts for the evolution of a pandemic that can judiciously provide information for policy and decision making. All codes and data accompanying this manuscript are available at https://github.com/PredictiveIntelligenceLab/DeepCOVID19.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Effects of the Grazing Incidence Geometry on X-ray Photon Correlation Spectroscopy Measurements

X-ray photon correlation spectroscopy (XPCS) is a versatile tool to measure dynamics on the nanometer to micrometer scale in bulk samples. XPCS has also been applied in grazing incidence (GI) geometry to examine the dynamics of surface layers. However, considering GI scattering experiments more universally, the GI geometry leads to a superposition of signals due to reflection and refraction effects, also known from the distorted-wave Born approximation (DWBA). In this paper, the impact of these reflection and refraction effects on the correlation analysis is determined experimentally by measuring grazing incidence transmission XPCS (GT-XPCS) and grazing incidence XPCS (GI-XPCS) simultaneously for a thin film sample, showing non-equilibrium dynamics. Here the results of the GI and GT geometry comparisons are combined within the framework of the standardly applied, simplified DWBA. These calculations allow identifying the main contributions of the detected signal from the leading scattering terms along the out-of-plane direction q z , which dominate the measured intensity pattern on the detector. In combination with the calculation of the non-linear effect of refraction in GTSAXS and GISAXS, it is possible to identify experimental conditions that can be chosen to run experiments and data analysis as close as possible to transmission XPCS and to explain which limitations for data interpretations are observed. Consequently, the beam exposure can be significantly reduced by using GI geometry only. Calculations of experimental settings prior to experiments are detailed to determine suitable q z regions for a variety of material systems measured in bulk-sensitive GI-XPCS experiments, allowing us to determine the scaling behavior of typical decay times as a function of q that is comparable to the scaling behavior obtained in distortion-free GT-XPCS or transmission XPCS experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Determination of proton PDF uncertainties with Markov chain Monte Carlo

We present an analysis of parton distribution functions (PDFs) of the proton using Markov chain Monte Carlo (MCMC) methods. The MCMC approach naturally implements Bayes’ theorem and, thus, provides a means to directly sample the underlying probability distribution—in this case, the probability distribution of the PDF parameters. This allows for a straightforward propagation of the resulting uncertainties into any PDF-dependent observable, preserving their simple probabilistic interpretation. In our analysis we include a broad set of deep inelastic scattering data from HERA, BCDMS and NMC experiments along with the Drell-Yan, 𝑊 and 𝑍 boson data from LHC and Tevatron experiments, which combined with theoretical calculations at next-to-next-to-leading order in QCD allow for realistic determination of PDFs. The main focus of this analysis is to explore alternative methods for PDF uncertainty estimation that are more firmly grounded in statistical principles. We show that the flexibility of the Bayes framework, allowing one, e.g., to account for non-Gaussianity or inconsistencies of datasets, is crucial to extract realistic uncertainties when such assumptions are not fulfilled. We also demonstrate that MCMC allows one to determine the Δ⁢𝜒 2 value corresponding to a given confidence level in the sample, which can, in turn, be used as a statistically well-founded tolerance criterion used in the Hessian method, thus addressing one of its main long-standing drawbacks.

Risse, Peter Clemens [Universität Münster (Germany↗
DOE OSTI2021

Understanding the dynamics of vacuole trafficking and secretion for enhanced plant storage reserves

Vacuoles are essential in plants and serve roles in both protein and organelle degradation for reutilization and storage of high energy molecules such as storage proteins which are critical for seed germination under normal and stress conditions and impact biomass production and nutrition. In addition, plants secrete many different molecules such as cell wall components, lipids, waxes and others that are important for plant protection, structure and biomass. Vacuole accumulation and secretion are closely linked and highly regulated pathways within the endomembrane trafficking network. We have discovered many drug-like synthetic small molecules that can suppress vacuole trafficking and enhance secretion; other chemicals can enhance vacuole trafficking and suppress secretion. The ability of our compounds to modulate these pathways indicates that these pathways are tightly coupled and highly regulated. Thus, to achieve the goals of enhanced energy reserves or enhanced secretion it is essential to understand how this dynamic is controlled. Bioactive chemicals have provided us unique insight into natural mechanisms controlling the dynamics of vacuole trafficking and secretion that was not possible previously. One such compound, ES2, negatively modulates secretion while enhancing vacuole trafficking. This is the first compound discovered that targets an essential secretory complex known as the exocyst and will reveal important regulatory domains controlling secretion and vacuole trafficking. During the current funding period, we determined the structure of the EXO70A1 in plants and identifying its binding site. To understand the dynamics of vacuole trafficking it is key to have compounds that inhibit vacuole trafficking as well as enhancing it (ES2). We have used sophisticated image analysis to identify another novel compound ES17 that enhances secretion while negatively modulating vacuole trafficking (opposite of ES2). We have already made rapid progress. Our preliminary data indicates that ES17 very specifically targets a late stage of vesicle fusion with the tonoplast which is highly regulated. Thus, ES17 is an ideal compound to dissect the regulation of vacuole targeting dynamics. In addition, we have now synthesized a biotinylated version of ES17 that is fully functional, permitting the use of chemical proteomics to identify the cognate target with a high probability of success. We propose to investigate the ES17 target which can negatively modulate vacuole trafficking and combine this with knowledge of the ES2 target which can promote vacuole trafficking. Dissecting the mechanisms controlling the dynamics of secretion vs vacuole accumulation and the use of synthetic modulators will permit us to understand and regulate trafficking to enhance secretion to the plasma membrane and cell wall or the accumulation of biological energy reserves for increased biomass, nutritional value, stress resistance, and the secretion of specialized plant products.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

​​Hydrogen Transportation and Distributed Energy Systems Seismic Risk Assessment for Cascadia Subduction Zone Airport Facilities​

Portland International Airport in Oregon is exploring operating a fleet of 28 fuel cell electric buses to support airport operations and possibly provide backup power during outages. In this report, we evaluate the risk of hydrogen deployment at an airport considering the potential for seismic activity. We present a seismic risk assessment for pressurized piping for a generic hydrogen refueling station at the airport. Only pressurized piping was considered in the risk analysis because of its vulnerability and availability of fragility information unlike the other components. Seismic capacity of other components should be included as more information becomes available. To characterize the seismic hazard, we considered the initiating event frequency of a Cascadia Subduction Zone earthquake event based on data from the literature. Seismic stressors that translate to the site from an earthquake event are expressed as peak ground acceleration. Site disturbance is characterized as a function of soil conditions, represented by different shear wave velocity values, which affect how seismic waves propagate and impact structures. Pipe fragility curves as a function of seismic capacity were used to correlate ground acceleration and failure probability. These fragility characteristics were combined with site-specific soil conditions to calculate the probability of pipe rupture during a seismic event. It should be noted that multiple simultaneous failures due to an earthquake event as a common cause are not considered. Some epistemic uncertainties such as self-exciting shocks, aftershocks, and the time variant nature of the ground motion are not considered either.

08 HYDROGEN↗