SEARCH · Engineering Papers

Results for “vertical velocities”

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 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Characterization of beam ion loss in high poloidal beta regime on EAST

A critical issue for achieving the integrated operation of steady-state long-pulse high-confinement (H-mode) plasmas on experimental advanced superconducting tokamak (EAST) is to improve beam ion population confinement during neutral beam injection (NBI). To study the characterization of beam ion loss and improve beam ion confinement, the steady-state long pulse scenario discharges were conducted on EAST (β p $\geqslant$ 2.0, β N $\geqslant$ 1.7, q 95 $\geqslant$ 6.7 and H 98y2 $\geqslant$ 1.1) with NBI heating. Based on neutron yield, the beam voltage and line-averaged electron density were adjusted from 50 kV to 60 kV and 4.4 × 10 19 m –3 to 5.0 × 10 19 m –3 , respectively. The results show that the dominant mechanisms of beam ion loss are shine-through loss, prompt loss, and stochastic ripple loss. The shine-through loss fraction is determined by initial velocity, flight time and entire beam path. The change in prompt loss fraction is caused by the change in the deposition of beam ions. The change in stochastic ripple loss fraction is caused by the change in the initial fraction of trapped-confined ions. Detailed physics shows that the prompt loss fraction during counter-Ip injections (~45%) is far larger than during co-Ip injections (~5%) due to the finite orbit width. The lost ions are mainly deposited on the lower divertor or below the midplane since the direction of magnetic drift is vertical down. The orbit types of prompt loss during counter-Ip injections are mainly trapped-lost and ctr-passing lost. To minimize the prompt loss fraction during counter-Ip injections, a reversed Ip configuration (rev-Ip) discharge #94758 was conducted. The result suggests that the beam ion wall load fraction during counter-Ip tangential injection (~3%) is far lower than that in normal Ip configuration (nor-Ip) discharge #94820. It is also found that the confinement of beam ion population in the counter-Ip injection #94758 was greatly improved when compared to #94820. Finally, this study can provide unique support for the improvement of beam ion population confinement and for the performance evaluation of the NBI system on EAST and future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Designing a Particle Flow Control Apparatus

Flow control within a particle-based Concentrated Solar Power (CSP) system is essential in determining the heat transfer coefficient, and therefore, the power generation capability of these systems. There are three areas where particle flow control is significant: the receivers, storage tanks, and particle-sCO2 heat exchangers. The focus of this work is on designing a new mechanism to control the flow in the particle-sCO2 heat exchangers due to the simplicity and potential cost savings when compared to the other areas of interest. The goal is for this new design to have quicker response times in terms of particle flowrate than a slide gate or flow control valve, which are designs currently used. The design resembles that of a chuck mechanism within a drill where a rotation of the sleeve elicits movement of the jaws both vertically and horizontally to close the outlet area of the nozzle. Additionally, this design will utilize the current actuator that is already used within these heat exchangers to reduce the complexity of implementation. The jaws are designed to be closed at an angle of 76˚ which is just slightly steeper than the hopper leading to the mechanism. Furthermore, this design can be tuned to limit particle bridging and other particle flow phenomena that result in blockages. The prototypes were 3D printed out of polylactic acid (PLA) and scaled up to 100%, 200%, and 400% to be able to observe the velocity profiles of the mechanism more clearly. Experiments are performed with this prototype to compare the inlet and outlet mass flow rates at different configurations of the jaws. The particles used in these experiments are 0.3mm HSP 40/70 that are commonly used in particle-based CSP systems.

14 SOLAR ENERGY↗

DOE OSTI2023

DESIGNING A PARTICLE FLOW CONTROL APPARATUS

Flow control within a particle-based Concentrated Solar Power (CSP) system is essential in determining the heat transfer coefficient, and therefore, the power generation capability of these systems. There are three areas where particle flow control is significant: the receivers, storage tanks, and particle-sCO2 heat exchangers. The focus of this paper is on designing a new mechanism to control the flow in the particle-sCO2 heat exchangers due to the simplicity and potential cost savings when compared to the other areas of interest. The goal is for this new design to have quicker response times in terms of particle flowrate than a slide gate or flow control valve, which are designs currently used. The design resembles that of a chuck mechanism within a drill where a rotation of the sleeve elicits movement of the jaws both vertically and horizontally to close the outlet area of the nozzle. Additionally, this design will utilize the current actuator that is already used within these heat exchangers to reduce the complexity of implementation. The jaws are designed to be closed at an angle of 76 ̊ which is just slightly steeper than the hopper leading to the mechanism. Furthermore, this design can be tuned to limit particle bridging and other particle flow phenomena that result in blockages. The prototypes were 3D printed out of polylactic acid (PLA) and scaled up to 100%, 200%, and 400% to be able to observe the velocity profiles of the mechanism more clearly. Experiments are performed with this prototype to compare the inlet and outlet mass flow rates at different configurations of the jaws. The particles used in these experiments are 0.3mm HSP 40/70 that are commonly used in particle-based CSP systems.

14 SOLAR ENERGY↗

DOE OSTI2026

W-Band ARM Scanning Cloud Radar (WSACR) 2nd Generation CF-Radial Spectral Data, Vertically-Pointing Scan, Cross Mode (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). Beamwidths at Ka and W bands are roughly matched at 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies that are unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the Ka-SACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. Measurements collected with the W-SACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, and linear depolarization ration.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

W-Band ARM Scanning Cloud Radar (WSACR) 2nd Generation CF-Radial Spectral Data, Vertically-Pointing Scan, Co-Polarization Mode (a1)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

W-Band ARM Scanning Cloud Radar (WSACR) 2nd Generation CF-Radial Spectral Data, Vertically-Pointing Scan, Cross-Polarization Model (a1)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗

DOE OSTI2021

Spectral Element Modeling of Acoustic to Seismic Coupling Over Topography

Acoustic waves in the atmosphere are commonly recorded on seismometers as they couple into the ground. These signals, here called ground coupled airwaves, are not commonly considered in numerical modeling of infrasound propagation, which often assumes a rigid unmeshed boundary. Starting from an analytically-tractable spherical wave model, we analyze how the coupling of an acoustic wave into a planar elastic halfspace changes over a wide range of scenarios. We use energy admittance to quantify acoustic to seismic coupling over both a planar elastic halfspace and meshed topography. Our spectral element and analytic calculations have different maxima as a function of incidence angle, with very high admittance values for near-vertical incidence (maximum ≈78%). Energy admittance calculations at shallow incidence angles are much smaller (less than 1%). In simulations over the complex topography of Sakurajima Volcano, we attribute the variable spatial pattern of energy admittance to changes in earth parameters between each model. The observed pressure difference over the simulated 15 km region appears to be < 2%. While this value is relatively small, the cumulative addition over 100s of km and multiple acoustic bounce points may be significant. Acoustic to seismic coupling along the propagation path may bias long distance yield estimates, particularly when infrasound propagates over regions with steep topography or particularly slow seismic velocities, such as alluvial planes.

58 GEOSCIENCES↗

DOE OSTI2024

Dark matter raining on DUNE and other large volume detectors

Direct detection is a powerful means of searching for particle physics evidence of dark matter (DM) heavier than about a GeV with o(kiloton) volume, low-threshold detectors. In many scenarios, some fraction of the DM may be boosted to large velocities enhancing and generally modifying possible detection signatures. We investigate the scenario where 100% of the DM is boosted at the Earth due to new attractive long-range forces. This leads to two main improvements in detection capabilities: (1) the large boost allows for detectable signatures of DM well below a GeV at large-volume neutrino detectors, such as DUNE, Super-K, Hyper-K, and JUNO, as possible DM detectors, and (2) the flux at the Earth’s surface is enhanced by a focusing effect. In addition, the model leads to a significant anisotropy in the signal with the DM flowing dominantly vertically at the Earth’s surface instead of the typical approximately isotropic DM signal. We develop the theory behind this model and also calculate realistic constraints using a detailed GENIE simulation of the signal inside detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Tidal Energy Resource Characterization, Velocity and Turbulence Measurements, Processed Data, Cook Inlet, AK, 2021

This submission contains processed datasets from a long-term deployment of 3 moorings and a transect survey of the proposed tidal energy site off the East Forelands in Cook Inlet, AK. The long-term mooring datasets were created from 8 instruments mounted on a Terrasond High Energy Oceanographic Mooring (THEOM) bottom lander and two Mid-Water Mooring (MWM) Stablemoor buoys from 1 July 2021 to 31 August 2021 (60 days). The west-most mooring (MWM1) was deployed at 60.720225 N, 151.436196 W in ~50 m of water. The middle mooring (THEOM) was deployed at 60.720703 N, 151.429500 W in ~52 m of water. The east-most buoy (MWM2) was deployed at 60.720081 N, 151.420896 W in ~50 m of water. Each Stablemoor carried three instruments: 1. A Nortek Vector acoustic Doppler velocimeter (ADV) mounted at the Stablemoor's nose. Data were recorded at 8 Hz on a 5 minute duty cycle every 20 minutes. Data was motion-corrected using the internal IMU and external ADCP bottom-track data and then bin-averaged into 4 minute bins and converted to the Principal (streamwise, cross-stream, vertical) coordinate system. (Note: 30 seconds were trimmed from the beginning and end of each 5 minute duty cycle to account for the filter end-effects from turning on and turning off the IMU.) 2. A down-looking Nortek Signature 1000 kHz acoustic Doppler current profiler (ADCP) mounted in the first Stablemoor instrument well. Data were recorded in 2 Hz with 5-beam burst and bottom-track enabled. Processed data has been averaged into 10 minute bins and converted into the Principal coordinate system. 3. An up-looking Nortek Signature 1000 kHz acoustic Doppler current profiler (ADCP) mounted in the second Stablemoor instrument well. Data were recorded at 4 Hz with 5 beam burst enabled. Processed data has been averaged into 10 minute bins and converted into the Principal coordinate system. Note: the down-facing ADCP on MWM1 failed on July 10th, 2021, only recording 9 days of data. Because ADV motion-correction required bottom track, the ADV from MWM1 also only has 9 days processed. Additionally, only 25 days of data were processed from the MWM2 ADV because it appeared to have been impacted by debris on 7/25. Two instruments were mounted on the THEOM (see MHKDR link further below for THEOM raw data): 4. A Nortek Vector acoustic Doppler velocimeter (ADV). Data were recorded at 8 Hz on a 5 minute duty cycle every 20 minutes. Data was bin-averaged into 5 minute bins, and converted to the Principal coordinate system. 5. A Nortek Signature 500 kHz acoustic Doppler current profiler (ADCP). Data were recorded in 4 Hz in the beam coordinate system from all 5 beams. Processed data has been averaged into 10 minutes bins and converted to the Principal coordinate system.

16 TIDAL AND WAVE POWER↗

DOE OSTI2023

A Moving-Wave Implementation in WRF to Study the Impact of Surface Water Waves on the Atmospheric Boundary Layer

Abstract While numerous modeling studies have focused on the interaction of ocean surface waves with the atmospheric boundary layer, most employ idealized waves that are either monochromatic or synthetically generated from a theoretical wave spectrum, and the atmospheric solvers are typically incompressible. To study wind–wave coupling in real-world scenarios, a model that can simulate both realistic meteorological and wave conditions is necessary. In this paper we describe the implementation of a moving bottom boundary condition into the Weather Research and Forecasting Model for large-eddy simulation applications. We first describe the moving bottom boundary conditions within WRF’s pressure-based vertical coordinate system. We then validate our code with idealized test cases that have analytical solutions, including flow over a monochromatic wave with and without viscosity. Finally, we present results from turbulent flows over a moving monochromatic wave with different wave ages, and demonstrate satisfactory agreement of the wave growth rate with results from the literature. We also compare atmospheric stress and wind parameters from two physically equivalent cases. The first specifies a wind moving in the same direction as a propagating wave, while the second involves a stationary wave with the wind adjusted such that the wind relative to the wave is the same as in the first case. Results indicate that the velocity and Reynolds stress profiles for the two cases match, further validating the moving bottom implementation.

17 WIND ENERGY↗

DOE OSTI2025

An intercomparison of wall fluxes in a turbulent thermal convection chamber: Direct numerical simulations and wall-modeled large-eddy simulations enhanced by machine learning

Thermal convection in a closed chamber is driven by a warm bottom, a cold top, and side walls at various temperatures. Although wall fluxes are the source of convection energy, accurately modeling these fluxes (i.e., the wall model) is challenging. In large-eddy simulations (LESs), many wall models are traditionally derived from the canonical boundary layer, which may be unsuitable for thermal convection bounded by both horizontal and vertical walls. This study conducts a model intercomparison of dry convection in a cubic-meter chamber using three direct numerical simulations (DNSs) and four LESs with different wall models. The LESs employ traditional wall models, a new wall model employing physics-aware neural networks, and a refined grid near the walls. The experiment involves four cases with varying sidewall temperatures. Our results show that LESs capture the main flow features and the trends of mean fluxes. The physics-aware neural networks and refined wall grids can improve the temporally averaged local fluxes when the large-scale circulation has a preferred direction. Even without the local improvement of wall fluxes, the LES flow quantities (temperature and velocities) can still largely match those in DNSs, provided the mean flux largely matches the DNSs. Additionally, DNSs reveal that a variation in corner treatments has minimal impacts on the flow quantities away from corners. Finally, LESs underestimate the mean fluxes of the entire wall due to their inability to resolve corner regions, but their mean flux away from the corner can better match DNS.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Distributed Acoustic Sensing (DAS) Seismic Monitoring of CO 2 Injected for Enhanced Oil Recovery in Northern Michigan

As part of the Midwest Regional Carbon Sequestration Partnership (MRCSP) Phase III project, a monitoring study was conducted to assess the effectiveness of DAS (Distributed Acoustic Sensing) -based VSP (Vertical Seismic Profiling) technology for delineating CO 2 injected into the Silurian-age pinnacle reefs in northern Michigan, the host rocks for the MRCSP Phase III demonstration project. The DAS VSP study was conducted in the Chester 16 reef, one of several reefs in Otsego County Michigan that is operated by Core Energy, LLC, of Traverse City, Michigan.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

On Sediment Mass Flux Directionality Induced by Yawed Permeable Vanes under Near-Critical Mobility Conditions

Not provided.

Engineering↗

DOE OSTI2024

A strainmeter array as the fulcrum of novel observatory sites along the Alto Tiberina Near Fault Observatory

Fault slip is a complex natural phenomenon involving multiple spatiotemporal scales from seconds to days to weeks. To understand the physical and chemical processes responsible for the full fault slip spectrum, a multidisciplinary approach is highly recommended. The Near Fault Observatories (NFOs) aim at providing high-precision and spatiotemporally dense multidisciplinary near-fault data, enabling the generation of new original observations and innovative scientific products. The Alto Tiberina Near Fault Observatory is a permanent monitoring infrastructure established around the Alto Tiberina fault (ATF), a 60 km long low-angle normal fault (mean dip 20°), located along a sector of the Northern Apennines (central Italy) undergoing an extension at a rate of about 3 mm yr –1 . The presence of repeating earthquakes on the ATF and a steep gradient in crustal velocities measured across the ATF by GNSS stations suggest large and deep (5–12 km) portions of the ATF undergoing aseismic creep. Both laboratory and theoretical studies indicate that any given patch of a fault can creep, nucleate slow earthquakes, and host large earthquakes, as also documented in nature for certain ruptures (e.g., Iquique in 2014, Tōhoku in 2011, and Parkfield in 2004). Nonetheless, how a fault patch switches from one mode of slip to another, as well as the interaction between creep, slow slip, and regular earthquakes, is still poorly documented by near-field observation. With the strainmeter array along the Alto Tiberina fault system (STAR) project, we build a series of six geophysical observatory sites consisting of 80–160 m deep vertical boreholes instrumented with strainmeters and seismometers as well as meteorological and GNSS antennas and additional seismometers at the surface. By covering the portions of the ATF that exhibits repeated earthquakes at shallow depth (above 4 km) with these new observatory sites, we aim to collect unique open-access data to answer fundamental questions about the relationship between creep, slow slip, dynamic earthquake rupture, and tectonic faulting.

58 GEOSCIENCES↗

DOE OSTI2023

Investigation of carbon incorporation in laser-assisted MOCVD of GaN

Background carbon (C) impurity incorporation in metalorganic chemical vapor deposition (MOCVD) grown gallium nitride (GaN) represents one of the major issues in further improving GaN vertical power device performance. This work presents a laser-assisted MOCVD (LA-MOCVD) technique to address the high-C issue in MOCVD homoepitaxial GaN under different growth rate (Rg) regimes and studies the correlations between [C] and Rg. [C] in LA-MOCVD GaN is reduced by 50%–90% as compared to the conventional MOCVD GaN for a wide growth rate range between 1 and 16 μm/h. A mass-transport based model is developed to understand the C incorporation at different Rg regimes. The results obtained from the developed model are in good agreement with experimental data. The model further reveals that LA-MOCVD effectively suppresses C incorporation by reducing the active C species in the gas phase. Moreover, high step velocity in step flow growth mode can facilitate C incorporation at fast Rg, exhibiting steeper C increase. The theoretical model indicates that [C] can be suppressed below 1016 cm−3 with a fast growth rate (Rg) of 10 μm/h by utilizing higher power LA-MOCVD and freestanding GaN substrates with larger off-cut angles.

Physics↗

DOE OSTI2020

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗