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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Multigrid calculation of internal flows in complex geometries

The development, validation, and application of a general purpose multigrid solution algorithm and computer program for the computation of elliptic flows in complex geometries is presented. This computer program combines several desirable features including a curvilinear coordinate system, collocated arrangement of the variables, and Full Multi-Grid/Full Approximation Scheme (FMG/FAS). Provisions are made for the inclusion of embedded obstacles and baffles inside the flow domain. The momentum and continuity equations are solved in a decoupled manner and a pressure corrective equation is used to update the pressures such that the fluxes at the cell faces satisfy local mass continuity. Despite the computational overhead required in the restriction and prolongation phases of the multigrid cycling, the superior convergence results in reduced overall CPU time. The numerical scheme and selected results of several validation flows are presented. Finally, the procedure is applied to study the flowfield in a side-inlet dump combustor and twin jet impingement from a simulated aircraft fuselage.

Smith, K. M.↗
NASA NTRS1992

Finite-difference algorithms for the time-domain Maxwell's equations - A numerical approach to RCS analysis

The applications of two CFD-based finite-difference methods to computational electromagnetics are investigated. In the first method, the time-domain Maxwell's equations are solved using the explicit Lax-Wendroff scheme and in the second method, the second-order wave equations satisfying the Maxwell's equations are solved using the implicit Crank-Nicolson scheme. The governing equations are transformed to a generalized curvilinear coordinate system and solved on a body-conforming mesh using the scattered-field formulation. The induced surface current and the bistatic radar cross section are computed and the results are validated for several two-dimensional test cases involving perfectly-conducting scatterers submerged in transverse-magnetic plane waves.

Vinh, Hoang↗
NASA NTRS1992

Stability investigations of airfoil flow by global analysis

As the result of global, non-parallel flow stability analysis the single value of the disturbance growth-rate and respective frequency is obtained. This complex value characterizes the stability of the whole flow configuration and is not referred to any particular flow pattern. The global analysis assures that all the flow elements (wake, boundary and shear layer) are taken into account. The physical phenomena connected with the wake instability are properly reproduced by the global analysis. This enhances the investigations of instability of any 2-D flows, including ones in which the boundary layer instability effects are known to be of dominating importance. Assuming fully 2-D disturbance form, the global linear stability problem is formulated. The system of partial differential equations is solved for the eigenvalues and eigenvectors. The equations, written in the pure stream function formulation, are discretized via FDM using a curvilinear coordinate system. The complex eigenvalues and corresponding eigenvectors are evaluated by an iterative method. The investigations performed for various Reynolds numbers emphasize that the wake instability develops into the Karman vortex street. This phenomenon is shown to be connected with the first mode obtained from the non-parallel flow stability analysis. The higher modes are reflecting different physical phenomena as for example Tollmien-Schlichting waves, originating in the boundary layer and having the tendency to emerge as instabilities for the growing Reynolds number. The investigations are carried out for a circular cylinder, oblong ellipsis and airfoil. It is shown that the onset of the wake instability, the waves in the boundary layer, the shear layer instability are different solutions of the same eigenvalue problem, formulated using the non-parallel theory. The analysis offers large potential possibilities as the generalization of methods used till now for the stability analysis.

Morzynski, Marek↗
NASA NTRS1993

A numerical study of confined turbulent jets

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Zhu, J.↗
NASA NTRS1992

Some Aspects of Bifurcation Structure of Laminar Flow in Curved Ducts

A bifurcation study is made of laminar flow in curved ducts. The problem is formulated in a curvilinear coordinate system, and the governing equations, after orthogonal mapping is applied, are solved numerically by an iterative finite-difference method. Many computer runs were made with various duct cross-sections ranging from a circle to a square, to learn the transition of bifurcation structure with this change in cross-section and to reconcile the differences between them. In addition, a simpler technique is proposed to generate symmetric four-cell solutions in a circular pipe and a means is put forward to stabilize four-vortex structures in a complete cross-section.

Kao, Hsiao C.↗
NASA NTRS1993

Numerical investigation of receptivity to freestream vorticity

The receptivity to freestream vorticity of the boundary layer over a flat plate with an elliptic leading edge is studied numerically. A finite-difference scheme that is second-order accurate in both time and space is utilized to solve the incompressible Navier-Stokes system in general curvilinear coordinates with the vorticity and stream function as dependent variables. Simulations show, for the conditions considered, a linear variation in the Tollmien-Schlichting response with forcing amplitude for perturbations of freestream velocity that are either symmetrical or asymmetrical with respect to the basic state stagnation streamline.

Buter, Thomas A.↗
NASA NTRS1993

On the generation of secondary motion in circular to rectangular transition ducts

A computer code has been developed for the calculation of three-dimensional turbulent flow in circular to rectangular transition ducts. The governing equations are written in a general curvilinear coordinate system so that all singly-connected geometrical configurations can be accommodated. Turbulence closure is through the standard k-epsilon model or a Reynolds stress closure model. Computed results show qualitative agreement with experimental data. Pressure-induced secondary flow is the dominant mode through the transition, but turbulence-induced secondary flow becomes prominent some distance downstream of the transition.

Demuren, A. O.↗
NASA NTRS1993

Waverider design for generalized shock geometries

A new design method for waverider configurations with generalized shock geometries is presented. A cross-stream Euler marching procedure is used to define the post-shock flow field, and an arbitrary 3D shock shape is specified as an input. This approach makes it possible to use nonaxisymmetric shock topologies with nonconstant shock strengths. The problem is reformulated in the proper curvilinear coordinate system to suppress its fundamental ill-posedness.

Jones, K. D.↗
NASA NTRS1992

A method of smooth bivariate interpolation for data given on a generalized curvilinear grid

A method of locally bicubic interpolation is presented for data given at the nodes of a two-dimensional generalized curvilinear grid. The physical domain is transformed to a computational domain in which the grid is uniform and rectangular by a generalized curvilinear coordinate transformation. The metrics of the transformation are obtained by finite differences in the computational domain. Metric derivatives are determined by repeated application of the chain rule for partial differentiation. Given the metrics and the metric derivatives, the partial derivatives required to determine a locally bicubic interpolant can be estimated at each data point using finite differences in the computational domain. A bilinear transformation is used to analytically transform the individual quadrilateral cells in the physical domain into unit squares, thus allowing the use of simple formulas for bicubic interpolation.

Zingg, David W.↗
NASA NTRS1993

Multigrid methods and high order finite difference for flow in transition - Effects of isolated and distributed roughness elements

The high order finite difference and multigrid methods have been successfully applied to direct numerical simulation (DNS) for flow transition in 3D channels and 3D boundary layers with 2D and 3D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semicoarsening multigrid method associated with line distributive relaxation scheme, and a new treatment of the outflow boundary condition, which needs only a very short buffer domain to damp all wave reflection, are developed. These approaches make the multigrid DNS code very accurate and efficient. This makes us not only able to do spatial DNS for the 3D channel and flat plate at low computational costs, but also able to do spatial DNS for transition in the 3D boundary layer with 3D single and multiple roughness elements. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments.

Liu, C.↗
NASA NTRS1993

Effect of curvature on stationary crossflow instability of a three-dimensional boundary layer

An incompressible three-dimensional laminar boundary-layer flow over a swept wing is used as a model to study both the wall-curvature and streamline-curvature effects on the stationary crossflow instability. The basic state is obtained by solving the full Navier-Stokes (N-S) equations numerically. The linear disturbance equations are cast on a fixed, body-intrinsic, curvilinear coordinate system. Those nonparallel terms which contribute mainly to the streamline-curvature effect are retained in the formulation of the disturbance equations and approximated by their local finite difference values. The resulting eigenvalue problem is solved by a Chebyshev collocation method. The present results indicate that the convex wall curvature has a stabilizing effect, whereas the streamline curvature has a destabilizing effect. A validation of these effects with an N-S solution for the linear disturbance flow is provided.

Lin, Ray-Sing↗
NASA NTRS1993

Analysis of direct numerical simulation data of a Mach 4.5 transitional boundary-layer flow

This paper describes the creation, by temporal direct numerical simulation and the analysis based on the Reynolds stress transport equations, of a high-quality data set that represents the laminar-turbulent transition of a high-speed boundary-layer flow. Following Pruett and Zang (1992), and with the help of algorithmic refinements, the evolution of an axial, Mach 4.5 boundary-layer flow along a hollow cylinder is simulated numerically. Favre-averaged Reynolds stress transport equations are derived in generalized curvilinear coordinates and are then specialized to the cylindrical geometry at hand. Reynolds stresses and various turbulence quantities, such as turbulent kinetic energy and turbulent Mach number, are calculated from the numerical data at various stages of the transition process. The kinetic energy 'budgets' are constructed from the transport equations. Various contributing terms for the evolution of kinetic energy, like the rates of production and dissipation, transport, and diffusion, are presented. The compressible dissipation rate is small in comparison with the solenoidal dissipation rate for all times. The pressure-dilatation term is of the same order of magnitude as the compressible dissipation rate.

Dinavahi, Surya P. G.↗
NASA NTRS1952

Generalization of Boundary-Layer Momentum-Integral Equations to Three-Dimensional Flows Including Those of Rotating System

The Navier-Stokes equations of motion and the equation of continuity are transformed so as to apply to an orthogonal curvilinear coordinate system rotating with a uniform angular velocity about an arbitrary axis in space. A usual simplification of these equations as consistent with the accepted boundary-layer theory and an integration of these equations through the boundary layer result in boundary-layer momentum-integral equations for three-dimensional flows that are applicable to either rotating or nonrotating fluid boundaries. These equations are simplified and an approximate solution in closed integral form is obtained for a generalized boundary-layer momentum-loss thickness and flow deflection at the wall in the turbulent case. A numerical evaluation of this solution carried out for data obtained in a curving nonrotating duct shows a fair quantitative agreement with the measures values. The form in which the equations are presented is readily adaptable to cases of steady, three-dimensional, incompressible boundary-layer flow like that over curved ducts or yawed wings; and it also may be used to describe the boundary-layer flow over various rotating surfaces, thus applying to turbomachinery, propellers, and helicopter blades.

Mager, Arthur↗
NASA NTRS1943

On a new method for calculating the potential flow past a body of revolution

A new method is presented for obtaining the velocity potential of the flow about a body of revolution moving uniformly in the direction of its axis of symmetry in a fluid otherwise at rest. This method is based essentially on the fact that the form of the differential equation for the velocity potential is invariant with regard to conformal transformation of the meridian plane. By means of the conformal transformation of the meridian profile into a circle a system of orthogonal curvilinear coordinates is obtained, the main feature of which is that one of the coordinate lines is the meridian profile itself. The use of this type of coordinate system yields a simple expression of the boundary condition at the surface of the solid and leads to a rational process of iteration for the solution of the differential equation for the velocity potential. It is shown that the velocity potential for an arbitrary body of revolution may be expressed in terms of universal functions which, although not normal, are obtainable by means of simple quadratures. The general results are applied to a body of revolution obtained by revolving a symmetrical Joukowski profile about its axis of symmetry. A numerical example further serves to illustrate the theory.

Kaplan, Carl↗
NASA NTRS1993

Multigrid methods for flow transition in three-dimensional boundary layers with surface roughness

The efficient multilevel adaptive method has been successfully applied to perform direct numerical simulations (DNS) of flow transition in 3-D channels and 3-D boundary layers with 2-D and 3-D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semi-coarsening multigrid method associated with line distributive relaxation scheme, and an improved outflow boundary-condition treatment, which needs only a very short buffer domain to damp all order-one wave reflections, are developed. These approaches make the multigrid DNS code very accurate and efficient. This allows us not only to be able to do spatial DNS for the 3-D channel and flat plate at low computational costs, but also to do spatial DNS for transition in the 3-D boundary layer with 3-D single and multiple roughness elements, which would have extremely high computational costs with conventional methods. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments. The contribution of isolated and distributed roughness to transition is analyzed.

Liu, Chaoqun↗
NASA NTRS1993

Multigrid properties of upwind-biased data reconstructions

The multigrid properties of two data reconstruction methods used for achieving second-order spatial accuracy when solving the two-dimensional Euler equations are examined. The data reconstruction methods are used with an implicit upwind algorithm which uses linearized backward-Euler time-differencing. The solution of the resulting linear system is performed by an iterative procedure. In the present study only regular quadrilateral grids are considered, so a red-black Gauss-Seidel iteration is used. Although the Jacobian is approximated by first-order upwind extrapolation, two alternative data reconstruction techniques for the flux integral that yield higher-order spatial accuracy at steady state are examined. The first method, probably most popular for structured quadrilateral grids, is based on estimating the cell gradients using one-dimensional reconstruction along curvilinear coordinates. The second method is based on Green's theorem. Analysis and numerical results for the two dimensional Euler equations show that data reconstruction based on Green's theorem has superior multigrid properties as compared to the one-dimensional data reconstruction method.

Warren, Gary P.↗
NASA NTRS1992

Streamline-curvature effect in three-dimensional boundary layers

The effect of including wall and streamline curvature terms in swept-wing boundary-layer stability calculations is studied. The linear disturbance equations are cast on a fixed, body-intrinsic, curvilinear coordinate system. Those nonparallel terms which contribute mainly to the streamline-curvature effect are retained in this formulation and approximated by their local finite-difference values. Convex-wall curvature has a stabilizing effect, while streamline curvature is destabilizing if the curvature exceeds a critical value.

Reed, Helen L.↗
NASA NTRS1994

Discrete-Layer Piezoelectric Plate and Shell Models for Active Tip-Clearance Control

The objectives of this work were to develop computational tools for the analysis of active-sensory composite structures with added or embedded piezoelectric layers. The targeted application for this class of smart composite laminates and the analytical development is the accomplishment of active tip-clearance control in turbomachinery components. Two distinct theories and analytical models were developed and explored under this contract: (1) a discrete-layer plate theory and corresponding computational models, and (2) a three dimensional general discrete-layer element generated in curvilinear coordinates for modeling laminated composite piezoelectric shells. Both models were developed from the complete electromechanical constitutive relations of piezoelectric materials, and incorporate both displacements and potentials as state variables. This report describes the development and results of these models. The discrete-layer theories imply that the displacement field and electrostatic potential through-the-thickness of the laminate are described over an individual layer rather than as a smeared function over the thickness of the entire plate or shell thickness. This is especially crucial for composites with embedded piezoelectric layers, as the actuating and sensing elements within these layers are poorly represented by effective or smeared properties. Linear Lagrange interpolation polynomials were used to describe the through-thickness laminate behavior. Both analytic and finite element approximations were used in the plane or surface of the structure. In this context, theoretical developments are presented for the discrete-layer plate theory, the discrete-layer shell theory, and the formulation of an exact solution for simply-supported piezoelectric plates. Finally, evaluations and results from a number of separate examples are presented for the static and dynamic analysis of the plate geometry. Comparisons between the different approaches are provided when possible, and initial conclusions regarding the accuracy and limitations of these models are given.

Heyliger, P. R.↗