Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave 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 1,369 records · Page 76

Evolution of interstellar pickup ions in the solar wind

A model is constructed for the evolution of an interstellar pickup ion distribution in the solar wind. The model assumes that the ions are immediately isotropized at ionization and follows the subsequent development of the distribution function as the particles are convected with the solar wind. The effects of energy diffusion in an ambient wave field with a power law spectrum, adiabatic deceleration in the expanding solar wind, and continual addition of newly ionized particles are all included in the model. An analytical expression describing the evolution of the distribution function in phase space velocity and heliocentric radius is obtained. The distribution quickly approaches an asymptotic shape in phase space which depends on the relative efficiency of the energy diffusion process compared to that of adiabatic deceleration. At large distances from the sun the density of pickup ions falls as 1/r in this model. An expression for the distribution function at large distances and for large particle speed is presented. The asymptotic shape should describe the distribution of pickup ions in the outer heliosphere and could be used as an input distribution for a model of the anomalous component of heliosphere and could be used as an input distribution for a model of the anomalous component of cosmic rays. Comparison of this work with the recent observation of He(+) at 1 AU implies that the energy diffusion process is very weak inside 1 AU.

Isenberg, Philip A.↗

Effective propagation constants in dense random media under effective medium approximation

The effective medium approximation is applied to study the effective propagation constants in a dense random medium. The dyadic Green's function is introduced to establish the effective medium approximation formalism for electromagnetic waves. The multiple scattering equations and the Lippmann-Schwinger equations for the transition operator are obtained in the configuration average form. The dispersion equations of multiple scattering is derived by using a standard method in quantum mechanics. To obtain an expression for the effective propagation constants the matrix elements of the configuration average dyadic transiton operator are calculated in momentum representation. Numerical illustations are carried out to demonstrate the difference in the effective propagation constants between the use of this approximation and the well-known quasi-crystalline approximation. A comparison is made with measured loss tangent in dry snow.

Zhu, P. Y.↗

Dynamics of atmospheric teleconnections during the northern summer

The mechanisms of northern summertime teleconnections are investigated using a barotropic model. In a series of numerical experiments the atmospheric response over the eastern Pacific-North America to an idealized local divergence source corresponding to the northward displacement of the ITCZ in the eastern Pacific is studied. It is found that the response is much stronger in June than in May and is strongest when the forcing is located north of about 10 deg N. This can be explained in terms of the refractive properties of the climatological summertime subtropical jet stream over North America. In another series of experiments the global response as a function of the longitudinal location of the tropical forcing is examined. A wave train emanating from the subtropics of the western Pacific near the Philippines, arching across the Aleutians and the Gulf of Alaska, and terminating with a high anomaly over the continental United States appears over a wide longitudinal range of local forcing, suggesting the existence of a normal mode for the northern summertime climatological flow. The normal-mode concept is supported by further experiments using extratropical forcings as well as free-mode integrations.

Lau, K.-M.↗

Core Losses In Three Alloys

Report describes experimental investigation of core losses and dynamic magnetic flux densities as functions of magnetomotive force of three soft magnetic materials under sine-wave voltage excitation at frequencies from 1 to 50 kHz and at temperatures from 23 to 300 degrees C. Prompted by need for data to design magnetic cores for efficient, lightweight power supplies operating in harsh environments.

Schwarze, G. E.↗

The helical decomposition and the instability assumption

Direct numerical simulations show that the triadic transfer function T(k,p,q) peaks sharply when q (or p) is much smaller than k. The triadic transfer function T(k,p,q) gives the rate of energy input into wave number k from all interactions with modes of wave number p and q, where k, p, q form a triangle. This observation was thought to suggest that energy is cascaded downscale through non-local interactions with local transfer and that there was a strong connection between large and small scales. Both suggestions were in contradiction with the classical Kolmogorov picture of the energy cascade. The helical decomposition was found useful in distinguishing between kinematically independent interactions. That analysis has gone beyond the question of non-local interaction with local transfer. In particular, an assumption about the statistical direction of triadic energy transfer in any kinematically independent interaction was introduced (the instability assumption). That assumption is not necessary for the conclusions about non-local interactions with local transfer recalled above. In the case of turbulence under rapid rotation, the instability assumption leads to the prediction that energy is transferred in spectral space from the poles of the rotation axis toward the equator. The instability assumption is thought to be of general validity for any type of triad interactions (e.g. internal waves). The helical decomposition and the instability assumption offer detailed information about the homogeneous statistical dynamics of the Navier-Stokes equations. The objective was to explore the validity of the instability assumption and to study the contributions of the various types of helical interactions to the energy cascade and the subgrid-scale eddy-viscosity. This was done in the context of spectral closures of the Direct Interaction or Quasi-Normal type.

Waleffe, Fabian A.↗

Model images of radio halos around supernova remnants

I present model calculations of profiles and two-dimensional images of the radio synchrotron emission of young supernova remnants, concentrating on observable effects of relativistic eletrons diffusing upstream of the shock wave. If the preshock electron scattering mean free path is sufficiently long, observable synchrotron halos outside the bulk of the radio emission can potentially result; their absence can constrain the mean free path from above. If scattering is primarily due, as expected, to Alfven waves with amplitude detla(B), the halo is expected to extend a distance of order r(sub g)c(delta(B)/B)(exp 2)/v(sub s) beyond the shock, where r(sub g) is the gyroradius of the electrons emitting at the observed frequency, B is the upstream magnetic field strength, v(sub s) is the shock velocity, and the amplitude delta(B) refers to wave with wavelength comparable to r(sub g), of order 10(exp 13) cm for typical supernova-remnant parameters. However, the detailed geometry of the halo varies with the assumptions about particle acceleration in the shock wave. I present an atlas of model profiles and images as a function of preshock diffusion length, of aspect angle between the magnetic field and the line of sight, and of other relevant parameters.

Reynolds, Stephen P.↗

Numerical Analysis of the Sea State Bias for Satellite Altimetry

Theoretical understanding of the dependence of sea state bias (SSB) on wind wave conditions has been achieved only for the case of a unidirectional wind-driven sea. Recent analysis of Geosat and TOPEX altimeter data showed that additional factors, such as swell, ocean currents, and complex directional properties of realistic wave fields, may influence SSB behavior. Here we investigate effects of two-dimensional multimodal wave spectra using a numerical model of radar reflection from a random, non-Gaussian surface. A recently proposed ocean wave spectrum is employed to describe sea surface statistics. The following findings appear to be of particular interest: (1) Sea swell has an appreciable effect in reducing the SSB coefficient compared with the pure wind sea case but has less effect on the actual SSB owing to the corresponding increase in significant wave height. (2) Hidden multimodal structure (the two-dimensional wavenumber spectrum contains separate peaks, for swell and wind seas, while the frequency spectrum looks unimodal) results in an appreciable change of SSB. (3) For unimodal, purely wind-driven seas, the influence of the angular spectral width is relatively unimportant; that is, a unidirectional sea provides a good qualitative model for SSB if the swell is absent. (4) The pseudo wave age is generally much better fo parametrization the SSB coefficient than the actual wave age (which is ill-defined for a multimodal sea) or wind speed. (5) SSB can be as high as 5% of the significant wave height, which is significantly greater than predicted by present empirical model functions tuned on global data sets. (6) Parameterization of SSB in terms of wind speed is likely to lead to errors due to the dependence on the (in practice, unknown) fetch.

Glazman, R. E.↗

Wind Induced Sediment Resuspension in a Microtidal Estuary

Bottom sediment resuspension frequency, duration and extent (% of bottom sediments affected) were characterized for the fifteen month period from September 1995 to January 1997 for the Barataria Basin, LA. An empirical model of sediment resuspension as a function of wind speed, direction, fetch and water depth was derived from wave theory. Water column turbidity was examined by processing remotely sensed radiance information from visible and near-IR AVHRR imagery. Based on model predictions, wind induced resuspension occurred during all seasons of this study. Seasonal characteristics for resuspension reveal that late fall, winter and early spring are the periods of most frequent and intense resuspension. Model predictions of the critical wind speed required to induce resuspension indicate that winds of 4 m/s (averaged over all wind directions resuspend approximately 50% of bottom sediments in the water bodies examined. Winds of this magnitude (4 m/s) occurred for 80% of the time during the late fall, winter and early spring and for approximately 30% of the time during the summer. More than 50% of the bottom sedimets are resuspended throughout the year, indicating the importance of resuspension as a process affecting sediment and biogeochemical fluxes in the Barataria Basin.

Booth, J. G.↗

Aeroelastic System Development Using Proper Orthogonal Decomposition and Volterra Theory

This research combines Volterra theory and proper orthogonal decomposition (POD) into a hybrid methodology for reduced-order modeling of aeroelastic systems. The out-come of the method is a set of linear ordinary differential equations (ODEs) describing the modal amplitudes associated with both the structural modes and the POD basis functions for the uid. For this research, the structural modes are sine waves of varying frequency, and the Volterra-POD approach is applied to the fluid dynamics equations. The structural modes are treated as forcing terms which are impulsed as part of the uid model realization. Using this approach, structural and uid operators are coupled into a single aeroelastic operator. This coupling converts a free boundary uid problem into an initial value problem, while preserving the parameter (or parameters) of interest for sensitivity analysis. The approach is applied to an elastic panel in supersonic cross ow. The hybrid Volterra-POD approach provides a low-order uid model in state-space form. The linear uid model is tightly coupled with a nonlinear panel model using an implicit integration scheme. The resulting aeroelastic model provides correct limit-cycle oscillation prediction over a wide range of panel dynamic pressure values. Time integration of the reduced-order aeroelastic model is four orders of magnitude faster than the high-order solution procedure developed for this research using traditional uid and structural solvers.

Lucia, David J.↗

Delaying Trains of Short Light Pulses in WGM Resonators

Suitably configured whispering-gallery-mode (WGM) optical resonators have been proposed as delay lines for trains of short light pulses. Until now, it has been common practice to implement an optical delay line as a coiled long optical fiber, which is bulky and tends to be noisy. An alternative has been to implement an optical delay line as a coupled-resonator optical waveguide (a chain of coupled optical resonators), which is compact but limits the width of the pulse spectrum to the width of an optical resonance and thereby places a lower limit on the duration of a pulse. In contrast, a delay line according to the proposal could be implemented as a single WGM resonator, and the pulses delayed by the resonator could be so short that their spectral widths could greatly exceed the spectral width of any single resonance. The proposal emerged from theoretical and experimental studies of the propagation of a pulse train in a WGM resonator. An important element of the theoretical study was recognition that the traditional definition of group velocity in effect, the velocity of a single pulse comprising a packet of waves propagating in a medium, the responsivity of which is a monotonous function of frequency does not necessarily apply in the case of a WGM resonator or other medium having a spectrum consisting of discrete resonance peaks at different frequencies. A new definition of group velocity, applicable to a train of pulses propagating in such a medium, was introduced and found to lead to the discovery of previously unknown features of propagation. Notably, it was found that in a micro - sphere optical resonator that supports a suitable combination of WGM modes, the group velocity for a train of light pulses could be positive, zero, or negative. A positive group velocity could be so small that the delay could be much longer than the ring-down time of the resonator; a delay of such great length is impossible for a single pulse interacting with either a linearly responding lossless resonator or a coupled- resonator optical waveguide.

Matsko, Andrey↗

Full-Circle Resolver-to-Linear-Analog Converter

A circuit generates sinusoidal excitation signals for a shaft-angle resolver and, like the arctangent circuit described in the preceding article, generates an analog voltage proportional to the shaft angle. The disadvantages of the circuit described in the preceding article arise from the fact that it must be made from precise analog subcircuits, including a functional block capable of implementing some trigonometric identities; this circuitry tends to be expensive, sensitive to noise, and susceptible to errors caused by temperature-induced drifts and imprecise matching of gains and phases. These disadvantages are overcome by the design of the present circuit. The present circuit (see figure) includes an excitation circuit, which generates signals Ksin(Omega(t)) and Kcos(Omega(t)) [where K is an amplitude, Omega denotes 2(pi)x a carrier frequency (the design value of which is 10 kHz), and t denotes time]. These signals are applied to the excitation terminals of a shaft-angle resolver, causing the resolver to put out signals C sin(Omega(t)-Theta) and C cos(Omega(t)-Theta). The cosine excitation signal and the cosine resolver output signal are processed through inverting comparator circuits, which are configured to function as inverting squarers, to obtain logic-level or square-wave signals .-LL[cos(Omega(t)] and -LL[cos(Omega(t)-Theta)], respectively. These signals are fed as inputs to a block containing digital logic circuits that effectively measure the phase difference (which equals Theta between the two logic-level signals). The output of this block is a pulse-width-modulated signal, PWM(Theta), the time-averaged value of which ranges from 0 to 5 VDC as Theta ranges from .180 to +180deg. PWM(Theta) is fed to a block of amplifying and level-shifting circuitry, which converts the input PWM waveform to an output waveform that switches between precise reference voltage levels of +10 and -10 V. This waveform is processed by a two-pole, low-pass filter, which removes the carrier-frequency component. The final output signal is a DC potential, proportional to Theta that ranges continuously from -10 V at Theta = -180deg to +10 V at Theta = +180deg..

Alhorn, Dean C.↗

An Empirical Temperature Variance Source Model in Heated Jets

An acoustic analogy approach is implemented that models the sources of jet noise in heated jets. The equivalent sources of turbulent mixing noise are recognized as the differences between the fluctuating and Favre-averaged Reynolds stresses and enthalpy fluxes. While in a conventional acoustic analogy only Reynolds stress components are scrutinized for their noise generation properties, it is now accepted that a comprehensive source model should include the additional entropy source term. Following Goldstein s generalized acoustic analogy, the set of Euler equations are divided into two sets of equations that govern a non-radiating base flow plus its residual components. When the base flow is considered as a locally parallel mean flow, the residual equations may be rearranged to form an inhomogeneous third-order wave equation. A general solution is written subsequently using a Green s function method while all non-linear terms are treated as the equivalent sources of aerodynamic sound and are modeled accordingly. In a previous study, a specialized Reynolds-averaged Navier-Stokes (RANS) solver was implemented to compute the variance of thermal fluctuations that determine the enthalpy flux source strength. The main objective here is to present an empirical model capable of providing a reasonable estimate of the stagnation temperature variance in a jet. Such a model is parameterized as a function of the mean stagnation temperature gradient in the jet, and is evaluated using commonly available RANS solvers. The ensuing thermal source distribution is compared with measurements as well as computational result from a dedicated RANS solver that employs an enthalpy variance and dissipation rate model. Turbulent mixing noise predictions are presented for a wide range of jet temperature ratios from 1.0 to 3.20.

Khavaran, Abbas↗

Application of a First-Principles Anomalous Transport Model for Electrons to Multiple Hall Thrusters and Operating Conditions

We have developed a physics-based model based on a pseudo-particle description of the electron cyclotron drift instability. A key improvement of the model with respect to previous work is that linear theory is not applied in the event of wave saturation and deviations of electrons or ions from a Maxwellian distribution function. In the acceleration region, the anomalous collision frequency is computed as the minimum value necessary to prevent the electron drift velocity from exceeding the thermal velocity. A functional based on the electron equilibration time is defined to control the transition from high to low resistivity regions. The model was previously applied to a single Hall thruster at its nominal operating condition, showing promising results that captured accurately the location of the thruster’s acceleration region. In this paper, we extend the use of this firstprinciples models to two additional thrusters, also considering multiple operating conditions for each of them. Numerical results are compared to experimental measurements obtained with non-invasive laser induced fluorescence. In general, the agreement between experiments and simulations is good. The model is able to predict the location of the acceleration region for all cases. We observe however that fine details, such as changes in the plasma potential gradient within the acceleration regions, are not captured. The model is also insensitive to changes in the magnetic field strength while experiments show that small shifts in location (of less than 5% of the acceleration channel length) occur. We plan to address the weaknesses of our method with the help of physical insight gained from kinetic simulations of the acceleration region.

Chaplin, Vernon H.↗

High-Order/Low-Dissipation Chain-Rule Flux Solution Reconstruction Schemes in FUN3D

In this paper, we report progress in the development of economically high-order flux-solution-reconstruction (FSR) schemes, which are second-order accurate on general unstructured grids but achieve high-order accuracy when a grid is regular, i.e., has the same stencil (with the same spacing) throughout the domain. Two variants of the FSR schemes are discussed: chain-rule-flux-solution reconstruction (CFSR) and quadratic-form-flux-solution reconstruction (QFSR), where the former is based on the chain rule and the latter on the flux reconstruction expressed as a function of solution variables. These schemes are tested for flows with shock waves with a limiter incorporated in the flux and solution reconstructions. Improved results, compared with second-order methods are demonstrated for inviscid and viscous flows with smooth grids.

high-order↗

Two-layer viscous instability in a rotating couette device

A novel experiment to study the interfacial shear instability between two liquids is described. Density-matched immiscible liquids are confined between concentric cylinders such that the interface is parallel to the cylinder walls. Interfacial waves that develop because of viscosity differences between the shearing fluids are studied as a function of rotation rate and depth ratio using optical techniques. Conditions neutral stability and the most unstable wavenumber agree reasonably well with predictions from linear stability analysis of the Navier-Stokes equations. Illumination using laser sheets allows precise measurement of the interface shape. Future experiments will verify the correctness of weakly nonlinear theories that describe energy transfer and saturation of wave growth by nonlinear effects. Measurements of solitary wave shapes, that occur far above neutral stability, will be compared to similar measurements for systems that have gravity as an important force to determine how gravity effects large disturbances. These results will be used to interpret slug and annular flow data that have been obtained in other mu g studies.

Gallagher, Christopher↗

Vector formulation for interferogram surface fitting

Interferometry is an optical testing technique that quantifies the optical path difference (OPD) between a reference wave front and a test wave front based on the interference of light. Fringes are formed when the OPD is an integral multiple of the illuminating wavelength. The resultant two-dimensional pattern is called an interferogram. The function of any interferogram analysis program is to extract this OPD and to produce a representation of the test wave front (or surface). This is accomplished through a three-step process of sampling, ordering, and fitting. We develop a generalized linear-algebra vector-notation model of the interferogram sampling and fitting process.

Fischer, David J.↗

SAVS: A Space and Atmospheric Visualization Science system

The research environment faced by space and atmospheric scientists in the 1990s is characterized by unprecedented volumes of new data, by ever-increasing repositories of unexploited mission files, and by the widespread use of empirical and large-scale computational models needed for the synthesis of understanding across data sets and discipline boundaries. The effective analysis and interpretation of such massive amounts of information have become the subjects of legitimate concern. With SAVS (a Space and Atmospheric Visualization Science System), we address these issues by creating a 'push-button' software environment that mimics the logical scientific processes in data acquisition, reduction, and analysis without requiring a detailed understanding of the methods, networks, and modules that link the tools and effectively execute the functions. SAVS provides (1) a customizable framework for accessing a powerful set of visualization tools based on the popular AVS visualization software with hooks to PV-Wave and access to Khoros modules, (2) a set of mathematical and statistical tools, (3) an extensible library of discipline-specific functions and models (e.g., MSIS, IRI, Feldstein Oval, IGRF, satellite tracking with CADRE-3, etc.), and (4) capabilities for local and remote data base access. The system treats scalar, vector, and image data, and runs on most common Unix workstations. We present a description of SAVS and its components, followed by several applications based on generic research interests in interplanetary and magnetospheric physics (IMP/ISTP), active experiments in space (CRRES), and mission planning focused on the Earth's thermospheric, ionospheric, and mesospheric domains (TIMED).

Szuszczewicz, E. P.↗

MODIS On-Board Blackbody Function and Performance

Two MODIS instruments are currently in orbit, making continuous global observations in visible to long-wave infrared wavelengths. Compared to heritage sensors, MODIS was built with an advanced set of on-board calibrators, providing sensor radiometric, spectral, and spatial calibration and characterization during on-orbit operation. For the thermal emissive bands (TEB) with wavelengths from 3.7 m to 14.4 m, a v-grooved blackbody (BB) is used as the primary calibration source. The BB temperature is accurately measured each scan (1.47s) using a set of 12 temperature sensors traceable to NIST temperature standards. The onboard BB is nominally operated at a fixed temperature, 290K for Terra MODIS and 285K for Aqua MODIS, to compute the TEB linear calibration coefficients. Periodically, its temperature is varied from 270K (instrument ambient) to 315K in order to evaluate and update the nonlinear calibration coefficients. This paper describes MODIS on-board BB functions with emphasis on on-orbit operation and performance. It examines the BB temperature uncertainties under different operational conditions and their impact on TEB calibration and data product quality. The temperature uniformity of the BB is also evaluated using TEB detector responses at different operating temperatures. On-orbit results demonstrate excellent short-term and long-term stability for both the Terra and Aqua MODIS on-board BB. The on-orbit BB temperature uncertainty is estimated to be 10mK for Terra MODIS at 290K and 5mK for Aqua MODIS at 285K, thus meeting the TEB design specifications. In addition, there has been no measurable BB temperature drift over the entire mission of both Terra and Aqua MODIS.

Xiaoxiong, Xiong↗