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 577 records · Page 32

Prediction of pressure drop in fluid tuned mounts using analytical and computational techniques

A simplified model for predicting pressure drop in fluid tuned isolator mounts was developed. The model is based on an exact solution to the Navier-Stokes equations and was made more general through the use of empirical coefficients. The values of these coefficients were determined by numerical simulation of the flow using the commercial computational fluid dynamics (CFD) package FIDAP.

Lasher, William C.↗

A mathematical model for simulating noise suppression of lined ejectors

A mathematical model containing the essential features embodied in the noise suppression of lined ejectors is presented. Although some simplification of the physics is necessary to render the model mathematically tractable, the current model is the most versatile and technologically advanced at the current time. A system of linearized equations and the boundary conditions governing the sound field are derived starting from the equations of fluid dynamics. A nonreflecting boundary condition is developed. In view of the complex nature of the equations, a parametric study requires the use of numerical techniques and modern computers. A finite element algorithm that solves the differential equations coupled with the boundary condition is then introduced. The numerical method results in a matrix equation with several hundred thousand degrees of freedom that is solved efficiently on a supercomputer. The model is validated by comparing results either with exact solutions or with approximate solutions from other works. In each case, excellent correlations are obtained. The usefulness of the model as an optimization tool and the importance of variable impedance liners as a mechanism for achieving broadband suppression within a lined ejector are demonstrated.

Watson, Willie R.↗

A flux splitting scheme with high-resolution and robustness for discontinuities

A flux splitting scheme is proposed for the general nonequilibrium flow equations with an aim at removing numerical dissipation of Van-Leer-type flux-vector splittings on a contact discontinuity. The scheme obtained is also recognized as an improved Advection Upwind Splitting Method (AUSM) where a slight numerical overshoot immediately behind the shock is eliminated. The proposed scheme has favorable properties: high-resolution for contact discontinuities; conservation of enthalpy for steady flows; numerical efficiency; applicability to chemically reacting flows. In fact, for a single contact discontinuity, even if it is moving, this scheme gives the numerical flux of the exact solution of the Riemann problem. Various numerical experiments including that of a thermo-chemical nonequilibrium flow were performed, which indicate no oscillation and robustness of the scheme for shock/expansion waves. A cure for carbuncle phenomenon is discussed as well.

Wada, Yasuhiro↗

Detection of Delaminations in Composite Beams Using Piezoelectric Sensors

This paper investigates the feasibility of a proposed technique for detecting delamination using piezoelectric layers or patches embedded or bonded to a composite structure. Variations in the voltage generated in the piezoelectric layers indicates the presence and location of delamination, when the structure is excited either externally or via actuators. The theoretical foundations of a method for predicting the dynamic response of delaminated composite beams with piezoelectric layers are described. The governing equations are presented for the case of external vibroacoustic excitation, as well as, for the case of locally induced vibrations by some of the embedded piezoelectric elements. An exact solution is developed within the limits of linear laminate theory. Applications illustrate the feasibility of delamination detection in cantilever beams. The results illustrate that the proposed technique may provide accurate detection of the presence, size, and location of a delamination.

Saravanos, Dimitris A.↗

Numerical results in the statistical theory of thermal turbulence

Several statistical theories of turbulence are applied to thermal convection problem, and are compared in a series of numerical calculations. Results are presented for both the case of very large and very small Prandtl numbers, and for moderate Rayleigh numbers. The theories compared are the direct-interaction approximation, the self-consistent field approximation, the quasi-linear approximation, the quasi-normal approximation, and a 'Markovianized' quasi-normal approximation. Results for these theories are then compared with some simple exact solutions to the statistical problem. The results indicate that the direct-interaction and the self-consistent field method give satisfactory results at all Rayleigh numbers investigated, whereas the procedures based on the quasi-normal approximation give sensible results only at small Rayleigh numbers.

Herring, J. R.↗

Calibration experiments with DC-8 SAR data

This paper is a summary of some research activities in calibrating DC-8 SAR data conducted in 1989/90 by the author and others. Firstly, by examining the process of 'symmetrization' in the production of DC-8 Stokes' matrix data, it is shown that an exact solution to the problem of calibrating this data exists, given certain assumptions. Secondly a new, more general, system model for radar polarimeters is proposed, together with a simple approach for verifying the model's validity. The importance of the orientation angle about the radar line of sight for some types of calibration devices is discussed and some results of an experiment to cross-calibrate the C-Band VV data from the DC-8 and DLR E-SAR sensors are presented.

Freeman, Anthony↗

Nonlinear stability of oscillatory core-annular flow: A generalized Kuramoto-Sivashinsky equation with time periodic coefficients

In this paper the nonlinear stability of two-phase core-annular flow in a pipe is examined when the acting pressure gradient is modulated by time harmonic oscillations and viscosity stratification and interfacial tension is present. An exact solution of the Navier-Stokes equations is used as the background state to develop an asymptotic theory valid for thin annular layers, which leads to a novel nonlinear evolution describing the spatio-temporal evolution of the interface. The evolution equation is an extension of the equation found for constant pressure gradients and generalizes the Kuramoto-Sivashinsky equation with dispersive effects found by Papageorgiou, Maldarelli & Rumschitzki, Phys. Fluids A 2(3), 1990, pp. 340-352, to a similar system with time periodic coefficients. The distinct regimes of slow and moderate flow are considered and the corresponding evolution is derived. Certain solutions are described analytically in the neighborhood of the first bifurcation point by use of multiple scales asymptotics. Extensive numerical experiments, using dynamical systems ideas, are carried out in order to evaluate the effect of the oscillatory pressure gradient on the solutions in the presence of a constant pressure gradient.

Coward, Adrian V.↗

Random vibration of nonlinear beams by the new stochastic linearization technique

In this paper, the beam under general time dependent stationary random excitation is investigated, when exact solution is unavailable. Numerical simulations are carried out to compare its results with those yielded by the conventional linearization techniques. It is found that the modified version of the stochastic linearization technique yields considerably more accurate results for the mean square displacement of the beam than the conventional equivalent linearization technique, especially in the case of large nonlinearity.

Fang, J.↗

Evolution of wave function in a dissipative system

For a dissipative system with Ohmic friction, we obtain a simple and exact solution for the wave function of the system plus the bath. It is described by the direct product in two independent Hilbert space. One of them is described by an effective Hamiltonian, the other represents the effect of the bath, i.e., the Brownian motion, thus clarifying the structure of the wave function of the system whose energy is dissipated by its interaction with the bath. No path integral technology is needed in this treatment. The derivation of the Weisskopf-Wigner line width theory follows easily.

Yu, Li-Hua↗

Discrete-Layer Piezoelectric Plate and Shell Models for Active Tip-Clearance Control

The objectives of this work were to develop computational tools for the analysis of active-sensory composite structures with added or embedded piezoelectric layers. The targeted application for this class of smart composite laminates and the analytical development is the accomplishment of active tip-clearance control in turbomachinery components. Two distinct theories and analytical models were developed and explored under this contract: (1) a discrete-layer plate theory and corresponding computational models, and (2) a three dimensional general discrete-layer element generated in curvilinear coordinates for modeling laminated composite piezoelectric shells. Both models were developed from the complete electromechanical constitutive relations of piezoelectric materials, and incorporate both displacements and potentials as state variables. This report describes the development and results of these models. The discrete-layer theories imply that the displacement field and electrostatic potential through-the-thickness of the laminate are described over an individual layer rather than as a smeared function over the thickness of the entire plate or shell thickness. This is especially crucial for composites with embedded piezoelectric layers, as the actuating and sensing elements within these layers are poorly represented by effective or smeared properties. Linear Lagrange interpolation polynomials were used to describe the through-thickness laminate behavior. Both analytic and finite element approximations were used in the plane or surface of the structure. In this context, theoretical developments are presented for the discrete-layer plate theory, the discrete-layer shell theory, and the formulation of an exact solution for simply-supported piezoelectric plates. Finally, evaluations and results from a number of separate examples are presented for the static and dynamic analysis of the plate geometry. Comparisons between the different approaches are provided when possible, and initial conclusions regarding the accuracy and limitations of these models are given.

Heyliger, P. R.↗

Harmonic oscillator representation in the theory of scattering and nuclear reactions

The following questions, concerning the application of the harmonic oscillator representation (HOR) in the theory of scattering and reactions, are discussed: the formulation of the scattering theory in HOR; exact solutions of the free motion Schroedinger equation in HOR; separable expansion of the short range potentials and the calculation of the phase shifts; 'isolated states' as generalization of the Wigner-von Neumann bound states embedded in continuum; a nuclear coupled channel problem in HOR; and the description of true three body scattering in HOR. As an illustration the soft dipole mode in the (11)Li nucleus is considered in a frame of the (9)Li+n+n cluster model taking into account of three body continuum effects.

Smirnov, Yuri F.↗

Effect of mesh distortion on the accuracy of transverse shear stresses and their sensitivity coefficients in multilayered composites

A study is made of the effect of mesh distortion on the accuracy of transverse shear stresses and their first-order and second-order sensitivity coefficients in multilayered composite panels subjected to mechanical and thermal loads. The panels are discretized by using a two-field degenerate solid element, with the fundamental unknowns consisting of both displacement and strain components, and the displacement components having a linear variation throughout the thickness of the laminate. A two-step computational procedure is used for evaluating the transverse shear stresses. In the first step, the in-plane stresses in the different layers are calculated at the numerical quadrature points for each element. In the second step, the transverse shear stresses are evaluated by using piecewise integration, in the thickness direction, of the three-dimensional equilibrium equations. The same procedure is used for evaluating the sensitivity coefficients of transverse shear stresses. Numerical results are presented showing no noticeable degradation in the accuracy of the in-plane stresses and their sensitivity coefficients with mesh distortion. However, such degradation is observed for the transverse shear stresses and their sensitivity coefficients. The standard of comparison is taken to be the exact solution of the three-dimensional thermoelasticity equations of the panel.

Noor, Ahmed K.↗

Acoustic scattering from ellipses by the modal element method

The modal element method is used to study acoustic scattering from ellipses, which may be acoustically soft (absorbing) or hard (reflecting). Because exact solutions are available, the results provide a benchmark for algorithm performance for scattering from airfoils and similar shapes. Numerical results for scattering from rigid ellipses are presented for a wide variety of eccentricities at moderate frequencies. These results indicate that the method is practical.

Kreider, Kevin L.↗

A comparative study of upwind and MacCormack schemes for CAA benchmark problems

In this study, upwind schemes and MacCormack schemes are evaluated as to their suitability for aeroacoustic applications. The governing equations are cast in a curvilinear coordinate system and discretized using finite volume concepts. A flux splitting procedure is used for the upwind schemes, where the signals crossing the cell faces are grouped into two categories: signals that bring information from outside into the cell, and signals that leave the cell. These signals may be computed in several ways, with the desired spatial and temporal accuracy achieved by choosing appropriate interpolating polynomials. The classical MacCormack schemes employed here are fourth order accurate in time and space. Results for categories 1, 4, and 6 of the workshop's benchmark problems are presented. Comparisons are also made with the exact solutions, where available. The main conclusions of this study are finally presented.

Viswanathan, K.↗

FORTRAN program for analyzing ground-based radar data: Usage and derivations, version 6.2

A postflight FORTRAN program called 'radar' reads and analyzes ground-based radar data. The output includes position, velocity, and acceleration parameters. Air data parameters are also provided if atmospheric characteristics are input. This program can read data from any radar in three formats. Geocentric Cartesian position can also be used as input, which may be from an inertial navigation or Global Positioning System. Options include spike removal, data filtering, and atmospheric refraction corrections. Atmospheric refraction can be corrected using the quick White Sands method or the gradient refraction method, which allows accurate analysis of very low elevation angle and long-range data. Refraction properties are extrapolated from surface conditions, or a measured profile may be input. Velocity is determined by differentiating position. Accelerations are determined by differentiating velocity. This paper describes the algorithms used, gives the operational details, and discusses the limitations and errors of the program. Appendices A through E contain the derivations for these algorithms. These derivations include an improvement in speed to the exact solution for geodetic altitude, an improved algorithm over earlier versions for determining scale height, a truncation algorithm for speeding up the gradient refraction method, and a refinement of the coefficients used in the White Sands method for Edwards AFB, California. Appendix G contains the nomenclature.

Haering, Edward A., Jr.↗

Finite-surface method for the Maxwell equations with corner singularities

The finite-surface method for the two-dimensional Maxwell equations in generalized coordinates is extended to treat perfect conductor boundaries with sharp corners. Known singular forms of the grid and the electromagnetic fields in the neighborhood of each corner are used to obtain accurate approximations to the surface and line integrals appearing in the method. Numerical results are presented for a harmonic plane wave incident on a finite flat plate. Comparisons with exact solutions show good agreement.

Vinokur, Marcel↗

Thin shell formation in radiative shocks. 1: Supernova remnants in low-density media

This paper explores the onset of thin-shell formation in interstellar shocks associated with supernova explosions. We outline a simple but useful scheme that indicates the time at which thin shell formation begins for supernova remnants (SNRs) evolving in a range of interstellar environments, extending the previous analytical models to arbitrary power-law density media. The result depends on the gas cooling properties and the shock velocity and radius. This is then applied to the specific case of SNRs in low-density media. The procedure for defining the time for the onset of shell formation, t(sub sf), equates the value of the adiabat, kappa = p/rho(exp gamma), to zero using the known time dependence of the shock radius and velocity. For the case of a power-law density ambient medium of the form rho(r) = Br(exp -omega), it is found that shell formation can be prevented when the ambient density drops faster than a critical rate. For a cooling function of the form Lambda = Lambda(sub 0) tau(exp beta), with beta = -0.5 (appropriate for line cooling), shell formation never occurs for omega greater than or equal to 9/5. The shell formation time is then computed for spherical shocks in a power-law density medium. For omega = 0, the onset of shell formation is found to be at t(sub sf) approx. equal to 2.87 x 10(exp 4) E51(exp 3/14) n(sub 0 exp -4.7) yr, which agrees well with previous estimates derived by other means. We compare the analytical shell formation time with the results of detailed numerical models for omega = 0 and three different ambient densities and find good agreement. The extension of the criterion for the onset of thin shell formation using the ratio of cooling to swept-up column density is also described. This method provides a useful approximation for cases when the exact solution is not known.

Franco, Jose↗

Nonreflective boundary conditions for high-order methods

A different approach to nonreflective boundary conditions for the Euler equations is presented. This work is motivated by a need for inflow and outflow boundary conditions that do not limit the useful accuracy of high-order accurate methods. The primary interest is in the propagation and convection of continuous acoustic and convective waves. This new approach employs the exact solution to finite waves to relate interior values and ambient conditions to boundary values. The method is first presented in one dimension and then generalized to multidimensions. Grid refinement studies are used to demonstrate high-order convergence for both one-dimensional and two-dimensional flows.

Atkins, H.↗