Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Methods”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Modeling of Convective-Stratiform Precipitation Processes: Sensitivity to Partitioning Methods and Numerical Advection Schemes

Six different convective-stratiform separation techniques, including a new technique that utilizes the ratio of vertical and terminal velocities, are compared and evaluated using two-dimensional numerical simulations of a tropical [Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional Experiment for STORM-Central (PRESTORM)] squall line. The simulations are made using two different numerical advection schemes: 4th order and positive definite advection. Comparisons are made in terms of rainfall, cloud coverage, mass fluxes, apparent heating and moistening, mean hydrometeor profiles, CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent heating retrieval. Overall, it was found that the different separation techniques produced results that qualitatively agreed. However, the quantitative differences were significant. Observational comparisons were unable to conclusively evaluate the performance of the techniques. Latent heating retrieval was shown to be sensitive to the use of separation technique mainly due to the stratiform region for methods that found very little stratiform rain. The midlatitude PRESTORM simulation was found to be nearly invariant with respect to advection type for most quantities while for TOGA COARE fourth order advection produced numerous shallow convective cores and positive definite advection fewer cells that were both broader and deeper penetrating above the freezing level.

Lang, Steve↗
NASA NTRS1992

Navier-Stokes analysis of three-dimensional unsteady flows inside turbine stages

This study presents a numerical method for solving the 3D Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully 3D Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms, while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flowfield of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flowfields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗
NASA NTRS2010

Faster and More Accurate Transport Procedures for HZETRN

Several aspects of code verification are examined for HZETRN. First, a detailed derivation of the numerical marching algorithms is given. Next, a new numerical method for light particle transport is presented, and improvements to the heavy ion transport algorithm are discussed. A summary of various coding errors is also given, and the impact of these errors on exposure quantities is shown. Finally, a coupled convergence study is conducted. From this study, it is shown that past efforts in quantifying the numerical error in HZETRN were hindered by single precision calculations and computational resources. It is also determined that almost all of the discretization error in HZETRN is caused by charged target fragments below 50 AMeV. Total discretization errors are given for the old and new algorithms, and the improved accuracy of the new numerical methods is demonstrated. Run time comparisons are given for three applications in which HZETRN is commonly used. The new algorithms are found to be almost 100 times faster for solar particle event simulations and almost 10 times faster for galactic cosmic ray simulations.

Slaba, Tony C.↗
NASA NTRS1993

Evaluation of advanced materials through experimental mechanics and modelling

Composite materials have been frequently used in aerospace vehicles. Very often defects are inherited during the manufacture and damages are inherited during the construction and services. It becomes critical to understand the mechanical behavior of such composite structure before it can be further used. One good example of these composite structures is the cylindrical bottle of solid rocket motor case with accidental impact damages. Since the replacement of this cylindrical bottle is expensive, it is valuable to know how the damages affects the material, and how it can be repaired. To reach this goal, the damage must be characterized and the stress/strain field must be carefully analyzed. First the damage area, due to impact, is surveyed and identified with a shearography technique which uses the principle of speckle shearing interferometry to measure displacement gradient. Within the damage area of a composite laminate, such as the bottle of solid rocket motor case, all layers are considered to be degraded. Once a lamina being degraded the stiffness as well as strength will be drastically decreased. It becomes a critical area of failure to the whole bottle. And hence the stress/strain field within and around a damage should be accurately evaluated for failure prediction. To investigate the stress/strain field around damages a Hybrid-Numerical method which combines experimental measurement and finite element analysis is used. It is known the stress or strain at the singular point can not be accurately measured by an experimental technique. Nevertheless, if the location is far away from the singular spot, the displacement can be found accurately. Since it reflects the true displacement field locally regardless of the boundary conditions, it is an excellent input data for a finite element analysis to replace the usually assumed boundary conditions. Therefore, the Hybrid-Numerical method is chosen to avoid the difficulty and to take advantage of both experimental technique and finite element analysis. Experimentally, the digital image correlation technique is employed to measure the displacement field. It is done by comparing two digitized images, before and after loading. Numerically, the finite element program, ABAQUS (version 5.2), is used to analyze the stress and strain field. It takes advantage of the high speed and huge memory size of modern supercomputer, CRAY Y-MP, at NASA Marshall Space Flight Center.

Yang, Yii-Ching↗
NASA NTRS1993

Navier-Stokes analysis of three-dimensional flow and heat transfer inside turbine blade rows

A numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow and heat transfer through multiple turbomachinery blade rows is presented. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow and heat transfer field of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗
NASA NTRS1984

Theoretical and experimental engine-inlet flow fields for fighter forebodies

The capability of two numerical methods, one for transonic and one for supersonic flows, to predict the flow fields about representative fighter aircraft forebodies in the vicinity of the engine inlets was examined. The Mach number range covered was 0.9 to 2.5 and the angle-of-attack range was 0 deg to 25 deg. The computer progams that implement each of the numerical methods are described as to their features and usage, and results are compared with comprehensive wind tunnel data. Although both prediction methods were inviscid, results show that the aerodynamic effects of the forebody, with and without a wing, can be simulated fairly well. Futher work is needed to include the effects of viscosity, including vortex shedding.

Yaros, S. F.↗
NASA NTRS1992

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗
NASA NTRS1976

The axisymmetric jet counterflow problem

A numerical method for computing the flow field of an underexpanded axisymmetric supersonic jet flowing counter to a uniform supersonic stream is described. The strong shock waves and shear layers which occur in the flow are treated as normal and tangential discontinuities bounding regions of inviscid flow. Within the regions, the unsteady gasdynamic equations are solved in integral form on a moving computational mesh. As the numerical flow field evolves, the mesh is deformed according to the jump relations applied at the discontinuities, thus determining their positions as part of the solution. Also presented are computed results obtained for a wide range of jet total pressure. The standoff distances of the bow shock and jet shock are in good agreement with those measured experimentally, thus demonstrating the accuracy of the numerical method.

Schiff, L. B.↗
NASA NTRS1977

Efficient numerical treatment of periodic systems with application to stability problems

Two efficient numerical methods for dealing with the stability of linear periodic systems are presented. Both methods combine the use of multivariable Floquet-Liapunov theory with an efficient numerical scheme for computing the transition matrix at the end of one period. The numerical properties of these methods are illustrated by applying them to the simple parametric excitation problem of a fixed end column. The practical value of these methods is shown by applying them to some helicopter rotor blade aeroelastic and structural dynamics problems. It is concluded that these methods are numerically efficient, general and practical for dealing with the stability of large periodic systems.

Friedmann, P.↗
NASA NTRS1994

A new flux-conserving numerical scheme for the steady, incompressible Navier-Stokes equations

This paper is concerned with the continued development of a new numerical method, the space-time solution element (STS) method, for solving conservation laws. The present work focuses on the two-dimensional, steady, incompressible Navier-Stokes equations. Using first an integral approach, and then a differential approach, the discrete flux conservation equations presented in a recent paper are rederived. Here a simpler method for determining the flux expressions at cell interfaces is given; a systematic and rigorous derivation of the conditions used to simulate the differential form of the governing conservation law(s) is provided; necessary and sufficient conditions for a discrete approximation to satisfy a conservation law in E2 are derived; and an estimate of the local truncation error is given. A specific scheme is then constructed for the solution of the thin airfoil boundary layer problem. Numerical results are presented which demonstrate the ability of the scheme to accurately resolve the developing boundary layer and wake regions using grids which are much coarser than those employed by other numerical methods. It is shown that ten cells in the cross-stream direction are sufficient to accurately resolve the developing airfoil boundary layer.

Scott, James R.↗
NASA NTRS2000

Polar Plate Theory for Orthogonal Anisotropy

Laminated fiber-reinforced (or filamentary) composites are used today for their high strength-to-weight and stiffness-to-weight ratios. However, because of the anisotropic behavior of composites, determining the response on a macroscopic scale is challenging. This is particularly evident in the evaluation of the governing differential equations of a circular disk with the fibers of the lamina oriented with rectilinear orthogonality. This includes any situation involving a composite plate of circular geometry in which out-of-plane displacements due to load are desired, such as fastener pull through loading of a composite plate. Current analysis techniques use numerical methods with rectilinear coordinate systems to solve problems with circular geometry. These analyses over predict plate stiffness by 20% and underpredict failure by 70%. Consequently, there is a need to transform classical composite plate theory to a polar coordinate system. In order to better analyze structures with circular geometries the classical composite plate equations are transformed into the plate equations for a rectilinearly anisotropic composite in polar coordinates. A composite plate is typically a laminate of fibers in rectilinear directions. Subsequent to the lay-tip the necessary geometry is cut out of a rectangular plate. In a similar manner, the derivation of the plate equation starts with the fundamental definitions of strain, displacement and curvature and incorporates the material property angular dependence into the equilibrium equations for a differential polar element. In the transformed state, the stiffness coefficients are no longer constant, adding to the complexity of the governing differential equations. This paper discusses the new derivation and evaluation of the plate equations for a circular composite disk with orthogonal rectilinear anisotropy. The resultant new three partial differential equations, which describe the circular anisotropic plate, can be used to evaluate out-of-plane displacements for given load conditions in design of composite Structures. Without the formation of these unique plate equations, design of such structures is not a precise engineering accomplishment due to the lack of precise design tools. Such structures are over designed to compensate and costly tests need to be performed after production to validate safety. Though impossible to present the mathematical solution within the limits of this paper, the application of the newly derived plate equation in its orthotropic form for a balanced symmetric laminate compared favorably with test results and surpassed the results of standard numerical methods. Within the limitations of the problem foundation the predictive model provides a mechanism for establishing out-of-plane deflection levels for a circular composite plate. Future work in this area should include incorporating the circular composite plate equations in Finite Element Models. The new anisotropic polar plate equations call be utilized as a design tool so that the design more accurately meets the requirements, thus taking full advantage of the weight savings of composites. These accomplishments have not previously been presented to the engineering community nor ever proposed to the IAF. Should the proposed paper be accepted for presentation, attendance is assured.

Bailey, Michelle D.↗
NASA NTRS1997

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗
NASA NTRS1970

A method of numerically controlled machine part programming

Computer program is designed for automatically programmed tools. Preprocessor computes desired tool path and postprocessor computes actual commands causing machine tool to follow specific path. It is used on a Cincinnati ATC-430 numerically controlled machine tool.

Source record↗
NASA NTRS1971

The calculation of electromagnetic fields in the Fresnel and Fraunhofer regions using numerical integration methods

Some results obtained with a digital computer program written at Goddard Space Flight Center to obtain electromagnetic fields scattered by perfectly reflecting surfaces are presented. For purposes of illustration a paraboloidal reflector was illuminated at radio frequencies in the simulation for both receiving and transmitting modes of operation. Fields were computed in the Fresnel and Fraunhofer regions. A dual-reflector system (Cassegrain) was also simulated for the transmitting case, and fields were computed in the Fraunhofer region. Appended results include derivations which show that the vector Kirchhoff-Kottler formulation has an equivalent form requiring only incident magnetic fields as a driving function. Satisfaction of the radiation conditions at infinity by the equivalent form is demonstrated by a conversion from Cartesian to spherical vector operators. A subsequent development presents the formulation by which Fresnel or Fraunhofer patterns are obtainable for dual-reflector systems. A discussion of the time-average Poynting vector is also appended.

Schmidt, R. F.↗
NASA NTRS1978

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗
NASA NTRS1981

The experimental verification of a streamline curvature numerical analysis method applied to the flow through an axial flow fan

The experimental verification of an inviscid, incompressible through-flow analysis method is presented. The primary component of this method is an axisymmetric streamline curvature technique which is used to compute the hub-to-tip flow field of a given turbomachine. To analyze the flow field in the blade-to-blade plane of the machine, the potential flow solution of an infinite cascade of airfoils is also computed using a source model technique. To verify the accuracy of such an analysis method an extensive experimental verification investigation was conducted using an axial flow research fan. Detailed surveys of the blade-free regions of the machine along with intra-blade surveys using rotating pressure sensing probes and blade surface static pressure taps provide a one-to-one relationship between measured and predicted data. The results of this investigation indicate the ability of this inviscid analysis method to predict the design flow field of the axial flow fan test rotor to within a few percent of the measured values.

Pierzga, M. J.↗