SEARCH · Search NASA
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.
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.
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.
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.
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.
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.
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.
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.
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 .
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.
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.
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.
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.
Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution
Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.
EGS Collab Experiment 1: TOUGH2-CSM Simulation of Embedded Natural Fractures and Chemical Tracer Transport and Sorption
The EGS Collab SIGMA-V project is a multi-lab and university collaborative research project that is being undertaken at the Sanford Underground Research Facility (SURF) in South Dakota. The project consists of studying stimulation, fluid-flow, and heat transfer processes at a scale of 10-20 m, which is readily amenable to detailed characterization and monitoring. One objective of the project is to establish circulation from injector to producer by hydraulically fracturing the injector. Data generated during these experiments is to be compared with predictions from coupled thermal, hydrological, mechanical, and chemical simulators. One such a simulator, TOUGH2-CSM, has been enhanced in order to simulate EGS Collab SIGMA-V project experiments. These modifications include adding tracers, the capability to model tracer sorption, and an embedded fracture formulation. A set of example problems validate our conservative tracer transport and sorption formulations. We then simulated tracer transport and thermal breakthrough for the first EGS Collab SIGMA-V experiment. This dataset includes the TOUGH2-CSM input and output files associated with the thermal and tracer simulations. A conference paper is included for additional context.
Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow
Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.
Tracer Gas Model Development and Verification in PFLOTRAN
Tracer gases, whether they are chemical or isotopic in nature, are useful tools in examining the flow and transport of gaseous or volatile species in the underground. One application is using detection of short-lived argon and xenon radionuclides to monitor for underground nuclear explosions. However, even chemically inert species, such as the noble gases, have bene observed to exhibit non-conservative behavior when flowing through porous media containing certain materials, such as zeolites, due to gas adsorption processes. This report details the model developed, implemented, and tested in the open source and massively parallel subsurface flow and transport simulator PFLOTRAN for future use in modeling the transport of adsorbing tracer gases.
Stress-Induced Anomalous Transport in Natural Fracture Networks
In this work we investigate the effects of geological stress on fluid flow and tracer transport in natural fracture networks. We show the emergence of non-Fickian (anomalous) transport from the interplay among fracture network geometry, aperture heterogeneity, and geological stress. In this study, we extract the fracture network geometry from the geological map of an actual rock outcrop, and we simulate the geomechanical behavior of fractured rock using a hybrid finite-discrete element method. We analyze the impact of stress on the aperture distribution, fluid flow field, and tracer transport properties. Both stress magnitude and orientation have strong effects on the fracture aperture field, which in turn affects fluid flow and tracer transport through the system. We observe that stress anisotropy may cause significant shear dilation along long, curved fractures that are preferentially oriented to the stress loading. This, in turn, induces preferential flow paths and anomalous early arrival of tracers. An increase in stress magnitude enhances aperture heterogeneity by introducing more small apertures, which exacerbates late-time tailing. This effect is stronger when there is higher heterogeneity in the initial aperture field. To honor the flow field with strong preferential flow paths, we extend the Bernoulli Continuous Time Random Walk model to incorporate dual velocity correlation length scales. The proposed upscaled transport model captures anomalous transport through stressed fracture networks and agrees quantitatively with the high-fidelity numerical simulations.