Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow and tracer simulation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

High-resolution simulation of free-surface flow and tracer retention over streambeds with ripples

This study presents a novel high-resolution simulation of free-surface flow and tracer retention over a streambed with ripples based on varying ripple morphologies, surface hydraulics and the transport of a tracer pulse from surface water to surface dead zone. For the simulations, the computational fluid dynamics (CFD) model OpenFOAM was used to solve the three-dimensional Navier-Stokes equations in combination with an implemented transport equation. Pressure gradients at the streambed were used to account for hyporheic exchange, assuming water flow from high pressure zones to low pressure zones. Flow velocities, ripple sizes and spacing showed to significantly affect these pressure gradients, but also the transport of a passive tracer at the streambed, which was not investigated so far. Due to the velocity field, large parts of the tracer mass were transported alongside the main stream above the ripples. Tracer mass reaching the space between the ripples was temporarily retained due to low velocities and recirculations. It was shown that the retention is depending on the ripple size and space between the ripples as well as on the flow velocity. Decreasing ripple sizes and higher flow velocities lead to a smaller tracer retention. Furthermore we showed that the ripple length to height ratio controls the generation of recirculation zones which affect the residence time of the tracer significantly. Ripple spacing leads to temporarily higher tracer concentration at the streambed, but smaller tracer retention. We conclude that the impact of the streambed morphology on the hydraulics in combination with tracer retention should be addressed for a comprehensive understanding of compound movement, exchange and transformation within the hyporheic zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Chemistry in CESM-SE: Evaluation, Performance and Optimization

The purpose of the proposed work focused on development of chemistry representation within the Spectral Element (SE) dynamical core as implemented in the Community Earth System Model (CESM). More specifically, a main focus was on the ability of SE to accurately represent tracer transport. The proposed approach was to incrementally increase the complexity of the problem, starting from specified two-dimensional flow and tracers to simulations using specified dynamics and full chemistry. As demonstrated below, we have successfully studied all aspects of the proposed work, although only part of the work has been published in the refereed literature so far. Furthermore, because the SE dynamical core has been found to have several deficiencies that are still being investigated for solution, not all proposed tasks were finalized. In addition to the tests for SE performance, in an effort to decrease the computational burden of interactive chemistry, especially in the case of a large number of chemical species and chemical reactions, development on a faster chemical solver and implementation on GPUs has been implemented in CESM under the leadership of John Drake (U. Tennessee).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Helium bubble flow visualization of the spanwise separation on a NACA 0012 with simulated glaze ice

Research has been performed to experimentally visualize and document the flow separation due to simulated glaze ice accretion on a NACA 0012 semispan with 30-deg sweep using helium bubbles as flow tracers. Results are compared to Navier-Stokes computational simulations for different angles of attack. Prior to acquiring data for the semispan model, a two-dimensional experiment was conducted to determine the accuracy of using the helium bubbles as flow tracers. Results from the 3D experiment compare well to the computational simulations.

Kerho, M.↗
NASA NTRS1992

Helium bubble flow visualization of the spanwise separation on a NACA 0012 with simulated glaze ice

Research was performed to experimentally visualize and document the flow separation due to simulated glaze ice accretion on a NACA 0012 semispan with 30 deg sweep using helium bubbles as flow tracers. Results are compared to Navier-Stokes computational simulations for different angles of attack. Prior to acquiring data for the semispan model, a two dimensional experiment was conducted to determine the accuracy of using the helium bubbles as flow tracers. Results from the three dimensional experiment compare well to the computational simulations.

Kerho, M.↗
DOE OSTI2024

Integrating a Water Tracer Model Into WRF‐Hydro for Characterizing the Effect of Lateral Flow in Hydrologic Simulations

Abstract Most current land models approximate terrestrial hydrological processes as one‐dimensional vertical flow, neglecting lateral water movement from ridges to valleys. Such lateral flow is fundamental at catchment scales and becomes crucial for finer‐scale land models. To test the effect of incorporating lateral flow toward three‐dimensional representations of hydrological processes in the next generation land models, we integrate a water tracer model into the WRF‐Hydro framework to track water movement from precipitation to discharge and evapotranspiration. This hydrologic‐tracer integrated system allows us to identify the key mechanisms by which lateral flow affects the flow paths and transit times in WRF‐Hydro. By comparing modeling experiments with and without lateral routing in two contrasting catchments, we determine the impacts of lateral flow on the transit times of precipitation event‐water. Results show that with limited hydrologic connectivity, lateral flow extends the transit times by reducing (increasing) event‐water drainage loss (accumulation) in ridges (valleys) and allowing reinfiltration of infiltration‐excess flow, which is missing in most land models. On the contrary with high hydrologic connectivity, lateral flow can effectively accelerate the water release to streams and reduce the transit time. However, the transit times are substantially underestimated by the model compared with isotope‐derived estimates, indicating model limitations in representing flow paths and transit times. This study provides some insights on the fundamental differences in terrestrial hydrology simulated by land models with and without lateral flow representation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Multiscale solute transport upscaling for a three‐dimensional hierarchical porous medium

Abstract A laboratory‐generated hierarchical, fully heterogeneous aquifer model (FHM) provides a reference for developing and testing an upscaling approach that integrates large‐scale connectivity mapping with flow and transport modeling. Based on the FHM, three hydrostratigraphic models (HSMs) that capture lithological (static) connectivity at different resolutions are created, each corresponding to a sedimentary hierarchy. Under increasing system ln K variances (0.1, 1.0, 4.5), flow upscaling is first conducted to calculate equivalent hydraulic conductivity for individual connectivity (or unit) of the HSMs. Given the computed flow fields, an instantaneous, conservative tracer test is simulated by all models. For the HSMs, two upscaling formulations are tested based on the advection‐dispersion equation (ADE), implementing space versus time‐dependent macrodispersivity. Comparing flow and transport predictions of the HSMs against those of the reference model, HSMs capturing connectivity at increasing resolutions are more accurate, although upscaling errors increase with system variance. Results suggest: (1) by explicitly modeling connectivity, an enhanced degree of freedom in representing dispersion can improve the ADE‐based upscaled models by capturing non‐Fickian transport of the FHM; (2) when connectivity is sufficiently resolved, the type of data conditioning used to model transport becomes less critical. Data conditioning, however, is influenced by the prediction goal; (3) when aquifer is weakly‐to‐moderately heterogeneous, the upscaled models adequately capture the transport simulation of the FHM, despite the existence of hierarchical heterogeneity at smaller scales. When aquifer is strongly heterogeneous, the upscaled models become less accurate because lithological connectivity cannot adequately capture preferential flows; (4) three‐dimensional transport connectivities of the hierarchical aquifer differ quantitatively from those analyzed for two‐dimensional systems.

Zhang, Mingkan↗
DOE OSTI2015

CFD Modeling and Image Analysis of Exhaled Aerosols due to a Growing Bronchial Tumor: towards Non-Invasive Diagnosis and Treatment of Respiratory Obstructive Diseases

Diagnosis and prognosis of tumorigenesis are generally performed with CT, PET, or biopsy. Such methods are accurate, but have the limitations of high cost and posing additional health risks to patients. In this study, we introduce an alternative computer aided diagnostic tool that can locate malignant sites caused by tumorigenesis in a non-invasive and low-cost way. Our hypothesis is that exhaled aerosol distribution is unique to lung structure and is sensitive to airway structure vari-ations. With appropriate approaches, it is possible to locate the disease site, determine the disease severity, and subsequently formulate a targeted drug delivery plan to treat the disease. This study numerically evaluated the feasibility of the proposed breath test in an image-based lung model with varying pathological stages of a bronchial squamous tumor. Large eddy simulations and a Lagran-gian tracking approach were used to model respiratory airflows and aerosol dynamics. Respira-tions of tracer aerosols of 1 µm at a flow rate of 20 L/min were simulated, with the distributions of exhaled aerosols recorded on a filter at the mouth exit. Aerosol patterns were quantified with multiple analytical techniques such as concentration disparity, spatial scanning and fractal analysis. We demonstrated that a growing bronchial tumor induced notable variations in both the airflow and exhaled aerosol distribution. These variations became more apparent with increasing tumor severity. The exhaled aerosols exhibited distinctive pattern parameters such as spatial probability, fractal dimension, and multifractal spectrum. Results of this study show that morphometric measures of the exhaled aerosol pattern can be used to detect and monitor the pathological states of respiratory diseases in the upper airway. The proposed breath test also has the potential to locate the site of the disease, which is critical in developing a personalized, site-specific drug de-livery protocol.

Aerosol breath test↗
DOE OSTI2018

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↗
NASA NTRS2001

MHD Instability and Turbulence in the Tachocline

In this quarter we have begun simulations on the Cray T3E at PSC and we are debugging our code on the TSC. The PSC simulations are examining stratified shear turbulence with a flow-aligned magnetic field and passive tracer particles. We have conducted analysis of neutral simulations to establish a firm basis of comparison. Second-order structure functions have been computed, fit, and compared to theoretical expressions relating the dissipation fields and the structure-function-fit parameters. Agreement with high-Reynolds number observations is excellent, giving us confidence that the lower-Re simulations are relevant to higher-Re flows. We have also evaluated the neutral layer anisotropy.

Werne, Joseph↗
DOE OSTI2022

Simulating Wind Around Isolated Buildings with the System for Atmospheric Modeling

A method to represent the flow around obstacles such as buildings is implemented in the System for Atmospheric Modeling (SAM). The method, which we call a Quasi-Solid Box Method (QSBM), introduces a simple modification to the anelastic equations that forces the flow to stagnate within the obstacle's boundaries. The performance of the modified SAM is evaluated using CEDVAL (Compilation of Experimental Data for Validation of Microscale Dispersion Models) wind tunnel measurements of the wind and tracer dispersion around a single rectangular building. All major observed features of the flow are quite well reproduced. The dispersion of the tracer released at the building base is also simulated. The results of simulation of a flow around a cubic building rotated by 45 degrees relative to the flow, and around a building in the form of a cylinder of aspect ratio of one, are also presented.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Colloidal State Machines as Smart Tracers for Chemical Reactor Analysis

A widely utilized tool in reactor analysis is passive tracers that report the residence time distribution, allowing estimation of the conversion and other properties of the system. Recently, advances in microrobotics have introduced powered and functional entities with sizes comparable to some traditional tracers. This has motivated the concept of Smart Tracers that could record the local chemical concentrations, temperature, or other conditions as they progress through reactors. Herein, the design constraints and advantages of Smart Tracers by simulating their operation in a laminar flow reactor model conducting chemical reactions of various orders are analyzed. It is noted that far fewer particles are necessary to completely map even the most complex concentration gradients compared with their conventional counterparts. Design criteria explored herein include sampling frequency, memory storage capacity, and ensemble number necessary to achieve the required accuracy to inform a reactor model. Cases of severe particle diffusion and sensor noise appear to bind the functional upper limit of such probes and require consideration for future design. The results of the study provide a starting framework for applying the new technology of microrobotics to the broad and impactful set of problems classified as chemical reactor analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Numerical Mixing Suppresses Submesoscale Baroclinic Instabilities Over Sloping Bathymetry

Abstract The impacts of spurious numerical salinity mixing on the larger‐scale flow and tracer fields are characterized using idealized simulations. The idealized model is motivated by realistic simulations of the Texas‐Louisiana shelf and features oscillatory near‐inertial wind forcing. can exceed the physical mixing from the turbulence closure in frontal zones and within the mixed layer. This suggests that simulated mixing processes in frontal zones are driven largely by . Near‐inertial alongshore wind stress amplitude is varied to identify a base case that maximizes the ratio of to in simulations with no prescribed horizontal mixing. We then test the sensitivity of the base case with three tracer advection schemes (MPDATA, U3HC4, and HSIMT) and conduct ensemble runs with perturbed bathymetry. Instability growth is evaluated using the volume‐integrated eddy kinetic energy and available potential energy . While all schemes have similar total mixing, the HSIMT simulations have over double the volume‐integrated and 20% less relative to other schemes, which suppresses the release of and reduces the by roughly 25%. This results in reduced isohaline variability and steeper isopycnals, evidence that enhanced suppresses instability growth. Differences in and between the MPDATA and U3HC4 simulations are marginal. However, the U3HC4 simulations have 25% more . Experiments with variable horizontal viscosity and diffusivity coefficients show that small amounts of prescribed horizontal mixing improve the representation of the ocean state for all advection schemes by reducing the and increasing the .

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

The meteorological environment of the tropospheric ozone maximum over the tropical South Atlantic Ocean

Atmospheric flow patterns are examined over the South Atlantic Ocean where a maximum of tropospheric ozone has been observed just west of southern Africa. We investigate the flow climatology during October and perform a case study for six days during October 1989. Horizontal and vertical motions are examined and used to prepare 3D backward trajectories from the region of greatest ozone. An initially zonally symmetric distribution of ozone is treated as a passive tracer and advected by 3D flows forecast by the global model. Results from the passive tracer simulation indicate that 3D advection alone can produce a maximum of tropospheric ozone in the observed location. In addition, the trajectories suggest that by-products of biomass burning could be transported to the area of maximum ozone. Low-level flow from commonly observed regions of burning in Africa streams westward to the area of interest. Over Brazil, if the burning by-products are carried into the upper troposphere by convective process, they then could be transported eastward to the ozone feature in approximately five days. There is considerable subsidence over the tropical southern Atlantic, such that stratospheric influences also are a factor in producing the ozone maximum. Both planetary-scale and transient synoptic-scale circulation features play major roles in the various transport processes that influence the region. In summary, the observed tropospheric ozone maximum appears to be caused by a complex set of horizontal and vertical advections, transport from regions of biomass burning, and stratospheric influences.

Krishnamurti, T. N.↗
NASA NTRS2001

A High-Resolution, Three-Dimensional Model of Jupiter's Great Red Spot

The turbulent flow at the periphery of the Great Red Spot (GRS) contains many fine-scale filamentary structures, while the more quiescent core, bounded by a narrow high- velocity ring, exhibits organized, possibly counterrotating, motion. Past studies have neither been able to capture this complexity nor adequately study the effect of vertical stratification L(sub R)(zeta) on the GRS. We present results from a series of high-resolution, three-dimensional simulations that advect the dynamical tracer, potential vorticity. The detailed flow is successfully captured with a characteristic value of L(sub R) approx. equals 2000 km, independent of the precise vertical stratification profile.

Cho, James Y.-K.↗
DOE OSTI2022

Evaluation of Turbulence and Dispersion in Multiscale Atmospheric Simulations over Complex Urban Terrain during the Joint Urban 2003 Field Campaign

Abstract This paper evaluates the representation of turbulence and its effect on transport and dispersion within multiscale and microscale-only simulations in an urban environment. These simulations, run using the Weather Research and Forecasting Model with the addition of an immersed boundary method, predict transport and mixing during a controlled tracer release from the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma. This work extends the results of a recent study through analysis of turbulence kinetic energy and turbulence spectra and their role in accurately simulating wind speed, direction, and tracer concentration. The significance and role of surface heat fluxes and use of the cell perturbation method in the numerical simulation setup are also examined. Our previous study detailed the model development necessary for our multiscale simulations, examined model skill at predicting wind speeds and tracer concentrations, and demonstrated that dynamic downscaling from mesoscale to microscale through a sequence of nested simulations can improve predictions of transport and dispersion relative to a microscale-only simulation forced by idealized meteorology. Here, predictions are compared with observations to assess qualitative agreement and statistical model skill at predicting wind speed, wind direction, tracer concentration, and turbulent kinetic energy at locations throughout the city. We also investigate the scale distribution of turbulence and the associated impact on model skill, particularly for predictions of transport and dispersion. Our results show that downscaled large-scale turbulence, which is unique to the multiscale simulations, significantly improves predictions of tracer concentrations in this complex urban environment. Significance Statement Simulations of atmospheric transport and mixing in urban environments have many applications, including pollution modeling for urban planning or informing emergency response following a hazardous release. These applications include phenomena with spatial scales spanning from millimeters to kilometers. Most simulations resolve flow only within the urban area of interest, omitting larger scales of turbulence and regional influences. This study examines a method that resolves both the small and large-scale flow features. We evaluate simulation accuracy by comparing predictions with observations from an experiment involving the release of a tracer gas in Oklahoma City, Oklahoma, with emphasis on correctly modeling turbulent fluctuations. Our results demonstrate the importance of resolving large-scale flow features when predicting transport and dispersion in urban environments.

42 ENGINEERING↗
DOE OSTI2018

On visualizing material flow and precipitate evolution during probeless friction stir spot welding of an Al-Li alloy

Highlights: • Severe thermal cycle and deformation refined grains and intermetallic particles in SZ. • T{sub 1} phases fully dissolved in SZ followed by re-precipitation of GP zone and δ′ phases. • The predominant flow direction remained in spiral-like pattern during P-FSSW. • Local abrasion existed in the lap-weld interface due to flow difference, forming joining. - Abstract: The precipitate evolution and material flow visualization during probeless friction stir spot welding (P-FSSW) of an Al-Li alloy has been investigated experimentally and numerically in this work. Due to severe heat input during welding, the main strengthening phases of T{sub 1} fully dissolved in the stir zone followed by the re-precipitation of GP zone and δ′ phases, assisted with the refinement of the intermetallic particles and grains. The history of material flow was three-dimensionally visualized with a tracer, studied with X-ray tomography. The experiments and the simulation results provided a consistent interpretation of a circumferential material flow which moved inwards radially near the top surface. Upward flow in the bottom sheet improved with dwell time, causing the interface to bend. Local abrasion between the two sheets due to flow difference led to the disruption of oxide layer, and formed a connection at the interface.

36 MATERIALS SCIENCE↗