Results for “Kinetic simulations” — Page 18

NASA NTRS1988

The effect of the near earth micrometeoroid environment on a highly reflective mirror surface

A resurgence of interest in placing large solar concentrator solar dynamic systems in space for power generation has brought up again a concern for maintaining the integrity of the optical properties of highly specular reflecting surfaces in the near earth space environment. One of the environmental hazards needing evaluation is the micrometeoroid environment. It has been shown that highly reflective polished metals and thin film coatings degrade when exposed to simulated micrometeoroids in the lab. At NASA-Lewis, a shock tube was used to simulate the phenomenon of micrometeoroid impact by accelerating micron sized particles to hypervelocities. Any changes in the optical properties of surfaces exposed to this impact were then evaluated. The degradation of optical properties of polished metals and thin metallic films after exposure to simulated micrometeoroids was determined as a function of impacting kinetic energy area of the particles. A calibrated sensor was developed to not only detect the micrometeoroid environment, but also to evaluate the degradation of the optical properties of thin aluminum films in space. Results of the simulation are presented and discussed.

Mirtich, Michael J.↗

NASA NTRS1988

The effect of the near earth micrometeoroid environment on a highly reflective mirror surface

A resurgence of interest in placing large solar concentrator solar dynamic systems in space for power generation has brought up again a concern for maintaining the integrity of the optical properties of highly specular reflecting surfaces in the near earth space environment. One of the environmental hazards needing evaluation is the micrometeoroid environment. It has been shown that highly reflective polished metals and thin film coatings degrade when exposed to simulated micrometeoroids in the lab. At NASA-Lewis, a shock tube was used to simulate the phenomenon of micrometeoroid impact by accelerating micron sized particles to hypervelocities. Any changes in the optical properties of surfaces exposed to this impact were then evaluated. The degradation of optical properties of polished metals and thin metallic films after exposure to simulated micrometeoroids was determined as a function of impacting kinetic energy area of the particles. A calibrated sensor was developed to not only detect the micrometeoroid environment, but also to evaluate the degradation of the optical properties of thin aluminum films in space. Results of the simulation are presented and discussed.

Mirtich, Michael J.↗

NASA NTRS1995

Bounce-Averaged Hamiltonian for Charged Particles in an Axisymmetric but Nondipolar Model Magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M, J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. The L value of a field line is (by definition) inversely proportional to the flux enclosed by the corresponding magnetic shell of equatorial radius r(sub 0), and the L value at the neutral line (r(sub 0) = b) is denoted L*. The azimuthal coordinate phi measures magnetic local time. The best functional representation found for the normalized difference (L(exp 3)a(exp 3)/mu(sub E))(B(sub m) - B(sub 0)) between mirror-point field B(sub m) and equatorial field B(sub 0) along any field line is a 5-term expansion in powers (2/3 through 6/3) of the quantity X equivalent to (La/mu(sub E))(exp 1/2)K, where K equivalent to (J(exp 2)/8m(sub 0)M)(exp 1/2) is an adiabatically conserved quantity independent of particle energy, m(sub 0) is the rest mass of the particle, and a is the radius of the Earth. This functional form is motivated by results for limiting cases in which particles mirror very near and very far from the magnetic equator. Expansion coefficients corresponding to various powers of X are obtained from least squares fits to numerically computed results for X as a function of L and B(sub m). These are accurately expressible as fourth-order polynomials in (r(sub 0)/b)(exp 3), hence indirectly as functions of L/L* = 3La/2b. This representation, which leads (except for a manageably small region of parameter space) to better than 1% accuracy in the specification of B(sub m) as a function of K and L, allows bounce-averaged guiding center simulations to be performed without actually tracing the bounce motions of individual particles. Bounce-averaged drifts L' (meridional) and phi' (azimuthal) are proportional to derivatives of the Hamiltonian H (sum of kinetic and potential energies) with respect to phi and L, respectively. Our formulation thus provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗

NASA NTRS1984

Chemical kinetic reaction mechanism for the combustion of propane

A detailed chemical kinetic reaction mechanism for the combustion of propane is presented and discussed. The mechanism consists of 27 chemical species and 83 elementary chemical reactions. Ignition and combustion data as determined in shock tube studies were used to evaluate the mechanism. Numerical simulation of the shock tube experiments showed that the kinetic behavior predicted by the mechanism for stoichiometric mixtures is in good agrement with the experimental results over the entire temperature range examined (1150-2600K). Sensitivity and theoretical studies carried out using the mechanism revealed that hydrocarbon reactions which are involved in the formation of the HO2 radical and the H2O2 molecule are very important in the mechanism and that the observed nonlinear behavior of ignition delay time with decreasing temperature can be interpreted in terms of the increased importance of the HO2 and H2O2 reactions at the lower temperatures.

Jachimowski, C. J.↗

NASA NTRS2002

Fundamentals of Plasma Sails Propulsion Concept

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.

Khazanov, G. V.↗

NASA NTRS1990

Filtering analysis of a direct numerical simulation of the turbulent Rayleigh-Benard problem

A filtering analysis of a turbulent flow was developed which provides details of the path of the kinetic energy of the flow from its creation via thermal production to its dissipation. A low-pass spatial filter is used to split the velocity and the temperature field into a filtered component (composed mainly of scales larger than a specific size, nominally the filter width) and a fluctuation component (scales smaller than a specific size). Variables derived from these fields can fall into one of the above two ranges or be composed of a mixture of scales dominated by scales near the specific size. The filter is used to split the kinetic energy equation into three equations corresponding to the three scale ranges described above. The data from a direct simulation of the Rayleigh-Benard problem for conditions where the flow is turbulent are used to calculate the individual terms in the three kinetic energy equations. This is done for a range of filter widths. These results are used to study the spatial location and the scale range of the thermal energy production, the cascading of kinetic energy, the diffusion of kinetic energy, and the energy dissipation. These results are used to evaluate two subgrid models typically used in large-eddy simulations of turbulence. Subgrid models attempt to model the energy below the filter width that is removed by a low-pass filter.

Eidson, T. M.↗

NASA NTRSDate not supplied

Integrated Process-Structure-Property Simulations for Additive Manufacturing Using the Open-Source Materialite Package

The microstructure and properties of additively manufactured (AM) metals are strongly dependent on process conditions. Therefore, process-structure-property (PSP) simulations are a useful tool for exploring process parameter space, studying process variations, and quantifying uncertainty in material properties. However, integrating process-structure and structure-property simulations often involves connecting multiple software packages. Each package may use unique data structures and require substantial domain knowledge. This presentation demonstrates PSP simulation capabilities of Materialite, an open-source package developed at NASA Langley Research Center. Materialite simplifies model linkages by using a common data structure and model interface, enabling straightforward simulation across a PSP model chain. Physics-based models, including kinetic Monte Carlo and crystal plasticity, are implemented within the package. The model interface is also intended to simplify implementation of new models and enable integration with external simulation tools. Example use cases include uncertainty quantification with PSP models and GPU-accelerated powder bed fusion AM process models.

additive manufacturing↗

NASA NTRSDate not supplied

A Green’s Function Sensor Fusion Approach for Evaluating Spacecraft Entry Heating From on-Board Thermal Instrumentation

During atmospheric entry, distributed thermal measurements are critical to enable the evaluation of heat loads on spacecraft thermal protection systems (TPS). In recent space exploration missions, Schmidt-Boelter-type heat flux gauges have been integrated into the TPS alongside conventional temperature measurement instrumentation to measure total (convective and radiative) and radiative heat transfer rates1,2. While direct heat flux sensors (HFS) are able to provide valuable information detailing the thermal loads experienced by spacecraft, the interpretation of these measurements in unsteady, convective environments requires a correction factor to account for local heating augmentations at the cold wall HFS surface2,3. Current efforts to estimate cold wall correction factors, and thus recover the hot wall TPS heat flux, rely on time-marching computational fluid dynamics (CFD) simulations2. Simulation-based methods are susceptible to large uncertainties, however, as they require estimations of vehicle trajectory, gas kinetics, wall catalysis models, and other flight conditions as input parameters2,4. Furthermore, CFD simulations are computationally expensive and cannot efficiently survey all possible entry scenarios, exacerbating the uncertainty of reconstructed hot wall heat flux values. These drawbacks motivate the development of alternative hot wall heat flux reconstruction methods that are not reliant on CFD-based correction factors.

Kenneth McAfee↗

NASA NTRSDate not supplied

A Green’s Function Sensor Fusion Approach for Evaluating Spacecraft Entry Heating From on-Board Thermal Instrumentation

During atmospheric entry, distributed thermal measurements are critical to enable the evaluation of heat loads on spacecraft thermal protection systems (TPS). In recent space exploration missions, Schmidt-Boelter-type heat flux gauges have been integrated into the TPS alongside conventional temperature measurement instrumentation to measure total (convective and radiative) and radiative heat transfer rates1,2. While direct heat flux sensors (HFS) are able to provide valuable information detailing the thermal loads experienced by spacecraft, the interpretation of these measurements in unsteady, convective environments requires a correction factor to account for local heating augmentations at the cold wall HFS surface2,3. Current efforts to estimate cold wall correction factors, and thus recover the hot wall TPS heat flux, rely on time-marching computational fluid dynamics (CFD) simulations2. Simulation-based methods are susceptible to large uncertainties, however, as they require estimations of vehicle trajectory, gas kinetics, wall catalysis models, and other flight conditions as input parameters2,4. Furthermore, CFD simulations are computationally expensive and cannot efficiently survey all possible entry scenarios, exacerbating the uncertainty of reconstructed hot wall heat flux values. These drawbacks motivate the development of alternative hot wall heat flux reconstruction methods that are not reliant on CFD-based correction factors.

Kenneth McAfee↗

NASA NTRSDate not supplied

A Green’s Function Sensor Fusion Approach for Evaluating Spacecraft Entry Heating From on-Board Thermal Instrumentation

During atmospheric entry, distributed thermal measurements are critical to enable the evaluation of heat loads on spacecraft thermal protection systems (TPS). In recent space exploration missions, Schmidt-Boelter-type heat flux gauges have been integrated into the TPS alongside conventional temperature measurement instrumentation to measure total (convective and radiative) and radiative heat transfer rates1,2. While direct heat flux sensors (HFS) are able to provide valuable information detailing the thermal loads experienced by spacecraft, the interpretation of these measurements in unsteady, convective environments requires a correction factor to account for local heating augmentations at the cold wall HFS surface2,3. Current efforts to estimate cold wall correction factors, and thus recover the hot wall TPS heat flux, rely on time-marching computational fluid dynamics (CFD) simulations2. Simulation-based methods are susceptible to large uncertainties, however, as they require estimations of vehicle trajectory, gas kinetics, wall catalysis models, and other flight conditions as input parameters2,4. Furthermore, CFD simulations are computationally expensive and cannot efficiently survey all possible entry scenarios, exacerbating the uncertainty of reconstructed hot wall heat flux values. These drawbacks motivate the development of alternative hot wall heat flux reconstruction methods that are not reliant on CFD-based correction factors.

Kenneth McAfee↗

NASA NTRS1991

Atmospheric pressure flow reactor: Gas phase chemical kinetics under tropospheric conditions without wall effects

A flow reactor for simulating the interaction in the troposphere is set forth. A first reactant mixed with a carrier gas is delivered from a pump and flows through a duct having louvers therein. The louvers straighten out the flow, reduce turbulence and provide laminar flow discharge from the duct. A second reactant delivered from a source through a pump is input into the flowing stream, the second reactant being diffused through a plurality of small diffusion tubes to avoid disturbing the laminar flow. The commingled first and second reactants in the carrier gas are then directed along an elongated duct where the walls are spaced away from the flow of reactants to avoid wall interference, disturbance or turbulence arising from the walls. A probe connected with a measuring device can be inserted through various sampling ports in the second duct to complete measurements of the first and second reactants and the product of their reaction at selected XYZ locations relative to the flowing system.

Koontz, Steven L.↗

NASA NTRS2000

A Study of Wake Development and Structure in Constant Pressure Gradients

Motivated by the application to high-lift aerodynamics for commercial transport aircraft, a systematic investigation into the response of symmetric/asymmetric planar turbulent wake development to constant adverse, zero, and favorable pressure gradients has been conducted. The experiments are performed at a Reynolds number of 2.4 million based on the chord of the wake generator. A unique feature of this wake study is that the pressure gradients imposed on the wake flow field are held constant. The experimental measurements involve both conventional LDV and hot wire flow field surveys of mean and turbulent quantities including the turbulent kinetic energy budget. In addition, similarity analysis and numerical simulation have also been conducted for this wake study. A focus of the research has been to isolate the effects of both pressure gradient and initial wake asymmetry on the wake development. Experimental results reveal that the pressure gradient has a tremendous influence on the wake development, despite the relatively modest pressure gradients imposed. For a given pressure gradient, the development of an initially asymmetric wake is different from the initially symmetric wake. An explicit similarity solution for the shape parameters of the symmetric wake is obtained and agrees with the experimental results. The turbulent kinetic energy budget measurements of the symmetric wake demonstrate that except for the convection term, the imposed pressure gradient does not change the fundamental flow physics of turbulent kinetic energy transport. Based on the turbulent kinetic energy budget measurements, an approach to correct the bias error associated with the notoriously difficult dissipation estimate is proposed and validated through the comparison of the experimental estimate with a direct numerical simulation result.

Thomas, Flint O.↗

NASA NTRS1992

Compressible homogeneous shear: Simulation and modeling

Compressibility effects were studied on turbulence by direct numerical simulation of homogeneous shear flow. A primary observation is that the growth of the turbulent kinetic energy decreases with increasing turbulent Mach number. The sinks provided by compressible dissipation and the pressure dilatation, along with reduced Reynolds shear stress, are shown to contribute to the reduced growth of kinetic energy. Models are proposed for these dilatational terms and verified by direct comparison with the simulations. The differences between the incompressible and compressible fields are brought out by the examination of spectra, statistical moments, and structure of the rate of strain tensor.

Sarkar, S.↗

NASA NTRS1991

Compressible homogeneous shear - Simulation and modeling

Compressibility effects were studied on turbulence by direct numerical simulation of homogeneous shear flow. A primary observation is that the growth of the turbulent kinetic energy decreases with increasing turbulent Mach number. The sinks provided by compressible dissipation and the pressure dilatation, along with reduced Reynolds shear stress, are shown to contribute to the reduced growth of kinetic energy. Models are proposed for these dilatational terms and verified by direct comparison with the simulations. The differences between the incompressible and compressible fields are brought out by the examination of spectra, statistical moments, and structure of the rate of strain tensor.

Sarkar, S.↗

NASA NTRS1988

The seasonal cycle of energetics from the GLAS/UMD climate GCM

The annual cycle of atmospheric energetics from a 2-year integration of the GLAS/UMD Climate GCM is computed and compared to results from the European Centre analyses of the GWE year, and to previously published results on a global basis. All calculations are done in the mixed space-time domain. The main conclusions are: (1) the seasonal cycle of today's eddy kinetic energy (in both hemispheres), and of the transient eddy available potential energy and the potential-to-kinetic energy conversions (mean and eddy) in the Northern Hemisphere are well simulated by the GCM; (2) the GCM's tendency to have anomalously large mean u-winds at upper levels in high latitudes leads to excessive wintertime values of mean kinetic and available potential energies, and causes distortions in the GCM latitude-height distribution of kinetic energy and of many of the conversions; (3) the eddy conversion of available potential-to-kinetic energy obtained from the ageostrophic wind in these analyses; and (4) the conversions in the Southern Hemisphere are not well simulated by the GCM, although the observations are somewhat questionable.

Straus, David M.↗

NASA NTRS1994

The asymptotic state of rotating homogeneous turbulence at high Reynolds numbers

The long-time, asymptotic state of rotating homogeneous turbulence at high Reynolds numbers has been examined using large-eddy simulation of the incompressible Navier-Stokes equations. The simulations were carried out using 128 x 128 x 512 collocation points in a computational domain that is four times longer along the rotation axis than in the other directions. Subgrid-scale motions in the simulations were parameterized using a spectral eddy viscosity modified for system rotation. Simulation results show that in the asymptotic state the turbulence kinetic energy undergoes a power-law decay with an exponent which is independent of rotation rate, depending only on the low-wavenumber form of the initial energy spectrum. Integral lengthscale growth in the simulations is also characterized by power-law growth; the correlation length of transverse velocities exhibiting much more rapid growth than observed in non-rotating turbulence.

Squires, Kyle D.↗

NASA NTRS1995

The making of an Alfvenic fluctuation: The resolution of a second-order analysis

Ulysses observations of the high speed polar streams show that they are largely occupied by very large amplitude Alfvenic fluctuations accompanied by many rotational discontinuities. These fluctuations have a nearly constant magnetic intensity or amplitude, and the magnetic field direction per wave cycle sweeps only through a limited arc, much as a car wiperblade would do. Barnes and Hollweg (JGR, 79, 2302, 1974) suggested that this unusual waveform could arise from an obliquely propagating and linearly polarized Alfven wave of finite amplitude. From a second-order analysis, they showed that the existence of a particular solution with a constant amplitude but could not resolve the outcome of the homogeneous solution which consisted of fast waves. They suggested that Landau damping of these fast waves may be needed to get the observed waveform. We present a 1 1/2 D hybrid simulation which is fully nonlinear and correctly describes the ion kinetics for an initially monochromatic and linearly polarized Alfven wave propagating obliquely to the background magnetic field. The wave has a large amplitude and a wavelength so long that it can be considered dispersionless for simulation times. At early times, the second harmonic in density and in magnetic field transverse to the initial wave magnetic field are generated and have more power than other harmonics. Steepening is observed with a weak fast shock emerging, but no rotational discontinuity is left behind, and instead a constant amplitude and an arc-shaped waveform is made. The compressional component which develops after the shocks have dissipated is to zeroth order better described as a pure acoustic wave than as a fast wave. This might be explained by the relaxing of the Alfven wave to a state where its ponderomotive force vanishes so that the compressional component can travel almost independently of it.

Vasquez, Bernard J.↗

NASA NTRS1987

Martin Marietta atomic oxygen Low Earth Orbit (LEO) simulation

An atomic oxygen beam apparatus that produces flux levels and atomic kinetic energy similar to that encountered by spacecraft in low Earth orbit (LEO) is described. The beam apparatus consists of an electric discharge ion source, mass filter, decelerator, and neutralizer. Specific design goals include a 1.3 cm beam diameter, a 5 eV beam energy, and a flux density on the order of 10(15) cm(-2) s(-1). The total fluence will be on the order of 10(19) cm(-2) for an 8 hour test. The neutral oxygen beam will expose various materials contained within a large target chamber. Within the chamber will be a rather complex suite of instrumentation that will allow real-time studies of material mass loss and reactant species spatial distribution. In addition, a UV solar simulator will aid in the understanding of various synergistic effects.

Sjolander, Gary W.↗