Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit solution”

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 325 records · Page 18

A 3D finite element multigrid solver for the Euler equations

A low storage, computationally efficient algorithm for the solution of the compressible Euler equations on unstructured tetrahedral meshes is developed. The algorithm takes the form of a centered scheme with the explicit addition of a high accuracy artificial viscosity and the solution is advanced to steady state by means of a multistage timestepping method. The side-based data structure which is employed enables a clear connection to be established between the proposed algorithm and upwind cell vertex schemes for unstructured meshes. The computational efficiency of the procedure is improved by incorporating an unstructured multigrid acceleration procedure. A number of flows of practical interest are analyzed to demonstrate the numerical performance of the proposed approach.

Peraire, J.↗

Fast and optimal solution to the Rankine-Hugoniot problem

A new, definitive, reliable and fast iterative method is described for determining the geometrical properties of a shock (i.e., theta sub Bn, yields N, V sub s and M sub A), the conservation constants and the self-consistent asymptotic magnetofluid variables, that uses the three dimensional magnetic field and plasma observations. The method is well conditioned and reliable at all theta sub Bn angles regardless of the shock strength or geometry. Explicit proof of uniqueness of the shock geometry solution by either analytical or graphical methods is given. The method is applied to synthetic and real shocks, including a bow shock event and the results are then compared with those determined by preaveraging methods and other iterative schemes. A complete analysis of the confidence region and error bounds of the solution is also presented.

Vinas, A. F.↗

Fast and optimal solution to the 'Rankine-Hugoniot problem'

A new, definitive, reliable and fast iterative method is described for determining the geometrical properties of a shock (i.e., theta sub Bn, yields N, V sub s and M sub A), the conservation constants and the self-consistent asymptotic magnetofluid variables, that uses the three dimensional magnetic field and plasma observations. The method is well conditioned and reliable at all theta sub Bn angles regardless of the shock strength or geometry. Explicit proof of uniqueness of the shock geometry solution by either analytical or graphical methods is given. The method is applied to synthetic and real shocks, including a bow shock event and the results are then compared with those determined by preaveraging methods and other iterative schemes. A complete analysis of the confidence region and error bounds of the solution is also presented.

Vinas, A. F.↗

A simple analytic approximation for dusty Stroemgren spheres

An analytic approximation is illustrated to Stromgren's solution for H II regions which permits explicit exhibition of the effects of internal dust on the ionization structure. Far infrared observations of H II regions are accounted for in terms of true absorption by internal dust of a significant fraction of the Lyman continuum photons.

Petrosian, V.↗

Algebraic Riccati equations in zero-sum differential games

The procedure for finding the closed-loop Nash equilibrium solution of two-player zero-sum linear time-invariant differential games with quadratic performance criteria and classical information pattern may be reduced in most cases to the solution of an algebraic Riccati equation. Based on the results obtained by Willems, necessary and sufficient conditions for existence of solutions to these equations are derived, and explicit conditions for a scalar example are given.

Johnson, T. L.↗

Finite difference methods for transient signal propagation in stratified dispersive media

Explicit difference equations are presented for the solution of a signal of arbitrary waveform propagating in an ohmic dielectric, a cold plasma, a Debye model dielectric, and a Lorentz model dielectric. These difference equations are derived from the governing time-dependent integro-differential equations for the electric fields by a finite difference method. A special difference equation is derived for the grid point at the boundary of two different media. Employing this difference equation, transient signal propagation in an inhomogeneous media can be solved provided that the medium is approximated in a step-wise fashion. The solutions are generated simply by marching on in time. It is concluded that while the classical transform methods will remain useful in certain cases, with the development of the finite difference methods described, an extensive class of problems of transient signal propagating in stratified dispersive media can be effectively solved by numerical methods.

Lam, D. H.↗

Transonic flow over the NACA 64A006 with an oscillating flap-calculations based on the Euler equations

Exploratory calculations of transonic flows over the airfoil with a quarter-chord oscillating flap were made using a program which obtains approximate solutions to the Euler equations with an explicit, shock-capturing, finite-difference scheme. The calculations, essentially inviscid and for the airfoil at zero angle-of-attack in a free-stream, are at Mach numbers and reduced frequencies which were tested in experiments by Tijdeman. The oscillatory lifts from analogus calculations by various investigators generally agree with one another better than they agree with Tijdeman's data. Inclusion in the calculations of an approximate modelling of boundary conditions expected at slotted wind tunnel walls tends to shift some of the results closer to the experimental values.

Magnus, R. J.↗

Effect of sidewall suction on flow in two-dimensional wind tunnels

A closed-form analysis of flow in a two-dimensional subsonic wind tunnel which uses sidewall suction around the model to reduce sidewall boundary-layer effects is presented. The model problem which is treated involves a flat plate airfoil in a tunnel with a suction window shaped to permit an analytic solution. This solution shows that the lift coefficient depends explicitly on the porosity parameter of the suction window and implicitly on the suction pressure differential. For a given sidewall displacement thickness, the lift coefficient increases as the suction-window porosity decreases.

Barnwell, R. W.↗

Adaptive grid refinement for the Euler and compressible Navier-Stokes equations

The incorporation of a simple a posteriori adaptive mesh method into an explicit finite element based procedure for the solution of compressible flow problems is described. The re-interpolation problem is discussed and results are presented that indicate the improvement in solution quality which can be obtained via such an adaptive mesh approach. The method may lead to the appearance of some badly deformed elements which are removed from the domain of computation.

Lohner, R.↗

On a large time-step high resolution scheme

This paper presents a class of new second-order accurate (2K + 3)-point explicit schemes for the computation of weak solutions of hyperbolic conservation laws, that are total-variation-diminishing under a Courant-Friedrichs-Lewy restriction of K. These highly nonlinear schemes are obtained by applying a nonoscillatory first-order accurate (2K + 1)-point scheme to a modified flux. The derived second-order accurate schemes achieve high resolutions, while retaining the robustness of the original first-order accurate scheme.

Harten, A.↗

A simple formulation of the delta-four-stream approximation for radiative transfer parameterizations

A systematic development of the delta-four-stream approximation for calculations of radiative fluxes in planetary atmosphere is presented. It is shown that an analytic solution for this approximation can be derived explicitly, with minimum computational effort for flux calculations. Relative accuracy checks for reflection, transmission, and absorption for numerous asymmetry factors, single-scattering albedos, optical depths, and solar zenith angles have been performed with respect to the 'exact' results computed from the adding method for radiative transfer. Overall, results from the delta-four-stream approximation yield relative accuracies within about 5 percent. This approximation is well suited to radiative transfer parameterizations involving flux and heating calculations in aerosol and cloudy atmospheres.

Liou, Kuo-Nan↗

The simulation of 2D compressible viscous high speed flow by the finite element method

Am implicit/explicit procedure for the solution of problems of two-dimensional steady compressible viscous high-speed flows is presented. In the vicinity of solid walls, a grid which need only exhibit structure in the normal direction is employed while, away from this region, the grid is totally unstructured. The implicit form of the algorithm is used near solid walls, with the grid structure being utilized in an equation solution approach, based upon line relaxation. The explicit form of the algorithm is used elsewhere. Grid adaptation is achieved by means of adaptive remeshing. To illustrate the performance of the proposed method, problems of shock-boundary layer interaction and flow over a simulated forebody at high Mach number are included.

Hassan, O.↗

Progress with multigrid schemes for hypersonic flow problems

Several multigrid schemes are considered for the numerical computation of viscous hypersonic flows. For each scheme, the basic solution algorithm uses upwind spatial discretization with explicit multistage time stepping. Two level versions of the various multigrid algorithms are applied to the two dimensional advection equation, and Fourier analysis is used to determine their damping properties. The capabilities of the multigrid methods are assessed by solving three different hypersonic flow problems. Some new multigrid schemes based on semicoarsening strategies are shown to be quite effective in relieving the stiffness caused by the high aspect ratio cells required to resolve high Reynolds number flows. These schemes exhibit good convergence rates for Reynolds numbers up to 200 x 10(exp 6) and Mach numbers up to 25.

Radespiel, R.↗

Effect of sidewall suction on flow in two-dimensional wind tunnels

A closed-form analysis of flow in a two-dimensional subsonic wind tunnel that uses sidewall suction around the model to reduce sidewall boundary-layer effects is presented. The model problem that is treated involves a flat plate airfoil in a tunnel with a suction window shaped to permit an analytic solution. This solution shows that the lift coefficient depends explicitly on the porosity parameter of the suction window and implicitly on the suction pressure differential. For a given sidewall displacement thickness, the lift coefficient increases as the suction-window porosity decreases.

Barnwell, Richard W.↗

Multigrid schemes for viscous hypersonic flows

Several multigrid schemes are considered for the numerical computation of viscous hypersonic flows. For each scheme, the basic solution algorithm employs upwind spatial discretization with explicit multistage time stepping. Two-level versions of the various multigrid algorithms are applied to the two-dimensional advection equation, and Fourier analysis is used to determine their damping properties. The capabilities of the multigrid methods are assessed by solving two different hypersonic flow problems. Some new multigrid schemes, based on semicoarsening strategies, are shown to be quite effective in relieving the stiffness caused by the high-aspect-ratio cells required to resolve high Reynolds number flows. These schemes exhibit good convergence rates for Reynolds numbers up to 200 x 10(exp 6).

Swanson, R. C.↗

Automation life-cycle cost model

The problem domain being addressed by this contractual effort can be summarized by the following list: Automation and Robotics (A&R) technologies appear to be viable alternatives to current, manual operations; Life-cycle cost models are typically judged with suspicion due to implicit assumptions and little associated documentation; and Uncertainty is a reality for increasingly complex problems and few models explicitly account for its affect on the solution space. The objectives for this effort range from the near-term (1-2 years) to far-term (3-5 years). In the near-term, the envisioned capabilities of the modeling tool are annotated. In addition, a framework is defined and developed in the Decision Modelling System (DEMOS) environment. Our approach is summarized as follows: Assess desirable capabilities (structure into near- and far-term); Identify useful existing models/data; Identify parameters for utility analysis; Define tool framework; Encode scenario thread for model validation; and Provide transition path for tool development. This report contains all relevant, technical progress made on this contractual effort.

Gathmann, Thomas P.↗

An Algebraic Multigrid Solver for Navier-Stokes Problems in the Discrete Second-Order Approximation

An algebraic multigrid scheme is presented for solving the discrete Navier-Stokes equations to second-order accuracy using the defect-correction method. Solutions have been obtained for problems involving both structured and unstructured meshes, with the resolution and resolution grading controlled by global and local mesh refinements. The solver is efficient and robust to the extent that no underrelaxation of variables has been required to ensure convergence, but rates of convergence can be improved with small amounts of underrelaxation of the velocity-pressure coupling. Provided that the computational mesh can resolve the flow field, convergence characteristics are almost mesh independent. Rates of convergence actually improve with refinement, asymptotically approaching mesh independent values. For extremely coarse meshes where dispersive truncation errors would be expected to prevent convergence (or even induce divergence), solutions can still be obtained by using explicit underrelaxation in the iterative cycle.

Webster, R.↗