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,549 records · Page 86

Electrical conductivity and velocity of highly ionized plasma flows - Theory and experiment.

Use of an immersible, three-coil, magnetic-induction probe, previously tested in a low-density supersonic argon jet, to measure electrical conductivity and velocity profiles of a highly ionized high-density nitrogen jet in the continuum flow regime where effects due to probe bow shocks and boundary layers might not be negligible. Measured centerline values of electrical conductivity and velocity were compared with predictions based on a theoretical analysis previously developed to study the gas as it expanded adiabatically and inviscidly from an equilibrium sonic state to the nozzle exit. The resulting numerical exit plane values for electron density and electron temperature were then substituted into the Spitzer-Haerm conductivity formula to compute a theoretical conductivity which agreed within 40% of the measured conductivity, while the calculated and experimental velocity values differed by as much as 50%. The lack of agreement was attributed to the possible use of invalid assumptions and boundary conditions in the computer analysis or to the unknown effects of shocks on the probe data.

Vendell, E. W.↗

Using PCDS to study the influence of UV flux variation on the middle atmosphere

The value of the Pilot Climate Data System (PCDS) in examining the influences of short-period solar flux variations on the middle atmosphere was demonstrated. Several Nimbus satellite data sets proposed for the study exist in the PCDS. Planned for retrieval through the PCDS are ozone mixing ratios and cumulative ozone profiles available from the Backscatter Ultraviolet Spectrophotometer on Nimbus-4. Also to be accessed are Nimbus-7 data sets that will provide ozone and nitrogen dioxide mixing ratios, temperature profiles from the Limb Infrared Monitor of the Stratosphere, and ozone mixing ratio profiles from the Solar Backscatter Ultraviolet Spectrophotometer. Bypassing the time-consuming process of reading raw data tapes, the researchers plan to transfer the PCDS processed data to an IBM PC at Iowa State University. The IBM PC will serve as an intermediate vehicle for transferring the data to the NOAA CYBER 840 at Boulder, Colorado, where research will continue on both an eight-layer radiative-photochemical numerical model and on a nonlinear dynamical model. in this research are the FGGE and ERB data sets. These relevant data sets also reside in the PCDS. Problems associated with remote access to the PCDS were discussed in regard to all the PCDS data sets.

Heilman, W.↗

Simulation of collisional transport processes and the stability of planetary rings

The utility of the phase-space fluid method for the study of planetary ring dynamics is presently demonstrated through the numerical solution of a model kinetic equation for a flattened Keplerian disk. Attention is given to ringlets composed of single-sized particles, as well as to ringlets composed of two different-sized particles; in the latter case, the ringlets evolve in such a way that the lighter particles are confined by the heavier ones. The results obtained indicate that some natural process may sharpen the optical depth profile of edges even without an external forcing mechanism, and that intermediate optical depths are dynamically preferred in some cases.

Brophy, Thomas G.↗

Integrated flow and structural modeling for rocket engine component test facility propellant systems

A set of PC-based computational Dynamic Fluid Flow Simulation models is presented for modeling facility gas and cryogenic systems. Data obtained provide important information regarding performance envelope parameters for the facility using different engine components; time-dependent valve setting for controlling steady-state, quasi-steady state, and transient profiles; optimum facility pipe and pipe component sizes and parameters; momentum transfer loads; and fluid conditions at critical points. A set of COSMIC NASTRAN-based finite element models is also presented to evaluate the loads and stresses on test facility piping systems from fluid and gaseous effects, thermal chill down, and occasional wind loads. The models are based on Apple Macintosh software which makes it possible to change numerous parameters.

Dequay, L.↗

On turbulent spots in plane Poiseuille flow

Turbulence characteristics inside a turbulent spot in plume Poiseuille flow are investigated by analyzing a data base obtained from a direct numerical simulation. The spot is found to consist of two distinct regions - a turbulent area and a wave area. The flow inside the turbulent area has a strong resemblance to that found in the fully developed turbulent channel. Suitably defined mean and rms fluctuations as well as the internal shear layer structures are found to be similar to the turbulent counterpart. In the wave area, the inflexional mean spanwise profiles cause a rapid growth of oblique waves, which break down to turbulence. The breakdown process of the oblique waves is reminiscent of the secondary instability observed during transition to turbulence in channel and boundary layer flows. Other detailed characteristics associated with the Poiseuille spot are presented and are compared with experimental results.

Henningson, Dan S.↗

Energy buildup in sheared force-free magnetic fields

Photospheric displacement of the footpoints of solar magnetic field lines results in shearing and twisting of the field, and consequently in the buildup of electric currents and magnetic free energy in the corona. The sudden release of this free energy may be the origin of eruptive events like coronal mass ejections, prominence eruptions, and flares. An important question is whether such an energy release may be accompanied by the opening of magnetic field lines that were previously closed, for such open field lines can provide a route for matter frozen into the field to escape the sun altogether. This paper presents the results of numerical calculations showing that opening of the magnetic field is permitted energetically, in that it is possible to build up more free energy in a sheared, closed, force-free magnetic field than is in a related magnetic configuration having both closed and open field lines. Whether or not the closed force-free field attains enough energy to become partially open depends on the form of the shear profile; the results presented compare the energy buildup for different shear profiles. Implications for solar activity are discussed briefly.

Wolfson, Richard↗

Metallographic cooling rates of L-group ordinary chondrites

Shock metamorphism appears to be a ubiquitous feature in L-group ordinary chondrites. Brecciation and heterogeneous melting obscure much of the early history of this meteorite group and have caused confusion as to whether L chondrites have undergone thermal metamorphism within onion-shell or rubble-pile parent bodies. Employing the most recent shock criteria, we have examined 55 Antarctic and 24 non-Antarctic L chondrites in order to identify those which have been least affected by post-accretional shock. Six low-shock samples (those with shock grade less than S4) of petrographic types L3-L5 were selected from both populations and metallographic cooling rates were obtained following the technique of Willis and Goldstein. All non-Antarctic L6 chondrites inspected were too heavily shocked to be included in this group. However, 4 shocked L6 chondrites were analyzed in order to determine what effects shock may impose on metallographic cooling rates. Metallographic cooling rates were derived by analyzing the cores of taenite grains and then measuring the distance to the nearest grain edge. Taenites were identified using backscatter imaging on a Cameca SX-50 electron microprobe. Using backscatter we were able to locate homogeneous, rust-free, nearly spherical grains. M-shaped profiles taken from grain traverses were also used to help locate the central portions of selected grains. All points which contained phosphorus above detection limits were discarded. Plots of cooling-rate data are summarized and data from the high-shock samples are presented. The lack of coherency of cooling rates for individual samples is indicative of heterogeneous cooling following shock. The data confirms the statement expressed by numerous workers that extreme care must be taken when selecting samples of L chondrites for cooling-rate studies. Data for the 6 non-Antarctic low-shock samples are also presented. The samples display a general trend in cooling rates. The lowest metamorphic grade yielded the slowest cooling rates and an increase in grade follows an increase in cooling rate. This is the opposite relationship to that predicted by the onion-shell model.

Bennett, Marvin E.↗

Hydrodynamic Analyses and Evaluation of New Fluid Film Bearing Concepts

Over the past several years, numerical and experimental investigations have been performed on a waved journal bearing. The research work was undertaken by Dr. Florin Dimofte, a Senior Research Associate in the Mechanical Engineering Department at the University of Toledo. Dr. Theo Keith, Distinguished University Professor in the Mechanical Engineering Department was the Technical Coordinator of the project. The wave journal bearing is a bearing with a slight but precise variation in its circular profile such that a waved profile is circumscribed on the inner bearing diameter. The profile has a wave amplitude that is equal to a fraction of the bearing clearance. Prior to this period of research on the wave bearing, computer codes were written and an experimental facility was established. During this period of research considerable effort was directed towards the study of the bearing's stability. The previously developed computer codes and the experimental facility were of critical importance in performing this stability research. A collection of papers and reports were written to describe the results of this work. The attached captures that effort and represents the research output during the grant period.

Keith, Theo G., Jr.↗

Terrestrial Planet Formation Constrained By Mars and the Structure of the Asteroid Belt

Reproducing the large Earth/Mars mass ratio requires a strong mass depletion in solids within the protoplanetary disc between 1 and 3 au. The Grand Tack model invokes a specific migration history of the giant planets to remove most of the mass initially beyond 1 au and to dynamically excite the asteroid belt. However, one could also invoke a steep density gradient created by inward drift and pile-up of small particles induced by gas drag, as has been proposed to explain the formation of close-in super-Earths. Here we show that the asteroid belt’s orbital excitation provides a crucial constraint against this scenario for the Solar system. We performed a series of simulations of terrestrial planet formation and asteroid belt evolution starting from discs of planetesimals and planetary embryos with various radial density gradients and including Jupiter and Saturn on nearly circular and coplanar orbits. Discs with shallow density gradients reproduce the dynamical excitation of the asteroid belt by gravitational self-stirring but form Mars analogues significantly more massive than the real planet. In contrast, a disc with a surface density gradient proportional to r(exp −5.5) reproduces the Earth/Mars mass ratio but leaves the asteroid belt in a dynamical state that is far colder than the real belt. We conclude that no disc profile can simultaneously explain the structure of the terrestrial planets and asteroid belt. The asteroid belt must have been depleted and dynamically excited by a different mechanism such as, for instance, in the Grand Tack scenario.

Numerical methods↗

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗

Numerical assessment of the impact of the guiding-centre approximation on fast ion simulations in NSTX

Guiding-centre (GC) and full-orbit (FO) simulations of the beam-injected fast ion distribution and the corresponding neutron emissivity have been carried out for magnetohydrodynamics-quiescent National Spherical Torus eXperiment (NSTX) plasmas, using a combination of ASCOT5 and DRESS, to assess the suitability of the GC approximation for fast ions in NSTX. It was found that GC and FO simulations predicted substantially different steady-state distributions in both position and velocity space and different neutron emissivity profiles, leading to a 15% reduction in the predicted global neutron rate for FO relative to GC. These changes accompany a higher magnetic moment in FO, and correspond to a change in particle orbits from co-passing to trapped and stagnation orbits. ASCOT5 was also benchmarked against TRANSP/NUBEAM with input loaded entirely from TRANSP/NUBEAM output files, with agreement found between the GC simulations when finite Larmor radius (FLR) corrections were omitted. ASCOT5 FO and TRANSP/NUBEAM with FLR produced fast ion distributions which differed in localised regions, but predicted global neutron rates which agree within 3%.

ASCOT↗

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗

Direct numerical simulation of single-species and binary-species boundary layers at high pressure

Direct numerical simulations of single-species and binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate flow physics in a binary-species boundary layer at high pressure, with special attention to mass diffusion. The working fluids are nitrogen in the single-species cases, and a mixture of nitrogen and methane in the binary-species cases. An investigation of mean profiles shows that velocity and temperature profiles have steep gradients near the wall, whereas the mass fraction profiles do not have a steep gradient. This result indicates that the similarity law between velocity, temperature and mass fraction does not hold true in the binary-species boundary layer. The comparison of turbulent fluctuations shows that the qualitative characteristics of velocity and temperature are similar to each other. In contrast, profiles of mass fraction fluctuations are largely different from those of velocity and temperature. This result indicates that mass diffusion is not similar to momentum and thermal diffusions in binary-species boundary layers.

Toki, Takahiko↗

Soliton Generation in Negative Thermal Expansion Materials

Strain solitons have been observed statically in several 2D materials and dynamically in substrate materials using ultrafast laser pulses. The latter case relies on lattice relaxation in response to ultrafast heating in a light-absorbing transducer material, a process which is sensitive to the thermal expansion coefficient. Here we consider an unusual case where the sign of the thermal expansion coefficient is negative, a scenario which is experimentally feasible in light of rapid and recent advances in the discovery of negative thermal expansion materials. We present numerical solutions to a nonlinear differential equation which has been repeatedly demonstrated to quantitatively model experimental data and discuss the salient results using realistic parameters for material linear and nonlinear elasticity. The solitons that emerge from the initial value problem with negative and positive thermal expansion are qualitatively different in several ways. The new case of negative thermal expansion gives rise to a nearly-periodic soliton train with chirped profile and free of an isolated shock front. We suggest this unanticipated result may be realized experimentally and assess the potential for certain applications of this generic effect.

36 MATERIALS SCIENCE↗

Computer program for calculating laminar, transitional, and turbulent boundary layers for a compressible axisymmetric flow

A finite-difference program is described for calculating the viscous compressible boundary layer flow over either planar or axisymmetric surfaces. The flow may be initially laminar and progress through a transitional zone to fully turbulent flow, or it may remain laminar, depending on the imposed boundary conditions, laws of viscosity, and numerical solution of the momentum and energy equations. The flow may also be forced into a turbulent flow at a chosen spot by the data input. The input may contain the factors of arbitrary Reynolds number, free-stream Mach number, free-stream turbulence, wall heating or cooling, longitudinal wall curvature, wall suction or blowing, and wall roughness. The solution may start from an initial Falkner-Skan similarity profile, an approximate equilibrium turbulent profile, or an initial arbitrary input profile.

Albers, J. A.↗

Simulator study of a pictorial display for general aviation instrument flight

A simulation study of a computer drawn pictorial display involved a flight task that included an en route segment, terminal area maneuvering, a final approach, a missed approach, and a hold. The pictorial display consists of the drawing of boxes which either move along the desired path or are fixed at designated way points. Two boxes may be shown at all times, one related to the active way point and the other related to the standby way point. Ground tracks and vertical profiles of the flights, time histories of the final approach, and comments were obtained from time pilots. The results demonstrate the accuracy and consistency with which the segments of the flight are executed. The pilots found that the display is easy to learn and to use; that it provides good situation awareness, and that it could improve the safety of flight. The small size of the display, the lack of numerical information on pitch, roll, and heading angles, and the lack of definition of the boundaries of the conventional glide slope and localizer areas were criticized.

Adams, J. J.↗