Search NASASearch

SEARCH · Search NASA

Results for “Numerical analysis”

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 199 records · Page 11

Numerical analysis of three-dimensional viscous internal flows

A 3-D Navier-Stokes code has been developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system rotating about the x-axis, and then mapped to a general body-fitted coordinate system. Streamwise viscous terms are neglected using the thin-layer assumption, and turbulence effects are modeled using the Baldwin-Lomax turbulence model. The equations are discretized using finite differences on stacked C-type grids and are solved using a multistage Runge-Kutta algorithm with a spatially-varying time step and implicit residual smoothing. Calculations have been made of a horseshoe vortex formed in front of a flat plate with a round leading edge standing in a turbulent endwall boundary layer. Comparisons are made with experimental data taken by Eckerle and Langston for a circular cylinder under similar conditions. Computer and measured results are compared in terms of endwall flow visualization pictures and total pressure loss contours and vector plots on the symmetry plane. Calculated details of the primary vortex show excellent agreement with the experimental data. The calculations also show a small secondary vortex not seen experimentally.

Chima, Rodrick V.

Numerical analysis of pulse facility transient flow through a supersonic combustion model

The supersonic high enthalpy flow from a shock tunnel facility flowing through a two-dimensional duct simulating a combustor has been computed for two test conditions for both steady and transient inflow. The nominal test conditions had stagnation enthalpies of 4.3 and 12.1 MJ/kg, with corresponding flight Mach numbers of 9 and 16, respectively. A full Navier-Stokes code (SPARK) was used to compute the flow development and the mixing of a hydrogen jet injected from the base of a central strut. Results of the calculations are presented and compared with wall pressure data. Additional comparisons of the transient and steady inflow results are presented to assess the effect of transient flow on the data. This paper presents a status of the progress in applying CFD codes to analyze the transient flow from a pulse facility.

Rogers, R. Clayton

Numerical analysis of three-dimensional viscous internal flows

A 3-D Navier-Stokes code has been developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system rotating about the x-axis, and then mapped to a general body-fitted coordinate system. Streamwise viscous terms are neglected using the thin-layer assumption, and turbulence effects are modeled using the Baldwin-Lomax turbulence model. The equations are discretized using finite differences on stacked C-type grids and are solved using a multistage Runge-Kutta algorithm with a spatially-varying time step and implicit residual smoothing. Calculations have been made of a horseshoe vortex formed in front of a flat plate with a round leading edge standing in a turbulent endwall boundary layer. Comparisons are made with experimental data taken by Eckerle and Langston for a circular cylinder under similar conditions. Computer and measured results are compared in terms of endwall flow visualization pictures and total pressure loss contours and vector plots on the symmetry plane. Calculated details of the primary vortex show excellent agreement with the experimental data. The calculations also show a small secondary vortex not seen experimentally.

Chima, Rodrick V.

Numerical analysis of flow about a total temperature sensor

The unsteady flowfield about an airfoil-shaped inlet temperature sensor has been investigated using the thin-layer and full Navier-Stokes equations. A finite-volume formulation of the governing equations was used in conjunction with a Runge-Kutta time stepping scheme to analyze the flow about the sensor. Flow characteristics for this configuration were established at Mach numbers of 0.5 and 0.8 for different Reynolds numbers. The results were obtained for configurations of increasing complexity; important physical phenomena such as shock formation, boundary-layer separation, and unsteady wake formation were noted. Based on the computational results, recommendations for further study and refinement of the inlet temperature sensor were made.

Von Lavante, Ernst

A numerical analysis of the effects of conjugate heat transfer, vapor compressibility, and viscous dissipation in heat pipes

The effects of conjugate heat transfer, vapor compressibility, and viscous dissipation in heat pipes are discussed. The accuracy of the partially parabolic versus the elliptic presentation of the governing equations is also examined. The results show that the axial wall conduction has a tendency to make the temperature distribution more uniform for heat pipes with large ratios of pipe wall to effective liquid-wick thermal conductivity. The compressible and incompressible models show very close agreement for the total pressure drop, while the local pressure variations along the heat pipe are quite different for these two models when the radial Reynolds number at the interface is high.

Faghri, Amir

Numerical analysis of ignition transients in solid rocket motors

A model to analyze the unstready multidimensional turbulent flow in the head-end, star slot section of a solid rocket motor during the ignition transient is developed. The present paper examines the fluid dynamic aspects of the starting transient. A finite-difference solution of the unsteady, compressible, full Navier-Stokes equations, together with a two-equation k-epsilon turbulence model, is obtained utilizing the MacCormack explicit, predictor-corrector technique. Computed results for the flow in the slot are compared with experimental data obtained in cold flow tests of a scaled model of the Space Shuttle SRM head-end, star grain section. The agreement between the calculated flowfield and existing experimental data is very good.

Ciucci, A.

Numerical analysis and simulation of an assured crew return vehicle flow field

A lifting body was proposed as a candidate for the Assured Crew Return Vehicle (ACRV) which will serve as a crew rescue vehicle for the Space Station Freedom. The focus is on body surface definition, both surface and volume grid definition, and the computation of inviscid flow fields about the vehicle at wind tunnel conditions. Very good agreement is shown between the computed aerodynamic characteristics of the vehicle at M(sub infinity) = 10 and those measured in wind tunnel tests at high Reynolds numbers.

Weilmuenster, K. James

Numerical analysis of the Iosipescu specimen for composite materials

A finite element analysis of the Iosipescu shear tests for unidirectional and cross-ply composites is presented. It is shown that an iterative analysis procedure must be used to model the fixture-specimen kinematics. The correction factors which are needed to compensate for the nonuniformity of stress distribution in calculating shear modulus are shown to be dependent on the material orthotropic ratio and the finite element loading models. Test section strain distributions representative of typical graphite-epoxy specimens are also presented.

Ho, H.

Numerical analysis of bipropellant combustion in liquid thrust chambers by an Eulerian-Eulerian approach

The liquid thrust chambers performance (LTCP) code is used for parametric studies of flow and combustion in liquid rocket engines. Multiphase flow equations are solved in an Eulerian-Eulerian framework, and multistep finite rate chemistry is incorporated. The discretization scheme is fully implicit and is based on the total variation diminishing (TVD) scheme, which is accurate, robust, very efficient and capable of handling steep gradients and stiff chemistry. Effects of injection velocity and chamber size have been considered, and the effect of group combustion on the evaporation rate has been studied for a dense spray.

Dang, A. L.

Computational and numerical analysis of hypersonic nozzle flows with comparisons to wind tunnel calibration data

Navier-Stokes (NS), parabolized NS codes, and classical method of characteristics prediction techniques are employed to compute hypersonic nozzle flowfields for the 15 in Mach 6 High Temperature Tunnel and the 16 in Mach 17 Nitrogen Tunnel, which are part of the Langley Research Center. The study focuses on defining the accuracy of these techniques by comparing the computational results to wind tunnel pitot pressure measurements conducted in the newly designed nozzle of the Mach 6 tunnel and in the original nozzle of the Mach 17 tunnel. The comparisons between predicted results and test section pitot surveys for the Mach 6 nozzle are in good agreement and show highly uniform flow over a range of reservoir pressures and temperatures.

Hackett, Charles M.

Numerical analysis of right-half plane zeros for a single-link manipulator

The purpose of this research is to further develop an understanding of how nonminimum phase zero location is affected by structural link design. As the demand for light-weight robots that can operate in a large workspace increases, the structural flexibility of the links become more of an issue in controls problems. When the objective is to accurately position the tip while the robot is actuated at the base, the system is nonminimum phase. One important characteristic of nonminimum phase systems is system zeros in the right half of the Laplace plane. The ability to pick the location of these nonminimum phase zeros would give the designer a new freedom similar to pole placement. The research targets a single-link manipulator operating in the horizontal plane and modeled as a Euler-Bernoulli beam with pinned-free end conditions. Using transfer matrix theory, one can consider link designs that have variable cross-sections along the length of the beam. A FORTRAN program was developed to determine the location of poles and zeros given the system model. The program was used to confirm previous research on nonminimum phase systems, and develop a relationship for designing linearly tapered links. The method allows the designer to choose the location of the first pole and zero and then defines the appropriate taper to match the desired locations. With the pole and zero location fixes, the designer can independently change the link's moment of inertia about its axis of rotation by adjusting the height of the beam. These results can be applied to inverse dynamic algorithms currently under development at Georgia Tech.

Girvin, Douglas Lynn

Numerical analysis of nonminimum phase zero for nonuniform link design

As the demand for light-weight robots that can operate in a large workspace increases, the structural flexibility of the links becomes more of an issue in control. When the objective is to accurately position the tip while the robot is actuated at the base, the system is nonminimum phase. One important characteristic of nonminimum phase systems is system zeros in the right half of the Laplace plane. The ability to pick the location of these nonminimum phase zeros would give the designer a new freedom similar to pole placement. This research targets a single-link manipulator operating in the horizontal plane and modeled as a Euler-Bernoulli beam with pinned-free end conditions. Using transfer matrix theory, one can consider link designs that have variable cross-sections along the length of the beam. A FORTRAN program was developed to determine the location of poles and zeros given the system model. The program was used to confirm previous research on nonminimum phase systems, and develop a relationship for designing linearly tapered links. The method allows the designer to choose the location of the first pole and zero and then defines the appropriate taper to match the desired locations. With the pole and zero location fixed, the designer can independently change the link's moment of inertia about its axis of rotation by adjusting the height of the beam. These results can be applied to the inverse dynamic algorithms that are currently under development.

Girvin, Douglas L.

Experimental and numerical analysis of the wing rock characteristics of a 'wing-body-tail' configuration

Free-to-roll wind tunnel tests were conducted and a computer simulation exercise was performed in an effort to investigate in detail the mechanism of wing rock on a configuration that consisted of a highly-slender forebody and a 78 deg swept delta wing. In the wind tunnel test, the roll angle and wing surface pressures were measured during the wing rock motion. A limit cycle oscillation was observed for angles of attack between 22 deg and 30 deg. In general, the wind tunnel test confirmed that the main flow phenomena responsible for the wing-body-tail wing rock are the interactions between the forebody and the wing vortices. The variation of roll acceleration (determined from the second derivative of the roll angle time history) with roll angle clearly showed the energy balance necessary to sustain the limit cycle oscillation. Pressure measurements on the wing revealed the hysteresis of the wing rock process. First, second and nth order models for the aerodynamic damping were developed and examined with a one degree of freedom computer simulation. Very good agreement with the observed behavior from the wind tunnel was obtained.

Suarez, Carlos J.

Numerical analysis of the Iosipescu specimen for composite materials

A finite element analysis of the Iosipescu shear tests for unidirectional and cross-ply composites is presented. It is shown that an iterative analysis procedure must be used to model the fixture-specimen kinematics. The correction factors which are needed to compensate for the nonuniformity of stress distribution in calculating shear modulus are shown to be dependent on the material orthotropic ratio and the finite element loading models. Test section strain distributions representative of typical graphite-epoxy specimens are also presented.

Ho, H.

Numerical analysis of the sensitivity of crystal growth experiments to spacecraft residual acceleration

An analysis is conducted of the sensitivity of the Bridgman-Stockbarger crystal growth method, using an idealized model for a range of operating and boundary conditions over a variety of accelerations. Attention is given to the dopant nonuniformity at the melt-crystal interface. The largest compositional nonuniformities are found to occur for disturbances whose amplitudes are greater than 10 exp 6 g, and frequencies below 0.1 Hz.

Alexander, J. I. D.