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 379 records · Page 21

The 1997-1999 Abrupt Change of the Upper Ocean Temperature in the North Central Pacific

The abrupt warming of the north central Pacific Ocean from 1997 to 1999 is studied using an ocean data assimilation product. During this period, the average mixed-layer temperature in the region of 170-210(deg)E, 25-40(deg)N rises by 1.8 K. The major contributors to the warming are surface heat flux (1.3 K), geostrophic advection (0.7 K), and entrainment (0.7 K). For the geostrophic advection, the contributions by the zonal, meridional, and vertical components are 0.4, -0.1 and 0.3 K, respectively. Mixing and meridional Ekman advection have cooling effect. The significance of the geostrophic advection indicates the importance of ocean dynamics in controlling the abrupt warming tendency during the 1997-99 period and the inadequacy of a slab-mixed-layer model in simulating such warming tendency.

climate↗

Seasonality of the Migrating Semidiurnal Tide in the Tropical Upper Mesosphere and Lower Thermosphere and its Thermodynamic and Momentum Budget

This work uses the Specified Dynamics-Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD-WACCM-X) to determine and explain the seasonality of the migrating semidiurnal tide (SW2) components of tropical upper mesosphere and lower thermosphere (UMLT) temperature, zonal wind, and meridional wind. This work also quantifies aliasing due to SW2 in satellite-based tidal estimates. Results show that during equinox seasons, the vertical profiles of tropical UMLT temperature SW2 and zonal-wind SW2’s amplitudes have a double-peak structure while they, along with meridional-wind SW2, have a single-peak structure in their amplitudes in June solstice. Hough mode reconstruction reveals that a linear combination of five SW2 Hough modes cannot fully reproduce these features. Tendency analysis reveals that for temperature, the adiabatic term, nonlinear advection term, and linear advection term are important. For the winds, the classical terms, nonlinear advection term, linear advection term, and gravity wave drag are important. Results of our alias analysis then indicate that SW2 can induce an ∼60% alias in zonal-mean and DW1 components calculated from sampling like that of the Thermosphere–Ionosphere–Mesosphere Energetics and Dynamics satellite and the Aura satellite. This work concludes that in situ generation by wave–wave interaction and/or by gravity waves plays significant roles in the seasonality of tropical UMLT temperature SW2, zonal-wind SW2, and meridional-wind SW2. The alias analysis further adds that one cannot simply assume that SW2 in the tropical UMLT is negligible.

Tides↗

Assessment of Oceanographic Conditions during the North Atlantic EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) Field Campaign

This manuscript presents an overview of NASA’s EXport Processes in the Ocean from Remote Sensing 2021 Field Campaign in the North Atlantic (EXPORTS NA) and provides quantitative and dynamical descriptions of the physical processes modulating water transformations during the study. A major programmatic goal was to conduct the sampling in a Lagrangian mode so that ocean ecological and biogeochemical changes can be observed independent from physical advective processes. To accomplish this goal, EXPORTS NA conducted a multi-ship, multi-asset field sampling program within a retentive, anticyclonic mode water eddy. Beneath depths of ∼100 m, Lagrangian sampling assets remained within the eddy core waters (ECWs) throughout the experiment, demonstrating that the ECWs within the mode water eddy were retentive. However, strong westerly winds from four storm events deepened the mixed layer (ML) of the surface core waters (SCWs) above the eddy’s mode water core by 25-40 m and exchanged some of the SCWs with surface waters outside of the eddy via Ekman transport. Estimates of flushing times ranged from 5-8 days, with surface exchange fractions ranging from 20-75% and were consistent with particle tracking advected by combined geostrophic and Ekman velocities. The relative contributions of horizontal and vertical advection on changes in ECW tracers depended on the horizontal and vertical gradients of that tracer. For example, horizontal advection played a large role in ECW salinity fluxes, while vertical entrainment played a larger role in the fluxes of nutrients into SCW ML. Each storm injected nutrients and low oxygen waters into the ML, after which the surface ocean ecosystem responded by reducing nutrient concentrations and increasing %O2 saturation levels. Overall, ECW values of chlorophyll and POC were the largest at the onset of the field program and decreased throughout the campaign. The analysis presented provides a physical oceanographic context for the many measurements made during the EXPORTS NA field campaign while illustrating the many challenges of conducting a production-flux experiment, even in a Lagrangian frame, and the inherent uncertainties of interpreting biological carbon pump observations that were collected in a Eulerian frame of reference.

oceanography↗

A hybrid incremental projection method for thermal-hydraulics applications

In this paper, a new second-order accurate, hybrid, incremental projection method for time-dependent incompressible viscous flow is introduced in this paper. The hybrid finite-element/finite-volume discretization circumvents the well-known Ladyzhenskaya–Babuška–Brezzi conditions for stability, and does not require special treatment to filter pressure modes by either Rhie–Chow interpolation or by using a Petrov–Galerkin finite element formulation. The use of a co-velocity with a high-resolution advection method and a linearly consistent edge-based treatment of viscous/diffusive terms yields a robust algorithm for a broad spectrum of incompressible flows. The high-resolution advection method is shown to deliver second-order spatial convergence on mixed element topology meshes, and the implicit advective treatment significantly increases the stable time-step size. The algorithm is robust and extensible, permitting the incorporation of features such as porous media flow, RANS and LES turbulence models, and semi-/fully-implicit time stepping. A series of verification and validation problems are used to illustrate the convergence properties of the algorithm. The temporal stability properties are demonstrated on a range of problems with 2 ≤ CFL ≤ 100. The new flow solver is built using the Hydra multiphysics toolkit. The Hydra toolkit is written in C++ and provides a rich suite of extensible and fully-parallel components that permit rapid application development, supports multiple discretization techniques, provides I/O interfaces, dynamic run-time load balancing and data migration, and interfaces to scalable popular linear solvers, e.g., in open-source packages such as HYPRE, PETSc, and Trilinos.

97 MATHEMATICS AND COMPUTING↗

Regional-scale analysis of karst underground flow deduced from tracing experiments: examples from carbonate aquifers in Malaga province, southern Spain

Tracer concentration data from field experiments conducted in several carbonate aquifers (Malaga province, southern Spain) were analyzed following a dual approach based on the graphical evaluation method (GEM) and solute transport modeling to decipher flow mechanisms in karst systems at regional scale. The results show that conduit system geometry and flow conditions are the principal factors influencing tracer migration through the examined karst flow routes. Solute transport is mainly controlled by longitudinal advection and dispersion throughout the conduit length, but also by flow partitioning between mobile and immobile fluid phases, while the matrix diffusion process appears to be less relevant. The simulation of tracer breakthrough curves (BTCs) suggests that diffuse and concentrated flow through the unsaturated zone can have equivalent transport properties under extreme recharge, with high flow velocities and efficient mixing due to the high hydraulic gradients generated. Tracer mobilization within the saturated zone under low flow conditions mainly depends on the hydrodynamics (rather than on the karst conduit development), which promote a lower longitudinal advection and retardation in the tracer migration, resulting in a marked tailing effect of BTCs. The analytical advection-dispersion equation better approximates the effective flow velocity and longitudinal dispersion estimations provided by the GEM, while the non-equilibrium transport model achieves a better adjustment of most asymmetric and long-tailed BTCs. The assessment of karst underground flow properties from tracing tests at regional scale can aid design of groundwater management and protection strategies, particularly in large hydrogeological systems (i.e. transboundary carbonate aquifers) and/or in poorly investigated ones.

58 GEOSCIENCES↗

ENSO regimes and the late 1970's climate shift: The role of synoptic weather and South Pacific ocean spiciness

South Pacific subtropical density compensated temperature and salinity (spiciness) anomalies are known to be associated with decadal equatorial variability, however, the mechanisms by which such disturbances are generated, advect and the degree to which they modulate the equatorial thermocline remains controversial. During the late 1970's a climate regime transition preceded a period of strong and sustained El Nino events. Using an ocean general circulation model forced by the constituent mechanical and thermodynamic components of the reanalysed atmosphere we show that the late 1970's transition coincided with the arrival of a large-scale, subsurface cold and fresh water anomaly in the central tropical Pacific. An ocean reanalysis for the period 1990–2007 that assimilates subsurface Argo, XBT and CTD data, reveals that disturbances occur due to the subduction of negative surface salinity anomalies from near 30° S, 100° W which are advected along the σ=25–26 kgm{sup −3} isopycnal surfaces. These anomalies take, on average, seven years to reach the central equatorial Pacific where they may substantially perturb the thermocline before the remnants ultimately ventilate in the region of the western Pacific warm pool. Positive (warm–salty) disturbances, known to occur due to late winter diapycnal mixing and isopycnal outcropping, arise due to both subduction of subtropical mode waters and subsurface injection. On reaching the equatorial band (10° S–0° S) these disturbances tend to deepen the thermocline reducing the model's ENSO. In contrast the emergence of negative (cold–fresh) disturbances at the equator are associated with a shoaling of the thermocline and El Nino events. Process studies are used to show that the generation and advection of anomalous density compensated thermocline disturbances critically depend on stochastic forcing of the intrinsic ocean by weather. We further show that in the absence of the inter-annual component of the atmosphere forcing Central Pacific El Nino events are manifest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An unstructured-mesh finite-volume MPDATA for compressible atmospheric dynamics

An advancement of the unstructured-mesh finite-volume MPDATA (Multidimensional Positive Definite Advection Transport Algorithm) is presented that formulates the error-compensative pseudo-velocity of the scheme to rely only on face-normal advective fluxes to the dual cells, in contrast to the full vector employed in previous implementations. This is essentially achieved by expressing the temporal truncation error underlying the pseudo-velocity in a form consistent with the flux-divergence of the governing conservation law. The development is especially important for integrating fluid dynamics equations on non-rectilinear meshes whenever face-normal advective mass fluxes are employed for transport compatible with mass continuity—the latter being essential for flux-form schemes. In particular, the proposed formulation enables large-time-step semi-implicit finite-volume integration of the compressible Euler equations using MPDATA on arbitrary hybrid computational meshes. Furthermore, it facilitates multiple error-compensative iterations of the finite-volume MPDATA and improved overall accuracy. The advancement combines straightforwardly with earlier developments, such as the nonoscillatory option, the infinite-gauge variant, and moving curvilinear meshes. A comprehensive description of the scheme is provided for a hybrid horizontally-unstructured vertically-structured computational mesh for efficient global atmospheric flow modelling. The proposed finite-volume MPDATA is verified using selected 3D global atmospheric benchmark simulations, representative of hydrostatic and non-hydrostatic flow regimes. Besides the added capabilities, the scheme retains fully the efficacy of established finite-volume MPDATA formulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Preserving Tracer Correlations in Moment-Based Atmospheric Transport Models

A linear non-diffusive algorithm for advective transport is developed that greatly improves the detail at which aerosols and clouds can be represented in atmospheric models. Linear advection schemes preserve tracer correlations but the most basic linear scheme is rarely used by atmospheric modelers on account of its excessive numerical diffusion. Higher-order schemes are in widespread use, but these present new problems as nonlinear adjustments are required to avoid occurrences of negative concentrations, spurious oscillations, and other non-physical effects. Generally successful at reducing numerical diffusion during the advection of individual tracers, for example, particle number or mass, the higher-order schemes fail to preserve even the simplest of correlations between interrelated tracers. As a result, important attributes of aerosol and cloud populations including radial moments of particle size distributions, molecular precursors related through chemical equilibria, aerosol mixing state, and distribution of cloud phase are poorly represented. We introduce a new transport scheme, minVAR, that is both non-diffusive and preservative of tracer correlations, thereby combining the best features of the basic and higher-order schemes while enabling new features such as the tracking of sub-grid information at arbitrarily fine scales with high computational efficiency.

54 ENVIRONMENTAL SCIENCES↗

Evolution of accretion disks in tidal disruption events

During a stellar tidal disruption event (TDE), an accretion disk forms as stellar debris returns to the disruption site and circularizes. Rather than being confined within the circularizing radius, the disk can spread to larger radii to conserve angular momentum. A spreading disk is a source of matter for re-accretion at rates that may exceed the later stellar fallback rate, although a disk wind can suppress its contribution to the central black hole accretion rate. A spreading disk is detectible through a break in the central accretion rate history or, at longer wavelengths, by its own emission. We model the evolution of TDE disk size and accretion rate by accounting for the time-dependent fallback rate, for the influence of wind losses in the early advective stage, and for the possibility of thermal instability for accretion rates intermediate between the advection-dominated and gas-pressure-dominated states. The model provides a dynamic basis for modeling TDE light curves. All or part of a young TDE disk will precess as a solid body because of the Lense-Thirring effect, and precession may manifest itself as a quasi-periodic modulation of the light curve. The precession period increases with time. Applying our results to the jetted TDE candidate Swift J1644+57, whose X-ray light curve shows numerous quasi-periodic dips, we argue that the data best fit a scenario in which a main-sequence star was fully disrupted by an intermediate mass black hole on an orbit significantly inclined from the black hole equator, with the apparent jet shutoff at t = 500 days corresponding to a disk transition from the advective state to the gas-pressure-dominated state.

79 ASTRONOMY AND ASTROPHYSICS↗

THE COSMIC BATTERY IN ASTROPHYSICAL ACCRETION DISKS

The aberrated radiation pressure at the inner edge of the accretion disk around an astrophysical black hole imparts a relative azimuthal velocity on the electrons with respect to the ions which gives rise to a ring electric current that generates large-scale poloidal magnetic field loops. This is the Cosmic Battery established by Contopoulos and Kazanas in 1998. In the present work we perform realistic numerical simulations of this important astrophysical mechanism in advection-dominated accretion flows, ADAFs. We confirm the original prediction that the inner parts of the loops are continuously advected toward the central black hole and contribute to the growth of the large-scale magnetic field, whereas the outer parts of the loops are continuously diffusing outward through the turbulent accretion flow. This process of inward advection of the axial field and outward diffusion of the return field proceeds all the way to equipartition, thus generating astrophysically significant magnetic fields on astrophysically relevant timescales. We confirm that there exists a critical value of the magnetic Prandtl number between unity and 10 in the outer disk above which the Cosmic Battery mechanism is suppressed.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical investigation of microbial quorum sensing under various flow conditions

Microorganisms efficiently coordinate phenotype expressions through a decision-making process known as quorum sensing (QS). We investigated QS amongst distinct, spatially distributed microbial aggregates under various flow conditions using a process-driven numerical model. Model simulations assess the conditions suitable for QS induction and quantify the importance of advective transport of signaling molecules. In addition, advection dilutes signaling molecules so that faster flow conditions require higher microbial densities, faster signal production rates, or higher sensitivities to signaling molecules to induce QS. However, autoinduction of signal production can substantially increase the transport distance of signaling molecules in both upstream and downstream directions. We present empirical approximations to the solutions of the advection–diffusion–reaction equation that describe the concentration profiles of signaling molecules for a wide range of flow and reaction rates. These empirical relationships, which predict the distribution of dissolved solutes along pore channels, allow to quantitatively estimate the effective communication distances amongst multiple microbial aggregates without further numerical simulations.

59 BASIC BIOLOGICAL SCIENCES↗

An Evaluation of Size-Resolved Cloud Microphysics Scheme Numerics for use with Radar Observations Part II: Condensation and Evaporation

Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part 1 addressed collision-coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision-coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.

Cloud microphysics↗

Linking Soil Moisture Variation and Abundance of Plants to Geomorphic Processes: A Generalized Model for Erosion-Uplifting Landscapes

American Geophysical Union. All Rights Reserved. The diffusive and advective erosion-created landscapes have similar structure (hillslopes and channels) across different scales regardless of variations in drivers and controls. The relative magnitude of diffusive erosion to advective erosion ( D/K ratio) in a landscape development model controls hillslope length, shape, and drainage density, which regulate soil moisture variation, one of the critical resources of plants, through the contributing area ( A ) and local slope ( S ) represented by a topographic index (TI). Here we explore the theoretical relation between geomorphic processes, TI, and the abundance and distribution of plants. We derived an analytical model that expresses the TI with D, K , and A . This gives us the relation between soil moisture variation and geomorphic processes. Plant tolerance curves are used to link plant performance to soil moisture. Using the hypothetical tolerance curves of three plants, we show that the abundance and distribution of xeric, mesic, and hydric plants on the landscape are regulated by the D/K ratio. Where diffusive erosion is the major erosion process (large D/K ratio), mesic plants have higher abundance relative to xeric and hydric plants and the landscape has longer and convex-upward hillslope and low channel density. Finally, increasing the dominance of advective erosion increases relative abundance of xeric and hydric plants dominance, and the landscape has short and concave hillslope and high channel density.

59 BASIC BIOLOGICAL SCIENCES↗

Effective Stochastic Model For Reactive Transport

We propose an effective stochastic advection-diffusion-reaction (SADR) model. Unlike traditional advection-dispersion-reaction models, the SADR model describes mechanical and diffusive mixing as two separate processes. In the SADR model, the mechanical mixing is driven by random advective velocity with the variance given by the coefficient of mechanical dispersion. The diffusive mixing is modeled as a fickian diffusion with the effective diffusion coefficient. Both coefficients are given in terms of Peclet number (Pe) and the coefficient of molecular diffusion. We use the experimental results of \cite{Gramling} to demonstrate that for transport and bimolecular reactions in porous media the SADR model is significantly more accurate than the traditional dispersion model, which overestimates the mass of the reaction product by as much as 25\%.

stochastic equations, mixing controlled reactions,↗

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗

Reactive fluid flow in CO 2 storage reservoirs: A 2‐D pore network model study

Abstract Water acidifies in the presence of CO 2 and prompts mineral dissolution. A 2‐D pore network model scheme is developed to investigate reactive fluid flow in CO 2 storage reservoirs during injection when advective transport prevails. Mineral dissolution satisfies kinetic rate laws and continues until thermodynamic equilibrium is reached. In advection‐dominant regimes, network simulation results show that species concentration, tube enlargement and flow rate can be summarized in terms of the dimensionless Damköhler number Da which is the ratio between advection time along a pore and the reaction time. Reservoirs will tend to experience localized enlargement near injection wells (before water drying) and compact dissolution in the far‐field. The Damköhler number couples with initial pore‐size variability to distort the relationship between mean tube diameter and either local or network‐average flow rates. Both the Damköhler number and pore‐size variability should be considered in field‐scale numerical simulators.

Kim, Seunghee↗

Understanding the vertical structure of potential vorticity in tropical depressions

Here, potential vorticity (PV) has been used to understand the intensification and motion of a variety of tropical vortices. In this work, atmospheric reanalyses and idealized models are used to understand how the vertical structures of moist convective heating and adiabatic advection jointly shape the vertical structures of PV in tropical depressions. Observationally based estimates reveal a top-heavy PV structure in tropical depressions, contrasting with bottom-heavy structures of absolute vorticity and diabatic PV generation. These distinct vertical structures are reproduced in an axisymmetric model which employs the weak temperature gradient approximation for conceptual simplicity and is forced by stratiform and deep convective heating. When provided in isolation, the stratiform and deep convective heatings produce PV maxima at 500 hPa and near the surface, respectively. When these two heatings are applied simultaneously, interactions between the stratiform and deep convective modes enhance the adiabatic advective tendencies produced by the transverse circulation, making the PV distribution more top-heavy. In the lower and middle troposphere, radial advection also greatly reduces the radius of the PV structure relative to that of the imposed heating, consistent with structures in observed tropical depressions; the implications of these differences in radial structures for using the flux form of the relevant conservation equations (e.g. for PV substance or absolute vorticity) are discussed.

54 ENVIRONMENTAL SCIENCES↗

Explaining persistent incomplete mixing in multicomponent reactive transport with Eulerian stochastic model

We present an Eulerian stochastic advection–diffusion–reaction (SADR) model and use it to explain incomplete mixing typically observed in transport experiments with bimolecular reactions. Unlike traditional advection– dispersion–reaction (ADR) models, the SADR model describes mechanical and diffusive mixing as two separate processes. In the SADR model, mechanical mixing is driven by random advective velocity whose variance is given by the coefficient of mechanical dispersion. The diffusive mixing is modeled as a Fickian diffusion process with the effective diffusion coefficient. We demonstrate that the sum of the two coefficients is equal to the dispersion coefficient, but only the effective diffusion coefficient contributes to the mixing-controlled reactions. We use experimental results of Gramling et al. * * to show that for transport and bimolecular reactions in porous media, the SADR model is significantly more accurate than the ADR model, which overestimates the concentration of the reaction product by as much as 60%. We also show that the SADR model predicts an experimentally observed bell- shaped spatial distribution of the reactive product concentration, while the ADR model results in a concentration distribution with an unphysical kink.

Bimolecular reactions, incomplete mixing, Stochast↗