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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1988

A transonic interactive boundary-layer theory for laminar and turbulent flow over swept wings

A 3-D laminar and turbulent boundary-layer method is developed for compressible flow over swept wings. The governing equations and curvature terms are derived in detail for a nonorthogonal, curvilinear coordinate system. Reynolds shear-stress terms are modeled by the Cebeci-Smith eddy-viscosity formulation. The governing equations are descretized using the second-order accurate, predictor-corrector finite-difference technique of Matsuno, which has the advantage that the crossflow difference formulas are formed independent of the sign of the crossflow velocity component. The method is coupled with a full potential wing/body inviscid code (FLO-30) and the inviscid-viscous interaction is performed by updating the original wing surface with the viscous displacement surface calculated by the boundary-layer code. The number of these global iterations ranged from five to twelve depending on Mach number, sweep angle, and angle of attack. Several test cases are computed by this method and the results are compared with another inviscid-viscous interaction method (TAWFIVE) and with experimental data.

Woodson, Shawn H.↗

NASA NTRS1989

Development of direct-inverse 3-D methods for applied transonic aerodynamic wing design and analysis

Progress in the direct-inverse wing design method in curvilinear coordinates has been made. This includes the remedying of a spanwise oscillation problem and the assessment of grid skewness, viscous interaction, and the initial airfoil section on the final design. It was found that, in response to the spanwise oscillation problem that designing at every other spanwise station produced the best results for the cases presented, a smoothly varying grid is especially needed for the accurate design at the wing tip, the boundary layer displacement thicknesses must be included in a successful wing design, the design of high and medium aspect ratio wings is possible with this code, and the final airfoil section designed is fairly independent of the initial section.

Carlson, Leland A.↗

NASA NTRS1989

A direct procedure for interpolation on a structured curvilinear two-dimensional grid

A direct procedure is presented for locally bicubic interpolation on a structured, curvilinear, two-dimensional grid. The physical (Cartesian) space is transformed to a computational space in which the grid is uniform and rectangular by a generalized curvilinear coordinate transformation. Required partial derivative information is obtained by finite differences in the computational space. The partial derivatives in physical space are determined by repeated application of the chain rule for partial differentiation. A bilinear transformation is used to analytically transform the individual quadrilateral cells in physical space into unit squares. The interpolation is performed within each unit square using a piecewise bicubic spline.

Zingg, David W.↗

NASA NTRS1989

An interactive three-dimensional boundary-layer method for transonic flow over swept wings

A three-dimensional laminar/turbulent boundary-layer method is developed for transonic flow over swept wings. The governing equations and curvature terms are written for a nonorthogonal body-oriented curvilinear coordinate system. The viscous method is coupled to the full-potential inviscid code, FLO-30, through a displacement-surface interaction. Typically, for transonic Mach numbers and moderate sweep angles, between 5 and 12 viscous updates are required for convergence. The method is applied to a variety of wing planforms, and the results are also compared to those obtained with the three-dimensional integral inviscid-viscous interaction method TAWFIVE (Streett and Melson, 1983).

Woodson, Shawn H.↗

NASA NTRS1989

Numerical solution of the incompressible Navier-Stokes equations for steady-state and time-dependent problems

An algorithm for the solution of the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. The algorithm is based on the method of artificial compressibility and uses a higher-order flux-difference splitting technique for the convective terms and a second-order central difference for the viscous terms. The steady-state solution of flow through a square duct with a 90 deg bend is computed and the results are compared with experimental data. Good agreement is observed. A comparison with an analytically known exact solution is then performed to verify the time accuracy of the algorithm. Finally, the flow through an artificial heart configuration with moving boundaries is calculated and presented.

Rogers, Stuart E.↗

NASA NTRS1989

Numerical simulation of axisymmetric turbulent flow in combustors and diffusers

Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Yung, Chain-Nan↗

NASA NTRS1989

An analysis of curvature effects for the control of wall-bounded shear flows

The Reynolds stress transport equations are used to predict the effects of simultaneous and sequential combinations of distortions on turbulent boundary layers. The equations are written in general orthogonal curvilinear coordinates, with the curvature terms expressed in terms of the principal radii of curvature of the respective coordinate surfaces. Results are obtained for the cases of two-dimensional and three-dimensional flows in the limit where production and pressure-strain redistribution dominate over diffusion effects.

Gatski, T. B.↗

NASA NTRS1988

A study of turbulent flow computations in an angled duct with a step

This paper presents a method for effective computation of turbulent flows in an angled duct which has a step at the angled section by using a grid generation technique. An effect of flow curvature was analyzed on the separating and reattaching flows. The computations are based on the Reynolds averaged transport equations which are transformed for use on a generalized curvilinear coordinate system. A special care was taken on the treatment of pressure corrections. Wall function boundary conditions for k-epsilon model were developed.

Maruszewski, J. P.↗

NASA NTRS1989

Flow of rarefied gases over two-dimensional bodies

A kinetic-theory analysis is made of the flow of rarefied gases over two-dimensional bodies of arbitrary curvature. The Boltzmann equation simplified by a model collision integral is written in an arbitrary orthogonal curvilinear coordinate system, and solved by means of finite-difference approximation with the discrete ordinate method. A numerical code is developed which can be applied to any two-dimensional submerged body of arbitrary curvature for the flow regimes from free-molecular to slip at transonic Mach numbers. Predictions are made for the case of a right circular cylinder.

Jeng, Duen-Ren↗

NASA NTRS1989

A time accurate finite volume high resolution scheme for three dimensional Navier-Stokes equations

A time accurate, three-dimensional, finite volume, high resolution scheme for solving the compressible full Navier-Stokes equations is presented. The present derivation is based on the upwind split formulas, specifically with the application of Roe's (1981) flux difference splitting. A high-order accurate (up to the third order) upwind interpolation formula for the inviscid terms is derived to account for nonuniform meshes. For the viscous terms, discretizations consistent with the finite volume concept are described. A variant of second-order time accurate method is proposed that utilizes identical procedures in both the predictor and corrector steps. Avoiding the definition of midpoint gives a consistent and easy procedure, in the framework of finite volume discretization, for treating viscous transport terms in the curvilinear coordinates. For the boundary cells, a new treatment is introduced that not only avoids the use of 'ghost cells' and the associated problems, but also satisfies the tangency conditions exactly and allows easy definition of viscous transport terms at the first interface next to the boundary cells. Numerical tests of steady and unsteady high speed flows show that the present scheme gives accurate solutions.

Liou, Meng-Sing↗

NASA NTRS1989

Curvature effects on the stability of three-dimensional laminar boundary layers

The linear stability equations for compressible, three-dimensional laminar boundary layer flow are derived in an orthogonal curvilinear coordinate system. The system of equations is solved using a finite difference scheme in order to study the effects of streamline and surface curvature and compressibility on the stability of the flow past a swept wing. It is known that convex surface curvature can have a stabilizing effect on the laminar boundary layer. Conversely, concave surface curvature can be destabilizing. The magnitude of these effects for swept wings is determined. Results indicate that amplification rates and hence, N-factors, for the flow over the convex upper surface of a swept wing can be reduced by about 15 to 45 percent when curvature effects are included in the linear stability analysis. The results of the calculations show that concave curvature destabilizes crossflow type disturbances with a significant increase in amplification rate. In addition, comparisons are made with some experimental results on a swept concave-convex surface. Calculated velocity vector plots show good agreement with observed disturbances in the laminar boundary layer over the concave surface.

Collier, F. S., Jr.↗

NASA NTRS1989

Compact finite volume methods for the diffusion equation

An approach to treating initial-boundary value problems by finite volume methods is described, in which the parallel between differential and difference arguments is closely maintained. By using intrinsic geometrical properties of the volume elements, it is possible to describe discrete versions of the div, curl, and grad operators which lead, using summation-by-parts techniques, to familiar energy equations as well as the div curl = 0 and curl grad = 0 identities. For the diffusion equation, these operators describe compact schemes whose convergence is assured by the energy equations and which yield both the potential and the flux vector with second order accuracy. A simplified potential form is especially useful for obtaining numerical results by multigrid and alternating direction implicit (ADI) methods. The treatment of general curvilinear coordinates is shown to result from a specialization of these general results.

Rose, Milton E.↗

NASA NTRS1989

Extension-torsion coupling behavior of advanced composite tilt-rotor blades

An analytic model was developed to study the extension-bend-twist coupling behavior of an advanced composite helicopter or tilt-rotor blade. The outer surface of the blade is defined by rotating an arbitrary cross section about an initial twist axis. The cross section can be nonhomogeneous and composed of generally anisotropic materials. The model is developed based upon a three dimensional elasticity approach that is recast as a coupled two-dimensional boundary value problem defined in a curvilinear coordinate system. Displacement solutions are written in terms of known functions that represent extension, bending, and twisting and unknown functions for local cross section deformations. The unknown local deformation functions are determined by applying the principle of minimum potential energy to the discretized two-dimensional cross section. This is an application of the Ritz method, where the trial function family is the displacement field associated with a finite element (8-node isoparametric quadrilaterals) representation of the section. A computer program was written where the cross section is discretized into 8-node quadrilateral subregions. Initially the program was verified using previously published results (both three-dimensional elasticity and technical beam theory) for pretwisted isotropic bars with an elliptical cross section. In addition, solid and thin-wall multi-cell NACA-0012 airfoil sections were analyzed to illustrate the pronounced effects that pretwist, initial twist axis location, and spar location has on coupled behavior. Currently, a series of advanced composite airfoils are being modeled in order to assess how the use of laminated composite materials interacts with pretwist to alter the coupling behavior of the blade. These studies will investigate the use of different ply angle orientations and the use of symmetric versus unsymmetric laminates.

Kosmatka, J. B.↗

NASA NTRS1990

Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls.

Casper, Jay↗

NASA NTRS1990

Numerical simulation of the incompressible internal flow through a tilting disk valve

A numerical simulation of the incompressible viscous flow through a prosthetic tilting disk heart valve is presented in order to demonstrate the current capability to model unsteady flows with moving boundaries. Both steady and unsteady flow calculations are performed by solving the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates. In order to handle the moving boundary problems, the chimera grid embedding scheme which decomposes a complex computational domain into several simple subdomains is used. An algebraic turbulence model for internal flows is incorporated to reach the physiological values of Reynolds number. Good agreement is obtained between the numerical results and experimental measurements. It is found that the tilting disk valve causes large regions of separated flow, and regions of high shear.

Chang, I-Dee↗

NASA NTRS1990

Algorithm development for Maxwell's equations for computational electromagnetism

A new algorithm has been developed for solving Maxwell's equations for the electromagnetic field. It solves the equations in the time domain with central, finite differences. The time advancement is performed implicitly, using an alternating direction implicit procedure. The space discretization is performed with finite volumes, using curvilinear coordinates with electromagnetic components along those directions. Sample calculations are presented of scattering from a metal pin, a square and a circle to demonstrate the capabilities of the new algorithm.

Goorjian, Peter M.↗

NASA NTRS1990

Splitting of inviscid fluxes for real gases

Flux-vector and flux-difference splittings for the inviscid terms of the compressible flow equations are derived under the assumption of a general equation of state for a real gas in equilibrium. No necessary assumptions, approximations for auxiliary quantities are introduced. The formulas derived include several particular cases known for ideal gases and readily apply to curvilinear coordinates. Applications of the formulas in a TVD algorithm to one-dimensional shock-tube and nozzle problems show their quality and robustness.

Liou, Meng-Sing↗

NASA NTRS1990

Preconditioned upwind methods to solve 3-D incompressible Navier-Stokes equations for viscous flows

A computational method for calculating low-speed viscous flowfields is developed. The method uses the implicit upwind-relaxation finite-difference algorithm with a nonsingular eigensystem to solve the preconditioned, three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. The technique of local time stepping is incorporated to accelerate the rate of convergence to a steady-state solution. An extensive study of optimizing the preconditioned system is carried out for two viscous flow problems. Computed results are compared with analytical solutions and experimental data.

Hsu, C.-H.↗