SEARCH · Search NASA
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.
A numerical procedure for empirically determining a 'continuum' level in low-dispersion spectra from the IUE is described. The technique was developed for the automated fitting of SWP-LO (1150 A - 2000 A) spectra of stars of late spectral type, but it is more widely applicable.
The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expression for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations. The numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the Courant-Friedrichs-Lewy (CFL) limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about six - eight grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 - 300 time steps. Three numerical examples of waves in compressible flow are included: (1) sound propagation in a duct with linear shear, (2) linear wave growth in a compressible free shear layer, and (3) vortex pairing in a compressible free shear layer perturbed at two frequencies.
The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expressions for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations, and the numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the CFL limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about 6 to 8 grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 to 300 time steps. Three numerical examples of waves in compressible flow are included.
It is well known that due to vortex shedding in high speed flow over cutouts, cavities, and gaps, intense noise may be generated. Strong tonal oscillations occur in a feedback cycle in which the vortices shed from the upstream edge of the cavity convect downstream and impinge on the cavity lip, generating acoustic waves that propagate upstream to excite new vortices. Numerical simulation of such a complicated process requires a scheme that can: (1) resolve acoustic waves with low dispersion and numerical dissipation, (2) handle nonlinear and discontinuous waves (e.g. shocks), and (3) have an effective (near field) nonreflecting boundary condition (NRBC). The new space time conservation element and solution element method, or CE/SE for short, is a numerical method that meets the above requirements.
Superalloys can be strengthened by oxide dispersions to potentially improve their mechanical properties at high temperatures. Additive manufacturing has opened new paths for introducing oxide particles into superalloys with increased processing flexibility, to potentially provide oxide dispersion strengthening. Numerically controlled electrical discharge machining (EDM) has often been used as a machining process for initial sectioning, detailed roughing, and even finishing operations of conventional superalloys. The objective of this study was to examine how additively manufactured superalloys containing oxide dispersions respond when sectioned using EDM. Consistent sections of several superalloys, two containing an oxide dispersion, were sliced using a typical EDM machine with standard brass wire and a typical roughing condition set. Roughness, microstructure, and residual stress at the EDM surfaces were compared for the varied superalloys.
The growth of electrostatic waves near the plasma frequency (fp) due to an unstable electron beam is investigated by solving the unmagnetized electrostatic dispersion equation numerically. The numerical solutions are compared with analytic theories for reactive (or fluid-type) and kinetic versions of the beam instability, and for the O'Neil/Malmberg connection of the beam and Langmuir modes. Conditions for growth significantly above or below fp are given. Three general results are found: (1) The unstable waves do not grow on a mode with Langmuir dispersion except in the limit of a very dilute beam with growth on O'Neil/Malmberg's connected mode; (2) The properties of the unstable mode depend strongly on beam parameters such as beam density, speed and temperature; and (3) The frequency of maximum growth frequently lies significantly above or below fp, and differs significantly from that predicted by the Langmuir dispersion relation. Results imply important consequences for theories of strong turbulence and nonlinear wave-wave processes, and observational identification of fp from observed wave frequencies.
The several wake vortices which originate from aircraft wingtips, flap edges, engine pylons, etc. usually merge, in the far field, to form a single pair whose structure determines the hazard posed to encountering aircraft. To gain an understanding of the process whereby vortices merge and disperse, a numerical study was made of the interaction of two-dimensional, time-dependent, inviscid vortical regions. It was found that discrete boundaries, which depend on the structure and spacing of the vortices, distinguish merging from nonmerging situations. Furthermore, certain arrays of finite vortex cores that alternate in sign were found to undergo division and merging that may be useful in alleviating the hazard posed by aircraft wakes.
A description is given of a number of numerical schemes to solve an evolution equation (Korteweg-deVries) that arises when modelling the propagation of water waves in a channel. The discussion also includes the results of numerical experiments made with each of the schemes. It is suggested, on the basis of these experiments, that one of the schemes may have (discrete) solitary-wave solutions.
Two main areas were addressed in support of an effort to understand mechanism responsible for the broadband electrostatic noise (BEN) observed in the magnetotail. The first area concerns the generation of BEN in the boundary layer region of the magnetotail whereas the second area concerns the occassional presence of BEN in the neutral sheet region. For the generation of BEN in the boundary layer region, a hybrid simulation code was developed to perform reliable longtime, quiet, highly resolved simulations of field aligned electron and ion beam flow. The result of the simulation shows that broadband emissions cannot be generated by beam-plasma instability if realistic values of the ion beam parameters are used. The waves generated from beam-plasma instability are highly discrete and are of high frequencies. For the plasma sheet boundary layer condition, the wave frequencies are in the kHz range, which is incompatible with the observation that the peak power in BEN occur in the 10's of Hz range. It was found that the BEN characteristics are more consistent with lower hybrid drift instability. For the occasional presence of BEN in the neutral sheet region, a linear analysis of the kinetic cross-field streaming instability appropriate to the neutral sheet condition just prior to onset of substorm expansion was performed. By solving numerically the dispersion relation, it was found that the instability has a growth time comparable to the onset time scale of substorm onset. The excited waves have a mixed polarization in the lower hybrid frequency range. The imposed drift driving the instability corresponds to unmagnetized ions undergoing current sheet acceleration in the presence of a cross-tail electric field. The required electric field strength is in the 10 mV/m range which is well within the observed electric field values detected in the neutral sheet during substorms. This finding can potentially account for the disruption of cross-tail current and its diversion to the ionosphere to form the substorm current wedge. Furthermore, a number of features associated with substorm expansion onset can be understood based on this substorm onset scenario.
The Landau damping of electrostatic Langmuir waves and ion-acoustic waves in a hot, isotropic, nonmagnetized, generalized Lorentzian plasma is analyzed using the modified plasma dispersion function. Numerical solutions for the real and imaginary parts of the wave frequency omega sub 0 - (i)(gamma) have been obtained as a function of the normalized wave number (k)(lambda sub D), where lambda sub D is the electron Debye length. For both particle distributions the electrostatic modes are found to be strongly damped at short wavelengths. At long wavelengths, this damping becomes less severe, but the attenuation of Langmuir waves is much stronger for a generalized Lorentzian plasma than for a Maxwellian plasma. It is concluded that Landau damping of ion-acoustic waves is only slightly affected by the presence of a high energy tail, but is strongly dependent on the ion temperature.
The results of a spatial stability analysis of a two-dimensional slab jet, in which optically thin radiative cooling is dynamically important, are presented. We study both magnetized and unmagnetized jets at external Mach numbers of 5 and 20. We model the cooling rate by using two different cooling curves: one appropriate to interstellar gas, and the other to photoionized gas of reduced metallicity. Thus, our results will be applicable to both protostellar (Herbig-Haro) jets and optical jets from active galactic nuclei. We present analytical solutions to the dispersion relations in useful limits and solve the dispersion relations numerically over a broad range of perturbation frequencies. We find that the growth rates and wavelengths of the unstable Kelvin-Helmholtz (K-H) modes are significantly different from the adiabatic limit, and that the form of the cooling function strongly affects the results. In particular, if the cooling curve is a steep function of temperature in the neighborhood of the equilibrium state, then the growth of K-H modes is reduced relative to the adiabatic jet. On the other hand, if the cooling curve is a shallow function of temperature, then the growth of K-H modes can be enhanced relative to the adiabatic jet by the increase in cooling relative to heating in overdense regions. Inclusion of a dynamically important magnetic field does not strongly modify the important differences between an adiabatic jet and a cooling jet, provided the jet is highly supermagnetosonic and not magnetic pressure-dominated. In the latter case, the unstable modes behave more like the transmagnetosonic magnetic pressure-dominated adiabatic limit. We also plot fluid displacement surfaces associated with the various waves in a cooling jet in order to predict the structures that might arise in the nonlinear regime. This analysis predicts that low-frequency surface waves and the lowest order body modes will be the most effective at producing observable features in the jet.
Between 27 June and 14 July 2019 aerosol layers were observed by the United Kingdom (UK) Raman lidar network in the upper troposphere and lower stratosphere. The arrival of these aerosol layers in late June caused some concern within the London Volcanic Ash Advisory Centre (VAAC) as according to dispersion simulations the volcanic plume from the 21 June 2019 eruption of Raikoke was not expected over the UK until early July. Using dispersion simulations from the Met Office Numerical Atmospheric-dispersion Modelling Environment (NAME), and supporting evidence from satellite and in situ aircraft observations, we show that the early arrival of the stratospheric layers was not due to aerosols from the explosive eruption of the Raikoke volcano but due to biomass burning smoke aerosols associated with intense forest fires in Alberta, Canada, that occurred 4 d prior to the Raikoke eruption. We use the observations and model simulations to describe the dispersion of both the volcanic and forest fire aerosol clouds and estimate that the initial Raikoke ash aerosol cloud contained around 15 Tg of volcanic ash and that the forest fires produced around 0.2 Tg of biomass burning aerosol. The operational monitoring of volcanic aerosol clouds is a vital capability in terms of aviation safety and the synergy of NAME dispersion simulations, and lidar data with depolarising capabilities allowed scientists at the Met Office to interpret the various aerosol layers over the UK and attribute the material to their sources. The use of NAME allowed the identification of the observed stratospheric layers that reached the UK on 27 June as biomass burning aerosol, characterised by a particle linear depolarisation ratio of 9 %, whereas with the lidar alone the latter could have been identified as the early arrival of a volcanic ash–sulfate mixed aerosol cloud. In the case under study, given the low concentration estimates, the exact identification of the aerosol layers would have made little substantive difference to the decision-making process within the London VAAC. However, our work shows how the use of dispersion modelling together with multiple observation sources enabled us to create a more complete description of atmospheric aerosol loading.
Results of a numerical investigation of the dispersion of solid particles in decaying isotropic turbulence are presented. The 3D time-dependent velocity field of a homogeneous nonstationary turbulence is computed using the method of direct numerical simulation (DNS). The dispersion characteristics of three different solid particles (corn, copper, and glass) injected in the flow are obtained by integrating the complete equation of particle motion along the instantaneous trajectories of 22-cubed particles for each particle type, and then by performing ensemble averaging. Good agreement was achieved between the present DNS results and the measured time development of the mean-square displacement of the particles. Questions of how and why the dispersion statistics of a solid particle differ from those of its corresponding fluid point and surrounding fluid and what influences inertia and gravity have on these statistics are also discussed.
Management of off-nominal situations on-board the International Space Station (ISS) is important to its continuous operation. One situation of concern is an accidental release of a chemical into the ISS atmosphere. In particular, introduction of ammonia into the cabin atmosphere can occur via the interface heat exchangers (IFHX) between the external thermal control system containing ammonia and internal thermal control system that uses water as a coolant to remove heat from ISS subsystems. Breach of the water/ammonia barrier of the IFHX can lead to a catastrophic rupture. Once the liquid water/ammonia mixture exits the ITCS, it instantly vaporizes and mixes with the U.S. Laboratory cabin atmosphere that results in rapid contamination of the cabin. The goal of the study is to assess the amount of ammonia in the Russian Segment by the time the crew is able to isolate the U.S. Segment. A Computational Fluid Dynamics (CFD) model for an accurate prediction of airflow and ammonia transport in the frozen flow field within the assembly complete ISS cabin was developed. The localized effects of ammonia dispersion are examined and discussed.
Release of ammonia into the International Space Station (ISS) cabin atmosphere can occur if the water/ammonia barrier breach of the active thermal control system (ATCS) interface heat exchanger (IFHX) happens. After IFHX breach liquid ammonia is introduced into the water-filled internal thermal control system (ITCS) and then to the cabin environment through a ruptured gas trap. Once the liquid water/ammonia mixture exits ITCS, it instantly vaporizes and mixes with the U.S. Laboratory cabin air that results in rapid deterioration of the cabin conditions. The goal of the study is to assess ammonia propagation in the Station after IFHX breach to plan the operation procedure. A Computational Fluid Dynamics (CFD) model for accurate prediction of airflow and ammonia transport within each of the modules in the ISS cabin was developed. CFD data on ammonia content in the cabin aisle way of the ISS and, in particular, in the Russian On- Orbit Segment during the period of 15 minutes after gas trap rupture are presented for four scenarios of rupture response. Localized effects of ammonia dispersion and risk mitigation are discussed.
A theoretical procedure to determine the heat transfer characteristics of the interline region (junction of liquid-solid-vapor) from the macroscopic optical and thermophysical properties of the system is outlined. The analysis is based on the premise that the interline transport processes are controlled by the London-van der Waals dispersion force between condensed phases (solid and liquid). Numerical values of the dispersion constant are presented. The procedure is used to compare the relative size of the interline heat sink of various systems using a constant heat flux mode. This solution demonstrates the importance of the interline heat flow number, which is evaluated for various systems.
In this paper, we study numerically the excitation of plasma waves by electron beams, in the auroral region above 2000 km of altitude. We have solved the fully kinetic dispersion relation, using numerical method and found the real frequency and the growth rate of the plasma wave modes. We have examined the instability properties of low-frequency waves such as the Electromagnetic Ion Cyclotron (EMIC) wave as well as Lower-Hybrid (LH) wave in the range of high-frequency. In all cases, the source of free energy are electron beams propagating parallel to the geomagnetic field. We present some features of the growth rate modes, when the cold plasma parameters are changed, such as background electrons and ions species (H(+) and O(+)) temperature, density or the electron beam density and/or drift velocity. These results can be used in a test-particle simulation code, to investigate the ion acceleration and their implication in the auroral acceleration processes, by wave-particle interaction.
The problem 2 in Category 3 of the 4th Computational Aeroacoustic(CAA) Workshop is solved using the space-time conservation element and solution element (CE/SE) method. This problem models rotor-stator interaction in a 2D cascade. It involves complex geometries and flow physics including vortex shedding and acoustic radiation. The parallel version of the 2D nonlinear Euler solver is used with an unstructured triangular mesh to solve this problem. The Giles approach is incorporated with the CE/SE method to handle non-equal pitches of the rotor and stator. Validation on the Giles approach is performed using Problem 3.1 in the 2nd CAA Workshop. The space-time CE/SE method is a finite volume method with second-order accuracy in both space and time. The flux conservation is enforced in both space and time instead of space only. It has low numerical dissipation and dispersion errors. It uses simple non-reflecting boundary conditions and is compatible with unstructured meshes. It is simple, flexible, and generate reasonably accurate solutions. The CE/SE method has been successfully applied to solve numerous practical problems, especially aeroacoustic problems. Some preliminary numerical results of the benchmark problem 3.2 of the 4th CAA Workshop are shown. The steady-state pressure contour is plotted. The mean pressure distribution on the blade surface is compared with Turbo solution showing a good agreement. The sound pressure level versus the rotor harmonic n at the six designated positions on the blade surface, three locations at inlet plane, and three locations at the outlet plane are plotted. It can be seen that the acoustic response exists only at the excitation frequencies (n = 1,2,3). On the blade surface, the acoustic wave at n = 1 is dominant, while at the inlet and outlet planes, the sound pressure level at n = 2 becomes the largest, which is similar to the results presented. The distribution of sound pressure level at different spatial modes along the z- direction is plotted for n = 1,2,3, respectively. It shows that the spatial modes m = -32 and 22 at n = 1 exponentially decay, and the spatial modes m = 10 at n = 2, m = -42 and 12 at n = 3 propagate both upstream and downstream, which agrees with the prediction based on the linearized theory. Some oscillations are observed, which needs to be investigated further. In the final paper, the numerical results will be compared with a frequency-domain solver LINFLUX solution if it is available.