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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Computation of turbulent boundary layers employing the defect wall-function method

In order to decrease overall computational time requirements of spatially-marching parabolized Navier-Stokes finite-difference computer code when applied to turbulent fluid flow, a wall-function methodology, originally proposed by R. Barnwell, was implemented. This numerical effort increases computational speed and calculates reasonably accurate wall shear stress spatial distributions and boundary-layer profiles. Since the wall shear stress is analytically determined from the wall-function model, the computational grid near the wall is not required to spatially resolve the laminar-viscous sublayer. Consequently, a substantially increased computational integration step size is achieved resulting in a considerable decrease in net computational time. This wall-function technique is demonstrated for adiabatic flat plate test cases from Mach 2 to Mach 8. These test cases are analytically verified employing: (1) Eckert reference method solutions, (2) experimental turbulent boundary-layer data of Mabey, and (3) finite-difference computational code solutions with fully resolved laminar-viscous sublayers. Additionally, results have been obtained for two pressure-gradient cases: (1) an adiabatic expansion corner and (2) an adiabatic compression corner.

Brown, Douglas L.↗
NASA NTRS1999

On the Role of Surface Friction in Tropical Intraseasonal Oscillation

The Madden-Julian oscillation (MJO), or the tropical intraseasonal oscillation, has attracted much attention, ever since its discovery in the early seventies for reasons of both scientific understanding and practical forecasts. Among the theoretical interpretations of the MJO, the wave-CISK (conditional instability of the second kind) mechanism is the most popular. The basic idea of the wave-CISK interpretation is that the cooperation between the low-level convergence associated with the eastward moving Kelvin wave and the cumulus convection generates an eastward moving Kelvin-wave-like mode. Later it was recognized that the MJO has an important Rossby-wave-like component. However linear analysis and numerical simulations based on it (even when conditional heating is used) have revealed two problems with the wave-CISK interpretation; i.e., excessive speed and the most preferred scale being zero or grid scale. Chao (1995) presented a discussion of these problems and attributed these problems to the particular type of expression for the cumulus heating used in the linear analyses and numerical studies (i.e., the convective heating is proportional to low-level convergence and a fixed vertical heating profile). It should be pointed out that in the relatively successful simulation of MJO with general circulation models the problem of grid scale being the most preferred scale does not appear and die problem of excessive speed is not as severe as in the linear analysis.

Chao, Winston C.↗
NASA NTRS1998

Guidance and Control Concepts for the X-33 Technology Demonstrator

The X-33 technology demonstrator is a suborbital precursor to the Reusable Launch Vehicle (RLV) with first flight planned for summer of 1999. The flight test program will include about 15 flights originating from Edwards Air Force Base, California, each with widely varying flight profiles in order to test new thermal protection system (TPS) materials, structures, and linear aerospike engines. The first flights will be relatively short range flights with about a 300 nmi range, maximum Mach number of 7, maximum altitude of 190,000 feet, whereas the latter flights will cover about 800 nmi range, with max altitude of about 260,000 feet and max Mach of about 15. The guidance algorithms must be flexible enough to accommodate these various profiles and to adapt to severe off-nominal dispersions, such as early engine failure (partial or total) where possibly more than half the thrust is lost. An onboard real-time performance monitor will be used to assess the viability of the nominal landing site as well as alternate landing sites that would potentially be used in extreme off-nominal conditions. During ascent, a single entry guidance-related parameter, which is easy to calculate, is used to assess the viability of the nominal landing site as well as alternate landing sites. Real-time adjustment of the stored ascent attitude profile will be performed, as required, to maximize the probability of making it to the nominal landing site. Numerical results are given for various engine-out cases to illustrate the adaptability of the performance monitor.

Dukeman, Gregory A.↗
NASA NTRS2003

Magnetohydrodynamic Jump Conditions for Oblique Relativistic Shocks with Gyrotropic Pressure

Shock jump conditions, i.e., the specification of the downstream parameters of the gas in terms of the upstream parameters, are obtained for steady-state, plane shocks with oblique magnetic fields and arbitrary flow speeds. This is done by combining the continuity of particle number flux and the electromagnetic boundary conditions at the shock with the magnetohydrodynamic conservation laws derived from the stress-energy tensor. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be solved analytically. For mildly relativistic shocks, analytic solutions are obtained for isotropic pressure using an approximation for the adiabatic index that is valid in high sonic Mach number cases. Examples assuming isotropic pressure illustrate how the shock compression ratio depends on the shock speed and obliquity. In the more general case of gyrotropic pressure, the jump conditions cannot be solved analytically with- out additional assumptions, and the effects of gyrotropic pressure are investigated by parameterizing the distribution of pressure parallel and perpendicular to the magnetic field. Our numerical solutions reveal that relatively small departures from isotropy (e.g., approximately 20%) produce significant changes in the shock compression ratio, r , at all shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with gyrotropic pressure is obtained. In particular, either dynamically important fields or significant pressure anisotropies can incur marked departures from the canonical gas dynamic value of r = 3 for a shocked ultrarelativistic flow and this may impact models of particle acceleration in gamma-ray bursts and other environments where relativistic shocks are inferred. The jump conditions presented apply directly to test-particle acceleration, and will facilitate future self-consistent numerical modeling of particle acceleration at oblique, relativistic shocks; such models include the modification of the fluid velocity profile due to the contribution of energetic particles to the momentum and energy fluxes.

Double, Glen P.↗
NASA NTRS2006

Validation of High-Speed Turbulent Boundary Layer and Shock-Boundary Layer Interaction Computations with the OVERFLOW Code

The capability of the OVERFLOW code to accurately compute high-speed turbulent boundary layers and turbulent shock-boundary layer interactions is being evaluated. Configurations being investigated include a Mach 2.87 flat plate to compare experimental velocity profiles and boundary layer growth, a Mach 6 flat plate to compare experimental surface heat transfer,a direct numerical simulation (DNS) at Mach 2.25 for turbulent quantities, and several Mach 3 compression ramps to compare computations of shock-boundary layer interactions to experimental laser doppler velocimetry (LDV) data and hot-wire data. The present paper describes outlines the study and presents preliminary results for two of the flat plate cases and two small-angle compression corner test cases.

Oliver, A. B.↗
NASA NTRS2016

A Study of Heat Transfer and Flow Characteristics of Rising Taylor Bubbles

Practical application of flow boiling to ground- and space-based thermal management systems hinges on the ability to predict the systems heat removal capabilities under expected operating conditions. Research in this field has shown that the heat transfer coefficient within two-phase heat exchangers can be largely dependent on the experienced flow regime. This finding has inspired an effort to develop mechanistic heat transfer models for each flow pattern which are likely to outperform traditional empirical correlations. As a contribution to the effort, this work aimed to identify the heat transfer mechanisms for the slug flow regime through analysis of individual Taylor bubbles.An experimental apparatus was developed to inject single vapor Taylor bubbles into co-currently flowing liquid HFE 7100. The heat transfer was measured as the bubble rose through a 6 mm inner diameter heated tube using an infrared thermography technique. High-speed flow visualization was obtained and the bubble film thickness measured in an adiabatic section. Experiments were conducted at various liquid mass fluxes (43-200 kgm2s) and gravity levels (0.01g-1.8g) to characterize the effect of bubble drift velocityon the heat transfer mechanisms. Variable gravity testing was conducted during a NASA parabolic flight campaign.Results from the experiments showed that the drift velocity strongly affects the hydrodynamics and heat transfer of single elongated bubbles. At low gravity levels, bubbles exhibited shapes characteristic of capillary flows and the heat transfer enhancement due to the bubble was dominated by conduction through the thin film. At moderate to high gravity, traditional Taylor bubbles provided small values of enhancement within the film, but large peaks in the wake heat transfer occurred due to turbulent vortices induced by the film plunging into the trailing liquid slug. Characteristics of the wake heat transfer profiles were analyzed and related to the predicted velocity field. Results were compared and shown to agree with numerical simulations of colleagues from EPFL, Switzerland.In addition, a preliminary study was completed on the effect of a Taylor bubble passing through nucleate flow boiling, showing that the thinning thermal boundary layer within the film suppressed nucleation, thereby decreasing the heat transfer coefficient.

capillary flow↗
NASA NTRS1993

Particle-Gas Dynamics in the Midplane of a Protoplanetary Nebula

In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.

Cuzzi, Jeffrey N.↗
NASA NTRS1994

Numerical and Experimental Study of Transport Phenomena in Directional Solidification of Succinonitrile

A numerical and experimental study of the growth of succinonitrile (SCN) using a horizontal Bridginan furnace and transparent glass ampoule was conducted. Two experiments were considered: one in which the temperature profile was fixed relative to the ampoule (no-growth case); and a second in which the thermal profile was translated at a constant rate (steady growth case). Measured temperature profiles on the outer surface of the ampoule were used as thermal boundary conditions for the modelling. The apparent heat capacity formulation combined with the variable viscositymeth~ was used to model the phase change in SeN. Both 2-D and 3-D models were studied and numerical solutions obtained using the commercial finite element code, FIDAP1. Comparison of the numerical results to experimental data showed excellent agreement. The complex 3-D shallow-cavity flow in the melt, differences between 2-D and 3-D models, effects of natural convection on the thermal gradient and shape of the solid/liquid interface, and the sensitivity of simulations to specific assumptions, are also discussed.

de Groh, Henry C., III↗
NASA NTRS1989

Comparison of LDV measurements and Navier-Stokes solutions in a two-dimensional 180-degree turn-around duct

Results from an experimental and numerical investigation of turbulent subsonic flow inside a two-dimensional, strongly curved, 180 deg turn-around duct are presented. Data measured with a two-component, two-color laser Doppler velocimeter include profiles of mean axial velocity and local flow angle. Static pressure distributions are also measured. Results are obtained at a Mach number of 0.1, and at Reynolds numbers of 100,000 and a million based on channel height. Numerical calculations are performed using the incompressible Navier-Stokes equations with a Prandtl mixing length zero-equation turbulence model modified for internal flows. Theory and experiment are compared to evaluate the ability of the turbulence model to predict this class of internal flows with strong curvature.

Monson, Daryl J.↗
NASA NTRS2015

A Numerical and Experimental Study of Coflow Laminar Diffusion Flames: Effects of Gravity and Inlet Velocity

In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both computationally and experimentally. A series of flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity conditions with the fuel stream CH4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth vorticity-velocity formulation of the governing equations was employed to describe the reactive gaseous mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional computational domain was employed, where the grid on the axisymmetric domain was a nonuniform tensor product mesh. The governing equations and boundary conditions were discretized on the mesh by a nine-point finite difference stencil, with the convective terms approximated by a monotonic upwind scheme and all other derivatives approximated by centered differences. The resulting set of fully coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton's method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off height, and soot temperature were determined by flame emission images recorded by a digital camera, and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first compared to assess the accuracy of the numerical model. After its validity was established, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of laminar coflow methane-air diffusion flames was explored further by examining quantities derived from the computational results.

microgravity↗
NASA NTRS2005

Analysis of plasmaspheric plumes: CLUSTER and IMAGE observations and numerical simulations

Plasmaspheric plumes have been routinely observed by CLUSTER and IMAGE. The CLUSTER mission provides high time resolution four-point measurements of the plasmasphere near perigee. Total electron density profiles can be derived from the plasma frequency and/or from the spacecraft potential (note that the electron spectrometer is usually not operating inside the plasmasphere); ion velocity is also measured onboard these satellites (but ion density is not reliable because of instrumental limitations). The EUV imager onboard the IMAGE spacecraft provides global images of the plasmasphere with a spatial resolution of 0.1 RE every 10 minutes; such images acquired near apogee from high above the pole show the geometry of plasmaspheric plumes, their evolution and motion. We present coordinated observations for 3 plume events and compare CLUSTER in-situ data (panel A) with global images of the plasmasphere obtained from IMAGE (panel B), and with numerical simulations for the formation of plumes based on a model that includes the interchange instability mechanism (panel C). In particular, we study the geometry and the orientation of plasmaspheric plumes by using a four-point analysis method, the spatial gradient. We also compare several aspects of their motion as determined by different methods: (i) inner and outer plume boundary velocity calculated from time delays of this boundary observed by the wave experiment WHISPER on the four spacecraft, (ii) ion velocity derived from the ion spectrometer CIS onboard CLUSTER, (iii) drift velocity measured by the electron drift instrument ED1 onboard CLUSTER and (iv) global velocity determined from successive EUV images. These different techniques consistently indicate that plasmaspheric plumes rotate around the Earth, with their foot fully co-rotating, but with their tip rotating slower and moving farther out.

Darouzet, Fabien↗
NASA NTRS2000

Transverse Mode Dynamics of VCSELs Undergoing Current Modulation

Transverse mode dynamics of a 20-micron-diameter vertical-cavity surface-emitting laser (VCSEL) undergoing gain switching by deep current modulation is studied numerically. The direct current (dc) level is set slightly below threshold and is modulated by a large alternating current (ac). The resulting optical pulse train and transverse-mode patterns are obtained numerically. The ac frequency is varied from 2.5 GHz to 10 GHz, and the ac amplitude is varied from one-half to four times that of the dc level. At high modulation frequencies, a regular pulse train is not generated unless the ac amplitude is large enough. At all modulation frequencies, the transverse spatial profile switches from single-mode to multiple-mode pattern as the ac pumping level is increased. Optical pulse widths vary in the range 5-30 ps. with the pulse width decreasing when either the frequency is increased or the ac amplitude is decreased. The numerical modeling uses an approximation form of the semiconductor Maxwell-Bloch equations. Temporal evolution of the spatial profiles of the laser (and of carrier density) is determined without any assumptions about the type or number of modes. Keywords: VCSELs, current modulation, gain switching, transverse mode dynamics, computational modeling

Goorjian, Peter M.↗
NASA NTRS1973

Global electromagnetic induction in the moon and planets

Experiments and analyses concerning electromagnetic induction in the moon and other extraterrestrial bodies are summarized. The theory of classical electromagnetic induction in a sphere is first considered, and this treatment is extended to the case of the moon, where poloidal eddy-current response has been found experimentally to dominate other induction modes. Analysis of lunar poloidal induction yields lunar internal electrical conductivity and temperature profiles. Two poloidal-induction analytical techniques are discussed: a transient-response method applied to time-series magnetometer data, and a harmonic-analysis method applied to data numerically Fourier-transformed to the frequency domain, with emphasis on the former technique. Attention is given to complicating effects of the solar wind interaction with both induced poloidal fields and remanent steady fields. The static magnetization field induction mode is described, from which are calculated bulk magnetic permeability profiles. Magnetic field measurements obtained from the moon and from fly-bys of Venus and Mars are studied to determine the feasibility of extending theoretical and experimental induction techniques to other bodies in the solar system.

Dyal, P.↗
NASA NTRS1990

Fluid dynamics and low gravity effects of chemical vapor deposition

Based on the comparison between experimental data and numerical results for the growth of GaAs from TMGa, it was shown that 3D simulations are necessary to simulate rectangular CVD reactors even when operated under subcritical (Ra) conditions. The important points found are summarized in the three attached reprints. The experimental studies of mixed convection in horizontal channels have shown three regimes of high Ra (22,220) number flows. At Re = 18.5, the rolls develop very quickly, significantly modulating the axial velocity even before it reaches the beginning of the hot plate. A few centimeters downstream, the velocities become asymmetric about the vertical centerplane and at x = 12 cm, become unsteady. These asymmetries were predicted theoretically, but experimental evidence has not been published prior to this work. At Re = 36, the axial velocity is only slightly modified at x = 0. Although the flow remains steady and symmetric about the vertical centerplane, there is a small spatial oscillation in the velocities over the length of the channel. The period of this oscillation was around 5 cm. At Re = 54, the longitudinal rolls developed smoothly over a length of 30 cm, with no asymmetries, unsteadiness, or spatial oscillations. Comparison of numerical simulations of these flows to experiments has revealed the importance and difficulty of setting proper thermal boundary conditions on the sidewalls. Calculated flows and experimentally measured flows showed very similar profiles, but at different axial locations, with the rolls developing more rapidly in the experiments. This is directly attributable to partially conducting sidewalls of the apparatus being hotter in the entrance section than the adiabatic walls of the simulations. A thorough comparison of the experimental data and numerical results for a variety of sidewall boundary conditions is in preparation.

Nyce, Thomas A.↗
NASA NTRS1984

Boundary-layer linear stability theory

Most fluid flows are turbulent rather than laminar and the reason for this was studied. One of the earliest explanations was that laminar flow is unstable, and the linear instability theory was first developed to explore this possibility. A series of early papers by Rayleigh produced many notable results concerning the instability of inviscid flows, such as the discovery of inflectional instability. Viscosity was commonly thought to act only to stabilize the flow, and flows with convex velocity profiles appeared to be stable. The investigations that led to a viscous theory of boundary layer instability was reported. The earliest application of linear stability theory to transition prediction calculated the amplitude ratio of the most amplified frequency as a function of Reynolds number for a Blasius boundary layer, and found that this quantity had values between five and nine at the observed Ret. The experiment of Schubauer and Skramstad (1947) completely reversed the prevailing option and fully vindicated the Gottingen proponents of the theory. This experiment demonstrated the existence of instability waves in a boundary layer, their connection with transition, and the quantitative description of their behavior by the theory of Tollmien and Schlichting. It is generally accepted that flow parameters such as pressure gradient, suction and heat transfer qualitatively affect transition in the manner predicted by the linear theory, and in particular that a flow predicted to be stable by the theory should remain laminar. The linear theory, in the form of the e9, or N-factor is today in routine use in engineering studies of laminar flow. The stability theory to boundary layers with pressure gradients and suction was applied. The only large body of numerical results for exact boundary layer solutions before the advent of the computer age by calculating the stability characteristics of the Falkner-Skan family of velocity profiles are given. When the digital computer reached a stage of development which permit the direct solution of the primary differential equations, numerical results were obtained from the linear theory during the next 10 years for many different boundary layer flows: three dimensional boundary layers; free convention boundary layers; compressible boundary layers; boundary layers on compliant walls; a recomputation of Falkner-Skan flows; unsteady boundary layers; and heated wall boundary layers.

Leslie M Mack↗
NASA NTRS1992

A quantitative measure of the structure of gamma-ray burst time profiles

A cursory examination of cosmic gamma-ray burst time profiles indicates an inhomogeneous distribution of structure. In the first approximation, there seem to be two types of profiles; smooth ones with little structure and highly variable ones with lots of structure. To put this observation to the test, we have examined the statistical nature of the profile derivative to choose which parameter might best be called the burst 'spikiness'. We have found that a good estimator is given by a count of the number of 'spikes' (defined by a specific numerical recipe) and not by the rms deviations from either a pre-burst background or any type of moving average background. The application of this parameter to 30 burst time histories shows it to be consistent over a wide range of profile types. The analysis also reveals a preferred average time between spikes of approximately 1.5 seconds.

Lestrade, John P.↗
NASA NTRS1979

Application of numerical optimization to the design of advanced supercritical airfoils

An application of numerical optimization to the design of advanced airfoils for transonic aircraft showed that low-drag sections can be developed for a given design Mach number without an accompanying drag increase at lower Mach numbers. This is achieved by imposing a constraint on the drag coefficient at an off-design Mach number while minimizing the drag coefficient at the design Mach number. This multiple design-point numerical optimization has been implemented with the use of airfoil shape functions which permit a wide range of attainable profiles during the optimization process. Analytical data for the starting airfoil shape, a single design-point optimized shape, and a double design-point optimized shape are presented. Experimental data obtained in the NASA Ames two-by two-foot wind tunnel are also presented and discussed.

Johnson, R. R.↗