SEARCH · Engineering Papers

Results for “complex fluids”

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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRS1992

A gridless Euler/Navier-Stokes solution algorithm for complex two-dimensional applications

The development of a gridless computational fluid dynamics (CFD) method for the solution of the two-dimensional Euler and Navier-Stokes equations is described. The method uses only clouds of points and does not require that the points be connected to form a grid as is necessary in conventional CFD algorithms. The gridless CFD approach appears to resolve the problems and inefficiencies encountered with structured or unstructured grid methods. As a result, the method offers the greatest potential for accurately and efficiently solving viscous flows about complex aircraft configurations. The method is described in detail, and calculations are presented for standard Euler and Navier-Stokes cases to assess the accuracy and efficiency of the capability.

Batina, John T.↗

NASA NTRS1996

Segmentation of Unstructured Datasets

Datasets generated by computer simulations and experiments in Computational Fluid Dynamics tend to be extremely large and complex. It is difficult to visualize these datasets using standard techniques like Volume Rendering and Ray Casting. Object Segmentation provides a technique to extract and quantify regions of interest within these massive datasets. This thesis explores basic algorithms to extract coherent amorphous regions from two-dimensional and three-dimensional scalar unstructured grids. The techniques are applied to datasets from Computational Fluid Dynamics and from Finite Element Analysis.

Bhat, Smitha↗

DOE OSTI2020

Nanofabrication of synthetic nanoporous geomaterials: from nanoscale-resolution 3D imaging to nano-3D-printed digital (shale) rock

Abstract Advances in imaging have made it possible to view nanometer and sub-nanometer structures that are either synthesized or that occur naturally. It is believed that fluid dynamic and thermodynamic behavior differ significantly at these scales from the bulk. From a materials perspective, it is important to be able to create complex structures at the nanometer scale, reproducibly, so that the fluid behavior may be studied. New advances in nanoscale-resolution 3D-printing offer opportunities to achieve this goal. In particular, additive manufacturing with two-photon polymerization allows creation of intricate structures. Using this technology, a creation of the first nano-3D-printed digital (shale) rock is reported. In this paper, focused ion beam-scanning electron microscopy (FIB-SEM) nano-tomography image dataset was used to reconstruct a high-resolution digital rock 3D model of a Marcellus Shale rock sample. Porosity of this 3D model has been characterized and its connected/effective pore system has been extracted and nano-3D-printed. The workflow of creating this novel nano-3D-printed digital rock 3D model is described in this paper.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2010

Multi-Fluid Block-Adaptive-Tree Solar Wind Roe-Type Upwind Scheme: Magnetospheric Composition and Dynamics During Geomagnetic Storms, Initial Results

The magnetosphere contains a significant amount of ionospheric O{+}, particularly during geomagnetically active times. The presence of ionospheric plasma in the magnetosphere has a notable impact on magnetospheric composition and processes. We present a new multifluid MHD version of the BATS-R-US model of the magnetosphere to track the fate and consequences of ionospheric outflow. The multi-fluid MHD equations are presented as are the novel techniques for overcoming the formidable challenges associated with solving them. Our new model is then applied to the May 4, 1998 and March 31, 2001 geomagnetic storms. The results are juxtaposed with traditional single- fluid MHD and multispecies MHD simulations from a previous study, thereby allowing us to assess the benefits of using a more complex model with additional physics. We find that our multi-fluid MHD model (with outflow) gives comparable results to the multi-species MHD model (with outflow), including a more strongly negative Dst, reduced CPCP, and a drastically improved magnetic field at geosynchronous orbit, as compared to single-fluid MHD with no outflow. Significant differences in composition and magnetic field are found between the multi-species and multi-fluid approach further away from the Earth. We further demonstrate the ability to explore pressure and bulk velocity differences between H{+} and O(+}, which is not possible when utilizing the other techniques considered.

Gkocer, A.↗

NASA NTRS2016

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗

DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗

DOE OSTI2025

An experimental study on gas-liquid phase fluid migration in hydrate-bearing sediments during hydrate dissociation

Natural gas hydrate production tests face problems such as severe sand blockage, poor gas-liquid phase separation, and significant land subsidence. This is because of the insufficient understanding of the complex phase transition and gas-liquid multi-phase fluid migration during hydrate dissociation. In hydrate-bearing sediment systems, hydrate phase transition couples with gas-liquid fluid migration. The phase transition causes changes in pore structure, which in turn modifies porous infiltration parameters and fluid flow capacity. Meanwhile, alterations in phase interfaces affect key parameters like surface tension and wettability. Gas-liquid fluid migration influences heat and mass transfer, thus affecting phase equilibrium and dissociation rates. To bridge the gap in describing gas-liquid fluid migration during hydrate dissociation in experiments, this research innovatively integrated an unsteady-state gas displacement by water and a quantitative hydrate dissociation process, independently developed a multi-phase seepage experimental system suitable for hydrate dissociation and determined the relationship between seepage parameters and hydrate saturation under different porosity. The results are as follows: a) Core samples with higher initial porosity show a greater recovery rate of fluid flow capacity. b) The retarding effect of multi-phase fluid has a more significant impact on the migration of the wetting phase fluid (water) than that of the non-wetting phase fluid (methane). c) During hydrate dissociation, the evolution of absolute permeability shows an “S-shaped” pattern, and the evolution of relative permeability shows a “wiring-harness” pattern. In conclusion, the findings can provide a theoretical basis for preventing geological disasters and for geotechnical engineering design during hydrate production.

58 GEOSCIENCES↗

NASA NTRS2005

Three Dimensional Multi-Fluid Simulation of Comet Halley

The Cometary Atmosphere Simulator (CASIM) is designed to simulate in 2-D and 3-D the complex interaction between the cometary atmosphere and the hypersonic solar wind using a multi-fluid approach. Our simulator is based on the solution of multi-fluid equations using an efficient adaptive Cartesian mesh. It is designed to use the capabilities of highly parallel super-cluster computers.

Benna, M.↗

NASA NTRS2004

Extraction of Organic Molecules from Terrestrial Material: Quantitative Yields from Heat and Water Extractions

In the robotic search for life on Mars, different proposed missions will analyze the chemical and biological signatures of life using different platforms. The analysis of samples via analytical instrumentation on the surface of Mars has thus far only been attempted by the two Viking missions. Robotic arms scooped relogith material into a pyrolysis oven attached to a GC/MS. No trace of organic material was found on any of the two different samples at either of the two different landing sites. This null result puts an upper limit on the amount of organics that might be present in Martian soil/rocks, although the level of detection for each individual molecular species is still debated. Determining the absolute limit of detection for each analytical instrument is essential so that null results can be understood. This includes investigating the trade off of using pyrolysis versus liquid solvent extraction to release organic materials (in terms of extraction efficiencies and the complexity of the sample extraction process.) Extraction of organics from field samples can be accomplished by a variety of methods such utilizing various solvents including HCl, pure water, supercritical fluid and Soxhelt extraction. Utilizing 6N HCl is one of the most commonly used method and frequently utilized for extraction of organics from meteorites but it is probably infeasible for robotic exploration due to difficulty of storage and transport. Extraction utilizing H2O is promising, but it could be less efficient than 6N HCl. Both supercritical fluid and Soxhelt extraction methods require bulky hardware and require complex steps, inappropriate for inclusion on rover spacecraft. This investigation reports the efficiencies of pyrolysis and solvent extraction methods for amino acids for different terrestrial samples. The samples studied here, initially created in aqueous environments, are sedimentary in nature. These particular samples were chosen because they possibly represent one of the best terrestrial analogs of Mars and they represent one of the absolute best case scenarios for finding organic molecules on the Martian surface.

Beegle, L. W.↗

NASA NTRS1993

Computational heat transfer analysis for oscillatory channel flows

An accurate finite-difference scheme has been utilized to investigate oscillatory, laminar and incompressible flow between two-parallel-plates and in circular tubes. The two-parallel-plates simulate the regenerator of a free-piston Stirling engine (foil type regenerator) and the channel wall was included in the analysis (conjugate heat transfer problem). The circular tubes simulate the cooler and heater of the engine with an isothermal wall. The study conducted covered a wide range for the maximum Reynolds number (from 75 to 60,000), Valensi number (from 2.5 to 700), and relative amplitude of fluid displacement (0.714 and 1.34). The computational results indicate a complex nature of the heat flux distribution with time and axial location in the channel. At the channel mid-plane we observed two thermal cycles (out of phase with the flow) per each flow cycle. At this axial location the wall heat flux mean value, amplitude and phase shift with the flow are dependent upon the maximum Reynolds number, Valensi number and relative amplitude of fluid displacement. At other axial locations, the wall heat flux distribution is more complex.

Ibrahim, Mounir↗

DOE OSTI2021

Internally cooled membrane-based absorber for dehumidification and water heating: Validated model and simulation study

This report presents a comprehensive model for advanced membrane-based absorber components accounting for three separate fluid streams, one of them providing internal cooling, with complex flow patterns. The model was implemented with various liquid desiccants and their properties (including LiCl, CaCl 2 , and [emim][OMS] ionic liquid [IL]). The model has been validated using data from two laboratory prototype absorbers, one with CaCl 2 as the working fluid and one with an IL as the working fluid. Entropy analysis was further carried out to verify the model and understand effect of internal cooling on the entropy variation of three fluid streams. The developed model and codes are expected to enable detailed configuration optimization and provide in-depth understanding of three-fluid heat and mass exchanger (HMX) performance. The paper also describes parametric studies of the HMX components in utilizing the latent heat removed in space cooling to heat water and numerically explores the unique benefit of its application in a semi-open absorption heat pump water heater. The results show that a single three-fluid HMX has the potential to achieve simultaneous dehumidification and water heating efficiently at cost effectiveness, particularly in hot and humid climates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Improving the reliability of machine learned potentials for modeling inhomogeneous liquids

The atomic-scale response of inhomogeneous fluids at interfaces and surrounding solute particles plays a critical role in governing chemical, electrochemical, and biological processes. Classical molecular dynamics simulations have been applied extensively to simulate the response of fluids to inhomogeneities directly, but are limited by the accuracy of the underlying interatomic potentials. Here, we use neural network potentials (NNPs) trained to ab initio simulations to accurately predict the inhomogeneous responses of two distinct fluids: liquid water and molten NaCl. Although NNPs can be readily trained to model complex bulk systems across a range of state points, we show that to appropriately model a fluid's response at an interface, relevant inhomogeneous configurations must be included in the training data. In order to sufficiently sample appropriate configurations of such inhomogeneous fluids, we develop protocols based on molecular dynamics simulations in the presence of external potentials. We demonstrate that NNPs trained on inhomogeneous fluid configurations can more accurately predict several key properties of fluids—including the density response, surface tension and size-dependent cavitation free energies—for liquid water and molten NaCl, compared to both empirical interatomic potentials and NNPs that are not trained on such inhomogeneous configurations. This work therefore provides a first demonstration and framework to extract the response of inhomogeneous fluids from first principles for classical density-functional treatment of fluids free from empirical potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Forest structural complexity and ignition pattern influence simulated prescribed fire effects

Background: Forest structural characteristics, the burning environment, and the choice of ignition pattern each influence prescribed fire behaviors and resulting fire effects; however, few studies examine the influences and interactions of these factors. Understanding how interactions among these drivers can influence prescribed fire behavior and effects is crucial for executing prescribed fires that can safely and effectively meet management objectives. To analyze the interactions between the fuels complex and ignition patterns, we used FIRETEC, a three-dimensional computational fluid dynamics fire behavior model, to simulate fire behavior and effects across a range of horizontal and vertical forest structural complexities. For each forest structure, we then simulated three different prescribed fires each with a unique ignition pattern: strip-head, dot, and alternating dot. Results: Forest structural complexity and ignition pattern affected the proportions of simulated crown scorch, consumption, and damage for prescribed fires in a dry, fire-prone ecosystem. Prescribed fires in forests with complex canopy structures resulted in increased crown consumption, scorch, and damage compared to less spatially complex forests. The choice of using a strip-head ignition pattern over either a dot or alternating-dot pattern increased the degree of crown foliage scorched and damaged, though did not affect the proportion of crown consumed. We found no evidence of an interaction between forest structural complexity and ignition pattern on canopy fuel consumption, scorch, or damage. Conclusions: We found that forest structure and ignition pattern, two powerful drivers of fire behavior that forest managers can readily account for or even manipulate, can be leveraged to influence fire behavior and the resultant fire effects of prescribed fire. These simulation findings have critical implications for how managers can plan and perform forest thinning and prescribed burn treatments to meet risk management or ecological objectives.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Surfactant-Aided Low-Salinity Waterflooding for Low-Temperature Carbonate Reservoirs

Carbonate reservoirs tend to be oil-wet/mixed-wet and heterogeneous because of mineralogy and diagenesis. The objective of this study is to improve oil recovery in low-temperature dolomite reservoirs using low-salinity and surfactant-aided spontaneous imbibition. The low-salinity brine composition was optimized using ζ-potential measurements, contact-angle (CA) experiments, and a novel wettability-alteration measure. Significant wettability alteration was observed on dolomite rocks at a salinity of 2,500 ppm. We evaluated 37 surfactants by performing CA, interfacial-tension (IFT), and spontaneous-imbibition experiments. Here, three (quaternary ammonium) cationic and one (sulfonate) anionic surfactants showed significant wettability alteration and produced 43–63% of original oil in place (OOIP) by spontaneous imbibition. At a low temperature (35°C), oil recovery by low-salinity effect is small compared with that by wettability-altering surfactants. Coreflood tests were performed with a selected low-salinity cationic surfactant solution. A novel coreflood was proposed that modeled heterogeneity and dynamic imbibition into low-permeability regions. The results of the “heterogeneous” coreflood were consistent with that of spontaneous-imbibition tests. These experiments demonstrated that a combination of low-salinity brine and surfactants can make originally oil-wet dolomite rocks more water-wet and improve oil recovery from regions bypassed by waterflood at a low temperature of 35°C.

02 PETROLEUM↗

DOE OSTI2023

A Deep Learning-Based Workflow for Fast Prediction of 3D State Variables in Geological Carbon Storage: A Dimension Reduction Approach

In this study, we used deep learning techniques, which are a form of artificial intelligence, to create fast and effective models for predicting how fluids flow in underground geological formations. This is important for managing geological carbon storage, a method used to fight climate change by storing carbon dioxide underground. The challenge lies in the complex nature of these underground spaces and the large amount of data needed to accurately simulate them. To overcome these issues, we developed a new workflow that reduces the data’s complexity before training the deep learning model and then reconstructs the predicted results in their original form. We also proposed a unique approach to handle the specific complexities found in 3D saturation fields, a crucial aspect of fluid flow prediction. We tested our method using real-world data from the Gulf of Mexico. Our results show that our approach not only accurately predicts fluid behavior but also significantly reduces computation time. This will greatly improve real-time decision-making and risk assessment in large-scale geological carbon storage operations.

Wang, Hongsheng↗

NASA NTRS1997

Constrained Multipoint Aerodynamic Shape Optimization Using an Adjoint Formulation and Parallel Computers

An aerodynamic shape optimization method that treats the design of complex aircraft configurations subject to high fidelity computational fluid dynamics (CFD), geometric constraints and multiple design points is described. The design process will be greatly accelerated through the use of both control theory and distributed memory computer architectures. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on a higher order CFD method. In order to facilitate the integration of these high fidelity CFD approaches into future multi-disciplinary optimization (NW) applications, new methods must be developed which are capable of simultaneously addressing complex geometries, multiple objective functions, and geometric design constraints. In our earlier studies, we coupled the adjoint based design formulations with unconstrained optimization algorithms and showed that the approach was effective for the aerodynamic design of airfoils, wings, wing-bodies, and complex aircraft configurations. In many of the results presented in these earlier works, geometric constraints were satisfied either by a projection into feasible space or by posing the design space parameterization such that it automatically satisfied constraints. Furthermore, with the exception of reference 9 where the second author initially explored the use of multipoint design in conjunction with adjoint formulations, our earlier works have focused on single point design efforts. Here we demonstrate that the same methodology may be extended to treat complete configuration designs subject to multiple design points and geometric constraints. Examples are presented for both transonic and supersonic configurations ranging from wing alone designs to complex configuration designs involving wing, fuselage, nacelles and pylons.

Reuther, James↗

DOE OSTI2021

Large Eddy Simulations of Turbulent and Buoyant Flows in Urban and Complex Terrain Areas Using the Aeolus Model

Fast and accurate predictions of the flow and transport of materials in urban and complex terrain areas are challenging because of the heterogeneity of buildings and land features of different shapes and sizes connected by canyons and channels, which results in complex patterns of turbulence that can enhance material concentrations in certain regions. To address this challenge, we have developed an efficient three-dimensional computational fluid dynamics (CFD) code called Aeolus that is based on first principles for predicting transport and dispersion of materials in complex terrain and urban areas. The model can be run in a very efficient Reynolds average Navier–Stokes (RANS) mode or a detailed large eddy simulation (LES) mode. The RANS version of Aeolus was previously validated against field data for tracer gas and radiological dispersal releases. As a part of this work, we have validated the Aeolus model in LES mode against two different sets of data: (1) turbulence quantities measured in complex terrain at Askervein Hill; and (2) wind and tracer data from the Joint Urban 2003 field campaign for urban topography. As a third set-up, we have applied Aeolus to simulate cloud rise dynamics for buoyant plumes from high-temperature explosions. For all three cases, Aeolus LES predictions compare well to observations and other models. These results indicate that Aeolus LES can be used to accurately simulate turbulent flow and transport for a wide range of applications and scales.

42 ENGINEERING↗