Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear implicitness”

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 37 records · Page 2

Nonlinear Attitude Filtering Methods

The extended Kalman filter (EKF) is the workhorse of real-time spacecraft attitude estimation. Since the group SO3 of rotation matrices has dimension three, most attitude determination EKFs use lower- dimensional attitude parameterizations than the nine-parameter attitude matrix itself. The fact that all three- parameter representations of SO3 are singular or discontinuous for certain attitudes has led to extended discussions of constraints and attitude representations in EKFs. The most successful EKF uses a nonsingular parameterization for the global attitude, which necessarily has more than three parameters, while employing a three-component representation for the attitude errors. This filter has become known as the Multiplicative Extended Kalman Filter. These issues are now well understood, however, and the EKF has performed admirably in the vast majority of attitude determination applications. Nevertheless, poor performance or even divergence arising from the linearization implicit in the EKF has led to the development of nonlinear filters, most recently sigma point or unscented filters and particle filters.

Markley, F. Landis↗

Development and Verification of the Charring Ablating Thermal Protection Implicit System Solver

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method with first and second order implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton's method, while the fully implicit linear system is solved with the Generalized Minimal Residual method. Verification results from exact solutions and the Method of Manufactured Solutions are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Amar, Adam J.↗

Implicit TVD schemes for hyperbolic conservation laws in curvilinear coordinates

The Harten (1983, 1984) total variation-diminishing (TVD) schemes, constituting a one-parameter explicit and implicit, second-order-accurate family, have the property of not generating spurious oscillations when applied to one-dimensional, nonlinear scalar hyperbolic conservation laws and constant coefficient hyperbolic systems. These methods are presently extended to the multidimensional hyperbolic conservation laws in curvilinear coordinates. Means by which to linearize the implicit operator and solution strategies, in order to improve the computation efficiency of the implicit algorithm, are discussed. Numerical experiments with steady state airfoil calculations indicate that the proposed linearized implicit TVD schemes are accurate and robust.

Yee, H. C.↗

Linear iterative solvers for implicit ODE methods

The numerical solution of stiff initial value problems, which lead to the problem of solving large systems of mildly nonlinear equations are considered. For many problems derived from engineering and science, a solution is possible only with methods derived from iterative linear equation solvers. A common approach to solving the nonlinear equations is to employ an approximate solution obtained from an explicit method. The error is examined to determine how it is distributed among the stiff and non-stiff components, which bears on the choice of an iterative method. The conclusion is that error is (roughly) uniformly distributed, a fact that suggests the Chebyshev method (and the accompanying Manteuffel adaptive parameter algorithm). This method is described, also commenting on Richardson's method and its advantages for large problems. Richardson's method and the Chebyshev method with the Mantueffel algorithm are applied to the solution of the nonlinear equations by Newton's method.

Saylor, Paul E.↗

An implicit adaptation algorithm for a linear model reference control system

This paper presents a stable implicit adaptation algorithm for model reference control. The constraints for stability are found using Lyapunov's second method and do not depend on perfect model following between the system and the reference model. Methods are proposed for satisfying these constraints without estimating the parameters on which the constraints depend.

Mabius, L.↗

Hyperbolic/parabolic development for the GIM-STAR code

Flow fields in supersonic inlet configurations were computed using the eliptic GIM code on the STAR computer. Spillage flow under the lower cowl was calculated to be 33% of the incoming stream. The shock/boundary layer interaction on the upper propulsive surface was computed including separation. All shocks produced by the flow system were captured. Linearized block implicit (LBI) schemes were examined to determine their application to the GIM code. Pure explicit methods have stability limitations and fully implicit schemes are inherently inefficient; however, LBI schemes show promise as an effective compromise. A quasiparabolic version of the GIM code was developed using elastical parabolized Navier-Stokes methods combined with quasitime relaxation. This scheme is referred to as quasiparabolic although it applies equally well to hyperbolic supersonic inviscid flows. Second order windward differences are used in the marching coordinate and either explicit or linear block implicit time relaxation can be incorporated.

Spradley, L. W.↗

Implicit finite-difference simulations of steady and unsteady transonic flows

Implicit methods for several fluid dynamic formulations have been developed and applied to steady-state and low-frequency transonic flows. The basic steps involved in the construction of implicit schemes include: selection of linearly stable accurate implicit difference operators, time-linearization of nonlinear terms, and approximate factorization of the implicit operators into easily solved systems of equations. The proposed schemes are found very efficient for the simpler formulations.

Ballhaus, W. F.↗

Analysis and design of steady transonic flow over airfoils by the method of parametric differentiation

The results reported here are based on applying the method of parametric differentiation (MPD) to transform the nonlinear differential equation governing small-disturbance transonic flow to a linear equation. Implicit approximate factorization and monotone methods were used to accelerate convergence of the linear problem by an order of magnitude over successive line over-relaxation. The relative merits of using MPD are discussed in comparison to conventional small-disturbance applications. Several MPD analyses are performed on an array of airfoils. A design procedure utilizing MPD is discussed and demonstrated for two nonlifting cases.

Halt, D. W.↗

Numerical Behaviour of a Smooth Local Correlation-based Transition Model in a Newton-Krylov Flow Solver

The numerical behaviour of transport-equation-based transition models, including both iterative and grid convergence, is influenced by the source terms. Transition models contain source terms that are large and highly nonlinear, and can be destabilizing in a strong implicit solver. Linearization strategies with varying levels of coupling are evaluated in conjunction with a source-term time step restriction to determine best-practices for solving the SA-sLM2015smooth local correlation-based transition model in an implicit Newton-Krylov flow solver. Achieving deep iterative convergence facilitates a detailed investigation of the grid convergence of these free-transition simulations, which are evaluated relative to fully-turbulent simulations performed using the Spalart-Allmaras turbulence model. Simulations of the NLF0416 general aviation airfoil, VA-2 supercritical airfoil, and NASA CRM-NLF wing-body geometry are performed over a range of grid levels. The results demonstrate that both a fully-coupled linearization strategy and a source-term time step restriction improve nonlinear convergence as the complexity of the free-transition simulations increases. In general, additional grid resolution is required for free-transition simulations relative to fully-turbulent simulations in order to achieve a similar level of accuracy, with the grid convergence of free-transition simulations sensitive to the streamwise grid spacings in the transition regions.

AATT↗

A compressible solution of the Navier-Stokes equations for turbulent flow about an airfoil

A compressible time dependent solution of the Navier-Stokes equations including a transition turbulence model is obtained for the isolated airfoil flow field problem. The equations are solved by a consistently split linearized block implicit scheme. A nonorthogonal body-fitted coordinate system is used which has maximum resolution near the airfoil surface and in the region of the airfoil leading edge. The transition turbulence model is based upon the turbulence kinetic energy equation and predicts regions of laminar, transitional, and turbulent flow. Mean flow field and turbulence field results are presented for an NACA 0012 airfoil at zero and nonzero incidence angles of Reynolds number up to one million and low subsonic Mach numbers.

Shamroth, S. J.↗

The prediction of the turbulent flow field about an isolated airfoil

A compressible time-dependent solution of the Navier-Stokes equations including a transition-turbulence model is obtained for the isolated airfoil flow field problem. The equations are solved by a consistently split linearized block implicit scheme due to Briley and McDonald. A nonorthogonal body fitted coordinate system is used which has maximum resolution near the airfoil surface and in the region of the airfoil leading edge. The transition-turbulence model is based upon the turbulence kinetic energy equation and predicts regions of laminar, transitional and turbulent flow. Mean flow field and turbulence field results are presented for an NACA 0012 airfoil at zero and nonzero incidence angles at Reynolds number up to one million and low subsonic Mach numbers.

Shamroth, S. J.↗

A viscous flow analysis of the tip vortex generation process

A viscous primary-secondary flow analysis designed for numerical solution by forward marching integration is applied to the tip vortex generation problem. Equations governing streamwise momentum, streamwise vorticity, and a gas law are solved in conjunction with equations governing scalar and vector surface potentials for the secondary velocities. The numerical method used combines a consistently split linearized block implicit (LBI) scheme for parabolic equations and a scalar iterative ADI scheme for elliptic equations. Computed results are presented for a wing of rectangular planform immersed in a high Reynolds number stream at 6 degrees incidence. The physical mechanism for generation of the tip vortex, as revealed by the analysis and computations, is discussed. Detailed flow results for the case of 6 degree incidence are included.

Shamroth, S. J.↗

Computation of three-dimensional viscous supersonic flow in inlets

A new approach has been developed for the computation of the three-dimensional viscous supersonic flow with embedded subsonic regions adjacent to solid boundaries and is applied to a mixed-compression supersonic inlet typical of current designs. The approach uses a reduced form of the three-dimensional Navier-Stokes equations so that the resultant equations can be treated as an initial boundary value problem and thus be solved by non-iterative forward marching in space. The numerical procedure utilizes an efficient consistently-split linearized block implicit technique to solve the finite difference analogues to the set of governing partial differential equations.

Buggeln, R. C.↗

Analysis of turbulent flow about an isolated airfoil using a time dependent Navier-Stokes procedure

The procedure solves the Navier-Stokes equations by the consistently split linearized block implicit method of Briley and McDonald in a body fitted coordinate system. The procedure is described and results are presented for flow about an airfoil whose incidence changes from 6 degrees to 19 degrees at a Reynolds number of one million and Mach number of 0.2. In addition, the unsteady flow about an airfoil held at a constant 19 degree incidence is examined and compared to data.

Shamroth, S. J.↗

Numerical solutions of Navier-Stokes equations for compressible turbulent two/three dimensional flows in terminal shock region of an inlet/diffuser

The multidimensional ensemble averaged compressible time dependent Navier Stokes equations in conjunction with mixing length turbulence model and shock capturing technique were used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations were carried out: (1) 2 D, steady, subsonic; (2) 2 D, steady, transonic with normal shock; (3) 2 D, steady, supersonic with terminal shock; (4) 2 D, transient process of shock development and (5) 3 D, steady, transonic with normal shock. The numerical results obtained for the 2 D and 3 D transonic shocked flows were compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data.

Liu, N. S.↗

Numerical solution of the Navier-Stokes equations for compressible turbulent two/three dimensional flows in the terminal shock region of an inlet/diffuser

The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique have been used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations have been carried out: (1) 2-D, steady, subsonic; (2) 2-D, steady, transonic with normal shock, (3) 2-D, steady, supersonic with terminal shock, (4) 2-D, transient process of shock development and (5) 3-D, steady, transonic with normal shock. The numerical results obtained for the 2-D and 3-D transonic shocked flows have been compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data. The predicted transient responses of the flow to externally applied impulsive disturbance are consistent with results obtained via asymptotic analysis.

Liu, N.-S.↗

Solution of 3-dimensional time-dependent viscous flows. Part 3: Application to turbulent and unsteady flows

A numerical scheme is developed for solving the time dependent, three dimensional compressible viscous flow equations to be used as an aid in the design of helicopter rotors. In order to further investigate the numerical procedure, the computer code developed to solve an approximate form of the three dimensional unsteady Navier-Stokes equations employing a linearized block implicit technique in conjunction with a QR operator scheme is tested. Results of calculations are presented for several two dimensional boundary layer flows including steady turbulent and unsteady laminar cases. A comparison of fourth order and second order solutions indicate that increased accuracy can be obtained without any significant increases in cost (run time). The results of the computations also indicate that the computer code can be applied to more complex flows such as those encountered on rotating airfoils. The geometry of a symmetric NACA four digit airfoil is considered and the appropriate geometrical properties are computed.

Weinberg, B. C.↗