Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “component 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,693 records · Page 94

Advancing Multi-Hazard Risk and Safety Considerations for Aging Nuclear Facilities (Revision 1)

This project will demonstrate a multi-hazard time-dependent probabilistic risk assessment (PRA) approach for nuclear facilities considering aging-related deterioration of structures. A generic pressurized water reactor (PWR) reactor subjected to seismic mainshock-aftershock sequences considering the aging of the containment structure will be used as a case study to demonstrate the multi-hazard PRA approach. Using advanced modeling and simulation, seismic mainshock-aftershock fragility functions will be simulated for the containment structure considering aging effects. A multi-hazard PRA model for a generic PWR reactor will be built to quantify the multi-hazard core damage frequency (CDF) and large early release frequency (LERF) with explicit time-dependent modeling of event sequences. To date, the cascading impacts of multi-hazards are not adequately accounted for in the PRA models for nuclear facilities. In addition, for both initial and periodic evaluation of facilities to withstand natural phenomena hazards (NPH), deterioration of the SSCs due to aging and other effects is not adequately considered. By advancing multi-hazard considerations accounting for aging effects, this project will contribute to improved understanding of the safety of aging nuclear facilities. Project outcomes such as the multi-hazard CDF and LERF will also allow facility owners to optimize upgrade and retrofit protocols. This project is divided into two components: (1) advanced modeling and simulations; and (2) multi-hazard PRA. For the first component, Idaho National Laboratory’s (INL) Multi-hazard Analysis for STOchastic time-DOmaiN phenomena (MASTODON) and BlackBear codes will be used to simulate the seismic response and damage of a containment structure under cascading mainshock-aftershock sequences with aging effects. Uncertainties related to the seismic inputs, material parameters, and environmental factors such as temperature and humidity will be identified and propagated to the fragility functions. These fragility functions, which consider aging effects, will be time dependent. The second component will take the fragility information from the first component to build a multi-hazard time-dependent PRA model starting from a generic PWR model to evaluate the multi-hazard CDF and LERF and characterize the associated consequences. For this component, OpenPRA Web Application and SAPHIRE code will be used. Events such as the Fukushima Daiichi accident have highlighted the importance of considering the cascading impacts of multi-hazards for PRA. Moreover, many of the reactors in the current nuclear fleet in the United States (US) are already operating well beyond their initially planned design life, and applications for further extensions to operating licenses are being considered. Therefore, considering multi-hazard effects and aging deterioration in the NPH risk assessment process will contribute to the safety of both existing and future nuclear facilities. Outcomes of this project will thus directly benefit the standards DOE-STD-1020-16 and DOE-HDBK-1224-18.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A model for simulation and processing of radar images

A model for recording, processing, presentation, and analysis of radar images in digital form is presented. The observed image is represented as having two random components, one which models the variation due to the coherent addition of electromagnetic energy scattered from different objects in the illuminated areas. This component is referred to as fading. The other component is a representation of the terrain variation which can be described as the actual signal which the radar is attempting to measure. The combination of these two components provides a description of radar images as being the output of a linear space-variant filter operating on the product of the fading and terrain random processes. In addition, the model is applied to a digital image processing problem using the design and implementation of enhancement scene. Finally, parallel approaches are being employed as possible means of solving other processing problems such as SAR image map-matching, data compression, and pattern recognition.

Stiles, J. A.↗

The multicomponent structure of bulges

The morphology of disk galaxies is usually described by two major components, the centrally concentrated spheroidal component, called the bulge, and an oblate disk. The ratio of there contribution to the total luminosity - the bulge-to-disk ratio - is one of the parameters characterizing the Hubble sequence. Following de Vaucouleurs (1948), for most galaxies the radial distribution of the outer spheroid is fairly well described by the r exp 1/4 law I(r)=I(sub 0) exp(-(alpha)r), whereas the radial luminosity distribution of the disk follows an exponential law: I(r)=I(sub 0) exp(-alpha(r exp 1/4)) (Freeman 1970), with r the radial distance from the center. I(sub 0) and alpha are characteristic constants for each individual galaxy. Parameters for the structural properties of these components give important constraints for models of the formation and evolution of galaxies. Therefore we have tried to decompose disk and bulge components from high S/N CCD observations of a sample of edge-on disk galaxies. A common procedure for the decomposition is to model one component in a region where it dominates and subtract it from the combined light distribution. This technique was successfully carried out e.g. by van der Kruit & Earl (1981, 1982) and Wakamatsu & Hamabe (1984, 1989). Here we present two more examples of bulge-dominated edge-on SO galaxies, namely ESO 506-G33 and NGC 7123, which show an additional small and concentrated central component besides disk and 'bulge'.

Dettmar, Ralf-Juergen↗

Synthesis of Correct Digital Controller Models from Specifications by Model Transformation (21-0320)

The design of high consequence controllers (in weapons systems, autonomy, etc.) that do what they are supposed to do is a significant challenge. Testing simply does not come close to meeting the requirements for assurance. Today circuit designers at Sandia (and elsewhere) typically capture the core behavior of their components using state models in tools such as STATEFLOW. They then check that their models meet certain requirements (e.g. “The system bus must not deadlock” or “both traffic lights at an intersection must not be green at the same time”) using tools called model checkers. If the model checker returns “yes” then the property is guaranteed to be satisfied by the model. However, there are several drawbacks to this industry practice: (1) there is a lot of detail to get right, this is particularly challenging when there are multiple components requiring complex coordination (2) any errors returned by the model checker have to be traced back through the design and fixed, necessitating rework, (3) there are severe scalability problems with this approach, particularly when dealing with concurrency. All this places high demands on the designers who now face not only an accelerated schedule but also controllers of increasing complexity. This report describes a new and fundamentally different approach to the construction of safety-critical digital controllers. Instead of directly constructing a complete model and then trying to verify it, the designer can start with an initial abstract (think “sketch”) model plus the requirements, from which a correct concrete model is automatically synthesized. There is no need for post-hoc verification of required functional properties. Having tool to carry this out will significantly impact the nation’s ability to ensure the safety of high-consequence digital systems. The approach has been implemented in a prototype tool, along with a suite of examples, including ones that reflect actual problems faced by designers. Our approach operates on a variant of Statecharts developed at Sandia called Qspecs. Statecharts are a widely used formalism for developing concurrent reactive systems, supporting scalability through allowing state models containing composite states, which are the serial or parallel composition of substates which can themselves contain statecharts. Statecharts enable an incremental style of development, in which states are progressively refined to incorporate greater detail in an incremental model of software development. Our approach formulates a set of constraints from the structure of the models and the requirements and propagates these constraints to a fixpoint. The solution to the constraints is an inductive invariant along with guards on the transitions. We also show how our approach extends to implementation refinement, decomposition, composition, and elaboration. We currently handle safety requirements written in LTL (Linear Temporal Logic)

42 ENGINEERING↗

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗

The Community Land Model Version 5: Description of New Features, Benchmarking, and Impact of Forcing Uncertainty

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

58 GEOSCIENCES↗

Apollo cryogenic integrated systems program

The integrated systems program is capable of simulating both nominal and anomalous operation of the Apollo cryogenics storage system (CSS). Two versions of the program exist; one for the Apollo 14 configuration and the other for J Type Mission configurations. The program consists of two mathematical models which are dynamically coupled. A model of the CSS components and lines determines the oxygen and hydrogen flowrate from each storage tank given the tank pressures and temperatures, and the electrical power subsystem and environmental control subsystem flow demands. Temperatures and pressures throughout the components and lines are also determined. A model of the CSS tankage determines the pressure and temperatures in the tanks given the flowrate from each tank and the thermal environment. The model accounts for tank stretch and includes simplified oxygen tank heater and stratification routines. The program is currently operational on the Univac 1108 computer.

Seto, R. K. M.↗

Quantitative risk assessment system (QRAS)

A quantitative risk assessment system (QRAS) builds a risk model of a system for which risk of failure is being assessed, then analyzes the risk of the system corresponding to the risk model. The QRAS performs sensitivity analysis of the risk model by altering fundamental components and quantifications built into the risk model, then re-analyzes the risk of the system using the modifications. More particularly, the risk model is built by building a hierarchy, creating a mission timeline, quantifying failure modes, and building/editing event sequence diagrams. Multiplicities, dependencies, and redundancies of the system are included in the risk model. For analysis runs, a fixed baseline is first constructed and stored. This baseline contains the lowest level scenarios, preserved in event tree structure. The analysis runs, at any level of the hierarchy and below, access this baseline for risk quantitative computation as well as ranking of particular risks. A standalone Tool Box capability exists, allowing the user to store application programs within QRAS.

Weinstock, Robert M↗

3D Material Response Analysis of PICA Pyrolysis Experiments

Primarily interested in improving ablation modeling for use in inverse reconstruction of flight environments on ablative heat shields. Ablation model is essentially a component of the heat flux sensor, so model uncertainties lead to measurement uncertainties. Non-equilibrium processes have been known to be significant in low density ablators for a long time, but increased accuracy requirements of the reconstruction process necessitates incorporating this physical effect. Attempting to develop a pyrolysis model for implementation in material response based on the PICA data produced by Bessire and Minton. Pyrolysis gas species molar yields as a function of temperature and heating rate. Several problems encountered while trying to fit Arrhenius models to the data led to further investigation of the experimental setup.

Oliver, Brandon A.↗

Developing High-Throughput Organ-on-a-Chip Models to Investigate the Effects of Ionizing Radiation on the Central Nervous System

One of the main health risks in human space exploration is central nervous system (CNS) damage by ionizing radiation due to exposure to the galactic cosmic rays (GCRs). In animal models, irradiation with simulated GCRs or their components has been shown to cause neuronal damage and neuroinflammation associated with cognitive and behavioral dysfunction. In general, the extent of CNS damage is partially regulated by the blood-brain barrier (BBB), which enables immune cells to enter the CNS. The main cellular regulators of BBB permeability are astrocytes, which also modulate neuronal death, immune responses and oxidative stress, and thus could serve as a robust CNS-specific target for countermeasure development. However, studies on BBB permeability and astrocyte functions in regulating CNS responses to ionizing radiation have been limited, especially in human tissue/organ analogs. Therefore, we established a high-throughput 3D organ-on-a-chip system to study human CNS and BBB impairments in response to ionizing radiation, based on commercially available OrganoPlates (Mimetas, Inc.) seeded with primary or induced pluripotent stem cell-derived human cells. We investigated both immediate and delayed CNS responses to major GCR components: 0.15-0.5 Gy 250MeV/n 4-He, and 0.3-0.8 Gy 600 MeV/n 56-Fe; as well as to 0.5-1 Gy X-rays. We observed ionizing radiation-mediated increases in BBB permeability that was exacerbated by astrocyte presence and accompanied by morphological changes in endothelial cells and tight junctions, altered cytokine profile including TNFa upregulation, and increased oxidative stress. We also quantified irradiation-mediated changes in astrocyte activation and neuronal functions, revealing major astrocyte damage mediated by 600MeV/n 56-Fe particles. Thus, we demonstrate that deep space radiation may contribute to CNS damage by disrupting both astrocyte and endothelial cell components of the blood-brain barrier. Our next steps include mapping and validating the transcriptomic changes induced by simulated GCRs and their components in human CNS models. Ultimately, we aim to uncover potential novel targets for countermeasure developments to mitigate CNS damage in long duration spaceflight.

Radiation↗

WHONDRS-GUI: a web application for global survey of surface water metabolites

Background The Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) is a consortium that aims to understand complex hydrologic, biogeochemical, and microbial connections within river corridors experiencing perturbations such as dam operations, floods, and droughts. For one ongoing WHONDRS sampling campaign, surface water metabolite and microbiome samples are collected through a global survey to generate knowledge across diverse river corridors. Metabolomics analysis and a suite of geochemical analyses have been performed for collected samples through the Environmental Molecular Sciences Laboratory (EMSL). The obtained knowledge and data package inform mechanistic and data-driven models to enhance predictions of outcomes of hydrologic perturbations and watershed function, one of the most critical components in model-data integration. To support efforts of the multi-domain integration and make the ever-growing data package more accessible for researchers across the world, a Shiny/R Graphical User Interface (GUI) called WHONDRS-GUI was created. Results The web application can be run on any modern web browser without any programming or operational system requirements, thus providing an open, well-structured, discoverable dataset for WHONDRS. Together with a context-aware dynamic user interface, the WHONDRS-GUI has functionality for searching, compiling, integrating, visualizing and exporting different data types that can easily be used by the community. The web application and data package are available at https://data.ess-dive.lbl.gov/view/doi:10.15485/1484811 , which enables users to simultaneously obtain access to the data and code and to subsequently run the web app locally. The WHONDRS-GUI is also available for online use at Shiny Server ( https://xmlin.shinyapps.io/whondrs/ ).

59 BASIC BIOLOGICAL SCIENCES↗

On the simulation of the oceanic general circulation

Two global ocean simulations based on the physics of a highly viscous ocean are discussed, one having realistic atmospheric functions for calculating the thermal forcing of the ocean. The velocity field of this model compares reasonably well with the relatively small amount of real data available. Temperature and heat budget components of the model reproduce best the equatorial band of heating. A simulation based upon the physics of a weakly viscous ocean is described, which produces the correct pattern of isotherms. This model has not yet been run to thermal equilibrium. The Gulf Stream and eddies and their influence on the oceanic and atmospheric heat budgets are discussed in the terms of their importance in coupled ocean-atmosphere models.

Mintz, Y.↗

Animal models for simulating weightlessness

NASA has developed a rat model to simulate on earth some aspects of the weightlessness alterations experienced in space, i.e., unloading and fluid shifts. Comparison of data collected from space flight and from the head-down rat suspension model suggests that this model system reproduces many of the physiological alterations induced by space flight. Data from various versions of the rat model are virtually identical for the same parameters; thus, modifications of the model for acute, chronic, or metabolic studies do not alter the results as long as the critical components of the model are maintained, i.e., a cephalad shift of fluids and/or unloading of the rear limbs.

Morey-Holton, E.↗

Modeling Bird Migration under Climate Change: A Mechanistic Approach

How will migrating birds respond to changes in the environment under climate change? What are the implications for migratory success under the various accelerated climate change scenarios as forecast by the Intergovernmental Panel on Climate Change? How will reductions or increased variability in the number or quality of wetland stop-over sites affect migratory bird species? The answers to these questions have important ramifications for conservation biology and wildlife management. Here, we describe the use of continental scale simulation modeling to explore how spatio-temporal changes along migratory flyways affect en-route migration success. We use an individually based, biophysical, mechanistic, bird migration model to simulate the movement of shorebirds in North America as a tool to study how such factors as drought and wetland loss may impact migratory success and modify migration patterns. Our model is driven by remote sensing and climate data and incorporates important landscape variables. The energy budget components of the model include resting, foraging, and flight, but presently predation is ignored. Results/Conclusions We illustrate our model by studying the spring migration of sandpipers through the Great Plains to their Arctic breeding grounds. Why many species of shorebirds have shown significant declines remains a puzzle. Shorebirds are sensitive to stop-over quality and spacing because of their need for frequent refueling stops and their opportunistic feeding patterns. We predict bird "hydrographs that is, stop-over frequency with latitude, that are in agreement with the literature. Mean stop-over durations predicted from our model for nominal cases also are consistent with the limited, but available data. For the shorebird species simulated, our model predicts that shorebirds exhibit significant plasticity and are able to shift their migration patterns in response to changing drought conditions. However, the question remains as to whether this behavior can be maintained over increasing and sustained environmental change. Also, the problem is much more complex than described by the current processes captured in our model. We have taken some important and interesting steps, and our model does demonstrate how local scale information about individual stop-over sites can be linked into the migratory flyway as a whole. We are incorporating additional, species specific, mechanistic processes to better reflect different climate change scenarios

Smith, James A.↗

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗

Microgravity acceleration modeling for orbital systems

In view of the decisive importance of a disturbance-free environment on the Space Station, and on other orbital systems, for materials processing experiments, a theoretical and semi-experimental analysis of the acceleration environment to be expected on large orbiting spacecraft was undertaken. A unified model of such spacecraft cannot be established; therefore, a number of sub-models representing major components of typical large spacecraft must be investigated. In order to obtain experimental data of forces, a typical spacecraft - an engineering model of the Spacelab - was suspended on long ropes in a high-bay hangar, and equipped with a number of accelerometers. Active components on the Spacelab (fans, pumps, air conditioners, valves, levers) were operated, and astronautics moved boxes, drawers, sleds, and their own bodies. Generally speaking, the response of the Spacelab structure was very similar to the environment measured on Spacelabs SL-1, SL-2, and D-1. At frequencies in the broad range between 1 and about 100 Hz, acceleration peaks reached values of 10(exp -3) and 10(exp -2) g sub o, and even higher.

Knabe, Walter↗

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗