Engineering Papers⌕ Search

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.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Boundary-fitted coordinate systems for numerical solution of partial differential equations - A review

A comprehensive review of methods of numerically generating curvilinear coordinate systems with coordinate lines coincident with all boundary segments is given. Some general mathematical framework and error analysis common to such coordinate systems is also included. The general categories of generating systems are those based on conformal mapping, orthogonal systems, nearly orthogonal systems, systems produced as the solution of elliptic and hyperbolic partial differential equations, and systems generated algebraically by interpolation among the boundaries. Also covered are the control of coordinate line spacing by functions embedded in the partial differential operators of the generating system and by subsequent stretching transformation. Dynamically adaptive coordinate systems, coupled with the physical solution, and time-dependent systems that follow moving boundaries are treated. References reporting experience using such coordinate systems are reviewed as well as those covering the system development.

Thompson, J. F.↗
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 NTRS1987

Nonlinear beam kinematics for small strains and finite rotations

A simple matrix expression is obtained for the strain components of a beam in which the magnitudes of neither beam displacements nor rotations are explicitly restricted. The only kinematical restrictions are that: (1) strains are small compared to unity and (2) components of local rotation, a newly identified kinematical quantity, are of the order of the strains to a fractional power equal to at least one half. Local rotations are defined as the change of orientation of material elements off the beam reference axis relative to those on the beam reference axis. Local rotations appear explicitly in the resulting strain expressions, facilitating the treatment of both open- and closed-section beams in applications of the theory. The resulting strain components are expressed in a local Cartesian coordinate system and can be calculated directly in that way. Thus, one can use a curvilinear coordinate system that is natural to the beam problem without the complications that usually surround such an approach. Examples show the simplicity and the generality of the present approach as well as why previously published results differ among themselves concerning tension-torsion coupling.

Hodges, Dewey H.↗
NASA NTRS1988

Incompressible and anelastic flow simulations on numerically generated grids

In the numerical simulation of incompressible and anelastic flows, it is necessary to solve an elliptic equation at each time step. When the boundaries of such flows are nonrectangular, it may be advantageous to solve the equations on a new, numerically generated coordinate grid, in which the property of orthogonality has been preserved. Flow equations in general curvilinear coordinates maintaining the conservative form are given for both anelastic models using the momentum equations, and for incompressible models, using the vorticity equation. The general problem of grid-generation in two dimensions is presented, and a quasi-conformal transformation technique is discussed in detail. Some examples of grids generated by this technique are exhibited. Three examples of the flow of a stratified fluid over obstacles are presented, in which the grid-generation permits some new results to be obtained.

Sharman, R. D.↗
NASA NTRS1988

A formulation for the boundary-layer equations in general coordinates

This is a working paper in which a formulation is given for solving the boundary-layer equations in general body-fitted curvilinear coordinates while retaining the original Cartesian dependent variables. The solution procedure does not require that any of the coordinates be orthogonal, and much of the software developed for many Navier-Stokes schemes can be readily used. A limited number of calculations has been undertaken to validate the approach.

Steger, Joseph L.↗
NASA NTRS1982

Reynolds averaged Navier-Stokes computations of transonic flows The-state-of-the-art

The state of the art in computer simulations for transonic flowfields requiring solutions for the Navier-Stokes equations is assessed. It is noted that current simulations of transonic flowfields require comparisons with experimental results because the simulations are not free from discretization errors. Cases of turbulence are treated with weighted variables in a time-averaged scheme to yield Reynolds averaged Navier-Stokes equations. The turbulence is modeled in a first-order approach with a Reynolds stress tensor or a second-order approach where the tensor is obtained from the Navier-Stokes equations. The choice of coordinate systems and dependent variables for simulating the flow around airfoils is discussed, with particular attention to curvilinear coordinates. The determinations of boundary conditions is examined, along with numerical methods related to physical phenomena

Mehta, U.↗
NASA NTRS1984

Generation of surface-fitted coordinate grids

An accurate numerical solution of the Navier-Stokes equations is strongly dependent on the description of the flow geometry. In the past, rectangular grids led to an ill description of curved boundaries and consequently hindered finite-difference solutions. Now, it is accepted that body-fitted coordinate systems give a much better description of complex flow geometries. The generation of body-fitted curvilinear coordinate grids in two dimensions by the solution of a nonlinear system of elliptic equations has been proposed Thompson. It has been extended to three dimensions and applied in practical configurations.

Garon, A.↗
NASA NTRS1986

A numerical method for three-dimensional incompressible flows using nonorthogonal body-fitted coordinate systems

In the present paper, a numerical method for solving the equations of motion of three-dimensional incompressible flows in nonorthogonal body-fitted coordinate (BEC) systems has been developed and evaluated. 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 and a central differencing plus artificial dissipation scheme are used to approximate the convection terms in the governing equations. Effects of these two schemes on the accuracy of numerical predictions are studied. 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 numerical method.

Chen, Y. S.↗
NASA NTRS1991

Choice of velocity variables for complex flow computation

The issue of adopting the velocity components as dependent velocity variables for the Navier-Stokes flow computations is investigated. The viewpoint advocated is that a numerical algorithm should preferably honor both the physical conservation law in differential form and the geometric conservation law in discrete form. With the use of Cartesian velocity vector, the momentum equations in curvilinear coordinates can retain the full conservation-law form and satisfy the physical conservation laws. With the curvilinear velocity components, source terms appear in differential equations and hence the full conservation law form can not be retained. In discrete expressions, algorithms based on the Cartesian components can satisfy the geometric conservation-law form for convection terms but not for viscous terms; those based on the curvilinear components, on the other hand, cannot satisfy the geometric conservation-law form for either convection or viscous terms. Several flow solutions for domain with 90 and 360 degree turnings are presented to illustrate the issues of using the Cartesian velocity components and the staggered grid arrangement.

Shyy, W.↗
NASA NTRS1981

Mesh generation by conformal and quasiconformal mappings

It is pointed out that many recent advances in the finite-difference solution of elliptic equations have been limited to regions whose boundary contours coincide with coordinate lines of the Cartesian coordinate system. The reason for this is related to the fact that in the case of an arbitrary curvilinear coordinate system the original equation becomes much more complex. However, there is no added complexity if an orthogonal coordinate system is generated from a conformal mapping. In the present investigation, a finite difference method developed for the construction of conformal mappings has been generalized to construct quasi-conformal mappings. It is expected that the use of more sophisticated numerical algorithms could lead to improvements in both speed and accuracy. Quasi-conformal mappings have applications not only in the solution of elliptic equations but also in other areas such as orthogonal mesh generation on surfaces and the solution of certain fluid flow problems.

Mastin, C. W.↗
NASA NTRS1996

Coupled mixed-field laminate theory and finite element for smart piezoelectric composite shell structures

Mechanics for the analysis of laminated composite shells with piezoelectric actuators and sensors are presented. A new mixed-field laminate theory for piezoelectric shells is formulated in curvilinear coordinates which combines single-layer assumptions for the displacements and a layerwise representation for the electric potential. The resultant coupled governing equations for curvilinear piezoelectric laminates are described. Structural mechanics are subsequently developed and an 8-node finite-element is formulated for the static and dynamic analysis of adaptive composite structures of general laminations containing piezoelectric layers. Evaluations of the method and comparisons with reported results are presented for laminated piezoelectric-composite plates, a closed cylindrical shell with a continuous piezoceramic layer and a laminated composite semi-circular cantilever shell with discrete cylindrical piezoelectric actuators and/or sensors.

Saravanos, Dimitris A.↗
DOE OSTI2020

Non-Cartesian Coordinate Systems in the Portage Library

In order to be more broadly useful, the Portage framework needs to be able to support non-Cartesian coordinate systems. This document will lay out the most important details of implementing non-Cartesian coordinate systems in Portage. The most immediate need for this feature is integration with the EAP code base. Towards that end, the focus will be on certain curvilinear coordinate systems. However, this will provide a framework that can be used to implement other coordinate systems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Numerical solutions for viscous and potential flow about arbitrary two-dimensional bodies using body-fitted coordinate systems

A procedure for numerical solution of the time-dependent, incompressible Navier-Stokes equations for the flow about arbitrarily shaped two-dimensional bodies is given. This solution is based on a technique of automatic numerical generation of a curvilinear coordinate system having a coordinate line coincident with the body contour regardless of its shape. The implicit solution utilizes the vorticity-stream function formulation with a false-position iterative adjustment of the surface vorticity in satisfaction of the no-slip boundary condition. Excellent agreement with the Blasius boundary layer solution is obtained for a semi-infinite flat plate. Results are presented for Reynolds numbers up to 2000 for several airfoils and a cambered rock.

Thames, F. C.↗
NASA NTRS1983

High resolution applications of the Osher upwind scheme for the Euler equations

The 'Osher' scheme was introduced by Osher (1981). It represents an upwind finite-difference method for hyperbolic systems of conservation laws, including the Euler equations. In studies conducted by Osher (1981) and Osher and Solomon (1982), the method was applied to the nonisentropic form of the Euler equations in one dimension and the isentropically restricted form in two spatial dimensions, both in Cartesian coordinates. Chakravarthy and Osher (1982) have shown an approach for extending the Osher scheme to the Euler equations written for general geometries, taking into account the use of mappings to arbitrary curvilinear coordinate systems. The present investigation is concerned with the high resolution extension of the Osher scheme to second-order accuracy. Results are presented for several example problems, giving attention to quasi-one-dimensional Laval nozzle flow, a one-dimensional shock tube problem, and supersonic flow over a cylinder.

Chakravarthy, S. R.↗
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 NTRS1982

Elliptic grid generation

Various types of generating systems for boundary-conforming coordinate systems based on the numerical solution of systems of elliptic partial differential equations are discussed. Particular emphasis is given to the determination of functions in these equations which control the distribution of the curvilinear coordinate lines in the field.

Thompson, J. F.↗
NASA NTRS1980

An approximate factorization solution of the Navier-Stokes equations for transonic flow using body-fitted coordinates with application to NACA 64A010 airfoils

The implementation of the approximate factorization algorithm and its ability to efficiently and accurately describe transonic flow about an NACA 64A010 airfoil section is examined. The approximate factorization algorithm is developed from the nondimensional, conservative, vectorized Navier-Stokes equations expressed in curvilinear coordinates. Equations of state and transport coefficient relations appropriate to atmospheric air are appended to close the system of partial differential equations. An algebraic turbulence model is also incorporated into the equation set. This algorithm was verified by investigating the flow about an NACA 64A010 airfoil at 0, 2, and 3.5 deg angle of attack for free-stream conditions of 2,000,000 Reynolds number and 0.8 Mach number. Overall results were in good qualitative agreement with wind tunnel data sets. However, while nondimensional times of six were attained, numerical difficulties prevented any case from reaching a true steady state.

Copper, G. K.↗
NASA NTRS1978

Numerical solution of the Navier-Stokes equations for arbitrary blunt bodies in supersonic flows

A time-dependent, two-dimensional Navier-Stokes code employing the body-fitted coordinate technique has been developed for supersonic flows past blunt bodies of arbitrary shape. The computer program is based on the finite-difference approximation of the compressible Navier-Stokes equations transformed to nonorthogonal curvilinear coordinates with the contravariant components of the velocity vector as dependent variables. The bow shock ahead of the body is obtained as part of the solution, by 'shock capturing'. Numerical solutions of the complete equations are presented in detail for free-stream Mach number 4.6, Reynolds number 10,000, and an isothermal wall temperature of 556 K for a circular cylinder with the free-stream outer boundaries forming a hyperbola in the front and a circular arc in the back.

Warsi, Z. U. A.↗