Search NASA⌕ Search

SEARCH · Search NASA

Results for “potential vorticity”

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.

126 records · Page 7

Localized Flexoelectric Effect Around Ba(CuNb) Nano‐Clusters in Epitaxial BiFeO 3 Films for Enhancement of Electric and Multiferroic Properties

Abstract Room‐temperature (RT) multiferroic materials have received significant research attention for various potential applications; however, their properties are not suitable for real‐world implementation. In this study, a nano‐scale localized flexoelectric effect is introduced to enhance the RT multiferroic performance of epitaxial bismuth iron oxide (BiFeO 3 ; BFO) thin films by embedding 10 mol% Ba(Cu 1/3 Nb 2/3 )O 3 (BCN) nano‐clusters into the host BFO film, which originally has a rhombohedral crystal structure. By utilizing nano‐clustering, a large out‐of‐plane coherent strain is localized around the nano‐clusters, resulting in a highly strained tetragonality of the BFO structure; subsequently, the films exhibit peculiar types of domains and domain walls, such as nano‐scale rotational vortices and antiparallel dipole configurations. These peculiar domain structures, which originate from the localized flexoelectric effect at the nano‐scale, enable excellent ferroelectric, ferromagnetic, and RT multiferroic magnetoelectric coupling. This study reveals that the local variation in the localized flexoelectric field around nano‐clusters considerably impacts the formation of unusual domain‐wall structures. This suggests that the controlled introduction of nano‐clusters with different crystal structures is promising for achieving the desired multiferroic properties.

Song, Hyunseok↗

Identifying Characteristics of Wildfire Towers and Troughs

Wildfire behavior is dictated by the complex interaction of numerous physical phenomena including dynamic ambient and fire-induced winds, heat transfer, aerodynamic drag on the wind by the fuel and combustion. These phenomena create complex feedback effects between the fire and its surroundings. In this study, we aim to study the mechanisms by which buoyant flame dynamics along with vortical motions and instabilities control wildfire propagation. Specifically, this study employs a suite of simulations conducted with the physics-based coupled fire-atmosphere behavior model (FIRETEC). The simulations are initialized with a fire line and the fires are allowed to propagate on a grass bed, where the fuel heights and wind conditions are varied systematically. Flow variables are extracted to identify the characteristics of the alternating counter-rotational vortices, called towers and troughs, that drive convective heat transfer and fire spread. These vortices have previously been observed in wildfires and laboratory fires, and have also been observed to arise spontaneously in FIRETEC due to the fundamental physics incorporated in the model. However, these past observations have been qualitative in nature and no quantitative studies can be found in the literature which connected these coherent structures fundamental to fire behavior with the constitutive flow variables. To that end, a variety of state variables are examined in the context of these coherent structures under various wind profile and grass height conditions. Identification of various correlated signatures and fire-atmosphere feedbacks in simulations provides a hypothesis that can be tested in future observational or experimental efforts, potentially assisting experimental design, and can aid in the interpretation of data from in situ detectors.

54 ENVIRONMENTAL SCIENCES↗

The emergence of three-dimensional chiral domain walls in polar vortices

Abstract Chirality or handedness of a material can be used as an order parameter to uncover the emergent electronic properties for quantum information science. Conventionally, chirality is found in naturally occurring biomolecules and magnetic materials. Chirality can be engineered in a topological polar vortex ferroelectric/dielectric system via atomic-scale symmetry-breaking operations. We use four-dimensional scanning transmission electron microscopy (4D-STEM) to map out the topology-driven three-dimensional domain walls, where the handedness of two neighbor topological domains change or remain the same. The nature of the domain walls is governed by the interplay of the local perpendicular (lateral) and parallel (axial) polarization with respect to the tubular vortex structures. Unique symmetry-breaking operations and the finite nature of domain walls result in a triple point formation at the junction of chiral and achiral domain walls. The unconventional nature of the domain walls with triple point pairs may result in unique electrostatic and magnetic properties potentially useful for quantum sensing applications.

36 MATERIALS SCIENCE↗

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗

Direct current magnetic Hall probe technique for measurement of field penetration in thin film superconductors for superconducting radio frequency resonators

Superconducting Radio Frequency (SRF) cavities used in particle accelerators are typically formed from or coated with superconducting materials. Currently, high purity niobium is the material of choice for SRF cavities that have been optimized to operate near their theoretical field limits. This brings about the need for significant R & D efforts to develop next generation superconducting materials that could outperform Nb and keep up with the demands of new accelerator facilities. To achieve high quality factors and accelerating gradients, the cavity material should be able to remain in the superconducting Meissner state under a high RF magnetic field without penetration of quantized magnetic vortices through the cavity wall. Therefore, the magnetic field at which vortices penetrate a superconductor is one of the key parameters of merit of SRF cavities. Techniques to measure the onset of magnetic field penetration on thin film samples need to be developed to mitigate the issues with the conventional magnetometry measurements that are strongly influenced by the film orientation and shape and edge effects. In this work, we report the development of an experimental setup to measure the field of full flux penetration through films and multi-layered superconductors. Our system combines a small superconducting solenoid that can generate a magnetic field of up to 500 mT at the sample surface and three Hall probes to detect the full flux penetration through the superconductor. We report this setup can be used to study alternative materials that could potentially outperform niobium, as well as superconductor–insulator–superconductor (SIS) multilayer coatings on niobium.

47 OTHER INSTRUMENTATION↗

A new era in ferroelectrics

Topological structures in ferroic materials have drawn great interest in recent years due to the richness of the underlying physics and the potential for applications in next generation electronics. Recent advances in atomically precise thin-film materials synthesis and characterization of structural/physical phenomena at unprecedented length/energy/time scales have enabled us to study exotic phases and their associated physics [Rößler et al., Nature 442, 797 (2006); S. Das, Nature 568, 368 (2019); Yadav et al., Nature 530, 198 (2016); and Stoica et al., Nat. Mater. 18, 377 (2019)]. It is appropriate that, in the second century of ferroelectrics, some dramatic discoveries are propelling the field into directions heretofore unimaginable. In this review, we explore the recent progress in ferroelectric-oxide superlattices in which researchers can control structure and physical properties through the application of epitaxial strain, layer thickness, temperature, electric field, etc. We provide a discussion of exotic topological structures (e.g., closure domains, vortices, polar skyrmions, and other exotic phases) and associated functionalities in ferroelectric/paraelectric superlattices. We conclude with a brief overview of and prospects for how the field may evolve in the coming years.

36 MATERIALS SCIENCE↗

Novel Air System for 300 kW Heavy Duty Fuel Cell

Hydrogen fuel cell powered vehicles for heavy duty trucks is a potential path for reducing heavy duty vehicle emissions in the future. The air handling system delivers the proper amount of air (oxygen) to react with fuel (hydrogen) in the fuel cell to produce power. Air delivery requires significant power and is the largest parasitic loss for a 300 kW fuel cell. Today’s systems use an electric motor to power an air compressor that supplies oxygen to the fuel cell stack. In addition to parasitic power loss, hydrogen fuel cell systems often have reliability issues associated with the air handling system. Reliability is of significant concern for heavy duty applications (especially long-haul applications). This project aims to improve both the electrical power consumption and reliability of hydrogen fuel cell air handling systems to meet the needs of heavy duty on-highway vehicle applications. The air handling is provided by a twin vortices series (TVS) compressor in addition to adding a TVS expander to recover waste heat energy back into the compressor. The final configuration includes a 600 V, 39 kW motor connected with a single shaft to the compressor and expander. This configuration reduced the total electrical power consumption from 48.6 kW to 37 kW at full load, 27 kW to 9 kW at half load and 0.44 kW to 0.22 kW at idle. The response time requirement was to be less than 2 seconds while the final demonstration yielded 0.62 sec. The study includes additional technologies including water dosing into the compressor, a recuperator and elimination of the intercooler.

Reich, Evan [Eaton Corporation, Southfield, MI (Un↗

Thermalization and criticality on an analogue–digital quantum simulator

Abstract Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators 1,2 . Unlocking the full potential of such systems towards this goal requires flexible initial state preparation, precise time evolution and extensive probes for final state characterization. Here we present a quantum simulator comprising 69 superconducting qubits that supports both universal quantum gates and high-fidelity analogue evolution, with performance beyond the reach of classical simulation in cross-entropy benchmarking experiments. This hybrid platform features more versatile measurement capabilities compared with analogue-only simulators, which we leverage here to reveal a coarsening-induced breakdown of Kibble–Zurek scaling predictions 3 in theXYmodel, as well as signatures of the classical Kosterlitz–Thouless phase transition 4 . Moreover, the digital gates enable precise energy control, allowing us to study the effects of the eigenstate thermalization hypothesis 5–7 in targeted parts of the eigenspectrum. We also demonstrate digital preparation of pairwise-entangled dimer states, and image the transport of energy and vorticity during subsequent thermalization in analogue evolution. These results establish the efficacy of superconducting analogue–digital quantum processors for preparing states across many-body spectra and unveiling their thermalization dynamics.

Science & Technology - Other Topics↗

Near-wake dynamics of a vertical-axis turbine

Cross-flow, or vertical-axis, turbines are a promising technology for capturing kinetic energy in wind or flowing water and their inherently unsteady fluid mechanics present unique opportunities for control optimization of individual rotors or arrays. To explore the potential for beneficial interactions between turbines in an array, as well as to characterize important cycle-to-cycle variations, coherent structures in the wake of a single two-bladed cross-flow turbine are examined using planar stereo particle image velocimetry in a water channel experiment. There are three main objectives in the present work. First, the mean wake structure of this high chord-to-radius ratio rotor is described, compared with previous studies, and a simple explanation for observed wake deflection is presented. Second, the unsteady flow is then analysed via the triple decomposition, with the periodic component extracted using a combination of traditional techniques and a novel implementation of the optimized dynamic mode decomposition. The latter method is shown to outperform conditional averaging and Fourier methods, as well as uncover frequencies suggesting a transition to bluff-body shedding in the far wake. Third, vorticity and finite-time Lyapunov exponents are then employed to further analyse the oscillatory wake component. Vortex streets on both sides of the wake are identified, and their formation mechanisms and effects on the mean flow are discussed. Strong axial (vertical) flow is observed in vortical structures shed on the retreating side of the rotor where the blades travel downstream. Time-resolved tracking of these vortices is performed, which demonstrates that vortex trajectories have significant rotation-to-rotation variation within one diameter downstream. This variability suggests it would be challenging to harness or avoid such structures at greater downstream distances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Monolayer Superconductivity and Tunable Topological Electronic Structure at the Fe(Te,Se)/Bi 2 Te 3 Interface

The interface between 2D topological Dirac states and an s-wave superconductor is expected to support Majorana-bound states (MBS) that can be used for quantum computing applications. Realizing these novel states of matter and their applications requires control over superconductivity and spin-orbit coupling to achieve spin-momentum-locked topological interface states (TIS) which are simultaneously superconducting. While signatures of MBS have been observed in the magnetic vortex cores of bulk FeTe 0.55 Se 0.45 , inhomogeneity and disorder from doping make these signatures unclear and inconsistent between vortices. In this work, superconductivity is reported in monolayer (ML) FeTe 1–y Se y (Fe(Te,Se)) grown on Bi 2 Te 3 by molecular beam epitaxy (MBE). Spin and angle-resolved photoemission spectroscopy (SARPES) directly resolve the interfacial spin and electronic structure of Fe(Te,Se)/Bi 2 Te 3 heterostructures. For y = 0.25, the Fe(Te,Se) electronic structure is found to overlap with the Bi 2 Te 3 TIS and the desired spin-momentum locking is not observed. In contrast, for y = 0.1, reduced inhomogeneity measured by scanning tunneling microscopy (STM) and a smaller Fe(Te,Se) Fermi surface with clear spin-momentum locking in the topological states are found. Hence, it is demonstrated that the Fe(Te,Se)/Bi 2 Te 3 system is a highly tunable platform for realizing MBS where reduced doping can improve characteristics important for Majorana interrogation and potential applications.

36 MATERIALS SCIENCE↗

Needle-Like Ferroelastic Domains in Individual Ferroelectric Nanoparticles

Superior structural, physical and electronic properties make ferroelectric nanocrystals essential in enabling a range of next-generation devices. Ferroelectric responses are determined by crystal structure and domain morphology. The ability to reversibly displace, create and annihilate elastic domains is critical to device applications. Whereas electric-field control has been demonstrated for ferroelectric 180° surface domain walls and vortices, similar control of ferroelastic domains and domain boundaries within individual nanocrystals remains challenging. Using controlled external compressive and tensile axial stress, deterministic and reversible control of highly mobile ferroelastic domains and axial polarization in three dimensions is demonstrated in this study. While many studies exist on ferroelastic domains in thin films and bulk, little is known about ferroelastic interactions at the single nanoparticle level, especially involving domain boundaries. Through combining Bragg coherent X-ray diffractive imaging and Landau theory, strain gradients in individual BaTiO3 nanocrystals are shown to stabilize needle-like ferroelastic twin domains. These domains are highly labile under applied axial stress, producing a locally enhanced electric polarization mediated by a ferroelectric phase transition. The efficacy of Bragg coherent X-ray diffractive imaging in studying in operando domains in three-dimensions is demonstrated, while synergy with theory provides a paradigm for domain boundary engineering and potential for nanoscale functional devices.

36 MATERIALS SCIENCE↗

Experimental investigation of turbulent wake flows in a helically wrapped rod bundle in presence of localized blockages

In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Vortex phase diagram of the kagome superconductor CsV 3 ⁢Sb 5

The screening response of the kagome superconductor CsV 3 ⁢Sb 5 was obtained from high-resolution ac mutual inductance measurements. At zero applied magnetic field and low temperatures, we observe no evidence for gapless quasiparticles, while near T c we find evidence for enhanced fluctuations. A rich vortex state appears above H c⁢1 ≈30 Oe, exhibiting successive emergence of vortex phases. For a fixed magnetic field, lowering the temperature below T c⁡ (H) leads to a wide range of vortex liquid state in a landscape of weak pinning potential, which gives rise to an irreversibility line. Further lowering the temperature, we identify the vortex melting line followed by the peak effect manifested in enhanced vortex pinning strength and critical current. Finally, we suggest that such an unusual behavior, where the peak effect region is fully contained within the vortex lattice state below the irreversibility line, is a consequence of the strong anisotropy and weak bulk pinning in CsV 3 ⁢Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurements of Magnetic Field Penetration of Materials for Superconducting Radiofrequency Cavities

Superconducting Radio Frequency (SRF) cavities used in particle accelerators are typically formed from or coated with superconducting materials. Currently high purity niobium is the material of choice for SRF cavities which have been optimized to operate near their theoretical field limits. This brings about the need for significant R&D efforts to develop next generation superconducting materials which could outperform Nb and keep up with the demands of new accelerator facilities. To achieve high quality factors and accelerating gradients, the cavity material should be able to remain in the superconducting Meissner state under high RF magnetic field without penetration of quantized magnetic vortices through the cavity wall. Therefore, the magnetic field at which vortices penetrate in a supercon- ductor is one of the key parameters of merit of SRF cavities. Techniques to measure the onset of magnetic field penetration on thin film samples need to be developed to mitigate the issues with the conventional magnetometry measurements which are strongly influenced by the film orientation and shape and edge effects. The applied magnetic field is also needed to be parallel to the one side of the superconductor to resemble the magnetic field profile at the surface of the SRF cavities operating at fundamental accelerating mode. In this work we report the development of an experimental setup called Magnetic Field Penetration (MFP) magnetometer to measure the field of full flux penetration through bulk, thin films and multilayered superconductors. Our system combines a small superconducting solenoid which can generate the magnetic field up to 500 mT at the sample surface and three Hall probes to detect the full flux penetration through the planner superconductor with 2 inch diameter. This setup was used to study alternative materials which could potentially outperform niobium, as well as SIS multilayer coatings on niobium.

Senevirathne, Iresha↗

Interfacial Spin-Orbit Coupling: A Platform for Superconducting Spintronics

Spin-orbit coupling (SOC) is a key interaction in spintronics, allowing electrical control of spin or magnetization and, vice versa, magnetic control of electrical current. However, recent advances have revealed much broader implications of SOC that is also central to the design of topological states with potential applications from low-energy dissipation and faster magnetization switching to high tolerance of disorder. SOC and the resulting emergent interfacial spin-orbit fields are simply realized in junctions through structural inversion asymmetry, while the anisotropy in magnetoresistance (MR) allows their experimental detection. Surprisingly, we demonstrate that an all-epitaxial ferromagnet/MgO/metal junction with a single ferromagnetic region and only negligible MR anisotropy undergoes a remarkable transformation below the superconducting transition temperature of the metal. The superconducting junction has a MR anisotropy 3 orders of magnitude higher and could enable novel applications in superconducting spintronics. In contrast to common realizations of MR effects that require a finite applied magnetic field, our system is designed to have two stable zero-field states with mutually orthogonal magnetizations: in plane and out of plane. This bistable magnetic anisotropy allows us to rule out orbital and vortex effects due to an applied magnetic field and identify the SOC origin of the observed MR. Such MR reaches approximately 20% without an applied magnetic field and could be further increased for large magnetic fields that support vortices. Our findings call for a revisit of the role of SOC, even when it seems negligible in the normal state, and suggest an alternative platform for superconducting spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Improved Subseasonal Forecasting of Extreme Polar Vortices Using Machine Learning

Our research was focused on forecasting the position and shape of the winter stratospheric polar vortex at a subseasonal timescale of 15 days in advance. To achieve this, we employed both statistical and neural network machine learning techniques. The analysis was performed on 42 winter seasons of reanalysis data provided by NASA giving us a total of 6,342 days of data. The state of the polar vortex for determined by using geometric moments to calculate the centroid latitude and the aspect ratio of an ellipse fit onto the vortex. Timeseries for thirty additional precursors were calculated to help improve the predictive capabilities of the algorithm. Feature importance of these precursors was performed using random forest to measure the predictive importance and the ideal number of precursors. Then, using the precursors identified as important, various statistical methods were tested for predictive accuracy with random forest and nearest neighbor performing the best. An echo state network, a type of recurrent neural network that features sparsely connected hidden layer and a reduced number of trainable parameters that allows for rapid training and testing, was also implemented for the forecasting problem. Hyperparameter tuning was performed for each methods using a subset of the training data. The algorithms were trained and tuned on the first 41 years of data, then tested for accuracy on the final year. In general, the centroid latitude of the polar vortex proved easier to predict than the aspect ratio across all algorithms. Random forest outperformed other statistical forecasting algorithms overall but struggled to predict extreme values. Forecasting from echo state network suggested a strong predictive capability past 15 days, but further work is required to fully realize the potential of recurrent neural network approaches.

54 ENVIRONMENTAL SCIENCES↗

The sound of a Martian dust devil

Dust devils (convective vortices loaded with dust) are common at the surface of Mars, particularly at Jezero crater, the landing site of the Perseverance rover. They are indicators of atmospheric turbulence and are an important lifting mechanism for the Martian dust cycle. Improving our understanding of dust lifting and atmospheric transport is key for accurate simulation of the dust cycle and for the prediction of dust storms, in addition to being important for future space exploration as grain impacts are implicated in the degradation of hardware on the surface of Mars. Here we describe the sound of a Martian dust devil as recorded by the SuperCam instrument on the Perseverance rover. The dust devil encounter was also simultaneously imaged by the Perseverance rover’s Navigation Camera and observed by several sensors in the Mars Environmental Dynamics Analyzer instrument. Combining these unique multi-sensorial data with modelling, we show that the dust devil was around 25 m large, at least 118 m tall, and passed directly over the rover travelling at approximately 5 m s –1 . Acoustic signals of grain impacts recorded during the vortex encounter provide quantitative information about the number density of particles in the vortex. The sound of a Martian dust devil was inaccessible until SuperCam microphone recordings. This chance dust devil encounter demonstrates the potential of acoustic data for resolving the rapid wind structure of the Martian atmosphere and for directly quantifying wind-blown grain fluxes on Mars.

54 ENVIRONMENTAL SCIENCES↗

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↗