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 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Frequency response in short thermocouple wires

Theoretical expressions are derived for the steady state frequency response of a thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for a nonuniform wire with unequal material properties and wire diameters across the junction. The amplitude ratio at low frequency omega approaches 0 agrees with the results of Scadron and Warshawsky (1952) for a steady state temperature distribution. Moreover, the frequency response for a nonuniform wire in the limit of infinite length l approaches infinity is shown to reduce to a simple expression that is analogous to the classic first order solution for a thermocouple wire with uniform properties. Theoretical expressions are also derived for the steady state frequency response of a supported thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for both a supported one material wire and a two material wire with unequal material properties across the junction. For the case of a one material supported wire, an exact solution is derived which compares favorably with an approximate expression that only matches temperatures at the support junction. Moreover, for the case of a two material supported wire, an analytical expression is derived that closely correlates numerical results. Experimental measurements are made for the steady state frequency response of a supported thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for both a supported one material wire (type K) and a two material wire (type T) with unequal material properties across the junction. The data for the amplitude ratio and phase angle are correlated to within 10 pct. with the theoretical predictions of Forney and Fralick (1991). This is accomplished by choosing a natural frequency omega sub n for the wire data to correlate the first order response at large gas temperature frequencies. It is found that a large bead size, however, will increase the amplitude ratio at low frequencies but decrease the natural frequency of the wire. The phase angle data are also distorted for imperfect junctions.

Forney, L. J.↗

NASA NTRS2016

Hyperbolic Method for Dispersive PDEs: Same High-Order of Accuracy for Solution, Gradient, and Hessian

In this paper, we introduce a new hyperbolic first-order system for general dispersive partial differential equations (PDEs). We then extend the proposed system to general advection-diffusion-dispersion PDEs. We apply the fourth-order RD scheme of Ref. 1 to the proposed hyperbolic system, and solve time-dependent dispersive equations, including the classical two-soliton KdV and a dispersive shock case. We demonstrate that the predicted results, including the gradient and Hessian (second derivative), are in a very good agreement with the exact solutions. We then show that the RD scheme applied to the proposed system accurately captures dispersive shocks without numerical oscillations. We also verify that the solution, gradient and Hessian are predicted with equal order of accuracy.

Mazaheri, Alireza↗

NASA NTRS2004

2004: Finite-Difference Time Domain Solution of Light Scattering by an Infinite Dielectric Column Immersed in an Absorbing Medium

The two-dimensional (2-D) finite-difference time domain (FDTD) method is applied to calculate light scattering and absorption by an arbitrarily shaped infinite column embedded in an absorbing dielectric medium. A uniaxial perfectly matched layer (UPML) absorbing boundary condition (ABC) is used to truncate the computational domain. The single-scattering properties of the infinite column embedded in the absorbing medium, including scattering phase functions, extinction and absorption efficiencies, are derived using an area integration of the internal field. An exact solution for light scattering and absorption by a circular cylinder in an absorbing medium is used to examine the accuracy of the 2-D UPML FDTD code. With use of a cell size of 1/120 incident wavelength in the FDTD calculations, the errors in the extinction and absorption efficiencies and asymmetry factors from the 2-D UPML FDTD are generally smaller than approx .1%. The errors in the scattering phase functions are typically smaller than approx .4%. Using the 2-D UPML FDTD technique, light scattering and absorption by long noncircular columns embedded in absorbing media can be accurately solved.

Sun, W.↗

NASA NTRS1953

Deflection and Stress Analysis of Thin Solid Wings of Arbitrary Plan Form with Particular Reference to Delta Wings

The structural analysis of arbitrary solid cantilever wings by small-deflection thin-plate theory is reduced to the solution of linear ordinary differential equations by the assumption that the chordwise deflections at any spanwise station may be expressed in the form of a power series in which the coefficients are functions of the spanwise coordinate. If the series is limited to the first two and three terms (that is, if linear and parabolic chordwise deflections, respectively, are assumed), the differential equations for the coefficients are solved exactly for uniformly loaded solid delta wings of constant thickness ratio. For cases for which exact solutions to the differential equations cannot be obtained, a numerical procedure is derived. Experimental deflection and stress data for constant-thickness delta-plate specimens of 45 degree and 60 degree sweep are presented and are found to compare favorably with the present theory.

Stein, Manuel↗

NASA NTRS2011

Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods

A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.

Atkins, Harold L.↗

NASA NTRS2011

Simulation of Benchmark Cases with the Terminal Area Simulation System (TASS)

The hydrodynamic core of the Terminal Area Simulation System (TASS) is evaluated against different benchmark cases. In the absence of closed form solutions for the equations governing atmospheric flows, the models are usually evaluated against idealized test cases. Over the years, various authors have suggested a suite of these idealized cases which have become standards for testing and evaluating the dynamics and thermodynamics of atmospheric flow models. In this paper, simulations of three such cases are described. In addition, the TASS model is evaluated against a test case that uses an exact solution of the Navier-Stokes equations. The TASS results are compared against previously reported simulations of these benchmark cases in the literature. It is demonstrated that the TASS model is highly accurate, stable and robust.

Ahmad, Nashat N.↗

NASA NTRS2011

Simulation of Benchmark Cases with the Terminal Area Simulation System (TASS)

The hydrodynamic core of the Terminal Area Simulation System (TASS) is evaluated against different benchmark cases. In the absence of closed form solutions for the equations governing atmospheric flows, the models are usually evaluated against idealized test cases. Over the years, various authors have suggested a suite of these idealized cases which have become standards for testing and evaluating the dynamics and thermodynamics of atmospheric flow models. In this paper, simulations of three such cases are described. In addition, the TASS model is evaluated against a test case that uses an exact solution of the Navier-Stokes equations. The TASS results are compared against previously reported simulations of these banchmark cases in the literature. It is demonstrated that the TASS model is highly accurate, stable and robust.

Ahmad, Nash'at↗

NASA NTRS1994

Stochastic solution to quantum dynamics

The quantum Liouville equation in the Wigner representation is solved numerically by using Monte Carlo methods. For incremental time steps, the propagation is implemented as a classical evolution in phase space modified by a quantum correction. The correction, which is a momentum jump function, is simulated in the quasi-classical approximation via a stochastic process. The technique, which is developed and validated in two- and three- dimensional momentum space, extends an earlier one-dimensional work. Also, by developing a new algorithm, the application to bound state motion in an anharmonic quartic potential shows better agreement with exact solutions in two-dimensional phase space.

John, Sarah↗

NASA NTRS1995

Evaluation of numerical schemes for the analysis of sound generation by blade-gust interaction

In this investigation three different numerical algorithms have been utilized to compute the flow about a flat plate in the presence of a transverse gust described by a sinusoidal disturbance. The three schemes include the MacCormack explicit finite difference scheme which is second order accurate in both time and space, the Gottlieb and Turkel modification of MacCormack's scheme which is fourth order accurate in space and second order accurate in time, (referred to as the 2-4 scheme), and a two step scheme developed by Bayliss et al. which has second order temporal accuracy and sixth order spatial accuracy (a 2-6 scheme). The flow field results are obtained with these schemes by using the same code with the only difference being the implementation of the respective solution algorithms. The problem is set up so that the sinusoidal disturbance is imposed at the surface of the flat plate as a surface boundary condition. Thus the problem is treated as scattering problem. The computed results include the time average of the acoustic pressure squared along grid lines five points away from the boundaries; distribution throughout the computational domain is monitored at various times. The numerical results are compared with an exact solution obtained by Atassi, Dusey, and Davis.

Scott, J. N.↗

NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗

NASA NTRS1983

Time domain response envelope for structural dynamic systems

A transient envelope solution method is developed for payload structural systems. This method requires that the external forcing functions are decomposed into a quasi-static part and a dynamic part. The amplitude of the peak envelope and the corresponding varying frequencies for the generalized forcing functions are required for obtaining the response envelopes. The proposed method is applied to a sample problem as well as a spacecraft structural system and the results are in good agreement with the exact solution.

Chen, J.-C.↗

NASA NTRS1983

Ultrasonic attenuation of a void-containing medium for very long wavelengths

Ultrasonic longitudinal through-thickness attenuation in an isotropic medium due to scattering by randomly distributed voids is considered analytically. The attenuation is evaluated on the assumption of no interaction between voids. The scattered power is assumed to be entirely lost, thus accounting for the ultrasonic attenuation. The scattered power due to the presence of a void is described in terms of the scattering cross section of the void. An exact solution exists for the scattering cross section of a spherical void. An approximate solution for the scattering cross section of an ellipsoidal void is developed based on the so-called Born approximation commonly used in quantum mechanics. This approximate solution is valid for k sub p a sub i 1, where k sub p is the wave number of the incident longitudinal wave and a sub i is the largest dimension of the void. It is found that the shape of the void has negligible effect on the scattering cross section and that only the volume of the void is important. Thus, it is noted that in cases where k sup p a sub i 1, the exact scattering cross section of a spherical void having the same volume as an arbitrarily shaped void can be used for evaluating ultrasonic attenuation.

Williams, J. H., Jr.↗

NASA NTRS1986

Numerical solutions of the complete Navier-Stokes equations

Using ideas from the kinetic theory, the Navier-Stokes equations are modified in such a way that they can be cast as a set of first order hyperbolic equations. This is achieved by incorporating time dependent terms into the definition of the stress tensor and the heat flux vectors. The boundary conditions are then determined from the theory of characteristics. Because the resulting equations reduce to the traditional Navier-Stokes equations when the steady state is reached, the present approach provides a straightforward scheme for the determination of inflow and outflow boundary conditions. The method is validated by comparing its predictions with known exact solutions of the steady Navier-Stokes equations.

Hassan, H. A.↗

NASA NTRS1975

Second-order small-disturbance solutions for hypersonic flow over power-law bodies

Similarity solutions were found which give the adiabatic flow of an ideal gas about two-dimensional and axisymmetric power-law bodies at infinite Mach number to second order in the body slenderness parameter. The flow variables were expressed as a sum of zero-order and perturbation similarity functions for which the axial variations in the flow equations separated out. The resulting similarity equations were integrated numerically. The solutions, which are universal functions, are presented in graphic and tabular form. To avoid a singularity in the calculations, the results are limited to body power-law exponents greater than about 0.85 for the two-dimensional case and 0.75 for the axisymmetric case. Because of the entropy layer induced by the nose bluntness (for power-law bodies other than cones and wedges), only the pressure function is valid at the body surface. The similarity results give excellent agreement with the exact solutions for inviscid flow over wedges and cones having half-angles up to about 20 deg. They give good agreement with experimental shock-wave shapes and surface-pressure distributions for 3/4-power axisymmetric bodies, considering that Mach number and boundary-layer displacement effects are not included in the theory.

Townsend, J. C.↗

NASA NTRS1983

Fast Euler solver for steady, 1-dimensional flows

A numerical technique to solve the Euler equations for steady, one dimensional flows is presented. The technique is essentially implicit, but is structured as a sequence of explicit solutions for each Riemann variable separately. Each solution is obtained by integrating in the direction prescribed by the propagation of the Riemann variables. The technique is second-order accurate. It requires very few steps for convergence, and each step requires a minimal number of operations. Therefore, it is three orders of magnitude more efficient than a standard time-dependent technique. The technique works very well for transonic flows and provides shock fitting with errors as small as 0.001. Results are presented for subsonic problems. Errors are evaluated by comparison with exact solutions.

Moretti, G.↗

NASA NTRS1984

Ultrasonic attenuation of a void-containing medium for very long wavelengths

Ultrasonic longitudinal through-thickness attenuation in an isotropic medium due to scattering by randomly distributed voids is considered analytically. The attenuation is evaluated on the assumption of no interaction between voids. The scattered power is assumed to be entirely lost, thus accounting for the ultrasonic attenuation. The scattered power due to the presence of a void is described in terms of the scattering cross section of the void. An exact solution exists for the scattering cross section of a spherical void. An approximate solution for the scattering cross section of an ellipsoidal void is developed based on the so-called Born approximation commonly used in quantum mechanics. This approximate solution is valid for k sub p a sub i much less than 1, where k sub p is the wave number of the incident longitudinal wave and a sub i is the largest dimension of the void. It is found that the shape of the void has negligible effect on the scattering cross section and that only the volume of the void is important. Thus, it is noted that in cases where k sub p a sub i is much less than 1, the exact scattering cross section of a spherical void having the same volume as an arbitrarily shaped void can be used for evaluating ultrasonic attenuation. Previously announced in STAR as N83-28466

Williams, J. H., Jr.↗

NASA NTRS1998

The Space-Time Conservation Element and Solution Element Method-A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: Numerical Simulation of Shock Waves and Contact Discontinuities - 2

Without resorting to special treatment for each individual test case, the 1D and 2D CE/SE shock-capturing schemes described previously (in Part I) are used to simulate flows involving phenomena such as shock waves, contact discontinuities, expansion waves and their interactions. Five 1D and six 2D problems are considered to examine the capability and robustness of these schemes. Despite their simple logical structures and low computational cost (for the 2D CE/SE shock-capturing scheme, the CPU time is about 2 micro-secs per mesh point per marching step on a Cray C90 machine), the numerical results, when compared with experimental data, exact solutions or numerical solutions by other methods, indicate that these schemes can accurately resolve shock and contact discontinuities consistently.

Wang, Xiao-Yen↗

NASA NTRS1993

H(sub infinity)-type filter for spacecraft attitude estimation

A nonlinear filtering theory from a deterministic point of view is presented and an application to attitude determination is considered. The approach that is taken in this paper is motivated largely by the H(sun infinity) control and estimation theory for linear systems which has been evolved within the last decade. Rather than formulating the estimation problem as a game played by two adversaries, as has been done in the linear case, we employ in this work some notions from the theory of dissipative systems as a vehicle for arriving at certain Hamilton-Jacobi inequality, which in turn, provides a solution to the filtering problem, whenever it is satisfied. Application of this method to a linear estimation problem and to the problem of estimating a spacecraft attitude quaternion and gyro drift bias vector are presented. In limiting cases, these give the Kalman filter and the extended Kalman filter, respectively. The main advantages of this approach over its probabilistic counterpart are that this approach does not require a prior knowledge of any statistics, and that in general it is more amenable to a quantitative assessment regarding approximations such as linearization, and that in certain cases this approach yields an exact solution to the nonlinear filtering.

Markley, F. Landis↗