SEARCH · Engineering Papers

Results for “EXACT SOLUTIONS”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1975

Numerical techniques for the solution of the compressible Navier-Stokes equations and implementation of turbulence models

The time-splitting explicit numerical method of MacCormack is applied to separated turbulent boundary layer flow problems. Modifications of this basic method are developed to counter difficulties associated with complicated geometry and severe numerical resolution requirements of turbulence model equations. The accuracy of solutions is investigated by comparison with exact solutions for several simple cases. Procedures are developed for modifying the basic method to improve the accuracy. Numerical solutions of high-Reynolds-number separated flows over an airfoil and shock-separated flows over a flat plate are obtained. A simple mixing length model of turbulence is used for the transonic flow past an airfoil. A nonorthogonal mesh of arbitrary configuration facilitates the description of the flow field. For the simpler geometry associated with the flat plate, a rectangular mesh is used, and solutions are obtained based on a two-equation differential model of turbulence.

Baldwin, B. S.↗

NASA NTRS1985

Grid-free simulation of diffusion using random walk methods

The simulation of the diffusion of a continuum field by the random walk (RW) displacement of a set of particles is considered. Elements of the gradients of the diffusive concentration are transported by computational particles. It is demonstrated that, by the use of concentration gradients in the RW process, statistical errors are reduced and each realization of the numerical solution is a representation of the exact solution. The algorithm is grid-free, and the computational elements move to follow the gradients; hence, the algorithm is self-adaptive, and uniform resolution is achieved for all times.

Ghoniem, A. F.↗

NASA NTRS1989

Solutions of differential equations with regular coefficients by the methods of Richmond and Runge-Kutta

Numerical solutions of the differential equation which describe the electric field within an inhomogeneous layer of permittivity, upon which a perpendicularly-polarized plane wave is incident, are considered. Richmond's method and the Runge-Kutta method are compared for linear and exponential profiles of permittivities. These two approximate solutions are also compared with the exact solutions.

Cockrell, C. R.↗

NASA NTRS1991

Vibration suppression using a constrained rate-feedback-threshold control strategy

A finite time, minimum force rate-feedback-threshold controller is developed to bring a system with or without known external disturbances back into an 'allowable' state bound in finite time. The disturbances are assumed to be expandable in terms of Fourier series. The optimal control is defined by a two-point boundary value problem coupled to a set of definite integral constraints. Quasi-closed form solutions are derived which replace the solution of the two-point boundary value problem and definite integral constraints with the solution of algebraic equations and the calculation of matrix exponentials. Examples are provided which demonstrate the threshold control technique and compare the quasi-closed form solutions with numerical and MACSYMA generated exact solutions.

Zimmerman, D. C.↗

NASA NTRS1991

Multiple shooting algorithms for jump-discontinuous problems in optimal control and estimation

Multiple shooting algorithms are developed for jump-discontinuous two-point boundary value problems arising in optimal control and optimal estimation. Examples illustrating the origin of such problems are given to motivate the development of the solution algorithms. The algorithms convert the necessary conditions, consisting of differential equations and transversality conditions, into algebraic equations. The solution of the algebraic equations provides exact solutions for linear problems. The existence and uniqueness of the solution are proved.

Mook, D. J.↗

NASA NTRS1958

Application of the Method of Coordinate Perturbation to Unsteady Duct Flow

The method of coordinate perturbation is applied to the unsteady flow of a compressible fluid in ducts of variable cross section. Solutions, in the form of perturbation series, are obtained for unsteady flows in ducts for which the logarithmic derivative of area variation with respect to the space coordinate is a function of the 'smallness' parameter of the perturbation series. This technique is applied to the problem of the interaction of a disturbance and a shock wave in a diffuser flow. It is found that, for a special choice of the function describing the disturbance, the path of the shock wave can be expressed in closed form to first order. The method is then applied to the determination of the flow field behind a shock wave moving on a prescribed path in the x,t-plane. Perturbation series solutions for quite general paths are developed. The perturbation series solutions are compared with the more exact solutions obtained by the application of the method of characteristics. The approximate solutions are shown to be in reasonably accurate agreement with the solutions obtained by the method of characteristics.

Himmel, Seymour C.↗

NASA NTRSDate not supplied

Intent Modeling and Intent Conflict Probability Calculation for Operations in Upper Class E Airspace

This work presents probabilistic intent models for two typical vehicles in ETM operational environment. Several methods, including analytical methods for approximated solutions and a numerical method for exact solution, are presented and applied to compute the intent conflict probability for ETM operations. A comparison of these methods will be conducted in the final paper. Furthermore, simulations will be performed to verify the probabilistic intent model and verify the results of intent conflict probability from different methods.

Operational intent↗

NASA NTRS1995

An accuracy assessment of Cartesian-mesh approaches for the Euler equations

A critical assessment of the accuracy of Cartesian-mesh approaches for steady, transonic solutions of the Euler equations of gas dynamics is made. An exact solution of the Euler equations (Ringleb's flow) is used not only to infer the order of the truncation error of the Cartesian-mesh approaches, but also to compare the magnitude of the discrete error directly to that obtained with a structured mesh approach. Uniformly and adaptively refined solutions using a Cartesian-mesh approach are obtained and compared to each other and to uniformly refined structured mesh results. The effect of cell merging is investigated as well as the use of two different K-exact reconstruction procedures. The solution methodology of the schemes is explained and tabulated results are presented to compare the solution accuracies.

Coirier, William J.↗

NASA NTRS1980

Comparison of two leading uniform theories of edge diffraction with the exact uniform asymptotic solution

The diffraction of an arbitrary cylindrical wave by a half-plane has been treated by Rahmat-Samii and Mittra who used a spectral domain approach. In this paper, their exact solution for the total field is expressed in terms of a new integral representation. For large wave number k, two rigorous procedures are described for the exact uniform asymptotic expansion of the total field solution. The uniform expansions obtained are valid in the entire space, including transition regions around the shadow boundaries. The final results are compared with the formulations of two leading uniform theories of edge diffraction, namely, the uniform asymptotic theory and the uniform theory of diffraction. Some unique observations and conclusions are made in relating the two theories.

Boersma, J.↗

NASA NTRS1983

Buckling of imperfect periodic lattice structures

A simplified buckling analysis is presented for a family of periodic lattice structures such as those proposed for large space structures. A transcendental 6 x 6 matrix of eigenvalues is shown to be sufficient for modeling buckling behavior because member stiffnesses are based on an exact solution of the beam-column equation. Exact stiffnesses are derived for a curved member, thus allowing modeling of imperfect lattice structures. Comparisons of predictions of the lattice model with those available from shell and beam theory underscore the inaccuracies introduced by treating the lattice structure as a continuum. Sample calculations are provided for an isogrid cylinder and a three element double-laced truss.

Anderson, M. S.↗

NASA NTRS1985

Construction of approximate analytical solutions to a new class of non-linear oscillator equation

The principle of harmonic balance is invoked in the development of an approximate analytic model for a class of nonlinear oscillators typified by a mass attached to a stretched wire. By assuming that harmonic balance will hold, solutions are devised for a steady state limit cycle and/or limit point motion. A method of slowly varying amplitudes then allows derivation of approximate solutions by determining the form of the exact solutions and substituting into them the lowest order terms of their respective Fourier expansions. The latter technique is actually a generalization of the method proposed by Kryloff and Bogoliuboff (1943).

Mickens, R. E.↗

NASA NTRS1986

An efficient method for solving the steady Euler equations

An efficient numerical procedure for solving a set of nonlinear partial differential equations, the steady Euler equations, using Newton's linearization procedure is presented. A theorem indicating quadratic convergence for the case of differential equations is demonstrated. A condition for the domain of quadratic convergence Omega(2) is obtained which indicates that whether an approximation lies in Omega(2) depends on the rate of change and the smoothness of the flow vectors, and hence is problem-dependent. The choice of spatial differencing, of particular importance for the present method, is discussed. The treatment of boundary conditions is addressed, and the system of equations resulting from the foregoing analysis is summarized and solution strategies are discussed. The convergence of calculated solutions is demonstrated by comparing them with exact solutions to one and two-dimensional problems.

Liou, M.-S.↗

NASA NTRS1986

Complex and singular solutions of KdV and MKdV equations

The Korteweg-de Vries (KdV) and the modified Korteweg-de Vries (MKdV) equations are shown to have, besides the regular real solutions, exact regular complex as well as singular solutions. The singular solution for the KdV is real but for the MKdV it is pure imaginary. Implications of the complex solutions are discussed.

Buti, B.↗

NASA NTRS1998

Parallel 3D Multi-Stage Simulation of a Turbofan Engine

A 3D multistage simulation of each component of a modern GE Turbofan engine has been made. An axisymmetric view of this engine is presented in the document. This includes a fan, booster rig, high pressure compressor rig, high pressure turbine rig and a low pressure turbine rig. In the near future, all components will be run in a single calculation for a solution of 49 blade rows. The simulation exploits the use of parallel computations by using two levels of parallelism. Each blade row is run in parallel and each blade row grid is decomposed into several domains and run in parallel. 20 processors are used for the 4 blade row analysis. The average passage approach developed by John Adamczyk at NASA Lewis Research Center has been further developed and parallelized. This is APNASA Version A. It is a Navier-Stokes solver using a 4-stage explicit Runge-Kutta time marching scheme with variable time steps and residual smoothing for convergence acceleration. It has an implicit K-E turbulence model which uses an ADI solver to factor the matrix. Between 50 and 100 explicit time steps are solved before a blade row body force is calculated and exchanged with the other blade rows. This outer iteration has been coined a "flip." Efforts have been made to make the solver linearly scaleable with the number of blade rows. Enough flips are run (between 50 and 200) so the solution in the entire machine is not changing. The K-E equations are generally solved every other explicit time step. One of the key requirements in the development of the parallel code was to make the parallel solution exactly (bit for bit) match the serial solution. This has helped isolate many small parallel bugs and guarantee the parallelization was done correctly. The domain decomposition is done only in the axial direction since the number of points axially is much larger than the other two directions. This code uses MPI for message passing. The parallel speed up of the solver portion (no 1/0 or body force calculation) for a grid which has 227 points axially.

Turner, Mark G.↗

NASA NTRS2005

Multi-Dimensional Spectral Difference Method for Unstructured Grids

A new, high-order, conservative, and efficient method for conservation laws on unstructured grids is developed. It combines the best features of structured and unstructured grid methods to attain computational efficiency and geometric flexibility; it utilizes the concept of discontinuous and high-order local representations to achieve conservation and high accuracy; and it is based on the finite-difference formulation for simplicity. Universal reconstructions are obtained by distributing unknowns in a geometrically similar manner for all unstructured cells. Placements of the unknown and flux points with various order of accuracy are given for the line, triangular and tetrahedral elements. The data structure of the new method permits an optimum use of cache memory, resulting in further computational efficiency on modern computers. A new pointer system is developed that reduces memory requirements and simplifies programming for any order of accuracy. Numerical solutions are presented and compared with the exact solutions for wave propagation problems in both two and three dimensions to demonstrate the capability of the method. Excellent agreement has been found. The method is simpler and more efficient than previous discontinuous Galerkin and spectral volume methods for unstructured grids.

Liu, Yen↗

NASA NTRS1971

An exact closed form solution for constant area compressible flow with friction and heat transfer

The well-known differential equation for the one-dimensional flow of a compressible fluid with heat transfer and wall friction has no known solution in closed form for the general case. This report presents a closed form solution for the special case of constant heat flux per unit length and constant specific heat. The solution was obtained by choosing the square of a dimensionless flow parameter as one of the independent variables to describe the flow. From this exact solution, an approximate simplified form is derived that is applicable for predicting subsonic flow performance characteristics for many types of constant area passages in internal flow. The data included in this report are considered sufficiently accurate for use as a guide in analyzing and designing internal gas flow systems.

Sturas, J. I.↗

NASA NTRS1974

A comparison of numerical solutions to the inviscid equations of fluid motion

The second-order MacCormack method and the third-order Rusanov and Kutler-Warming-Lomax methods are applied to the inviscid 1-d Burgers' equation, wedge flow, and the problem of shock reflection from a rigid boundary. The numerical solution in each case is compared to the exact solution and the quantitative estimates of accuracy are obtained. Results of this study show that the third-order Kutler-Warming-Lomax or the Rusanov methods tuned for minimum dissipation or minimum dispersion provide the most accurate solution in each of the examples considered.

Anderson, D. A.↗

NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗