Search NASA⌕ Search

SEARCH · Search NASA

Results for “Euler characteristic”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

A fast efficient implicit scheme for the gasdynamic equations using a matrix reduction technique

An efficient implicit finite-difference algorithm for the gasdynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved without loss of the stability characteristics generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows and used to show that the method offers computational advantages over the conventional Beam and Warming scheme. Existing Beam and Warming codes can be retrofit with minimal effort. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations.

Barth, T. J.↗
NASA NTRS1993

Energy transfer and constrained simulations in isotropic turbulence

The defining characteristic of turbulent flows is their ability to dissipate energy, even in the limit of zero viscosity. The Euler equations, if constrained in such a way that the velocity derivatives remain bounded, conserve energy. But when they arise as the limit of the Navier-Stokes (NS) equations, when the Reynolds number goes to infinity, there is persuasive empirical evidence that the gradients become singular as just the right function of Re for the dissipation to remain non-zero and to approach a well defined limit. It is generally believed that this limiting value of the dissipation is a property of the Euler equations themselves, independent of the particular dissipative mechanism involved, and that it can be normalized with the large scale properties of the turbulent flow (e.g. the kinetic energy per unit volume u'(exp 2)/2, and the integral scale L) without reference to the Reynolds number or to other dissipative quantities. This is usually taken to imply that the low wave number end of the energy spectrum, far from the dissipative range, is also independent of the particular mechanism chosen to dispose of the energy transfer. In the following sections, we present some numerical experiments on the effect of substituting different dissipation models into the truncated Euler equations. We will see that the effect is mainly felt in the 'near dissipation' range of the energy spectrum, but that this range can be quite wide in some cases, contaminating a substantial range of wave numbers. In the process, we will develop a 'practical' approximation to the subgrid energy transfer in isotropic turbulence, and we will gain insight into the structure of the nonlinear interactions among turbulent scales of comparable size, and into the nature of energy backscatter. Some considerations on future research directions are offered at the end.

Jimenez, Javier↗
NASA NTRS1995

Implicit schemes and parallel computing in unstructured grid CFD

The development of implicit schemes for obtaining steady state solutions to the Euler and Navier-Stokes equations on unstructured grids is outlined. Applications are presented that compare the convergence characteristics of various implicit methods. Next, the development of explicit and implicit schemes to compute unsteady flows on unstructured grids is discussed. Next, the issues involved in parallelizing finite volume schemes on unstructured meshes in an MIMD (multiple instruction/multiple data stream) fashion are outlined. Techniques for partitioning unstructured grids among processors and for extracting parallelism in explicit and implicit solvers are discussed. Finally, some dynamic load balancing ideas, which are useful in adaptive transient computations, are presented.

Venkatakrishnam, V.↗
NASA NTRS2004

Computed Linear/Nonlinear Acoustic Response of a Cascade for Single/Multi Frequency Excitation

This paper examines mode generation and propagation characteristics of a 2-D cascade due to incident vortical disturbances using a time domain approach. Full nonlinear Euler equations are solved employing high order accurate spatial differencing and time marching techniques. The solutions show the generation and propagation of mode orders that are expected from theory. Single frequency excitations show linear response over a wide range of amplitudes. The response for multi-frequency excitations tend to become nonlinear due to interaction between frequencies and self interaction.

Nallasamy, M.↗
NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗
NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗
NASA NTRS1984

An efficient approximate factorization implicit scheme for the equations of gasdynamics

An efficient implicit finite-difference algorithm for the gas dynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved while maintaining the same favorable stability characteristics and generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows about a NACA 0012 airfoil. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations. It is shown that the method offers advantages over the conventional Beam and Warming scheme and can retrofit existing Beam and Warming codes with minimal effort.

Barth, T. J.↗
NASA NTRS1992

CFD predictions of near-field pressure signatures of a low-boom aircraft

A three dimensional Euler marching code has been utilized to predict near-field pressure signatures of an aircraft with low boom characteristics. Computations were extended to approximately six body lengths aft of the aircraft in order to obtain pressure data at three body lengths below the aircraft for a cruise Mach number of 1.6. The near-field pressure data were extrapolated to the ground using a Whitham based method. The distance below the aircraft where the pressure data are attained is defined in this paper as the 'separation distance.' The influences of separation distances and the still highly three-dimensional flow field on the predicted ground pressure signatures and boom loudness are presented in this paper.

Fouladi, Kamran↗
NASA NTRS1993

An aerodynamic model for one and two degree of freedom wing rock of slender delta wings

The unsteady aerodynamic effects due to the separated flow around slender delta wings in motion were analyzed. By combining the unsteady flow field solution with the rigid body Euler equations of motion, self-induced wing rock motion is simulated. The aerodynamic model successfully captures the qualitative characteristics of wing rock observed in experiments. For the one degree of freedom in roll case, the model is used to look into the mechanisms of wing rock and to investigate the effects of various parameters, like angle of attack, yaw angle, displacement of the separation point, and wing inertia. To investigate the roll and yaw coupling for the delta wing, an additional degree of freedom is added. However, no limit cycle was observed in the two degree of freedom case. Nonetheless, the model can be used to apply various control laws to actively control wing rock using, for example, the displacement of the leading edge vortex separation point by inboard span wise blowing.

Hong, John↗
NASA NTRS1994

Variations in the modal characteristics of a telescopically deploying beam

The equations of motion for a two-segment deploying telescopic beam are derived through application of Lagrange's equation. The outer tube of the beam is fixed at one end and the inner tube slides freely relative to the fixed segment. The resulting nonlinear, non-autonomous set of equations is linearized and simplified to the standard Euler-Bernoulli partial differential equations for an elastic beam by freezing the deployment process at various stages of deployment, and examining the small amplitude and natural modes of vibration of the resulting configuration. Application of the natural boundary conditions and compatibility of motion relations for the two segments in their common region of overlap leads to a transcendental characteristic equation in the frequency parameter Beta(L). Numerical solution of the equation for the characteristic roots determines the modal frequencies, and the corresponding mode shapes are obtained from the general solution of the Euler-Bernoulli equation tailored to the natural boundary conditions. Sample results of modal frequencies and shapes are presented for various stages of deployment and discussed. It is shown that for all intermediate stages of deployment (between 0 and 100 percent) the spectral distribution is drastically altered by the appearance of regions of very closely spaced modal frequencies. The sources of this modal agglomeration are explored.

Amos, Anthony K.↗
NASA NTRS2013

Combined Experimental and Numerical Simulations of Thermal Barrier Coated Turbine Blades Erosion

A combined experimental and computational study was conducted to investigate the erosion of thermal barrier coated (TBC) blade surfaces by alumina particles ingestion in a single stage turbine. In the experimental investigation, tests of particle surface interactions were performed in specially designed tunnels to determine the erosion rates and particle restitution characteristics under different impact conditions. The experimental results show that the erosion rates increase with increased impingement angle, impact velocity and temperature. In the computational simulations, an Euler-Lagrangian two stage approach is used in obtaining numerical solutions to the three-dimensional compressible Reynolds Averaged Navier-Stokes equations and the particles equations of motion in each blade passage reference frame. User defined functions (UDF) were developed to represent experimentally-based correlations for particle surface interaction models which were employed in the three-dimensional particle trajectory simulations to determine the particle rebound characteristics after each surface impact. The experimentally based erosion UDF model was used to predict the TBC erosion rates on the turbine blade surfaces based on the computed statistical data of the particles impact locations, velocities and angles relative to the blade surface. Computational results are presented for the predicted TBC blade erosion in a single stage commercial APU turbine, for a NASA designed automotive turbine, and for the NASA turbine scaled for modern rotorcraft operating conditions. The erosion patterns in the turbines are discussed for uniform particle ingestion and for particle ingestion concentrated in the inner and outer 5 percent of the stator blade span representing the flow cooling the combustor liner.

Hamed, Awate↗
NASA NTRS1989

A natural flow wing design employing 3-D nonlinear analysis applied at supersonic speeds

A wing-design study has been conducted on a 65-deg-swept leading-edge delta wing in which a near-conical geometry was employed to take advantage of the naturally occurring conical flow which arises over such a wing in a supersonic flow field. Three-dimensional nonlinear analysis methods were used in the study. In preliminary design, wing planform, design conditions, and near-conical concept were derived and a baseline standard wing (conventional airfoil distribution) and a baseline near-conical wing were chosen. During the initial analysis, a full-potential solver was employed to determine the aerodynamic characteristics of the baseline standard delta wing and the near-conical delta wing. Modifications due to airfoil thickness, leading-edge radius, and camber were then applied to the baseline near-conical wing. The final design employed a Euler solver to analyze the best wing configurations found in the initial design, and to extend this study to develop a more refined wing. Benefits due to each modification are discussed, and a final natural flow wing geometry is chosen and its aerodynamic characteristics are compared with the baseline wings.

Bauer, Steven X. S.↗
NASA NTRS1992

Computational fluid dynamics application to hypersonic flow over a Martian entry vehicle (MEV) - A correlation with experiment

A class of vehicles for a mission to Mars are analyzed for aerodynamic characteristics using advanced Computational Fluid Dynamics (CFD). The general configuration is a modified cone-conical-frustum geometry where the nose radius has a large influence on the flowfield. Inviscid-compressible flow using the Euler equations and viscous-compressible flow using the thin-layer Navier-Stokes equations is applied to the configuration. The surface modeling, grid generation and application of state-of-the-art CFD software are described. The effects of nose radius, angle of attack, and hypersonic velocity on the flight characteristics of the vehicle are discussed. The numerical simulations demonstrate the merits of the inviscid and viscous software. Results are compared with wind tunnel experiments.

Stewart, John E.↗
NASA NTRSDate not supplied

The Initial Orbit Determination (IOD) Problem With Range, Range-Rate and Angles

The Initial Orbit Determination (IOD) problem has been tacked since Laplace, Euler and Gauss in the service of determining the orbits of planets, asteroids, and comets. The measurements that established the orbit characteristics were angles over long arcs.1 The Space Age brought two new measurement types to bear: range and range-rate.2, 3 Initially these measurements were provided by radars through skin-tracking. Soon, however, transponders allowed better measurement accuracy and availability. To date, there are three types of IOD problems: angles only, range and angles, and range and range-rate. The Orion program exposed the need for an additional capability in the case of an extended loss of communication: when the measurements included antenna angles (azimuth and elevation or right-ascension and declination), range and range-rate. This paper contains the details of this new methodology. We will begin with state determination with range, range-rate and angle measurements. First the state estimate, given a sequence of these four time-synchronized measurements, will be analytically developed. Next the covariance matrix associated with this state estimate will be analytically developed. Finally, these will be verified by Monte Carlo simulations.

Christopher D’Souza↗
NASA NTRS1990

Experimental and predicted aerodynamic characteristics of a proposed Assured Crew Return Vehicle (ACRV) lifting-body configuration at Mach 6 and 10

The effects of Mach number, Reynolds number, and ratio of specific heats on the aerodynamic characteristics of a proposed Assured Crew Return Vehicle (ACRV) lifting-body configuration were examined for a range of angles of attack from -5 deg to 50 deg. Predictions made with a Langley-developed, three-dimensional Navier Stokes solver known as LAURA, which was exercised as an Euler solver for the present study, are compared with the experimental results. Unlike the Shuttle Orbiter, which experienced a significant nose-up increment in pitching moment with decreasing specific heat ratio (i.e., real gas effects), the aerodynamic characteristics of this lifting-body configuration are insensitive to changes in specific heat ratio. The maximum trimmed lift-to-drag ratio achieved was about 1.5. Predicted inviscid values of aerodynamic coefficients were generally in good agreement with measurement.

Micol, John R.↗
NASA NTRS1982

A natural conservative flux difference splitting for the hyperbolic systems of gasdynamics

Gasdynamic hyperbolic systems are treated by a novel conservative flux difference splitting upwind method which is applicable to explicit or implicit and iterative or direct schemes, for marching in time or space with the Euler or compressible Navier-Stokes equations. The method is able to capture sharp shocks correctly by maintaining global conservation. The embedded characteristics formulation is cast in the volumetric internal energy primitive variable, which is effective for real as well as perfect gases. Higher order upwind formulae are constructed from the distribution of pieces of the simple internal flux difference splitting, and the method reverts to first order at the appropriate points in order not to violate the domain of dependence by differencing across discontinuities. At insignificant computational overhead, switching is done by simple algebraic extensions of the truth functions that determine upwind direction for the first order scheme.

Lombard, C. K.↗
NASA NTRS1990

Near-field noise of a single-rotation propfan at an angle of attack

The near field noise characteristics of a propfan operating at an angle of attack are examined utilizing the unsteady pressure field obtained from a 3-D Euler simulation of the propfan flowfield. The near field noise is calculated employing three different procedures: a direct computation method in which the noise field is extracted directly from the Euler solution, and two acoustic-analogy-based frequency domain methods which utilize the computed unsteady pressure distribution on the propfan blades as the source term. The inflow angles considered are -0.4, 1.6, and 4.6 degrees. The results of the direct computation method and one of the frequency domain methods show qualitative agreement with measurements. They show that an increase in the inflow angle is accompanied by an increase in the sound pressure level at the outboard wing boom locations and a decrease in the sound pressure level at the (inboard) fuselage locations. The trends in the computed azimuthal directivities of the noise field also conform to the measured and expected results.

Nallasamy, M.↗
NASA NTRS2005

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗