Self-similar Imploding Solutions of the 1D Compressible Euler Equations with a Far Field Cutoff
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Engineering Papers
Results for “Euler flow”
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.
Pressure measurements in subcritical pure superfluid helium flow.
Pressure measurements in pure superfluid helium 2 flow at subcritical velocities, showing quantitative agreement with Euler equation with zero curl
Corke, H. E.↗
Implicit finite-difference procedures for the computation of vortex wakes
Implicit finite-difference procedures for the primitive form of the incompressible Navier-Stokes and the compressible Euler equations are used to compute vortex wake flows. The partial differential equations in strong conservation-law form are transformed to cluster grid points in regions with large changes in vorticity. In addition to clustering, fourth-order accurate, spatial difference operators are used to help resolve the flow-field gradients. The use of implicit time-differencing permits large time steps to be taken since temporal variations are typically small. Computational efficiency is achieved by approximate factorization. Both two-dimensional and preliminary three-dimensional calculations are described and qualitatively compared with existing experimental data.
Steger, J. L.↗
Multiple-grid acceleration of Lax-Wendroff algorithms
A technique for accelerating the convergence of a one-step Lax-Wendroff method to a steady-state solution is discussed and its applicability extended to the more general class of two-step Lax-Wendroff methods. Several two-step methods which lead to quite efficient multiple grid algorithms are discussed. Computational results are presented using the full two dimensional Euler equations for both subcritical and shocked supercritical flows. Extensions and generalizations are mentioned.
Johnson, G. M.↗
Multiple-gridding of the Euler equations with an implicit scheme
The multiple-grid scheme of Ni (1981) for the solution of the unsteady Euler equations for quasi-one-dimensional transonic flow problems is analyzed according to its ability to accelerate convergence to a steady state solution, its applicability to an implicit scheme, its accuracy, and its stability limits. This scheme is applied successfully to McCormack's (1981) implicit method of solving quasi-one-dimensional transonic flow problems for CFL numbers significantly greater than unity. It is determined that the second order accuracy of this implicit method is preserved by employing the multiple gridding. However, stability limitations that are flow dependent are introduced into the implicit method by the use of the multiple-grid scheme which will curtail its use under certain flow conditions. A scaling modification for the multiple grid scheme is developed which further improves its performance when used with explicit and artificially-damped implicit methods without any increase in computational work.
Stubbs, R. M.↗
Recent developments in multigrid methods for the steady Euler equations
The solution by multigrid techniques of the steady inviscid compressible equations of gas dynamics, the Euler equations is investigated. Steady two dimensional transonic flow over an airfoil section is studied intensively. Most of the material is applicable to three dimensional flow problems of aerodynamic interest.
Jespersen, D. C.↗
Parametric effects of CFL number and artificial smoothing on numerical solutions using implicit approximate factorization algorithm
An implicit approximate factorization algorithm is employed to quantify the parametric effects of Courant number and artificial smoothing on numerical solutions of the unsteady 3-D Euler equations for a windmilling propeller (low speed) flow field. The results show that propeller global or performance chracteristics vary strongly with Courant number and artificial dissipation parameters, though the variation is such less severe at high Courant numbers. Candidate sets of Courant number and dissipation parameters could result in parameter-dependent solutions. Parameter-independent numerical solutions can be obtained if low values of the dissipation parameter-time step ratio are used in the computations. Furthermore, it is realized that too much artificial damping can degrade numerical stability. Finally, it is demonstrated that highly resolved meshes may, in some cases, delay convergence, thereby suggesting some optimum cell size for a given flow solution. It is suspected that improper boundary treatment may account for the cell size constraint.
Daso, E. O.↗
Coupled aerodynamic and acoustical predictions for turboprops
To predict the noise fields for proposed turboprop airplanes, an existing turboprop noise code by Farassat has been modified to accept blade pressure inputs from a three-dimensional aerodynamic code. A Euler-type code can handle the nonlinear transonic flow of these high-speed, highly swept blades. This turbofan code was modified to allow the calculation mesh to extend to about twice the blade radius and to apply circumferential periodicity rather than solid-wall boundary conditions on the blade in the region between the blade tip and the outer shroud. Outputs were added for input to the noise prediction program and for color contour plots of various flow variables. The Farassat input subroutines were modified to read files of blade coordinates and predicted surface pressures. Aerodynamic and acoustic results are shown for the SR-3 model blade. Comparison of the acoustic predicted results with measured data show good agreement.
Clark, Bruce J.↗
A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach
A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.
Koenig, Herbert A.↗
Adaptive computational methods for aerothermal heating analysis
The development of adaptive gridding techniques for finite-element analysis of fluid dynamics equations is described. The developmental work was done with the Euler equations with concentration on shock and inviscid flow field capturing. Ultimately this methodology is to be applied to a viscous analysis for the purpose of predicting accurate aerothermal loads on complex shapes subjected to high speed flow environments. The development of local error estimate strategies as a basis for refinement strategies is discussed, as well as the refinement strategies themselves. The application of the strategies to triangular elements and a finite-element flux-corrected-transport numerical scheme are presented. The implementation of these strategies in the GIM/PAGE code for 2-D and 3-D applications is documented and demonstrated.
Price, John M.↗
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↗
A genuinely multi-dimensional upwind cell-vertex scheme for the Euler equations
The solution of the two-dimensional Euler equations is based on the two-dimensional linear convection equation and the Euler-equation decomposition developed by Hirsch et al. The scheme is genuinely two-dimensional. At each iteration, the data are locally decomposed into four variables, allowing convection in appropriate directions. This is done via a cell-vertex scheme with a downwind-weighted distribution step. The scheme is conservative, and third-order accurate in space. The derivation and stability analysis of the scheme for the convection equation, and the derivation of the extension to the Euler equations are given. Preconditioning techniques based on local values of the convection speeds are discussed. The scheme for the Euler equations is applied to two channel-flow problems. It is shown to converge rapidly to a solution that agrees well with that of a third-order upwind solver.
Powell, Kenneth G.↗
A genuinely multi-dimensional upwind cell-vertex scheme for the Euler equations
A scheme of solving the two-dimensional Euler equations is developed. The scheme is genuinely two-dimensional. At each iteration, the data are locally decomposed into four variables, allowing convection in appropriate directions. This is done via a cell-vertex scheme with a downwind-weighted distribution step. The scheme is conservative and third-order accurate in space. The derivation and stability analysis of the scheme for the convection equation, and the derivation of the extension to the Euler equations are given. Preconditioning techniques based on local values of the convection speeds are discussed. The scheme for the Euler equations is applied to two channel-flow problems. It is shown to converge rapidly to a solution that agrees well with that of a third-order upwind solver.
Powell, Kenneth G.↗
Unsteady blade-surface pressures on a large-scale advanced propeller: Prediction and data
An unsteady 3-D Euler analysis technique is employed to compute the flow field of an advanced propeller operating at an angle of attack. The predicted blade pressure waveforms are compared with wind tunnel data at two Mach numbers, 0.5 and 0.2. The inflow angle is three degrees. For an inflow Mach number of 0.5, the predicted pressure response is in fair agreement with data: the predicted phases of the waveforms are in close agreement with data while the magnitudes are underpredicted. At the low Mach number of 0.2 (takeoff), the numerical solution shows the formation of a leading edge vortex which is in qualitative agreement with measurements. However, the highly nonlinear pressure response measured on the blade suction surface is not captured in the present inviscid analysis.
Nallasamy, M.↗
Vorticity dynamics of inviscid shear layers
The inviscid evolution of a two-dimensional shear layer is simulated numerically by a scheme based on a kinematic decomposition of the unsteady flow. Lagrangian and Weber transformations of the incompressible Euler equations result in a Clebsch representation that separates the flowfield into rotational and irrotational components. These transformations produce the initial construction of the flowfield and define its subsequent evolution.
Yokota, Jeffrey W.↗
Computational analysis of methods for reduction of induced drag
The purpose of this effort was to perform a computational flow analysis of a design concept centered around induced drag reduction and tip-vortex energy recovery. The flow model solves the unsteady three-dimensional Euler equations, discretized as a finite-volume method, utilizing a high-resolution approximate Riemann solver for cell interface flux definitions. The numerical scheme is an approximately-factored block LU implicit Newton iterative-refinement method. Multiblock domain decomposition is used to partition the field into an ordered arrangement of blocks. Three configurations are analyzed: a baseline fuselage-wing, a fuselage-wing-nacelle, and a fuselage-wing-nacelle-propfan. Aerodynamic force coefficients, propfan performance coefficients, and flowfield maps are used to qualitatively access design efficacy. Where appropriate, comparisons are made with available experimental data.
Janus, J. M.↗
The block adaptive multigrid method applied to the solution of the Euler equations
In the present study, a scheme capable of solving very fast and robust complex nonlinear systems of equations is presented. The Block Adaptive Multigrid (BAM) solution method offers multigrid acceleration and adaptive grid refinement based on the prediction of the solution error. The proposed solution method was used with an implicit upwind Euler solver for the solution of complex transonic flows around airfoils. Very fast results were obtained (18-fold acceleration of the solution) using one fourth of the volumes of a global grid with the same solution accuracy for two test cases.
Pantelelis, Nikos↗
Euler Technology Assessment for Preliminary Aircraft Design: Compressibility Predictions by Employing the Cartesian Unstructured Grid SPLITFLOW Code
The objective of the second phase of the Euler Technology Assessment program was to evaluate the ability of Euler computational fluid dynamics codes to predict compressible flow effects over a generic fighter wind tunnel model. This portion of the study was conducted by Lockheed Martin Tactical Aircraft Systems, using an in-house Cartesian-grid code called SPLITFLOW. The Cartesian grid technique offers several advantages, including ease of volume grid generation and reduced number of cells compared to other grid schemes. SPLITFLOW also includes grid adaption of the volume grid during the solution to resolve high-gradient regions. The SPLITFLOW code predictions of configuration forces and moments are shown to be adequate for preliminary design, including predictions of sideslip effects and the effects of geometry variations at low and high angles-of-attack. The transonic pressure prediction capabilities of SPLITFLOW are shown to be improved over subsonic comparisons. The time required to generate the results from initial surface data is on the order of several hours, including grid generation, which is compatible with the needs of the design environment.
Finley, Dennis B.↗