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 829 records · Page 46

On the Convergence of Numerical Computations for Both Exact and Approximate Solutions for Electromagnetic Scattering by Nonspherical Dielectric Particles

We summarize the size parameter range of the applicability of four lightscattering computational methods for nonspheric dielectric particles. These methods include two exact methods - the extended boundary condition method (EBCM) and the invariant imbedding T-matrix method (II-TM) and two approximate approaches - the physical-geometric optics method (PGOM) and the improved geometric optics method (IGOM). For spheroids, the single-scattering properties computed by EBCM and II-TM agree for size parameters up to 150, and the comparison gives us confidence in using IITM as a benchmark for size parameters up to 150 for other geometries (e.g., hexagonal columns) because the applicability of II-TM with respect to particle shape is generic, as demonstrated in our previous studies involving a complex aggregate. This study demonstrates the convergence of the exact II-TM and approximate PGOM solutions for the complete set of single-scattering properties of a nonspherical shape other than spheroids and circular cylinders with particle sizes of ~48λ, specifically a hexagonal column with a size parameter of length as kL=300, where k=2π/λ and L is the column length. IGOM is also quite accurate except near the exact 180º backscattering direction. This study demonstrates that a synergetic combination of the numerically-exact II-TM and the approximate PGOM can seamlessly cover the entire size parameter range of practical interest. To demonstrate the applicability of the approach, we compute the optical properties of dust particles with a downstream application to the retrieval of dust aerosol optical thickness and effective particle size from polarimetric observations.

Ping Yang↗

Polarized Bidirectional Reflectance of Optically Thick Sparse Particulate Layers: an Efficient Numerically Exact Radiative-Transfer Solution

We describe a simple yet efficient numerical algorithm for computing polarized bidirectional reflectance of an optically thick (semi-infinite), macroscopically flat layer composed of statistically isotropic and mirror symmetric random particles. The spatial distribution of the particles is assumed to be sparse, random, and statistically uniform. The 44 Stokes reflection matrix is calculated by iterating the Ambartsumian's vector nonlinear integral equation. The result is a numerically exact solution of the vector radiative transfer equation and as such fully satisfies the energy conservation law and the fundamental reciprocity relation. Since this technique bypasses the computation of the internal radiation field, it is very fast and highly accurate. The FORTRAN implementation of the technique is publicly available on the World Wide Web at http://www.giss.nasa.gov/staff/ mmishchenko/brf. It can be combined with several existing computer programs providing the requisite single-scattering properties of spherical or morphologically complex particles and applied to a wide range of optical characterization problems. Benchmark results obtained with this program can be used for testing alternative solvers of the vector radiative transfer equation.

radiative transfer↗

Solution of the exact equations for three-dimensional atmospheric entry using directly matched asymptotic expansions

The problem of determining the trajectories, partially or wholly contained in the atmosphere of a spherical, nonrotating planet, is considered. The exact equations of motion for three-dimensional, aerodynamically affected flight are derived. Modified Chapman variables are introduced and the equations are transformed into a set suitable for analytic integration using asymptotic expansions. The trajectory is solved in two regions: the outer region, where the force may be considered a gravitational field with aerodynamic perturbations, and the inner region, where the force is predominantly aerodynamic, with gravity as a perturbation. The two solutions are matched directly. A composite solution, valid everywhere, is constructed by additive composition. This approach of directly matched asymptotic expansions applied to the exact equations of motion couched in terms of modified Chapman variables yields an analytical solution which should prove to be a powerful tool for aerodynamic orbit calculations.

Busemann, A.↗

The role of multiple scattering in one-dimensional radiative transfer

The usual methods of solving the radiative transfer equation yield answers which embrace all orders of scattering and thus shed little light on the underlying physical process. The present analysis examines the contributions of the various orders of scattering to the one dimensional transfer of radiation. In the one dimensional case an exact analytical solution exists and the problem reduces to that of expanding these exact solutions in powers of the albedo for single scattering. Formulas are given which permit the calculation of any order of scattering in an atmosphere of arbitrary optical thickness, particle albedo, and asymmetry parameter. The results should aid in identifying those physical situations where only the lowest orders of scattering play a significant role and where appropriate approximate methods might provide results of acceptable accuracy.

Adamson, D.↗

A Method for Computing Turbulent Heat Transfer in the Presence of a Streamwise Pressure Gradient for Bodies in High-Speed Flow

A modified Reynolds analogy between skin friction and heat transfer which depends upon local pressure gradient is derived. Exact and approximate solutions are derived from the differential equations; the exact solution is applicable for arbitrary initial (transition) conditions and the approximate solution requires fully developed turbulent flow from stagnation point or leading edge. The exact solution (restricted to stagnation initial conditions) and the approximate solutions are shown to agree with one another within 5 percent when applied to several blunt shapes. The present solutions generally predict the measured heating rates on these bodies within the accuracy of the measurements provided transition began upstream of the peak heating region. The present solutions appear to be sufficiently accurate for design purposes.

Cohen, Nathaniel B.↗

The extended auxiliary equation mapping method to determine novel exact solitary wave solutions of the nonlinear fractional PDEs

Abstract In this paper, some new nonlinear fractional partial differential equations (PDEs) have been considered.Three models are including the space-time fractional-order Boussinesq equation, space-time (2 + 1)-dimensional breaking soliton equations, and space-time fractional-order SRLW equation describe the behavior of these equations in the diverse applications. Meanwhile, the fractional derivatives in the sense of β -derivative are defined. Some fractional PDEs will convert to the considered ordinary differential equations by the help of transformation of β -derivative. These equations are analyzed utilizing an integration scheme, namely, the extended auxiliary equation mapping method. The different kinds of traveling wave solutions, solitary, topological, dark soliton, periodic, kink, and rational, fall out as a by-product of this scheme. Finally, the existence of the solutions for the constraint conditions is also shown. The outcome indicates that some fractional PDEs are used as a growing finding in the engineering sciences, mathematical physics, and so forth.

Engineering↗

Some examples of exact and approximate solutions in small particle scattering - A progress report

The formulation of basic equations from which the scattering of radiation by a localized variation in a medium is discussed. These equations are developed in both the differential and the integral form. Primary interest is in the scattering of electromagnetic waves for which the solution of the vector wave equation with appropriate boundary conditions must be considered. Scalar scattering by an infinite homogeneous isotropic circular cylinder, and scattering of electromagnetic waves by infinite circular cylinders are treated, and the case of the finite circular cylinder is considered. A procedure is given for obtaining angular scattering distributions from spheroids.

Greenberg, J. M.↗

Minimization theory of induced drag subject to constraint conditions

Exact analytical solutions in terms of induced drag influence coefficients can be attained which define the spanwise loading with minimized induced drag, subject to specified constraint conditions, for any nonplanar wing shape or number of lift plus wing bending moment about a given wing span station. Example applications of the theory are made to a biplane, a wing in ground effect, a cruciform wing, a V-wing, a planar-wing winglet, and linked wingtips in formation flying. For minimal induced drag, the spanwise loading, relative to elliptic, is outboard for the biplane and is inboard for the wing in ground effect and for the planar-wing winglet. A spinoff of the triplane solution provides mathematically exact equations for downwash and sidewash about a planar vorticity sheet having an arbitrary loading distribution.

Deyoung, J.↗

Validation Tests of TASS for Application to 3-D Vortex Simulations

Direct analytical solutions can be useful in validating the Core formulation of numerical systems. In this document an exact analytical solution to the nonlinear Navier-Stokes equation is compared to the numerical results from the three-dimensional Terminal Area Simulation System (TASS). This exact solution, of which the derivation is included, is for Beltrami type flow. Direct comparison of TASS to the analytical Beltrami solution is then used in evaluating the accuracy of TASS.

Switzer, George F.↗

Numerical solution of laminar jet mixing with and without free stream.

Systematic numerical solutions of two-dimensional and axisymmetrical laminar jet of an incompressible fluid with and without free stream have been obtained. For the case without free stream, the exact numerical solutions have been obtained for different initial velocity profiles and compared with experimental results and similarity solutions. For the case with free stream, the numerical solutions have been compared with the linearized analytical solutions. An approximate numerical solution for the axial velocity distribution of the three-dimensional laminar jet is proposed. The accuracy of the method has been determined from the exact solutions of the two limiting cases of three-dimensional jets, i.e., the two-dimensional and axisymmetric cases.

Pai, S. I.↗

The Principle of Energetic Consistency: Application to the Shallow-Water Equations

If the complete state of the earth's atmosphere (e.g., pressure, temperature, winds and humidity, everywhere throughout the atmosphere) were known at any particular initial time, then solving the equations that govern the dynamical behavior of the atmosphere would give the complete state at all subsequent times. Part of the difficulty of weather prediction is that the governing equations can only be solved approximately, which is what weather prediction models do. But weather forecasts would still be far from perfect even if the equations could be solved exactly, because the atmospheric state is not and cannot be known completely at any initial forecast time. Rather, the initial state for a weather forecast can only be estimated from incomplete observations taken near the initial time, through a process known as data assimilation. Weather prediction models carry out their computations on a grid of points covering the earth's atmosphere. The formulation of these models is guided by a mathematical convergence theory which guarantees that, given the exact initial state, the model solution approaches the exact solution of the governing equations as the computational grid is made more fine. For the data assimilation process, however, there does not yet exist a convergence theory. This book chapter represents an effort to begin establishing a convergence theory for data assimilation methods. The main result, which is called the principle of energetic consistency, provides a necessary condition that a convergent method must satisfy. Current methods violate this principle, as shown in earlier work of the author, and therefore are not convergent. The principle is illustrated by showing how to apply it as a simple test of convergence for proposed methods.

Cohn, Stephen E.↗

Improved Data Reduction Algorithm for the Needle Probe Method Applied to In-Situ Thermal Conductivity Measurements of Lunar and Planetary Regoliths

The needle probe method (also known as the' hot wire' or 'line heat source' method) is widely used for in-situ thermal conductivity measurements on soils and marine sediments on the earth. Variants of this method have also been used (or planned) for measuring regolith on the surfaces of extra-terrestrial bodies (e.g., the Moon, Mars, and comets). In the near-vacuum condition on the lunar and planetary surfaces, the measurement method used on the earth cannot be simply duplicated, because thermal conductivity of the regolith can be approximately 2 orders of magnitude lower. In addition, the planetary probes have much greater diameters, due to engineering requirements associated with the robotic deployment on extra-terrestrial bodies. All of these factors contribute to the planetary probes requiring much longer time of measurement, several tens of (if not over a hundred) hours, while a conventional terrestrial needle probe needs only 1 to 2 minutes. The long measurement time complicates the surface operation logistics of the lander. It also negatively affects accuracy of the thermal conductivity measurement, because the cumulative heat loss along the probe is no longer negligible. The present study improves the data reduction algorithm of the needle probe method by shortening the measurement time on planetary surfaces by an order of magnitude. The main difference between the new scheme and the conventional one is that the former uses the exact mathematical solution to the thermal model on which the needle probe measurement theory is based, while the latter uses an approximate solution that is valid only for large times. The present study demonstrates the benefit of the new data reduction technique by applying it to data from a series of needle probe experiments carried out in a vacuum chamber on JSC-1A lunar regolith stimulant. The use of the exact solution has some disadvantage, however, in requiring three additional parameters, but two of them (the diameter and the volumetric heat capacity of the probe) can be measured and the other (the volumetric heat capacity of the regolith/stimulant) may be estimated from the surface geologic observation and temperature measurements. Therefore, overall, the new data reduction scheme would make in-situ thermal conductivity measurement more practical on planetary missions.

JSC-1A↗

Crack-mouth displacements for semielliptical surface cracks subjected to remote tension and bending loads

The exact analytical solution for an embedded elliptical crack in an infinite body subjected to arbitrary loading was used in conjunction with the finite element alternating method to obtain crack-mouth-opening displacements (CMOD) for surface cracks in finite plates subjected to remote tension. Identical surface-crack configurations were also analyzed with the finite element method using 20-noded element for plates subjected to both remote tension and bending. The CMODs from these two methods generally agreed within a few percent of each other. Comparisons made with experimental results obtained from surface cracks in welded aluminum alloy specimens subjected to tension also showed good agreement. Empirical equations were developed for CMOD for a wide range of surface-crack shapes and sizes subjected to tension and bending loads. These equations were obtained by modifying the Green-Sneddon exact solution for an elliptical crack in an infinite body to account for finite boundary effects. These equations should be useful in monitoring surface-crack growth in tests and in developing complete crack-face-displacement equations for use in three-dimensional weight-function methods.

Raju, Ivatury S.↗

An Analysis of Exact and Approximate Equations for the Temperature Distribution in an Insulated Thick Skin Subjected to Aerodynamic Heating

The problem of calculating the temperature distribution in an insulated slab is investigated. Exact and approximate solutions are obtained, and the results are compared to determine the ranges of applicability of the approximations. The approximations are found to be within 5 percent of the exact solution when the ratio of the thermal capacitance of the metal to that of the insulation and the ratio of the conductance of the metal to that of the insulation are sufficiently large. The roots of the characteristic equation of the exact solution are generally applicable to the two-slab heat-transfer problem and are tabulated up to the first nine roots.

Harris, Robert S., Jr.↗

On a class of unsteady three-dimensional Navier Stokes solutions relevant to rotating disc flows: Threshold amplitudes and finite time singularities

A class of exact steady and unsteady solutions of the Navier Stokes equations in cylindrical polar coordinates is given. The flows correspond to the motion induced by an infinite disc rotating with constant angular velocity about the z-axis in a fluid occupying a semi-infinite region which, at large distances from the disc, has velocity field proportional to (x,-y,O) with respect to a Cartesian coordinate system. It is shown that when the rate of rotation is large, Karman's exact solution for a disc rotating in an otherwise motionless fluid is recovered. In the limit of zero rotation rate a particular form of Howarth's exact solution for three-dimensional stagnation point flow is obtained. The unsteady form of the partial differential system describing this class of flow may be generalized to time-periodic equilibrium flows. In addition the unsteady equations are shown to describe a strongly nonlinear instability of Karman's rotating disc flow. It is shown that sufficiently large perturbations lead to a finite time breakdown of that flow whilst smaller disturbances decay to zero. If the stagnation point flow at infinity is sufficiently strong, the steady basic states become linearly unstable. In fact there is then a continuous spectrum of unstable eigenvalues of the stability equations but, if the initial value problem is considered, it is found that, at large values of time, the continuous spectrum leads to a velocity field growing exponentially in time with an amplitude decaying algebraically in time.

Hall, Philip↗

Details of Exact Low Prandtl Number Boundary-Layer Solutions for Forced and For Free Convection

A detailed report is given of exact (numerical) solutions of the laminar-boundary-layer equations for the Prandtl number range appropriate to liquid metals (0.003 to 0.03). Consideration is given to the following situations: (1) forced convection over a flat plate for the conditions of uniform wall temperature and uniform wall heat flux, and (2) free convection over an isothermal vertical plate. Tabulations of the new solutions are given in detail. Results are presented for the heat-transfer and shear-stress characteristics; temperature and velocity distributions are also shown. The heat-transfer results are correlated in terms of dimensionless parameters that vary only slightly over the entire liquid-metal range. Previous analytical and experimental work on low Prandtl number boundary layers is surveyed and compared with the new exact solutions.

Sparrow, E. M.↗

A two-step iterative method for evolving nonlinear acoustic systems to a steady-state

A new approach for evolving two-dimensional nonlinear acoustic systems with flow to a steady state is presented. The approach is a two-step iterative method which is tested on a benchmark acoustic problem for which an exact analytical solution is available. Results are also calculated for a nonlinear acoustic problem for which an exact analytical solution is not known. Results indicate that the two-step method represents a powerful, efficient, and stable method for evolving two-dimensional acoustic systems to a steady state, and that the method is applicable to any number of spatial dimensions and to other hyperbolic systems. It is noted that for the benchmark problem only a single iteration on the method is required when the transient and steady-state field are of the same order of magnitude; however, four iterations are required when the steady-state field is several orders of magnitude smaller than the transient field. This method requires six iterations before achieving a steady state for the nonlinear test problem.

Watson, Willie R.↗