SEARCH · Engineering Papers

Results for “curvilinear coordinates”

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 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1978

Forebody and afterbody solutions of the Navier-Stokes equations for supersonic flow over blunt bodies in a generalized orthogonal coordinate system

A coordinate transformation, which can approximate many different two-dimensional and axisymmetric body shapes with an analytic function, is used as a basis for solving the Navier-Stokes equations for the purpose of predicting 0 deg angle of attack supersonic flow fields. The transformation defines a curvilinear, orthogonal coordinate system in which coordinate lines are perpendicular to the body and the body is defined by one coordinate line. This system is mapped in to a rectangular computational domain in which the governing flow field equations are solved numerically. Advantages of this technique are that the specification of boundary conditions are simplified and, most importantly, the entire flow field can be obtained, including flow in the wake. Good agreement has been obtained with experimental data for pressure distributions, density distributions, and heat transfer over spheres and cylinders in supersonic flow. Approximations to the Viking aeroshell and to a candidate Jupiter probe are presented and flow fields over these shapes are calculated.

Gnoffo, P. A.↗

DOE OSTI2020

Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

In this work, we develop a numerical method for solving the acoustic wave equation in covariant form on staggered curvilinear grids in an energy conserving manner. The use of a covariant basis decomposition leads to a rotationally invariant scheme that outperforms a Cartesian basis decomposition on rotated grids. The discretization is based on high order Summation-By-Parts (SBP) operators and preserves both symmetry and positive definiteness of the contravariant metric tensor. To improve accuracy and decrease computational cost, we also derive a modified discretization of the metric tensor that leads to a conditionally stable discretization. Bounds are derived that yield a point-wise condition that can be evaluated to check for stability of the modified discretization. This condition shows that the interpolation operators should be constructed such that their norm is close to one.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1972

Large amplitude flexural vibration of thin elastic flat plates and shells

The general equations governing the large amplitude flexural vibration of any thin elastic shell using curvilinear orthogonal coordinates are derived and consist of two coupled, nonlinear, partial differential equations in the normal displacement w and the stress function F. From these equations, the governing equations for the case of shells of revolution or flat plates can be readily obtained as special cases. The material of the shell or plate is isotropic and homogeneous and Hooke's law for the two-dimensional case is valid. It is suggested that the difference between the hardening type of nonlinearity in the case of flat plates and straight beams and the softening type of nonlinearity in the case of shells and rings can, in general, be traced to the amount of curvature present in the underformed median surface of the structure concerned.

Pandalia, K. A. V.↗

NASA NTRS1980

Nonlinear finite element analysis: An alternative formulation

A geometrical nonlinear analysis based on an alternative definition of strain is presented. Expressions for strain are obtained by computing the change in length of the base vectors in the curvilinear element coordinate system. The isoparametric element formulation is assumed in the global Cartesian coordinate system. The approach is based on the minimization of the strain energy, and the resulting nonlinear equations are solved by the modified Newton method. Integration of the first and second variation of the strain energy is performed numerically in the case of two and three dimensional elements. Application is made to a simple long cantilever beam.

Merazzi, S.↗

NASA NTRS1983

Multi-grid solution of Neumann pressure problem for viscous flows using primitive variables

The multi-grid (MG) technique has been advanced for use with the Neumann boundary-value problem in clustered curvilinear orthogonal coordinates. This comprises an important step in the analysis of viscous flows using the velocity-pressure formulation of the Navier-Stokes equations. With successive over-relaxation (SOR) as the smoothing operator and with suitably formulated restriction and coarse-grid-correction operators, a 4-grid procedure enhances the efficiency of fine-grid solutions of the Neumann problem by a factor of 3 to 5, depending on the problem parameters. Thy influence of the smoothing operator is also examined by employing the alternating-direction implicit and the strongly implicit techniques instead of SOR.

Ghia, U.↗

NASA NTRS1983

A direct method for the solution of unsteady two-dimensional incompressible Navier-Stokes equations

The unsteady incompressible Navier-Stokes equations are formulated in terms of vorticity and stream function in generalized curvilinear orthogonal coordinates to facilitiate analysis of flow configurations with general geometries. The numerical method developed solves the conservative form of the transport equation using the alternating-direction implicit method, whereas the stream-function equation is solved by direct block Gaussian elimination. The method is applied to a model problem of flow over a back-step in a doubly infinite channel, using clustered conformal coordinates. One-dimensional stretching functions, dependent on the Reynolds number and the asymptotic behavior of the flow, are used to provide suitable grid distribution in the separation and reattachment regions, as well as in the inflow and outflow regions. The optimum grid distribution selected attempts to honor the multiple length scales of the separated-flow model problem. The asymptotic behavior of the finite-differenced transport equation near infinity is examined and the numerical method is carefully developed so as to lead to spatially second-order accurate wiggle-free solutions, i.e., with minimum dispersive error. Results have been obtained in the entire laminar range for the backstep channel and are in good agreement with the available experimental data for this flow problem.

Ghia, K. N.↗

NASA NTRS1985

A 3-dimensional Navier-Stokes-Euler code for blunt-body flow computations

The shock-layer flowfield is obtained with or without viscous and heat-conducting dissipations from the conservative laws of fluid dynamics equations using a shock-fitting implicity finite-difference technique. The governing equations are cast in curvilinear-orthogonal coordinates and transformed to the domain between the shock and the body. Another set of equations is used for the singular coordinate axis, which, together with a cone generator away from the stagnation point, encloses the computation domain. A time-dependent alternating direction implicit factorization technique is applied to integrate the equations with local-time increment until a steady solution is reached. The shock location is updated after the flowfield computation, but the wall conditions are implemented into the implicit procedure. Innovative procedures are introduced to define the initial flowfield, to treat both perfect and equilibrium gases, to advance the solution on a coarse-to-fine grid sequence, and to start viscous flow computations from their corresponding inviscid solutions. The results are obtained from a grid no greater than 28 by 18 by 7 and converged within 300 integration steps. They are of sufficient accuracy to start parabolized Navier-Stokes or Euler calculations beyond the nose region, to compare with flight and wind-tunnel data, and to evaluate conceptual designs of reentry spacecraft.

Li, C. P.↗

NASA NTRS1984

Multi-grid technique for the solution of incompressible Navier-Stokes equations

A computationally efficient method is proposed for obtaining fine-mesh solutions of the Navier-Stokes equations for high-Re flow, including separation. The method involves implementing a multi-grid solution procedure with suitably chosen elements, such that the resulting overall solution technique is efficient as well as robust. The robustness and efficiency of the solution technique are demonstrated by applying it to three model problems: flow in a driven cavity, downstream asymptotic flow in curved ducts of square and polar sections, and Newmann boundary-value problem in clustered curvilinear orthogonal coordinates.

Ghia, K. N.↗

NASA NTRS1985

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗

NASA NTRS1990

A three-dimensional, compressible, laminar boundary-layer method for general fuselages. Volume 1: Numerical method

A procedure for calculating 3-D, compressible laminar boundary layer flow on general fuselage shapes is described. The boundary layer solutions can be obtained in either nonorthogonal 'body oriented' coordinates or orthogonal streamline coordinates. The numerical procedure is 'second order' accurate, efficient and independent of the cross flow velocity direction. Numerical results are presented for several test cases, including a sharp cone, an ellipsoid of revolution, and a general aircraft fuselage at angle of attack. Comparisons are made between numerical results obtained using nonorthogonal curvilinear 'body oriented' coordinates and streamline coordinates.

Wie, Yong-Sun↗

NASA NTRS1993

Stability and transition on swept wings

This paper describes the extension and application of the Parabolized Stability Equations (PSE) to the stability and transition of the supersonic three-dimensional laminar boundary layer on a swept wing. The problem formulation uses a general coordinate transformation for arbitrary curvilinear body-fitted computational grids. Some testing using these coordinates is briefly described to help validate the software used for the investigation. The disturbance amplitude ratios as a function of chord position for supersonic (Mach 1.5) boundary layers on untapered, untwisted wings of different sweep angles are then presented and compared with those obtained from local parallel analyses.

Stuckert, Greg↗

NASA NTRS1982

Stability of Goertler vortices in boundary layers

A formal analysis of Goertler-type instability is presented. The boundary-layer and disturbance equations are formulated in a general, orthogonal, curvilinear system of coordinates constructed from the inviscid flow over a curved surface. Effects of curvature on the boundary-layer flow are analyzed. The basic approximation for the disturbance equations is presented and solved numerically. Previous analyses are discussed and compared with our analysis. It is shown that the general system of coordinates developed in this analysis and the correct order-of-magnitude analysis of the disturbance velocities with two velocity scales leads to a rational foundation for future work in Goertler vortices.

Floryan, J. M.↗

NASA NTRS1984

Goertler instability of compressible boundary layers

The instability of the laminar compressible boundary-layer flows along concave surfaces is investigated. The linearized disturbance equations for the three-dimensional, counter-rotating, longitudinal-type vortices in two-dimensional boundary layers are presented in an orthogonal curvilinear system of coordinates. The basic approximation of the disturbance equations, which includes the effect of the growth of the boundary layer, is considered and solved numerically.

El-Hady, N. M.↗

NASA NTRS1990

Numerical solution of the boundary-layer equations for a general aviation fuselage

Numerical solutions of the three-dimensional, compressible laminar boundary-layer equations for a general aviation fuselage are presented. The numerical procedure is second-order accurate and independent of the cross-flow velocity direction. Numerical results are presented for a Mach number and unit Reynolds number of 0.3 and 7 x 10 to the 6th/m, respectively, for angles of attack of 0 and 3 deg. Comparisons are made between results obtained using a general nonorthogonal curvilinear body oriented coordinate system and a streamline coordinate system. Axisymmetric analogue results are also compared with the three-dimensional solutions.

Wie, Yong-Sun↗

NASA NTRS1986

An implicit flux-difference splitting scheme for three-dimensional, incompressible Navier-Stokes solutions to leading edge vortex flows

A new, implicit finite-difference scheme designed to solve the conservative, flux-difference split Navier-Stokes equations is used to compute incompressible vortex flows around delta wings. The completely vectorizable hybrid algorithm is constructed in delta form for steady state solutions independent of the time-step sizes. The scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. The governing equations are solved in curvilinear, body-fitted coordinates for treating complex geometries. The computed flow field results are compared with other theoretical and experimental data.

Hartwich, P.-M.↗

NASA NTRS1981

Principal stress formulation of continua motion

Principal stress coordinates are formulated for a material without tension, with the equations of motion projected onto a selected coordinate system. The translations are achieved by forming the local basis vector in correspondence with the local stress tensor. The principal stress coordinates are described as orthogonal and curvilinear, and the metric is defined. The kinematics and dynamics of the principal stress coordinates are analyzed, the latter in orthogonal coordinates, and an analytic formulation of the dynamics of a no-tension material is provided. Finally, elastic granular material such as sand, soil, etc., which are more complicated versions of no-tension material, are explored.

Zak, M.↗

NASA NTRS1979

Logarithmic spiral grids for image processing

A picture digitization grid based on logarithmic spirals rather than Cartesian coordinates is presented. Expressing this curvilinear grid as a conformal exponential mapping reveals useful image processing properties. The mapping induces a computational simplification that suggests parallel architectures in which most geometric transformations are effected by data shifting in memory rather than arithmetic on coordinates. These include fast, parallel noise-free rotation, scaling, and some projective transformations of pixel defined images. Conformality of the mapping preserves local picture-processing operations such as edge detection.

Weiman, C. F. R.↗

NASA NTRS1979

Numerical analysis of turbulent wakes of axial flow turbomachinery rotor blades

The wakes of rotor blades are turbulent, three-dimensional, and are subjected to curvature and rotation effects. The objective of this study is to predict the development of such wakes and compare the predictions with the existing data. A finite difference procedure is employed in the numerical analysis of the wake utilizing the continuity, momentum, and turbulence closure equations in the rotating curvilinear and non-orthogonal coordinate system. The turbulence closure is affected by the two equations model and the simplified Reynolds stress model. The predictions from both the models agree well with the mean velocity and the turbulence data.

Hah, C.↗