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,117 records · Page 62

A unified viscous theory of lift and drag of 2-D thin airfoils and 3-D thin wings

A unified viscous theory of 2-D thin airfoils and 3-D thin wings is developed with numerical examples. The viscous theory of the load distribution is unique and tends to the classical inviscid result with Kutta condition in the high Reynolds number limit. A new theory of 2-D section induced drag is introduced with specific applications to three cases of interest: (1) constant angle of attack; (2) parabolic camber; and (3) a flapped airfoil. The first case is also extended to a profiled leading edge foil. The well-known drag due to absence of leading edge suction is derived from the viscous theory. It is independent of Reynolds number for zero thickness and varies inversely with the square root of the Reynolds number based on the leading edge radius for profiled sections. The role of turbulence in the section induced drag problem is discussed. A theory of minimum section induced drag is derived and applied. For low Reynolds number the minimum drag load tends to the constant angle of attack solution and for high Reynolds number to an approximation of the parabolic camber solution. The parabolic camber section induced drag is about 4 percent greater than the ideal minimum at high Reynolds number. Two new concepts, the viscous induced drag angle and the viscous induced separation potential are introduced. The separation potential is calculated for three 2-D cases and for a 3-D rectangular wing. The potential is calculated with input from a standard doublet lattice wing code without recourse to any boundary layer calculations. Separation is indicated in regions where it is observed experimentally. The classical induced drag is recovered in the 3-D high Reynolds number limit with an additional contribution that is Reynold number dependent. The 3-D viscous theory of minimum induced drag yields an equation for the optimal spanwise and chordwise load distribution. The design of optimal wing tip planforms and camber distributions is possible with the viscous 3-D wing theory.

Yates, John E.↗

A numerical model including PID control of a multizone crystal growth furnace

This paper presents a 2D axisymmetric combined conduction and radiation model of a multizone crystal growth furnace. The model is based on a programmable multizone furnace (PMZF) designed and built at NASA Lewis Research Center for growing high quality semiconductor crystals. A novel feature of this model is a control algorithm which automatically adjusts the power in any number of independently controlled heaters to establish the desired crystal temperatures in the furnace model. The control algorithm eliminates the need for numerous trial and error runs previously required to obtain the same results. The finite element code, FIDAP, used to develop the furnace model, was modified to directly incorporate the control algorithm. This algorithm, which presently uses PID control, and the associated heat transfer model are briefly discussed. Together, they have been used to predict the heater power distributions for a variety of furnace configurations and desired temperature profiles. Examples are included to demonstrate the effectiveness of the PID controlled model in establishing isothermal, Bridgman, and other complicated temperature profies in the sample. Finally, an example is given to show how the algorithm can be used to change the desired profile with time according to a prescribed temperature-time evolution.

Panzarella, Charles H.↗

Electro-Thermal Simulation Studies of SiC Junction Diodes Containing Screw Dislocations Under High Reverse-Bias Operation

The objective of this work was to conduct a modeling study of SiC P-N junction diodes operating under high reverse biased conditions. Analytical models and numerical simulation capabilities were to be developed for self-consistent electro-thermal analysis of the diode current-voltage (I-V) characteristics. Data from GRC indicate that screw dislocations are unavoidable in large area SiC devices, and lead to changes in the SiC diode electrical response characteristics under high field conditions. For example, device instability and failures linked to internal current filamentation have been observed. The physical origin of these processes is not well understood, and quantitative projections of the electrical behavior under high field and temperature conditions are lacking. Thermal calculations for SiC devices have not been reported in the literature either. So estimates or projections of peak device temperatures and power limitations do not exist. This numerical study and simulation analysis was aimed at resolving some of the above issues. The following tasks were successfully accomplished: (1) Development of physically based models using one- and two-dimensional drift-diffusion theory for the transport behavior and I-V characteristics; (2) One- and two-dimensional heat flow to account for internal device heating. This led to calculations of the internal temperature profiles, which in turn, were used to update the electrical transport parameters for a self-consistent analysis. The temperature profiles and the peak values were thus obtainable for a given device operating condition; (3) Inclusion of traps assumed to model the presence of internal screw dislocations running along the longitudinal direction; (4) Predictions of the operating characteristics with and without heating as a function of applied bias with and without traps. Both one and two-dimensional cases were implemented; (5) Assessment of device stability based on the operating characteristics. The presence of internal non-uniformities, particularly filamentary structures, was probed and demonstrated; (6) Cause and physical origins of filamentary behavior and unstable I-V characteristics were made transparent; (7) It was demonstrated that diodes containing defects would be more prone to thermal breakdown associated with the temperature dependent decrease in the thermal conductivity; and (8) Finally, negative differential resistance (S-shaped NDR) which can potential lead to device instability and filamentary behavior was shown to occur for diodes containing a line of defects such as could be associated with a screw dislocation line.

Joshi, R. P.↗

Modified Involute Helical Gears: Computerized Design, Simulation of Meshing, and Stress Analysis

The computerized design, methods for generation, simulation of meshing, and enhanced stress analysis of modified involute helical gears is presented. The approaches proposed for modification of conventional involute helical gears are based on conjugation of double-crowned pinion with a conventional helical involute gear. Double-crowning of the pinion means deviation of cross-profile from an involute one and deviation in longitudinal direction from a helicoid surface. Using the method developed, the pinion-gear tooth surfaces are in point-contact, the bearing contact is localized and oriented longitudinally, and edge contact is avoided. Also, the influence of errors of aligment on the shift of bearing contact, vibration, and noise are reduced substantially. The theory developed is illustrated with numerical examples that confirm the advantages of the gear drives of the modified geometry in comparison with conventional helical involute gears.

Handschuh, Robert↗

Modified Involute Helical Gears: Computerized Design, Simulation of Meshing and Stress Analysis

The computerized design, methods for generation, simulation of meshing, and enhanced stress analysis of modified involute helical gears is presented. The approaches proposed for modification of conventional involute helical gears are based on conjugation of double-crowned pinion with a conventional helical involute gear. Double-crowning of the pinion means deviation of cross-profile from an involute one and deviation in longitudinal direction from a helicoid surface. Using the method developed, the pinion-gear tooth surfaces are in point-contact, the bearing contact is localized and oriented longitudinally, and edge contact is avoided. Also, the influence of errors of alignment on the shift of bearing contact, vibration, and noise are reduced substantially. The theory developed is illustrated with numerical examples that confirm the advantages of the gear drives of the modified geometry in comparison with conventional helical involute gears.

COMPUTER AIDED DESIGN↗

Velocity Measurements in the Wake of the Swept Wing Flow Test Model at the National Transonic Facility

Femtosecond laser electronic excitation tagging (FLEET) was applied to obtain flowfield velocity data for the Swept Wing Flow Test (SWiFT) at the NASA Langley Research Center National Transonic Facility (NTF). Despite numerous challenges associated with performing flow velocimetry measurements within a large-scale cryogenic wind tunnel facility, the experimental campaign was conducted under various conditions, including dry air at 320 K, cool nitrogen at 240 K, and Mach numbers of 0.2 and 0.8. FLEET velocimetry measurements were performed in the downstream wake of the SWiFT model, providing a quantitative dataset. This dataset includes one-dimensional velocity profiles and two-dimensional velocity maps acquired at different angles of attack and Reynolds numbers. The two-dimensional single component velocity maps indicate relatively uniform flow across the 120-mm wide wake flow survey suggesting that the main flow features can be represented by single-position linear velocity profiles. The measured instantaneous velocity profiles at variable angles of attack are compared to model lift coefficient information obtained at Mach 0.2. At high Reynolds numbers, a sudden velocity decrease was observed in the FLEET measurements at the same time as stall in lift coefficient was observed. At a low Reynolds number, both the velocity profiles and the lift coefficient show a smoother transition, without a sudden stall. Furthermore, the two-dimensional, one component velocity map reveals a velocity deficit region at Mach 0.8 at various angles of attack. Both single shot and mean velocity measurements were acquired allowing assessment of flowfield fluctuations and measurement precisions. The uncertainties are within 4 m/s in mean measurements based on repeatability data and about 5 m/s in instantaneous single-shot measurements. Measurements are reported with ~4.8 mm spatial resolution with 32-pixel averaging used to reduce measurement errors.

femtosecond↗

Velocity Measurements in the Wake of the Swept Wing Flow Test (SWIFT) Model at the National Transonic Facility

Femtosecond laser electronic excitation tagging (FLEET) was applied to obtain flowfield velocity data for the Swept Wing Flow Test (SWiFT) at the NASA Langley Research Center National Transonic Facility (NTF). Despite numerous challenges associated with performing flow velocimetry measurements within a large-scale cryogenic wind tunnel facility, the experimental campaign was conducted under various conditions, including dry air at 320 K, cool nitrogen at 240 K, and Mach numbers of 0.2 and 0.8. FLEET velocimetry measurements were performed in the downstream wake of the SWiFT model, providing a quantitative dataset. This dataset includes one-dimensional velocity profiles and two-dimensional velocity maps acquired at different angles of attack and Reynolds numbers. The two-dimensional single component velocity maps indicate relatively uniform flow across the 120-mm wide wake flow survey suggesting that the main flow features can be represented by single-position linear velocity profiles. The measured instantaneous velocity profiles at variable angles of attack are compared to model lift coefficient information obtained at Mach 0.2. At high Reynolds numbers, a sudden velocity decrease was observed in the FLEET measurements at the same time as stall in lift coefficient was observed. At a low Reynolds number, both the velocity profiles and the lift coefficient show a smoother transition, without a sudden stall. Furthermore, the two-dimensional, one component velocity map reveals a velocity deficit region at Mach 0.8 at various angles of attack. Both single shot and mean velocity measurements were acquired allowing assessment of flowfield fluctuations and measurement precisions. The uncertainties are within 4 m/s in mean measurements based on repeatability data and about 5 m/s in instantaneous single-shot measurements. Measurements are reported with ~4.8 mm spatial resolution with 32-pixel averaging used to reduce measurement errors.

Transonic↗

Electrostatic analysis of charge-coupled structures

An analysis is presented which is based upon a numerical solution of Gauss's law for the multiple dielectric semiconductor and electrode structure. The formulation resulted in a large set of equations, usually nonlinear, which were solved by an iterative scheme based upon the Gauss-Seidel method employing a relaxation parameter. Proper consideration of the nonlinearity of the space-charge in the semiconductor allowed convergent solutions which can give reasonable approximations to the channel potential profiles though not as good for estimating the space charge itself. The results of this work, which are three computer programs, are listed. None of the programs give the electric field explicitly as output data. However, when this information is desired, the programs can be expanded to include field estimations based on a polynominal fit to the potential field, or they may serve as subprograms for main programs dealing with general aspects of charge transfer.

Gassaway, J. D.↗

A model of active prominences

A one-dimensional numerical model of active loop prominences is investigated. The model includes the effects of gravity, the geometry of the magnetic field, and conduction losses to the chromosphere. Calculations indicate that, as originally proposed by Goldsmith, the thermal instability mechanism is, by itself, sufficient to account for the appearance of bright H-alpha knots in postflare loops. Under certain conditions, initial perturbations in the loop temperature and density profiles of small, but finite, amplitude (approximately 5%) and large size scale (greater than or approximately equal to 10 to the 9th cm) can grow into condensations with temperature and density differences of over an order of magnitude and size scales of less than 10 to the 8th cm. In agreement with observations, the conditions that must be satisfied are such that loop prominence systems are likely to occur only in large flares and such that knots preferentially form at the tops of loops. The velocities, densities, and lifetimes calculated for the loop material are also in agreement with observations. It is concluded that in order for H-alpha knots to occur, heating of some form must continue into the late cooling phase of a flare loop, and that this heating is more intense near the loop base than near the apex.

Antiochos, S. K.↗

Numerical calculations relevant to the initial expansion of the polar wind

The results of a numerical simulation of the expansion of an H(+)-0(+) electron plasma into a vacuum are reported. Variations were introduced in the ionic density ratios, the scale lengths for the density gradient at the plasma-vacuum interface, and the initial electron-ion temperature ratio. Conditions of a gravitationally bound O(+) plasma were also considered. It was found that if H(+) is a major ion, the H(+) density profile in the expansion region is always concave, with a phase-space divided into three regions, i.e., an undisturbed plasma, a rarefaction region, and an expansion region in which ions are accelerated. When the H(+) is a minor ion, the O(+) ions accelerate the H(+) ions, the H(+) phase-space is divided into five regions, and the increased electron temperatures produce enhanced electrostatic potentials and H(+) drifts in the expansion region.

Singh, N.↗

Numerical simulation of the atmosphere during a CAT encounter

In an attempt to determine whether clear air turbulence (CAT) which caused a brief, intense turbulence encounter to a jetliner might have been related to gravity waves generated by a line of thunderstorms, local sounding data taken within two hours of the incident have been used as input to a numerical model designed to simulate stratified flow over obstacles. In the model a rigid obstacle at the lower boundary acted as the source of the gravity waves. The model results show a large amplitude disturbance localized over the obstacle at and above flight level, in regions where CAT would be expected. Stratification studies have been conducted, varying the height and shape of the obstacle and the input wind profile. The success of the experiment leads to the conclusion that atmospheric gravity wave simulation can be useful in understanding and possibly avoiding CAT.

Keller, T. L.↗

Sloshing dynamics on rotating helium dewar tank

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by both the gravity gradient and jitter accelerations applicable to scientific spacecraft which is eligible to carry out spinning motion and/or slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with Gravity Probe-B (GP-B) spacecraft which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics were based on the non-inertia frame spacecraft bound coordinate, and solve time dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers were derived. Results were widely published in the open journals.

Hung, R. J.↗

Onboard Atmospheric Modeling and Prediction for Autonomous Aerobraking Missions

Aerobraking has proven to be an effective means of increasing the science payload for planetary orbiting missions and/or for enabling the use of less expensive launch vehicles. Though aerobraking has numerous benefits, large operations cost have been required to maintain the aerobraking time line without violating aerodynamic heating or other constraints. Two operations functions have been performed on an orbit by orbit basis to estimate atmospheric properties relevant to aerobraking. The Navigation team typically solves for an atmospheric density scale factor using DSN tracking data and the atmospheric modeling team uses telemetric accelerometer data to recover atmospheric density profiles. After some effort, decisions are made about the need for orbit trim maneuvers to adjust periapsis altitude to stay within the aerobraking corridor. Autonomous aerobraking would reduce the need for many ground based tasks. To be successful, atmospheric modeling must be performed on the vehicle in near real time. This paper discusses the issues associated with estimating the planetary atmosphere onboard and evaluates a number of the options for Mars, Venus and Titan aerobraking missions.

Tolson, Robert H.↗

Evaluating the Contribution of NASA Remotely-Sensed Data Sets on a Convection-Allowing Forecast Model

The Short-term Prediction Research and Transition (SPoRT) Center is a collaborative partnership between NASA and operational forecasting partners, including a number of National Weather Service forecast offices. SPoRT provides real-time NASA products and capabilities to help its partners address specific operational forecast challenges. One challenge that forecasters face is using guidance from local and regional deterministic numerical models configured at convection-allowing resolution to help assess a variety of mesoscale/convective-scale phenomena such as sea-breezes, local wind circulations, and mesoscale convective weather potential on a given day. While guidance from convection-allowing models has proven valuable in many circumstances, the potential exists for model improvements by incorporating more representative land-water surface datasets, and by assimilating retrieved temperature and moisture profiles from hyper-spectral sounders. In order to help increase the accuracy of deterministic convection-allowing models, SPoRT produces real-time, 4-km CONUS forecasts using a configuration of the Weather Research and Forecasting (WRF) model (hereafter SPoRT-WRF) that includes unique NASA products and capabilities including 4-km resolution soil initialization data from the Land Information System (LIS), 2-km resolution SPoRT SST composites over oceans and large water bodies, high-resolution real-time Green Vegetation Fraction (GVF) composites derived from the Moderate-resolution Imaging Spectroradiometer (MODIS) instrument, and retrieved temperature and moisture profiles from the Atmospheric Infrared Sounder (AIRS) and Infrared Atmospheric Sounding Interferometer (IASI). NCAR's Model Evaluation Tools (MET) verification package is used to generate statistics of model performance compared to in situ observations and rainfall analyses for three months during the summer of 2012 (June-August). Detailed analyses of specific severe weather outbreaks during the summer will be presented to assess the potential added-value of the SPoRT datasets and data assimilation methodology compared to a WRF configuration without the unique datasets and data assimilation.

Zavodsky, Bradley T.↗

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.↗