SEARCH · Engineering Papers
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.
A multi-stage Runge-Kutta method is analyzed for solving the Euler equations exterior to an airfoil. Highly subsonic, transonic and supersonic flows are evaluated. Various techniques for accelerating the convergence to a steady state are introduced and analyzed.
A modified Runge-Kutta integration technique is applied to condensational cloud droplet growth calculations, using the device of limiting dr, the growth of a droplet of radius r over the time step dt to some function of the absolute value of r(eq)-r, where r(eq) denotes the equilibrium value. The use of variable time steps to obtain an efficient choice of dt is also discussed. Results using these techniques are compared with those obtained with the standard fourth-order Runge-Kutta method, showing that the modified technique reduces the requisite computer time by a factor of about one hundred.
A multi-stage Runge-Kutta method is analyzed for solving the two-dimensional Euler equations for external and internal flow problems. Subsonic, supersonic and, highly supersonic flows are studied. Various techniques for accelerating the convergence to a steady state are described and analyzed. Effects of the grid aspect ratio on the rate of convergence are evaluated. An enthalpy damping technique applicable to supersonic flows is described in detail. Numerical results for supersonic flows containing both oblique and normal shocks are presented confirming the efficiency of the method.
A multi-stage Runge-Kutta method is analyzed for solving the two-dimensional Euler equations for external and internal flow problems. Subsonic, supersonic and, highly supersonic flows are studied. Various techniques for accelerating the convergence to a steady state are described and analyzed. Effects of the grid aspect ratio on the rate of convergence are evaluated. An enthalpy damping technique applicable to supersonic flows is described in detail. Numerical results for supersonic flows containing both oblique and normal shocks are presented confirming the efficiency of the method.
Rotordynamic studies of complex systems utilizing multiple degree-of-freedom analysis have been performed to understand response, loads, and stability. In order to understand the fundamental nature of rotordynamic response, the Jeffcott rotor model has received wide attention. The purpose of this paper is to provide a generic rotordynamic analysis of a stiffness cross-coupled Jeffcott rotor model to illustrate characteristics of a second order stiffness-coupled linear system. The particular characteristics investigated were forced response, force vector diagrams, response orbits, and stability. Numerical results were achieved through a fourth order Runge-Kutta method for solving differential equations and the Routh Hurwitz stability criterion. The numerical results were verified to an exact mathematical solution for the steady state response.
A multistage Runge-Kutta method is analyzed for solving the Euler equations exterior to an airfoil. Highly subsonic, transonic and supersonic flows are evaluated. Various techniques for accelerating the convergence to a steady state are introduced and analyzed.
The effects of numerical dissipation upon solutions to the Euler equations are considered, and results for transonic flows past airfoils are presented to demonstrate the effects of the dissipative terms. The equations are approximated using a finite-volume spatial approximation with added dissipation provided by an adaptive mixture of second and fourth differences. The resulting difference equations are solved using either an explicit multistage Runge-Kutta method or a diagonalized implicit method. It is found that errors in surface values can be introduced by the averaging required to calculate derived quantities of interest.
The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.
An improvement on the enthalpy damping procedure currently in use in solving supersonic flow fields is described. A correction based on entropy values is shown to produce a very efficient scheme for simulation of high Mach number three-dimensional flows. Substantial improvements in convergence rates have been achieved by incorporating this enthalpy damping scheme in a finite-volume Runge-Kutta method for solving the Euler equations. Results obtained for blended wing-body geometries at very high Mach numbers are presented.
Finite-volume applications of high-order accurate ENO schemes to two-dimensional boundary-value problems are studied. These schemes achieve high-order spatial accuracy, in smooth regions, by a piecewise polynomial approximation of the solution from cell averages. In addition, this spatial operation involves an adaptive stencil algorithm in order to avoid the oscillatory behavior that is associated with interpolation across steep gradients. High-order TVD Runge-Kutta methods are employed for time integration, thus making these schemes best suited for unsteady problems. Fifth- and sixth-order accurate applications are validated through a grid refinement study involving the solutions of scalar hyperbolic equations. A previously proposed extension for the Euler equations of gas dynamics is tested, including its application to solutions of boundary-value problems involving solid walls and curvilinear coordinates.
The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.
The steady-state transonic viscous flows over two wing-fuselage configurations are solved numerically. The 3D compressible thin-layer Navier-Stokes equations are solved by means of multistage timestepping similar to the Runge-Kutta method. The timestepping technique is applied with a multigrid scheme to make the transonic Navier-Stokes computations using the Baldwin-Lomax turbulence for closure and converging the solution on a grid of 1.52 million grid points with a C-O topology. Both generic and fighter wing-fuselage geometries are solved, and the results are compared to experimental studies. Solution convergence is found to be fast, and the fine grid generates accurate results for both configurations by means of the Baldwin-Lomax turbulence model. Results from the Johnson-King model calculations produce results that are more consistent with experimental results for the fighter configuration.
Numerical experiments for the SONIC schemes on 2D inviscid, compressible, steady, and unsteady problems are presented. These schemes belong to a new class of uniformly second-order accurate nonoscillatory schemes introduced by Huynh, with the well known UNO2 scheme of Harten and Osher being the most 'diffusive' in this class. The SONIC schemes can also be considered as uniformly second order accurate extensions of the popular TVD schemes. For simplicity, a MUSCL approach for spatial discretization and a Runge-Kutta method for time integration are used. Test problems include steady oblique shock reflection and the well known unsteady double Mach reflection problem. Results confirm that the SONIC schemes are more accurate than their TVD counterparts.
The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expressions for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations, and the numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the CFL limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about 6 to 8 grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 to 300 time steps. Three numerical examples of waves in compressible flow are included.
Several explicit integration algorithms with self-adative time integration strategies are developed and investigated for efficiency and accuracy. These algorithms involve the Runge-Kutta second order, the lower Runge-Kutta method of orders one and two, and the exponential integration method. The algorithms are applied to viscoplastic models put forth by Freed and Verrilli and Bodner and Partom for thermal/mechanical loadings (including tensile, relaxation, and cyclic loadings). The large amount of computations performed showed that, for comparable accuracy, the efficiency of an integration algorithm depends significantly on the type of application (loading). However, in general, for the aforementioned loadings and viscoplastic models, the exponential integration algorithm with the proposed self-adaptive time integration strategy worked more (or comparably) efficiently and accurately than the other integration algorithms. Using this strategy for integrating viscoplastic models may lead to considerable savings in computer time (better efficiency) without adversely affecting the accuracy of the results. This conclusion should encourage the utilization of viscoplastic models in the stress analysis and design of structural components.
The objective of this research is to develop turbulence models to predict the flow and heat transfer fields dominated by the curvature effect such as those encountered in turbine cascades and turn-around ducts. A Navier-Stokes code has been developed using an explicit Runge-Kutta method with a two layer k-epsilon/ARSM (Algebraic Reynolds Stress Model), Chien's Low Reynolds Number (LRN) k-epsilon model and Coakley's LRN q-omega model. The near wall pressure strain correlation term was included in the ARSM. The formulation is applied to Favre-averaged N-S equations and no thin-layer approximations are made in either the mean flow or turbulence transport equations. Anisotropic scaling of artificial dissipation terms was used. Locally variable timestep was also used to improve convergence. Detailed comparisons were made between computations and data measured in a turbine cascade by Arts et al. at Von Karman Institute. The surface pressure distributions and wake profiles were predicted well by all the models. The blade heat transfer is predicted well by k-epsilon/ARSM model, as well as the k-epsilon model. It's found that the onset of boundary layer transition on both surfaces is highly dependent upon the level of local freestream turbulence intensity, which is strongly influenced by the streamline curvature. Detailed computation of the flow in the turn around duct has been carried out and validated against the data by Monson as well as Sandborn. The computed results at various streamwise locations both on the concave and convex sides are compared with flow and turbulence data including the separation zone on the inner well. The k-epsilon/ARSM model yielded relatively better results than the two-equation turbulence models. A detailed assessment of the turbulence models has been made with regard to their applicability to curved flows.
The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expression for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations. The numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the Courant-Friedrichs-Lewy (CFL) limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about six - eight grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 - 300 time steps. Three numerical examples of waves in compressible flow are included: (1) sound propagation in a duct with linear shear, (2) linear wave growth in a compressible free shear layer, and (3) vortex pairing in a compressible free shear layer perturbed at two frequencies.
The plotting of streamlines is an effective way of visualizing fluid motion in steady flows. Additional information about the flowfield, such as local rotation and expansion, can be shown by drawing in the form of a ribbon or tube. In this paper, we present efficient algorithms for the construction of streamlines, streamribbons and streamtubes on unstructured grids. A specialized version of the Runge-Kutta method has been developed to speed up the integration of particle paths. We have also derived closed-form solutions for calculating angular rotation rate and radius to construct streamribbons and streamtubes, respectively. According to our analysis and test results, these formulations are two to four times better in performance than previous numerical methods. As a large number of traces are calculated, the improved performance could be significant.