Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “asymptotic analysis”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Pseudospectral collocation methods for fourth order differential equations

Collocation schemes are presented for solving linear fourth order differential equations in one and two dimensions. The variational formulation of the model fourth order problem is discretized by approximating the integrals by a Gaussian quadrature rule generalized to include the values of the derivative of the integrand at the boundary points. Collocation schemes are derived which are equivalent to this discrete variational problem. An efficient preconditioner based on a low-order finite difference approximation to the same differential operator is presented. The corresponding multidomain problem is also considered and interface conditions are derived. Pseudospectral approximations which are C1 continuous at the interfaces are used in each subdomain to approximate the solution. The approximations are also shown to be C3 continuous at the interfaces asymptotically. A complete analysis of the collocation scheme for the multidomain problem is provided. The extension of the method to the biharmonic equation in two dimensions is discussed and results are presented for a problem defined in a nonrectangular domain.

Malek, Alaeddin↗
DOE OSTI2020

Steady detonation propagation in thin channels with strong confinement

We examine asymptotically the dynamics of two-dimensional, steady detonation wave propagation and failure for a strongly confined high explosive (HE), in which the width of the explosive is small relative to the reaction zone length. An energy balance equation is derived, which shows how the longitudinal acceleration of subsonic flow behind the detonation shock is influenced both by chemical reaction and by the effects of HE boundary streamline deflection, specifically via the induced rate of change of mass flux through the detonation wave. The latter serves to either counteract or reinforce the acceleration of longitudinal flow, depending on the sign of the gradient of the boundary streamline deflection at the detonation shock. The analysis is valid for general equations of state and chemical reaction rates in the HE. The asymptotically derived form of the energy equation represents an eigenvalue problem for the determination of the steady detonation propagation speed, solved via a shooting method. We explore specific results for ideal and stiffened equations of state, along with a pressure-dependent reaction rate for which changes in the pressure exponent and reaction order are also studied. We consider the influences of both straight and curved HE boundary streamline shapes. The asymptotic analysis reveals significant physical insights into how detonation propagation and failure are affected by strong confinement.

42 ENGINEERING↗
NASA NTRS1998

Asymptotic Spreading Rate of Initially Compressible Jets: Experiment and Analysis

Experimental results for the spreading and centerline velocity decay rates for round, compressible jets, from a convergent and a convergent-divergent nozzle, are presented. The spreading rate is determined from the variation of streamwise mass flux obtained from Pitot probe surveys. Results for the far asymptotic region show that both spreading and centerline velocity decay rates, when nondimensionalized by parameters at the nozzle exit, decrease with increasing "jet Mach number" M(sub J). Dimensional analysis with the assumption of momentum conservation, together with compressible flow calculations for the conditions at the nozzle exit, predict this Mach number dependence well. The analysis also demonstrates that an increase in the "potential core length" of the jet occurring with increasing M(sub J), a commonly observed trend, is largely accounted for simply by the variations in the density and static pressure at the nozzle exit. The effect of decreasing mixing efficiency with increasing compressibility is inferred to contribute only partially to the latter trend.

Zaman, K. B. M. Q.↗
NASA NTRS1979

A numerical comparison with an exact solution for the transient response of a cylinder immersed in a fluid

The transient response of an elastic cylindrical shell immersed in an acoustic media that is engulfed by a plane wave is determined numerically. The method applies to the USA-STAGS code which utilizes the finite element method for the structural analysis and the doubly asymptotic approximation for the fluid-structure interaction. The calculations are compared to an exact analysis for two separate loading cases: a plane step wave and an exponentially decaying plane wave.

Giltrud, M. E.↗
NASA NTRS1991

Asymptotic structure of hydromagnetically driven relativistic winds

A fully relativistic analysis has been performed of the asymptotic structure of stationary axisymmetric hydromagnetic winds. If a flow fills the region containing the rotation axis, then the flux surfaces in the flow must collimate to a set of current-carrying cylindrical surface extending to infinite transverse radius, collimate to a set of cylindrical surfaces extending to a finite radius, or collimate to a current-free paraboloidal field configuration which fills up the entire space. If an asymptotically cylindrical flow carries a finite current at radii well beyond the light cylinder, then the Lorentz factor of the terminal flow speed on a given flux surface is proportional to the total current enclosed within this flux surface. If a flow is of type II paraboloidal, then its asymptotic energy flux is carried entirely by the gas motion rather than the electromagnetic fields.

Chiueh, Tzihong↗
NASA NTRS1992

Nature of buoyancy-driven flows in a reduced-gravity environment

The role of buoyancy-driven convection in reduced-gravity environments has been emphasized. It is shown that for some materials-processing experiments the values of the fundamental dimensionless parameters such as Gr or Ra are shifted from the very large (ground-based values) to moderately large (space-based values). As a consequence, in cases where approximate analytical solutions are desired, the accuracy of the usual asymptotic (boundary-layer) analysis in which infinitely large parametric ranges are assumed is reduced. Approximate analytical techniques for refinement of the asymptotic solutions are identified that will extend the accuracy to the moderately large parametric ranges associated with space-based processing.

Kassemi, Siavash A.↗
NASA NTRS1983

Helicopter rotor loads using discretized matched asymptotic expansions

The numerical practicality of a matched asymptotic expansion approach for the computation of unsteady three dimensional airloads on a helicopter rotor was improved. This effort utilizes a discretized repesentation of the doublet strength distribution and helical streamlines to decrease the computational requirements of the original analysis. The continuous variation of the doublet strength was approximated by piecewise constant or piecewise quadratic distributions, and the helical trajectory of a fluid particle was approximated by connected straight line segments. As a direct result of these simplified representations the computational time required for the execution of a typical flight condition was reduced by an order of magnitude with respect to the requirements of the original analysis. Airloads which were computed using the discretized method for a two bladed model rotor and a full scale four bladed rotor are in close agreement with measured results and airloads from the original asymptotic analysis. For conditions characterized by significant rotor/wake interaction the piecewise constant representation requires a reduced azimuth spacing to maintain acceptable accuracy.

Pierce, G. A.↗
NASA NTRS1998

Asymptotic Spreading Rate of Initially Compressible Jets-Experiment and Analysis

Experimental results for the spreading and centerline velocity decay rates for round, compressible jets, from a convergent and a convergent-divergent nozzle, are presented. The spreading rate is determined from the variation of streamwise mass flux obtained from Pitot probe surveys. Results for the far asymptotic region show that both spreading and centerline velocity decay rates, when nondimensionalized by parameters at the nozzle exit, decrease with increasing "jet Mach number" M(sub j). Dimensional analysis with the assumption of momentum conservation, together with compressible flow calculations for the conditions at the nozzle exit, predict this Mach number dependence well. The analysis also demonstrates that an increase in the "potential core length" of the jet occurring with increasing M(sub j), a commonly observed trend, is largely accounted for simply by the variations in the density and static pressure at the nozzle exit. The effect of decreasing mixing efficiency with increasing compressibility is inferred to contribute only partially to the latter trend.

Zaman, K. B. M. Q.↗
NASA NTRS1983

The evolution of Tollmien-Schlichting waves near a leading edge

The method of matched asymptotic expansions is used to study the generation of Tollmien-Schlichting waves by free-stream disturbances incident on a flat-plate boundary layer. Near the leading edge, the motion is governed by the unsteady boundary-layer equation, while farther downstream it is governed (to lowest order) by the Orr-Sommerfeld equation with slowly varying coefficients. It is shown that there is an overlap domain where the Tollmien-Schlichting wave solutions to the Orr-Sommerfeld equation and appropriate asymptotic solutions of the unsteady boundary-layer equation match, in the matched-asymptotic-expansion sense. The analysis explains how long-wavelength free-stream disturbances can generate Tollmien-Schlichting waves of much shorter wavelength. It also leads to a set of scaling laws for the asymptotic structure of the unsteady boundary layer.

Goldstein, M. E.↗
DOE OSTI2021

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗
NASA NTRS1985

An asymptotic theory for the interference of large aspect ratio swept wings and multiple propeller slipstreams

This paper presents an asymptotic method for the analysis of the interference of multiple tractor propeller slipstream with large aspect ratio swept wings. It is assumed that the height of the slipstream is of the order of the wing chord and its spanwise extent is of the order of the wing span. Three different flow regions are identified by employing different stretching transformations. Asymptotic expansions are made in each of the three regions, using the chord to span ratio as the small expansion parameter. The details of the nonuniform flow in the slipstream enter into the wing sectional analysis. In the outer limit, the wing shrinks to a swept lifting line, and the slipstream reduces to a thin sheet of jet carrying the momentum gain from the propeller. The curvature of this jet sheet results in a pressure difference which is represented by a vortex sheet. The governing equations are solved by discretization. Several examples are considered for which experimental data are available. Comparison of the present results with the experimental data as well as other numerical solutions showed generally good agreement.

Prabhu, R. K.↗
NASA NTRS1994

A uniform geometrical optics and an extended uniform geometrical theory of diffraction for evaluating high frequency EM fields near smooth caustics and composite shadow boundaries

A uniform geometrical optics (UGO) and an extended uniform geometrical theory of diffraction (EUTD) are developed for evaluating high frequency electromagnetic (EM) fields within transition regions associated with a two and three dimensional smooth caustic of reflected rays and a composite shadow boundary formed by the caustic termination or the confluence of the caustic with the reflection shadow boundary (RSB). The UGO is a uniform version of the classic geometrical optics (GO). It retains the simple ray optical expressions of classic GO and employs a new set of uniform reflection coefficients. The UGO also includes a uniform version of the complex GO ray field that exists on the dark side of the smooth caustic. The EUTD is an extension of the classic uniform geometrical theory of diffraction (UTD) and accounts for the non-ray optical behavior of the UGO reflected field near caustics by using a two-variable transition function in the expressions for the edge diffraction coefficients. It also uniformly recovers the classic UTD behavior of the edge diffracted field outside the composite shadow boundary transition region. The approach employed for constructing the UGO/EUTD solution is based on a spatial domain physical optics (PO) radiation integral representation for the fields which is then reduced using uniform asymptotic procedures. The UGO/EUTD analysis is also employed to investigate the far-zone RCS problem of plane wave scattering from two and three dimensional polynomial defined surfaces, and uniform reflection, zero-curvature, and edge diffraction coefficients are derived. Numerical results for the scattering and diffraction from cubic and fourth order polynomial strips are also shown and the UGO/EUTD solution is validated by comparison to an independent moment method (MM) solution. The UGO/EUTD solution is also compared with the classic GO/UTD solution. The failure of the classic techniques near caustics and composite shadow boundaries is clearly demonstrated and it is shown that the UGO/EUTD results remain valid and uniformly reduce to the classic results away from the transition regions. Mathematical details on the asymptotic properties and efficient numerical evaluation of the canonical functions involved in the UGO/EUTD expressions are also provided.

Constantinides, E. D.↗
NASA NTRS1990

The inviscid stability of supersonic flow past heated or cooled axisymmetric bodies

The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions, are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.

Shaw, Stephen J.↗
NASA NTRS1992

The inviscid stability of supersonic flow past heated or cooled axisymmetric bodies

The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.

Shaw, Stephen J.↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1974

Aerodynamic influence coefficient method using singularity splines

A numerical lifting surface formulation, including computed results for planar wing cases is presented. This formulation, referred to as the vortex spline scheme, combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of loading function methods. The formulation employes a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfied all the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise. Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral. The current formulation uses the elementary horseshoe vortex as the basic singularity and is therefore restricted to linearized potential flow. As part of the study, a non planar development was considered, but the numerical evaluation of the lifting surface concept was restricted to planar configurations. Also, a second order sideslip analysis based on an asymptotic expansion was investigated using the singularity spline formulation.

Mercer, J. E.↗
NASA NTRS1984

Progress in wind tunnel wall interference assessment/correction procedures at the NAE

Wall corrections as a function of wall porosity in the transonic wall interference problem was assessed. Effective porosities primarily for the two dimensional case were established as follows: (1) comparison of experimental data for two geometrically similar models of different chord/height ratio, an overall value of wall porosity could be deduced; (2) theoretical development which allows for unequal porosity for the floor and ceiling and wall boundary pressure measurements, porosities for floor and ceiling could be deduced; (3) a scheme was developed which allowed unequal porosity of floor and ceiling and streamwise varying porosity. The boundary layer development along the perforated floor and ceiling under the influence of the model pressure field, variations in boundary layer thickness underlining the difficulties in deducing meaningful values of wall porosity were determined. Wall boundary pressure measurement, in combination with singularity modelling of the airfoil, was sufficient to yield required information on the wall interference flow without having to establish some value for wall porosity. The singularity modelling of the airfoil initially covered only lift and volume but was extended to include drag and pitching moment, and second order volume term. It is shown by asymptotic transonic small disturbance analysis, that the derived corrections to angle of attack and free stream Mach number are correct to the first order.

Ohman, L. H.↗
NASA NTRS1987

Thermochemical erosion

A quasi-steady multidimensional Stefan problem is treated by means of similarity, asymptotic methods, and numerical analysis. The continuous wearing away of a solid 'erodee' is effected by gradual consumption of a solid 'eroder'. In the process under study, the erodee melts and the (liquid) melt diffuses to the eroder, where an exothermic, indefinitely rapid, diffusion-controlled surface reaction occurs. The products of reaction convey some of the chemically-generated heat to the surface of the erodee, so that further melting may occur, and the process may proceed. The analysis seeks to estimate the quasi-steady rates of consumption of the erodee and eroder, for the (convenient) self-similar case of confocal parabolas defining the boundaries of the melt layer. The analysis is simplified by observing that the flow is of lubrication-theory type: creeping, primarily unidirectional motion owing to the thinness of the melt layer between the erodee and eroder, with diffusion of heat and mass predominantly perpendicular to the principal motion. The analysis is completed by consideration of the force balance on the eroder.

Butler, Gerald↗