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 703 records · Page 39

Approximate Solution Methods for Spectral Radiative Transfer in High Refractive Index Layers

Some ceramic materials for high temperature applications are partially transparent for radiative transfer. The refractive indices of these materials can be substantially greater than one which influences internal radiative emission and reflections. Heat transfer behavior of single and laminated layers has been obtained in the literature by numerical solutions of the radiative transfer equations coupled with heat conduction and heating at the boundaries by convection and radiation. Two-flux and diffusion methods are investigated here to obtain approximate solutions using a simpler formulation than required for exact numerical solutions. Isotropic scattering is included. The two-flux method for a single layer yields excellent results for gray and two band spectral calculations. The diffusion method yields a good approximation for spectral behavior in laminated multiple layers if the overall optical thickness is larger than about ten. A hybrid spectral model is developed using the two-flux method in the optically thin bands, and radiative diffusion in bands that are optically thick.

Siegel, R.↗

Micromechanical model of crack growth in fiber reinforced ceramics

A model based on the micromechanical mechanism of crack growth resistance in fiber reinforced ceramics is presented. The formulation of the model is based on a small scale geometry of a macrocrack with a bridging zone, the process zone, which governs the resistance mechanism. The effect of high toughness of the fibers in retardation of the crack advance, and the significance of the fiber pullout mechanism on the crack growth resistance, are reflected in this model. The model allows one to address issues such as influence of fiber spacing, fiber flexibility, and fiber matrix friction. Two approaches were used. One represents the fracture initiation and concentrated on the development of the first microcracks between fibers. An exact closed form solution was obtained for this case. The second case deals with the development of an array of microcracks between fibers forming the bridging zone. An implicit exact solution is formed for this case. In both cases, a discrete fiber distribution is incorporated into the solution.

Rubinstein, Asher A.↗

Micromechanical model of crack growth in fiber reinforced brittle materials

A model based on the micromechanical mechanism of crack growth resistance in fiber reinforced ceramics is presented. The formulation of the model is based on a small scale geometry of a macrocrack with a bridging zone, the process zone, which governs the resistance mechanism. The effect of high toughness of the fibers in retardation of the crack advance, and the significance of the fiber pullout mechanism on the crack growth resistance, are reflected in this model. The model allows one to address issues such as influence of fiber spacing, fiber flexibility, and fiber matrix friction. Two approaches were used. One represents the fracture initiation and concentrated on the development of the first microcracks between fibers. An exact closed form solution was obtained for this case. The second case deals with the development of an array of microcracks between fibers forming the bridging zone. An implicit exact solution is formed for this case. In both cases, a discrete fiber distribution is incorporated into the solution.

Rubinstein, Asher A.↗

Interior impedance wedge diffraction with surface waves

The exact impedance wedge solution is evaluated asymptotically using the method of steepest descents for plane wave illumination at normal incidence. Uniform but different impedances on each face are considered for both soft and hard polarizations. The asymptotic solution isolates the incident, singly reflected, multiply reflected, diffracted, and surface wave fields. Multiply reflected fields of any order are permitted. The multiply reflected fields from the exact solution are written as ratios of auxiliary Maliuzhinets functions, whereas a geometrical analysis gives the reflected fields as products of reflection coefficients. These two representations are shown to be identical in magnitude, phase and the angular range over which they exist. The diffracted field includes four Fresnel transition functions as in the perfect conductor case, and the expressions for the appropriate discontinuities at the shadow boundaries are presented. The surface wave exists over a finite angular range and only for certain surface impedances. A surface wave transition field is included to retain continuity. Computations are presented for interior wedge diffractions although the formulation is valid for both exterior and interior wedges.

Balanis, Constantine A.↗

Analytic solutions for single and multiple cylinders of gravitating polytropes in magnetostatic equilibrium

Exact analytic solutions for the static equilibrium of a gravitating plasma polytrope in the presence of magnetic fields are presented. The means of generating various equilibrium configurations to illustrate directly the complex physical relationships between pressure, magnetic fields, and gravity in self-gravitating systems is demonstrated. One of the solutions is used to model interstellar clouds suspended by magnetic fields against the galactic gravity such as may be formed by the Parker (1966) instability. It is concluded that the pinching effect of closed loops of magnetic fields in the clouds may be a dominant agent in further collapsing the clouds following their formation.

Lerche, I.↗

Light Scattering by Coated Sphere Immersed in Absorbing Medium: A Comparison between the FDTD and Analytic Solutions

A recently developed finite-difference time domain scheme is examined using the exact analytic solutions for light scattering by a coated sphere immersed in an absorbing medium. The relative differences are less than 1% in the extinction, scattering, and absorption efficiencies and less than 5% in the scattering phase functions. The definition of apparent single-scattering properties is also discussed. (C) 2003 Elsevier Ltd. All rights reserved.

Sun, W.↗

Singular perturbation analysis of the atmospheric orbital plane change problem

A three-state model is presented for the aeroassisted orbital plane change problem. A further model order reduction to a single state model is examined using singular perturbation theory. The optimal solution for this single state model compares favorably with the exact numerical solution using a four-state model; however, a separate boundary layer solution is required to satisfy the terminal constraint on altitude. This, in general, involves the solution of a two-point boundary value problem, but for a two-state model. An approximation is introduced to obtain an analytical control solution for lift and bank angle. Included are numerical simulation results of a guidance law derived from this analysis, along with comparison to earlier work by other researchers.

Calise, A. J.↗

On traveling waves in beams

The basic equations of Timoshenko for the motion of vibrating nonuniform beams, which allow for effects of transverse shear deformation and rotary inertia, are presented in several forms, including one in which the equations are written in the directions of the characteristics. The propagation of discontinuities in moment and shear, as governed by these equations, is discussed. Numerical traveling-wave solutions are obtained for some elementary problems of finite uniform beams for which the propagation velocities of bending and shear discontinuities are taken to be equal. These solutions are compared with modal solutions of Timoshenko's equations and, in some cases, with exact closed solutions. (author)

Leonard, Robert W↗

Exact, E = 0, classical and quantum solutions for general power-law oscillators

For zero energy, E = 0, we derive exact, classical and quantum solutions for all power-law oscillators with potentials V(r) = -gamma/r(exp nu), gamma greater than 0 and -infinity less than nu less than infinity. When the angular momentum is non-zero, these solutions lead to the classical orbits (p(t) = (cos mu(phi(t) - phi(sub 0)t))(exp 1/mu) with mu = nu/2 - 1 does not equal 0. For nu greater than 2, the orbits are bound and go through the origin. We calculate the periods and precessions of these bound orbits, and graph a number of specific examples. The unbound orbits are also discussed in detail. Quantum mechanically, this system is also exactly solvable. We find that when nu is greater than 2 the solutions are normalizable (bound), as in the classical case. Further, there are normalizable discrete, yet unbound, states. They correspond to unbound classical particles which reach infinity in a finite time. Finally, the number of space dimensions of the system can determine whether or not an E = 0 state is bound. These and other interesting comparisons to the classical system will be discussed.

Nieto, Michael Martin↗

Improved tangent-cone method for the Aerodynamic Preliminary Analysis System (APAS) version of the hypersonic arbitrary-body program

The Aerodynamic Preliminary Analysis System (APAS) utilizes a modified version of the Hypersonic Arbitrary-Body Program (HABP) Mark 3 code in its analysis rationale. Four methods are considered for incorporation into the code as the tangent-cone method. The combination of second-order slender body theory and the approximate solution of Hammitt and Murthy shows the best agreement with the exact numerical solutions and is thus included in the APAS production version of the HABP code.

Cruz, Christopher I.↗

Analytical description of the breakup of liquid jets in air

A viscous or inviscid cylindrical jet with surface tension in a vacuum tends to pinch due to the mechanism of capillary instability. Similarity solutions are constructed which describe this phenomenon as a critical time is encountered, for two physically distinct cases: inviscid jets governed by the Euler equations and highly viscous jets governed by the Stokes equations. In both cases the only assumption imposed is that at the time of pinching the jet shape has a radial length scale which is smaller than the axial length scale. For the inviscid case, we show that our solution corresponds exactly to one member of the one-parameter family of solutions obtained from slender jet theories and the shape of the jet is locally concave at breakup. For highly viscous jets our theory predicts local shapes which are monotonic increasing or decreasing indicating the formation of a mother drop connected to the jet by a thin fluid tube. This qualitative behavior is in complete agreement with both direct numerical simulations and experimental observations.

Papageorgiou, Demetrios T.↗

Benchmark problems and solutions

The scientific committee, after careful consideration, adopted six categories of benchmark problems for the workshop. These problems do not cover all the important computational issues relevant to Computational Aeroacoustics (CAA). The deciding factor to limit the number of categories to six was the amount of effort needed to solve these problems. For reference purpose, the benchmark problems are provided here. They are followed by the exact or approximate analytical solutions. At present, an exact solution for the Category 6 problem is not available.

Tam, Christopher K. W.↗

Relaxation of magnetotail plasmas

A quasi-thermodynamic model is presented for the relaxation of magnetotail plasmas during substorms, followed by quiet times. It is proposed that the plasma relaxes to a state of low-potential energy subject to a small number of global constraints. The constraints are exactly preserved by all ideal motions and, approximately, by a wide class of motions of the plasma undergoing magnetic reconnection. A variational principle which minimizes the free energy predicts the relaxed state. Exact, two-dimensional solutions of the relaxed state are obtained. A universal feature of the exact solutions is a chain of magnetic islands along the tail axis. Sufficient conditions for the stability of relaxed states are obtained from the second variation of the free-energy functional.

Bhattacharjee, A.↗

A fast numerical solution of scattering by a cylinder: Spectral method for the boundary integral equations

It is known that the exact analytic solutions of wave scattering by a circular cylinder, when they exist, are not in a closed form but in infinite series which converges slowly for high frequency waves. In this paper, we present a fast number solution for the scattering problem in which the boundary integral equations, reformulated from the Helmholtz equation, are solved using a Fourier spectral method. It is shown that the special geometry considered here allows the implementation of the spectral method to be simple and very efficient. The present method differs from previous approaches in that the singularities of the integral kernels are removed and dealt with accurately. The proposed method preserves the spectral accuracy and is shown to have an exponential rate of convergence. Aspects of efficient implementation using FFT are discussed. Moreover, the boundary integral equations of combined single and double-layer representation are used in the present paper. This ensures the uniqueness of the numerical solution for the scattering problem at all frequencies. Although a strongly singular kernel is encountered for the Neumann boundary conditions, we show that the hypersingularity can be handled easily in the spectral method. Numerical examples that demonstrate the validity of the method are also presented.

Hu, Fang Q.↗

A {3,2}-Order Bending Theory for Laminated Composite and Sandwich Beams

A higher-order bending theory is derived for laminated composite and sandwich beams thus extending the recent {1,2}-order theory to include third-order axial effect without introducing additional kinematic variables. The present theory is of order {3,2} and includes both transverse shear and transverse normal deformations. A closed-form solution to the cylindrical bending problem is derived and compared with the corresponding exact elasticity solution. The numerical comparisons are focused on the most challenging material systems and beam aspect ratios which include moderate-to-thick unsymmetric composite and sandwich laminates. Advantages and limitations of the theory are discussed.

Cook, Geoffrey M.↗

Semi-Analytic Reconstruction of Flux in Finite Volume Formulations

Semi-analytic reconstruction uses the analytic solution to a second-order, steady, ordinary differential equation (ODE) to simultaneously evaluate the convective and diffusive flux at all interfaces of a finite volume formulation. The second-order ODE is itself a linearized approximation to the governing first- and second- order partial differential equation conservation laws. Thus, semi-analytic reconstruction defines a family of formulations for finite volume interface fluxes using analytic solutions to approximating equations. Limiters are not applied in a conventional sense; rather, diffusivity is adjusted in the vicinity of changes in sign of eigenvalues in order to achieve a sufficiently small cell Reynolds number in the analytic formulation across critical points. Several approaches for application of semi-analytic reconstruction for the solution of one-dimensional scalar equations are introduced. Results are compared with exact analytic solutions to Burger s Equation as well as a conventional, upwind discretization using Roe s method. One approach, the end-point wave speed (EPWS) approximation, is further developed for more complex applications. One-dimensional vector equations are tested on a quasi one-dimensional nozzle application. The EPWS algorithm has a more compact difference stencil than Roe s algorithm but reconstruction time is approximately a factor of four larger than for Roe. Though both are second-order accurate schemes, Roe s method approaches a grid converged solution with fewer grid points. Reconstruction of flux in the context of multi-dimensional, vector conservation laws including effects of thermochemical nonequilibrium in the Navier-Stokes equations is developed.

Gnoffo, Peter A.↗

Rendezvous Guidance Technology

The literature on rendezvous guidance technology is reviewed and a bibliography of about 160 documents is presented. The importance of the proper selection of the ascent trajectory and of the orbit for a manned space station or an orbiting launch complex is emphasized. The way the launch delay problem affects the selection of the ascent trajectory and the target orbit is discussed. A specific ascent trajectory and an orbit for a manned space station or an orbiting launch complex is recommended on the premise that launch delay problems will continue to outweigh most of the other orbit selection criteria. A new method for determining the velocity corrections required for mid-course guidance is proposed . The proposed methods consist of an "exact-numerical" solution of the relative equations of motion in the "Shell-coordinate system." The proposed "exact-numerical" method extends the applicability of the mid- course equations to a much greater range than that of the currently available linear methods .

Trajectory↗

Determination of Nonlinear Stiffness Coefficients for Finite Element Models with Application to the Random Vibration Problem

In this paper, a method for obtaining nonlinear stiffness coefficients in modal coordinates for geometrically nonlinear finite-element models is developed. The method requires application of a finite-element program with a geometrically non- linear static capability. The MSC/NASTRAN code is employed for this purpose. The equations of motion of a MDOF system are formulated in modal coordinates. A set of linear eigenvectors is used to approximate the solution of the nonlinear problem. The random vibration problem of the MDOF nonlinear system is then considered. The solutions obtained by application of two different versions of a stochastic linearization technique are compared with linear and exact (analytical) solutions in terms of root-mean-square (RMS) displacements and strains for a beam structure.

Muravyov, Alexander A.↗