SEARCH · Engineering Papers

Results for “numerical computation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Characteristics of Non-Premixed Turbulent Flames in Microgravity

This project is concerned with the characteristics of turbulent hydrocarbon (primarily propane) gas-jet diffusion flames in microgravity. A microgravity environment provides the opportunity to study the structure of turbulent diffusion flames under momentum-dominated conditions (large Froude number) at moderate Reynolds number which is a combination not achievable in normal gravity. This paper summarizes progress made since the last workshop. Primarily, the features of flame radiation from microgravity turbulent jet diffusion flames in a reduced gravity environment are described. Tests were conducted for non-premixed, nitrogen diluted propane flames burning in quiescent air in the NASA Glenn 5.18 Second Zero Gravity Facility. Measured flame radiation from wedge-shaped, axial slices of the flame are compared for microgravity and normal gravity flames. Results from numerical computations of the flame using a k-e model for the turbulence are also presented to show the effects of flame radiation on the thermal field. Flame radiation is an important quantity that is impacted by buoyancy as has been shown in previous studies by the authors and also by Urban et al. It was found that jet diffusion flames burning under microgravity conditions have significantly higher radiative loss (about five to seven times higher) compared to their normal gravity counterparts because of larger flame size in microgravity and larger convective heat loss fraction from the flame in normal gravity. These studies, however, were confined to laminar flames. For the case of turbulent flames, the flame radiation is a function of time and both the time-averaged and time-dependent components are of interest. In this paper, attention is focused primarily on the time-averaged level of the radiation but the turbulent structure of the flame is also assessed from considerations of the radiation power spectra.

Hegde, U.↗

NASA NTRS2001

Outwardly Propagating Flames at Elevated Pressures

Spherical, outwardly-propagating flames of CH4-O2-inert and H2-O2-inert mixtures were experimentally studied in a high pressure apparatus. Stretch-free flame speeds and Markstein lengths were extracted for a wide range of pressures and equivalence ratios for spherically-symmetric, smooth flamefronts and compared to numerical computations with detailed chemistry and transport, as well as existing data in the literature. Wrinkle development was examined for propagating flames that were unstable under our experimental conditions. Hydrodynamic cells developed for most H2-air and CH4-air flames at elevated pressures, while thermal-diffusive instabilities were also observed for lean and near-stoichiometric hydrogen flames at pressures above atmospheric. Strategies in suppressing or delaying the onset of cell formation have been assessed. Buoyancy effects affected sufficiently off-stoichiometric CH4 mixtures at high pressures.

Law, C. K.↗

NASA NTRS2001

The Capture of Interstellar Dust: The Pure Poynting-Robertson Case

Ulysses and Galileo spacecraft have discovered interstellar dust particles entering the solar system. In general, particles trajectories not altered by Lorentz forces or radiation pressure should encounter the sun on open orbits. Under Newtonian forces alone these particles return to the interstellar medium. Dissipative forces, such as Poynting Robertson (PR) and corpuscular drag and non-dissipative Lorentz forces can modify open orbits to become closed. In particular, it is possible for the orbits of particles that pass close to the Sun to become closed due to PR drag. Further, solar irradiation will cause modification of the size of the dust particle by evaporation. The combination of these processes gives rise a class of capture orbits and bound orbits with evaporation. Considering only the case of pure PR drag a minimum impact parameter is derived for initial capture by Poynting-Robertson drag. Orbits in the solar radiation field are computed numerically accounting for evaporation with optical and material properties for ideal interstellar particles modeled. The properties of this kind of particle capture are discussed for the Sun but is applicable to other stars.

Jackson, A. A.↗

NASA NTRS1995

On the Calculation of Anisotropic Extinction Coefficients for Rigid Fibrous Ceramic Insulations

The specific anisotropic extinction coefficient e* couples the effective Radiative properties of a fibrous insulation into the radiation diffusion equation. This coefficient can be calculated using Mie scattering theory if fiber diameters, refractive indices and fiber orientation distributions are known. In general, fiber orientation distributions are not readily accessible and past calculations have considered fibers as either randomly distributed or normal to the heat flow direction. In certain rigid fibrous ceramic insulations neither of these cases apply well and a simple procedure is described for approximating e* from values calculated for the random and normal orientation cases. The intrinsic error associated with this scaling procedure is investigated. Numerical computations for several test structures and fiber materials show the average error to be less than 5% for net heat flux and radiation conductivity calculations.

Marshall, Jochen↗

NASA NTRS1997

Simulation of Ultra-Small Electronic Devices: The Classical-Quantum Transition Region

Concern is increasing about how quantum effects will impact electronic device operation as down-scaling continues along the SIA Roadmap through 2010. This document describes part of a new semiconductor device modeling (SDM) program at NAS to investigate these concerns by utilizing advanced NAS and third-party numerical computation software to rapidly implement and investigate electronic device models including quantum effects. This SDM project will investigate quantum effects in devices in the classical-quantum transition region, especially sub-0.1 mm MOSFETs. Specific tasks planned for this project include the use of quantum corrections to the classical drift-diffusion and hydrodynamic models of electron transport, arid the use of nominally quantum models including significant scattering.

Biegel, Bryan A.↗

NASA NTRS2002

Surface Enhancement Improves Crack Resistance

The low plasticity burnishing (LPB) process produces a deep layer of surface compression in a quick and affordable manner to produce metal surfaces free of scratches, nicks, and gouges. The process, designed for easy inclusion in the manufacturing environment, can be performed with conventional Computer Numerical Control machine tools. This allows parts to be processed during manufacturing, rather than as a post process in a separate facility. A smooth, free-rolling spherical ball suspended in a fluid allows for single-point contact. The ball comes into mechanical contact only with the surface to be burnished, and can be moved in any direction. LPB can be applied to all types of carbon and alloy steel, stainless steel, cast iron, aluminum, titanium, and nickel- based super alloys. In addition to improving a surface's resistance to fatigue and damage, treatment stops the growth of shallow cracks. The LPB process is used on the leading edges of turbine blades to improve resistance to foreign object damage and crack growth. This means significant savings for aircraft owners, since maintenance requirements to inspect for fatigue damage, replace parts, and remove corrosion damage increase the cost of operation.

Source record↗

NASA NTRS1999

An O(Nm(sup 2)) Plane Solver for the Compressible Navier-Stokes Equations

A hierarchical multigrid algorithm for efficient steady solutions to the two-dimensional compressible Navier-Stokes equations is developed and demonstrated. The algorithm applies multigrid in two ways: a Full Approximation Scheme (FAS) for a nonlinear residual equation and a Correction Scheme (CS) for a linearized defect correction implicit equation. Multigrid analyses which include the effect of boundary conditions in one direction are used to estimate the convergence rate of the algorithm for a model convection equation. Three alternating-line- implicit algorithms are compared in terms of efficiency. The analyses indicate that full multigrid efficiency is not attained in the general case; the number of cycles to attain convergence is dependent on the mesh density for high-frequency cross-stream variations. However, the dependence is reasonably small and fast convergence is eventually attained for any given frequency with either the FAS or the CS scheme alone. The paper summarizes numerical computations for which convergence has been attained to within truncation error in a few multigrid cycles for both inviscid and viscous ow simulations on highly stretched meshes.

Thomas, J. L.↗

NASA NTRS2004

Numerical Simulation of 3-D Supersonic Viscous Flow in an Experimental MHD Channel

The 3-D supersonic viscous flow in an experimental MHD channel has been numerically simulated. The experimental MHD channel is currently in operation at NASA Ames Research Center. The channel contains a nozzle section, a center section, and an accelerator section where magnetic and electric fields can be imposed on the flow. In recent tests, velocity increases of up to 40% have been achieved in the accelerator section. The flow in the channel is numerically computed using a new 3-D parabolized Navier-Stokes (PNS) algorithm that has been developed to efficiently compute MHD flows in the low magnetic Reynolds number regime. The MHD effects are modeled by introducing source terms into the PNS equations which can then be solved in a very e5uent manner. To account for upstream (elliptic) effects, the flowfield can be computed using multiple streamwise sweeps with an iterated PNS algorithm. The new algorithm has been used to compute two test cases that match the experimental conditions. In both cases, magnetic and electric fields are applied to the flow. The computed results are in good agreement with the available experimental data.

Kato, Hiromasa↗

NASA NTRS2004

Simulation of 3-D Nonequilibrium Seeded Air Flow in the NASA-Ames MHD Channel

The 3-D nonequilibrium seeded air flow in the NASA-Ames experimental MHD channel has been numerically simulated. The channel contains a nozzle section, a center section, and an accelerator section where magnetic and electric fields can be imposed on the flow. In recent tests, velocity increases of up to 40% have been achieved in the accelerator section. The flow in the channel is numerically computed us ing a 3-D parabolized Navier-Stokes (PNS) algorithm that has been developed to efficiently compute MHD flows in the low magnetic Reynolds number regime: The MHD effects are modeled by introducing source terms into the PNS equations which can then be solved in a very efficient manner. The algorithm has been extended in the present study to account for nonequilibrium seeded air flows. The electrical conductivity of the flow is determined using the program of Park. The new algorithm has been used to compute two test cases that match the experimental conditions. In both cases, magnetic and electric fields are applied to the seeded flow. The computed results are in good agreement with the experimental data.

Gupta, Sumeet↗

NASA NTRS1995

Effect of Far-Field Boundary Conditions on Boundary-Layer Transition

The effect of far-field boundary conditions on the evolution of a finite-amplitude two-dimensional wave in the Blasius boundary layer is assessed. With the use of the parabolized stability equations (PSE) theory for the numerical computations, either asymptotic, Dirichlet, Neumann or mixed boundary conditions are imposed at various distances from the wall. The results indicate that asymptotic and mixed boundary conditions yield the most accurate mean-flow distortion and unsteady instability modes in comparison with the results obtained with either Dirichlet or Neumann conditions.

Bertolotti, Fabio P.↗

NASA NTRS1994

Effect of Far-Field Boundary Conditions on Boundary-Layer Transition

Bertolotti, Fabio P.↗

NASA NTRS2004

Evaluation of CFD to Determine Two-Dimensional Airfoil Characteristics for Rotorcraft Applications

The efficient prediction of helicopter rotor performance, vibratory loads, and aeroelastic properties still relies heavily on the use of comprehensive analysis codes by the rotorcraft industry. These comprehensive codes utilize look-up tables to provide two-dimensional aerodynamic characteristics. Typically these tables are comprised of a combination of wind tunnel data, empirical data and numerical analyses. The potential to rely more heavily on numerical computations based on Computational Fluid Dynamics (CFD) simulations has become more of a reality with the advent of faster computers and more sophisticated physical models. The ability of five different CFD codes applied independently to predict the lift, drag and pitching moments of rotor airfoils is examined for the SC1095 airfoil, which is utilized in the UH-60A main rotor. Extensive comparisons with the results of ten wind tunnel tests are performed. These CFD computations are found to be as good as experimental data in predicting many of the aerodynamic performance characteristics. Four turbulence models were examined (Baldwin-Lomax, Spalart-Allmaras, Menter SST, and k-omega).

Smith, Marilyn J.↗

NASA NTRS2004

A Survey Of Earth-Moon Libration Orbits: Stationkeeping Strategies And Intra-Orbit Transfers

Cislunar space is a readily accessible region that may well develop into a prime staging area in the effort to colonize space near Earth or to colonize the Moon. While there have been statements made by various NASA programs regarding placement of resources in orbit about the Earth-Moon Lagrangian locations, there is no survey of the total cost associated with attaining and maintaining these unique orbits in an operational fashion. Transfer trajectories between these orbits required for assembly, servicing, and positioning of these resources have not been extensively investigated. These orbits are dynamically similar to those used for the Sun-Earth missions, but differences in governing gravitational ratios and perturbation sources result in unique characteristics. We implement numerical computations using high fidelity models and linear and nonlinear targeting techniques to compute the various maneuver (Delta)V and temporal costs associated with orbits about each of the Earth-Moon Lagrangian locations (L1, L2, L3, L4, and L5). From a dynamical system standpoint, we speak to the nature of these orbits and their stability. We address the cost of transfers between each pair of Lagrangian locations.

Folta, David↗

NASA NTRS2004

Comparisons of the Solutions for the Category 3--Problem 2: Cascade-Gust Interaction

Six different solutions were submitted for this benchmark problem. These were obtained using a variety of methods that can be conveniently categorized in two main groups, a nonlinear time-domain group and a linearized frequency-domain group. The first includes solutions submitted by (1) Hixon, (2) Nallasamy et. al, (3) Shieh et. al, and (4) Wang et. al, and the second includes solutions submitted by (5) Coupland and (6) Serrano et. al. Methods (1) and (2) use sixth order compact differencing schemes and the rest are essentially second order in space. With the exception of the solution submitted by Shieh et. al, all are individually discussed in great detail in the workshop proceedings. Comparisons of the submitted solutions with the benchmark solution are presented below. Due to differences in the level of solution detail provided to the author by the participants, the comparisons do not always include results from all submissions. It should be noted at the outset that, since the benchmark solution itself was numerically computed, the comparisons are somewhat subjective. In order to provide maximum latitude for the participants of the workshop, no restrictions were placed on the type of method that could be used to solve the problem. Neither was there were any stipulations to use a particular grid topology or grid density. Therefore, without a detailed study of the critical features of the computed solutions, it is not possible to make concrete statements about the relative merits of one method over another. Such a study is beyond the scope of the current exercise, especially since complete flowfield details were not provided to the author by all participants. Instead a package, containing the information about the benchmark solution (both the steady and unsteady parts of it), is included on the proceedings CD should the authors who submitted solutions for this problem wish to examine in detail the benchmark solution and compare their results to it.

Envia, Edmane↗

NASA NTRS2005

Fully Associative, Nonisothermal, Potential-Based Unified Viscoplastic Model for Titanium-Based Matrices

A number of titanium matrix composite (TMC) systems are currently being investigated for high-temperature air frame and propulsion system applications. As a result, numerous computational methodologies for predicting both deformation and life for this class of materials are under development. An integral part of these methodologies is an accurate and computationally efficient constitutive model for the metallic matrix constituent. Furthermore, because these systems are designed to operate at elevated temperatures, the required constitutive models must account for both time-dependent and time-independent deformations. To accomplish this, the NASA Lewis Research Center is employing a recently developed, complete, potential-based framework. This framework, which utilizes internal state variables, was put forth for the derivation of reversible and irreversible constitutive equations. The framework, and consequently the resulting constitutive model, is termed complete because the existence of the total (integrated) form of the Gibbs complementary free energy and complementary dissipation potentials are assumed a priori. The specific forms selected here for both the Gibbs and complementary dissipation potentials result in a fully associative, multiaxial, nonisothermal, unified viscoplastic model with nonlinear kinematic hardening. This model constitutes one of many models in the Generalized Viscoplasticity with Potential Structure (GVIPS) class of inelastic constitutive equations.

Source record↗

NASA NTRS2002

Comparison Made of Operating Characteristics of Spiral Bevel Gears Manufactured Using Different Methods

Spiral bevel gears are important components on all current rotorcraft drive systems. These components are required to operate at high speeds, high loads, and for an extremely large number of load cycles. In this application, spiral bevel gears are used to redirect the shaft from the horizontal gas turbine engine to the vertical rotor. Because of the high expense of manufacturing these gears, methods that can achieve the same level of performance at reduced cost are highly desirable to aerospace gear manufacturers. Gears manufactured for aerospace applications use high-quality materials and are manufactured to tight tolerances. Special manufacturing machine tools and computer numerically controlled coordinate measurement systems have enabled rotorcraft drive system manufacturers to produce extremely high-quality gears during their normal production. Because of low production rates for rotorcraft, these gears are manufactured in small batches, and thus are unable to benefit from the economics of high production numbers as in other industries. In this investigation, two different manufacturing methods, face-milled and face-hobbed, were used to fabricate spiral bevel gears. For face-milled spiral bevel gears, grinding of the contacting surfaces is the final manufacturing step. At least two different specialty machines are needed to generate the teeth for face-milled spiral bevel gears. For face-hobbed gears, hard cutting is the final manufacturing process. The same machine is used to rough cut and finish cut the gears. This study compared the operational behavior of face-milled spiral bevel gears with that of face-hobbed spiral bevel gears. Test hardware was manufactured to fit within NASA Glenn Research Center's Spiral Bevel Test Facility and to aerospace quality standards. Tests were conducted for stress, vibration, and noise. A comparison of the results attained indicated that the face-hobbed gears had a lower alternating stress level with a more even distribution of loading across the teeth, and slightly reduced levels of vibration and noise. Results of this study show that the face-hobbed method is a viable and lower-cost alternative for producing aerospace-quality spiral-bevel gears.

Handschuh, Robert F.↗

NASA NTRS2003

Novel Engineering and Fabrication Techniques Tested in Low-Noise- Research Fan Blades

A major source of fan noise in commercial turbofan engines is the interaction of the wake from the fan blades with the stationary vanes (stators) directly behind them. The Trailing Edge Blowing (TEB) project team at the NASA Glenn Research Center designed and fabricated new fan blades to study the effects of fan trailing edge blowing as a potential noise-reduction concept. The intent is to fill the rotor wake by supplying air to the rotor blade trailing edge at the proper conditions to minimize the wake deficit, and thus generate less noise. The TEB hardware is designed for the Active Noise Control Fan (ANCF) test rig in Glenn's Aeroacoustic Propulsion Laboratory. For this test, the air is fed from an external supply through the shaft of the rig. It is distributed to the base of each blade through an impeller, where it is forced into a plenum at the core of each blade. In actual engine configuration, air would most likely be bled from the compressor, but only at times when noise is an issue, such as takeoffs and landings. Glenn researchers designed and manufactured the blades in-house, using new techniques and concepts. The skins, which were designed for maximum strength in the directions of highest stress, were molded from multiple layers of carbon fiber. Considerable use was made of rapid prototyping techniques, such as laser sintering. The core was sintered from a lightweight polymer, and the retainer was CNC-machined (computer numerical control machined) from aluminum. All the components were joined with a cold-cure aerospace adhesive. These techniques and processes reduced the overall cost and allowed the new concept to be studied much sooner than would be possible using traditional fabrication methods. Since this test rig did not support the use of blade-monitoring techniques such as strain gauges, extensive bench testing was required to qualify the design. The blades were examined using a variety of methods including holography, pull tests (cyclic and failure), shake tests, rap tests, and nondestructive inspection. Acoustic testing of the ANCF fan using TEB has been ongoing since January of 2001. The fan has completed about 100 hr of testing with no structural, vibrational, or fatigue problems. Far-field acoustic measurements, in-duct mode measurements, precise hot wire surveys, and detailed performance measurements are providing data for evaluating the concept. The far-field noise data show that tone noise was reduced significantly with the initial ANCF TEB fan design. In addition, a significant reduction in unsteady stator loading has been measured, indicating the potential for stator broadband noise reduction. The acoustic benefits will be assessed and important design parameters identified to improve the ability to fully exploit any benefit provided by this technique. On the basis of the success of trailing edge blowing, Glenn plans to continue this research with a higher speed, higher pressure ratio fan operating in an acoustic wind tunnel to simulate flight conditions.

Cunningham, Cameron C.↗

NASA NTRS2006

Large Array Channel Capacity in the Presence of Interference

We develop a model for a large array ground receiver system for use in deep-space communications, and analyze the resulting array channel capacity. The model includes effects of array geometry, time-dependent spacecraft orbital trajectory, point and extended interference sources, and elevation-dependent noise and tropospheric channel variations. Channel capacity is expressed as the ratio of determinants of covariance matrices characterizing source, interference, and additive noise, and then reduced to a simpler quadratic form more amenable to analysis and numerical computation. This formulation facilitates inclusion of array and channel characteristics into the model, as well as comparison of optimal, suboptimal, and equivalent single antenna configurations on achievable throughput. Realistic examples of ground array channel capacity calculations are presented, demonstrating the impact of array geometry, planetary interference sources, and array combining algorithm design upon the achievable data throughput.

Vilnrotter, V.↗