Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “higher order beam dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 37 records · Page 2

Development of a composite tailoring procedure for airplane wing

The development of a composite wing box section using a higher order-theory is proposed for accurate and efficient estimation of both static and dynamic responses. The theory includes the effect of through-the-thickness transverse shear deformations which is important in laminated composites and is ignored in the classical approach. The box beam analysis is integrated with an aeroelastic analysis to investigate the effect of composite tailoring using a formal design optimization technique. A hybrid optimization procedure is proposed for addressing both continuous and discrete design variables.

Chattopadhyay, Aditi↗

Higher Order Multipole Analysis of Beam Crabbing Mode at 197 MHz for Electron-Ion Collider

In the local crabbing scheme for Electron Ion Collider interaction region, beam rotation is initiated by one set of crab cavities at a high beta region with near 90 deg phase advance to the interaction point (IP), and ceased by an identical set of cavities on the symmetrically opposite side. The fundamental mode of all crab cavities is designed to apply a momentum kick to the particle bunches to complete this half cycle of rotation in the horizontal plane, meanwhile the electromagnetic (EM) field can be expanded into a series of harmonics. Despite the lowest dipole mode, higher harmonics cause instabilities during beam operation and loss of particles. We used SimTrack, which is the 6-d symplectic element-by-element tracking code developed at Brookhaven National Laboratory in 2015, to study the threshold of each higher-order harmonic multipolar mode with the criteria of maintaining a dynamic aperture of 6 σ or greater for the Hadron Storage Ring (HSR). In this paper, we showed the threshold for each higher harmonic multipole and their comparison with the multipole strength for the current crab cavity design.

43 PARTICLE ACCELERATORS↗

Transverse BBU Suppression with Feedback for Energy Recovery Linacs

In this study, we present beam dynamics simulations of the Strong Hadron Cooling (SHC) Energy Recovery Linac (ERL) that incorporate a feedback (FB) system designed to mitigate Beam Breakup (BBU) instabilities. The FB effectively suppresses BBU and increases the threshold current by roughly an order of magnitude. Our analysis revealed the aliasing behavior of the dominant higher-order mode (HOM) and its suppression by the FB, as evident from the power spectral density (PSD) of the beam centroid’s horizontal motion at the exit of the final linac section. Although the FB system— even when affected by noise—successfully reduced the BBU peak, it also amplified the PSD at other frequencies. These results provide important insights into BBU suppression and threshold current enhancement in ERLs.

Setiniyaz, S. [Thomas Jefferson National Accelerat↗

Evaluation of several modal methods for transient structural analysis

The convergence properties of three modal methods for solving structural transient response (STR) problems are examined by comparing the accuracy of the total number of modes used to represent the responses. The goal is to identify a method which involves as few degrees of freedom as possible. Two mode superposition methods, the mode displacement method (MDM) and the mode acceleration method (MAM), are fitted with equivalent expressions for damped response. An error norm is employed to compare the accuracies of the methods in calculating displacements, moments, and shear forces for a viscously-damped cantilevered beam experiencing several dynamic loading conditions and levels of damping. The MAM converges quicker and uses fewer dof, mainly because it has a pseudo-static response term which accounts for some of the flexibility of higher modes. A third method, which integrates the convolution integral iteratively over time, generates a higher order approximation to the STR with a reduced number of modes. It converges fastest among the three, but requires definition of an appropriate forcing function. Finally, the MAM and MDM models are applied to study the effects of the forcing frequency and the initial conditions.

Camarda, Charles J.↗

Enhanced beam-beam modeling to include longitudinal variation during weak-strong simulation

Beam-beam interactions pose substantial challenges in the design and operation of circular colliders, significantly affecting their performance. In particular, the weak-strong simulation approach is pivotal for investigating single-particle dynamics during the collider design phase. This paper evaluates the limitations of existing models in weak-strong simulations, noting that while they accurately account for energy changes due to slingshot effects, they fail to incorporate longitudinal coordinate changes ( z variation). To address this gap, we introduce two novel transformations that enhance Hirata’s original framework by including both z variation and slingshot effect-induced energy changes. Through rigorous mathematical analysis and extensive weak-strong simulation studies, we validate the efficacy of these enhancements in achieving a more precise simulation of beam-beam interactions. Our results reveal that although z variation constitutes a higher-order effect and does not substantially affect the emittance growth rate within the specific design parameters of the Electron-Ion Collider, the refined model offers improved accuracy, particularly in scenarios involving the interaction between beam-beam effects and other random diffusion processes, as well as in simulations incorporating realistic lattice models. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A new pre-loaded beam geometric stiffness matrix with full rigid body capabilities

Space structures, such as the Space Station solar arrays, must be extremely light-weight, flexible structures. Accurate prediction of the natural frequencies and mode shapes is essential for determining the structural adequacy of components, and designing a controls system. The tension pre-load in the 'blanket' of photovoltaic solar collectors, and the free/free boundary conditions of a structure in space, causes serious reservations on the use of standard finite element techniques of solution. In particular, a phenomenon known as 'grounding', or false stiffening, of the stiffness matrix occurs during rigid body rotation. The authors have previously shown that the grounding phenomenon is caused by a lack of rigid body rotational capability, and is typical in beam geometric stiffness matrices formulated by others, including those which contain higher order effects. The cause of the problem was identified as the force imbalance inherent in the formulations. In this paper, the authors develop a beam geometric stiffness matrix for a directed force problem, and show that the resultant global stiffness matrix contains complete rigid body mode capabilities, and performs very well in the diagonalization methodology customarily used in dynamic analysis.

Bosela, P. A.↗

Very-High Dynamic Range, 10,000 Frames/Second Pixel Array Detector for Electron Microscopy

Precision and accuracy of quantitative scanning transmission electron microscopy (STEM) methods such as ptychography, and the mapping of electric, magnetic, and strain fields depend on the dose. Reasonable acquisition time requires high beam current and the ability to quantitatively detect both large and minute changes in signal. A new hybrid pixel array detector (PAD), the second-generation Electron Microscope Pixel Array Detector (EMPAD-G2), addresses this challenge by advancing the technology of a previous generation PAD, the EMPAD. The EMPAD-G2 images continuously at a frame-rates up to 10 kHz with a dynamic range that spans from low-noise detection of single electrons to electron beam currents exceeding 180 pA per pixel, even at electron energies of 300 keV. The EMPAD-G2 enables rapid collection of high-quality STEM data that simultaneously contain full diffraction information from unsaturated bright-field disks to usable Kikuchi bands and higher-order Laue zones. Test results from 80 to 300 keV are presented, as are first experimental results demonstrating ptychographic reconstructions, strain and polarization maps. We introduce a new information metric, the maximum usable imaging speed (MUIS), to identify when a detector becomes electron-starved, saturated or its pixel count is mismatched with the beam current.

36 MATERIALS SCIENCE↗

Beam Dynamics of the Strong Hadron Cooler ERL at the Electron-Ion Collider

The Strong Hadron Cooler Energy Recovery Linac (SHC-ERL) for the Electron-Ion Collider (EIC) is required to produce an electron beam with a bunch charge of 1 nC, an average current of 100 mA, and a beam energy of both 150 and 55 MeV, with strict requirements for the transverse emittance, slice energy spread, and longitudinal distribution. This paper covers the design in detail, including preliminary considerations of higher order effects and diagnostics.

Deitrick, K.↗

Structural dynamic model obtained from flight use with piloted simulation and handling qualities analysis

The ability to use flight data to determine an aircraft model with structural dynamic effects suitable for piloted simulation. and handling qualities analysis has been developed. This technique was demonstrated using SR-71 flight test data. For the SR-71 aircraft, the most significant structural response is the longitudinal first-bending mode. This mode was modeled as a second-order system, and the other higher order modes were modeled as a time delay. The distribution of the modal response at various fuselage locations was developed using a uniform beam solution, which can be calibrated using flight data. This approach was compared to the mode shape obtained from the ground vibration test, and the general form of the uniform beam solution was found to be a good representation of the mode shape in the areas of interest. To calibrate the solution, pitch-rate and normal-acceleration instrumentation is required for at least two locations. With the resulting structural model incorporated into the simulation, a good representation of the flight characteristics was provided for handling qualities analysis and piloted simulation.

Powers, Bruce G.↗

Third integer resonant extraction transit time simulation studies

In this work, we present the investigation of transit time of particles in the non-linear third-integer resonant extraction process. Transit time is defined as the number of turns a particle takes to get extracted once it is in the unstable region in the phase space, i.e., outside the triangular separatrix in case of third-integer resonance. The study of transit time is important because transit time directly contributes to the beam response time during resonant extraction and thus knowing it apriori would be practically useful in designing of the extraction system. In this work, we shall investigate the analytical derivation of the transit time of particles (to the first order Kobayashi Hamiltonian) in different parts of the phase space distribution and compare against the analytical results. We also compare the simulation result of the transit time of particles (with higher statistics) for the static as well as dynamic extraction conditions cases, particularly in the context of resonant extraction parameters for Mu2e experiment at Fermilab.

Narayanan, Aakaash [Fermilab]↗

Third Integer Resonant Extraction Transit Time Simulation Studies

In this work, we present the investigation of transit time of particles in the non-linear third-integer resonant extraction process. Transit time is defined as the number of turns a particle takes to get extracted once it is in the unstable region in the phase space, i.e., outside the triangular separatrix in case of third-integer resonance. The study of transit time is important because transit time directly contributes to the beam response time during resonant extraction and thus knowing it apriori would be practically useful in designing of the extraction system. In this work, we shall investigate the analytical derivation of the transit time of particles (to the first order Kobayashi Hamiltonian) in different parts of the phase space distribution and compare against the analytical results. We also compare the simulation result of the transit time of particles (with higher statistics) for the static as well as dynamic extraction conditions cases, particularly in the context of resonant extraction parameters for Mu2e experiment at Fermilab.

Narayanan, Aakaash [Fermilab] (ORCID:0000000157944↗

Derivation of a Tappered p-Version Beam Finite Element

A tapered p-version beam finite element suitable for dynamic applications is derived. The taper in the element is represented by allowing the area moments of inertia to vary as quartic polynomials along the length of the beam, and the cross-sectional area to vary as a quadratic polynomial. The p-version finite-element characteristics are implemented through a set of polynomial shape functions. The lower-order shape functions are identical to the classical cubic and linear shape functions normally associated with a beam element. The higher-order shape functions are a hierarchical set of polynomials that are integrals of orthogonal polynomials. Explicit expressions for the mass and stiffness matrices are presented for an arbitrary value of p. The element has been verified to be numerically stable using shape functions through 22nd order.

Hinnant, Howard E.↗

The 2022 magneto-optics roadmap

Abstract Magneto-optical (MO) effects, viz. magnetically induced changes in light intensity or polarization upon reflection from or transmission through a magnetic sample, were discovered over a century and a half ago. Initially they played a crucially relevant role in unveiling the fundamentals of electromagnetism and quantum mechanics. A more broad-based relevance and wide-spread use of MO methods, however, remained quite limited until the 1960s due to a lack of suitable, reliable and easy-to-operate light sources. The advent of Laser technology and the availability of other novel light sources led to an enormous expansion of MO measurement techniques and applications that continues to this day (see section 1). The here-assembled roadmap article is intended to provide a meaningful survey over many of the most relevant recent developments, advances, and emerging research directions in a rather condensed form, so that readers can easily access a significant overview about this very dynamic research field. While light source technology and other experimental developments were crucial in the establishment of today’s magneto-optics, progress also relies on an ever-increasing theoretical understanding of MO effects from a quantum mechanical perspective (see section 2), as well as using electromagnetic theory and modelling approaches (see section 3) to enable quantitatively reliable predictions for ever more complex materials, metamaterials, and device geometries. The latest advances in established MO methodologies and especially the utilization of the MO Kerr effect (MOKE) are presented in sections 4 (MOKE spectroscopy), 5 (higher order MOKE effects), 6 (MOKE microscopy), 8 (high sensitivity MOKE), 9 (generalized MO ellipsometry), and 20 (Cotton–Mouton effect in two-dimensional materials). In addition, MO effects are now being investigated and utilized in spectral ranges, to which they originally seemed completely foreign, as those of synchrotron radiation x-rays (see section 14 on three-dimensional magnetic characterization and section 16 on light beams carrying orbital angular momentum) and, very recently, the terahertz (THz) regime (see section 18 on THz MOKE and section 19 on THz ellipsometry for electron paramagnetic resonance detection). Magneto-optics also demonstrates its strength in a unique way when combined with femtosecond laser pulses (see section 10 on ultrafast MOKE and section 15 on magneto-optics using x-ray free electron lasers), facilitating the very active field of time-resolved MO spectroscopy that enables investigations of phenomena like spin relaxation of non-equilibrium photoexcited carriers, transient modifications of ferromagnetic order, and photo-induced dynamic phase transitions, to name a few. Recent progress in nanoscience and nanotechnology, which is intimately linked to the achieved impressive ability to reliably fabricate materials and functional structures at the nanoscale, now enables the exploitation of strongly enhanced MO effects induced by light–matter interaction at the nanoscale (see section 12 on magnetoplasmonics and section 13 on MO metasurfaces). MO effects are also at the very heart of powerful magnetic characterization techniques like Brillouin light scattering and time-resolved pump-probe measurements for the study of spin waves (see section 7), their interactions with acoustic waves (see section 11), and ultra-sensitive magnetic field sensing applications based on nitrogen-vacancy centres in diamond (see section 17). Despite our best attempt to represent the field of magneto-optics accurately and do justice to all its novel developments and its diversity, the research area is so extensive and active that there remains great latitude in deciding what to include in an article of this sort, which in turn means that some areas might not be adequately represented here. However, we feel that the 20 sections that form this 2022 magneto-optics roadmap article, each written by experts in the field and addressing a specific subject on only two pages, provide an accurate snapshot of where this research field stands today. Correspondingly, it should act as a valuable reference point and guideline for emerging research directions in modern magneto-optics, as well as illustrate the directions this research field might take in the foreseeable future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Rarefied gas dynamics; Proceedings of the Eighth International Symposium, Stanford University, Stanford, Calif., July 10-14, 1972

Recent studies in the theory of rarefied gases (kinetic theory, shock structure, molecular forces and interactions, evaporation and condensation, expansions and jets, flow past bodies, gas-solid interactions) are presented along with descriptions of numerous experiments in these areas. Some of the topics covered include: singular limits of solutions of Boltzmann's equations, the Boltzmann and Klimontovich formalisms with reference to turbulent resolution, higher-order kinetic theory for the Lorentz gas, intermolecular potentials from shock structure experiments, cooling of particles in aerosol beams, hypersonic-subsonic transition in spherically symmetric expanding flows, absolute velocity measurements in a rarefied gas flow by an ion time-of-flight technique, the product distribution of gas-solid reactions from absolute rate theory, and relaxation of a gas with internal degrees of freedom. Individual items are announced in this issue.

Karamcheti, K.↗

Space Charge Simulations of High Intensity Proton Beams in the AGS Booster

Computer simulation studies have been performed to understand the beam behavior and to explore intensity limitations of proton beams in the AGS Booster at higher beam intensities. During the 100 GeV polarized proton operations of RHIC Run 2024, sPHENIX operated in modes with a crossing angle at collisions in order to mitigate beam-beam effects. Three different running modes were employed: (a) sPHENIX operated with a negative (-2 mrad) crossing angle, and STAR operated with 0 mrad. Both experiments were brought into collisions at the start of the store. (b) sPHENIX was brought into collisions with 0 mrad first. Then STAR was brought into collisions after the beam-beam parameter from sPHENIX reduced to below $10 \times 10^{-3}$. (c) sPHENIX operated with a positive (+1.5 mrad) crossing angle, and STAR operated with 0 mrad. Both experiments were brought into collisions at the start of the store. The collisions with a crossing angle of up to $\pm 2$ mrad, as in running modes (a) and (c), lead to large Piwinski angle in the new sPHENIX detector, which reduces luminosity if other parameters are unchanged. There are two ways to compensate the reduction in luminosity: squeeze $\beta^{*}$ if there is sufficient dynamic aperture, or increase the injected beam intensity. The first part of polarized proton operation during RHIC Run 2024 was dedicated to increasing the intensity. Different configurations were tested with crossing angle and lattice adjustments on RHIC. At the same time, new injector configurations were developed and tested in an effort to push for both higher intensity and better quality of the beam injected into RHIC. When the beam intensity is increased, space charge is a concern particularly in the lower energy stages of acceleration, such as during the injection and the early part of the Booster cycle, which could become a dominant effect in limiting the intensity of the beam that can be delivered to RHIC.

43 PARTICLE ACCELERATORS↗

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗

Novel deuterium–tritium (DT) neutron activation diagnostics for increased neutron yield operations at Sandia’s Z-facility (invited)

As inertial confinement fusion experiments at Sandia National Laboratories’ Z-facility begin using deuterium–tritium (DT) fuels, the expected DT-neutron yields are expected to increase many orders of magnitude. This brings about the need for diagnostics that have a significantly higher dynamic range of operation and ability to obtain yields where there will be significant radiation holds. Using the current activation diagnostic infrastructure at the Z-facility, a pair of activation diagnostics are being developed. Here, the first is fluorine activation, which will be used in a similar manner to current copper activation diagnostics that are fielded in a rabbit for quick removal from the experiments without opening the vacuum experiment containment. The second is a prompt DT-neutron detector that exploits the activation of oxygen with similar operation characteristics as the current Be-probes fielded at Z. We will discuss the design of this pair of instruments, along with calibration work performed at Sandia’s Ion Beam Laboratory in support of the development and anticipated instrument performance.

Mangan, Michael Antony↗

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗