Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ADVECTION”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Hydrothermal discharge from the El Tatio basin, Atacama, Chile

El Tatio in northern Chile is one of the best-studied geothermal fields in South America. However, there remain open questions about the mass and energy budgets, water recharge rates and residence time in the subsurface, origin of dissolved solutes, and processes affecting the phase and chemical composition of groundwater and surface water. Here we measured and sampled surface manifestations of the geothermal system (geysers, perpetual spouters, mud pools/volcanoes, and non-eruptive hot springs) and meteoric water. From the isotopic composition we infer that the thermal water has a meteoric origin that is different from the composition of local meteoric water. The absence of detectable tritium in thermal waters indicates that most of the recharge occurred pre-1950. Boiling and steam separation from the deep reservoir appear to be the main subsurface processes affecting the thermal fluids. A large amount of heat is lost to the atmosphere by evaporation from surface water and by steam emitted from erupting geysers. Finally, using the chloride inventory method, we estimate thermal water discharge to be 218 to 234 L/s, and the advective heat flow to be 120 to 170 MW.

15 GEOTHERMAL ENERGY↗
DOE OSTI2017

Benthic storms, nepheloid layers, and linkage with upper ocean dynamics in the western North Atlantic

Benthic storms are episodic periods of strong abyssal currents and intense, benthic nepheloid (turbid) layer development. In order to interpret the driving forces that create and sustain these storms, we synthesize measurements of deep ocean currents, nephelometer-based particulate matter (PM) concentrations, and seafloor time-series photographs collected during several science programs that spanned two decades in the western North Atlantic. Benthic storms occurred in areas with high sea-surface eddy kinetic energy, and they most frequently occurred beneath the meandering Gulf Stream or its associated rings, which generate deep cyclones, anticyclones, and/or topographic waves; these create currents with sufficient bed-shear stress to erode and resuspend sediment, thus initiating or enhancing benthic storms. Occasionally, strong currents do not correspond with large increases in PM concentrations, suggesting that easily erodible sediment was previously swept away. Periods of moderate to low currents associated with high PM concentrations are also observed; these are interpreted as advection of PM delivered as storm tails from distal storm events. Outside of areas with high surface and deep eddy kinetic energy, benthic nepheloid layers are weak to non-existent, indicating that benthic storms are necessary to create and maintain strong nepheloid layers. Origins and intensities of benthic storms are best identified using a combination of time-series measurements of bottom currents, PM concentration, and bottom photographs, and these should be coupled with water-column and surface-circulation data to better interpret the specific relations between shallow and deep circulation patterns. Understanding the generation of benthic nepheloid layers is necessary in order to properly interpret PM distribution and its influence on global biogeochemistry.

58 GEOSCIENCES↗
DOE OSTI2015

Permeability and porosity of hydrate-bearing sediments in the northern Gulf of Mexico

Hydrate-bearing sands are being actively explored because they contain the highest concentrations of hydrate and are the most economically recoverable hydrate resource. However, relatively little is known about the mechanisms or timescales of hydrate formation, which are related to methane supply, fluid flux, and host sediment properties such as permeability. We used logging-while-drilling data from locations in the northern Gulf of Mexico to develop an effective medium theory-based model for predicting permeability based on clay-sized sediment fraction. The model considers permeability varying between sand and clay endpoint permeabilities that are defined from laboratory data. We verified the model using permeability measurements on core samples from three boreholes, and then used the model to predict permeability in two wells drilled in Walker Ridge Block 313 during the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II expedition in 2009. We found that the cleanest sands (clay-sized fraction <0.05) had intrinsic (hydrate-free) permeability contrasts of 5-6 orders of magnitude with the surrounding clays, which is sufficient to provide focused hydrate formation due to advection of methane from a deep source or diffusion of microbial methane from nearby clay layers. In sands where the clay-sized fraction exceeds 0.05, the permeability reduces significantly and focused flow is less pronounced. In these cases, diffusion of dissolved microbial methane is most likely the preferred mode of methane supply for hydrate formation. In conclusion, our results provide important constraints on methane supply mechanisms in the Walker Ridge area and have global implications for evaluating rates of methane migration and hydrate formation in hydrate-bearing sands.

58 GEOSCIENCES↗
DOE OSTI2019

FVCOM-plume – A Three-dimensional Lagrangian Outfall Plume Dilution and Transport Model for Dynamic Tidal Environments: Model Development

Although numerous steady-state effluent plume dilution models are in use, their application in tidal environments remains a challenge. Three-dimensional dynamic circulation models are also inadequate, often due to the lack of required resolution and simplifying assumptions. To overcome these limitations, we present FVCOM-plume—an outfall plume dilution and transport model that operates within the Finite Volume Coastal Ocean Model (FVCOM) framework. It provides simultaneous inclusion of near-field dilution and far-field plume transport processes. The near-field is based on UM3 model using Lagrangian Control Volume approach to compute buoyant plume trajectory and dilution from multiport diffusers. The far-field uses neutrally buoyant particles with point masses and the random walk method to solve unsteady advection-diffusion processes. A density kernel approach is used to compute concentrations at point locations and analyze far-field plume characteristics. The results demonstrate the ability of FVCOM-plume to simultaneously capture near-field and far-field effluent plume dynamics in tidal environments.

Plume Dilution Models↗
DOE OSTI2025

Nonlinear simulation of under-resolved flows with shocks

Here, we consider the numerical simulation of advection-dominated flows whose wide range of physical length scales exceed the memory capacity of finite computers. Simulating flows with shocks and turbulence presented challenges for the earliest computers that were quickly overcome by the development of new numerical methodology. Principal among those new ideas were artificial viscosity and finite volume methods, concepts that remain in common use today. We begin by describing the history of those methods, the innovators and their motivations. We then describe the development of finite scale theory, a reformulation of Navier–Stokes theory that exposes the physical principles on which artificial viscosity is based. We discuss the essential properties of the finite scale equations, the observer, unresolved kinetic energy and inviscid energy dissipation. We briefly consider the implementation of the finite scale equations on the computer from the point of view of Gisin’s conjectures about finite information.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Electro-osmo-dialysis through nanoporous layers physically conjugated to micro-perforated ion-exchange membranes: Highly selective accumulation of trace coions

This study explores theoretically and numerically the phenomenon of counter-flow accumulation of trace coions during electro-osmo-dialysis through composite membranes containing nanoporous and perforated ion-exchange layers. The combination of charged nanopores and perforations allows electroosmosis, and ion concentration polarization appears near the ion-exchange layer. Accumulation of trace coions occurs because their advective and electromigration flow velocity components have opposite signs and are locally compensating in depletion regions where electric fields are high. This local compensation of flow components can give rise to high accumulation factors (>1000 in some cases) for trace coions at some distance from the interface between the nanoporous and ion-exchange layers. Moreover, the extent of this accumulation depends on the trace-coion diffusion coefficient. The resulting selectivity between different coions might serve as a basis for environmentally friendly separation of ionic species such as synthetic peptides or rare-earth elements. For example, a periodic process might include accumulation followed by species release into smaller volumes of solution. For both rare-earth ions and peptide-like species, the accumulation selectivity between ions with small differences in mobilities depends on system parameters such as nanoporous-layer thickness, distance between the perforations, and the zeta potential of the nanopore surface. Here, the use of optimal values of these membrane parameters may enable tuning of the electro-osmo-dialysis system for separation of specific ions. Modelling reveals that despite inhomogeneous patterns of fluid and ion flows close to the nanoporous/perforated ion-exchange interface, at relatively short distances from this interface (less than half the distance between the perforations) all the flows become practically one-dimensional (1D). Similarly, the concentration profiles of trace coions also become 1D. The flow one-dimensionality at short distances from the current-polarized interface enables a simple 1D approximation that is useful for an approximate analysis of trends in this multi-parameter system.

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

Highly selective ion separations based on counter-flow electromigration in nanoporous membranes

Highly selective ion separations are vital for recovering important salts, and membrane-based techniques may enable environmentally friendly ion separations that operate continuously. This study demonstrates that simply opposing convective flow with electromigration in track-etched membranes leads to remarkable selectivities among monovalent ions. Because the membrane-pore diameters are orders of magnitude larger than the ion diameters, advection moves all ions at the same velocity, whereas electrophoretic velocities are proportional to ion mobilities. Thus during counter-flow cation electromigration, less mobile cations have lower electromigration velocity components and pass through porous membranes more rapidly than more mobile cations. Using membranes with 400 nm pores, counter-flow electromigration gives Li + /K + selectivities of 100 and Li + /Na + selectivities around 30 in mixed-salt studies. In this work, numerical simulations based on the extended Nernst-Planck equation agree with trends in experimental data and highlight the importance of high Péclet numbers, high current-to-flow ratios, and uniform current distributions for achieving high selectivities. Importantly, experimental separations give high selectivities at 0.3 M ionic strength and low Li + /K + ratios, which suggest these methods may also work in brine solutions. Nevertheless, energy costs for such separations are high.

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

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗
DOE OSTI2020

Enhanced detection of volatile organic compounds (VOCs) by caffeine modified carbon nanotube junctions

Detection of vapor phase analytes and gas molecules is primarily essential for numerous applications including medical diagnostic, environmental monitoring, chemical threat detection, food quality control, and chemical warfare agents. Here, we developed Caffeine functionalized carbon nanotubes (CNTs) random network based vapor quantum resistive sensors (vQRSs) for volatile organic compounds (VOCs) detection. The localized transmission electron microscopy (TEM) image and energy dispersive X-ray spectroscopy (EDS) confirm migration of caffeine to CNT–CNT junction (CNTj). We have recorded more than one order of magnitude higher sensitivity and unique selectivity towards selective VOCs (e.g., acetone, toluene, ethanol, methanol, etc.) with short response time (few seconds) for caffeine modified CNT junction (Caf-CNTj) vQRSs. Molecular dynamics (MD) calculations suggest that the sensitivity of the sensor is diffusion limited, with an advective contribution for ethanol. Furthermore, the designed Caf-CNTj sensors can potentially be important candidate to be a part of electronic nose used for various application such as breath analysis and food quality monitoring.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Active-noise-induced dynamic clustering of passive colloidal particles

Active fluids generate spontaneous, often chaotic mesoscale flows. Harnessing these flows to drive soft materials embedded within an active fluid into structures with controlled length scales and lifetimes is a key challenge at the interface between the fields of active matter and nonequilibrium self-assembly. Here, we present a simple and efficient computational approach to model soft materials advected by active fluids, by simulating particles moving in a spatiotemporally correlated noise field. To illustrate our approach, we simulate the dynamical self-organization of repulsive colloids within such an active noise field. The colloids form structures whose sizes and dynamics can be tuned by the correlation time and length of the active fluid, and range from small rotating droplets to clusters with internal flows and system-spanning sizes that vastly exceed the active correlation length. Our results explain how the interplay between active fluid time and length scales and emergent driven assembly can be used to rationally design functional assemblies. More broadly, our approach can be used to efficiently simulate diverse active fluids and other systems with spatiotemporally correlated noise.

Brownian dynamics↗
DOE OSTI2019

Simulation of edge localized mode heat pulse using drift-kinetic ions and Boltzmann electrons

This article reports on 1D+2V heat pulse propagation studies using the COGENT guiding center kinetic code. The model uses magnetized kinetic ions and a simple Boltzmann electron model. Results agree with previous kinetic and fluid modeling benchmark studies that correspond to the parameters of edge localized modes (ELMs) observed on the JET tokamak. The plasma parameters for the edge pedestal and ensuing ELM dynamics are in the low collisionality regime. Hence, the dominant balance between the assumed Maxwellian ELM source and collisionless parallel advection causes the ion PDF to develop a significantly anisotropic velocity distribution. Adding nonlinear Coulomb ion-ion collisions to the model acts to smooth the sharp features of the ion distribution function, but the anisotropy remains robust due to the low collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Development and Implementation of a CTF Code Verification Suite.

CTF is a thermal hydraulic subchannel code developed to predict light water reactor (LWR) core behavior. It is a version of Coolant Boiling in Rod Arrays (COBRA) developed by Oak Ridge National Laboratory (ORNL) and North Carolina State University (NCSU) and used in the Consortium for the Advanced Simulation of LWRs (CASL). Here, the existing CTF code verification matrix is expanded, which ensures that the code is a faithful representation of the underlying mathematical model. The suite of code verification tests are mapped to the underlying conservation equations of CTF and significant gaps are addressed. As such, five new problems are incorporated: isokinetic advection, conduction, pressure drop, convection, and pipe boiling. Convergence behavior and numerical errors are quantified for each of the tests and all tests converge at the correct rate to their corresponding analytic solution. A new verification utility that generalizes the code verification process is used to incorporate these problems into the CTF automated test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Alteration in micro-mechanical characteristics of wellbore cement fracture surfaces due to fluid exposure

The stiffness of wellbore cement fracture surfaces was measured after exposing to the advective flow of nitrogen, silicone oil, and medium sweet dead crude oil for different exposure periods. The test specimens were extracted from fractured cement cylinders, where the cement fracture surfaces were exposed to the different fluids up to 15 weeks. A nanoindenter with a Berkovich indenter tip was used to measure load-indentation depth data, which was used to extract the elastic modulus (E) and nano-hardness (H) of the cement fracture surfaces. A reduction in the elastic modulus compared with an unexposed specimen were observed in all the specimens. Both elastic modulus and nano-hardness for the specimens exposed to silicone oil were lower than specimens exposed to nitrogen gas and varied with the period of exposure. The elastic modulus and nano-hardness of the specimens exposed to crude oil were the lowest with a significant decrement with the exposure period. The frequency distribution of the nanoindentation measurements shows that the volume-fraction ratio of the two types of cement hydrated nanocomposites for both the unexposed and test specimens is about 70:30%. Phase transformation beneath the indenter is observed for all of the specimens, with more obvious plastic deformation in specimens exposed to crude oil. Analytical measurements (SEM, EDS, FT-IR, and XRD) on exposed cement fracture surfaces reveal different levels of physical and chemical alteration that are consistent with the reduction in stiffness measured by nanoindentation. Furthermore, the study suggests that cement stiffness will decrease due to crude oil exposure, and the fracture will be sensitive to stress and pore pressure with time.

42 ENGINEERING↗
DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗
DOE OSTI2020

A passive model for the evolution of subgrid-scale instabilities in turbulent flow regimes

In this work, a zero-dimensional model for the development of turbulent mixing layers driven by Rayleigh–Taylor (RT) or Richtmyer–Meshkov (RM) instabilities is presented, for application in Reynolds-Averaged Navier–Stokes (RANS) simulations. The model approximates the behavior of the turbulent mixing layer at a passively advected tracer particle that sits on the material interface of interest in the simulation, and is used to provide initial conditions for the RANS model once the mixing layer width has developed to a size that may be resolved on the computational mesh. The model relies on approximations to the Besnard–Harlow–Rauenzahn (BHR) RANS closure equations and is thus consistent the turbulent development of the mixing layer, unlike previous models of this type which have focused on the early-time laminar development of RT or RM instabilities. These instabilities will often grow through an initially laminar state before developing into turbulence, and criteria for switching from an existing laminar model to the proposed turbulent model are considered. The proposed model reproduces the mixing layer development of RT and RM instabilities seen in mixed width resolving BHR RANS simulations to a reasonable degree of accuracy, and relaxes the need for the RANS simulation to resolve the mixed width as soon as the mixing layer becomes turbulent, which may be impractical in cases where the mixing layer grows through a wide range of scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Latent-space time evolution of non-intrusive reduced-order models using Gaussian process emulation

Non-intrusive reduced-order models (ROMs) have recently generated considerable interest for constructing computationally efficient counterparts of nonlinear dynamical systems emerging from various domain sciences. They provide a low-dimensional emulation framework for systems that may be intrinsically high-dimensional. This is accomplished by utilizing a construction algorithm that is purely data-driven. It is no surprise, therefore, that the algorithmic advances of machine learning have led to non-intrusive ROMs with greater accuracy and computational gains. However, in bypassing the utilization of an equation-based evolution, it is often seen that the interpretability of the ROM framework suffers. This becomes more problematic when black-box deep learning methods are used which are notorious for lacking robustness outside the physical regime of the observed data. In this article, we propose the use of a novel latent-space interpolation algorithm based on Gaussian process regression. Notably, this reduced-order evolution of the system is parameterized by control parameters to allow for interpolation in space. The use of this procedure also allows for a continuous interpretation of time which allows for temporal interpolation. The latter aspect provides information, with quantified uncertainty, about full-state evolution at a finer resolution than that utilized for training the ROMs. This research assesses the viability of this algorithm for an advection-dominated system given by the inviscid shallow water equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Return flows from beaver ponds enhance floodplain-to-river metals exchange in alluvial mountain catchments

River to floodplain hydrologic connectivity is strongly enhanced by beaver- (Castor canadensis) engineered channel water diversions. The hydroecological impacts are wide ranging and generally positive, however, the hydrogeochemical characteristics of beaver-induced flowpaths have not been thoroughly examined. Using a suite of complementary ground- and drone-based heat tracing and remote sensing methodology we characterized the physical template of beaver-induced floodplain exchange for two alluvial mountain streams near Crested Butte, Colorado, USA. A flowpath-oriented perspective to water quality sampling allowed characterization of the chemical evolution of channel water diverted through floodplain beaver ponds and ultimately back to the channel in 'beaver pond return flows'. Subsurface return flow seepages were universally suboxic, while ponds and surface return flows showed a range of oxygen concentration due to in-situ photosynthesis and atmospheric mixing. Median concentrations of reduced metals: manganese (Mn), iron (Fe), aluminum (Al), and arsenic (As) were substantially higher along beaver-induced flowpaths than in geologically controlled seepages and upstream main channel locations. The areal footprint of reduced return seepage flowpaths were imaged with surface electromagnetic methods, indicating extensive zones of high-conductivity shallow groundwater flowing back toward the main channels and emerging at relatively warm bank seepage zones observed with infrared. Multiple-depth redox dynamics within one focused seepage zone showed coupled variation over time, likely driven by observed changes in seepage rate that may be controlled by pond stage. High-resolution times series of dissolved Mn and Fe collected downstream of the beaver-impacted reaches demonstrated seasonal dynamics in mixed river metal concentrations. Al time series concentrations showed proportional change to Fe at the smaller stream location, indicating chemically reduced flowpaths were sourcing Al to the channel. Overall our results indicated beaver-induced floodplain exchanges create important, and perhaps dominant, transport pathways for floodplain metals by expanding chemically-reduced zones paired with strong advective exchange.

58 GEOSCIENCES↗
DOE OSTI2023

Contrasting changes in ozone during 2019–2021 between eastern and the other regions of China attributed to anthropogenic emissions and meteorological conditions

Ozone pollution is one of the most severe air quality issues in China that poses a serious threat to human health and ecosystems. During 2019–2021, the maximum daily 8-h average ozone concentrations in eastern China (110–122.5°E, 26–42°N) and the rest of China (ROC) show different decreasing patterns, with ozone concentrations in eastern China decreasing by 14.9 μg/m 3 , which is much larger than 4.8 μg/m 3 in ROC. Here, based on two independent methods, the atmospheric chemical transport model (GEOS-Chem) simulations and the machine learning (ML) model (LightGBM) predictions, the reasons for the differences in ozone changes between eastern China and ROC during the warm season (April to September) are investigated. According to the GEOS-Chem (LightGBM) results, changes in the meteorological conditions contributed to an ozone decrease by 7.3 (6.8) μg/m 3 in eastern China due to decreased chemical production and an ozone decrease by 6.8 (7.0) μg/m 3 in ROC attributed to the weakened horizontal and vertical advection. With the influence of meteorological factors excluded, the observations show that changes in anthropogenic emissions resulted in an ozone decrease by 7.6 (8.1) μg/m 3 in eastern China and an ozone increase by 2.0 (2.2) μg/m 3 in ROC, which is primarily induced by the changes in NO x emissions. The surface measurements and satellite retrievals also indicate that the reduction in NO x emissions in ROC is less efficient than that in the more developed eastern China, leading to contrasting changes in ozone concentrations between eastern China and ROC during 2019–2021. In conclusion, our results highlight the critical need to reduce ozone precursor emissions in the rest regions of China apart from eastern China.

54 ENVIRONMENTAL SCIENCES↗