Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic coupling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
DOE OSTI2026

Mixed-anion electrolytes: from bulk speciation to interfacial dynamics in divalent metal electrodeposition

Mixing anions is emerging as a promising strategy for multivalent electrolyte design, allowing for adjustment of the solvation structure of bulk cations and enhancing the efficiency of electrochemical processes (e.g. metal deposition for batteries and catalysis). Further progress in electrolyte development requires a fundamental understanding of how tailored electrolyte speciation in mixed anion systems can modify the dynamic electrochemical interface during metal cycling. In this study, we present an anode-focused mechanistic study of exemplar Mg electrolytes containing three different secondary anions, correlating electrochemical behavior with bulk speciation and operando interfacial dynamics. Electrospray Ionization-Mass Spectrometry (ESI-MS) results reveal a general trend of forming mixed anion contact ion pairs (CIPs) across various anions, with the extent of ion pairing influenced by the association strength of the secondary anion. Operando multiharmonic electrochemical quartz crystal microbalance with dissipation (EQCM-D) reveals how these bulk species influence interfacial mass uptake, viscoelasticity, and solvent-coupled hydrodynamic behavior during deposition and stripping. The results indicate that Mg-containing ion pairs and solvated complexes shape adsorption, nucleation, and deposit growth, leading to distinct anion-dependent interphases ranging from more permeable, solvent-coupled layers to relatively compact and rigid deposits. This work establishes a quantitative link between bulk speciation and interfacial dynamics in divalent metal electrodeposition and provides mechanistic guidance for electrolyte design.

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

Simulation study of chemo-mechanical impacts of CO 2 injection in morrow b sandstone reservoir

Here, this study investigates the impacts of geomechanical and geochemical changes on carbon storage in a partially depleted oil reservoir, using results from four different coupled simulation models. Models were used to examine the relative importance of storage mechanisms, and how changing reservoir parameters might affect these mechanisms through time. The study uses data from a Morrowan sandstone reservoir in the Farnsworth Unit (FWU), Ochiltree County, Texas which is currently undergoing CO 2 enhanced oil recovery (EOR). Partially depleted oil reservoirs such as the FWU offer attractive carbon utilization and/or storage targets because of existing infrastructure and economic benefits from incremental oil recovery as well as tax credits. However, prediction of storage capacity or long-term fluid migration in these fields can be difficult because of the wide variation in formation fluids and operational histories that may have undergone. CO 2 injection can cause complex geomechanical and geochemical responses in a reservoir as a result of interplay between dynamic changes in pore pressure, reservoir temperature, fluid composition, and interactions between formation fluids, CO 2 , and reservoir rock. Thus, multiple coupled numerical simulation models must be developed and used to more precisely understand what CO 2 storage mechanisms are most significant, as well as the long-term fate of the stored CO 2 . Our study used results from hydrodynamic, coupled hydro-geomechanical, coupled hydro-geochemical, and coupled hydro-geomechanical-geochemical models to examine how changes in geomechanical and geochemical properties can impact the injectivity or storage capacity of CO 2 . Models simulated historical field operations and then forward-modeled a water-alternate-gas (WAG) operation for 20 years, followed by a 1000-year post-injection monitoring. The work demonstrates that in this specific reservoir, geomechanical impacts appear to be more significant than any geochemical processes, which showed minimal impact on reservoir properties. Furthermore, these simulation scenarios corroborated other studies showing more than 90% of CO 2 injected are being stored and will remain stored in the reservoir. The study was able to draw on a rich set of data from this currently operational project for the work. Because the Morrow is a widespread oil-bearing reservoir in the midcontinent of the United States and is fairly typical of many sandstone reservoirs elsewhere, the study has significance beyond the boundaries of our field area.

58 GEOSCIENCES↗
NASA NTRS1989

Remote sensing and numerical modeling of suspended sediment in Laguna de Terminos, Campeche, Mexico

It is necessary to understand the complex physical processes at work in coastal lagoons in order to manage them effectively. Improved methods of data collection and analysis must be found to provide synoptic, timely hydrodynamic information because of the sheer size of some lagoons and the difficulty of acquiring in situ data (particularly in the tropics). This paper summarizes research to model salinity and suspended sediment distributions in Laguna de Terminos, Mexico, using (1) a coupled hydrodynamic and dispersion model and (2) analysis of two Landsat Thematic Mapper images collected on November 25, 1984 and April 24, 1987. Atmospherically corrected chromaticity data derived from Thermatic Mapper data were significantly correlated with modeled total suspended sediment concentrations for the two dates. Comparison between numerically modeled and remotely sensed suspended sediment maps at 1.5 x 1.5 km resolution yielded a covariation map useful for identifying areas of discrepancy between the remotely sensed data and model output.

Jensen, John R.↗
DOE OSTI2020

Parallel Implicit Hydrodynamics for High Explosive Burn Calculations

High explosives in hostile environments will require calculational capabilities that model processes, which evolve on timescales from minutes to nanoseconds. Eventually the HE will begin to move metal. To handle this temporal evolution an implicit hydrodynamics coupled to the chemical release of the HE energy is required. In addition, the use of chemical kinetics to model the transition from, the initially, slow heating of a confined high explosive through to deflagration and on to detonation requires many computational zones to model high explosive engineering systems. This requirement means that a fully parallel implicit hydrodynamics is essential. In this paper we present the calculation of a nonlinear matrix equation for the advanced particle pressure that has been made parallel and implemented in our AMR code, BABBO. This new parallel implicit hydrodynamics has been applied to a cookoff problem, as well as, one and two dimensional shock problems. Results are presented and discussed.

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

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗
NASA NTRS2022

Modeling the Influence of Moulin Shape on Subglacial Hydrology

Subglacial models represent moulins as cylinders or cones, but field observations suggest that the upper part of moulins in the Greenland Ice Sheet has more complex shapes. These more complex shapes should cause englacial water storage within moulins to vary as a function of depth, a relationship not currently accounted for in models. Here, we use a coupled englacial-subglacial channel model to explore how moulin shape affects depth-dependent moulin water storage and water pressure dynamics within a subglacial channel. We simulate seven different moulin shapes across a range of moulin sizes. We find that the englacial storage capacity at the water level is the main control over the daily water level oscillation range and that depth-varying changes in englacial water storage control the temporal shape of this oscillation. Further, the cross-sectional area of the moulin within the daily oscillation range, but not above or below this range, controls pressures within the connected subglacial channel. Specifically, large cross-sectional areas can dampen daily to weekly oscillations that occur in the surface meltwater supply. Our findings suggest that further knowledge of the shape of moulins around the equilibrium water level would improve englacial storage parameterization in subglacial hydrological models and aid predictions of hydrodynamic coupling.

Moulin↗
DOE OSTI2023

Electrochemical Flow Reactor Design Allows Tunable Mass Transport Conditions for Operando Surface Enhanced Infrared Absorption Spectroscopy

Abstract In situ attenuated total reflection surface enhanced infrared absorption spectroscopy (ATR‐SEIRAS) is often used to investigate the near‐surface electrocatalytic reaction environment. However, there is a gap in directly correlating the near‐surface reaction environment with electrocatalytic reaction rates. To that end, we designed an electrochemical flow reactor for operando electrochemical ATR‐SEIRAS and demonstrate its capability with the CO 2 reduction reaction (CO 2 RR). Roughened gold catalyst thin films are prepared on ATR silicon crystals as a model system to probe local species under CO 2 RR conditions in 0.1 M KHCO 3 . We measured changes in the interfacial CO 2 concentration as a function of applied potential and electrolyte flow rate in operando , allowing us to correlate the changes in reaction rates with the observed CO 2 concentration. Including the choice of the catalyst and electrolyte, coupling hydrodynamic control with ATR‐SEIRAS in this platform enables investigations of how the local microenvironment affects the activity and selectivity of electrochemical reactions.

Avilés Acosta, Jaime E.↗
DOE OSTI2020

Nature of Non-Adiabatic Electron-Ion Forces in Liquid Metals

We report that an accurate description of electron-ion interactions in materials is crucial for our understanding of their equilibrium and non-equilibrium properties. Here, we assess the properties of frictional forces experienced by ions in non-crystalline metallic systems, including liquid metals and warm dense plasmas, that arise from electronic excitations driven by the nuclear motion due to the presence of a continuum of low-lying electronic states. To this end, we perform detailed ab-initio calculations of the full friction tensor that characterizes the set of friction forces. The nonadiabatic electron-ion interactions introduce hydrodynamic couplings between the ionic degrees of freedom, which are sizeable between nearest neighbors. The friction tensor is generally inhomogeneous, anisotropic and non-diagonal, especially at lower densities.

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

Laser material interactions in tamped materials on picosecond time scales in aluminum

Here, a 100 ps laser is used to probe the pressure generation, depth of the non-solid ablator, and the non-linear optical effects through tamper materials. Samples consisted of an aluminum ablator with tampers of sapphire and coverslip glass. In general, the sapphire tamped sample achieves higher pressures at lower laser intensities as compared to the coverslip glass tamped sample. Attempts to model the details of this set of experimental data with standard available radiation coupled hydrodynamic codes make clear that more physics is needed in these simulations to accurately predict the impact of the tamper material on the pressure generation and the depth of non-solid aluminum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Characterization of laser-accelerated proton beams from a 0.5 kJ sub-picosecond laser for radiography applications

Laser-accelerated ion beams show promise for many applications, including high-resolution flash imaging of static or dynamic objects in next-generation radiography to probe materials and plasmas in extreme environments and inertial confinement fusion. To scale up ion beam production for radiography applications, we conducted experiments using sub-picosecond lasers up to 0.5 kJ at the OMEGA-EP facility to characterize proton beams from solid targets, primarily CH/CD sub-micron thin films from which ion beams were also used for static and dynamic radiography for the first time. For standalone sub-micron thin CH films, the highest detected proton energy is in the range of 72–97 MeV. Proton beams with highest energy near or above 60 MeV at full laser energy and similar beam profiles are also measured from low-density CD foams and flat CH foil target of micrometer-scale thickness. The ~ 700–800 nm CH/CD foils achieve the highest ion yield among the targets tested. For sub-micron thin films, the laser prepulse can expand the target and lead to complex interactions, which is simulated using coupled hydrodynamic and two-step kinetic models. Simulations suggest the presence of a micrometer-scale preplasma plateau with near-critical density and further indicate that target normal sheath acceleration, electron heating from Relativistic transparency in the preplasma plateau, and background proton reflection from carbon ion front at the rear side contribute to the resulting proton spectrum from these sub-micron thin targets at various stages. These proton beams show strong potential for radiography and for production of secondary sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A High Resolution Simulation of a Single Shock-Accelerated Particle

We report particle drag models, which capture macroviscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a two-dimensional axisymmetric Navier–Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted using the hydrocode FLAG. FLAGs capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated—a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAGs Navier–Stokes and heat diffusion solutions and to provide a rationale for recent experimental particle drag measurements.

42 ENGINEERING↗
DOE OSTI2023

Nanomotor-enhanced transport of passive Brownian particles in porous media

Artificial micro/nanomotors are expected to perform tasks in interface-rich and species-rich environments for biomedical and environmental applications. In these highly confined and interconnected pore spaces, active species may influence the motion of coexisting passive participants in unexpected ways. Using three-dimensional super-resolution single-nanoparticle tracking, we observed enhanced motion of passive nanoparticles due to the presence of dilute well-separated nanomotors in an interconnected pore space. This enhancement acted at distances that are large compared to the sizes of the particles and cavities, in contrast with the insignificant effect on the passive particles with the same dilute concentration of nanomotors in an unconfined liquid. Experiments and simulations suggested an amplification of hydrodynamic coupling between self-propelled and passive nanoparticles in the interconnected confined environment, which enhanced the effective energy for passive particles to escape cavities through small holes. This finding represents an emergent behavior of confined nanomotors and suggests new strategies for the development of antifouling membranes and drug delivery systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Practical Effects of Integrating Temperature with Strang Split Reactions

For astrophysical reacting flows, operator splitting is commonly used to couple hydrodynamics and reactions. Each process operates independent of one another, but by staggering the updates in a symmetric fashion (via Strang splitting) second order accuracy in time can be achieved. However, approximations are often made to the reacting system, including the choice of whether or not to integrate temperature with the species. Here we demonstrate through a simple convergence test that integrating an energy equation together with reactions achieves the best convergence when modeling reactive flows with Strang splitting. Additionally, second order convergence cannot be achieved without integrating an energy or temperature equation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Ionospheric D and E regions

Some areas of D and E region research are examined, with particular reference to the photochemistry and transport of minor neutral constitutents in the mesosphere and lower thermosphere. The review shows that the D and E regions become more complex physical and chemical systems as research progresses, so that some previous views on these regions have to be revised. In particular, the lower ionosphere appears to be hydrodynamically coupled to the lower atmosphere and electrodynamically coupled to the magnetosphere.

Sechrist, C. F., Jr.↗
NASA NTRS1978

Modulation of short waves by long waves

Wave-tank experiments were performed to investigate the cyclic short-wave energy changes, related in phase to an underlying long wave, which occur during active generation of the short-wave field by wind. Measurements of time series of the short-wave slope were made by a laser-optical system, where the basic long-wave parameters were controlled and wind speeds were accurately reproducible. The short-wave slope variances were found to exhibit cyclic variations that are related to the phase of the long wave. The variations result from two combined effects: (1) the short wave frequency is varied by the long-wave orbital velocity; (2) the energy of the short waves is modulated by the actions of aerodynamic and hydrodynamic couplings that operate on the short waves in a manner related to the long-wave phase.

Reece, A. M., Jr.↗
NASA NTRS1980

The Eddington limit and supercritical accretion. I - Time-independent calculations

Spherically symmetric, steady state accretion of an ionized hydrogen plasma onto a neutron star is considered for accretion rates which exceed a critical rate at which the Eddington limiting luminosity is produced. The coupled hydrodynamic and frequency integrated, radiative transfer equations are solved for accretion rates up to 10 times the nominal limit. Steady state solutions are presented that imply a multiplicity of different luminosity solutions for a single accretion rate in this 'supercritical' regime.

Burger, H. L.↗