Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time-accurate”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Numerical simulation and comparison of symmetrical/supercritical airfoils for the near tip region of a helicopter in forward flight

Aerodynamic loads on a multi-bladed helicopter rotor in forward flight at transonic tip conditions are calculated. The unsteady, three-dimensional, time-accurate compressible Reynolds-averaged thin layer Navier-Stokes equations are solved in a rotating coordinate system on a body-conformed, curvilinear grid of C-H topology. Detailed boundary layer and global numerical comparisons of NACA-0012 symmetrical and CAST7-158 supercritical airfoils are made under identical forward flight conditions. The rotor wake effects are modeled by applying a correction to the geometric angle of attack of the blade. This correction is obtained by computing the local induced downwash velocity with a free wake analysis program. The calculations are performed on the Numerical Aerodynamic Simulation Cray 2 and the VPS32 (a derivative of a Cyber 205 at the Langley Research Center) for a model helicopter rotor in forward flight.

Badavi, F. F.↗
NASA NTRS1989

A dual potential formulation of the Navier-Stokes equations

A dual potential formulation for numerically solving the Navier-Stokes equations is developed and presented. The velocity field is decomposed using a scalar and vector potential. Vorticity and dilatation are used as the dependent variables in the momentum equations. Test cases in two dimensions verify the capability to solve flows using approximations from potential flow to full Navier-Stokes simulations. A three-dimensional incompressible flow formulation is also described. An interesting feature of this approach to solving the Navier-Stokes equations is the decomposition of the velocity field into a rotational part (vector potential) and an irrotational part (scalar potential). The Helmholtz decomposition theorem allows this splitting of the velocity field. This approach has had only limited use since it increases the number of dependent variables in the solution. However, it has often been used for incompressible flows where the solution scheme is known to be fast and accurate. This research extends the usage of this method to fully compressible Navier-Stokes simulations by using the dilatation variable along with vorticity. A time-accurate, iterative algorithm is used for the uncoupled solution of the governing equations. Several levels of flow approximation are available within the framework of this method. Potential flow, Euler and full Navier-Stokes solutions are possible using the dual potential formulation. Solution efficiency can be enhanced in a straightforward way. For some flows, the vorticity and/or dilatation may be negligible in certain regions (e.g., far from a viscous boundary in an external flow). It is possible to drop the calculation of these variables then and optimize the solution speed. Also, efficient Poisson solvers are available for the potentials. The relative merits of non-primitive variables versus primitive variables for solution of the Navier-Stokes equations are also discussed.

Gegg, S. G.↗
NASA NTRS1988

Computations of unsteady dissociating nitrogen flows

A time-accurate computer program has been developed to solve flows involving nonequilibrium chemical reactions and vibrational relaxation. The program is based on a multidomain approach and an explicit MacCormack scheme. Results are presented for three simple test cases, involving nitrogen dissociation in a closed box.

Marraffa, Lionel↗
NASA NTRS1988

Explicit and implicit solution of the Navier-Stokes equations on a massively parallel computer

The design, implementation, and performance of a two-dimensional time-accurate Navier-Stokes solver for the CM2 supercomputer are described. The program uses a single processor for each grid point. Two different time-stepping methods have so far been implemented: an explicit third-order Runge-Kutta method and an implicit approximation-factorization method. The CM2 results are checked against those of a mature well-vectorized Cray 2 program, both for correctness and performance. The code is found to be correct, and the performance in some cases is up to several times that of the Cray 2.

Levit, Creon↗
NASA NTRS1989

Unsteady Navier-Stokes computations past oscillating delta wing at high incidence

The unsteady, thin-layer, compressible Navier-Stokes equations, written in the moving frame of reference for the flow relative motion, is solved for the steady and unsteady supersonic flow around a round-edged delta wing. For supersonic flow, local conical flow solution has been obtained from the three-dimensional equations. Pseudotime stepping is used for the steady flow, while time-accurate stepping is used for the unsteady flow. The computational scheme is an implicit approximately-factored finite volume scheme which uses explicit and implicit dissipation terms. The scheme is verified for the steady flow solution. The scheme is then applied to a delta wing undergoing rolling oscillation at a reduced frequency of 1.137 with 15- deg maximum amplitude about a mean angle of attack of 10 deg for a Mach number of 2 and a Reynolds number of 500,000.

Kandil, Osama A.↗
NASA NTRS1989

Numerical simulation of vortex unsteadiness on slender bodies of revolution at large incidence

Time-accurate, fine-grid Navier-Stokes solutions were obtained for flow over a slender ogive-cylinder body of revolution at angles of attack ranging from 10 deg to 40 deg. The results indicate the progressive growth of crossflow separation and the development of the leeward side vortex pattern with increasing incidence. The computed flows show good agreement with experimental measurements. As the angle of attack was increased, the flows become less damped, and at 40 deg a nonsteady flow exhibiting self-sustained fluctuations was observed. The nonsteadiness was linked to the presence of small-scale three-dimensional vortices moving along the primary surfaces of crossflow separation. The behavior of the fluctuations with incidence parallels the trend observed in experiments.

Schiff, Lewis B.↗
NASA NTRS1989

Grid refinement studies of turbine rotor-stator interaction

Results from a three-dimensional, time-accurate Navier-Stokes simulation of rotor-stator interaction in an axial turbine stage are presented. The present study uses a fine grid in the spanwise direction to better resolve endwall and tip cllearance effects and complements coarse-grid calculations that were reported earlier. A realistic turbine stage with 22 stator vanes and 28 rotor blades is simulated as a single-stator, single-rotor airfoil combination with the stator geometry modified to properly account for blockage effects. This is in contrast to the earlier coarse-grid calculations where the rotor geometry was modified. The improved grid resolution and the unmodified rotor geometry result in a more accurate simulation of the flow field, particularly in the rotor channel where the interaction effects are more severe. The numerical results are compared to experimental data wherever possible and to earlier calculations.

Madavan, N. K.↗
NASA NTRS1989

Numerical solution of unsteady incompressible viscous flows in generalized moving coordinate systems

A solution method of the time-accurate, incompressible Navier-Stokes equations in generalized curvilinear moving coordinate systems is presented in this paper. Accuracy is achieved by a conservative finite-volume discretization which satisfies the geometric conservation laws in generalized moving coordinate systems. The solution method is second-order accurate in space and first-order accurate in time. A fractional step solution method is used to efficiently solve the discrete equations. The unknowns, namely the pressure and the volume-fluxes, are chosen to facilitate the formulation of a consistent Poisson equation and to obtain a robust Poisson solver with favorable convergence properties. The method is validated by comparisons to other numerical and experimental solutions. The comparisons show good agreement.

Rosenfeld, Moshe↗
NASA NTRS1989

Numerical simulation of vortical flows on flexible wings

A procedure to simultaneously solve the Navier-Stokes equations and modal structural equations of motion is presented for computing aeroelastic responses of wings. The Navier-Stokes flow equations are solved by a finite-difference scheme with dynamic grids. The coupled aeroelastic equations of motion are solved using the linear-acceleration method. The aeroelastic, configuration-adaptive dynamic grids are time-accurately generated using the aeroelastically deformed shape of the wing. The unsteady flow calculations are validated with the experiment. Present development is demonstrated for computing vortical flows over flexible wings.

Guruswamy, Guru P.↗
NASA NTRS1988

Navier-Stokes simulations of tip vortices for fixed and rotating helicopter blades

Flowfield and tip vortex results are presented for a hovering rotor blade at subcritical and supercritical flow conditions for both nonlifting and lifting configurations. These results are calculated numerically by solving, in a time-accurate fashion, the unsteady thin layer Navier-Stokes equations written in rotor coordinates. The lifting calculations use an induced downwash correction, estimated from a free wake analysis, to the geometric angle of attack of the blades to account for the wake effect. Comparison of numerical results with the experimental data shows very good agreement for all cases considered. Alternate methods of calculating hovering rotor flowfield as steady state flowfield on isolated fixed-blade that have the same circulation distribution as that of rotor in hover are also explored.

Srinivasan, G. R.↗
NASA NTRS1988

A numerical simulation of finite-length Taylor-Couette flow

Results from numerical simulations of finite-length Taylor-Couette flow are presented. Included are time-accurate and steady-state studies of the change in the nature of the symmetric two-cell/asymmetric one-cell bifurcation with varying aspect ratio and of the Reynolds number/aspect ratio locus of the two-cell/four-cell bifurcation. Preliminary results from wavy-vortex simulations at low aspect ratios are also presented.

Streett, C. L.↗
NASA NTRS1989

Vortical flow computations on swept flexible wings using Navier-Stokes equations

A procedure to couple the Navier-Stokes solutions with modal structural equations of motion is presented for computing aeroelastic responses of swept flexible wings. The Navier-Stokes flow equations are solved by a finite-difference scheme with dynamic grids. The coupled aeroelastic equations of motion are solved using the linear-acceleration method. The configuration-adaptive dynamic grids are time-accurately generated using the aeroelastically deformed shape of the wing. The calculations are compared with the experiment when available. Effects of flexibility, sweep angle, and pitch rate are demonstrated for flows with vortices.

Guruswamy, Guru P.↗
NASA NTRS1989

Calculation of unsteady flows in turbomachinery using the linearized Euler equations

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗
NASA NTRS1989

A time accurate finite volume method for propulsion chamber flows

An implicit three-dimensional time-accurate method for propulsion chamber flows is proposed which uses line Gauss-Seidel relaxation and multiple axial sweeps for the convergence of each time step. The general time-integration algorithm employed includes such schemes as the Euler implicit method. The results of spatial and temporal accuracy tests reveal that Roe's (1981) flux difference splitting provides excellent tracking of acoustic wave speeds. In comparison with other methods, no low mean flow Mach number convergence limitation or Courant number stabilization restriction is observed.

Beddini, R. A.↗
NASA NTRS1989

Application of a self-adaptive grid method to complex flows

A directional-split, modular, user-friendly grid point distribution code is applied to several test problems. The code is self-adaptive in the sense that grid point spacing is determined by user-specified constants denoting maximum and minimum grid spacings and constants relating the relative influence of smoothness and orthogonality. Estimates of truncation error, in terms of flow-field gradients and/or geometric features, are used to determine the point distribution. Points are redistributed along grid lines in a specified direction in an elliptic manner over a user-specified subdomain, while orthogonality and smoothness are controlled in a parabolic (marching) manner in the remaining directions. Multidirectional adaption is achieved by sequential application of the method in each coordinate direction. The flow-field solution is redistributed onto the newly distributed grid points after each unidirectional adaption by a simple one-dimensional interpolation scheme. For time-accurate schemes such interpolation is not necessary and time-dependent metrics are carried in the fluid dynamic equations to account for grid movement.

Deiwert, G. S.↗
NASA NTRS1989

Unsteady aerodynamic simulation of multiple bodies in relative motion: A prototype method

A prototype method for time-accurate simulation of multiple aerodynamic bodies in relative motion is presented. The method is general and features unsteady chimera domain decomposition techniques and an implicit approximately factored finite-difference procedure to solve the time-dependent thin-layer Navier-Stokes equations. The method is applied to a set of two- and three- dimensional test problems to establish spatial and temporal accuracy, quantify computational efficiency, and begin to test overall code robustness.

Meakin, Robert L.↗
NASA NTRS1988

Numerical investigations in three-dimensional internal flows

In order to develop an understanding of flow fields that exist in sidewall compression inlet models, an explicit, time-accurate full Navier-Stokes code was used to model them. Both laminar and turbulent boundary layers were numerically simulated. These simulations were examined as they were influenced by Mach number, changes in sidewall compression angle, hypersonic viscous interaction effects, and thick entering boundary layers. Cause and effect relationships were established. The numerical simulations were compared with the limited data available to aid in understanding both the simulations and the experiments.

Rose, William C.↗
NASA NTRS1989

Numerical modeling of enclosure convection

A numerical study on the steady and unsteady natural convection in two-dimensional rectangular enclosures has been performed by a time-accurate ADI finite difference scheme. The study covered a range of Rayleigh numbers between 1000 and 10 to the 7th, aspect ratios between 0.2 and 10.0, and tilt angles between -90 (heating from bottom) and +90 deg (heating from top). Various Prandtl numbers have been studied, but only the results of water (Pr = 7.0) are reported here due to space limitations. The physics revealed, however, includes the convection phenomena and the Rayleigh-Benard stability, as well as the combined mechanism of these two. The onset of secondary cells is determined by using a velocity map, which is simpler and cleaner, instead of a streamline plot. The critical Ra number for the occurrence of these secondary cells is shown to be lower than can be detected by experimental studies. On the Rayleigh-Benard stability part, a second transition from stable single-cell convection to periodic multicellular convection is disclosed.

Duh, J. C.↗