Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Instability control”

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

Synthesis, Properties, and Electrochemical Proton Reduction of a Homoleptic Tetrathiolato Ni-Site Model of [NiFe]-Hydrogenase

[NiFe]-hydrogenase enzymes process H 2 at a nonplanar tetracysteinato-Ni site, the sole participator in proton binding/redox chemistry during turnover. With the objective of assessing whether a simple tetrahedral/ tetrathiolato-Ni 2+ core could promote H 2 evolution reaction (HER), we synthesized (Et 4 N) 2 [Ni(S-p-CF 3 −Ph) 4 ] (1) employing para-trifluoromethylbenzenethiolate ( − S-p-CF 3 −Ph) as a Ni-site analog of [NiFe]-hydrogenase. Spectroscopic measurements and X-ray crystallography confirm the distorted tetrahedral geometry of 1. Dissolution of 1 results in partial thiolate dissociation and formation of S,S-bridged complexes such as (Et 4 N) 2 [Ni 2 (S-p-CF 3 −Ph) 6 ] (3) among other ill-defined species. Dissociation is further accelerated in the presence of Brønsted acids, complicating the assessment of 1 for proton reduction. However, this dissociation/proton instability is suppressed in the presence of additional thiolate ligand to ensure tetrahedral/tetrathiolato 1 persists in solution. Electrochemical HER activity was evaluated by monitoring the current response of an MeCN solution of 1/ excess thiolate after sequential titration with a weak Brønsted acid (acetic acid). The results suggest that 1 is a modest electrocatalyst for the HER with a turnover frequency of 14.5 ± 3.6 s −1 and an overpotential of 0.72 ± 0.02 V. Control experiments and supplementary DFT computations indicate that 1, or a species derived from 1, is responsible for the HER and suggest an ECCE-type mechanism.

Evolution reactions↗
DOE OSTI2019

Discovery of a chemical probe for PRDM9

PRDM9 is a PR domain containing protein which trimethylates histone 3 on lysine 4 and 36. Its normal expression is restricted to germ cells and attenuation of its activity results in altered meiotic gene transcription, impairment of double-stranded breaks and pairing between homologous chromosomes. There is growing evidence for a role of aberrant expression of PRDM9 in oncogenesis and genome instability. Here we report the discovery of MRK-740, a potent (IC 50 : 80 ± 16 nM), selective and cell-active PRDM9 inhibitor (Chemical Probe). MRK-740 binds in the substrate-binding pocket, with unusually extensive interactions with the cofactor S-adenosylmethionine (SAM), conferring SAM-dependent substrate-competitive inhibition. In cells, MRK-740 specifically and directly inhibits H3K4 methylation at endogenous PRDM9 target loci, whereas the closely related inactive control compound, MRK-740-NC, does not. The discovery of MRK-740 as a chemical probe for the PRDM subfamily of methyltransferases highlights the potential for exploiting SAM in targeting SAM-dependent methyltransferases.

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

Origin and regulation of oxygen redox instability in high-voltage battery cathodes

Oxygen redox at high-voltage has emerged as a transformative paradigm for high-energy battery cathodes such as layered transition-metal oxides by offering extra capacity beyond conventional transition-metal redox. However, these cathodes suffer from voltage hysteresis, voltage fade, and capacity drop upon cycling; single-crystalline cathodes have recently shown some improvements but these challenges still remain. We reveal the fundamental origin of oxygen redox instability to be originated from the domain boundaries that are present in single-crystalline cathode particles. By synthesizing single-crystalline cathodes free of domain boundaries, we show that the elimination of domain boundaries enhances the reversible lattice oxygen redox while inhibiting the irreversible oxygen release. Further, this leads to significantly suppressed structural degradation and improved mechanical integrity during battery cycling and abuse heating. The robust oxygen redox enabled through domain boundary control provides practical opportunities towards high-energy, long-cycling, and safe batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

3D Printed Carbon Aerogels for Polymer-Electrolyte Fuel Cells

Carbon aerogels (CA) are known for ultra-low density, very high surface area, fine porous structure, and excellent electrical conductivity. Recent advances in high resolution 3D printing at LLNL enable the manufacturing of complex structured CA electrodes which have been demonstrated in electrolysis, flow battery, and supercapacitor applications. Polymer-electrolyte fuel cells (PEFC) may be an excellent candidate for these materials/methods as designers seek to optimize membrane-electrode assemblies to achieve higher efficiencies. A major challenge to PEFC design is two-phase flow due to water production which causes pore flooding and leads to cell instability and degradation. To mitigate these effects, groups have investigated the creation of liquid specific pathways in the porous layers to alleviate liquid buildup and enhance gas diffusion. While some progress has been made, traditionally manufactured GDLs, which consist of randomly distributed carbon fibers, are limiting. 3D printed CAs are a promising solution as complex macroporous electrode structures can be manufactured out of micro and nano porous filaments with a high level of control. Additionally, the process may allow for the consolidation of multiple layers (MPL, GDL, & flow channels) into a single printed architecture. In this work various CA lattice structures were printed using a direct ink write (DIW) method and liquid breakthrough pressures were measured to demonstrate the ability to tune the water pathways to achieve specific pressure thresholds. Furthermore, a discussion of the prospects for printed CA electrodes in PEFC will focus on material properties, synthesis techniques, and current challenges.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Numerical Simulations of Precipitation Processes, Microphysics, and Microwave Radiative Properties of flood Producing Storms in Mediterranean & Adriatic Basins

A comprehensive understanding of the meteorological and microphysical nature of Mediterranean storms requires a combination of in situ data analysis, radar data analysis, and satellite data analysis, effectively integrated with numerical modeling studies at various scales. An important aspect of understanding microphysical controls of severe storms, is first understanding the meteorological controls under which a storm has evolved, and then using that information to help characterize the dominant microphysical processes. For hazardous Mediterranean storms, highlighted by the October 5-6, 1998 Friuli flood event in northern Italy, a comprehensive microphysical interpretation requires an understanding of the multiple phases of storm evolution. This involves intense convective development, Sratiform decay, orographic lifting, and sloped frontal lifting processes, as well as the associated vertical motions and thermodynamical instabilities governing physical processes that effect details of the size distributions and fall rates of the various types of hydrometeors found within the storm environment. This talk overviews the microphysical elements of a severe Mediterranean storm in such a context, investigated with the aid of TRMM satellite and other remote sensing measurements, but guided by a nonhydrostatic mesoscale model simulation of the Friuli flood event. The data analysis for this paper was conducted by my research groups at the Global Hydrology and Climate Center in Huntsville, AL and Florida State University in Tallahassee, and in collaboration with Dr. Alberto Mugnai's research group at the Institute of Atmospheric Physics in Rome. The numerical modeling was conducted by Professor Oreg Tripoli and Ms. Giulia Panegrossi at the University of Wisconsin in Madison, using Professor Tripoli's nonhydrostatic modeling system (NMS). This is a scalable, fully nested mesoscale model capable of resolving nonhydrostatic circulations from regional scale down to cloud scale and below.

Smith, Eric A.↗
NASA NTRS2009

Application of Stereo PIV on a Supersonic Parachute Model

The Mars Science Laboratory (MSL) is the next step in NASA's Mars Exploration Program, currently scheduled for 2011. The spacecraft's descent into the Martian atmosphere will be slowed from Mach 2 to subsonic speeds via a large parachute system with final landing under propulsive control. A Disk-Band-Gap (DBG) parachute will be used on MSL similar to the designs that have been used on previous missions, however; the DBG parachute used by MSL will be larger (21.5 m) than in any of the previous missions due to the weight of the payload and landing site requirements. The MSL parachute will also deploy at higher Mach number (M 2) than previous parachutes, which can lead to instabilities in canopy performance. Both the increased size of the DBG above previous demonstrated configurations and deployment at higher Mach numbers add uncertainty to the deployment, structural integrity and performance of the parachute. In order to verify the performance of the DBG on MSL, experimental testing, including acquisition of Stereo Particle Imaging Velocimetry (PIV) measurements were required for validating CFD predictions of the parachute performance. A rigid model of the DBG parachute was tested in the 10x10 foot wind tunnel at GRC. Prior to the MSL tests, a PIV system had never been used in the 10x10 wind tunnel. In this paper we discuss some of the technical challenges overcome in implementing a Stereo PIV system with a 750x400 mm field-of-view in the 10x10 wind tunnel facility and results from the MSL hardshell canopy tests.

Wernet, Mark P.↗
NASA NTRS2018

Application of a First-Principles Anomalous Transport Model for Electrons to Multiple Hall Thrusters and Operating Conditions

We have developed a physics-based model based on a pseudo-particle description of the electron cyclotron drift instability. A key improvement of the model with respect to previous work is that linear theory is not applied in the event of wave saturation and deviations of electrons or ions from a Maxwellian distribution function. In the acceleration region, the anomalous collision frequency is computed as the minimum value necessary to prevent the electron drift velocity from exceeding the thermal velocity. A functional based on the electron equilibration time is defined to control the transition from high to low resistivity regions. The model was previously applied to a single Hall thruster at its nominal operating condition, showing promising results that captured accurately the location of the thruster’s acceleration region. In this paper, we extend the use of this firstprinciples models to two additional thrusters, also considering multiple operating conditions for each of them. Numerical results are compared to experimental measurements obtained with non-invasive laser induced fluorescence. In general, the agreement between experiments and simulations is good. The model is able to predict the location of the acceleration region for all cases. We observe however that fine details, such as changes in the plasma potential gradient within the acceleration regions, are not captured. The model is also insensitive to changes in the magnetic field strength while experiments show that small shifts in location (of less than 5% of the acceleration channel length) occur. We plan to address the weaknesses of our method with the help of physical insight gained from kinetic simulations of the acceleration region.

Chaplin, Vernon H.↗
NASA NTRS1994

Low-freqency instabilities and the resulting velocity distributions of pickup ions at comet Halley

The interaction between the solar wind and newborn cometary ions is studied using a new analytical theory as well as one- and two-dimensional hybrid simulations. Using the observed parameters upstream of the comet Halley, a detailed study of wave excitation and the resulting particle distributions is presentd. Linear theory as well as simulations show that a variety of modes such as the fast magnetosonic mode, high frequency whistlers and obliquely propagating Alfven ion cyclotron waves can be excited. However, parallel propagating waves are found to be dominant in the wave spectrum and to control the scattering of the pickup ions. Several features of the observed distributions of pickup protons are explained. In particular, it is shown that the observed asymmetric pitch distribution for the pickup protons is due to the small saturation amplitude of the waves for the given parameters. Water group associated waves can lead to energy diffusion and further pitch angle scattering of protons. This effect is most likely to be important in the vicinity of the bow shock of comet Halley where the density of water group ions becomes comparable to that of protons. It is shown that the observed increase in the radius of the proton velocity shell just outside the bow shock can be due to water group waves. The nearby isotropic proton pitch angle distribution observed by Neugebauer et al. (1989) just outside the bow shock may, however, be related to the presence of a rotational discontinuity which has been identified in the magnetic field data. Just outside the bow shock, simulations show that parallel propagating water group waves can steepen with attached whistler wave packets. The steepening process at parallel propagation is a transient effect, in an important contrast to the case of steepening at oblique angles. The smaller beam densitites at comet Halley appears to be the main reason not only why waves at comet Halley have smaller amplitudes but also why oblique, steepening magnetosonic waves have not been detected at comet Halley, whereas they have been seen at comet Giacobini-Zinner. Halley, whereas they have been seen at comet Giacobini-Zinner.

Karimabadi, H.↗
NASA NTRS1987

Arcjet thruster research and technology, phase 1

The objectives of Phase 1 were to evaluate analytically and experimentally the operation, performance, and lifetime of arcjet thrusters operating between 0.5 and 3.0 kW with catalytically decomposed hydrazine (N2H4) and to begin development of the requisite power control unit (PCU) technology. Fundamental analyses were performed of the arcjet nozzle, the gas kinetic reaction effects, the thermal environment, and the arc stabilizing vortex. The VNAP2 flow code was used to analyze arcjet nozzle performance with non-uniform entrance profiles. Viscous losses become dominant beyond expansion ratios of 50:1 because of the low Reynolds numbers. A survey of vortex phenomena and analysis techniques identified viscous dissipation and vortex breakdown as two flow instabilities that could affect arcjet operation. The gas kinetics code CREK1D was used to study the gas kinetics of high temperature N2H4 decomposition products. The arc/gas energy transfer is a non-equilibrium process because of the reaction rate constants and the short gas residence times. A thermal analysis code was used to guide design work and to provide a means to back out power losses at the anode fall based on test thermocouple data. The low flow rate and large thermal masses made optimization of a regenerative heating scheme unnecessary.

Knowles, Steven C.↗
DOE OSTI2022

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1988

Effect of manipulator and feel system characteristics on pilot performance in roll tracking

The control system features of force vs displacement stick command, feel/force displacement gradient, 'feel' system vs command prefilter dynamic lag, and flight control system effective time delay, have been investigated by a fixed-base simulation employing both joysticks and side-sticks in order to ascertain their effect on pilot opinion, pilot performance, and their hand/arm neuromuscular subsystem interaction. The results obtained indicate that the use of displacement sensing added the feel system dynamic lag to the command path and increased effective system flight delay, thereby lowering pilot opinion and performance but reducing the possibility of roll instability; the use of force sensing as stick command to exclude feel system dynamic lag had the obverse effect.

Aponso, Bimal L.↗
NASA NTRS2000

Effects of G-Jitter on Interfacial Dynamics of Two Miscible Liquids: Application of MIM

We designed an experiment to examine the effects of g-jitter on mixing of two miscible liquids using the Microgravity Vibration Isolation Mount (MIM). The global bifurcation of the interface was observed with the MIM operating alternatively to either transmit the g-jitter, isolate from the g-jitter or to provide controlled vibration levels with well defined amplitude and frequency content. With the MIM in isolation mode, the interface remains stationary indicating buoyancy induced convection is negligibly small such that mixing occurs via intrinsic mass diffusion without the masking effect of vibration driven convection. Analytical and computational results are in agreement with the experimental findings. Operation of the MIM in forced mode with conditions typical of g-jitter shows that vibration induced convective flows can excite instability mechanisms such as Kelvin-Helmholtz to generate large amplitude quasi-stationary waves oriented vertically for various cases with Stokes-Reynolds number in the range of 0.003 to 0.5. The two and four mode quasi-stationary waves are also predicted with a mathematical model. Though unplanned, the effect of a primary thruster filing was captured and shown to cause a catastrophic bifurcation, enhancing local mass transport. In light of the findings, experiments planned for the International Space Station should consider the potential effects of g-jitter.

Duval, Walter M. B.↗
DOE OSTI2022

Efficacy of Stabilizing Calcium Battery Electrolytes through Salt-Directed Coordination Change

Achieving practical, high-energy-density calcium batteries requires controlling the stability of Ca 2+ electrolytes during calcium metal cycling. Because of the highly reactive nature of calcium, most typical electrolyte constituents are unstable, leading to electrode passivation and low Coulombic efficiency. Among various commercially available salts, calcium bis(trifluoromethylsulfonyl)imide (Ca(TFSI) 2 ) is attractive because of its oxidative stability and high solubility in a variety of solvents. However, this salt does not allow for calcium metal plating, and it has been proposed that TFSI – instability induced by Ca 2+ coordination is to blame. In this work, we test the ability of strongly coordinating Ca 2+ cosalts such as halides and borohydrides to displace TFSI – from the first coordination shell of Ca 2+ and thereby stabilize TFSI-based electrolytes to enable calcium plating. Through spectroscopic analysis, we find that the effectiveness of these cosalts at displacing the TFSI – anion is dependent on the solvent’s coordination strength toward Ca 2+ . Surprisingly, electrochemical calcium deposition behavior is not correlated to the population of bound or free TFSI – . Instead, the nature of the coordination interaction between Ca 2+ and the cosalt anion is more important for determining stability. Our findings indicate that TFSI – anions are inherently unstable during calcium deposition even in the nominally free state. Furthermore, strategies aimed at eliminating the interactions of these anions with the electrode surface via interface/interphase design are required.

25 ENERGY STORAGE↗
DOE OSTI2025

Asymptotic scaling laws for the stagnation conditions of Z-pinch implosions

Implosions of magnetically driven annular shells (Z pinches) are studied in the laboratory to produce high-energy-density plasmas. Such plasmas have a wide-range of applications including x-ray generation, controlled thermonuclear fusion, and astrophysics studies. In this work, we theoretically investigate the in-flight dynamics of a magnetically driven, imploding cylindrical shell that stagnates onto itself upon collision on axis. The converging flow of the Z-pinch is analyzed by considering the implosion trajectory in the (A, M) parametric plane, where A is the in-flight aspect ratio and M is the implosion Mach number. For an ideal implosion in the absence of instabilities and in the limit of A ≫ 1, we derive asymptotic scaling laws for hydrodynamic quantities evaluated at stagnation (e.g., density, temperature, and pressure) and for performance metrics (e.g., soft x-ray emission, K-shell x-ray emission, and neutron yield) as functions of target-design parameters.

Ruiz, D. E. [Sandia National Laboratories (SNL-NM)↗
DOE OSTI2025

Rippled metamaterials with scale-dependent tailorable elasticity

Thermally induced ripples are intrinsic features of nanometer-thick films, atomically thin materials, and cell membranes, significantly affecting their elastic properties. Despite decades of theoretical studies on the mechanics of suspended thermalized sheets, controversy still exists over the impact of these ripples, with conflicting predictions about whether elasticity is scale-dependent or scale-independent. Experimental progress has been hindered so far by the inability to have a platform capable of fully isolating and characterizing the effects of ripples. This knowledge gap limits the fundamental understanding of thin materials and their practical applications. Here, we show that thermal-like static ripples shape thin films into a class of metamaterials with scale-dependent, customizable elasticity. Utilizing a scalable semiconductor manufacturing process, we engineered nanometer-thick films with precisely controlled frozen random ripples, resembling snapshots of thermally fluctuating membranes. Resonant frequency measurements of rippled cantilevers reveal that random ripples effectively renormalize and enhance the average bending rigidity and sample-to-sample variations in a scale-dependent manner, consistent with recent theoretical estimations. The predictive power of the theoretical model, combined with the scalability of the fabrication process, was further exploited to create kirigami architectures with tailored bending rigidity and mechanical metamaterials with delayed buckling instability.

Applied Physical Sciences↗
NASA NTRS2003

Minimizing Segregation during the Controlled Directional Solidification of Dendric Alloys

Gravity-driven convection induced in the liquid by density gradients of temperature or composition disrupts uniform dendritic growth during controlled directional solidification and promotes severe macrosegregation. The solute-rich region about the dendrite tip appears to play a pivotal role in channel initiation. Allen and Hunt referred to this region as an "initial transient" or dynamic region constituting steep concentration gradients. Experimental investigation also point to the role the tip region plays in developing microstructure. Hellawell and co-workers showed that flow-through dendritic channels could be effectively disrupted, and segregation minimized, during the gradient freezing of bulk castings by rotating the melt through a slight angle with respect to Earth's gravity vector. Adapting this principle to controlled directional solidification, it has been shown" that segregation in dendritic alloys can be minimized, and properties improved, by processing the sample near horizontal in conjunction with a slow axial rotation of the crucible. It is postulated that the observed microstructural uniformity arises by maintaining the developing solute field about the dendrite tip. Solute rejected during vertical directional solidification will rise or sink parallel to the primary dendrite arms during axial rotation setting the stage for accumulation, instabilities, and segregation. In contrast, during horizontal growth, the rejected solute will sink or rise perpendicular to the primary dendrite. Now, in the presence of a slight axial rotation, solute that was initially sinking (or rising) will find itself above (or below) its parent dendrite, i.e., still about the tip region. The following is intended to experimentally demonstrate the viability of this concept in coordination with a model that gives predictive insight regarding solute distribution about growing dendrites. Alloys based on the lead-tin eutectic system were used in this study. The system is well characterized, the constituent metals are available in a very pure form, and the thermophysical properties are well known. During solidification of hypoeutectic alloys, e.g., 55 wt pct Pb, the primary dendrites reject the less dense tin, and for the hypereutectic alloys, e.g., 75 wt pct Sn, the primary dendrites reject denser lead. Alloys were prepared by melting appropriate amounts of lead and tin in a glass crucible after which the homogeneous liquid was sucked directly into 5-mm i.d. glass tubes. The sample tube, containing approximately 30 cm of alloy, was then mechanically driven into the directional solidification furnace assembly and positioned such that approx. 20 cm of the sample was remelted. Subsequently, directional solidification was initiated by withdrawing the sample through a water-cooled jacket at a constant growth velocity of 2 ,microns/s. After 5 to 6 cm of growth, the sample was quickly removed from the furnace and quenched in a water bath to preserve the solid-liquid interface. Samples were directionally solidified vertically upward, nearly horizontally, and some in conjunction with an applied axial rotation of the crucible. Temperature gradients at the solid-liquid interface were measured with an in-siru K-type thermocouple. Solidified samples were cut perpendicular and parallel to the growth direction and conventionally prepared for microscopic examination.

Grugel, Richard N.↗
NASA NTRS1999

Control of Flowing Liquid Films by Electrostatic Fields in Space

The interaction of a spacially varying electric field and a flowing thin liquid film is investigated experimentally for the design of a proposed light weight space radiator. Electrodes are utilized to create a negative pressure at the bottom of a fluid film and suppress leaks if a micrometeorite punctures the radiator surface. Experimental pressure profiles under a vertical falling film, which passes under a finite electrode, show that fields of sufficient strength can be used safely in such a device. Leak stopping experiments demonstrate that leaks can be stopped with an electric field in earth gravity. A new type of electrohydrodynamic instability causes waves in the fluid film to develop into 3D cones and touch the electrode at a critical voltage. Methods previously used to calculate critical voltages for non moving films are shown to be inappropriate for this situation. The instability determines a maximum field which may be utilized in design, so the possible dependence of critical voltage on electrode length, height above the film, and fluid Reynolds number is discussed.

Griffing, E. M.↗
DOE OSTI2023

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗