SEARCH · Search NASA

Results for “Compressible Euler equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1988

On the coupling of hyperbolic and parabolic systems: Analytical and numerical approach

The coupling of hyperbolic and parabolic systems is discussed in a domain Omega divided into two distinct subdomains omega(+) and omega(-). The main concern is to find the proper interface conditions to be fulfilled at the surface separating the two domains. Next, they are used in the numerical approximation of the problem. The justification of the interface conditions is based on a singular perturbation analysis, i.e., the hyperbolic system is rendered parabolic by adding a small artifical viscosity. As this goes to zero, the coupled parabolic-parabolic problem degenerates into the original one, yielding some conditions at the interface. These are taken as interface conditions for the hyperbolic-parabolic problem. Actually, two alternative sets of interface conditions are discussed according to whether the regularization procedure is variational or nonvariational. It is shown how these conditions can be used in the frame of a numerical approximation to the given problem. Furthermore, a method of resolution is discussed which alternates the resolution of the hyperbolic problem within omega(-) and of the parabolic one within omega(+). The spectral collocation method is proposed, as an example of space discretization (different methods could be used as well); both explicit and implicit time-advancing schemes are considered. The present study is a preliminary step toward the analysis of the coupling between Euler and Navier-Stokes equations for compressible flows.

Gastaldi, Fabio↗

NASA NTRS1992

LEWICE/E: An Euler based ice accretion code

A new version of the LEWICE ice accretion computer code was developed which calculates the ice growth on two dimensional surfaces, incorporating the effects of compressibility through the solution of the Euler equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, and plot results. This code increases the applicability of ice accretion predictions by allowing calculations at higher Mach numbers. The new elements of the code are described. Calculated results are compared to experiment for several cases, including a LEWICE example case and a thin airfoil section at a Mach number of 0.58.

Potapczuk, Mark G.↗

NASA NTRS1992

LEWICE/E - An Euler based ice accretion code

Potapczuk, Mark G.↗

NASA NTRS1995

A Hybrid Parallel Preconditioning Algorithm For CFD

A new hybrid preconditioning algorithm will be presented which combines the favorable attributes of incomplete lower-upper (ILU) factorization with the favorable attributes of the approximate inverse method recently advocated by numerous researchers. The quality of the preconditioner is adjustable and can be increased at the cost of additional computation while at the same time the storage required is roughly constant and approximately equal to the storage required for the original matrix. In addition, the preconditioning algorithm suggests an efficient and natural parallel implementation with reduced communication. Sample calculations will be presented for the numerical solution of multi-dimensional advection-diffusion equations. The matrix solver has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The full paper will show numerous examples in CFD to demonstrate the efficiency and robustness of the method.

Barth,Timothy J.↗

NASA NTRS1996

Some Experiences with Nonoverlapping Schur Complement Parallel Preconditioning for CFD Calculations

In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Barth, Timothy J.↗

NASA NTRS1983

Compressive behavior of unidirectional fibrous composites

The longitudinal compressive behavior of unidirectional fiber composites was investigated by using the Illinois Institute of Technology Research Institute (IITRI) test method with thick and thin test specimens. The test data obtained are interpreted by means of stress/strain curves from back-to-back strain gages, examination of fracture surfaces by scanning electron microscope, and predictive equations for distinct failure modes including fiber compression failure. Euler buckling, delamination, and flexure. The results show that longitudinal compressive fracture is induced by a combination of delamination, flexure, and fiber tier breaks. No distinct fracture surface characteristics can be associated with unique failure modes. An equation is described that can be used to extract the longitudinal compressive strength from the longitudinal tensile and flexural strengths of the same composite system.

Sinclair, J. H.↗

NASA NTRS1982

Compression behavior of unidirectional fibrous composite

The longitudinal compression behavior of unidirectional fiber composites is investigated using a modified Celanese test method with thick and thin test specimens. The test data obtained are interpreted using the stress/strain curves from back-to-back strain gages, examination of fracture surfaces by scanning electron microscope, and predictive equations for distinct failure modes including fiber compression failure, Euler buckling, delamination, and flexure. The results show that the longitudinal compression fracture is induced by a combination of delamination, flexure, and fiber tier breaks. No distinct fracture surface characteristics can be associated with unique failure modes. An equation is described which can be used to extract the longitudinal compression strength knowing the longitudinal tensile and flexural strengths of the same composite system.

Sinclair, J. H.↗

DOE OSTI2025

Characterizing Artificial Viscosity Parameters with Approximate Symmetries

We are often faced with trying to capture the physics of compressible shocks, which are governed by the Euler equations. However, Euler shocks are formally discontinuous at the shock front (translating to a step-function behavior of rel evant flow variables). This poses a practical problem for codes with finite-sized grid elements. As a result, one must make a concession in simulating the be havior of shocks within a discrete framework. In particular, we must blur, or ‘regularize’ Euler shocks so that they may be captured on a finite grid.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1984

Transonic viscous-inviscid interaction using Euler and inverse boundary-layer equations

The paper is concerned with the use of a zonal method for the computation of transonic viscous-inviscid interacting flow about airfoils. The inviscid portion of the flow is treated by using an Euler equation solution method, while an inverse integral compressible turbulent boundary-layer solution method is used for the viscous portion of the flow. The matching of the viscous and inviscid solutions is discussed, and some numerical results as well as comparisons with experimental data are presented.

Whitfield, D. L.↗

NASA NTRS1994

Multiblock, Multigrid Solution Of Euler Equations

Method of numerical solution of Euler equations of three-dimensional flow of compressible fluid involves combination of multiblock and multigrid strategies. In multiblock strategy, flow field divided, into multiple smaller, more computationally-convenient zones and computational grid fitted to applicable flow boundaries generated in each block. In multigrid strategy used here, different quantities computed, variously, on finer or coarser grids. Minimizing cost of computation by using fewest grid points yielding acceptably accurate values of affected variable. Multigrid strategy found effective in accelerating convergence to steady state, while multiblock strategy provides geometric flexibility.

Melson, N. Duane↗

NASA NTRS1994

Numerical Simulation of the Instability of an Inviscid Liquid Jet in a Coflowing Compressible Airstream

The nonlinear interfacial instability of a liquid jet in a coflowing compressible airstream is studied numerically. A high-resolution scheme which has second-order accuracy in space and time is coupled with a Lagrangian marker particle algorithm to visualize the large-scale motion of the interfaces in compressible flow. A numerical algorithm based on an approximate equation of state of a compressible liquid is developed to allow this two-fluid system to be governed by the nonlinear unsteady Euler equations in conservative form. The initial growth of small disturbances given by the simulations agrees well with linear theory. The process of jet disruption in compressible flow is demonstrated to consist of the formation of liquid spikes, interweaving of the gas and liquid and stretching and detachment of the liquid main center core.

Li, H. -S.↗

NASA NTRS1988

Accurate multigrid solution of the Euler equations on unstructured and adaptive meshes

A method for accurately solving inviscid compressible flow in the subcritical and supercritical regimes about complex configurations is presented. The method is based on the use of unstructured triangular meshes in two dimensions, and special emphasis is placed on the accuracy and efficiency of the solutions. High accuracy is achieved by careful scaling of the artificial dissipation terms, and by reformulating the inner and outer boundary conditions for both the convective and dissipative operators. An adaptive grid refinement strategy is presented which enhances the solution accuracy for complex flows. When coupled with an unstructured multigrid algorithm, this method is shown to produce an efficient solver for flows about arbitrary configurations.

Mavriplis, Dimitri J.↗

NASA NTRS1988

Accurate multigrid solution of the Euler equations on unstructured and adaptive meshes

Mavriplis, Dimitri J.↗

NASA NTRS1978

Lagrangian computation of inviscid compressible flows

A Lagrangian method is developed to solve the Euler equations of gas dynamics. The solution of the equations is obtained by a numerical computation with the well-known Flux-Corrected-Transport (FCT) numerical method. This procedure is modified so that the boundary treatment is accurate and relatively simple. Shock waves and other flow discontinuities are captured monotonically without any type of fitting procedures. The Lagrangian method is employed so that the problem of mesh generation is completely avoided. The method is applicable to all Mach numbers except the low subsonic range where compressibility effects are small. The method is applied to a one-dimensional Riemann problem (shock tube) and to a two-dimensional supersonic channel flow with reflecting shock waves.

Klopfer, G. H.↗

NASA NTRS1988

Three-dimensional unsteady transonic viscous-inviscid interaction using the Euler and boundary-layer equations

The objective of this study is the development of a numerical technique which can provide three-dimensional, time-accurate, compressible, turbulent flow solutions in a practical and relatively economical way. The approach taken is that of the method of viscous-inviscid interaction. The Euler equations are assumed to govern the outer inviscid portion of the flow, and the viscous layer close to the solid wall is described by a set of integral boundary-layer equations. The viscous solutions are obtained in a direct fashion with a weighted-average phase error scheme. The method of equivalent sources is used for viscous-inviscid coupling. Steady-state and unsteady computations for an AGARD airfoil and a wing show that satisfactory engineering solutions can be obtained for attached, high Reynolds number flows using this method. Quasi-unsteady interactions are shown to produce similar results to those provided by true-unsteady interactions. Considerable computer resources can be saved for some cases by using quasi-unsteady interactions.

Whitfield, David L.↗

NASA NTRS1984

Time dependent computation of the Euler equations for designing 2-D cascades, including the case of transonic shock free design

A numerical method to solve the inverse problem for airfoil cascades is described. The case of compressible, inviscid flow is considered. The proposed method is based on the finite difference approximation of the time-dependent Euler equations. Numerical examples show the capability of the method to design subsonic or transonic shockless blades.

Zannetti, L.↗

NASA NTRS1993

Directionally adaptive finite element method for multidimensional Euler and Navier-Stokes equations

A directionally adaptive finite element method for multidimensional compressible flows is presented. Quadrilateral and hexahedral elements are used because they have several advantages over triangular and tetrahedral elements. Unlike traditional methods that use quadrilateral/hexahedral elements, our method allows an element to be divided in each of the three directions in 3D and two directions in 2D. Some restrictions on mesh structure are found to be necessary, especially in 3D. The refining and coarsening procedures, and the treatment of constraints are given. A new implementation of upwind schemes in the constrained finite element system is presented. Some example problems, including a Mach 10 shock interaction with the walls of a 2D channel, a 2D viscous compression corner flow, and inviscid and viscous 3D flows in square channels, are also shown.

Tan, Zhiqiang↗

NASA NTRS1991

h-p adaptive finite element methods in computational fluid dynamics

The principal ideas of h-p adaptive finite element methods for fluid dynamics problems are discussed. Applications include acoustics, compressible Euler and both compressible and incompressible Navier-Stokes equations. Several numerical examples illustrate the presented concepts.

Oden, J. T.↗