Engineering PapersSearch

SEARCH · Engineering Papers

Results for “plume 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 19 records

Python Based Plume Dynamics Estimation Tool (PyPDET) Rapid Plume Strike Analysis for RPOD Maneuvers in Deep Space Operations

I worked as a NASA Intern during the Summer 2023 term in the DS-00 division under the supervision of my mentor, Dr. Jonathan Pitt. Our goal was to build on our previous work from 2022 to develop a plume strike estimation tool using a prescribed physics methodology and model plume impingement effects while considering the dynamics of a rendezvous, operations, proximity, and docking (RPOD) maneuver. This tool supports previously configured CFD-DSMC calculations by allowing for rapid analysis of initial designs using a low-fidelity source flow model. Engineers can then use the high-fidelity CFD-DSMC tool to consolidate results as they work towards finalizing a design. This year’s project was focused on developing a software application that other engineers would be using in their analysis. Thus, the user’s experience was considered in the development of this application. Proper documentation, testability, and modularity of the codebase was our priority. For example, the project included auto documentation procedures to start building towards a User Manual, while also including dedicated demonstration cases for more explicit communication of functionality. Also, this project included a framework for testing the source code for future developments. Additionally, care was taken to develop the code using an Object-Oriented Programming approach. Thus, allowing for a modular extensibility of functionality in anticipation of future developments. The core work of this project was developing an algorithm that would transform the visiting vehicle and associated thruster data according to the kinematics described in the jet firing history. It would then calculate the estimated plume strikes on two of the target vehicles and write data accordingly into a VTK file. Summer work is to conclude by developing and presenting a PowerPoint slide deck at the intern exit briefing on August 11 th , 2023. Once the model for simple plume strike calculations is developed and tested there are several avenues to explore to continue development of this tool. These are also discussed in this report.

Plume Impingement

The onset of plume dynamics in the spectral space

Time-dependent convection for infinite Prandtl number fluids has been investigated hitherto in the physical domain. Dynamics of plumelike structures resulting from boundary-layer instabilities can be interpreted in the spectral domain as having both direct and inverse cascades of energy taking place. In the physical domain only a small part of the spectrum is discernible. But the rest of the spectrum is needed for properly describing the nonlinear process. Only a relatively few modes, fewer than ten, are actually required for describing the essential features associated with the onset of plume dynamics.

Vincent, Alain P.

Chesapeake Bay plume dynamics from LANDSAT

LANDSAT images with enhancement and density slicing show that the Chesapeake Bay plume usually frequents the Virginia coast south of the Bay mouth. Southwestern (compared to northern) winds spread the plume easterly over a large area. Ebb tide images (compared to flood tide images) show a more dispersed plume. Flooding waters produce high turbidity levels over the shallow northern portion of the Bay mouth.

Munday, J. C., Jr.

Energy balance and plume dynamics in Triton's lower atmosphere

The present study of the thermal balance-affecting relationships among Triton lower atmosphere thermal conduction, eddy mixing, condensation, and radiative heating indicates that, while the temperature gradient is negative in the lower atmosphere, it becomes positive at higher altitudes due to the downward conduction of ionospheric heat. This temperature profile is essentially consistent with radio-occultation experiment data; the geyser-like plumes observed by Voyager suggest that the Trioton atmosphere's convective and conductive regions join near 10-km altitude, and that the values inferred for the eddy diffusion and heat-transport coefficients indicate a profile reminiscent of the earth's.

Yelle, Roger V.

Space Station flexible dynamics under plume impingement

Assembly of the Space Station requires numerous construction flights by the Space Shuttle. A particularly challenging problem is that of control of each intermediate station configuration when the shuttle orbiter is approaching it to deliver the next component. The necessary braking maneuvers cause orbiter thruster plumes to impinge on the station, especially its solar arrays. This in turn causes both overall attitude errors and excitation of flexible-body vibration modes. These plume loads are predicted to lead to CMG saturation during the approach of the orbiter to the SC-5 station configuration, necessitating the use of the station RCS jets for desaturation. They are also expected to lead to significant excitation of solar array vibrations. It is therefore of great practical importance to investigate the effects of plume loads on the flexible dynamics of station configuration SC-5 as accurately as possible. However, this system possesses a great many flexible modes (89 below 5 rad/s), making analysis time-consuming and complicated. Model reduction techniques can be used to overcome this problem, reducing the system model to one which retains only the significant dynamics, i.e. those which are strongly excited by the control inputs or plume disturbance forces and which strongly couple with the measured outputs. The particular technique to be used in this study is the subsystem balancing approach which was previously developed by the present investigator. This method is very efficient computationally. Furthermore, it gives accurate results even for the difficult case where the structure has many closed-spaced natural frequencies, when standard modal truncation can give misleading results. Station configuration SC-5 is a good example of such a structure.

Williams, Trevor

Fiber-Based Multi-Resolution Imaging and Multidimensional X-ray Diagnostics for Ejecta Dynamics in Plume Surface Interactions

The interaction of a rocket exhaust plume with a particulate-laden surface creates a complex, multiphase flow field that can destabilize the vehicle and damage nearby equipment. This study investigates two methods for capturing the particulate dynamics of plume-surface interaction (PSI): optical fiber-based multi-resolution Mie scattering and multi-dimensional X-ray radiography. Mie scattering was used to track PSI-interacted particles, trace their paths, and measure velocities. While effective for the jet periphery and early PSI stages, the technique becomes limited as the scattering cross-section increases over time due to particle displacement from the soil bed, causing the core to become optically dense and appear as a luminous, opaque region. To address this, X-ray radiography was explored as a complementary method for visualizing the optically dense core. PSI experiments were conducted with both reacting and non-reacting jets to evaluate these approaches across a range of optical and flow parameters. The results demonstrated the capability of the fiber-based multi-resolution imaging system to capture simultaneous fields of view at varying magnifications (1x, 2x, 4x), and the ability of X-ray imaging to penetrate the optically dense plume, revealing flow structures that would otherwise be obscured in scattering-based methods. Data were collected for various PSI parameters, including three different heights above the surface, to analyze the ejecta properties and the plume’s temporal evolution. These results provide the first imaging strategy capable of resolving flow structures over a wide spatial dynamic range while also offering the first visualization of the optically dense core.

Ejecta

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.

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.

Greendyke, R. B.

Mantle plumes - A boundary layer approach for Newtonian and non-Newtonian temperature-dependent rheologies

Stress is placed on the temperature dependence of both a linear Newtonian rheology and a nonlinear olivine rheology in accounting for narrow mantle flow structures. The boundary-layer theory developed incorporates an arbitrary temperature-dependent power-law rheology for the medium, in order to facilitate the study of mantle plume dynamics under real conditions. Thermal, kinematic, and dynamic structures of mantle plumes are modelled by a two-dimensional natural-convection boundary layer rising in a fluid with a temperature-dependent power-law relationship between shear stress and strain rate. An analytic similarity solution is arrived at for upwelling adjacent to a vertical isothermal stress-free plane. Newtonian creep as a deformation mechanism, thermal anomalies resulting from chemical heterogeneity, the behavior of plumes in non-Newtonian (olivine) mantles, and differences in the dynamics of wet and dry olivine are discussed.

Yuen, D. A.

Quantification of Asian Dust Plume Seasonal Dynamics and Regional Features

Dust is but one of many aerosols that are analyzed at the Jet Propulsion Laboratory in Pasadena. The purpose of this paper is to describe the process in analyzing and digitizing dust within a source region to better explain the work achieved by my internship. This paper will go over how to view collected data by Multi-angle Imaging SpectroRadiometer (MISR) [1] and the procedure of downloading data to be analyzed. With this data, one can digitize dust plumes using two methods called plume lines and plume polygons with the help of the software MISR INteractive eXplorer (MINX)[3]; thus, the theory of MINX's[3] algorithm and these methods are discussed in detail. Research was gathered from these techniques and emphasis is also focused on the obtained data and results.

aerosols

Aerosol Source Plume Physical Characteristics from Space-based Multiangle Imaging

Models that assess aerosol effects on regional air quality and global climate parameterize aerosol sources in terms of amount, type, and injection height. The multiangle imaging spectroradiometer (MISR) aboard NASA's Terra satellite retrieves total column aerosol optical thickness (AOT), and aerosol type over cloud-free land and water. A stereo-matching algorithm automatically retrieves reflecting-layer altitude wherever clouds or aerosol plumes have discernable spatial contrast, with about 500-m accuracy, at 1.1-km horizontal resolution. Near-source biomass burning smoke, volcanic effluent, and desert dust plumes are observed routinely, providing information about aerosol amount, particle type, and injection height useful for modeling applications. Compared to background aerosols, the plumes sampled have higher AOT, contain particles having expected differences in Angstrom exponent, size, single-scattering albedo, and for volcanic plume and dust cloud cases, particle shape. As basic thermodynamics predicts, thin aerosol plumes lifted only by regional winds or less intense heat sources are confined to the boundary layer. However, when sources have sufficient buoyancy, the representative plumes studied tend to concentrate within discrete, high-elevation layers of local stability; the aerosol is not uniformly distributed up to a peak altitude, as is sometimes assumed in modeling. MISR-derived plume heights, along with meteorological profile data from other sources, make it possible to relate radiant energy flux observed by the moderate resolution imaging spectroradiometer (MODIS), also aboard the Terra spacecraft, to convective heat flux that plays a major role in buoyant plume dynamics. A MISR climatology of plume behavior based on these results is being developed.

aerosol plumes

Tectonic plates, D (double prime) thermal structure, and the nature of mantle plumes

It is proposed that subducting tectonic plates can affect the nature of thermal mantle plumes by determining the temperature drop across a plume source layer. The temperature drop affects source layer stability and the morphology of plumes emitted from it. Numerical models are presented to demonstrate how introduction of platelike behavior in a convecting temperature dependent medium, driven by a combination of internal and basal heating, can increase the temperature drop across the lower boundary layer. The temperature drop increases dramatically following introduction of platelike behavior due to formation of a cold temperature inversion above the lower boundary layer. This thermal inversion, induced by deposition of upper boundary layer material to the system base, decays in time, but the temperature drop across the lower boundary layer always remains considerably higher than in models lacking platelike behavior. On the basis of model-inferred boundary layer temperature drops and previous studies of plume dynamics, we argue that generally accepted notions as to the nature of mantle plumes on Earth may hinge on the presence of plates. The implication for Mars and Venus, planets apparently lacking plate tectonics, is that mantle plumes of these planets may differ morphologically from those of Earth. A corollary model-based argument is that as a result of slab-induced thermal inversions above the core mantle boundary the lower most mantle may be subadiabatic, on average (in space and time), if major plate reorganization timescales are less than those acquired to diffuse newly deposited slab material.

Lenardic, A.

Simulation of Vacuum Chamber Pressure Distribution with Surrogate Modeling and Uncertainty Quantification

A major challenge in understanding differences in electric propulsion performance in ground tests and in space operations concerns the pressure distribution within the test vacuum chamber. The chamber backpressure is much higher than experienced in space, modifying thruster performance and plume dynamics. Numerical simulation is a key element to determining the background conditions in non-ideal vacuum chamber environments. An important parameter for the accurate simulation of chamber backpressure is the probability that an atom will stick to a cryogenic panel used to pump away the plume gases. This quantity can be used to model vacuum pumps in particle-based kinetic numerical methods. In this work, a three-dimensional direct simulation Monte Carlo code is used to model neutral xenon atoms flowing from the anode of the H9 Hall Effect Thruster within the University of Michigan’s Large Vacuum Test Facility. Simulated pressures are compared with ion gauge pressure measurements to infer the effective sticking coefficient of the chamber’s vacuum pumps. A pressure predicting surrogate model is developed for inference of pump sticking coefficients and for uncertainty quantification. This information enables accurate and useful kinetic simulations of electric propulsion thruster plasma plumes in vacuum chambers.

DSMC

MISR Interactive Explorer (MINX) : Production Digitizing to Retrieve Smoke Plume Heights and Validating Heights Against Lidar Data

The height at which smoke from a wildfire is injected into the atmosphere is an important parameter for climatology, because it determines how far the smoke can be transported. Using the MINX program to analyze MISR (Multi-angle Imaging Spectro-Radiometer) data, I digitized wildfire smoke plumes to add to an existing database of these heights for use by scientists studying smoke transport and plume dynamics. In addition to using MINX to do production digitizing of heights, I assisted in gathering lidar data for an ongoing validation of MINX and helped evaluate those data.

smoke plume

Dynamics of volcanic plumes on Io

Ballistic and aerodynamic models are proposed to explain the dynamics of the volcanic plumes on Io. A ballistic model, which assumes that eruptive gas and solid-particle flows become completely and permanently decoupled over a small distance as compared with the maximum particle height, is shown to be unable to account for the observed shape of Plume 3. An aerodynamic model which assumes the entrainment and transportation of solid particles by the gas flow into the central portion of the top of the plume, with ballistic particle trajectories only after passage through the first part of their descent, is then developed. The model accounts for plume shape by the formation of a weak decompression shock near the edge of the plume stalk, which constrains horizontal expansion until a sufficient altitude is reached. Calculated plume dimensions are then compared with observations of Ionian Plumes 1 and 3.

Cook, A. F.

Numerical modeling of the hybrid plume plasma rocket

This paper reports on ongoing research in the area of numerical simulations of hybrid plume dynamics, and describes the basic mechanisms for energy and particle transport. The solution approach assumes cylindrical geometry and includes a multiplicity of atomic reactions, as well as the presence of a strong magnetic field. The principal reactions are electron impact ionization, and resonant charge exchange between hot ions and cold neutrals. Radial particle and energy transport is mainly by diffusion. The present description includes a modified Bohm diffusion model for plasma in the core of the plume, and classical neutral particle diffusion in the cooler regions of the flow. Neutrals are allowed to 'free stream' in the low density regions, where the collision mean-free-path becomes comparable, or even larger than the characteristic dimensions of the system. A parameteric scan of the various regimes of interest has been completed and the most significant results are presented.

Chang, F. R.