SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Combustion hot section technology

The overall objective of the Turbine Engine Hot Section Technology Combustion Project is to develop and verify improved and more accurate analysis methods for increasing the ability to design with confidence the combustion system for advanced aircraft turbine engines. The analysis methods developed will be generically applicable to combustion systems and not restricted to one specific engine or manufacturer. This project's approach was to first assess and evaluate existing combustor aerothermal analysis models by means of a contracted effort initiated during FY 1982. This evaluation effort has assessed and quantified known models' strengths and deficiencies. During FY 1984 the Aerothermal Modeling Program, Phase 2 will be initiated, which is expected to have contracted model development efforts in the areas of improved numerical methods for turbulent viscous flows, flow interactions, and fuel spray flow foekd interactions. A Phase 3 effort is planned to address remaining model deficiencies. The primary inhouse effort in this area will be the determination of high pressure flame radiation characteristics in a full annular combustor. This experiment will be conducted in the NASA LeRC High Pressure Facility with the results compiled into a comprehensive flame radiation and liner heat flux model.

Ercegovic, D. B.↗

NASA NTRS1987

Guidance and control strategies for aerospace vehicles

A simplified method of matched asymptotic expansions was developed where the common part in composite solution is generated as a polynomial in stretched variable instead of actually evaluating the same from the outer solution. This methodology was applied to the solution of the exact equations for three dimensional atmospheric entry problems. Compared to previous works, the present simplified methodology yields explicit analytical expressions for various components of the composite solution without resorting to any type of transcendental equations to be solved only by numerical methods. The optimal control problem arising in the noncoplanar orbital transfer employing aeroassist was also addressed.

Hibey, Joseph L.↗

NASA NTRS1996

Development of the Tensoral Computer Language

The research scientist or engineer wishing to perform large scale simulations or to extract useful information from existing databases is required to have expertise in the details of the particular database, the numerical methods and the computer architecture to be used. This poses a significant practical barrier to the use of simulation data. The goal of this research was to develop a high-level computer language called Tensoral, designed to remove this barrier. The Tensoral language provides a framework in which efficient generic data manipulations can be easily coded and implemented. First of all, Tensoral is general. The fundamental objects in Tensoral represent tensor fields and the operators that act on them. The numerical implementation of these tensors and operators is completely and flexibly programmable. New mathematical constructs and operators can be easily added to the Tensoral system. Tensoral is compatible with existing languages. Tensoral tensor operations co-exist in a natural way with a host language, which may be any sufficiently powerful computer language such as Fortran, C, or Vectoral. Tensoral is very-high-level. Tensor operations in Tensoral typically act on entire databases (i.e., arrays) at one time and may, therefore, correspond to many lines of code in a conventional language. Tensoral is efficient. Tensoral is a compiled language. Database manipulations are simplified optimized and scheduled by the compiler eventually resulting in efficient machine code to implement them.

Ferziger, Joel↗

NASA NTRS2004

Computational Challenges of Viscous Incompressible Flows

Over the past thirty years, numerical methods and simulation tools for incompressible flows have been advanced as a subset of the computational fluid dynamics (CFD) discipline. Although incompressible flows are encountered in many areas of engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to the rather stringent requirements for predicting aerodynamic performance characteristics of flight vehicles, while flow devices involving low-speed or incompressible flow could be reasonably well designed without resorting to accurate numerical simulations. As flow devices are required to be more sophisticated and highly efficient CFD took become increasingly important in fluid engineering for incompressible and low-speed flow. This paper reviews some of the successes made possible by advances in computational technologies during the same period, and discusses some of the current challenges faced in computing incompressible flows.

Kwak, Dochan↗

NASA NTRS1976

Simplified solution for point contact deformation between two elastic solids

A linear-regression by the method of least squares is made on the geometric variables that occur in the equation for point contact deformation. The ellipticity and the complete eliptic integrals of the first and second kind are expressed as a function of the x, y-plane principal radii. The ellipticity was varied from 1 (circular contact) to 10 (a configuration approaching line contact). These simplified equations enable one to calculate easily the point-contact deformation to within 3 percent without resorting to charts or numerical methods.

Brewe, D. E.↗

NASA NTRS1980

Numerical calculation of transonic axial turbomachinery flows

A numerical method and the results of a computer program are presented for solving an exact, three dimensional, full potential equation that models rotating and nonrotating inviscid, absolutely irrotational, homentropic flows. Besides calculating the flows through an arbitrarily shaped rotor or stator blade row mounted on an axisymmetric hub and confined in an axisymmetric duct, the computer program is also capable of analyzing flow fields about arbitrarily shaped wing body combinations, propellers, helicopter rotors in hover, and wind turbine rotors. The governing equation is solved numerically in a fully conservative form by using an artificial time concept, a finite volume technique, rotated type dependent differencing, successive line overrelaxation, and sequential boundary conforming grid refinement. An artificial viscosity is added in fully conservative form, and an initial guess for the potential field is applied, as determined by a two dimensional cascade analysis.

Dulikravich, D. S.↗

NASA NTRS1982

Combustion hot section technology

The Turbine Engine Hot Section Technology Combustion Program is briefly described. The overall objective of the project is to develop and verify improved and more accurate analysis methods for increasing the ability to design with confidence the combustion system for advanced aircraft turbine engines. The approach is to first assess and evaluate existing combustor aerothermal analysis models by means of a contracted effort initiated during FY-82. The program also includes both analytical and experimental research efforts in the areas of aerothermal modeling and liner cyclic life. It is expected that the combustor model development effort will generate improved understanding in the areas of high pressure flame radiation characteristics, model numerical methods and solution schemes, complex geometrical boundary conditions, fuel spray - flow field interactions, combustion kinetics, flow and mixing of dilution jets, turbulence and heat transfer, and soot and carbon formation.

Ercegovic, D. B.↗

NASA NTRS1980

Numerical calculation of transonic axial turbomachinery flows

This paper presents a numerical method and the results of a computer program for solving an exact, three-dimensional, full-potential equation that models rotating and nonrotating inviscid, absolutely irrotational, homentropic flows. Besides calculating the flows through an arbitrarily shaped rotor or stator blade row mounted on an axisymmetric hub and confined in an axisymmetric duct, the computer program is also capable of analysing flow fields about arbitrarily shaped wing-body combinations, propellers, helicopter rotors in hover, and wind turbine rotors. The governing equation is solved numerically in a fully conservative form by using an artificial time concept, a finite volume technique, rotated type-dependent differencing, successive line overrelaxation, and sequential boundary-conforming grid refinement. An artificial viscosity is added in fully conservative form; and an initial guess for the potential field is applied, as determined by a two-dimensional cascade analysis.

Dulikravich, D. S.↗

NASA NTRS1983

A new explicit method for the numerical solution of parabolic differential equations

A new method is derived for solving parabolic partial differential equations arising in transient heat conduction or in boundary-layer flows. The method is based on a combination of the modified differential quadrature (MDQ) method with the rational Runge-Kutta time-integration scheme. It is fully explicit, requires no matrix inversion, and is stable for any time-step for the heat equations. Burgers equation and the one- and two-dimensional heat equations are solved to demonstrate the accuracy and efficiency of the proposed algorithm. The present method is found to be very accurate and efficient when results are compared with analytic solutions.

Satofuka, N.↗

NASA NTRS1983

Wave scattering and guidance by dielectric waveguides with periodic surfaces

The extended-boundary-condition method is used to solve the wave scattering and guidance by a stratified medium with a periodic surface. First, the transition matrix for two sets of Flocquet waves incident from above and below a periodic surface separating the two media is derived. The results are applied to the diffraction of an incident plane wave by a multilayered medium having a periodic surface. The theoretical results can easily be applied to the guidance problem for which the complex guiding constant of the periodic structure is obtained. Numerical results for both the scattering and the guidance problems are obtained. The theory is compared with experiment for metallic gratings with surface-plasmon excitations. The calculated guiding constants of a thin film with a periodic sinusoidal surface are also compared with other numerical methods.

Chuang, S. L.↗

NASA NTRS1986

A computer code for three-dimensional incompressible flows using nonorthogonal body-fitted coordinate systems

In this report, a numerical method for solving the equations of motion of three-dimensional incompressible flows in nonorthogonal body-fitted coordinate (BFC) systems has been developed. The equations of motion are transformed to a generalized curvilinear coordinate system from which the transformed equations are discretized using finite difference approximations in the transformed domain. The hybrid scheme is used to approximate the convection terms in the governing equations. Solutions of the finite difference equations are obtained iteratively by using a pressure-velocity correction algorithm (SIMPLE-C). Numerical examples of two- and three-dimensional, laminar and turbulent flow problems are employed to evaluate the accuracy and efficiency of the present computer code. The user's guide and computer program listing of the present code are also included.

Chen, Y. S.↗

NASA NTRS1999

Comparison of AGE and Spectral Methods for the Simulation of Far-Wakes

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

Bisset, D. K.↗

NASA NTRS1986

Detailed near-wake flowfield surveys with comparison to an Euler method of an aspect ratio 4 rectangular wing

An experimental investigation of the flowfield in the near-wake of an aspect ratio 4 rectangular wing was conducted, providing a complete detailed set of data for use in the validation of computational methods. An angle of attack of 8 degrees and two Reynolds numbers 530,000 and 391,000 were investigated using pitot and six-hole probes. In addition, two types of flow visualization were employed. The data presented includes contours of total pressure, mean velocity, flow angularity, and vorticity distribution data at five chordwise stations of the near-wake ranging from 0.167 to 5.00 chord lengths aft of the trailing edge. The experimental results were compared to the predicted results of a 2-D Euler numerical method. The results predicted by an Euler method failed to accurately define the flowfield. Tangential velocities remained relatively constant over the range of X/C considered though increased in angle of attack and Reynolds number did bring about corresponding increases. Axial velocities also increased with angle of attack and Reynolds number but showed greater sensitivity to increases in X/C. Graphic displays and contours of the total pressure data indicate that roll-up of the wing tip vortex is essentially complete one and one half chords downstream of the trailing edge.

Klinge, M. D.↗

NASA NTRS1994

Large-eddy simulation of a boundary layer with concave streamwise curvature

Turbulence modeling continues to be one of the most difficult problems in fluid mechanics. Existing prediction methods are well developed for certain classes of simple equilibrium flows, but are still not entirely satisfactory for a large category of complex non-equilibrium flows found in engineering practice. Direct and large-eddy simulation (LES) approaches have long been believed to have great potential for the accurate prediction of difficult turbulent flows, but the associated computational cost has been prohibitive for practical problems. This remains true for direct simulation but is no longer clear for large-eddy simulation. Advances in computer hardware, numerical methods, and subgrid-scale modeling have made it possible to conduct LES for flows or practical interest at Reynolds numbers in the range of laboratory experiments. The objective of this work is to apply ES and the dynamic subgrid-scale model to the flow of a boundary layer over a concave surface.

Lund, Thomas S.↗

NASA NTRS1988

Three-dimensional atmospheric entry problem using method of matched asymptotic expansions

The analysis of a three-dimensional atmospheric entry problem using the method of matched asymptotic expansions is considered. A composite solution is formed in terms of an outer solution, an inner solution and a common solution. The outer solution is obtained from the gravitationally dominant region, whereas the aerodynamically dominant region contributes to the inner solution. The common solution accounts for the overlap between the outer and inner regions. In comparison to previous work, this simplified methodology yields explicit analytical expressions for various components of the composite solution without resorting to any type of transcendental equations to be solved only by numerical methods.

Naidu, D. S.↗

NASA NTRS1991

Simulations of curved turbulent boundary layers

The objective of this work is to develop a space-time accurate numerical method for the solution of incompressible Navier-Stokes equations in generalized coordinates. The resulting code is to be used for direct and large-eddy simulation of turbulence in complex geometries. In a previous paper, the system of Navier-Stokes equations in general curvilinear coordinates was solved by a second-order accurate finite-difference scheme. Satisfactory results were obtained for several flows in two and three dimensions. The system of Navier-Stokes for the fluxes are given in Orlandi (1989). The main deficiency of the numerical scheme was the large CPU time required for the solution of the Poisson equation for the 'pressure' field. The point SOR relaxation, in conjunction with a multigrid scheme, was used for the Poisson equation. In some cases, particularly with very fine grids, it was impossible to obtain a divergent-free flow. A preliminary attempt is made to compute the spatially evolving flow of Swearingen & Blackwelder. To reduce the streamwise distance, the inflow was at a distance x = 60 cm from the leading edge.

Orlandi, Paolo↗

NASA NTRSDate not supplied

Evolution of Astrodynamics for Multibody Environments, Numerical Precision and Dynamical Understanding

The innovative nature of dynamical systems theory and its applications in trajectory design enabled many recent missions by opening up new transfer and capture opportunities to meet mission requirements, reduce fuel costs, and address various launch constraints and contingencies. With a combination of proven numerical methods and an advanced understanding of the dynamics, an ‘experienced’ user can rapidly design multibody and cislunar missions such as Lunar IceCube and Artemis/Themis and relatively quickly conceive trajectory concepts for Sun-Earth mission orbits, such as the Space Weather Follow On-L1 Lissajous orbit—one that can be achieved and maintained.

Dynamical Systems↗

NASA NTRS1989

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

Unsteady rotor wake interactions with the empennage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies was developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interacting vortex wake is represented by an array of discrete vortices which, in turn, are represented by a finite-core model. The evolution of the interacting vortex wake is calculated by Lagrangian techniques. The viscous flow field of the two-dimensional body is calculated on an Eulerian grid. The flow around circular and elliptic cylinders in the absence of an interacting vortex wake was calculated. These results compare very well with other numerical results and with results obtained from experiment and thereby demonstrate the accuracy of the viscous solution. The interaction of a rotor wake with the flow about a 4 to 1 elliptic cylinder at 45 degree incidence was calculated for a Reynolds number of 3000. The results demonstrate the significant variations in the lift and drag on the elliptic cylinder in the presence of the interacting rotor wake.

Stremel, Paul M.↗