Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fractional operators”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Identification of Safety Metrics for Airport Surface Operations

A large fraction of safety incidents occurs on the ground during airport surface operations. Although these incidents are mostly non-fatal with a few exceptions, they are high profile incidents that remain a source of concern for the National Transportation Safety Board (NTSB), the Federal Aviation Administration (FAA), major airlines, and other stakeholders of the National Airspace System (NAS). These incidents have historically been mitigated by implementing changes to regulations, policies, and procedures over time. This approach has minimized but not eliminated the risk of occurrence of safety incidents. It is thus important to develop integrated techniques to assess, model, and prevent these incidents by analyzing the risk and likelihood of occurrence and communicating results of the analysis to decision-making personnel who can mitigate and prevent incidents in real time. The work presented in this paper builds on a previously developed architecture for safety, Real-Time Safety Monitoring (RTSM), to enable monitoring and prediction of the safety of the NAS. In the RTSM framework, hazards to flight are translated to safety metrics such as wake vortex encounters or loss of separation, that can be modeled and analyzed offline and also predicted and monitored in real time (online). The intent of this paper is to integrate predictable incidents that occur during surface and ground operations into the safety portfolio of the RTSM project by (i) identifying suitable information sources from which ground incidents can be studied, (ii) developing safety metrics correlated with surface operations, and (iii) recommending suitable data sources that can be quantified and used for the computation of pertinent safety metrics.

safety↗
DOE OSTI2021

Identifiability and predictability of integer- and fractional-order epidemiological models using physics-informed neural networks

Here we analyze a plurality of epidemiological models through the lens of physics-informed neural networks (PINNs) that enable us to identify time-dependent parameters and data-driven fractional differential operators. In particular, we consider several variations of the classical susceptible-infectious-removed (SIR) model by introducing more compartments and fractional-order and time-delay models. We report the results for the spread of COVID-19 in New York City, Rhode Island and Michigan states and Italy, by simultaneously inferring the unknown parameters and the unobserved dynamics. For integer-order and time-delay models, we fit the available data by identifying time-dependent parameters, which are represented by neural networks. In contrast, for fractional differential models, we fit the data by determining different time-dependent derivative orders for each compartment, which we represent by neural networks. We investigate the structural and practical identifiability of these unknown functions for different datasets, and quantify the uncertainty associated with neural networks and with control measures in forecasting the pandemic.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1986

Finite element analysis of steady and transiently moving/rolling nonlinear viscoelastic structure. Part 1: Theory

In a three part series of papers, a generalized finite element analysis scheme is developed to handle the steady and transient response of moving/rolling nonlinear viscoelastic structure. This paper considers the development of the moving/rolling element strategy, including the effects of large deformation kinematics and viscoelasticity modelled by fractional integro-differential operators. To improve the solution strategy, a special hierarchical constraint procedure is developed for the case of steady rolling/translating as well as a transient scheme involving the use of a Grunwaldian representation of the fractional operator. In the second and third parts of the paper, 3-D extensions are developed along with transient contact strategies enabling the handling of impacts with obstructions. Overall, the various developments are benchmarked via comprehensive 2- and 3-D simulations. These are correlated with experimental data to define modelling capabilities.

Padovan, Joe↗
DOE OSTI2022

Learning functional priors and posteriors from data and physics

In this work, we develop a new Bayesian framework based on deep neural networks to be able to extrapolate in space-time using historical data and to quantify uncertainties arising from both noisy and gappy data in physical problems. Specifically, the proposed approach has two stages: (1) prior learning and (2) posterior estimation. At the first stage, we employ the physics-informed Generative Adversarial Networks (PI-GAN) to learn a functional prior either from a prescribed function distribution, e.g., Gaussian process, or from historical data and physics. At the second stage, we employ the Hamiltonian Monte Carlo (HMC) method to estimate the posterior in the latent space of PI-GANs. In addition, we use two different approaches to encode the physics: (1) automatic differentiation, used in the physicsinformed neural networks (PINNs) for scenarios with explicitly known partial differential equations (PDEs), and (2) operator regression using the deep operator network (DeepONet) for PDE-agnostic scenarios. We then test the proposed method for (1) meta-learning for one-dimensional regression, and forward/inverse PDE problems (combined with PINNs); (2) PDE-agnostic physical problems (combined with DeepONet), e.g., fractional diffusion as well as saturated stochastic (100-dimensional) flows in heterogeneous porous media; and (3) spatial-temporal regression problems, i.e., inference of a marine riser displacement field using experimental data from the Norwegian Deepwater Programme (NDP). The results demonstrate that the proposed approach can provide accurate predictions as well as uncertainty quantification given very limited scattered and noisy data, since historical data could be available to provide informative priors. In summary, the proposed method is capable of learning flexible functional priors, e.g., both Gaussian and non-Gaussian process, and can be readily extended to big data problems by enabling mini-batch training using stochastic HMC or normalizing flows since the latent space is generally characterized as low dimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

Two hundred passage three-way valve: Fraction collector

This paper describes the design and operation of a fraction collector used to direct flow of separated biological materials from 197 capillary tubes to either a collection tray or to a waste tank. This mechanism uses a 28-volt dc gear motor driving twin cams to force 197 needles through a self-sealing silicone rubber septum, where they inject the material in 197 separate pockets in a collection tray. The position of the collector tray is sensed by two optical limit switches. The time sequences are controlled automatically by an electronics control monitoring module.

Keffer, J. L.↗
DOE OSTI2024

On Road vs. Off Road Low Load Cycle Comparison

Reducing criteria pollutants while reducing greenhouse gases is an active area of research for commercial on-road vehicles as well as for off-road machines. The heavy duty on-road sector has moved to reducing NOx by 82.5% compared to 2010 regulations while increasing the engine useful life from 435,000 to 650,000 miles by 2027 in the United States (US). An additional certification cycle, the Low Load Cycle (LLC), has been added focusing on part load operation having tight NOx emissions levels. In addition to NOx, the total CO2 emissions from the vehicle will also be reduced for various model years. The off-road market is following with a 90% NOx reduction target compared to Tier 4 Final for 130-560 kW engines along with greenhouse gas targets that are still being established. The off-road market will also need to certify with a Low Load Application Cycle (LLAC), a version of which was proposed for evaluation in 2021. Since the LLAC has not been finalized, this study is being conducted to compare and contrast the LLC for on-road with the LLAC for off-road as there might be some shared learnings. A US off-road production 2023 Fiat Powertrain 13L engine and aftertreatment system was chosen for this work. This engine is used in production for both off-road and on-road products, so it is a good choice for this study. The associated off-road aftertreatment system was aged for more relevant comparisons. The engine calibration was not altered for either of the low load cycles. This study shows that the cycles are quite different in nature as the market needs are different. The LLC includes a large fraction of operation at idle and lower speeds, representing products that use the engine primarily for motive power, where lower vehicle speed means a lower engine speed and load. The LLAC has more time and load spent at high speeds and slightly higher loads. The off-road products represented by this cycle often use the engine to drive auxiliary equipment which means higher parasitic loads and hand/fixed throttle. The comparison will include the use profiles, tailpipe NOx and greenhouse gas emissions (CO2, N2O).

McCarthy, James↗
DOE OSTI2024

Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence

We consider a nearly collisionless plasma consisting of a species of “test particles” in one spatial and one velocity dimension, stirred by an externally imposed stochastic electric field—a kinetic analog of the Kraichnan model of passive advection. The mean effect on the particle distribution function is turbulent diffusion in velocity space—known as stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant C 2 ∝ ∫ ∫ d x d v 〈 f 2 〉 —a quantity that here plays the role of (negative) entropy of the distribution function. We find that C 2 is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of C 2 in Fourier space, with k and s denoting spatial and velocity wave numbers, respectively. In our model, the nonlinearity in the evolution equation for the spectrum turns into a fractional Laplacian operator in k space, leading to anomalous diffusion. Whereas even the linear phase mixing alone would lead to a constant flux of C 2 to high s (towards the collisional dissipation range) at every k , the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of C 2 is to both high k and high s simultaneously. Integrating over velocity (spatial) wave numbers, the k -space ( s -space) flux of C 2 is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the ( k , s ) plane, with power-law asymptotics at large k and s . Our model is fully analytically solvable, but the asymptotic scalings of the spectrum can also be found via a simple phenomenological theory whose key assumption is that the cascade is governed by a “critical balance” in phase space between the linear and nonlinear timescales. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Decoupled Fréchet kernels based on a fractional viscoacoustic wave equation

We formulate the Fréchet kernel computation using the adjoint-state method based on a fractional viscoacoustic wave equation. We first numerically prove that both the 1/2- and the 3/2-order fractional Laplacian operators are self-adjoint. Using this property, we show that the adjoint wave propagator preserves the dispersion and compensates the amplitude, while the time-reversed adjoint wave propagator behaves identically as the forward propagator with the same dispersion and dissipation characters. Without introducing rheological mechanisms, this formulation adopts an explicit Q parameterization, which avoids the implicit Q in the conventional viscoacoustic/viscoelastic full waveform inversion (Q-FWI). In addition, because of the decoupling of operators in the wave equation, the viscoacoustic Fréchet kernel is separated into three distinct contributions with clear physical meanings: lossless propagation, dispersion, and dissipation. Here, we find that the lossless propagation kernel dominates the velocity kernel, while the dissipation kernel dominates the attenuation kernel over the dispersion kernel.

58 GEOSCIENCES↗
NASA NTRS1988

The geology and petrogenesis of the southern closepet granite

The Archaean Closepet Granite is a polyphase body intruding the Peninsular Gneiss Complex and the associated supracrustal rocks. The granite out-crop runs for nearly 500 km with an approximate width of 20 to 25 km and cut across the regional metamorphic structure passing from granulite facies in the South and green schist facies in the north. In the amphibolite-granulite facies transition zone the granite is intimately mixed with migmatites and charnockite. Field observations suggests that anatexis of Peninsular gneisses led to the formation of granite melt, and there is a space relationship between migmatite formation, charnockite development and production and emplacement of granite magma. Based on texture and cross cutting relationships four major granite phases are recognized: (1) Pyroxene bearing dark grey granite; (2) Porphyritec granite; (3) Equigranular grey granite; and (4) Equigranular pink granite. The granite is medium to coarse grained and exhibit hypidiomorphic granular to porphyritic texture. The modal composition varies from granite granodiorite to quartz monzonite. Geochemical variation of the granite suite is consistent with either fractional crystallization or partial melting, but in both the cases biotite plus feldspar must be involved as fractionating or residual phases during melting to account trace element chemistry. The trace element data has been plotted on discriminant diagrams, where majority of samples plot in volcanic arc and within plate, tectonic environments. The granite show distinct REE patterns with variable total REE content. The REE patterns and overall abundances suggests that the granite suite represents a product of partial melting of crustal source in which fractional crystallization operated in a limited number of cases.

Jayananda, M.↗
NASA NTRS1998

Insights Into the Fractional Order Initial Value Problem via Semi-Infinite Systems

This paper considers various aspects of the initial value problem for fractional order differential equations. The main contribution of this paper is to use the solutions to known spatially distributed systems to demonstrate that fractional differintegral operators require an initial condition term that is time-varying due to past distributed storage of information.

Hartley, Tom T.↗
NASA NTRS2002

Experiments and Spectral Studies of Martian Volcanic Rocks: Implications for the Origin of Pathfinder Rocks and Soils

The composition and spectral properties of the Mars Pathfinder rocks and soils together with the identification of basaltic and andesitic Mars terrains based on Thermal Emission Spectrometer (TES) data raised interesting questions regarding the nature and origin of Mars surface rocks. We have investigated the following questions: (1) are the Pathfinder rocks igneous and is it possible these rocks could have formed by known igneous processes, such as equilibrium or fractional crystallization, operating within SNC magmas known to exist on Mars? If it is possible, what P (depth) and PH2O conditions are required? (2) whether TES-based interpretations of plagioclase-rich basalt and andesitic terrains in the south and north regions of Mars respectively are unique. Are the surface compositions of these regions plagioclase-rich, possibly indicating the presence of old AI-rich crust of Mars, or are the spectra being affected by something like surface weathering processes that might determine the spectral pyroxene to plagioclase ratio?

Rutherford, Malcolm J.↗
DOE OSTI2022

BETO 2021 Peer Review - FCIC Task 5 - Preprocessing

In the study of feedstock variability, this presentation specifically looks at preprocessing processes that will created well-defined and homogenous feedstock from variable biomass resources, in the quest to produce homogeneous intermediates that will be converted into market-ready products. The objective of this stage of the project is to develop science-based design and operation principles informed by TEA/LCA that result in predictable, reliable, and scalable performance of preprocessing unit operations (comminution, fractionation, deconstruction, and real-time imaging).

biomass↗
DOE OSTI2021

Connections between nonlocal operators: From vector calculus identities to a fractional Helmholtz decomposition

Nonlocal vector calculus, which is based on the nonlocal forms of gradient, divergence, and Laplace operators in multiple dimensions, has shown promising applications in fields such as hydrology, mechanics, and image processing. In this work, we study the analytical underpinnings of these operators. We rigorously treat compositions of nonlocal operators, prove nonlocal vector calculus identities, and connect weighted and unweighted variational frameworks. We combine these results to obtain a weighted fractional Helmholtz decomposition which is valid for sufficiently smooth vector fields. Our approach identifies the function spaces in which the stated identities and decompositions hold, providing a rigorous foundation to the nonlocal vector calculus identities that can serve as tools for nonlocal modeling in higher dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Self-consistent investigation of density fueling needs on ITER and CFETR utilizing the new Pellet Ablation Module

Abstract Self-consistent modeling using the stability, transport, equilibrium, and pedestal (STEP) workflow in the OMFIT integrated modeling framework (predicting pedestal with EPED, core profiles with TGYRO, current profile with ONETWO, and EFIT for equilibrium) suggests ITER and future devices such as China Fusion Engineering Test Reactor (CFETR) Zhuang et al (2019 Nucl. Fusion 59 112010) will benefit from high-density operation (Greenwald limit fraction f g w ≈ 0.7−1.3). Regimes with an operational density near the Greenwald limit will likely need peaked density profiles so that the pedestal density remains below the Greenwald limit. Peaked density profiles can be achieved with the help of pellet injection. A flexible Pellet Ablation Module (PAM), which predicts the density source based on a comprehensive analytical pellet ablation model, has been developed for predicting pellet fueling for transport studies, and has been incorporated into the STEP workflow for predictive modeling. This workflow is applied to DIII-D and finds good agreement with experiments. On ITER the effect of pellet fueling is examined in an advanced inductive scenario, where a fusion gain of up to Q = 9 is predicted with strong central pellet fueling. On CFETR, with a mid-radius density source, an average of 1.5 × 10 22 electrons s −1 are required to achieve the density and temperature profiles necessary for the 1000 MW advanced scenario with a tritium burn-up fraction of ∼ 3 % .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

CERF: IM3 Projected Western US Power Plant Locations

Overview The Capacity Expansion Regional Feasibility (CERF) model is an open-source geospatial python package that provides new power plant locations at a 1km resolution. The model ingests U.S. state or regional-scale electricity system capacity expansion plans, such as those produced by the Global Change Analysis Model (GCAM-USA), and identifies feasible, site-specific locations for individual new power plants (renewable and non-renewable). CERF combines high-resolution geospatial suitability analyses with an economic algorithm that selects individual plant siting locations based on grid interconnection costs and the locational marginal value of new generation. The model incorporates a wide range of dynamic constraints and opportunities, such as protected lands, population density, existing infrastructure, and water availability. This dataset provides CERF power plant siting results for IM3 Phase 2 simulations across eight different scenarios for the Western US through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 CERF siting results in this dataset correspond to capacity expansion plans in the GCAM-USA IM3 Phase 2 simulation data and are available for each of the above scenarios. Data Details Temporal Range: 2015-2055 in 5-year timesteps. Note that 2015 is the experiment base year and 2020 and beyond represent model simulation years. Spatial Range: Plant locations are provided for the eleven states in the Western US including Arizona, California, Colorado, Idaho, Montana, New Mexico, Nevada, Oregon, Utah, Washington, and Wyoming. Spatial Resolution: 1 km-squared, provided in x and y coordinates Geospatial Projection: Albers Equal Area Conic (ESRI:102003) File Type: csv The dataset contains subdirectories for each of the eight scenarios described in the overview. Each scenario folder contains two subfolders with the following information: 1. Power Plant Data This directory contains a single .csv file of power plant locations for both pre-existing (non-CERF sited plants in operation in 2015) and new (CERF-sited) power plants across the temporal range along with additional CERF model output parameters for CERF-sited plants. Plant with a siting year earlier than 2020 correspond to facilities that are operational leading into the first timestep CERF simulation. For a more detailed description of CERF model output parameters, see the CERF model documentation. Note that the cerf_plant_id parameter is unique within each scenario file but not across scenario files. Parameter Descriptions scenario - Name of scenario cerf_plant_id - Unique siting identifier cerf_sited - If True, indicates that plant was sited by CERF model. If False, indicates pre-existing facility region_name - Name of region (state) tech_id - Technology ID tech_name - Full generation technology name inclusive of cooling type (if applicable) and additional characteristics tech_simple - Simplified generation technology type unit_size_mw - Power plant unit size (MW) xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) index - Index position in the flattend 2D array buffer_in_km - Exclusion buffer around site (km) sited_year - Year of siting retirement_year - Year of retirement lmp_zone - Locational marginal price (LMP) zone ID locational_marginal_price_usd_per_mwh - Locational marginal price ($/MWh) generation_mwh_per_year - Generation output (MWh/yr) operating_cost_usd_per_year - Cost of plant operations ($/yr) net_operational_value - Net operational value based on LMP and and operating costs ($/yr) interconnection_cost - Cost of interconnection for transmission & gas pipeline (if applicable) net_locational_cost -- Difference of interconnection cost and operating value ($/yr) capacity_factor_fraction - Capacity factor (fraction) carbon_capture_rate_fraction - Carbon capture rate (fraction) fuel_co2_content_tons_per_btu - Fuel CO2 content (tons/Btu) fuel_price_usd_per_mmbtu - Fuel price ($/MMBtu) fuel_price_esc_rate_fraction - Fuel price escalation rate (fraction) heat_rate_btu_per_kWh - Heat rate (Btu/kWh) lifetime_yrs - Technology lifetime for annuity (years) operational_life_yrs - Operational lifetime for retirement (years) variable_om_usd_per_mwh - Variable operation and maintenance costs of yearly capacity use ($/MWh) variable_om_esc_rate_fraction - Variable operation and maintenance costs escalation rate (fraction) carbon_tax_usd_per_ton - Carbon tax ($/ton) carbon_tax_esc_rate_fraction - Carbon tax escalation rate (fraction) 2. Storage Data This directory contains information on new and pre-existing energy storage facilities operational in each timestep along with various storage operational parameters. The 2015 timestep provides pre-existing energy storage data and corresponds with facilities that are operational leading into the first model simulation timestep. Note that coordinates in the storage files correspond to the interconnection point on the grid (substation location), not individual energy storage locations. Energy storage is added in a cumulative process at each given interconnection point. That is, each individual file provides the total operational storage capacity interconnected to the specified substation for the given timestep, inclusive of previously installed storage at that location and new storage installed in that timestep at that location. Parameters scenario - Name of scenario timestep - Simulation timestep name - Unique storage identifier s_typ - Type of energy storage technology (battery or pumped storage hydro) s_node - Node ID of interconnecting substation xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) charge_rate - Maximum charge rate (power capacity) of storage system (MW) discharge_rate - Maximum discharge rate (power capacity) of storage system (MW) duration - Duration of storage system (hours) max_SoC - Allowed maximum state of charge (energy capacity) of storage system (MWh) min_SoC -Allowed minimum state of charge (energy capacity) of storage system (MWh) charge_eff - Efficiency of charge (fraction between 0 and 1) discharge_eff - Efficiency of discharge (fraction between 0 and 1) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

CERF↗
NASA NTRS2014

Investigation of the Effects of Cathode Flow Fraction and Position on the Performance and Operation of the High Voltage Hall Accelerator

The National Aeronautics and Space Administration (NASA) Science Mission Directorate In- Space Propulsion Technology office is sponsoring NASA Glenn Research Center (GRC) to develop a 4 kW-class Hall thruster propulsion system for implementation in NASA science missions. Tests were performed within NASA GRC Vacuum Facility 5 at background pressure levels that were six times lower than what has previously been attained in other vacuum facilities. A study was conducted to assess the impact of varying the cathode-to-anode flow fraction and cathode position on the performance and operational characteristics of the High Voltage Hall Accelerator (HiVHAc) thruster. In addition, the impact of injecting additional xenon propellant in the vicinity of the cathode was also assessed. Cathode-to-anode flow fraction sensitivity tests were performed for power levels between 1.0 and 3.9 kW. It was found that varying the cathode flow fraction from 5 to approximately 10% of the anode flow resulted in the cathode-to-ground voltage becoming more positive. For an operating condition of 3.8 kW and 500 V, varying the cathode position from a distance of closest approach to 600 mm away did not result in any substantial variation in thrust but resulted in the cathode-to-ground changing from -17 to -4 V. The change in the cathode-to-ground voltage along with visual observations indicated a change in how the cathode plume was coupling to the thruster discharge. Finally, the injection of secondary xenon flow in the vicinity of the cathode had an impact similar to increasing the cathode-to-anode flow fraction, where the cathode-to-ground voltage became more positive and discharge current and thrust increased slightly. Future tests of the HiVHAc thruster are planned with a centrally mounted cathode in order to further assess the impact of cathode position on thruster performance.

Cathode Coupling↗
DOE OSTI2025

Gas-Gap Calorimeter Sizing and Rating Tools for Heat Pipe Experimentation

Effective gas-gap calorimeter sizing and rating tools are required to design calorimeters that enable well-defined operating conditions, high cooling powers, and accurate calorimetry. The present report describes two tools utilizing Microsoft Excel and MathWorks MATLAB that enable the sizing and rating of gas-gap calorimeters with a He-Ar binary gas mixture in the gap. The Excel tools are two easy-to-use spreadsheets for calculating heat pipe operating temperature or He-Ar mole fractions given the required inputs. The MATLAB tool includes similar models with more robust solution algorithms along with sub-options for full-vacuum, static gas-gap, and forced circulation in the gas-gap. The accuracy of the developed tools were demonstrated by comparing with existing work in literature. The MATLAB tool was used for a parametric study to determine the gas-gap thickness, coolant gap thickness, coolant flow rate, and coolant inlet temperatures for the testing of a 3/4 in outer diameter heat pipe up to powers of 10 kW and temperatures of 1,000°C. The effects of heat pipe and calorimeter inner shell surface emissivities were investigated, considering the inability to control or accurately measure surface emissivities in most applications. It was found that controlling the He-Ar mole fractions provides a wide operating range for gas-gap thicknesses around ~ 0.042--0.090 in (~ 1.07--2.29 mm) for a surface emissivity range of 0.4--0.8, since conduction heat transfer is a significant fraction of the heat transfer across the gas-gap. In addition, it was found that turbulent flow in the coolant gap is needed at high input powers to prevent shell temperatures from approaching the boiling temperature of the water coolant. The parametric study resulted in choosing stainless steel tubes with an outer diameter and thickness of 1 x 0.065 in as the inner shell, and a 1-1/2 x 0.156 in as the outer shell of the calorimeter. Overall, this report presents the necessary information for the sizing and rating of gas-gap calorimeters with He-Ar mixtures for the testing of high-temperature heat pipes (~ 500--1,000°C).

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Methane foam performance evaluation in water-wet unconsolidated porous media: A systematic experimental investigation at elevated pressure and temperature conditions

In the petroleum industry, enhanced oil recovery (EOR) techniques employ foam extensively to establish conformance control in heterogeneous and fractured reservoirs in order to increase the sweep efficiency. In such applications, foam performance evaluation under complex subsurface conditions is pivotal for the effective and optimized deployment of foam treatment. However, there is a scarcity of hydrocarbon gas foam generation and evaluation studies that examine the relationships between foam performance and the critical foam parameters at high-pressure and high-temperature conditions. Here, this study aims at methodically investigating the effects of several foam parameters on methane foam performance in water-wet proppant packs under harsh operating conditions. This is in relation to the technical needs of hydrocarbon foam injection into hydraulically-induced, propped fractures in unconventional oil reservoirs. To this end, a state-of-the-art experimental foam generation apparatus was designed, fabricated, and commissioned. We performed a large number of foam flow experiments on proppant packs using methane gas and different foaming agents at 3500 psi and 115°C. Anionic and amphoteric surfactants were employed to probe the effect of their ionic nature on foam performance. Foam performance sensitivities to various foam generation parameters and operating conditions, such as surfactant concentration, gas fraction, total injection rate, operating pressure, salinity, and proppant pack length were investigated. To this end, steady-state pressure drops across the proppant packs during foam generation and foam's apparent viscosity were measured to quantify the foam performance of surfactants. The results were then analyzed to determine optimum values of the foam parameters and the interplay between these parameters are discussed here. The systematic results achieved from this work are in agreement with the trends available in the literature and provide new insights into complexities of in situ foam generation in water-wet, unconsolidated porous media at extreme reservoir conditions.

03 NATURAL GAS↗