Engineering Papers⌕ Search

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.

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

At least 1,045 records · Page 58

A quartic dispersion equation for internal gravity waves in the thermosphere

A new quartic dispersion equation in the square of the complex vertical wave number is derived by employing the 'shallow atmosphere' approximation and an ion drag approximation. These approximations allow the coefficients of the quartic equation to be given in terms of the corresponding cubic equation (which neglects the Coriolis force and the zonal ion drag component), but modified to take into account these neglected effects. Coupling between the extraordinary viscosity wave mode and the other three wave modes is highlighted and numerical solutions are compared for this quartic equation, an exact eighth order equation and the cubic equation. For the first time the validity of using the 'shallow atmosphere' approximation to describe internal gravity wave motions is demonstrated.

Hickey, M. P.↗

Study of three-dimensional effects on vortex breakdown

The incompressible axisymmetric steady Navier-Stokes equations in primitive variables are used to simulate vortex breakdown. The equations, discretized using a second-order, central-difference scheme, are linearized and then solved using an exact LU decomposition, Gaussian elimination, and Newton iteration. Solutions are presented for Reynolds numbers, based on vortex-core radius, as high as 1500. An attempt to study the stability of the axisymmetric solutions against three-dimensional perturbations is discussed.

Salas, M. D.↗

Vortex breakdown simulation - A circumspect study of the steady, laminar, axisymmetric model

The incompressible axisymmetric steady Navier-Stokes equations are written using the streamfunction-vorticity formulation. The resulting equations are discretized using a second-order central-difference scheme. The discretized equations are linearized and then solved using an exact LU decomposition, Gaussian elimination, and Newton iteration. Solutions are presented for Reynolds numbers (based on vortex core radius) 100-1800 and swirl parameter 0.9-1.1. The effects of inflow boundary conditions, the location of farfield and outflow boundaries, and mesh refinement are examined. Finally, the stability of the steady solutions is investigated by solving the time-dependent equations.

Salas, M. D.↗

Thin-plate spline quadrature of geodetic integrals

Thin-plate spline functions (known for their flexibility and fidelity in representing experimental data) are especially well-suited for the numerical integration of geodetic integrals in the area where the integration is most sensitive to the data, i.e., in the immediate vicinity of the evaluation point. Spline quadrature rules are derived for the contribution of a circular innermost zone to Stoke's formula, to the formulae of Vening Meinesz, and to the recursively evaluated operator L(n) in the analytical continuation solution of Molodensky's problem. These rules are exact for interpolating thin-plate splines. In cases where the integration data are distributed irregularly, a system of linear equations needs to be solved for the quadrature coefficients. Formulae are given for the terms appearing in these equations. In case the data are regularly distributed, the coefficients may be determined once-and-for-all. Examples are given of some fixed-point rules. With such rules successive evaluation, within a circular disk, of the terms in Molodensky's series becomes relatively easy. The spline quadrature technique presented complements other techniques such as ring integration for intermediate integration zones.

Vangysen, Herman↗

A numerical method for the solution of three dimensional, incompressible, viscous flows over slender bodies

A marching iterative method for the solution of the three dimensional, incompressibhle, steady and parabolized Navier-Stokes equations is described. The equations are written in primitive variables and discretized in general axisymmetric orthogonal coordinate systems. The coupled set of finite-difference equations are solved without any splitting or factorization errors. Moreover, the continuity equation and the two crossflow momentum equations are exactly satisfied at every step of the iterative process. The solution scheme is equivalent to the solution of one Poisson equation by the Successive Plane Over Relaxation method and has good convergence properties. Other existing solution methods resemble a Jacobi-type iterative scheme and therefore are less efficient. Numerical experiments include the laminar, incompressible flow over prolate spheroids at incidence.

Rosenfeld, Moshe↗

A solution to the problem of SAR range curvature

When synthetic aperture radar systems are pushed to attain finer resolution at larger ranges than was previously the case for remote sensing purposes, the geometric signal aberration known as range curvature arises. Known techniques for correcting range curvature are exact at only one selected range, thus forcing neighboring ranges to use the same correction as an approximation. A solution to the problem is proposed that is exact at all ranges, thus simplifying and improving the image processing for such systems.

Raney, R. K.↗

Hypersonic flows generated by parabolic and paraboloidal shock waves

A computer algorithm has been developed to determine the blunt-body flowfields supporting symmetric parabolic and paraboloidal shock waves at infinite free-stream Mach number. Solutions are expressed in an analytic form as high-order power series, in the coordinate normal to the shock, whose coefficients can be determined exactly. Analytic continuation is provided by the use of Pade approximations. Test cases provide solutions of very high accuracy. In the axisymmetric case for gamma equals 715 the solution has been found far downstream, where it agrees with the modified blast-wave results. For plane flow, on the other hand, a limit line appears within the shock layer, a short distance past the sonic line, suggesting the presence of an imbedded shock. Local solutions in the downstream limit are discussed.

Schwartz, L. W.↗

Time domain modelling and stability analysis of an integral pulse frequency modulated dc to dc power converter

Using state variable representation a nonlinear, discrete-time system is derived that models the converter exactly. This system is linearized about its steady state solution, and converter stability, transient response and audio susceptibility are studied. The steady state solution of the converter is stable if and only if all the roots of the linearized system are absolutely less than unity. Excellent agreement with laboratory test data has been observed.

Iwens, R. P.↗

Oscillatory conductive heat transfer for a fiber in an ideal gas

A description of the thermal effects created by placing a cylindrical fiber in an inviscid, ideal gas, through which an acoustic wave propagates, is presented. The fibers and the gas have finite heat capacities and thermal conductivities. Expressions for the temperature distribution in the gas and in the material are determined. The temperature distribution is caused by pressure oscillations in the gas which, in turn, are caused by the passage of an acoustic wave. The relative value of a dimensionless parameter is found to be indicative of whether the exact or approximate equations should be used in the solution. This parameter is a function of the thermal conductivities and heat capacities of the fiber and gas, the acoustic frequency, and the fiber diameter.

Kuntz, H. L.↗

A Comparison of Filter-based Approaches for Model-based Prognostics

Model-based prognostics approaches use domain knowledge about a system and its failure modes through the use of physics-based models. Model-based prognosis is generally divided into two sequential problems: a joint state-parameter estimation problem, in which, using the model, the health of a system or component is determined based on the observations; and a prediction problem, in which, using the model, the stateparameter distribution is simulated forward in time to compute end of life and remaining useful life. The first problem is typically solved through the use of a state observer, or filter. The choice of filter depends on the assumptions that may be made about the system, and on the desired algorithm performance. In this paper, we review three separate filters for the solution to the first problem: the Daum filter, an exact nonlinear filter; the unscented Kalman filter, which approximates nonlinearities through the use of a deterministic sampling method known as the unscented transform; and the particle filter, which approximates the state distribution using a finite set of discrete, weighted samples, called particles. Using a centrifugal pump as a case study, we conduct a number of simulation-based experiments investigating the performance of the different algorithms as applied to prognostics.

prognostics↗

Numerical integration of nearly-Hamiltonian systems

The reported investigation is concerned with the solution of systems of differential equations which are derived from a Hamiltonian function in the extended phase space. The problem selected involves a one-dimensional perturbed harmonic oscillator. The van der Pol equation considered has an exact asymptotic value for its amplitude. Comparisons are made between a numerical solution and a known analytical solution. In addition to the van der Pol problem, known solutions regarding the restricted problem of three bodies are used as examples for perturbed Keplerian motion. The extended phase space Hamiltonian discussed by Stiefel and Scheifele (1971) is considered. A description is presented of two canonical formulations of the perturbed harmonic oscillator.

Bond, V. R.↗

Calculation of the transient motion of elastic airfoils forced by control surface motion and gusts

The time-domain equations of motion of elastic airfoil sections forced by control surface motions and gusts were developed for the case of incompressible flow. Extensive use was made of special functions related to the inverse transform of Theodorsen's function. Approximations for the special cases of zero stream velocity, small time, large and time are given. A numerical solution technique for the solution of the general case is given. Examples of the exact transient response of an airfoil are presented.

Balakrishnan, A. V.↗

Buckling loads for stiffened panels subjected to combined longitudinal compression and shear loadings: Results obtained with PASCO, EAL, and STAGS computer

The shear buckling analyses used in PASCO are summarized. The PASCO analyses include the basic VIPASA analysis, which is essentially exact for longitudinal and transverse loads, and a smeared orthotropic solution which was added to alleviate a shortcoming in the VIPASA analysis. Buckling results are presented for six stiffened panels loaded by combinations of longitudinal compression and shear. The buckling results were obtained with the PASCO, EAL, and STAGS computer programs. The EAL and STAGS solutions were obtained with a fine finite element mesh and provide calculations for the entire range of combinations of longitudinal compression and shear loadings.

Stroud, W. J.↗

Acoustic transmission in non-uniform ducts with mean flow. I - The method of weighted residuals. II - The finite element method

The problem of acoustic transmission through nonuniform ducts containing a high-speed subsonic flow is studied by means of the method of weighted residuals in the form of a modified Galerkin method and a Galerkin formulation of the finite element method. The method of weighted residuals is shown to employ the basis functions generated from eigenvalue calculations for the case of no flow, and is verified by comparison with exact eigenvalue calculations in the uniform duct case and numerical solutions of the one-dimensional form of the equations in the nonuniform duct case. The finite element scheme based on both the Galerkin method and the residual least squares method and employing eight-noded isoparametric elements is presented and used to investigate multimodal propagation by the coupling of the solution in the duct nonuniform section to modal expansions in uniform sections. Comparison of the results of the two methods reveals them to be in substantial agreement, and predicts the importance of multimodal interactions at high Mach numbers.

Eversman, W.↗

A kinetic cross-field streaming instability

In a high-beta plasma the so-called modified-two-stream instability, which results from strongly magnetized electrons drifting relative to unmagnetized ions across a homogeneous magnetic field, is misnamed because the mode is highly kinetic, particularly when the relative streaming velocity exceeds the Alfven speed of the plasma. This kinetic cross-field streaming instability is investigated in detail, examining the effect of the electromagnetic terms and the stability boundaries in both low- and high-beta plasmas. An approximate dispersion relation showing the relation of this mode to the whistler is derived and solutions of it are compared with those obtained from the exact dispersion relation. The kinetic mode, unlike the usual modified-two-stream instability, is not stabilized by electromagnetic effects when the relative electron-ion drift speed exceeds the Alfven speed.

Wu, C. S.↗

Buckling loads of stiffened panels subjected to combined longitudinal compression and shear: Results obtained with PASCO, EAL, and STAGS computer programs

Buckling analyses used in PASCO are summarized with emphasis placed on the shear buckling analyses. The PASCO buckling analyses include the basic VIPASA analysis, which is essentially exact for longitudinal and transverse loads, and a smeared stiffener solution, which treats a stiffened panel as an orthotropic plate. Buckling results are then presented for seven stiffened panels loaded by combinations of longitudinal compression and shear. The buckling results were obtained with the PASCO, EAL, and STAGS computer programs. The EAL and STAGS solutions were obtained with a fine finite element mesh and are very accurate. These finite element solutions together with the PASCO results for pure longitudinal compression provide benchmark calculations to evaluate other analysis procedures.

Stroud, W. J.↗