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 73 records · Page 4

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion

The classical Landau/Levich models of liquid propellant combustion, which serve as seminal examples of hydrodynamic instability in reactive systems, have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and/or temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity, surface tension and viscosity on the hydrodynamic instability of the propagating liquid/gas interface. In particular, a composite asymptotic expression, spanning three distinguished wavenumber regimes, is derived for both cellular and pulsating hydrodynamic neutral stability boundaries A(sub p)(k), where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. For the case of cellular (Landau) instability, the results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limiting case of weak gravity, it is shown that cellular hydrodynamic instability in this context is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(l) wavenumber disturbances. It is also demonstrated that, in the large wavenumber regime, surface tension and both liquid and gas viscosity all produce comparable stabilizing effects in the large-wavenumber regime, thereby providing significant modifications to previous analyses of Landau instability in which one or more of these effects were neglected. In contrast, the pulsating hydrodynamic stability boundary is found to be insensitive to gravitational and surface-tension effects, but is more sensitive to the effects of liquid viscosity, which is a significant stabilizing effect for O(l) and higher wavenumbers. Liquid-propellant combustion is predicted to be stable (i.e., steady and planar) only for a range of negative pressure sensitivities that lie between the two types of hydrodynamic stability boundaries.

Margolis, Stephen B.↗

Development of a near-wall Reynolds-stress closure based on the SSG model for the pressure strain

In this research, a near-wall second-order closure based on the Speziable et al.(1991) or SSG model for the pressure-strain term is proposed. Unlike the LRR model, the SSG model is quasi-nonlinear and yields better results when applied to calculate rotating homogeneous turbulent flows. An asymptotic analysis near the wall is applied to both the exact and modeled, equations so that appropriate near-wall corrections to the SSG model and the modeled dissipation-rate equation can be derived to satisfy the physical wall boundary conditions as well as the asymptotic near-wall behavior of the exact equations. Two additional model constants are introduced and they are determined by calibrating against one set of near-wall channel flow data. Once determined, their values are found to remain constant irrespective of the type of flow examined. The resultant model is used to calculate simple turbulent flows, near separating turbulent flows, complex turbulent flows and compressible turbulent flows with a freestream Mach number as high as 10. In all the flow cases investigated, the calculated results are in good agreement with data. This new near-wall model is less ad hoc, physically and mathematically more sound and eliminates the empiricism introduced by Zhang. Therefore, it is quite general, as demonstrated by the good agreement achieved with measurements covering a wide range of Reynolds numbers and Mach numbers.

So, R. M. C.↗

Effects of Nose Radius and Aerodynamic Loading on Leading Edge Receptivity

An analysis is presented of the effects of airfoil thickness and mean aerodynamic loading on boundary-layer receptivity in the leading-edge region. The case of acoustic free-stream disturbances, incident on a thin cambered airfoil with a parabolic leading edge in a low Mach number flow, is considered. An asymptotic analysis based on large Reynolds number is developed, supplemented by numerical results. The airfoil thickness distribution enters the theory through a Strouhal number based on the nose radius of the airfoil, S = (omega)tau(sub n)/U, where omega is the frequency of the acoustic wave and U is the mean flow speed. The influence of mean aerodynamic loading enters through an effective angle-of-attack parameter ti, related to flow around the leading edge from the lower surface to the upper. The variation of the receptivity level is analyzed as a function of S, mu, and characteristics of the free-stream acoustic wave. For an unloaded leading edge, a finite nose radius dramatically reduces the receptivity level compared to that for a flat plate, the amplitude of the instability waves in the boundary layer being decreased by an order of magnitude when S = 0.3. Modest levels of aerodynamic loading are found to further decrease the receptivity level for the upper surface of the airfoil, while an increase in receptivity level occurs for the lower surface. For larger angles of attack close to the critical angle for boundary layer separation, a local rise in the receptivity level occurs for the upper surface, while for the lower surface the receptivity decreases. The effects of aerodynamic loading are more pronounced at larger values of S. Oblique acoustic waves produce much higher receptivity levels than acoustic waves propagating downstream parallel to the airfoil chord.

Hammerton, P. W.↗

Analysis of a turning point problem in flight trajectory optimization

The optimal control policy for the aeroglide portion of the minimum fuel, orbital plane change problem for maneuvering entry vehicles is reduced to the solution of a turning point problem for the bank angle control. For this problem a turning point occurs at the minimum altitude of the flight, when the flight path angle equals zero. The turning point separates the bank angle control into two outer solutions that are valid away from the turning point. In a neighborhood of the turning point, where the bank angle changes rapidly, an inner solution is developed and matched with the two outer solutions. An asymptotic analysis of the turning point problem is given, and an analytic example is provided to illustrate the construction of the bank angle control.

Gracey, C.↗

Boundary Conditions for Unsteady Compressible Flows

This paper explores solutions to the spherically symmetric Euler equations. Motivated by the work of Hagstrom and Hariharan and Geer and Pope, we modeled the effect of a pulsating sphere in a compressible medium. The literature available on this suggests that an accurate numerical solution requires artificial boundary conditions which simulate the propagation of nonlinear waves in open domains. Until recently, the boundary conditions available were in general linear and based on nonreflection. Exceptions to this are the nonlinear nonreflective conditions of Thompson, and the nonlinear reflective conditions of Hagstrom and Hariharan. The former are based on the rate of change of the incoming characteristics; the latter rely on asymptotic analysis and the method of characteristics and account for the coupling of incoming and outgoing characteristics. Furthermore, Hagstrom and Hariharan have shown that, in a test situation in which the flow would reach a steady state over a long time, Thompson's method could lead to an incorrect steady state. The current study considers periodic flows and includes all possible types and techniques of boundary conditions. The technique recommended by Hagstrom and Hariharan proved superior to all others considered and matched the results of asymptotic methods that are valid for low subsonic Mach numbers.

Hariharan, S. I.↗

The parallel diffusion of cosmic rays in a random magnetic field.

Within the quasi-linear approximation, the existence of the parallel diffusion coefficient for cosmic rays in a random magnetic field (homogeneous, isotropic), despite the slow decay of the interaction between particles and random field, is demonstrated. As an example, the results of a numerical calculation of the parallel diffusion coefficient for a Gaussian random-field correlation function are presented. The numerical results are corroborated by asymptotic analysis and are compared to those of other theories.

Klimas, A.↗

Laminar flow at a three-dimensional stagnation point with large rates of injection

Exact calculations of the titled flow are presented and compared to the predictions of an asymptotic analysis for large rates of injection. The inner layer of the boundary layer is found to involve outflow in both orthogonal directions whether the external flow along the y axis is inward or outward. As a result, the flow at a nearly two-dimensional stagnation point involves drastic changes as a weak outflow changes to a weak inflow. It is also found that the velocity profiles in the two directions in the inner layer are quite different.

Libby, P. A.↗

Observational dispersion relations

Pulse propagation in plasmas has been widely studied in relation to instabilities and their classification as convective or absolute. Such analyses emphasize the importance of considering the summation of a spectrum of modes, rather than simply making deductions from Brillouin dispersion relations plotted either for real frequency and complex wavenumber, or for complex frequency and real wavenumber. In particular, erroneous conclusions may be drawn about the magnitude and direction of propagation of a transient pulse. In this paper, the point of view is taken that an observer remote from the source would measure an 'observational dispersion relation'. The nature of this dispersion relation is examined for absorbing or convectively unstable plasmas by asymptotic analysis, and for different source waveforms. Illustrative results are presented for propagation in cold plasmas and for beam-plasma interaction.

Crawford, F. W.↗

Stationarity of magnetohydrodynamic fluctuations in the solar wind

Solar wind research and studies of charged particle propagation often assume that the interplanetary magnetic field represents a stationary random process. The extent to which ensemble averages of the solar wind magnetic fields follow the asymptotic behavior predicted by the ergodic theorem was investigated. Several time periods, including a span of nearly two years, are analyzed. Data intervals which span many solar rotations satisfy the conditions of weak stationarity if the effects of solar rotation are included in the asymptotic analysis. Shorter intervals which include a small integral number of interplanetary sectors also satisfy weak stationarity. The results are illustrated using magnetometer data from the ISEE-3, Voyager and IMP spacecraft.

Matthaeus, W. H.↗

Flow in a differentially rotated cylindrical drop at low Reynolds number

A liquid drop held captive between parallel disks that are differentially rotated is a model for the swirling flows induced by crystal rotation in the floating-zone process for growing semiconductor materials. An asymptotic analysis for a cylindrical drop is presented that elucidates the structure of the axisymmetric cellular motions caused by disk rotation at low Reynolds number. Variations of meniscus shape induced by these flows are described in the limit of small capillary number. Most cellular flow fields break the bifurcation point that corresponds to the Plateau-Rayleigh limit for the length of a static drop into two disjoint shape families and lower the maximum stable drop length. This effect is studied by a singular bifurcation analysis.

Harriott, G. M.↗

Numerical solution of the Navier-Stokes equations for compressible turbulent two/three dimensional flows in the terminal shock region of an inlet/diffuser

The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique have been used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations have been carried out: (1) 2-D, steady, subsonic; (2) 2-D, steady, transonic with normal shock, (3) 2-D, steady, supersonic with terminal shock, (4) 2-D, transient process of shock development and (5) 3-D, steady, transonic with normal shock. The numerical results obtained for the 2-D and 3-D transonic shocked flows have been compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data. The predicted transient responses of the flow to externally applied impulsive disturbance are consistent with results obtained via asymptotic analysis.

Liu, N.-S.↗

Flow in a differentially rotated cylindrical drop at moderate Reynolds number

Galerkin finite-element approximations are combined with computer-implemented perturbation methods for tracking families of solutions to calculate the steady axisymmetric flows in a differentially rotated cylindrical drop as a function of Reynolds number Re, drop aspect ratio and the rotation ratio between the two end disks. The flows for Reynolds numbers below 100 are primarily viscous and reasonably described by an asymptotic analysis. When the disks are exactly counter-rotated, multiple steady flows are calculated that bifurcate to higher values of Re from the expected solution with two identical secondary cells stacked symmetrically about the axial midplane. The new flows have two cells of different size and are stable beyond the critical value Re sub c. The slope of the locus of Re sub c for drops with aspect ratio up to 3 disagrees with the result for two disks of infinite radius computed assuming the similarity form of the velocity field. Changing the rotation ratio for exact counter-rotation ruptures the junction of the multiple flow fields into two separated flow families.

Harriott, G. M.↗

Numerical studies of interacting vortices

To get a basic understanding of the physics of flowfields modeled by vortex filaments with finite vortical cores, systematic numerical studies of the interactions of two dimensional vortices and pairs of coaxial axisymmetric circular vortex rings were made. Finite difference solutions of the unsteady incompressible Navier-Stokes equations were carried out using vorticity and stream function as primary variables. Special emphasis was placed on the formulation of appropriate boundary conditions necessary for the calculations in a finite computational domain. Numerical results illustrate the interaction of vortex filaments, demonstrate when and how they merge with each other, and establish the region of validity for an asymptotic analysis.

Liu, G. C.↗

A panel method study of vortex sheets with special emphasis on sheets of axisymmetric geometry

The self induced evolution of a vortex sheet was simulated by modeling the sheet using an integration of discrete elements of vorticity. Replacing small sections of a vortex sheet by flat panels of constant vorticity is found to reproduce more accurately the initial conditions for the Lagrangian simulation technique than replacement by point vortices. The flat panel method for the vortex sheet was then extended to model axisymmetric vortex sheets. The local and far field velocities induced by the axisymmetric panels were obtained using matched asymptotic analysis, and some of the uncertainties involved in other models of the axisymmetric vortex sheet have been eliminated. One important result of this analysis is the determination of the proper choice of core size for a circular vortex filament which may replace a section of an axisymmetric vortex sheet. Roll-up of both two dimensional and axisymmetric vortex sheets was computed using the panel methods developed in the report.

Sugioka, I.↗

Heat and mass transfer in combustion - Fundamental concepts and analytical techniques

Fundamental combustion phenomena and the associated flame structures in laminar gaseous flows are discussed on physical bases within the framework of the three nondimensional parameters of interest to heat and mass transfer in chemically-reacting flows, namely the Damkoehler number, the Lewis number, and the Arrhenius number which is the ratio of the reaction activation energy to the characteristic thermal energy. The model problems selected for illustration are droplet combustion, boundary layer combustion, and the propagation, flammability, and stability of premixed flames. Fundamental concepts discussed include the flame structures for large activation energy reactions, S-curve interpretation of the ignition and extinctin states, reaction-induced local-similarity and non-similarity in boundary layer flows, the origin and removal of the cold boundary difficulty in modeling flame propagation, and effects of flame stretch and preferential diffusion on flame extinction and stability. Analytical techniques introduced include the Shvab-Zeldovich formulation, the local Shvab-Zeldovich formulation, flame-sheet approximation and the associated jump formulation, and large activation energy matched asymptotic analysis. Potentially promising research areas are suggested.

Law, C. K.↗

Damage growth in composite laminates with interleaves

The influence of placing interleaves between fiber reinforced plies in multilayered composite laminates is investigated. The geometry of the composite is idealized as two dimensional, isotropic, linearly elastic media made of a damaged layer bonded between two half planes and separated by thin interleaves of low extensional and shear moduli. The damage in the layer is taken in the form of a symmetric crack perpendicular to the interface and may extend up to the interface. The case of an H-shaped crack in the form of a broken layer with delamination along the interface is also analyzed. The interleaves are modeled as distributed shear and tension springs. Fourier integral transform techniques are used to develop solutions in terms of singular integral equations. An asymptotic analysis of the integral equations based on Muskhelishvili's techniques reveals logarithmically singular axial stresses in the half plane at the crack tips for the broken layer. For the H shaped crack, similar singularities are found to exist in the axial stresses at the interface crack tips in the layer and the half plane. The solution of the equations is found numerically for the stresses and displacements by using the Hadamard's concept of direct differentiation of Cauchy integrals as well as Gaussian integration techniques.

Goree, James G.↗

Numerical studies of interacting vortices

To get a basic understanding of the physics of flowfields modeled by vortex filaments with finite vortical cores, systematic numerical studies of the interactions of two vortices and pairs of coaxial axisymmetric circular vortex rings were made. Finite difference solutions of the unsteady incompressible Navier-Stokes equations were carried out using vorticity and stream function as primary variables. Special emphasis was placed on the formulation of appropriate boundary conditions necessary for the calculations in a finite computational domain. Numerical results illustrate the interaction of vortex filaments, demonstrate when and how they merge with each other, and establish the region of validity for an asymptotic analysis.

Liu, G. C.↗