Search NASA⌕ Search

SEARCH · Search NASA

Results for “implicit method”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

The solution of the Navier-Stokes equations using Gauss-Seidel line relaxation

The Navier-Stokes equations in an implicit flux-split difference formulation are solved numerically using a Gauss-Seidel line-relaxation procedure. Particular attention is given to the selection of flux-vector splitting method and flux splitting in boundary layers. Results for sample problems involving (1) turbulent supersonic flow over a cone and (2) the viscous hypersonic flow of a chemically reacting gas in thermal nonequilibrium past a blunted cone are presented in extensive graphs and briefly characterized. The present flux-split procedures are shown to provide accurate shear-layer calculations.

Maccormack, Robert W.↗
NASA NTRS1992

Life extending control for rocket engines

The concept of life extending control is defined. A brief discussion of current fatigue life prediction methods is given and the need for an alternative life prediction model based on a continuous functional relationship is established. Two approaches to life extending control are considered: (1) the implicit approach which uses cyclic fatigue life prediction as a basis for control design; and (2) the continuous life prediction approach which requires a continuous damage law. Progress on an initial formulation of a continuous (in time) fatigue model is presented. Finally, nonlinear programming is used to develop initial results for life extension for a simplified rocket engine (model).

Lorenzo, C. F.↗
NASA NTRS1993

Investigation of high-speed free shear flows using improved pressure-strain correlated Reynolds stress turbulence model

A high-speed shear layer is studied using compressibility corrected Reynolds stress turbulence model which employs newly developed model for pressure-strain correlation. MacCormack explicit prediction-corrector method is used for solving the governing equations and the turbulence transport equations. The stiffness arising due to source terms in the turbulence equations is handled by a semi-implicit numerical technique. Results obtained using the new model show a sharper reduction in growth rate with increasing convective Mach number. Some improvements were also noted in the prediction of the normalized streamwise stress and Reynolds shear stress. The computed results are in good agreement with the experimental data.

Tiwari, S. N.↗
NASA NTRS1994

Accuracy enhancements for overset grids using a defect correction approach

A defect-correction approach is investigated as a means of enhancing the accuracy of flow computations on overset grids. Typically, overset-grid techniques process and pass information only at grid boundaries. In the current approach, error corrections at all overlapped interior points are injected between grids by using a defect-correction scheme. In some cases this is found to enhance the overall accuracy of the overset-grid method. Locally refined overset grids can be used to provide an efficient solution-adaptation method. The defect correction can also be ultilized as an error-correction technique for a coarse grid by evaluating the residual using a fine base grid, but solving the implicit equations only on the coarse grid. Numerical examples include an accuracy and dissipation study of an unsteady decaying vortex flow, the flow over a NACA 0012 airfoil, and the flow over a mulit-element high-lift airfoil.

Rogers, Stuart E.↗
NASA NTRS1995

Computational Study of Flow Establishment in Hypersonic Pulse Facilities

This paper presents a study of the temporal evolution of the combustion flowfield established by the interaction of ram-accelerator-type projectiles with an explosive gas mixture accelerated to hypersonic speeds in an expansion tube. The Navier-Stokes equations for a chemically reacting gas are solved in a fully coupled manner using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order, total variation diminishing (TVD) scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of expansion tube flowfields are then presented for two different configurations. In particular, the development of the shock-induced combustion process is followed. In one case, designed to ensure ignition only in the boundary layer, the lateral extent of the combustion front during the initial transient phase was surprisingly large. The time histories of the calculated thrust and drag forces on the ram accelerator projectile are also presented.

Yungster, S.↗
NASA NTRS1973

Preliminary study of a wall stabilized constructed arc

An iterative, implicit, finite-difference numerical technique is described which is suitable for obtaining solutions to the governing equations for a gas flowing in an axially symmetric constricted-arc heater. The method is shown to provide adequate solutions for three cases of simple pipe flows found in the literature, and for flow in a constricted-arc heater by direct comparison with experimental data. The comparison with arc-heater data includes static pressure, arc voltage, and wall heat flux, all as a function of axial location, and a radial temperature profile at one axial station. The arc-heater data were taken with air as the test gas at a heater inlet pressure approximately 0.40 atm and two flow rates of 2.2 and 4.8 g/sec. The arc currents investigated were between 377 and 584 amperes.

Graves, R. A.↗
NASA NTRS1983

Numerical simulation of electrothermal de-icing systems

Transient simulations of de-icing of composite aircraft components by electrothermal heating have been computed for both one and two-dimensional rectangular geometries. The implicit Crank-Nicolson formulation is used to insure stability of the finite-differenced heat conduction equations and the phase change in the ice layer is simulated using the Enthalpy method. Numerical solutions illustrating de-icer performance for various composite aircraft blades and environmental conditions are presented. Comparisons are made with previous studies and with available experimental data. Initial results using a coordinate mapping technique to describe the actual blade geometry are discussed.

De Witt, K. J.↗
NASA NTRS1983

Numerical simulation of two-dimensional heat transfer in composite bodies with application to de-icing of aircraft components

Transient, numerical simulations of the de-icing of composite aircraft components by electrothermal heating were performed for a two dimensional rectangular geometry. The implicit Crank-Nicolson formulation was used to insure stability of the finite-difference heat conduction equations and the phase change in the ice layer was simulated using the Enthalpy method. The Gauss-Seidel point iterative method was used to solve the system of difference equations. Numerical solutions illustrating de-icer performance for various composite aircraft structures and environmental conditions are presented. Comparisons are made with previous studies. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Chao, D. F. K.↗
NASA NTRS1989

Low-frequency sound propagation modeling over a locally-reacting boundary using the parabolic approximation

There is substantial interest in the analytical and numerical modeling of low-frequency, long-range atmospheric acoustic propagation. Ray-based models, because of frequency limitations, do not always give an adequate prediction of quantities such as sound pressure or intensity levels. However, the parabolic approximation method, widely used in ocean acoustics, and often more accurate than ray models for lower frequencies of interest, can be applied to acoustic propagation in the atmosphere. Modifications of an existing implicit finite-difference implementation for computing solutions to the parabolic approximation are discussed. A locally-reacting boundary is used together with a one-parameter impedance model. Intensity calculations are performed for a number of flow resistivity values in both quiescent and windy atmospheres. Variations in the value of this parameter are shown to have substantial effects on the spatial variation of the acoustic signal.

Robertson, J. S.↗
NASA NTRS1992

Active control of asymmetric vortical flows around cones using injection and heating

The effectiveness of certain active-control methods for asymmetric flows around circular cones is investigated by using computational solution of the unsteady, compressible full Navier-Stokes equations. Two main methods of active control which include flow injection and surface heating are used. For the flow-injection-control method, flow injection is used either in the normal direction to the surface or in the tangential direction to the surface. For the surface-heating-control method, the temperature of the cone surface is increased. The effectiveness of a hybrid method of flow control which combines normal injection with surface heating has also been studied. The Navier-Stokes equations, subjected to various surface boundary conditions, are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme for locally-conical flow solutions.

Kandil, Osama A.↗
NASA NTRS1995

The coupling of fluids, dynamics, and controls on advanced architecture computers

This grant provided for the demonstration of coupled controls, body dynamics, and fluids computations in a workstation cluster environment; and an investigation of the impact of peer-peer communication on flow solver performance and robustness. The findings of these investigations were documented in the conference articles.The attached publication, 'Towards Distributed Fluids/Controls Simulations', documents the solution and scaling of the coupled Navier-Stokes, Euler rigid-body dynamics, and state feedback control equations for a two-dimensional canard-wing. The poor scaling shown was due to serialized grid connectivity computation and Ethernet bandwidth limits. The scaling of a peer-to-peer communication flow code on an IBM SP-2 was also shown. The scaling of the code on the switched fabric-linked nodes was good, with a 2.4 percent loss due to communication of intergrid boundary point information. The code performance on 30 worker nodes was 1.7 (mu)s/point/iteration, or a factor of three over a Cray C-90 head. The attached paper, 'Nonlinear Fluid Computations in a Distributed Environment', documents the effect of several computational rate enhancing methods on convergence. For the cases shown, the highest throughput was achieved using boundary updates at each step, with the manager process performing communication tasks only. Constrained domain decomposition of the implicit fluid equations did not degrade the convergence rate or final solution. The scaling of a coupled body/fluid dynamics problem on an Ethernet-linked cluster was also shown.

Atwood, Christopher↗
NASA NTRS1994

Design of 3-D Nacelle near Flat-Plate Wing Using Multiblock Sensitivity Analysis (ADOS)

One of the major design tasks involved in reducing aircraft drag is the integration of the engine nacelles and airframe. With this impetus, nacelle shapes with and without the presence of a flat-plate wing nearby were optimized. This also served as a demonstration of the 3-D version of the recently developed aerodynamic design optimization methodology using sensitivity analysis, ADOS. The required flow analyses were obtained by solving the three-dimensional, compressible, thin-layer Navier-Stokes equations using an implicit, upwind-biased, finite volume scheme. The sensitivity analyses were performed using the preconditioned version of the SADD scheme (sensitivity analysis on domain decomposition). In addition to demonstrating the present method's capability for automatic optimization, the results offered some insight into two important issues related to optimizing the shapes of multicomponent configurations in close proximity. First, inclusion of the mutual interference between the components resulted in a different shape as opposed to shaping an isolated component. Secondly, exclusion of the viscous effects compromised not only the flow physics but also the optimized shapes even for isolated components.

Eleshaky, Mohamed E.↗
NASA NTRS1987

Evaluation of a research circulation control airfoil using Navier-Stokes methods

The compressible Reynolds time averaged Navier-Stokes equations were used to obtain solutions for flows about a two dimensional circulation control airfoil. The governing equations were written in conservation form for a body-fitted coordinate system and solved using an Alternating Direction Implicit (ADI) procedure. A modified algebraic eddy viscosity model was used to define the turbulent characteristics of the flow, including the wall jet flow over the Coanda surface at the trailing edge. Numerical results are compared to experimental data obtained for a research circulation control airfoil geometry. Excellent agreement with the experimental results was obtained.

Shrewsbury, George D.↗
NASA NTRS1998

Numerical Experiment with Time and Spatial Accuracy of Navier-Stokes Computation For Helicopter Problems

Helicopter flowfields are highly unsteady, nonlinear and three-dimensional. In forward flight and in hover, the rotor blades interact with the tip vortex and wake sheet developed by either itself or the other blades. This interaction, known as blade-vortex interactions (BVI), results in unsteady loading of the blades and can cause a distinctive acoustic signature. Accurate and cost-effective computational fluid dynamic solutions that capture blade-vortex interactions can help rotor designers and engineers to predict rotor performance and to develop designs for low acoustic signature. Such a predictive method must preserve a blade's shed vortex for several blade revolutions before being dissipated. A number of researchers have explored the requirements for this task. This paper will outline some new capabilities that have been added to the NASA Ames' OVERFLOW code to improve its overall accuracy for both vortex capturing and unsteady flows. To highlight these improvements, a number of case studies will be presented. These case studies consist of free convection of a 2-dimensional vortex, dynamically pitching 2-D airfoil including light-stall, and a full 3-D unsteady viscous solution of a helicopter rotor in forward flight In this study both central and upwind difference schemes are modified to be more accurate. Central difference scheme is chosen for this simulation because the flowfield is not dominated by strong shocks. The feature of shock-vortex interaction in such a flow is less important than the dominant blade-vortex interaction. The scheme is second-order accurate in time and solves the thin-layer Navier-Stokes equations in fully-implicit manner at each time-step. The spatial accuracy is either second and fourth-order central difference or third-order upwind difference using Roe-flux and MUSCLE scheme. This paper will highlight and demonstrate the methods for several sample cases and for a helicopter rotor. Preliminary computations on a rotor were performed by using this method and are in the process of documentation.

Ahmad, Jasim↗
NASA NTRS1988

Newton solution of inviscid and viscous problems

The application of Newton iteration to inviscid and viscous airfoil calculations is examined. Spatial discretization is performed using upwind differences with split fluxes. The system of linear equations which arises as a result of linearization in time is solved directly using either a banded matrix solver or a sparse matrix solver. In the latter case, the solver is used in conjunction with the nested dissection strategy, whose implementation for airfoil calculations is discussed. The boundary conditions are also implemented in a fully implicit manner, thus yielding quadratic convergence. Complexities such as the ordering of cell nodes and the use of a far field vortex to correct freestream for a lifting airfoil are addressed. Various methods to accelerate convergence and improve computational efficiency while using Newton iteration are discussed. Results are presented for inviscid, transonic nonlifting and lifting airfoils and also for laminar viscous cases.

Venkatakrishnan, V.↗
NASA NTRS1990

Applicability and evaluation of an implicit self-starting unconditionally stable methodology for the dynamics of structures

The applicability and evaluation of a new self-starting, unconditionally stable, implicit methodology of computation for the dynamics of structures is described. The methodology offers different perspectives and architecture for structural dynamics compared with the traditional (widely advocated and commonly used) time integration methods. It is based on velocity representations and architecture and uses finite elements as the principal analysis tool for structural dynamic modeling/analysis. In particular, the dynamics of beam-type flexural models are considered, and comparative results validate and support the proposed use of the self-starting methodology of computation for the dynamics of linear/nonlinear structures. The overall effectiveness and elegance strongly support its use in most existing commercial codes.

Tamma, Kumar K.↗
NASA NTRS1991

An h-p Taylor-Galerkin finite element method for compressible Euler equations

An extension of the familiar Taylor-Galerkin method to arbitrary h-p spatial approximations is proposed. Boundary conditions are analyzed, and a linear stability result for arbitrary meshes is given, showing the unconditional stability for the parameter of implicitness alpha not less than 0.5. The wedge and blunt body problems are solved with both linear, quadratic, and cubic elements and h-adaptivity, showing the feasibility of higher orders of approximation for problems with shocks.

Demkowicz, L.↗