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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Higher order turbulence closure models

Theoretical models are developed and numerical studies conducted on various types of flows including both elliptic and parabolic. The purpose of this study is to find better higher order closure models for the computations of complex flows. This report summarizes three new achievements: (1) completion of the Reynolds-stress closure by developing a new pressure-strain correlation; (2) development of a parabolic code to compute jets and wakes; and, (3) application to a flow through a 180 deg turnaround duct by adopting a boundary fitted coordinate system. In the above mentioned models near-wall models are developed for pressure-strain correlation and third-moment, and incorporated into the transport equations. This addition improved the results considerably and is recommended for future computations. A new parabolic code to solve shear flows without coordinate tranformations is developed and incorporated in this study. This code uses the structure of the finite volume method to solve the governing equations implicitly. The code was validated with the experimental results available in the literature.

Amano, Ryoichi S.↗
NASA NTRS1989

Calculation of steady and unsteady pressures at supersonic speeds with CAP-TSD

A finite difference technique is used to solve the transonic small disturbance flow equation making use of shock capturing to treat wave discontinuities. Thus the nonlinear effects of thickness and angle of attack are considered. Such an approach is made feasible by the development of a new code called CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance), and is based on a fully implicit approximate factorization (AF) finite difference method to solve the time dependent transonic small disturbance equation. The application of the CAP-TSD code to the calculation of low to moderate supersonic steady and unsteady flows is presented. In particular, comparisons with exact linear theory solutions are made for steady and unsteady cases to evaluate shock capturing and other features of the current method. In addition, steady solutions obtained from an Euler code are used to evaluate the small disturbance aspects of the code. Steady and unsteady pressure comparisons are made with measurements for an F-15 wing model and for the RAE tailplane model.

Bennett, Robert M.↗
NASA NTRS1989

Helicopter flight control system design and evaluation for NOE operations using controller inversion techniques

A flight control system design technique is described that allows rapid evaluation of candidate control laws using computer simulated trajectories. The two-step technique first computes the optimum trajectory for the basic aircraft, and then reconstructs the pilot's control displacements necessary to fly this trajectory, given a particular control system architecture. Since the flight control system itself does not add to the acceleration potential of the helicopter, one may evaluate several proposed control laws in terms of the resulting stick inputs for the same maneuver. The method is illustrated through application to implicit and explicit model-following control laws designed for a battlefield helicopter with and without an auxiliary propulsor. It is shown that the technique is a valuable aid for selection of promising control laws for further simulator studies and flight investigations.

Mckillip, Robert M., Jr.↗
NASA NTRS1990

Solution techniques for the vorticity-streamfunction formulation of two-dimensional unsteady incompressible flows

A review of the solution techniques of Tezduyar et al. (1988 and 1989) for the vorticity-streamfunction formulation of two-dimensional incompressible flows is presented. While both the viscous and inviscid cases are considered, the derivation of the proper finite-element formulations for multiply connected domains is emphasized. In all formulations associated with the vorticity transport equation, the streamline upwind/Petrov-Galerkin method is used. The adaptive implicit-explicit and grouped element-by-element solution strategies are employed to maximize the computational efficiency. The solutions obtained in all test cases compare well with solutions from previously published investigations. The convergence and benchmark studies performed in this paper show that the solution techniques presented are accurate, reliable, and efficient.

Tezduyar, T. E.↗
NASA NTRS1992

Using Derivatives Of Higher Order In Sensitivity Analysis

"What if" questions answered exactly, without using "increment-and-reevaluate" approach, and without finite differencing of system analysis, by use of sensitivity analysis to calculate first derivatives of behaviors of coupled systems with respect to design variables. If problem strongly nonlinear, efficiency improved by making derivatives of higher orders available to search algorithms. Extension effected by recursive application of same implicit function theorem that underlies previous method. Useful in design calculations for almost any system, including structures, machines, and distribution networks. Feedback control of processes improved by inclusion of derivative terms of higher order.

Sobieski, Jaroslaw↗
NASA NTRS1992

Upwind differencing and LU factorization for chemical non-equilibrium Navier-Stokes equations

By means of either the Roe or the Van Leer flux-splittings for inviscid terms, in conjunction with central differencing for viscous terms in the explicit operator and the Steger-Warming splitting and lower-upper approximate factorization for the implicit operator, the present, robust upwind method for solving the chemical nonequilibrium Navier-Stokes equations yields formulas for finite-volume discretization in general coordinates. Numerical tests in the illustrative cases of a hypersonic blunt body, a ramped duct, divergent nozzle flows, and shock wave/boundary layer interactions, establish the method's efficiency.

Shuen, Jian-Shun↗
NASA NTRS1992

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗
NASA NTRS1993

Application of fast Fourier transforms to the direct solution of a class of two-dimensional separable elliptic equations on the sphere

An efficient, direct, second-order solver for the discrete solution of a class of two-dimensional separable elliptic equations on the sphere (which generally arise in implicit and semi-implicit atmospheric models) is presented. The method involves a Fourier transformation in longitude and a direct solution of the resulting coupled second-order finite-difference equations in latitude. The solver is made efficient by vectorizing over longitudinal wave-number and by using a vectorized fast Fourier transform routine. It is evaluated using a prescribed solution method and compared with a multigrid solver and the standard direct solver from FISHPAK.

Moorthi, Shrinivas↗
NASA NTRS1993

Incompressible Navier-Stokes Calculations in Pump Flows

Flow through pump components, such as the SSME-HPFTP Impeller and an advanced rocket pump impeller, is efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. The resulting computer code is applied to the flow analysis inside an 11-inch SSME High Pressure Fuel Turbopump impeller, and an advanced rocket pump impeller. Numerical results of SSME-HPFTP impeller flow are compared with experimental measurements. In the advanced pump impeller, the effects of exit and shroud cavities are investigated. Flow analyses at design conditions will be presented.

Kiris, Cetin↗
NASA NTRS1994

Incompressible Navier-Stokes Computations with Heat Transfer

The existing pseudocompressibility method for the system of incompressible Navier-Stokes equations is extended to heat transfer problems by including the energy equation. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. Both forced and natural convection problems are examined. Numerical results from turbulent reattaching flow behind a backward-facing step will be compared against experimental measurements for the forced convection case. The validity of Boussinesq approximation to simplify the buoyancy force term will be investigated. The natural convective flow structure generated by heat transfer in a vertical rectangular cavity will be studied. The numerical results will be compared by experimental measurements by Morrison and Tran.

Kiris, Cetin↗
NASA NTRS1994

Progress In Incompressible Pump Flow Calculations

Steady and unsteady flows for propulsion systems are efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard kappa - epsilon model equations. The resulting computer code is applied to the flow analysis inside an advanced rocket pump impeller in steadily rotating reference frames. Numerical results are compared with experimental measurements. The effects of exit and shroud cavities with the leak-age flow are investigated. Time-accurate incompressible Navier-Stokes formulation with the overlapped grid scheme capability was evaluated by using MIT flapping foil experiment. The grid dependency, turbulence model effects, and the effect of order of differencing were investigated. Numerical results were compared against experimental data. The resulting procedure were applied to unsteady flapping foil calculations. Two upstream NACA 0025 foils perform high frequency synchronized motion and generate unsteady flow conditions to the downstream larger stationary foil. Comparison between unsteady experimental data and numerical results from two different moving boundary procedures will be presented.

Kiris, Cetin↗
NASA NTRS1997

Incompressible Navier-Stokes Computations for the Development of a Ventricular Assist Device

An incompressible flow analysis code, INS3D, has been applied to the development of a mechanical heart assist device. The solution method is based on the artificial compressibility approach and uses an implicit-upwind differencing scheme together with a Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal and the Coriolis force terms are included as source terms. The resulting computational procedure is validated for liquid rocket engine analysis and applied subsequently to analyze a Ventricular Assist Device (VAD). A new design configuration is developed which includes an inducer upstream of the impeller main blades, and substantial improvement is observed in the performance of the VAD.

Kiris, Cetin↗
NASA NTRS2019

An Initial Exploration of Improved Numerics within the Guidelines of the Negative Spalart-Allmaras Turbulence Model

A simple modification to the negative Spalart-Allmaras turbulence model is suggested so that when the turbulence working variable, ~v, is negative, diagonal dominance is increased, as is the tendency for the time-advancement scheme to push e toward positive values. Owing to the fact that the modification is only active when ~v is less than zero, the physical model is left unchanged. Using the proposed modification with a strong implicit solver based on Newton's method, convergence rates can be somewhat improved, with typical reductions in iterations and computer time on the order of 15-50%. The benefits are realized primarily when second- or higher-order accuracy is used for discretizing the convective terms in the turbulence model because of large overshoots that can occur with these schemes at the edges of boundary layers and wakes. For flowfields with few regions of negative ~v, or on very fine meshes where ~v is always greater than zero, little or no benefits should be expected.

Anderson, William K.↗
DOE OSTI2015

Multirate Particle-in-Cell Time Integration Techniques of Vlasov-Maxwell Equations for Collisionless Kinetic Plasma Simulations

A multi-rate PIC formulation was developed that employs large timesteps for slow field evolution, and small (adaptive) timesteps for particle orbit integrations. Implementation is based on a JFNK solver with nonlinear elimination and moment preconditioning. The approach is free of numerical instabilities (ω pe Δt >>1, and Δx >> λ D ), and requires many fewer dofs (vs. explicit PIC) for comparable accuracy in challenging problems. Significant gains (vs. conventional explicit PIC) may be possible for large scale simulations. The paper is organized as follows: Vlasov-Maxwell Particle-in-cell (PIC) methods for plasmas; Explicit, semi-implicit, and implicit time integrations; Implicit PIC formulation (Jacobian-Free Newton-Krylov (JFNK) with nonlinear elimination allows different treatments of disparate scales, discrete conservation properties (energy, charge, canonical momentum, etc.)); Some numerical examples; and Summary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

High-order, linearly stable, partitioned solvers for general multiphysics problems based on implicit–explicit Runge–Kutta schemes

This work introduces a general framework for constructing high-order, linearly stable, partitioned solvers for multiphysics problems from a monolithic implicit–explicit Runge–Kutta (IMEX-RK) discretization of the semi-discrete equations. The generic multiphysics problem is modeled as a system of n systems of partial differential equations where the i th subsystem is coupled to the other subsystems through a coupling term that can depend on the state of all the other subsystems. This coupled system of partial differential equations reduces to a coupled system of ordinary differential equations via the method of lines where an appropriate spatial discretization is applied to each subsystem. The coupled system of ordinary differential equations is taken as a monolithic system and discretized using an IMEX-RK discretization with a specific implicit–explicit decomposition that introduces the concept of a predictor for the coupling term. We propose four coupling predictors that enable the monolithic system to be solved in a partitioned manner, i.e., subsystem-by-subsystem, and preserve the IMEX-RK structure and therefore the design order of accuracy of the monolithic scheme. The four partitioned solvers that result from these predictors are high-order accurate, allow for maximum re-use of existing single-physics software, and two of the four solvers allow the subsystems to be solved in parallel at a given stage and time step. We also analyze the stability of a coupled, linear model problem with a specific coupling structure and show that one of the partitioned solvers achieves unconditional linear stability for this problem, while the others are unconditionally stable only for certain values of the coupling strength. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems including a simple linear system of ODEs, advection–diffusion–reaction systems, fluid–structure interaction problems, and particle-laden flows, where we verify the design order of the IMEX schemes and study various stability properties.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Improving the efficiency of aerodynamic shape optimization

The computational efficiency of an aerodynamic shape optimization procedure that is based on discrete sensitivity analysis is increased through the implementation of two improvements. The first improvement involves replacing a grid-point-based approach for surface representation with a Bezier-Bernstein polynomial parameterization of the surface. Explicit analytical expressions for the grid sensitivity terms are developed for both approaches. The second improvement proposes the use of Newton's method in lieu of an alternating direction implicit methodology to calculate the highly converged flow solutions that are required to compute the sensitivity coefficients. The modified design procedure is demonstrated by optimizing the shape of an internal-external nozzle configuration. Practically identical optimization results are obtained that are independent of the method used to represent the surface. A substantial factor of 8 decrease in computational time for the optimization process is achieved by implementing both of the design procedure improvements.

Burgreen, Greg W.↗
NASA NTRS1997

Mechanical System Reliability and Cost Integration Using a Sequential Linear Approximation Method

The development of new products is dependent on product designs that incorporate high levels of reliability along with a design that meets predetermined levels of system cost. Additional constraints on the product include explicit and implicit performance requirements. Existing reliability and cost prediction methods result in no direct linkage between variables affecting these two dominant product attributes. A methodology to integrate reliability and cost estimates using a sequential linear approximation method is proposed. The sequential linear approximation method utilizes probability of failure sensitivities determined from probabilistic reliability methods as well a manufacturing cost sensitivities. The application of the sequential linear approximation method to a mechanical system is demonstrated.

Kowal, Michael T.↗
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.↗