Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ASYMPTOTIC FUNCTION”

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 109 records · Page 6

Navier-Stokes calculations for the vortex of a rotor in hover

An efficient finite-difference scheme for the solution of the incompressible Navier-Stokes equation is used to study the vortex wake of a rotor in hover. The solution Procedure uses a vorticity-stream function formulation and incorporates an asymptotic far-field boundary condition enabling the size of the computational domain to be reduced in comparison to other methods. The results from the present method are compared with experimental data obtained by smoke flow visualization and hot-wire measurements for several rotor blade configurations.

Liu, C. H.↗

A wave envelope finite element scheme for acoustical radiation

The aeroacoustic problem associated with the radiation of fan noise from the inlet of a turbofan aircraft engine, the dimensions of which are generally many times larger than the acoustical wavelengths of the major energy-carrying frequencies, is considered. In the present approach, a conventional finite element solution in the inner region is compatibly matched to a 'wave envelope' finite element solution in a large but finite outer region. The inclusion of a wavelike variation with the correct asymptotic decay in the shape functions for the outer region preserves the correct behavior of the solution at large distances. The method is initially presented for a simple one-dimensional model based on the solution of Webster's horn equation. Results are presented for the specific case of a uniform cylindrical section joined to a conical expansion, and also for a simple axisymmetric test case of the calculation of acoustical pressure generated by a vibrating circular piston located at the center of an infinite rigid baffle.

Astley, R. J.↗

Stability of an oscillated fluid with a uniform density gradient

Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Jacqmin, David↗

Analysis of the two-point velocity correlations in turbulent boundary layer flows

The general objective of the present work is to explore the use of Rapid Distortion Theory (RDT) in analysis of the two-point statistics of the log-layer. RDT is applicable only to unsteady flows where the non-linear turbulence-turbulence interaction can be neglected in comparison to linear turbulence-mean interactions. Here we propose to use RDT to examine the structure of the large energy-containing scales and their interaction with the mean flow in the log-region. The contents of the work are twofold: First, two-point analysis methods will be used to derive the law-of-the-wall for the special case of zero mean pressure gradient. The basic assumptions needed are one-dimensionality in the mean flow and homogeneity of the fluctuations. It will be shown that a formal solution of the two-point correlation equation can be obtained as a power series in the von Karman constant, known to be on the order of 0.4. In the second part, a detailed analysis of the two-point correlation function in the log-layer will be given. The fundamental set of equations and a functional relation for the two-point correlation function will be derived. An asymptotic expansion procedure will be used in the log-layer to match Kolmogorov's universal range and the one-point correlations to the inviscid outer region valid for large correlation distances.

Oberlack, M.↗

Global-Local Finite Element Analysis for Thermo-Mechanical Stresses in Bonded Joints

An analysis of adhesively bonded joints using conventional finite elements does not capture the singular behavior of the stress field in regions where two or three dissimilar materials form a junction with or without free edges. However, these regions are characteristic of the bonded joints and are prone to failure initiation. This study presents a method to capture the singular stress field arising from the geometric and material discontinuities in bonded composites. It is achieved by coupling the local (conventional) elements with global (special) elements whose interpolation functions are constructed from the asymptotic solution.

Shkarayev, S.↗

CCKT Calculation of e-H Total Cross Sections

We are in the process of carrying out calculations of e-H total cross sections using the 'complex-correlation Kohn-T' (CCKT) method. In a later paper, we described the methodology more completely, but confined calculations to the elastic scattering region, with definitive, precision results for S-wave phase shifts. Here we extend the calculations to the (low) continuum (1 much less than k(exp 2) much less than 3) using a Green's function formulation. This avoids having to solve integro-differential equations; rather we evaluate indefinite integrals involving appropriate Green's functions and the (complex) optical potential to find the scattering function u(r). From the asymptotic form of u(r) we extract a T(sub L) which is a complex number. From T(sub L), elastic sigma(sub L)(elastic) = 4pi(2L+1)((absolute value of T(sub L))(exp 2)), and total sigma (sub L)(total) = 4pi/k(2L+1)Im(T(sub L)) cross sections follow.

Bhatia, Aaron K.↗

Human Mars Mission: Launch Window from Earth Orbit

The determination of orbital window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a DELTA V penalty. Usually, because of the DELTA V penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: (1) One impulsive maneuver from a Highly Elliptical Orbit (HEO) (2) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO) (3) One impulsive maneuver from a Low Earth Orbit (LEO) (4) Two impulsive maneuvers from LEO (5) Three impulsive maneuvers from LEO.

Young, Archie↗

Human Exploration Missions Study Launch Window from Earth Orbit

The determination of orbital launch window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a Delta(V) penalty. Usually, because of the Delta(V) penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: (1) One impulsive maneuver from a Low Earth Orbit (LEO), (2) Two impulsive maneuvers from LEO, (3) Three impulsive maneuvers from LEO, (4) One impulsive maneuvers from a Highly Elliptical Orbit (HEO), (5) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO) The formulation of these five different launch window modes provides a rapid means of generating realistic parametric data for space exploration studies. Also the formulation provides vector and geometrical data sufficient for use as a good starting point in detail trajectory analysis based on calculus of variations, steepest descent, or parameter optimization program techniques.

Young, Archie↗

Locating the Discontinuities of a Bounded Function by the Partial Sums of its Fourier Series I: Periodical Case

A key step for some methods dealing with the reconstruction of a function with jump discontinuities is the accurate approximation of the jumps and their locations. Various methods have been suggested in the literature to obtain this valuable information. In the present paper, we develop an algorithm based on identities which determine the jumps of a 2(pi)-periodic bounded not-too-highly oscillating function by the partial sums of its differentiated Fourier series. The algorithm enables one to approximate the locations of discontinuities and the magnitudes of jumps of a bounded function. We study the accuracy of approximation and establish asymptotic expansions for the approximations of a 27(pi)-periodic piecewise smooth function with one discontinuity. By an appropriate linear combination, obtained via derivatives of different order, we significantly improve the accuracy. Next, we use Richardson's extrapolation method to enhance the accuracy even more. For a function with multiple discontinuities we establish simple formulae which "eliminate" all discontinuities of the function but one. Then we treat the function as if it had one singularity following the method described above.

Kvernadze, George↗

The effect of sensor spacing on wind measurements at the Shuttle Landing Facility

This document presents results of a field study of the effect of sensor spacing on the validity of wind measurements at the Space Shuttle landing Facility (SLF). Standard measurements are made at one second intervals from 30 foot (9.1m) towers located 500 feet (152m) from the SLF centerline. The centerline winds are not exactly the same as those measured by the towers. This study quantifies the differences as a function of statistics of the observed winds and distance between the measurements and points of interest. The field program used logarithmically spaced portable wind towers to measure wind speed and direction over a range of conditions. Correlations, spectra, moments, and structure functions were computed. A universal normalization for structure functions was devised. The normalized structure functions increase as the 2/3 power of separation distance until an asymptotic value is approached. This occurs at spacings of several hundred feet (about 100m). At larger spacings, the structure functions are bounded by the asymptote. This enables quantitative estimates of the expected differences between the winds at the measurement point and the points of interest to be made from the measured wind statistics. A procedure is provided for making these estimates.

Merceret, Francis J.↗

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quasar evolution and the growth of black holes

A 'minimalist' model of AGN evolution is analyzed that links the measured luminosity function to an elementary description of black hole accretion. The observed luminosity function of bright AGN is extrapolated and simple prescriptions for the growth and luminosity of black holes are introduced to infer quasar birth rates, mean fueling rates, and relict black hole distribution functions. It is deduced that the mean accretion rate scales as (M exp -1./5)(t exp -6.7) and that, for the most conservative model used, the number of relict black holes per decade declines only as M exp -0.4 for black hole masses between 3 x 10 exp 7 and 3 x 10 exp 9 solar masses. If all sufficiently massive galaxies pass through a quasar phase with asymptotic black hole mass a monotonic function of the galaxy mass, then it is possible to compare the space density of galaxies with estimated central masses to that of distant quasars.

Small, Todd A.↗

Human Mars Mission: Launch Window from Earth Orbit

The determination of orbital window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to the earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a delta V penalty. Usually, because of the delta V penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: 1) One impulsive maneuver from a Highly Elliptical Orbit (HEO); 2) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO); 3) One impulsive maneuver from a Low Earth Orbit (LEO); 4) Two impulsive maneuvers form LEO; and 5) Three impulsive maneuvers form LEO. The formulation of these five different launch window modes provides a rapid means of generating realistic parametric data for space exploration studies. Also the formulation provides vector and geometrical data sufficient for use as a good starting point in detail trajectory analysis based on calculus of variations, steepest descent, or parameter optimization program techniques.

Young, Archie↗

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pulsed electromagnetic gas acceleration

Terminal voltage measurements with long cathodes in a high power, quasi-steady MPD discharge show that the critical current for the onset of voltage fluctuations, which was previously shown to be a function of cathode area, approaches an asymptote for cathodes of very large surface area. Floating potential measurements and photographs of the discharge luminosity indicate that the fluctuations are confined to the vicinity of the cathode and hence reflect a cathode emission process rather than a fundamental limit on MPD performance. Photoelectric measurements of particular argon neutral and ion transitions show that the higher electronic states are populated more heavily than would be calculated on the basis of Saha-Boltzmann equilibrium at the local electron temperature and number density. Preliminary optical depth measurements show that for a current of 4 kA and an argon mass flow of 12 g/sec, a population inversion exists between the upper and lower states of the 4880 A argon ion transition.

Jahn, R. G.↗