SEARCH · Engineering Papers

Results for “hydrostaticity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Large-scale paleo water-table rise in a deep desert aquifer recorded by dissolved noble gases

Illuminating past hydrological and climatological conditions in arid regions may provide insights into future trends in groundwater availability. The main objective of this study is to explore paleorecharge processes in the deep regional Nubian Sandstone Aquifer (NSA), which stretches below the hyperarid deserts of the Sinai Peninsula (Egypt) and the Negev (Israel), using dissolved noble gas data. Extremely high amounts of dissolved excess air (Ne concentrations up to 4 × solubility equilibrium) were observed in the ancient ( 81 Kr-depleted) groundwater that was recharged during past pluvial epochs in the basin. Here, the observed unique excess air signal and the clear spatiotemporal structure of the noble gas compositions in the aquifer are hypothesized to reflect two major characteristics of the groundwater system: (a) large-scale, long-term rises in the water table that facilitate the entrapment and dissolution of air bubbles under increased hydrostatic pressure, and (b) the geological settings in the southern recharge area in Sinai, including an almost horizontal position of the layers and intercalation of low-permeability rock formations within the aquifer, which inefficiently reject air bubbles from groundwater, thereby allowing for substantial entrapment of bubbles that dissolve as the water table rises. Enhanced dissolution of entrapped air was also observed in other paleo-groundwater-containing regional sandstone aquifers across North Africa, which, together with the recent findings from the NSA, suggest that largescale water table fluctuations have likely extensively occurred during past pluvial periods over these (and possibly other) arid regions. A preliminary assessment of noble gas recharge temperatures (NGTs) indicates an apparent decoupling of surface and water table temperatures in the case of deep aquifers with a thick (hundreds of m) unsaturated zone. This observed decoupling calls for a re-evaluation of previously obtained NGT records and a need for future work to consider the modification of NGTs due to soil air fractionation and geothermal heating of the underlying deep unsaturated zones.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Valence instability and collapse of ferromagnetism in EuB 6 at high pressures

Despite the simplicity of their cubic crystal lattice, rare-earth hexaborides display complex physical properties including a (long debated) onset of metallization via magnetic polaron formation at T c1 approximate to 15 K preceding ferromagnetic ordering at T c2 approximate to 12 K. Here, in this work, we used applied pressure to tune the interplay between electronic structure and magnetism in EuB6. We probed the magnetism, valence, and structure of EuB 6 under quasi-hydrostatic pressures up to 30 GPa using X-ray techniques. Our findings show evidence for collapse of ferromagnetism above 20 GPa following a monotonic increase of mean Eu valence. While X-ray diffraction measurements in the paramagnetic state at room temperature show that the lattice retains cubic symmetry, a measurable quadrupole interaction seen by time-domain synchrotron Mössbauer spectroscopy suggests a lowering of symmetry associated with magnetic ordering, becoming more prominent across the magnetic transition. The interplay between conduction band electron count and magnetism observed under applied pressure in EuB 6 opens possibilities for fine-tuning metallization and magnetic properties of similar Eu-based semi-metal systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Unconventional pressure-induced magnetic transitions and mechanical properties in non-magnetic LaFeSi

Here, this work demonstrates that two consecutive magnetic transitions occur in non-magnetic LaFeSi compound under negative hydrostatic pressure (small volume expansion). Electronic structure properties of LaFeSi were calculated using density functional theory (DFT) to understand the origin of these magnetic transitions. Mechanical properties of LaFeSi, especially elastic properties (shear and Young’s modulus, and Poisson’s ratio), were determined by DFT distortion calculations to elucidate lattice anisotropy and stability. The antiferromagnetic transition at –7.52 GPa, predicted in our work for LaFeSi, is consistent with the reported experimentally observed antiferromagnetic ground states in CeMnSi and LaMnSi compounds having similar crystal geometry. This finding of unexpected magnetic ordering predicts that layered non-magnetic materials, like LaFeSi, may become magnetically active with lattice expansion.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A cavitation and dynamic void growth model for a general class of strain-softening amorphous materials

In this work, we derive a homogenized framework for cavitation and dynamic void growth in a general class of strain-softening materials with particular emphasis on amorphous materials. Analytic solutions are derived by considering an assemblage of spherical shells subject to dynamic hydrostatic tensile pressure. The isochoric response of the matrix material is captured by a free volume theory based amorphous viscoplasticity model. The framework accounts for the coupled rate-dependent, strain-softening behavior typically exhibited in amorphous materials. Validation of the homogenized theory is carried out against numerous direct numerical simulations with excellent agreement. Furthermore, we find that a reduction in cavitation strength and significant acceleration of dynamic void growth rates is induced as a consequence of strain-softening in the amorphous matrix material. Lastly, the model is utilized to interpret some non-intuitive experimental observations of amorphous materials subject to shock compression and subsequent spall failure.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A model of thermal creep and annealing in finite domains based on coupled dislocation climb and vacancy diffusion

Here, we develop a framework to investigate thermal creep and annealing in finite domains, where the climb motion of discrete dislocations is coupled to the diffusion of a continuum vacancy field. The model is first formulated in a continuum finite-deformation setting. All governing equations and boundary conditions are obtained from a unified irreversible thermodynamics principle. The resulting model couples a mechanical boundary value problem (BVP), a vacancy diffusion BVP, and the climb and glide motion of the discrete dislocation network within the crystal. The framework is then linearized for implementation in three-dimensional (3D) discrete dislocation dynamics (DDD) simulations for arbitrary anisotropic crystals. A solution scheme is developed based on the superposition principle, which is imposed weakly on the dislocation network to obtain a Galerkin solution for the nodal climb velocities. The framework includes diffusional (Nabarro–Herring) creep deformation as well as dislocation creep by climb-assisted-glide. The method is applied to simulate the annealing of vacancy loops in Al, with good agreement to experimental measurements by Silcox and Hirsch. We further consider the effects of annealing under stress, and of the proximity of the vacancy loops to loaded and free boundaries Simulations in polycrystalline materials are carried out to highlight the effects of the grain size on dislocation climb and vacancy loop annealing. The method is also applied to estimate the creep rate due to climb-assisted glide of jogged-screw dislocations in γ-TiAl, and results are compared to experiments by Viswanathan et al. Finally, we discuss the effects of uniaxial and hydrostatic stresses on the two diffusive deformation pathways of the material, namely Nabarro–Herring creep and dislocation climb.

36 MATERIALS SCIENCE↗

DOE OSTI2022

High-pressure deformation of metallic glass nanoparticles

Here, the atomistic structure of metallic glasses is closely related to properties such as strength and ductility. Here, Ni 1-x B x metallic glass nanoparticles of two different sizes are compressed under quasi-hydrostatic high-pressure conditions in order to understand structural changes under stress. The structural changes in the nanoparticles were tracked using in situ high-pressure X-ray diffraction (XRD). The ambient pressure pair distribution functions generated from XRD showed that the smaller sized nanoparticles had a more compact amorphous structure with lower coordination number. XRD showed that the amorphous structure was stable up to the maximum pressures achieved. The bulk modulus of the smaller and larger sized nanoparticles was found to be 208 GPa and 178 GPa, respectively. This size-dependent high-pressure behavior was related to compositional differences between the nanoparticles. These results show that Ni 1-x B x metallic glass nanoparticles are highly stable under pressure, which could enable their use as inclusions in metal or ceramic matrix composites.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Directional amorphization of covalently-bonded solids: A generalized deformation mechanism in extreme loading

Shock compression subjects materials to a unique regime of high quasi-hydrostatic pressure and coupled shear stresses for durations on the order of 1–10 nanoseconds for laser-driven loading of samples. There is, additionally, an attendant temperature increase due to the shock and the mechanisms of plastic deformation in metals whereby dislocations, twins, and phase transitions nucleate and propagate at velocities near the sound speed. Covalently bonded materials have, by virtue of the directionality of their bonds, great difficulty in responding by conventional plastic deformation to this extreme regime of shock compression. Here, we propose that the shear from shock compression induces amorphization, as observed in Si, Ge, B4C, SiC, and olivine ((Mg, Fe) 2 SO 4 ) and that this is a general deformation mechanism in a broad class of covalently bonded materials. The crystalline structure transforms to amorphous along regions of maximum shear stress, forming nanoscale bands, and thereby relaxing the shear component of the imposed shock stress. This process is usually preceded by the emission and propagation of a critical concentration of dislocations.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Horn Location Sensors for LBNF

The Horn Location Sensor (HLS) system—named for the magnetic focusing horns it monitors, along with other critical beamline elements – is a high-precision alignment system developed for the Long-Baseline Neutrino Facility (LBNF) to support the Deep Underground Neutrino Experiment (DUNE). With minimal maintenance, the HLS can operate reliably in environments with high radiation, and it ensures that key components such as the protective baffle, focusing horns, and beam position monitors are aligned correctly – each essential for maintaining one of the world’s most intense muon-neutrino beams. A high-precision hydrostatic level sensor, a linear variable differential transformer, and INVAR rods are all used in the system to monitor vertical motion and tilt with sub-millimeter accuracy. A novel Sweep Tracker interferometer enhances calibration fidelity by correcting for non-linearities in laser wavelength and scan rate in real time. A critical part of DUNE’s precision alignment and flux prediction requirements, the HLS system initially supports beam power of up to 1.2 MW and can be upgraded to 2.4 MW.

Frequency scanning interferometry↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

3D Poisson equation↗

DOE OSTI2023

Turbulence modeling to aid tidal energy resource characterization in the Western Passage, Maine, USA

Numerical models combined with field measurements are regularly used to characterize tidal energy resources at potential energetic sites. However, most existing works only focus on the tidal hydrodynamic characteristics, and turbulence parameters are often not reported because of the lack of high-quality turbulence measurements and the limitations of numerical models in resolving turbulent eddies. In this study, we used FVCOM - a hydrostatic primitive equation (HPE) model - to characterize the tidal energy resource in the Western Passage, Maine, USA, by taking care of the essential macro-scale turbulence properties. We observed an excellent model performance using the Mellor-Yamada Level 2.5 Turbulence Model; estimating the spatial and vertical distribution of the turbulent kinetic energy and intensity added a new perspective to the site ranking for tidal energy converter (TEC) deployments. In addition, we also examined the role of channel geometry and bathymetry, such as headlands and underwater sills, in enhancing turbulent eddies around potential TEC siting locations. Ultimately, the detailed analysis of the turbulent flow characteristics has changed the site-ranking results and demonstrated that the regional-scale HPE models could be used for the relative understanding of more or less turbulent sites for a refined resource assessment.

16 TIDAL AND WAVE POWER↗

DOE OSTI2020

Turbulent exchanges between near-inertial waves and balanced flows

Wind generated near-inertial waves are ubiquitous in the upper ocean. An improved understanding of near-inertial wave dynamics following their excitation in the ocean and their subsequent interaction with mesoscale geostrophic balanced flows is key to decoding oceanic energy flow pathways. In this regard, multiple oceanic data sets accumulated over the past few decades reveal that the relative strength of near-inertial waves and geostrophic balanced eddy fields is highly variable, both geographically and seasonally. Inspired by these observations, we investigate turbulent interactions and energy exchanges between near-inertial waves and balanced flows using freely evolving numerical simulations of the non-hydrostatic Boussinesq equations. Here, we find accelerated vertical propagation and dissipation of the waves in regimes where balanced and wave fields have comparable strengths. In such regimes we also find that near-inertial waves directly extract energy from balanced flows, with O(10%) being the amount of balanced energy loss. In contrast, we find that near-inertial waves transfer energy to balanced flows in regimes where balance-to-wave energy ratio is small, with the gain in balanced energy being dependent on the relative strength of waves. Furthermore, these regimes are characterized by relatively weaker vertical propagation and dissipation of the near-inertial wave field. One of the key outcomes of this study is the demonstration of the lack of a unique direction for near-inertial wave-balanced flow energy transfers. Depending on the balance-to-wave energy ratio, near-inertial waves can act as an energy sink or energy source for the geostrophic balanced flow.

58 GEOSCIENCES↗

DOE OSTI2020

Small-angle Neutron Scattering (SANS) Characterization of Clay- and Carbonate-rich Shale at Elevated Pressures

Unconventional oil and gas from shale formations have emerged as some of the fastest growing energy resources in the United States, providing both cleaner energy to consumers and reducing the nation’s reliance on energy imports. To properly harness these important natural resources, the nanopore structure of associated shales must be fully understood, particularly under hydraulic fracturing conditions, where they are exposed to both overburden compressive and hydrostatic fluid pressures. The current study uses small-angle neutron scattering (SANS) to characterize pore structure, including porosity, pore accessibility, and pore size distribution, in the 1–100 nm regime at elevated pressures for mineralogically distinct clay- and carbonate-rich shales from the Permian Basin. Unlike typical porosity measurement techniques, SANS is uniquely capable of characterizing both open and closed porosity, allowing measurement of how pore accessibility changes with pressure and determination of the size range of accessible versus inaccessible pores. The porosity of the clay-rich shale was 7.7%, compared to 0.51% for the carbonate-rich shale. However, only 2.6% of the nanopores in the carbonate-rich shale were inaccessible to water at 8 kPSI (55.1 MPa) compared to 7.8% for the clay-rich shale. Furthermore, our analyses indicated that the closed pores fall within distinct size ranges, likely corresponding with the chemical nature of the pore host material. These results provide valuable insight into the effects of shale petrophysical properties on hydrocarbon extraction from unconventional reservoirs.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗

DOE OSTI2020

Modeling Molecules under Pressure with Gaussian Potentials

The computational modeling of molecules under high pressure is a growing research area that augments experimental high-pressure chemistry. Here, a new electronic structure method for modeling atoms and molecules under pressure, Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) approach, is introduced. In this method, a set of Gaussian potentials is distributed evenly on the van der Waals surface of the investigated chemical system, leading to a compression of the electron density and the atomic scaffold. Since no parameters other than pressure need to be specified, GOSTSHYP allows straightforward geometry optimizations and ab initio molecular dynamics simulations of chemical systems under pressure for nonexpert users. In this work, calculated energies, bond lengths, and dipole moments under pressure fall within the range of established computational methods for high-pressure chemistry. A Diels-Alder reaction and the cyclotrimerization of acetylene showcase the ability of GOSTSHYP to model pressure-induced chemical reactions. The connection to mechanochemistry is pointed out.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Phase Transition of Heptane under Pressure

A solid-to-solid phase transition was discovered in n-heptane with an onset at a pressure of 1.7 ± 0.1 GPa, evidenced by a discontinuity in the pressure–volume equation of state obtained from synchrotron X-ray diffraction experiments conducted at room temperature. The phase transition also resulted in discontinuous shifts of vibrational peaks in the infrared spectra. No additional phase transitions were detected up to 10 GPa, contradicting an earlier claim of a phase transition at 3 GPa. Furthermore, the hydrostatic limit of heptane in a diamond anvil cell was estimated to be 3.5 GPa. Above that pressure, the stress state is not well-defined, which may explain prior reports of a phase transition at 7.5 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Effect of Pressure on Crystal Structure and Phonon Density of States of FeSi

The strongly correlated material FeSi displays several unusual thermal, magnetic, and structural properties under varying pressure-temperature (P-T) conditions. It is a potential thermoelectric alloy and a material with several geochemical implications as a possible constituent at the Earth's core-mantle boundary (CMB). Previous theoretical studies predicted a pressure-induced B20-B2 transition at ambient temperature below 40 GPa; however, experimentally, the structural transition is observed only at high P-T conditions. In this study, we have performed high-pressure powder X-ray diffraction (XRD) up to 90 GPa and Nuclear Resonant Inelastic X-ray Scattering (NRIXS) measurements up to 120 GPa to understand the phase stability and lattice dynamics. Our study provides evidence for a non-hydrostatic stress-induced B20-B2 transition in FeSi at around 36 GPa. We deduced the Fe partial phonon density of states (PDOS) and thermal parameters from NRIXS measurements up to 120 GPa and compared them with density functional theory (DFT) calculations. Furthermore, the computations show pressure-induced metallization and the band gap closing around 12 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗