Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system analysis module”

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 505 records · Page 28

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

GaN-based W-band receiver chip development for fusion plasma diagnostics

Millimeter-wave diagnostics have proven effective on various magnetic fusion devices worldwide, yet the formidable challenges posed by the harsh environments of future burning plasma devices, characterized by extreme temperatures, pressures, and radiation levels, remain a significant hurdle. To address these challenges, the utilization of wide bandgap Gallium Nitride (GaN)-based millimeter-wave diagnostics is a most promising solution for fusion reactor safety monitoring and control. A noteworthy W-band GaN-based system-on-chip receiver has been the demonstrated by employing HRL T3 40 nm GaN technology. This receiver chip, compactly designed with dimensions of 3 × 5 mm 2 , incorporates essential components such as the 75–110 GHz RF Low-Noise Amplifier (LNA), mixer, Intermediate Frequency (IF) amplifier, and Local Oscillator (LO) chain. This receiver chip will be packaged as a millimeter-wave receiver module and applied on the DIII-D National Fusion Facility, for fusion plasma edge shape monitoring for operational safety and dangerous disruption prediction. The laboratory measurement results have demonstrated suitable performance. Furthermore, this advancement is pivotal for accurate analysis of plasma behavior in the extreme conditions of burning plasma devices, driving progress in fusion research and technology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Automated prediction of lattice parameters from X-ray powder diffraction patterns

A key step in the analysis of powder X-ray diffraction (PXRD) data is the accurate determination of unit-cell lattice parameters. This step often requires significant human intervention and is a bottleneck that hinders efforts towards automated analysis. This work develops a series of one-dimensional convolutional neural networks (1D-CNNs) trained to provide lattice parameter estimates for each crystal system. A mean absolute percentage error of approximately 10% is achieved for each crystal system, which corresponds to a 100- to 1000-fold reduction in lattice parameter search space volume. The models learn from nearly one million crystal structures contained within the Inorganic Crystal Structure Database and the Cambridge Structural Database and, due to the nature of these two complimentary databases, the models generalize well across chemistries. A key component of this work is a systematic analysis of the effect of different realistic experimental non-idealities on model performance. It is found that the addition of impurity phases, baseline noise and peak broadening present the greatest challenges to learning, while zero-offset error and random intensity modulations have little effect. However, appropriate data modification schemes can be used to bolster model performance and yield reasonable predictions, even for data which simulate realistic experimental non-idealities. In order to obtain accurate results, a new approach is introduced which uses the initial machine learning estimates with existing iterative whole-pattern refinement schemes to tackle automated unit-cell solution.

42 ENGINEERING↗

The Power and Energy Storage Systems Toolbox -- PSTess v.1.0

The Power and Energy Storage Systems Toolbox (PSTess) is a MATLAB-based computing package for dynamic simulation and analysis of utility-scale battery storage systems. This codebase is a fork of the Power Systems Toolbox Version 3.0, developed at Rensselaer Polytechnic Institute (RPI) and Cherry Tree Scientific Software. While PSTess shares a common lineage with PST, it is a substantially different application. As the name implies, the main distinguishing characteristic of PSTess is its ability to model inverter-based energy storage systems (ESS). The model that enables this is called ess.m, and it serves the dual role of representing ESS operational constraints and the generator/converter interface. With PSTess, the generator/converter interface is modeled as a controllable current source with the ability to modulate both active and reactive current. The model ess.m permits four-quadrant modulation, which allows it to represent a wide variety of inverter-based resources beyond energy storage when paired with an appropriate supplemental control model. Examples include utility-scale photovoltaic (PV) power plants, type 4 wind plants, and static synchronous compensators (STATCOM). This capability is especially useful for modeling hybrid plants that combine energy storage with renewable resources. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-8322 O

Elliott, RyanT.↗

Bioinformatic teaching resources - for educators, by educators - using KBase, a free, user-friendly, open source platform

Over the past year, biology educators and staff at the Department of Energy Systems Biology Knowledgebase (KBase) initiated a collaborative effort to develop a curriculum for bioinformatics education. KBase is a free and easily accessible data science platform that integrates many bioinformatics resources into a graphical user interface built upon reproducible analysis notebooks. KBase held conversations with college and high school instructors to understand how KBase could potentially support their educational goals. These conversations morphed into a working group of biological and data science instructors that adapted the KBase platform to their curriculum needs, specifically around concepts in Genomics, Metagenomics, Pangenomics, and Phylogenetics. The KBase Educators Working Group developed modular, adaptable, and customizable instructional units. Each instructional module contains teaching resources, publicly available data, analysis tools, and markdown capability to tailor instructions and learning goals for each class. The online user interface enables students to conduct hands-on data science research and analyses without requiring programming skills or their own computational resources (these are provided by KBase). Alongside these resources, KBase continues to work with instructors, supporting the development of additional curriculum modules. For anyone new to the platform, KBase, and the growing KBase Educators Organization, provides a community network, accompanied by community-sourced guidelines, instructional templates, and peer support to use KBase within a classroom whether virtual or in-person.

59 BASIC BIOLOGICAL SCIENCES↗

Integrated Technology for Cost-Effective CO2 Capture and Formic Acid Production: Modeling, Optimization, and Economic Analysis

A novel reactive technology is being investigated that electrochemically converts CO2 into valuable chemicals, particularly formic acid. This work focuses on identifying the optimal design and operation of an integrated membrane-based CO2 capture unit with the electrochemical conversion process. In this setup, the CO2 in the flue gas permeates through a CO2-selective membrane and enters an electrolyzer to produce formic acid, creating an integrated reaction module. To refine the chemical product, gas products from the electrolyzer are directed to a pressure swing adsorption unit, while the liquid product undergoes refinement to achieve commercial-grade formic acid using reactive distillation. A membrane CO2 capture model and an electrochemical conversion model have been developed using the IDAES Integrated Platform (Institute for the Design of Advanced Energy System), facilitating rigorous flowsheet modeling and process design and optimization.

Wang, Maojian↗

Geothermal Power Systems Analysis: Outcome of Industry Stakeholders Workshop

Geothermal cost and performance evaluation implemented via techno-economic assessment (TEA) modeling is critical for the U.S. Department of Energy (DOE) and other geothermal industry stakeholders in assessing the current state of geothermal technologies and to identify existing hurdles to commercially viable geothermal development. The Geothermal Electricity Technology Evaluation Model (GETEM) is a major TEA tool used in estimating the economic feasibility and levelized cost of energy (LCOE) of conventional hydrothermal systems and enhanced geothermal systems (EGS). Since 2021, GETEM has been transitioning from an intricate spreadsheet model to a user-friendly tool within the System Advisor Model (SAM) developed by the National Renewable Energy Laboratory (NREL). Apart from enabling an expanded visibility of the geothermal model among other renewable resources, having GETEM in SAM has the advantage of simulation automation, better usability, updates tracking, active user inputs/feedback, and extended financial modeling. GETEM is used in developing supply curves for NREL's Annual Technology Baseline (ATB), which provides inputs to the Renewable Energy Potential (reV) and the Regional Energy Deployment System (ReEDS) models. The geothermal module in NREL's reV model assesses the geothermal energy potential in the conterminous United States by defining the geospatial intersection of geothermal resources with existing grid infrastructure within the constraint of land use characteristics. The ReEDS model is a capacity expansion model used for simulating the long-term build-out and operation of the U.S. generation and transmission system based on current energy costs and policies. To ensure enhanced representation of current industry trends in our model transitions and development, we organized a two-day virtual workshop to elicit geothermal industry stakeholder input and recommendations on our current approaches and assumptions on techno-economic, resource assessment, and deployment scenarios modeling of geothermal technologies. Participants included developers, operators, investors, regulatory agencies, system modelers, national laboratory researchers, consultants, and other stakeholders. In this workshop, we gained stakeholder insights on current geothermal plant performance (i.e., capacity factors), updated drilling costs and learning curves, and next-generation technologies such as closed-loop and superhot rock geothermal. Other outcomes from this workshop and its impact on future geothermal development feasibility, resource availability, and capacity expansion studies are compiled and discussed.

annual technology baseline↗

Up‐And‐Coming Advances in Optical and Microwave Nonreciprocity: From Classical to Quantum Realm

Reciprocity is a fundamental physical principle that roots in the time‐reversal symmetry of physical laws. It allows making predictions on any arbitrary complex system's response and operation and hence simplifies the analysis. However, there are many practical situations in which it is advantageous to break reciprocity, e.g., isolators preventing wave scattering back to lasers and generators, full‐duplex systems for multiplexing transmission and receiving in the same channel, nonreciprocal cavity excitation, and protection of fragile states of superconductor quantum computers from thermal noise. The most widespread approach to time‐reversal symmetry breaking and nonreciprocity based on magnetic field biasing suffers from bulkiness, cost ineffectiveness, and loss, motivating researchers and engineers to search for more practical approaches. Herein, the up‐and‐coming advances in optical nonreciprocity, including new materials (Weyl semimetals, topological insulators, metasurfaces), active structures, time‐modulation, parity‐time (PT)‐symmetry breaking, nonlinearity combined with a structural asymmetry, quantum nonlinearity, unidirectional gain and loss, chiral quantum states and valley polarization are overviewed. A general description of nonreciprocal systems is provided and the pros and cons of the mentioned approaches to nonreciprocity are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Surface triggered stabilization of metastable charge-ordered phase in SrTiO 3

Charge ordering (CO), characterized by a periodic modulation of electron density and lattice distortion, has been a fundamental topic in condensed matter physics, serving as a potential platform for inducing novel functional properties. The charge-ordered phase is known to occur in a doped system with high d-electron occupancy, rather than low occupancy. Here, we report the realization of the charge-ordered phase in electron-doped (100) SrTiO 3 epitaxial thin films that have the lowest d-electron occupancy i.e., d 1 -d 0 . Theoretical calculation predicts the presence of a metastable CO state in the bulk state of electron-doped SrTiO 3 . Atomic scale analysis reveals that (100) surface distortion favors electron-lattice coupling for the charge-ordered state, and triggering the stabilization of the CO phase from a correlated metal state. This stabilization extends up to six unit cells from the top surface to the interior. Our approach offers an insight into the means of stabilizing a new phase of matter, extending CO phase to the lowest electron occupancy and encompassing a wide range of 3d transition metal oxides.

36 MATERIALS SCIENCE↗

The AGCVIII kinase Dw2 modulates cell proliferation, endomembrane trafficking, and MLG/xylan cell wall localization in elongating stem internodes of Sorghum bicolor

Stems of bioenergy sorghum (Sorghum bicolor L. Moench.), a drought–tolerant C4 grass, contain up to 50 nodes and internodes of varying length that span 4–5 m and account for approximately 84% of harvested biomass. Stem internode growth impacts plant height and biomass accumulation and is regulated by brassinosteroid signaling, auxin transport, and gibberellin biosynthesis. In addition, an AGCVIII kinase (Dw2) regulates sorghum stem internode growth, but the underlying mechanism and signaling network are unknown. Here we provide evidence that mutation of Dw2 reduces cell proliferation in internode intercalary meristems, inhibits endocytosis, and alters the distribution of heteroxylan and mixed linkage glucan in cell walls. Phosphoproteomic analysis showed that Dw2 signaling influences the phosphorylation of proteins involved in lipid signaling (PLDδ), endomembrane trafficking, hormone, light, and receptor signaling, and photosynthesis. Altogether, our results show that Dw2 modulates endomembrane function and cell division during sorghum internode growth, providing insight into the regulation of monocot stem development.

59 BASIC BIOLOGICAL SCIENCES↗

Clear Sky Toolkit: Decision Support User Guide

The Clear Sky Decision Support Template is a Microsoft Excel-based resource that enables users to screen and prioritize critical facilities for the addition of solar + storage. The Decision Support Template consolidates and streamlines the most relevant resources from the emergency management discipline, the electric utility sector, and national energy research institutions to assist in further site analysis and decision-making. The Decision Support Template is intended for use by local government and private sector decision-makers who aim to mitigate the effects of power outages on critical infrastructure, understand the potential for deploying photovoltaic (PV) and storage solutions on critical facilities, and align emergency management planning processes and utility needs with community resilience objectives. The Decision Support Template is divided into five modules that provide a structured, step-by-step approach users can follow and apply to unique circumstances. This process guides a more thorough and accurate assessment of a site's potential to accommodate a solar + storage system, provide grid benefits to utilities, and support community resilience and disaster mitigation for the broader region. After identifying potentially viable and preferred sites for solar + storage using the Decision Support Template, users should conduct a more detailed evaluation and consultation with a solar professional.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

High-throughput investigation of structural evolution upon solid-state in Cu–Cr–Co combinatorial multilayer thin-film

Cu-Cr-Co combinatorial multilayer thin-films were prepared by a high-throughput ion beam sputtering system. Based on the thickness ratio among the individual nanoscale monolayers (Cu, Cr, Co), the resulting stoichiometry covered the entire phase diagram. The chemical composition and structure of Cu-Cr-Co combinatorial chip upon solid-state reaction were studied by lab-based micro-X-ray fluorescence (mu-XRF) and high-throughput synchrotron X-ray diffraction (XRD), respectively. A composition-structure map for Cu-Cr-Co combinatorial chip was developed through automated data analysis employing hierarchical clustering techniques. The structural evolution of Cu-Cr-Co combinatorial chip as a function of heat-treatment temperature, time, and modulation period were studied systematically. Furthermore, the effect of the elemental distribution in-depth direction was investigated to gain more insights regarding phase transformation. This work provides an efficient method and new perspectives for the design and optimization of the composition and structure of high-performance thin-films.

36 MATERIALS SCIENCE↗

Modeling Tool Development and Validation for Solar Industry Process Heat Using Particle Thermal Energy Storage

U.S. industry sectors used 26.2 quadrillion Btu and accounted for 33% of total energy consumption in 2021 according to the Energy Information Agency. Industrial process heat accounts for 70% of industrial energy use with application temperatures ranging from 60 degrees -1100 degrees C. Industry processes, heavily relying on fossil fuels of cheap coal or natural gas, differ widely in operating conditions and load requirements which makes them difficult to standardize and imposes great challenges in decarbonization. Industry processes require reliable energy supply and vary widely in temperature ranges. Storing energy from renewable sources is necessary to improve reliability and to mitigate renewable intermittency when replacing carbon fuel-based heat supplies to achieve energy savings and reduce emissions. To this end, we have developed a particle-based thermal energy storage (TES) technology using low-cost and highly stable silica sand as a storage medium. The economic and performance-based analysis is key for renewable energy sources to reliably supply industry process heat and ultimately displace fossil fuels for decarbonization. The diversified industrial processes need case-by-case analysis and design. Therefore, an adaptive modeling tool is key for renewable power with energy storage to meet industry demands. Thus, a modeling tool to simulate a solar industry process heat system using the particle TES has been developed using the object-oriented equation-based language Modelica and the commercial platform of Modelon Impact. The Modelica-based software tool provides a general simulation environment for the design of reliable solar energy sources integrated with TES for various industrial process applications at different temperatures for economic competence with fossil fuels such as coal and natural gases. It uses both customized and standard component modeling modules in Modelon libraries for the flexibility to be adapted to a specific energy demand application. The particle TES system establishes a uniform energy supply platform with an efficient heat exchanger and particle thermal energy reservoir integrated with renewable powers. The particle TES system can provide a wide temperature range and can have a large storage temperature difference that increases storage energy density; therefore, it can be an adaptable energy storage system integrated with renewable power to supply 24/7 heat for industry decarbonization.

concentrated solar thermal↗

Optimization of Membrane-based Carbon Capture using Dimensional Analysis, CFD and Process System Engineering

Carbon capture is a promising option to mitigate CO2 emissions from existing coal-fired power plants, cement and steel industries, and petrochemical complexes. Among the available technologies, membrane-based carbon capture presents the lowest energy consumption, operating costs, and carbon footprint. In addition, membrane processes have important operational flexibil-ity and response times. On the other hand, the major challenges to widespread application of this technology are related to reducing capital costs and improving membrane stability and durability.To upscale the technology into stacked flat sheet configurations, high fidelity computational fluid dynamics (CFD) that describes the separation process accurately are required. High fidelity simulations have been shown to be effective in studying the complex transport phenomena in membrane systems. In addition, obtaining high CO2 recovery percentages and product purity requires a multi-stage membrane process, where the optimal network configuration of the membrane modules must be studied in a systematic way. In order to address the design problem at process scale, we formulate a superstructure for the membrane-based carbon capture, including up to three separation stages. In the formulation of the optimization problem, we include reduced models, based on rigorous CFD simulations of the membrane modules. Numerical results indicate that the optimal design includes three membrane stages, and the capture cost is 45.4 $/t-CO2.

Pedrozo, Hector A.↗

Microphysics of liquid water in sub-10 nm ultrafine aerosol particles

Abstract. Ultrafine aerosol particles with sizes smaller than 50 nm have been shown in recent studies to serve as a large source of cloud condensation nuclei that can promote additional cloud droplet formation under supersaturation conditions. Knowledge of the microphysics of liquid water in these droplets remains limited, particularly in the sub-10 nm particle size range, due to experimental and theoretical challenges associated with the complexity of aerosol components and the small length scales of interest (e.g., difficulty of precisely sampling the liquid–air interface, questionable validity of mean-field theoretical representations). Here, we carried out molecular dynamics (MD) simulations of aerosol particles with diameters between 1 and 10 nm and characterized atomistic-level structure and water dynamics in well-mixed and phase-separated systems with different particle sizes, NaCl salinities, and pimelic acid (PML) organic surface loadings as a function of distance from the time-averaged Gibbs dividing interface or instantaneous water–air interface. We define a sphericity factor (ϕ) that can shed light on the phase-mixing state of nanodroplets, and we reveal an unexpected dependence of mixing state on droplet size. Our results also evidence an ion concentration enhancement in ultrafine aerosols, which should modulate salt nucleation kinetics in sub-10 nm droplets, and provide detailed characterization of the influence of droplet size on surface tension and on water self-diffusivity near the interface. Analysis of water evaporation free energy and water activity demonstrates the validity of the Kelvin equation and Köhler theory at droplet sizes larger than 4 nm under moderate salinities and organic loadings and the need for further extension to account for ion concentration enhancement in sub-10 nm aerosols, droplet-size-dependent phase separation effects, and a sharp decrease in the cohesiveness of liquid water in sub-4 nm droplets. Finally, we show that an idealized fractional surface coating factor (fs) can be used to categorize and reconcile water accommodation coefficients (α*) observed in MD simulations and experimental results in the presence of organic coatings, and we resolve the droplet size dependence of α*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Protein3D: Enabling analysis and extraction of metal‐containing sites from the Protein Data Bank with molSimplify

Abstract Metalloenzymes catalyze a wide range of chemical transformations, with the active site residues playing a key role in modulating chemical reactivity and selectivity. Unlike smaller synthetic catalysts, a metalloenzyme active site is embedded in a larger protein, which makes interrogation of electronic properties and geometric features with quantum mechanical calculations challenging. Here we implement the ability to fetch crystallographic structures from the Protein Data Bank and analyze the metal binding sites in the program molSimplify. We show the usefulness of the newly created protein3D class to extract the local environment around non‐heme iron enzymes containing a two histidine motif and prepare 372 structures for quantum mechanical calculations. Our implementation of protein3D serves to expand the range of systems molSimplify can be used to analyze and will enable high‐throughput study of metal‐containing active sites in proteins.

Edholm, Freya↗

Online and Offline Data Quality Monitoring for the Mu2e Calorimeter

This thesis presents the design, implementation, and validation of a calorimeter Data Quality Monitoring (DQM) toolchain for the Mu2e experiment at Fermilab. Mu2e searches for charged lepton flavor violation via coherent muon-to-electron conversion in the field of an aluminum nucleus, $\mu^- Al \rightarrow e^-Al$, a process whose observation would constitute clear evidence of physics beyond the Standard Model. Achieving target sensitivity requires stringent control of detector performance and data integrity during acquisition, as subtle issues in readout configuration, data formatting, or electronics behavior can compromise reconstruction and bias downstream analyzes. To address these challenges, this work develops a multi-layer DQM approach spanning both raw data validation and reconstructed digi-level diagnostics. At the low level, a fragment analysis component performs word- and bit-field decoding of calorimeter readout blocks, enabling sanity checks of the expected structure and producing detailed error and integrity statistics useful for commissioning and troubleshooting. At the digi level, the CaloDigiDQM analyzer is implemented within the art framework and transforms each CaloDigiCollection into a structured hierarchy of ROOT histograms designed for fast drill-down diagnostics. The module generates coherent monitoring views at global, disk, board, and channel granularity, including occupancy, waveform-derived features (baseline, RMS, peak amplitude and position), and left-right sensor consistency metrics. Detector-aware channel-to-electronics mapping is performed through the conditions system (CaloDAQMap), ensuring that diagnostics remain aligned with hardware identifiers used in operations. For end-to-end testing without reliance on live DAQ data, a synthetic CaloDigi producer is developed to generate realistic waveforms with controlled noise and pulse shapes. The resulting system supports both offline ROOT-file production and online operation, including optional histogram streaming through otsdaq via ots::HistoSender. This toolchain provides a practical and scalable foundation for calorimeter commissioning and stable data collection, enabling early detection of anomalies and reducing operational risk for Mu2e.

Vakulenko, Mark [Drew U.] (ORCID:0009000276197818)↗

Finite Element Analysis (FEA) for Water-Foam Fracturing of Granite Rock

In addition to the foam data that were obtained from literature and that were collected from the current study, simulation data was also generated from finite element analysis (FEA) conducted in this study using COMSOL Multiphysics software. The FEA models were built to simulate the experiments conducted at Oak Ridge National Laboratory (ORNL) on cement and granite samples. In these FEA models, temperature was kept at ambient while the pressure profile resembled the loading conditions during the ORNL experiments, where pressure was either monotonically increased or applied cyclically. The cement material was used as a model material and was used to study Von Mises stress and tensile stress distribution for different bore hole length geometry using a parametric sweep with water as fracturing fluid using solid-fluid interaction module. For the granite material, FEA models were developed for stress analysis of cylindrical samples with water or foam fluids. The solid mechanics module in COMSOL was implemented to solve for Von Mises stress and tensile stress. The fluid-structure interaction module was implemented to solve for water-foam interaction on granite cylinder with addition of fluid-loading on structure, i.e., large deformation in solid mechanics with no impact on fluid deformation. Foam was considered as a pseudo single-phase compressible fluid for which material properties were calculated from water and gas (nitrogen) phases. The density of foam is calculated as a function of the densities of water and nitrogen, while viscosity is a function of temperature. Four types of FEA analyses were modelled: 1. Monotonic injection with water 2. Monotonic injection with foam 3. Cyclic injection with water 4. Cyclic injection with foam All the COMSOL files are converted to a zip file which is save in .mph.

15 GEOTHERMAL ENERGY↗