Search NASA⌕ Search

SEARCH · Search NASA

Results for “THEORIES - FLOW, LINEAR”

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 271 records · Page 15

Physics of pitch angle scattering and velocity diffusion. I - Theory

A general theory for the pitch angle scattering and velocity diffusion of particles in the field of a spectrum of waves in a magnetized plasma is presented. The test particle theory is used to analyze the particle motion. The form of diffusion surfaces is examined, and analytical expressions are given for the resonance width and bounce frequency. The resonance widths are found to vary strongly as a function of harmonic number. The resulting diffusion can be quite asymmetric with respect to pitch angle of 90 deg. The conditions for the onset of pitch angle scattering and energy diffusion are explained in detail. Some of the known shortcomings of the standard quasi-linear theory are also addressed, and ways to overcome them are shown. In particular, the often stated quasi-linear gap at 90 deg is found to exist only under very special cases. For instance, oblique wave propagation can easily remove the gap. The conditions for the existence of the gap are described in great detail. A new diffusion equation which takes into account the finite resonance widths is also discussed. The differences between this new theory and the standard resonance broadening theory is explained.

Karimabadi, H.↗

Low-freqency instabilities and the resulting velocity distributions of pickup ions at comet Halley

The interaction between the solar wind and newborn cometary ions is studied using a new analytical theory as well as one- and two-dimensional hybrid simulations. Using the observed parameters upstream of the comet Halley, a detailed study of wave excitation and the resulting particle distributions is presentd. Linear theory as well as simulations show that a variety of modes such as the fast magnetosonic mode, high frequency whistlers and obliquely propagating Alfven ion cyclotron waves can be excited. However, parallel propagating waves are found to be dominant in the wave spectrum and to control the scattering of the pickup ions. Several features of the observed distributions of pickup protons are explained. In particular, it is shown that the observed asymmetric pitch distribution for the pickup protons is due to the small saturation amplitude of the waves for the given parameters. Water group associated waves can lead to energy diffusion and further pitch angle scattering of protons. This effect is most likely to be important in the vicinity of the bow shock of comet Halley where the density of water group ions becomes comparable to that of protons. It is shown that the observed increase in the radius of the proton velocity shell just outside the bow shock can be due to water group waves. The nearby isotropic proton pitch angle distribution observed by Neugebauer et al. (1989) just outside the bow shock may, however, be related to the presence of a rotational discontinuity which has been identified in the magnetic field data. Just outside the bow shock, simulations show that parallel propagating water group waves can steepen with attached whistler wave packets. The steepening process at parallel propagation is a transient effect, in an important contrast to the case of steepening at oblique angles. The smaller beam densitites at comet Halley appears to be the main reason not only why waves at comet Halley have smaller amplitudes but also why oblique, steepening magnetosonic waves have not been detected at comet Halley, whereas they have been seen at comet Giacobini-Zinner. Halley, whereas they have been seen at comet Giacobini-Zinner.

Karimabadi, H.↗

Investigation at Mach Numbers of 0.20 to 3.50 of a Blended Diamond Wing and Body Combination of Sonic Design but with Low Wave-Drag Increase with Increasing Mach Number

A diamond wing and body combination was designed to have an area distribution which would result in near optimum zero-lift wave-drag coefficients at a Mach number of 1.00, and decreasing wave-drag coefficient with increasing Mach number up to near sonic leading-edge conditions for the wing. The airfoil section were computed by varying their shape along with the body radii (blending process) to match the selected area distribution and the given plan form. The exposed wing section had an average maximum thickness of about 3 percent of the local chords, and the maximum thickness of the center-line chord was 5.49 percent. The wing had an aspect ratio of 2 and a leading-edge sweep of 45 deg. Test data were obtained throughout the Mach number range from 0.20 to 3.50 at Reynolds numbers based on the mean aerodynamic chord of roughly 6,000,000 to 9,000,000. The zero-lift wave-drag coefficients of the diamond model satisfied the design objectives and were equal to the low values for the Mach number 1.00 equivalent body up to the limit of the transonic tests. From the peak drag coefficient near M = 1.00 there was a gradual decrease in wave-drag coefficient up to M = 1.20. Above sonic leading-edge conditions of the wing there was a rise in the wave-drag coefficient which was attributed in part to the body contouring as well as to the wing geometry. The diamond model had good lift characteristics, in spite of the prediction from low-aspect-ratio theory that the rear half of the diamond wing would carry little lift. The experimental lift-curve slope obtained at supersonic speeds were equal to or greater than the values predicted by linear theory. Similarly the other basic aerodynamic parameters, aerodynamic center position, and maximum lift-drag ratios were satisfactorily predicted at supersonic speeds.

Holdaway, George H.↗

Improving initialization and evolution accuracy of cosmological neutrino simulations

Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations — the impact of massive neutrinos should therefore be included in these simulations in a realistic, computationally tractable, and controlled manner. Here, a recent proposal to reduce the related computational cost employs a symmetric neutrino momentum sampling strategy in the initial conditions. We implement a modified version of this strategy into the Hardware/Hybrid Accelerated Cosmology Code (HACC) and perform convergence tests on its internal parameters. We illustrate that this method can impart Script $\mathcal{O}$(1%) numerical artifacts on the total matter field on small scales, similar to previous findings, and present a method to remove these artifacts using Fourier-space filtering of the neutrino density field. Moreover, we show that the converged neutrino power spectrum does not follow linear theory predictions on relatively large scales at early times at the 15% level, prompting a more careful study of systematics in particle-based neutrino simulations. We also present an improved method for backscaling linear transfer functions for initial conditions in massive neutrino cosmologies that is based on achieving the same relative neutrino growth as computed with Boltzmann solvers. Our self-consistent backscaling method yields sub-percent accuracy in the total matter growth function. Comparisons for the non-linear power spectrum with the Mira-Titan emulator at a neutrino mass of m ν = 0.15 eV are in very good agreement with the expected level of errors in the emulator and in the direct N-body simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear gas oscillations in pipes. I - Theory.

The problem of forced acoustic oscillations in a pipe is studied theoretically. The oscillations are produced by a moving piston in one end of the pipe, while a variety of boundary conditions ranging from a completely closed to a completely open mouth at the other end are considered. The linear theory predicts large amplitudes near resonance and that nonlinear effects become crucially important. By expanding the equations of motion in a series in the Mach number, both the amplitude and waveform of the oscillation are predicted there. In both the open- and closed-end cases the need for shock waves in some range of parameters is found. The amplitude of the oscillation is different for the two cases, however, being proportional to the square root of the piston amplitude in the closed-end case and to the cube root for the open end.

Jimenez, J.↗

Differential geometric methods in system theory.

Discussion of certain problems in system theory which have been or might be solved using some basic concepts from differential geometry. The problems considered involve differential equations, controllability, optimal control, qualitative behavior, stochastic processes, and bilinear systems. The main goal is to extend the essentials of linear theory to some nonlinear classes of problems.

Brockett, R. W.↗

Resonant triad in boundary-layer stability. Part 1: Fully nonlinear interaction

A first principles theory is developed to study the nonlinear spatial evolution of a near-resonance triad of instability waves in boundary layer transition. This triad consists of a plane wave at fundamental frequency and a pair of symmetrical, oblique waves at the subharmonic frequency. A low frequency, high Reynolds number asymptotic scaling leads to a distinct critical layer where nonlinearity first becomes important; the development of the triad's waves is determined by the critical layer's nonlinear, viscous dynamics. The resulting theory is fully nonlinear in that all nonlinearly generated oscillatory and nonoscillatory components are accounted for. The presence of the plane wave initially causes exponential of exponential growth of the oblique waves. However, the plane wave continues to follow the linear theory, even when the oblique waves' amplitude attains the same order of magnitude as that of the plane wave. A fully interactive stage then comes into effect when the oblique waves exceed a certain level compared to that of the plane wave. The oblique waves react back on the fundamental, slowing its growth rate. The oblique waves' saturation results from their self-interaction - a mechanism that does not require the presence of the plane wave. The oblique waves' saturation level is independent of their initial level, but decreases as the obliqueness angle increases.

Mankbadi, Reda R.↗

Diffuse Reflection of Laser Light From Clouds

Laser light reflected from an aqueous suspension of particles or "cloud" with known thickness and particle size distribution defines the "cloud radiative Green's function", G. G is sensitive to cloud thickness, allowing retrieval of that important quantity. We describe a laboratory simulation of G, useful in design of an offbeam Lidar instrument for remote sensing of cloud thickness. Clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The size distribution extends from 0.5 microns to 25 microns, roughly lognormal, similar to real clouds. Density of suspended spheres is adjusted so photon mean-free-path is about 10 cm, 1000 times smaller than in real clouds. The light source is a Nd:YAG laser at 530 nm. Detectors are flux and photon-counting PMTs, with a glass probe for precise positioning. A Labview 5 VI controls position and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider,and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns. Step size is selectable. Far from the beam, the rate of exponential increase in the beam direction scales as expected from diffusion theory, linearly with cloud thickness, and inversely as the square root of the reduced optical thickness, independent of particle size. Nearer the beam the signal increases faster than exponential and depends on particle size. Results verify 3D Monte Carlo simulations that demonstrate detectability of remotely sensed offbeam returns, without filters at night, with narrow bandpass filter in day.

Cahalan, R. F.↗

Anomalous heating of the polar E region by unstable plasma waves. II - Theory

It is found that anomalous electron temperatures in the disturbed high-latitude E region can be quantitatively explained in terms of heating by unstable plasma waves. The electron temperatures at 110 km have been measured to be as high as 1500 K instead of the expected value of about 300 K. It is shown that by using quasi-linear theory there is an ample source of heat in the unstable waves and that the measured electron temperature profiles have a shape very similar to what is expected from plasma wave heating by the modified two-stream instability. It is found that there is even more heating going to the ion gas, but that the resulting effect on the ion temperature may be difficult to measure. The best estimate of the wave heating rates leads to the conclusion that wave heating can be as much as 50% of the Joule heating for dc electric field strengths of the order of 45 mV/m or greater.

St.-Maurice, J. P.↗

A general aerodynamic approach to the problem of decaying or growing vibrations of thin, flexible wings with supersonic leading and trailing edges and no side edges

Indicial aerodynamic influence coefficients were evaluated from potential theory for a thin, flexible wing with supersonic leading and trailing edges only. The analysis is based on the use of small surface areas in which the downwash is assumed uniform. Within this limitation, the results are exact except for the restriction of linearized theory. The areas are not restricted either to square boxes or Mach boxes. A given area may be any rectangle or square which may or may not be cut by the Mach forecone, and any area can be used anywhere in the forecone without loss of accuracy.

Warner, R. W.↗

Numerical simulations of turbulent trapping in the weak beam-plasma instability

Numerical simulations of the weak beam-plasma instability were done in the turbulent regime where small-scale trapping is a dominant feature of the instability, a regime with behavior not predicted by quasi-linear theory. The results of the simulations were compared with those of a specific model of the turbulence, the so-called 'turbulent trapping' model, which gives precise formulas for the particle correlation functions, and predicts a growth rate well enhanced over the quasi-linear value. It was found that the model gives accurate predictions for the correlation functions. On the other hand, while growth rates were enhanced over the quasi-linear values, the enhancements observed are smaller than expected from the quantitative predictions of the model.

Theilhaber, K.↗

The elastic layer with a cylindrical hole subjected to a nonuniform axisymmetric radial displacement.

A problem in the linear theory of elasticity is considered wherein a layer with a circular cylindrical hole is subjected to a nonuniform axisymmetric radial displacement. The solution utilizes Navier's equations of elasticity which are solved by means of extended Hankel transforms. A special case in which the radial displacement is a linear function of the axial coordinate is presented. Numerical results are given in graphical form for the case when hole radius and layer thickness are equal. The inversion integrals were evaluated numerically using Longman's technique for computing infinite integrals of oscillatory functions.

Grissom, D. S.↗

Finite state modeling of aeroelastic systems

A general theory of finite state modeling of aerodynamic loads on thin airfoils and lifting surfaces performing completely arbitrary, small, time-dependent motions in an airstream is developed and presented. The nature of the behavior of the unsteady airloads in the frequency domain is explained, using as raw materials any of the unsteady linearized theories that have been mechanized for simple harmonic oscillations. Each desired aerodynamic transfer function is approximated by means of an appropriate Pade approximant, that is, a rational function of finite degree polynomials in the Laplace transform variable. The modeling technique is applied to several two dimensional and three dimensional airfoils. Circular, elliptic, rectangular and tapered planforms are considered as examples. Identical functions are also obtained for control surfaces for two and three dimensional airfoils.

Vepa, R.↗

Electromagnetic electron beam instabilities - Hot, isotropic beams

This paper considers the linear theory of electromagnetic instabilities driven by an electron beam in a homogeneous, nonrelativistic, Vlasov plasma. The beam is relatively hot, isotropic in its own frame, and streams parallel or antiparallel to a magnetic field B. Numerical solutions of the full dispersion equation for propagation parallel or antiparallel to B are presented, and the linear properties of the whistler heat flux and electron beam firehose instabilities are exhibited and compared. Under a broad range of parmameters the former mode has the lower beam speed threshold, and the larger maximum growth rate. In addition, it is demonstrated that, for a sufficiently large relative beam density, relative beam temperature, and plasma beta the whistler heat flux instability has a much lower beam speed threshold than the electrostatic electron beam instability. The application of these instabilities to first-order Fermi acceleration of electrons at space plasma shocks is discussed.

Gary, S. P.↗

Cloud Thickness from Diffusion of Lidar Pulses in Clouds

Measurements of the distribution of reflected light from a laser beam incident on an aqueous suspension of particles or "cloud" with known thickness and particle size distribution are reported. The distribution is referred to as the "cloud radiative Green's function", G. In the diffusion domain, G is sensitive to cloud thickness, allowing that important quantity to be retrieved. The goal of the laboratory simulation is to provide preliminary estimates of sensitivity of G to cloud thickness,for use in the optimal design of an offbeam Lidar instrument for remote sensing of cloud thickness (THOR, Thickness from Offbeam Returns). These clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The microsphere size distribution is roughly lognormal, from 0.5 microns to 25 microns, similar to real clouds. Density of suspended spheres is adjusted so mean-free-path of visible photons is about 10 cm, approximately 1000 times smaller than in real clouds. The light source is a ND:YAG laser at 530 nm. Detectors are flux and photon-counting Photomultiplier Tube (PMTS), with a glass probe for precise positioning. A Labview 5 VI controls positioning, and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider, and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns, and step size is selectable from the VI software. Far from the incident beam, the rate of exponential increase as the direction of the incident beam is approached scales as expected from diffusion theory, linearly with the cloud thickness, and inversely as the square root of the reduced optical thickness, and is independent of particle size. Near the beam the signal begins to increase faster than exponential, due to single and low-order scattering near the backward direction, and here the distribution depends on particle size. Results are being used to verify 3D Monte Carlo radiative transfer simulations, used to estimate signal-to-noise ratios for remotely sensed off beam returns, for both homogeneous and inhomogeneous clouds. Signal-to-noise estimates show that unfiltered observations are straight forward at night, while narrow band pass filters are being studied for day.

Cahalan, Robert F.↗

Using the HHT to Search for Gravitational Waves

Gravitational waves are a consequence of Einstein's theory of general relativity applied to the motion of very dense and massive objects such as black holes and neutron stars. Their detection will reveal a wealth of information about these mysterious objects that cannot be obtained with electromagnetic probes. Two projects are underway to attempt the detection of gravitational waves: NASA's Laser Interferometer Space Antenna (LISA), a space based mission being designed to search for waves from supermassive black holes at the centers of galaxies, and the NSF's Laser Interferometer Gravitational Wave Observatory (LIGO), a ground based facility that is now searching for waves from supernovae. pulsars, and the coalescence of black hole and neutron star systems. Because general relativity is an inherently non-linear theory, many of the predicted source waveforms show strong frequency modulation. In addition, the LIGO and LISA detectors are highly sensitive devices that produce a variety of non-linear transient noise features. Thus the unique capabilities of the HHT. the extraction of intrawave modulation and the characterization of non-linear and non-stationary signals, have a natural application to both signal detection and experimental characterization of the detectors. In this talk I will give an overview of the status of the field. including some of the expected sources of gravitational waves, and I will also describe the LISA and LIGO detectors. Then I will describe some applications of the HHT to waveform detection and detector noise characterization.

Camp, Jordan↗

Line-vortex theory for calculation of supersonic downwash

The perturbation field induced by a line vortex in a supersonic stream and the downwash behind a supersonic lifting surface are examined to establish approximate methods for determining the downwash behind supersonic wings. Lifting-lines methods are presented for calculating supersonic downwash. A bent lifting-line method is proposed for computing the downwash field behind swept wings. When applied to triangular wings with subsonic leading edges, this method gives results that, in general, are in good agreement with the exact linearized solution. An unbent lifting-line method (horseshoe-vortex system) is proposed for unswept wings. This method is applied to determine downwash behind rectangular wings with aspect ratios of 2 and 4. Excellent agreement with exact linearized theory is obtained for both aspect ratios by placing the lifting line at the 1/2-chord point. The use of lifting-lines therefore appears promising for obtaining estimates of the downwash behind supersonic wings.

Mirels, Harold↗

Magnetotail dynamics under isobaric constraints

Using linear theory and nonlinear MHD simulations, we investigate the resistive and ideal MHD stability of two-dimensional plasma configurations under the isobaric constraint dP/dt = 0, which in ideal MHD is equivalent to conserving the pressure function P = P(A), where A denotes the magnetic flux. This constraint is satisfied for incompressible modes, such as Alfven waves, and for systems undergoing energy losses. The linear stability analysis leads to a Schroedinger equation, which can be investigated by standard quantum mechanics procedures. We present an application to a typical stretched magnetotail configuration. For a one-dimensional sheet equilibrium characteristic properties of tearing instability are rediscovered. However, the maximum growth rate scales with the 1/7 power of the resistivity, which implies much faster growth than for the standard tearing mode (assuming that the resistivity is small). The same basic eigen-mode is found also for weakly two-dimensional equilibria, even in the ideal MHD limit. In this case the growth rate scales with the 1/4 power of the normal magnetic field. The results of the linear stability analysis are confirmed qualitatively by nonlinear dynamic MHD simulations. These results suggest the interesting possibility that substorm onset, or the thinning in the late growth phase, is caused by the release of a thermodynamic constraint without the (immediate) necessity of releasing the ideal MHD constraint. In the nonlinear regime the resistive and ideal developments differ in that the ideal mode does not lead to neutral line formation without the further release of the ideal MHD constraint; instead a thin current sheet forms. The isobaric constraint is critically discussed. Under perhaps more realistic adiabatic conditions the ideal mode appears to be stable but could be driven by external perturbations and thus generate the thin current sheet in the late growth phase, before a nonideal instability sets in.

Birn, Joachim↗