Nonlinear elastodynamic behavior of intact and fractured rock under in-situ stress and saturation conditions
Not provided.
SEARCH · Engineering Papers
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.
Not provided.
Electron-hole asymmetry is a fundamental property in solids that can determine the nature of quantum phase transitions and the regime of operation for devices. The observation of electron-hole asymmetry in graphene and recently in twisted graphene and moiré heterostructures has spurred interest into whether it stems from single-particle effects or from correlations, which are core to the emergence of intriguing phases in moiré systems. Here, we report an effective way to access electron-hole asymmetry in 2D materials by directly measuring the quasiparticle self-energy in graphene/Boron Nitride field-effect devices. As the chemical potential moves from the hole to the electron-doped side, we see an increased strength of electronic correlations manifested by an increase in the band velocity and inverse quasiparticle lifetime. These results suggest that electronic correlations intrinsically drive the electron-hole asymmetry in graphene and by leveraging this asymmetry can provide alternative avenues to generate exotic phases in twisted moiré heterostructures.
While most performance metrics of high–explosive (HE) based devices like detonation velocity, detonation pressure, and energy output are expected to degrade over time, the evolution of initiation threshold appears less clear, with claims of both increasing and decreasing trends in threshold having been made in the literature. This work analyzes D–optimally designed sequential binary test data for a few thermally conditioned porous–powder and polymer–bonded HE initiator systems using a Bayesian likelihood method employing the probit regression model. Here we find that in most cases the initiation threshold decreases (i.e., sensitivity increases) upon accelerated thermal conditioning. However, such results are nuanced and influenced by factors like the contact area of initiating stimulus, HE characteristics like density and specific surface area, as well as possible thermally induced changes to other materials and interfaces involved.
Abstract CH 4 emissions from inland water are highly uncertain in the current global CH 4 budget, especially for rivers and streams due to sparse measurements and the uncertainty of measurements caused by turbulent water flow. A previous study has revealed that vertical hydrologic exchange flow (VHEF) is the main regulator of CH 4 emissions from riverbed sediments. However, to what extent the understanding obtained from the plot‐scale can be extended to the reach scale and basin scale remains unknown. To address this challenge, we developed a process‐based model to estimate CH 4 flux at the air‐water interface using the attributes available in the national hydrography data set. It calculates the annual mean flux of VHEF, CH 4 production in sediments, and CH 4 transport in the river channel in a sequential manner. Model performance is evaluated by CH 4 efflux observed at the Hanford reach of the Columbia River. We show that reach‐wise sediment hydrologic and biogeochemical conditions estimated from the national hydrography data set could serve as a good indicator of CH 4 emissions from rivers. Aerobic methane oxidation and export to the downstream are the dominant ways of total CH 4 loss for the large lowland river. The hotspots of CH 4 emissions are likely to be at the reaches with fine sediments and slow channel velocity. This study demonstrates the possibility of quantifying CH 4 emissions at the reach scale and the modeling framework has the potential to be extended to the basin scale to improve estimates of CH 4 emissions from lotic inland water.
Bimetallic, reactive multilayers are uniformly structured materials composed of alternating sputter-deposited layers that may be ignited to produce self-propagating mixing and formation reactions. These nanolaminates are most commonly used as rapid-release heat sources. The specific chemical composition at each metal/metal interface determines the rate of mass transport in a mixing and formation reaction. The inclusion of engineered diffusion barriers at each interface will not only inhibit solid-state mixing but also may impede the self-propagating reactions by introducing instabilities to wavefront morphology. Here, this work examines the effect of adding diffusion barriers on the propagation of reaction waves in Co/Al multilayers. The Co/Al system has been shown to exhibit a reaction propagation instability that is dependent on the bilayer thickness, which allows for the occurrence of unstable modes in otherwise stable designs from the inclusion of diffusion barriers. Based on the known stability criteria in the Co/Al multilayer system, the way in which the inclusion of diffusion barriers changes a multilayer's heat of reaction, thermal conductivity, and material mixing mechanisms can be determined. These factors, in aggregate, lead to changes in the wavefront velocity and stability.
The permeability of mesh wicks is important for various applications, including two–phase heat transfer. However, the understanding of the permeability of single–layer, free–standing mesh wicks, with liquid–gas interfaces on both sides, is limited. A novel and simpler method is presented to determine the permeability of a free–standing wick and apply it to a representative mesh. This method involves modifying the capillary pressure via elevation and simultaneously measuring the permeability to determine the permeability–capillary pressure relationship. When applied to a copper mesh with plain weave having undergone surface cleaning, the permeability is found to decrease as capillary pressure for deionized water increases. A dimensional analysis is presented to generalize this data for other mesh sizes with similar weaves and fluids. The behavior of mesh in application is modeled, based on the integration of Darcy's law with an analytic function fit to measured data, and parametric studies are conducted to investigate the superficial velocity of liquids through the mesh under varying driving pressures, transport lengths, and liquid viscosity, based on the obtained capillary pressure–permeability relationship. This study provides valuable insights into the transport properties of mesh wicks, with potential applications in fields such as electronics cooling, electrochemical devices, and fluid purification technologies.
This work investigates the influence of sodium promotion on Ag/m-ZrO 2 catalysts for methanol steam reforming (MSR), focusing on activity, selectivity, surface chemistry, and mechanistic pathways. Temperature programmed reduction (TPR), XANES/EXAFS, CO 2 TPD, DRIFTS, and temperature programmed surface reaction methods were combined with fixed bed MSR testing to develop an integrated structure–function understanding of Na-modified Ag-ZrO 2 interfaces. Na addition systematically increases surface basicity, stabilizes strongly basic O 2− sites, and weakens the ν(CH) vibrational mode of surface formate, thereby facilitating C–H bond scission and accelerating decarboxylation to CO 2 . At moderate promoter levels (0.5–1.0 wt.% Na), the catalysts show significantly enhanced CO 2 selectivity and increased conversion relative to unpromoted Ag/m-ZrO 2 , while CH 4 formation remains negligible. Excessive Na (≥1.8 wt.%) leads to slower formate decomposition, greater carbonate stabilization, and suppressed conversion, revealing a narrow optimum around 1 wt.% Na. Short-term stability testing demonstrates steady conversion and product selectivity for both unpromoted and Na-promoted catalysts, with the latter maintaining markedly higher CO 2 selectivity. Although Pt/YSZ retains far superior intrinsic activity at ~10× higher space velocity, Ag offers a cost-advantaged alternative where lower cost metals are desirable. Collectively, these findings show that Na promotion enables tunable MSR selectivity on Ag/m-ZrO 2 by directing formate decomposition toward the CO 2 -forming pathway.
Three-dimensional non-rotating odd viscous liquids give rise to Taylor columns and support axisymmetric inertial-like waves ( J. Fluid Mech ., vol. 973, 2023, A30). When an odd viscous liquid is subjected to rigid-body rotation however, there arise in addition a plethora of other phenomena that need to be clarified. In this paper, we show that three-dimensional incompressible or two-dimensional compressible odd viscous liquids, rotating rigidly with angular velocity 𝛺, give rise to both oscillatory and evanescent inertial-like waves or a combination thereof (which we call of mixed type) that can be non-axisymmetric. By evanescent, we mean that along the radial direction, typically when moving away from a solid boundary, the velocity field decreases exponentially. These waves precess in a prograde or retrograde manner with respect to the rotating frame. The oscillatory and evanescent waves resemble respectively the body and wall-modes observed in (non-odd) rotating Rayleigh–Bénard convection ( J. Fluid Mech ., vol. 248, 1993, pp. 583–604). We show that the three types of waves (wall, body or mixed) can be classified with respect to pairs of planar wavenumbers 𝜅 which are complex, real or a combination, respectively. Experimentally, by observing the precession rate of the patterns, it would be possible to determine the largely unknown values of the odd viscosity coefficients. This formulation recovers as special cases recent studies of equatorial or topological waves in two-dimensional odd viscous liquids which provided examples of the bulk–interface correspondence at frequencies 𝜔 < 2𝛺. We finally point out that the two- and three-dimensional problems are formally equivalent. Their difference then lies in the way data propagate along characteristic rays in three dimensions, which we demonstrate by classifying the resulting Poincaré–Cartan equations.
We study the linear stability of a planar interface separating two fluids in relative motion, focusing on conditions appropriate for the boundaries of relativistic jets. The jet is magnetically dominated, whereas the ambient wind is gas-pressure-dominated. We derive the most general form of the dispersion relation and provide an analytical approximation of its solution for an ambient sound speed much smaller than the jet Alfvén speed v A , as appropriate for realistic systems. The stability properties are chiefly determined by the angle ψ between the wavevector and the jet magnetic field. For ψ = π/2, magnetic tension plays no role, and our solution resembles the one of a gas-pressure-dominated jet. Here, only sub-Alfvénic jets are unstable ($0\lt {M}_{e}\equiv (v/{v}_{{\rm{A}}})\cos \theta \lt 1$ , where v is the shear velocity and θ the angle between the velocity and the wavevector). For ψ = 0, the free energy in the velocity shear needs to overcome the magnetic tension, and only super-Alfvénic jets are unstable ($1\lt {M}_{e}\lt \sqrt{{(1+{{\rm{\Gamma }}}_{w}^{2})/[1+({v}_{{\rm{A}}}/c)}^{2}{{\rm{\Gamma }}}_{w}^{2}]}$, with Γ w the wind adiabatic index). Our results have important implications for the propagation and emission of relativistic magnetized jets.
This project developed a molten salt pump design for small modular reactors, concentrated solar, and Gen IV nuclear reactors with a magnetically levitated rotor that can operate at temperatures of up to 700 °C. This eliminates the need for rotating seals and roller element bearings which require maintenance and are prone to failure. The pump design is also more compact than existing molten salt pump designs which require long shafts to thermally isolate the electric motor. The aims of the project are to: 1) design and test the molten salt pump, 2) develop and test a high temperature position sensor needed for magnetic bearing control, and 3) develop and test high temperature coil fabrication methods for electric motor and magnetic bearing fabrication. To develop the high temperature coils, conductor materials were analyzed for oxidation, diffusion, cladding, and resistivity. Of the candidate conductor materials that met the requirements, silver was chosen. The long-term impact of radiation on silver resistivity was analyzed by developing a model of a generic molten salt fast reactor to model the conductor irradiation in worst case scenarios. The analysis showed that increased resistivity due to transmutation of silver decreases exponentially with the reactor reflector thickness. Different electrical insulation materials were tested and the final high temperature coils developed were able to withstand conductor to ground voltages up to 1800 V before insulation breakdown. Magnetic bearings require high-precision high-speed measurements of the rotor position to operate. Current commercial position sensors can only operate up to 550 °C so a novel position sensor was developed that can operate up to 800 °C. These sensors were also designed to measure the position of the metallic rotor through a thin metallic corrosion barrier that protects that coils, insulation, and magnetic cores so that the salt containment barrier does not have any penetrations. The high-temperature sensor was tested and had comparable performance to commercially available inductive position sensors. The sensor design and signal processing has led to several invention disclosures. The magnetically levitated molten salt pump design was completed and analyzed. The chloride salt test loop at Oak Ridge National Laboratory was chosen to test the pump and the pump hydraulics and power were designed for the test loop flow and pressure requirements. Classes of motor designs that do not utilize permanent magnetic were considered for the high-temperature motor. After extensive analysis of their performance, a synchronous reluctance motor was chosen for the pump. The design was refined and analyzed to optimize performance. The magnetic bearings were designed to meet the force and frequency requirements necessary to levitate the rotor and reject disturbances from rotor imbalance, motor forces, and hydraulic forces. Compliant mounts for the position sensors were designed that would not damage the ceramic sensor components and maintain sensor alignment over a wide temperature range. A thermal hydraulic analysis of the molten salt coolant flow used to cool the rotor and stator was performed to ensure that the maximum temperature in the rotor and stator will not exceed 725 °C. The corrosion barrier and casing were designed to maintain the salt fluid boundary by using compliant features to absorb large expansion mismatches due to different coefficients of thermal expansion (CTE). Finally, the salt drain and fill systems were designed along with the pump support structure, integration with the test loop, and thermal insulation and heating. Corrosion barrier fabrication experiments were performed to assess using laser cladding to apply a sub-millimeter corrosion barrier to the surface of a Fe-Co-V magnetic core material. The cladding method showed minimal mixing between materials at their interface and withstood temperature cycling without damage to the cladding. The technology commercialization is focused on the high-temperature position sensors and there has been significant interest from the nuclear industry in utilizing these sensors along with some interest from the aerospace industry. The project also led to a research collaboration with a fusion energy industrial partner to study barriers to scaling the pump design from the kW size to the MW size. This research collaboration will also study methods to create high-temperature coils that have conductor to ground breakdown voltages above 5 kV. The successful development of high-temperature electromagnetic coils and high-temperature position sensors greatly increases the temperature limits for electromagnetic devices that don’t rely on permanent magnets. This includes electric motors, magnetic bearings, linear actuators, rotary position and velocity sensing, actuated valves, and generators to name a few. This will fill a current need in molten salt reactors, concentrated solar, and fusion energy for these critical peripheral devices needed for practical reactor designs that do not currently have commercial solutions.
Additive manufacturing (AM) is rapidly emerging in high performance applications such as army ground vehicles, automotive and transportation. However, the response of AM parts/components to extreme loading such as high velocity impacts is less studied. In this work, the performance under ballistic impact of AM panels is evaluated using a medium velocity gas gun, generating projectile velocities up to 400 m/s. The preferential print orientation properties are considered in order to evaluate whether the panels exhibit isotropic or anisotropic behavior under impact. Surface morphology is investigated by milling the beads smooth on samples and comparing the impact on as-printed samples to those that are smoothed. The effect of nickel chromium micron (nichrome) wire embedded in the AM panels (during print) of polycarbonate-carbon fiber (PC-CF) and polycarbonate-glass fiber (PC-GF) are explored. Thermoplastic polyurethane-acrylonitrile butadiene styrene/carbon fiber (TPU-ABS/CF), Acrylonitrile butadiene styrene-carbon-fiber (ABS-CF) AM samples absorbed >50% of the impact energy. The ballistic performance was noted to be in the following order – ABS-CF > TPU-ABS/CF > PC. Scanning electron microscopy (SEM) was conducted to study the interface between the nichrome wire and the polymer-fiber matrix. This work is the first of its kind exploring into the capabilities of AM panels as ballistic materials. This study leads the way for developing AM panels that are easily manufactured and exhibit superior ballistic resistance.
Ductile materials subject to impact loading conditions can undergo spall fracture when an incoming compressive stress wave reflects off interfaces and free surfaces as a tensile stress wave. Experimental observations suggest that in ductile materials, spall fracture is driven by the evolution of porosity. However, the presence of initial porosity in ductile materials also introduces plastic compressibility, which can attenuate the incoming compressive stress wave and, as a result, reduce the amplitude of the reflected tensile stress wave. This, in turn, can mitigate or delay spall fracture. In this work, we report on finite deformation finite element calculations that analyze the response of porous ductile materials subjected to impact loading conditions. Two sets of calculations are carried out, in the first set the material contains initial porosity values ranging from 0% to 5% while in the second set the material also undergoes stress-controlled porosity nucleation. Both sets of calculations are carried out for a wide range of imposed impact velocities. Here, our results show that porosity in ductile materials can, under certain circumstances, mitigate spall fracture by attenuating stress waves. Results correlating the effects of impact velocity, initial porosity, and porosity nucleation on spall fracture are presented and the underlying mechanisms are discussed.
Hydrological exchange flows (HEF) across the river-aquifer interface and the associated residence time of river water in the aquifer have important implications for contaminant plume migration and biogeochemical processes in the river corridor. HEFs and residence time are influenced by both subsurface physical features and hydrologic forcing related to the transport process, which can exhibit complex spatial and temporal variations. In this study, we used a massively parallel subsurface flow model and a particle-tracking model to study the influences of different control factors on spatial variability of HEFs and residence time distributions (RTD) in the Hanford Reach of the Columbia River in Washington State. A total number of 100M particles were randomly injected in time and space and then tracked in a model domain that covers a 51-km 2 area (15.1M model cells). We used hourly river stages and groundwater levels to drive the model to provide dynamic velocity fields for the particle tracking in the simulation period that was longer than 2 years. The groundwater flow simulation and particle-tracking results provide the first comprehensive assessment of the spatial distribution of HEFs and residence time in large complex river corridors. Overall, our results show that the aquifer hydrogeological structure has the strongest correlation with the extent and magnitude of exchange flux. The residence time exhibits complex patterns that are impacted by all the river geomorphologic, hydrodynamic, and hydrogeologic factors and are strongly correlated with the downwelling ratio of exchange flux. The new insights gained through this study can be used to support the development of reduced-order models of HEFs and RTDs for large complex river systems.
The effect of fiber sizing (i.e., surface treatment) concentration (0 %, 1.36 %, 1.52 %, 1.94 %, and 2.13 %) on the mechanical properties (tensile, flexural, interlaminar shear strength (ILSS), and low velocity impact) of textile grade carbon fiber (TCF)-epoxy composite is examined. An inverse relationship between the strength and stiffness of the composite is observed with increased sizing concentration. The root mean square (RMS) roughness of the fiber surface increased from 17.8 nm (unsized) to 22.7 nm with 2.13 % sizing concentration. It was found that the tensile strength increased by 131 % from 221.4 ± 18.5 MPa (unsized) to 510.8 ± 28.05 MPa (for 1.36 % sizing) and further by 155 %–563.7 ± 14.95 MPa at 2.13 % sizing. On the contrary, the stiffness is initially increased by 126 % from 33.52 ± 7.80 GPa (unsized) to 75.9 ± 3.21 GPa (for 1.36 % sizing) but reduced with further increase in the sizing concentration. A single fiber pull-out test is simulated using the finite element method to validate the reverse trend in strength and stiffness. The varying sizing concentration is simulated by introducing an interface of varying thickness between fibers and matrix. Simulation results confirm that a thicker interface, corresponding to a higher sizing concentration, decreases interfacial shear stress, enhancing material strength while decreasing stiffness. The reverse trend in strength and stiffness with the sizing concentration aligns with experimental observations. In conclusion, the present study emphasizes the importance of sizing concentration for mechanical properties and provide a design criterion for customized high-strength and high-stiffness applications.
The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.
Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.
Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.
Seismic full wave inversion (FWI) is a widely used non-linear seismic imaging method used to reconstruct subsurface velocity images, however it is time consuming, has high computational cost and depend heavily on human interaction. Recently, deep learning has accelerated it’s use in several data-driven techniques, however most deep learning techniques suffer from overfitting and stability issues. In this work, we propose an edge computing-based data-driven inversion technique based on supervised deep convolutional neural network to accurately reconstruct the subsurface velocities. Deep learning based data-driven technique depends mostly on bulk data training. In this work, we train our deep convolutional neural network (DCN) (UNet and InversionNet) on the raw seismic data and their corresponding velocity models during the training phase to learn the non-linear mapping between the seismic data and velocity models. The trained network is then used to estimate the velocity models from new input seismic data during the prediction phase. The prediction phase is performed on a resource-constrained edge device such as Raspberry Pi. Raspberry Pi provides real-time and on-device computational power to execute the inference process. In addition, we demonstrate robustness of our models to perform inversion in the presence on noise by performing both noise-aware and no-noise training and feeding the resulting trained models with noise at different signal-to-noise (SNR) ratio values. We make great efforts to achieve very feasible inference times on the Raspberry Pi for both models. Specifically, the inference times per prediction for UNet and InversionNet models on Raspberry Pi were 22 and 4 s respectively whilst inference times for both models on the GPU were 2 and 18 s which are very comparable. Finally, we have designed a user-friendly interactive graphical user interface (GUI) to automate the model execution and inversion process on the Raspberry Pi.