Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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 1,045 records · Page 58

The stability of a compressible stratified shear layer

The stability of a shear layer under the effect of gravity is investigated using the compressible magnetohydrodynamic (MHD) equations, including an effective gravity term to represent the curvature effects of the flow and magnetic field line geometry. A general eigenmode equation is derived for a two-dimensional MHD fluid, and an energy-principle analysis to explain the effect of compressibility on the critical Richardson number is presented. For the case of a hyperbolic tangent shear flow and exponential density profile, it was found that, in the Boussinesq approximation, the compressibility raises the critical Richardson number from 1/4 to as much as 1/2, with the exact value depending on the value of the magnetic field at infinity. Under approximation of a strong asymptotic magnetic field, without invoking the Boussinesq approximation, it is shown both analytically and numerically that the density gradient terms cause the shear instability to be dispersive. The long-wavelength stability boundary for the Richardson number J = 0 is characterized by a normalized phase velocity c =

Wang, Z.↗

A 2-D oscillating flow analysis in Stirling engine heat exchangers

A two dimensional oscillating flow analysis was conducted, simulating the gas flow inside Stirling heat exchangers. Both laminar and turbulent oscillating pipe flow were investigated numerically for Re(max) = 1920 (Va = 80), 10800 (Va = 272), 19300 (Va = 272), and 60800 (Va = 126). The results are compared with experimental results of previous investigators. Also, predictions of the flow regime on present oscillating flow conditions were checked by comparing velocity amplitudes and phase differences with those from laminar theory and quasi-steady profile. A high Reynolds number k-epsilon turbulence model was used for turbulent oscillating pipe flow. Finally, performance evaluation of the K-epsilon model was made to explore the applicability of quasi-steady turbulent models to unsteady oscillating flow analysis.

Ahn, Kyung H.↗

A 2-D oscillating flow analysis in Stirling engine heat exchangers

A two-dimensional oscillating flow analysis was conducted, simulating the gas flow inside Stirling heat exchangers. Both laminar and turbulent oscillating pipe flow were investigated numerically for Re(max) = 1920 (Va = 80), 10800 (Va = 272), 19300 (Va = 272), and 60800 (Va = 126). The results are compared with experimental results of previous investigators. Also, predictions of the flow regime on present oscillating flow conditions were checked by comparing velocity amplitudes and phase differences with those from laminar theory and quasi-steady profile. A high Reynolds number k-epsilon turbulence model was used for turbulent oscillating pipe flow. Finally, performance evaluation of the K-epsilon model was made to explore the applicability of quasi-steady turbulent models to unsteady oscillating flow analysis.

Ahn, Kyung H.↗

The Effect of Strain Rate on the Evolution of Plane Wakes Subjected to Irrotational Strains

Direct numerical simulations of time-evolving turbulent plane wakes developing in the presence of irrotational plane strain applied at three different strain rates have been generated. The strain geometry is such that the flow is compressed in the streamwise direction and expanded in the cross-stream direction with the spanwise direction being unstrained. This geometry is the temporally evolving analogue of a spatially evolving wake in an adverse pressure gradient. A pseudospectral numerical method with up to 16 million modes is used to solve the equations in a reference frame moving with the irrotational strain. The initial condition for each simulation is taken from a previous turbulent self-similar plane wake direct numerical simulation at a velocity deficit Reynolds number, Re, of about 2,000. Although the evolutions of many statistics are nearly collapsed when plotted against total strain, there are some differences owing to the different strain rate histories. The impact of strain-rate on the wake spreading rate, the peak velocity deficit, the Reynolds stress profiles, and the flow structure is examined.

Rogers, Michael M.↗

Wear Mechanisms in Electron Sources for Ion Propulsion, 2: Discharge Hollow Cathode

The wear of the keeper electrode in discharge hollow cathodes is a major impediment to the implementation of ion propulsion onboard long-duration space science missions. The development of a predictive theoretical model for hollow cathode keeper life has long been sought, but its realization has been hindered by the complexities associated with the physics of the partially ionized gas and the associated erosion mechanisms in these devices. Thus, although several wear mechanisms have been hypothesized, a quantitative explanation of life test erosion profiles has remained incomplete. A two-dimensional model of the partially ionized gas in a discharge cathode has been developed and applied to understand the mechanisms that drove the erosion of the keeper in two long-duration life tests of a 30-cm ion thruster. An extensive set of comparisons between predictions by the numerical simulations and measurements of the plasma properties and of the erosion patterns is presented. It is found that the near-plume plasma oscillations, predicted by theory and observed by experiment, effectively enhance the resistivity of the plasma as well as the energy of ions striking the keeper.

hollow cathode↗

P161 Improved Impact of Atmospheric Infrared Sounder (AIRS) Radiance Assimilation in Numerical Weather Prediction

For over 6 years, AIRS radiances have been assimilated operationally into National (e.g. Environmental Modeling Center (EMC)) and International (e.g. European Centre for Medium-Range Weather Forecasts (ECMWF)), operational centers; assimilated in the North American Mesoscale (NAM) since 2008. Due partly to data latency and operational constraints, hyperspectral radiance assimilation has had less impact on the Gridpoint Statistical Interpolation (GSI) system used in the NAM and GFS. Objective of this project is to use AIRS retrieved profiles as a proxy for the AIRS radiances in situations where AIRS radiances are unable to be assimilated in the current operational system by evaluating location and magnitude of analysis increments.

Zavodsky, Bradley T.↗

NASA EOSDIS: Enabling Science by Improving User Knowledge

Lessons learned and impacts of applying these newer methods are explained and include several examples from our current efforts such as the interactive, on-line webinars focusing on data discovery and access including tool usage, informal and informative data chats with data experts across our EOSDIS community, data user profile interviews with scientists actively using EOSDIS data in their research, and improved conference and meeting interactions via EOSDIS data interactively used during hyper-wall talks and Worldview application. The suite of internet-based, interactive capabilities and technologies has allowed our project to expand our user community by making the data and applications from numerous Earth science missions more engaging, approachable and meaningful.

EOSDIS↗

Numerical Studies of Hall Thruster Acceleration Region Electron Transport

Azimuthal direction 1-D and azimuthal and axial direction 2-D Particle in Cell (PIC) simulations of the Hall thruster acceleration and near-plume regions are presented. The axial extent of the 2-D computational domain is based on the ExB drift parameter, 𝑫𝑬𝒙𝑩, the number of times around the channel an electron ExB drifts before being scattered. With only high energy main beam ions, the simulations produced a broad potential profile in the acceleration region. The addition of a small amount of ionization in the acceleration domain at a rate consistent with Hall thruster parameters causes the potential profile to develop two distinct regions. In the upstream portion of the computational domain, the potential profile is very steep, typical of Hall thruster acceleration regions. In the downstream portion, the potential profile is almost flat, like potentials in Hall thruster near plume regions. This bifurcation happens naturally; there is no other change in code parameters between the two results. Without the addition of a small amount of ionization in the acceleration region, the axial electron transport is very low. With ionization, the calculated ratio of electron current to ion current is the same order of those measured in Hall thrusters. This suggests that in regions with lower energy electrons and ions, plasma waves are effective in allowing crossfield transport, and axial electric fields are small. In the acceleration region, the plasma waves are less effective at scattering the energetic electrons, and high axial electric fields are needed to produce the observed cross-field electron transport.

Ortega, Alejandro Lopez↗

Simple models of the thermal structure of the Venusian ionosphere

Analytical and numerical models of plasma temperatures in the Venusian ionosphere are proposed. The magnitudes of plasma thermal parameters are calculated using thermal-structure data obtained by the Pioneer Venus Orbiter. The simple models are found to be in good agreement with the more detailed models of thermal balance. Daytime and nighttime temperature data along with corresponding temperature profiles are provided.

Whitten, R. C.↗

The stability to two-dimensional wakes and shear layers at high Mach numbers

This study is concerned with the stability properties of laminar free-shear-layer flows, and in particular symmetric two-dimensional wakes, for the supersonic through the hypersonic regimes. Emphasis is given to the use of proper wake profiles that satisfy the equations of motion at high Reynolds numbers. In particular the inviscid stability of a developing two-dimensional wake is studied as it accelerates at the trailing edge of a splitter plate. The nonparallelism of the flow is a leading-order effect in the calculation of the basic state, which is obtained numerically. Neutral stability characteristics are computed and the hypersonic stability is obtained by increasing the Mach number. It is found that the stability characteristics are altered significantly as the wake develops. Multiple modes (secondary modes) are found in the near wake that are closely related to the corresponding Blasius ones, but as the wake develops mode multiplicity is delayed to higher and higher Mach numbers. At a distance of about one plate length from the trailing edge, there is only one mode in a Mach number range of 0-20. The dominant mode emerging at all wake stations, and for high enough Mach numbers, is the so-called vorticity mode that is centered around the generalized inflection point layer. The structure of the dominant mode is also obtained analytically for all streamwise wake locations and it is shown how the far-wake limit is approached. Asymptotic results for the hypersonic mixing layer given by a tanh and a Lock distribution are also given.

Papageorgiou, Demetrios T.↗

Assessment of Radiation Scattering by Dust in Mars Entry

The impact of gas phase radiation scattering by dust particle presence in Mars entry scenarios is assessed using a suite of numerical tools. A hypersonic gas and dilute particle multiphase flow is simulated using a mixed continuum model and an efficient size distribution quadrature scheme. Radiative transport is simulated with a hybrid ray-tracing and P1 spherical harmonics model. The NASA HARA radiation code and a Lorentz-Mie scattering model are used to predict the gas-particle mixture radiative properties. Application of the numerical suite to a representative case based on the Mars 2020 entry shows that, even at conditions in excess of the strongest recorded dust storm, the augmentation in surface radiative heat flux due to scattering alone is lower than 1%, and therefore negligible. Predictions with even higher dust loading conditions confirm that even highly uncertain dust mass profiles are unlikely to be a concern for the scattering augmentation mechanism.

Andrew Hinkle↗

Numerical study of strong plasma shock waves produced in an electromagnetic shock tube.

A set of two-fluid Navier-Stokes equations with classical physical transport coefficients is used to compute the evolution and structure of collisional plasma shock waves in an electromagnetic shock tube. Shock speeds up to 200 cm/microsec and shocked-plasma temperatures of the order of kilovolts are studied. A strong transverse bias magnetic field is employed, which significantly alters the size and shape of the shock profiles, when compared with the zero bias field case. The wave structure is different from that of a two-fluid gasdynamic shock without transverse magnetic field. Especially significant is the effect of the small ion Larmor radius in reducing the gasdynamic shock thickness by at least an order of magnitude in the transverse shock case, which permits collisional shocks to have thicknesses much smaller than the post-shock mean free path. These collisional shock waves produce very hot ions, the ion temperature increasing with shock speed, although the formation time and distances also increase substantially with shock speed.

Schneider, S. H.↗

Photodissociation of O2 and H2O in the middle atmosphere: Comparison of numerical methods and impact on model O3 and OH

We have compared three photochemical diurnal models of O3 and OH in the upper stratosphere and mesosphere which use different techniques for calculating the absorption of solar ultraviolet radiation by the O2 Schumann-Runge bands. One model uses a detailed line-by-line representation of the O2 cross section from 1750-2050 A, while the two others use lower resolution, parameterized cross sections. Using the parameterized cross sections, the calculated O3 profiles for both day and night agree with those obtained from the line-by-line model to within 6%. This appears to eliminate inaccuracies in the parameterized O2 cross section as a major cause of previously reported model O3 deficits. A portion of the residual differences from the line-by-line model are attributed to inaccuracies in the calculated H2O photolysis rate. A parameterized H2O cross section is offered which improves the accuracy of this calculation.

Siskind, David E.↗

Shock-wave-induced turbulent boundary-layer separation at hypersonic speeds

An experiment is described that tests and guides computations of the interaction of a shock wave with a turbulent boundary layer. Numerical solutions of the time-averaged Navier-Stokes equations for the entire flow field employing algebraic eddy viscosity and turbulent Prandtl number models for shear stress and heat flux are presented and used to illustrate the dependence of the computations on the particulars of the turbulence models. To guide modifications in the models, the mean flow profiles and surface measurements of pressure, shear, and heat flux are analyzed critically. The results show that the models of eddy viscosity require substantial modifications in the interaction region. Improved solutions employing the experimentally modified models are presented.

Horstman, C. C.↗

Numerical analyses for treating diffusion in single-, two- and three-phase binary alloy systems

Numerical solutions were applicable for planar, cylindrical, or spherical geometries with any diffusion-zone size and any continuous variation of the diffusion coefficient with concentration. Special techniques were included in the analyses to account for differences in molal volumes, initiation and growth of an intermediate phase, disappearance of a phase, and the presence of an initial composition profile in the specimen. A major improvement in solution accuracy was achieved in the two phase analysis by employing a mass conservation criterion to establish the location of the interface rather than the conventional interface-flux-balance criterion. In the three phase analysis, computation time was minimized without sacrificing solution accuracy by treating the three phase problem as a two phase problem when the thickness of the intermediate phase was less than a preset small value. Three computer codes were developed to perform these analyses.

Tenney, D. R.↗

Refraction of microwave signals by water vapor

Tropospheric water vapor causes a refractive path length effect which is typically 5-10% of the 'dry' tropospheric effect and as large as several meters at elevation angles below 5 deg. The vertical water vapor profile is quite variable, and measurements of intensive atmospheric parameters such as temperature and humidity limited to the surface do not adequately predict the refractive effect. It is suggested that a water vapor refraction model that is a function of the amount of precipitable water alone can be successful at low elevation angles. From an extensive study of numerical ray tracings through radiosonde balloon data, such a model has been constructed. The model predicts the effect at all latitudes and elevation angles between 2 and 10 deg to an accuracy of better than 4% (11 cm at 3 deg elevation angle).

Goldfinger, A. D.↗

Algorithms and Libraries

This exploratory study initiated our inquiry into algorithms and applications that would benefit by latency tolerant approach to algorithm building, including the construction of new algorithms where appropriate. In a multithreaded execution, when a processor reaches a point where remote memory access is necessary, the request is sent out on the network and a context--switch occurs to a new thread of computation. This effectively masks a long and unpredictable latency due to remote loads, thereby providing tolerance to remote access latency. We began to develop standards to profile various algorithm and application parameters, such as the degree of parallelism, granularity, precision, instruction set mix, interprocessor communication, latency etc. These tools will continue to develop and evolve as the Information Power Grid environment matures. To provide a richer context for this research, the project also focused on issues of fault-tolerance and computation migration of numerical algorithms and software. During the initial phase we tried to increase our understanding of the bottlenecks in single processor performance. Our work began by developing an approach for the automatic generation and optimization of numerical software for processors with deep memory hierarchies and pipelined functional units. Based on the results we achieved in this study we are planning to study other architectures of interest, including development of cost models, and developing code generators appropriate to these architectures.

Dongarra, Jack↗

Lunar and Planetary Science XXXV: Meteorites

The session "Meteorites" included the following reports:Description of a New Stony Meteorite Find from Bulloch County, Georgia; Meteorite Ablation Derived from Cosmic Ray Track Data Dhofar 732: A Mg-rich Orthopyroxenitic Achondrite Halogens, Carbon and Sulfur in the Tagish Lake Meteorite: Implications for Classification and Terrestrial Alteration; Electromagnetic Scrape of Meteorites and Probably Columbia Tiles; Pre-Atmospheric Sizes and Orbits of Several Chondrites; Research of Shock-Thermal History of the Enstatite Chondrites by Track, Thermoluminescence and Neutron-Activation (NAA) Methods; Radiation and Shock-thermal History of the Kaidun CR2 Chondrite Glass Inclusions; On the Problem of Search for Super-Heavy Element Traces in the Meteorites: Probability of Their Discovery by Three-Prong Tracks due to Nuclear Spontaneous Fission Trace Element Abundances in Separated Phases of Pesyanoe, Enstatite Achondrite; Evaluation of Cooling Rate Calculated by Diffusional Modification of Chemical Zoning: Different Initial Profiles for Diffusion Calculation; Mineralogical Features and REE Distribution in Ortho- and Clinopyroxenes of the HaH 317 Enstatite Chondrite Dhofar 311, 730 and 731: New Lunar Meteorites from Oman; The Deuterium Content of Individual Murchison Amino Acids; Clues to the Formation of PV1, an Enigmatic Carbon-rich Chondritic Clast from the Plainview H-Chondrite Regolith Breccia ;Numerical Simulations of the Production of Extinct Radionuclides and ProtoCAIs by Magnetic Flaring.

Source record↗