Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solid dynamics”

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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Analysis of AC-DC Converters for Grid-tied High Temperature Steam Electrolysis Systems

Grid-tied HTSE systems have the prospects to produce clean hydrogen enabling power and broad energy systems decarbonization. Most commercially available power electronic converter systems (PECS) are designed for batteries, solar PVs, wind, and other well-established renewable energy resources. Standards (such as IEEE 1547, UL 1741, CA Rule-21, HI Rule14) exist for PECS used for renewable energy systems such as battery storage and solar PVs. However, those that consider the dynamic behavior of HTSEs and that can be used for large-scale H2 systems are yet to be developed. This paper investigates the performance of these PECS for the HTSE application that are set up at the Idaho National Laboratory for hydrogen production testing, research and development. In particular, the performance analysis of two grid-tied PECS (A and B) is conducted for a 100 kW solid oxide HTSE system. System A consists of 6 units of 30kW MOSFET-switched bidirectional AC-DC rectifier while system B has a single unit of 150kW thyristor-switched AC-DC rectifier. Both systems are connected to the HTSE stacks via a DC-DC converter. Different operational conditions of the HTSE system are tested to analyze the dynamic response of the HTSE’s PECS. The experimental results show the need to develop advanced control strategies for PECS that incorporates the dynamics of HTSE systems for improved performance.

High temperature steam electrolysis↗
DOE OSTI2020

Understanding Superionic Conductivity in Lithium and Sodium Salts of Weakly Coordinating Closo -Hexahalocarbaborate Anions

Solid-state ion conductors based on closo-polyborate anions combine high ionic conductivity with a rich array of tunable properties. Cation mobility in these systems is intimately related to the strength of the interaction with the neighboring anionic network and the energy for reorganizing the coordination polyhedra. In this paper, we explore such factors in solid electrolytes with two anions of the weakest coordinating ability, [HCB 11 H 5 Cl 6 ] – and [HCB 11 H 5 Br 6 ] – , and a total of 11 polymorphs are identified for their lithium and sodium salts. Our approach combines ab initio molecular dynamics, synchrotron X-ray powder diffraction, differential scanning calorimetry, and AC impedance measurements to investigate their structures, phase-transition behavior, anion orientational mobilities, and ionic conductivities. We find that M(HCB 11 H 5 X 6 ) (M = Li, Na, X = Cl, Br) compounds exhibit order–disorder polymorphic transitions between 203 and 305 °C and display Li and Na superionic conductivity in the disordered state. Through detailed analysis, we illustrate how cation disordering in these compounds originates from a competitive interplay among the lattice symmetry, the anion reorientational mobility, the geometric and electronic asymmetry of the anion, and the polarizability of the halogen atoms. These factors are compared to other closo-polyborate-based ion conductors to suggest guidelines for optimizing the cation–anion interaction for fast ion mobility. This study expands the known solid-state poly(carba)borate-based materials capable of liquid-like ionic conductivities, unravels the mechanisms responsible for fast ion transport, and provides insights into the development of practical superionic solid electrolytes.

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

A Transformational Natural Gas Fueled Dynamic SOFC for Data Center in-Rack Power (Final Technical Report)

Data centers consume 250TWh or 1% of global electricity use in 2019. This electricity demand is expected to grow exponentially in the future. The current data center power supply relies on electric grids with diesel engines and batteries as backup; the latter is critically important because a sudden loss of power could erase crucial data, causing catastrophic economic loss. However, such a double protection architecture is costly, inefficient and polluting. In search for better in-rack power system, fuel-cell based generators have recently emerged as a promising alternative to diesel engines and batteries because of its potential to be more cost-effective, more efficient and less polluting. Among many types of fuel cells, solid oxide fuel cells stand out to be a front runner because of their ability to directly operate on natural gas that has an existing infrastructure. However, there are several critical challenges facing SOFC generators to meet the rigorous power requirements of data centers: high cost, poor reliability and more importantly, slow response to the datacenters’ dynamic load change. The unexpected overload conditions in data centers can cause fuel starvation; the latter can initiate microcracks in the anode structure, thus jeopardizing the lifetime of the SOFC generator. The current efforts to increase transient power capabilities of traditional SOFCs and therefore protect SOFC stacks from internal damage during transient loading, rely primarily upon mathematical algorithms. One of the challenges to these control methodologies is the physical limitations of the fuel delivery system such as slow response of mass flow controllers to fast power demand.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS2014

Spectroscopic Measurement Techniques for Aerospace Flows

The conditions that characterize aerospace flows are so varied, that a single diagnostic technique is not sufficient for its measurement. Fluid dynamists use knowledge of similarity to help categorize and focus on different flow conditions. For example, the Reynolds number represents the ratio of inertial to viscous forces in a flow. When the velocity scales, length scales, and gas density are large and the magnitude of the molecular viscosity is low, the Reynolds number becomes large. This corresponds to large scale vehicles (e.g Airbus A380), fast moving objects (e.g. artillery projectiles), vehicles in dense fluids (e.g. submarine in water), or flows with low dynamic viscosity (e.g. skydiver in air). In each of these cases, the inertial forces dominate viscous forces, and unsteady turbulent fluctuations in the flow variables are observed. In contrast, flows with small length scales (e.g. dispersion of micro-particles in a solid rocket nozzle), slow moving objects (e.g. micro aerial vehicles), flows with low density gases (e.g. atmospheric re-entry), or fluids with a large magnitude of viscosity (e.g. engine coolant flow), all have low Reynolds numbers. In these cases, viscous forces become very important and often the flows can be steady and laminar. The Mach number, which is the ratio of the velocity to the speed of sound in the medium, also helps to differentiate types of flows. At very low Mach numbers, acoustic waves travel much faster than the object, and the flow can be assumed to be incompressible (e.g. Cessna 172 aircraft). As the object speed approaches the speed of sound, the gas density can become variable (e.g. flow over wing of Learjet 85). When the object speed is higher than the speed of sound (Ma > 1), the presences of shock waves and other gas dynamic features can become important to the vehicle performance (e.g. SR-71 Blackbird). In the hypersonic flow regime (Ma > 5), large changes in temperature begin to affect flow properties, causing real-gas effects to occur (e.g. X-43 Scramjet). At even higher Mach numbers, chemistry and nonequilibrium effects come into play (e.g. Startdust re-entry capsule), further complicating the measurement. These limits can be predicted by calculating the ratio of chemical and thermal relaxation time to the flow time scales. Other non-dimensional numbers can be used to further differentiate types of aerospace flows.

Danehy, Paul M.↗
NASA NTRS1996

CFD Analysis of the 24-inch JIRAD Hybrid Rocket Motor

A series of multispecies, multiphase computational fluid dynamics (CFD) analyses of the 24-inch diameter joint government industry industrial research and development (JIRAD) hybrid rocket motor is described. The 24-inch JIRAD hybrid motor operates by injection of liquid oxygen (LOX) into a vaporization plenum chamber upstream of ports in the hydroxyl-terminated polybutadiene (HTPB) solid fuel. The injector spray pattern had a strong influence on combustion stability of the JIRAD motor so a CFD study was initiated to define the injector end flow field under different oxidizer spray patterns and operating conditions. By using CFD to gain a clear picture of the flow field and temperature distribution within the JIRAD motor, it is hoped that the fundamental mechanisms of hybrid combustion instability may be identified and then suppressed by simple alterations to the oxidizer injection parameters such as injection angle and velocity. The simulations in this study were carried out using the General Algorithm for Analysis of Combustion SYstems (GALACSY) multiphase combustion codes. GALACSY consists of a comprehensive set of droplet dynamic submodels (atomization, evaporation, etc.) and a computationally efficient hydrocarbon chemistry package built around a robust Navier-Stokes solver optimized for low Mach number flows. Lagrangian tracking of dispersed particles describes a closely coupled spray phase. The CFD cases described in this paper represent various levels of simplification of the problem. They include: (A) gaseous oxygen with combusting fuel vapor blowing off the walls at various oxidizer injection angles and velocities, (B) gaseous oxygen with combusting fuel vapor blowing off the walls, and (C) liquid oxygen with combusting fuel vapor blowing off the walls. The study used an axisymmetric model and the results indicate that the injector design significantly effects the flow field in the injector end of the motor. Markedly different recirculation patterns are observed in the vaporization chamber as the oxygen velocity and/or spray pattern is varied. The ability of these recirculation patterns to stabilize the diffusion flame above the surface of the solid fuel gives a plausible explanation for the experimentally determined combustion stability characteristics of the JIRAD motor, and suggests how combustion stability can be assured by modifications to the injector design.

Liang, Pak-Yan↗
DOE OSTI2025

Computational Fluid Dynamics Simulations of Glass Vitrification Refractory Coupon Tests

The Waste Treatment and Immobilization Plant (WTP) at the Hanford site is nearing the start of the Direct-Feed Low-Activity Waste (DFLAW) operations. DFLAW is destined to convert a pretreated low activity waste portion of the 56 million gallons of tank waste into a stable solid glass. In the subsequent decade completion of the high-level waste (HLW) facility is anticipated. Sustained operational missions of both LAW and HLW melter facilities are expected over multiple decades. In high-temperature glass melters, the refractory lining corrodes over time, which could potentially be an issue for longer term operations, this refractory corrosion is higher at the level of the glass-air interface due to surface tension driven flow. The glass viscosity, melt pool temperature, and glass chemical composition can impact the rate at which the refractory corrodes. This rate is important to quantify for the various waste glasses to be produced at the WTP since the integrity of the refractory should not be a limiting factor affecting the lifetime of the melter. To this end, a series of glasses representative of the first batches of waste glass produced by the WTP will be melted in small-scale crucibles with Monofrax® K-3 coupons inserted. The corrosion of the K-3 will be measured in the melt and at the meltline (or neckline). A model for the corrosion rate will be constructed and implemented into a previously developed framework for a computational fluid dynamics (CFD) model of the full-scale WTP. To assist with experimental design and validate the implementation of the model in the full-scale melter, CFD simulations of the small-scale crucible tests were performed. The bubbling that occurs in the small-scale crucible is initially validated here with a model that uses silicone oil at room temperature. The viscosity of the oil ranges from 1 to 100 Pa•s, which corresponds to operating glass pool temperatures near 1150 °C down to idling temperatures near 950 °C. The simulation results show good agreement with the bubble sizes that form during experiments. CFD modeling of the crucible setup was used to determine bubbling characteristics to match the range of near-wall velocities expected in the full-scale WTP. This study presents the initial CFD modeling results, corrosion testing plan, and some preliminary corrosion samples with an outline for the next steps for the development of the corrosion model.

Abboud, Alexander W. [Idaho National Lab]↗
DOE OSTI2024

Tuning collective anion motion enables superionic conductivity in solid-state halide electrolytes

Halides of the family Li 3 MX 6 (M = Y, In, Sc and so on, X = halogen) are emerging solid electrolyte materials for all-solid-state Li-ion batteries. They show greater chemical stability and wider electrochemical stability windows than existing sulfide solid electrolytes, but have lower room-temperature ionic conductivities. Here, in this work, we report the discovery that the superionic transition in Li 3 YCl 6 is triggered by the collective motion of anions, as evidenced by synchrotron X-ray and neutron scattering characterizations and ab initio molecular dynamics simulations. Based on this finding, we used a rational design strategy to lower the transition temperature and thus improve the room-temperature ionic conductivity of this family of compounds. We accordingly synthesized Li 3 YCl x Br 6-x and Li 3 GdCl 3 Br 3 and achieved very high room-temperature conductivities of 6.1 and 11 mS cm -1 for Li 3 YCl 4.5 Br 1.5 and Li 3 GdCl 3 Br 3 , respectively. These findings open new routes to the design of room-temperature superionic conductors for high-performance solid batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Development of Advanced Solid Sorbents for Direct Air Capture (Final Report)

RTI International is partnering with Creare, and Mohammed VI Polytechnic University (UM6P) successfully met all major technical objectives as outlined below: Identified one MOF adsorbent for DAC; Identified one amine-P-dendrimer adsorbent for DAC; Performed Computational fluid dynamics simulations of the MOF and amine-P-dendrimer adsorbents and validate them with experimental data; Selected one adsorbent for DAC; Demonstrated the scale-up of selected candidate and perform cost review; Performed a preliminary process design; Performed technology maturation plan and environmental health and safety risk assessment. After a thorough comparison between both adsorbents (MOFs and P-dendrimers) and based on their technical merits, the amine-based P-dendrimers sorbent was selected as the best performing adsorbent for DAC. In fact, the amine-based P-dendrimers have demonstrated superior performance over the MOF sorbents under relevant DAC conditions. The selected sorbent has been used in subsequent projects (DE-FE0032099 and DE-AR0001412) and has demonstrated similar superior performance with no loss of activity after long-term sorption-desorption testing (e.g., 100 cycles). This sorbent has demonstrated a very high and stable sorption performance, excellent kinetics, and fast regeneration at 80 °C, specifically at 75% relative humidity (RH). The incorporation of RTI’s high-performance, high-durability amine-based P-dendrimers sorbent and Creare’s hybrid additive manufacturing (H-AM) technology that will produce high performance, compact heat and mass exchange structures at low cost using methods, will lead to a contactor that is optimized for wind-driven operation. This contactor technology provides a technically viable pathway for reducing the cost of DAC to <$100/tonne of CO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Supported Organoiridium-Pincer Catalysts for the Nonoxidative Dehydrogenation of High-Density Polyethylene

The iridium-pincer complex {p-OP( t Bu) 2 -C 6 H 2 -2,6-OP( t Bu) 2 ] 2 }Ir(C 2 H 4 ) ( P [Ir]) has been reported as a stable and active catalyst toward alkane dehydrogenation in homogeneous and supported heterogeneous systems. Dehydrogenation has been shown as a practical method toward functional polyolefins, with dehydrogenated high-density polyethylene (deHDPE) demonstrated as a valuable synthon for upcycling, as orthogonal C-H strategies are key to end-of-life upcycling. The heterogenization of P [Ir] on oxides (SiO 2 , Al 2 O 3 , and TiO 2 ; P [Ir]/E y O x ) yields a mixture of organometallic Ir-fragments whose catalytic nonoxidative dehydrogenation activity is modulated by the binding modes of the active metal on the surface. The binding mode was elucidated by a combination of solid-state NMR and XAFS analyses and supported by DFT calculations. Surface binding through the ligand enables active organoiridium that catalyzes internal olefination of deHDPE up to 1.23 mol % at 200 °C under dynamic vacuum. Alternatively, when the organoiridium is bonded though the metal center (Ir-O SiO ), catalyst activity is negligible. Furthermore, the catalytic activity of P [Ir]/SiO 2 showed comparable reactivity with the homogeneous analogue under the same catalytic conditions, and the heterogenized catalyst can be reused up to three cycles. In conclusion, this work highlights the importance of understanding how organometallic precursors react with hydroxylated metal oxide surfaces to establish structure-property relationships.

Hunt, Samuel B. [Argonne National Laboratory (ANL)↗
DOE OSTI2020

Thermomechanical response of thickly tamped targets and diamond anvil cells under pulsed hard x-ray irradiation

In the laboratory study of extreme conditions of temperature and density, the exposure of matter to high intensity radiation sources has been of central importance. Here, we interrogate the performance of multi-layered targets in experiments involving high intensity, hard x-ray irradiation, motivated by the advent of extremely high brightness hard x-ray sources, such as free electron lasers and 4th-generation synchrotron facilities. Intense hard x-ray beams can deliver significant energy in targets having thick x-ray transparent layers (tampers) around samples of interest for the study of novel states of matter and materials’ dynamics. Heated-state lifetimes in such targets can approach the microsecond level, regardless of radiation pulse duration, enabling the exploration of conditions of local thermal and thermodynamic equilibrium at extreme temperature in solid density matter. The thermal and mechanical responses of such thick layered targets following x-ray heating, including hydrodynamic relaxation and heat flow on picosecond to millisecond timescales, are modeled using radiation hydrocode simulation, finite element analysis, and thermodynamic calculations. Assessing the potential for target survival over one or more exposures and resistance to damage arising from heating and resulting mechanical stresses, this study doubles as an investigation into the performance of diamond anvil high pressure cells under high x-ray fluences. Long used in conjunction with synchrotron x-ray radiation and high power optical lasers, the strong confinement afforded by such cells suggests novel applications at emerging high intensity x-ray facilities and new routes to studying thermodynamic equilibrium states of warm, very dense matter.

Meza-Galvez, J. (ORCID:000000032699849X)↗
DOE OSTI2020

Compromise of Personal Protective Clothing from Liquid Exposure

Following critiques of multiple personal contamination events from entries into the Oak Ridge National Laboratory’s Spallation Neutron Source transfer bay, it was considered that the most likely causes for contamination were personal protective clothing doffing errors or that moisture (sweat) allowed contamination to wick through the protective clothing. To pursue this issue, radiological protection staff looked more closely, however, at the specific area of the clothing where contamination was highest; under enhanced lighting and photochromic manipulation, there appeared to have been some type of moisture in the area. Recognizing the possibility that moisture may have allowed for migration of contamination through the clothing, further experiments were undertaken to determine under which conditions this transport might have occurred. The objective for this work was to identify the susceptibility of different types of personal protective clothing to various liquids encountered in the workplace.Method: Several tests were performed to determine if perspiration had enabled migration of contamination, and to identify what other liquids might have affected contamination transport. Two layers of personal protective clothing were subjected to static conditions and dynamic conditions, to include active rubbing of the materials while wet. Food dye added to each of the liquids tested enabled visual indications of liquid breakthrough. Additional tests were conducted to see if solid contamination could be transported through the materials along with the liquids. As a result, it was found that all but one type of non-rubberized personal protective clothing in use at Oak Ridge National Laboratory were permanently compromised to some extent by the solvents used for decontamination. It was determined that most common cleaning agents immediately and permanently destroyed the hydrophobic nature of several of the splash-resistant protective clothing materials, allowing for radioactive contamination to penetrate through the material to the worker. Work around wet surfaces or performing wet decontamination will only be performed in protective clothing known to prevent transport of the wetting agent.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2026

Reverse phonon thermal flux from an applied electric field

In the phenomenon known as electron drag, a phonon thermal flux is established in a conducting crystal by an electric field applied under isothermal conditions through directed transfer of quasimomentum from the electronic charge current to the phonon subsystem. Prior understanding of this phenomenon involves a thermal current composed of low-frequency acoustic phonons that is in the same direction as the charge current. This results in an increase of the Peltier thermopower. Here, we show that it is also possible to establish a phonon thermal current that is in the opposite direction to the charge current, which we refer to as a reverse phonon thermal flux. We demonstrate this behavior both through a simple qualitative model and through first-principles calculations performed for three materials: p-type 𝜃-TaN, n-type BAs, and n-type Si. The reverse phonon flux is shown to arise through intervalley electron-phonon scattering processes involving high-frequency acoustic phonons. Unlike n-BAs and n-Si, the combined features in the band structure and phonon dispersions of p-type 𝜃-TaN promote a particularly large reverse phonon thermal flux comparable in magnitude to the oppositely directed thermal flux from low-frequency phonons. Finally, this work highlights a previously unrecognized behavior in the coupled electron-phonon system and advances our understanding of the rich physics of transport in solids.

lattice dynamics↗
DOE OSTI2023

Convection-Induced Compositional Patterning at Grain Boundaries in Irradiated Alloys

We consider the stability of precipitates formed at grain boundaries (GBs) by radiation-induced segregation in dilute alloys subjected to irradiation. Here, the effects of grain size and misorientation of symmetric-tilt GBs are quantified using phase field modeling. A novel regime is identified where, at long times, GBs are decorated by precipitate patterns that resist coarsening. Maps of the chemical Péclet number indicate that arrested coarsening takes place when solute advection dominates over thermal diffusion right up to the precipitate/matrix interface, preventing interfacial local equilibrium and overriding capillary effects. This contrasts with liquid-solid mixtures where convection always accelerates coarsening.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Attached Jet Velocity Profiles in Mixing Tanks

Free jets have been studied in detail over much of the last century, but the theory for offset and attached jets remains incomplete. Attached jets differ from free jets in that they lose momentum to nearby surfaces, attenuating their velocities. The velocity profiles of free circular jets are nearly Gaussian, with quantitative mathematical descriptions derived from first principles by Goertler and Tollmien (Rajaratnam, 1976). In contrast, mathematical descriptions of three-dimensional attached jets from circular nozzles remain much less mature. Agelin-Chaab and Tachie (2011) used particle imaging velocimetry of a three-dimensional attached jet to show that the scaled velocity decays with scaled distance from the nozzle with a power law exponent between -1.15 and -1.21, which is larger in magnitude than that of a free jet. However, quantitative analytical expressions for the velocity profiles of attached jets similar to those of free jets remain elusive. This paper addresses this critical gap. Here we evaluate the velocity profiles of three-dimensional offset jets emerging from circular nozzles that that become attached jets. These jets lose momentum due to interactions with nearby surfaces and are important to evaluating flows in mixing vessels and to suspending solids and trapped gases in radioactive waste tanks. Despite the importance of attached jets, prior insight has been purely experimental, limited to overly simplistic analytical models, or restricted to computationally expensive computational fluid dynamics case studies. We compare the expression of Verhoff (1963) to experimental results to find reasonable quantitative agreement. As stated by Agelin-Chaab and Tachie (2011), “detailed velocity measurements of 3D offset jets are rare.” Such remains the case. This study adds to the literature by providing information at two additional Reynolds numbers (1.43 106 and 1.87 106) and evaluating simple but accurate expressions for velocity profiles. These Reynolds numbers are an order of magnitude higher than prior reports. The semi-empirical stream wise velocity profile perpendicular to the surface proposed by Verhoff (1963) is in approximate agreement with these velocity profiles, which is surprising because these attached jets are three-dimensional instead of two-dimensional as evaluated by Verhoff. However, additional work is necessary to fully describe these profiles quantitatively.

jet, attached jet, offset jet↗
DOE OSTI2021

Deterministic Calibration of MFiX-PIC, Part 1: Settling Bed

The Particle-in-cell (PIC) numerical approach for modeling granular solids in fluid flow has gained significant interest in recent years. Valued for its often shorter time-to-solution, the PIC formulation relies on modeling statistical groupings of particles called parcels in cooperation with a solids stress model to affect local solids velocity. This is in contrast to the discrete element model (DEM) where every particle in a system is modelled individually and directly coupled to local solids velocity through Newtonian mechanics. The U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) develops and maintains Multiphase Flow with Interphase eXchanges (MFiX), a collection of open-source computational fluid dynamics (CFD) solvers. Included in the MFiX suite are traditional two-fluid model (TFM) and DEM solvers, and a recently added PIC solver (NETL, 2021). In general, PIC methodologies offer an accuracy trade-off in lieu of computational speed; and therefore, it is important to assess the credibility of MFiX-PIC simulations. For this purpose, a systematic verification, validation and uncertainty quantification (VVUQ) effort was initiated at NETL to assess the new PIC solver

42 ENGINEERING↗
DOE OSTI2024

Large-Volume Stimulation of Rock for Greatly Enhanced Fluids Recovery using Targeted Seismic-Assisted Hydraulic Fracturing (Final Technical Report)

This project has developed and demonstrated a new technology for large-volume and targeted comminution of rock in low permeability formations to enhance recovery from unconventional oil and gas (UOG) resources. The technology is based on a strategically designed interaction of multiple induced seismic pulses that assist the hydraulic fracturing process to enhance shear and multi-planar crack formation. This greatly increased rock stimulation, through bulk comminution, is expected to cause significant increase in permeability leading to enhancement of recovery factors for sub- surface fluids. The proposed technology is especially applicable for enhanced recovery in emerging UOG plays, such as ductile shales that are resistant to opening-mode fracturing by conventional hydraulic fracturing processes. The project combines an integrated experimental and computation approach to develop and demonstrate a modular technology that can be easily implement in the field to augment current practices. The effort integrates a fundamental scientific understanding of dynamic material response under constraint and damage-induced permeability and porosity enhancements at multiple length scales, along with models of comminution due to the local release of kinetic energy associated with high shear strain rate of dynamic deformation. The results will be validated through small-scale experiments and then implemented in a lab-scale field test to demonstrate the developed technology. The project involves collaboration between a solid mechanician, a materials scientist, and a petroleum engineering geo-mechanician. It builds on a successful history of success and collaboration that includes both fundamental scientific exploration and technology development.

04 OIL SHALES AND TAR SANDS↗