Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlocality”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Ground-based microwave observations of ozone in the upper stratosphere and mesosphere

A 9-month-long series of mesurements of ozone in the upper stratosphere and mesosphere is reported. The measurements are presented as monthly averages of profiles in blocks of roughly 20 min local time and as night-to-day ratios. An error analysis predicts accuracies of 5-26% for the monthly profiles and 2.5-9% for the ratios. The data are compared to historical data from Solar Mesosphere Explorer (SME) and limb infrared monitor of the stratosphere (LIMS), and it is shown how to remove the effect of different vertical resolution from the comparisons. The microwave data typically agree to better than 10% with SMF and nighttime LIMS ozone at all altitudes below the 0.1-mbar surface. Comparison of the microwave night-to-day ratio with the corresponding ratio from LIMS suggests that nonlocal thermodynamic equilibrium effects in the LIMS daytime data exceed 10% at all pressures less than or equal to 1 mbar.

Connor, Brian J.↗

Auroral activity associated with Kelvin-Helmholtz instability at the inner edge of the low-latitude boundary layer

Auroral activity occurred in the late afternoon sector (approx. 16 MLT) in the northern hemisphere during the passage at Earth of an interplanetary magnetic cloud on January 14, 1988. The auroral activity consisted of a very dynamic display which was preceded and followed by quiet auroral displays. During the quiet displays, discrete rayed arcs aligned along the geomagnetic L shells were observed. In the active stage, rapidly evolving spiral forms centered on magnetic zenith were evident. The activity persisted for many minutes and was characterized by the absence of directed motion. They were strongly suggestive of intense filaments of upward field-aligned currents embedded in the large-scale region 1 current system. Distortions of the flux ropes as they connect from the equatorial magnetosphere to the ionosphere were witnessed. We assess as possible generating mechanisms three nonlocal sources known to be associated with field-aligned currents. Of these, partial compressions of the magnetosphere due to variations of solar wind dynamic pressure seem an unlikely source. The possibility that the auroral forms are due to reconnection is investigated but is excluded because the active aurora were observed on the closed field line region just equatorward of the convection reversal boundary. To support this conclusion further, we apply recent results on the mapping of ionospheric regions to the equatorial plane based on the Tsyganenko 1989 model (Kaufmann et al., 1993). We find that for comparable magnetic activity the aurora map to the equatorial plane at X(sub GSM) = approx. 3 R(sub E) and approx. 2 R(sub E) inward of the magnetopause, that is, the inner edge of the boundary layer close to dusk. Since the auroral forms are manifestly associated with magnetic field shear, a vortical motion at the equatorial end of the flux rope is indicated, making the Kelvin-Helmholtz instability acting at the inner edge of the low-latitude boundary layer the most probable generating source.

Farrugia, C. J.↗

Sounding rocket observation of a hot atomic oxygen geocorona

A sounding rocket measurement of the ultraviolet, atomic oxygen dayglow reveals an excess of emission compared to standard thermospheric model calculations at exospheric altitudes. We explore two explanations for this discrepancy: a breakdown of the radiative transfer model due to nonlocal thermal equilibrium (non-LTE) conditions above the exobase and a hot atomic oxygen geocorona. In particular, the effects of non-LTE on the 3P2, 1, 0 sublevel populations are modeled, and a hot O component in the upper thermosphere and lower exosphere is added to investigate the effects on the modeled emissions. For both cases, the data are reanalyzed and compared with the results using a standard LTE model. A hot O geocorona having a peak density of 10(exp 6)/cc at 550 km and a temperature of 4000 K is consistent with the data and appears to be the most reasonable explanation of the high-altitude enhanced emissions observed in the data.

Cotton, Daniel M.↗

Localized tearing modes in the magnetotail driven by curvature effects

The stability of collisionless tearing modes is examined in the presence of curvature drift resonances and the trapped particle effects. A kinetic description for both electrons and ions is employed to investigate the stability of a two-dimensional equilibrium model. The main features of the study are to treat the ion dynamics properly by incorporating effects associated with particle trajectories in the tail fields and to include the linear coupling of trapped particle modes. Generalized dispersion relations are derived in several parameter regimes by considering two important sublayers of the reconnecting region. For a typical choice of parameters appropriate to the current sheet region, we demonstrate that localized tearing modes driven by ion curvature drift resonance effects are excited in the current sheet region with growth time of the order of a few seconds. Also, we examine nonlocal characteristics of tearing modes driven by curvature effects and show that modes growing in a fraction of a second arise when mode widths are larger than the current sheet width. Further, we show that trapped particle effects, in an interesting frequency regime, significantly enhance the growth rate of the tearing mode. The relevance of this theory for substorm onset phase and other features of the substorms is briefly discussed.

Sundaram, A. K.↗

Kinetic simulation of the Kelvin-Helmholtz instability in a finite sized jet

Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the Kelvin-Helmholtz (K-H) instability for a jet of material moving transverse to the magnetic field in the middle of a field reversal region. The instability shows two phases. This first phase is the development of a short scale length instability on the sides of the jet, which might be characterized as a local K-H instability. This instability saturates at low amplitude by modest broadening of the jet profile. Then a nonlocal K-H instability sets in where the jet exhibits a kinklike behavior which spreads and disrupts the jet in a dramatic fashion. The simulation results are discussed within the context of possible relevance to the magnetopause.

Thomas, V. A.↗

Non-local transport in turbulent MHD convection

The nonlocal non-diffusive transport of passive scalars in turbulent magnetohydrodynamic (MHD) convection is investigated using transilient matrices. These matrices describe the probability that a tracer particle beginning at one position in a flow will be advected to another position after some time. A method for the calculation of these matrices from simulation data which involves following the trajectories of passive tracer particles and calculating their transport statistics, is presented. The method is applied to study the transport in several simulations of turbulent, rotating, three dimensional compressible, penetrative MDH convection. Transport coefficients and other diagnostics are used to quantify the transport, which is found to resemble advection more closely than diffusion. Some of the results are found to have direct relevance to other physical problems, such as the light element depletion in sun-type stars. The large kurtosis found for downward moving particles at the base of the convection zone implies several extreme events.

Miesch, Mark↗

Multi-Particle Interferometry Based on Double Entangled States

A method for producing a 4-photon entangled state based on the use of two independent pair sources is discussed. Of particular interest is that each of the pair sources produces a two-photon state which is simultaneously entangled in both polarization and space-time variables. Performing certain measurements which exploit this double entanglement provides an opportunity for verifying the recent demonstration of nonlocality by Greenberger, Horne, and Zeilinger.

Pittman, Todd B.↗

Two-Photon Entanglement and EPR Experiments Using Type-2 Spontaneous Parametric Down Conversion

Simultaneous entanglement in spin and space-time of a two-photon quantum state generated in type-2 spontaneous parametric down-conversion is demonstrated by the observation of quantum interference with 98% visibility in a simple beam-splitter (Hanburry Brown-Twiss) anticorrelation experiment. The nonlocal cancellation of two-photon probability amplitudes as a result of this double entanglement allows us to demonstrate two different types of Bell's inequality violations in one experimental setup.

Sergienko, A. V.↗

Analysis of Subgrid Eddy Diffusivity in Turbulence Using Direct Numerical Simulation Databases

Direct numerical simulation data in isotropic turbulence at Taylor-scale Reynolds number 90 are used to compute subgrid transfer and spectral eddy diffusivities for passive scalars at Schmidt numbers 1 and 1/8. The results are consistent with spectral transfer processes dominated by a local forward cascade in wavenumber space which is a consequence of moderately nonlocal interactions in the scalar fields. Data decompositions show that resolvable-subgrid interactions make (in general) substantially stronger contributions than fully-subgrid interactions to the eddy diffusivity, which has a non-negligible inverse cascade, or backscatter, component. The eddy diffusivity has a strong cusp at wavenumbers close to the cutoff value separating resolvable and subgrid ranges. At low wavenumbers the eddy diffusivity is generally small, but there the details depend on Schmidt number. The behavior of the subgrid turbulent Schmidt number (the ratio of effective eddy viscosity to eddy diffusivity) is observed to depend strongly on Schmidt number. In interpreting the results at a given spectral cutoff it is important to note that the spectrum of less diffusive scalars (of higher Schmidt number) has more high-wavenumber content.

Yeung, P. K.↗

Artificial Boundary Conditions Based on the Difference Potentials Method

While numerically solving a problem initially formulated on an unbounded domain, one typically truncates this domain, which necessitates setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The issue of setting the ABC's appears to be most significant in many areas of scientific computing, for example, in problems originating from acoustics, electrodynamics, solid mechanics, and fluid dynamics. In particular, in computational fluid dynamics (where external problems present a wide class of practically important formulations) the proper treatment of external boundaries may have a profound impact on the overall quality and performance of numerical algorithms. Most of the currently used techniques for setting the ABC's can basically be classified into two groups. The methods from the first group (global ABC's) usually provide high accuracy and robustness of the numerical procedure but often appear to be fairly cumbersome and (computationally) expensive. The methods from the second group (local ABC's) are, as a rule, algorithmically simple, numerically cheap, and geometrically universal; however, they usually lack accuracy of computations. In this paper we first present a survey and provide a comparative assessment of different existing methods for constructing the ABC's. Then, we describe a relatively new ABC's technique of ours and review the corresponding results. This new technique, in our opinion, is currently one of the most promising in the field. It enables one to construct such ABC's that combine the advantages relevant to the two aforementioned classes of existing methods. Our approach is based on application of the difference potentials method attributable to V. S. Ryaben'kii. This approach allows us to obtain highly accurate ABC's in the form of certain (nonlocal) boundary operator equations. The operators involved are analogous to the pseudodifferential boundary projections first introduced by A. P. Calderon and then also studied by R. T. Seeley. The apparatus of the boundary pseudodifferential equations, which has formerly been used mostly in the qualitative theory of integral equations and PDE'S, is now effectively employed for developing numerical methods in the different fields of scientific computing.

Tsynkov, Semyon V.↗

External Boundary Conditions for Three-Dimensional Problems of Computational Aerodynamics

We consider an unbounded steady-state flow of viscous fluid over a three-dimensional finite body or configuration of bodies. For the purpose of solving this flow problem numerically, we discretize the governing equations (Navier-Stokes) on a finite-difference grid. The grid obviously cannot stretch from the body up to infinity, because the number of the discrete variables in that case would not be finite. Therefore, prior to the discretization we truncate the original unbounded flow domain by introducing some artificial computational boundary at a finite distance of the body. Typically, the artificial boundary is introduced in a natural way as the external boundary of the domain covered by the grid. The flow problem formulated only on the finite computational domain rather than on the original infinite domain is clearly subdefinite unless some artificial boundary conditions (ABC's) are specified at the external computational boundary. Similarly, the discretized flow problem is subdefinite (i.e., lacks equations with respect to unknowns) unless a special closing procedure is implemented at this artificial boundary. The closing procedure in the discrete case is called the ABC's as well. In this paper, we present an innovative approach to constructing highly accurate ABC's for three-dimensional flow computations. The approach extends our previous technique developed for the two-dimensional case; it employs the finite-difference counterparts to Calderon's pseudodifferential boundary projections calculated in the framework of the difference potentials method (DPM) by Ryaben'kii. The resulting ABC's appear spatially nonlocal but particularly easy to implement along with the existing solvers. The new boundary conditions have been successfully combined with the NASA-developed production code TLNS3D and used for the analysis of wing-shaped configurations in subsonic (including incompressible limit) and transonic flow regimes. As demonstrated by the computational experiments and comparisons with the standard (local) methods, the DPM-based ABC's allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable increase of the convergence rate of multigrid iterations.

Tsynkov, Semyon V.↗

An Application of the Difference Potentials Method to Solving External Problems in CFD

Numerical solution of infinite-domain boundary-value problems requires some special techniques that would make the problem available for treatment on the computer. Indeed, the problem must be discretized in a way that the computer operates with only finite amount of information. Therefore, the original infinite-domain formulation must be altered and/or augmented so that on one hand the solution is not changed (or changed slightly) and on the other hand the finite discrete formulation becomes available. One widely used approach to constructing such discretizations consists of truncating the unbounded original domain and then setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The role of the ABC's is to close the truncated problem and at the same time to ensure that the solution found inside the finite computational domain would be maximally close to (in the ideal case, exactly the same as) the corresponding fragment of the original infinite-domain solution. Let us emphasize that the proper treatment of artificial boundaries may have a profound impact on the overall quality and performance of numerical algorithms. The latter statement is corroborated by the numerous computational experiments and especially concerns the area of CFD, in which external problems present a wide class of practically important formulations. In this paper, we review some work that has been done over the recent years on constructing highly accurate nonlocal ABC's for calculation of compressible external flows. The approach is based on implementation of the generalized potentials and pseudodifferential boundary projection operators analogous to those proposed first by Calderon. The difference potentials method (DPM) by Ryaben'kii is used for the effective computation of the generalized potentials and projections. The resulting ABC's clearly outperform the existing methods from the standpoints of accuracy and robustness, in many cases noticeably speed up the multigrid convergence, and at the same time are quite comparable to other methods from the standpoints of geometric universality and simplicity of implementation.

Ryaben 'Kii, Victor S.↗

Cross-field Current Instability for Substorm Expansions

The funding provided by the above-referenced NASA grant has enabled us: (1) to investigate the quasi-linear evolution of the IWI [Lui et al., 1993] and that of the generalized MTSI/IWI [Yoon and Lui, 1993], (2) to carry out the linear analysis of the LHDI to elucidate the difference between it and the MTSI/PM instability [Yoon et al., 1994], (3) to conduct some preliminary nonlocal analyses of the MTSI [Lui et al., 1995] and the IWI [Yoon and Lui, 1996] modes, (4) to study low-frequency shear-driven instability and its nonlinear evolution, which might compete with the CCI [Yoon et al., 1996], and (5) to study the evolution of current sheet during late substorm growth phase by means of 2-D Hall-MHD simulation in order to obtain a better understanding of the current sheet equilibrium crucial for CCI theory [Yoon and Lui, 1997].

Lui, Anthony↗

Physics of Boundaries and Their Interactions in Space Plasmas

This report describes the work done by SciberNet, Inc. during the month of January. During this time, we primarily worked on further analysis of the results presented at the AGU as well as writing them up for publication. Using large scale simulations, we showed that the magnetopause during the southward IMF case is quite irregular with varying thickness, and has a complex flow pattern owing to the nonlinear effects of the convective flow superimposed on the flows generated in the reconnection layer. We used inflow-outflow boundary conditions to examine the kinetic nature of the discontinuities that are formed in the reconnection layer and concluded that nonlocal effects play a major role in the formation of such discontinuities and can alter their properties from the usual structures expected from 1-D simulations or from fluid theories. Finally, we used our 3-D simulations to examine the nonlinear interaction of the tearing mode with the Kelvin-Helmholtz instability. We showed that this interaction leads to the generation of a large core field which is observed both in the magnetotail as well as the magnetopause.

Omidi, Nojan↗

Long-Time Numerical Integration of the Three-Dimensional Wave Equation in the Vicinity of a Moving Source

We propose a family of algorithms for solving numerically a Cauchy problem for the three-dimensional wave equation. The sources that drive the equation (i.e., the right-hand side) are compactly supported in space for any given time; they, however, may actually move in space with a subsonic speed. The solution is calculated inside a finite domain (e.g., sphere) that also moves with a subsonic speed and always contains the support of the right-hand side. The algorithms employ a standard consistent and stable explicit finite-difference scheme for the wave equation. They allow one to calculate tile solution for arbitrarily long time intervals without error accumulation and with the fixed non-growing amount of tile CPU time and memory required for advancing one time step. The algorithms are inherently three-dimensional; they rely on the presence of lacunae in the solutions of the wave equation in oddly dimensional spaces. The methodology presented in the paper is, in fact, a building block for constructing the nonlocal highly accurate unsteady artificial boundary conditions to be used for the numerical simulation of waves propagating with finite speed over unbounded domains.

Ryabenkii, V. S.↗

A New Normalized Difference Cloud Retrieval Technique Applied to Landsat Radiances Over the Oklahoma ARM Site

We suggest a new approach to cloud retrieval, using a normalized difference of nadir reflectivities (NDNR) constructed from a non-absorbing and absorbing (with respect to liquid water) wavelength. Using Monte Carlo simulations we show that this quantity has the potential of removing first order scattering effects caused by cloud side illumination and shadowing at oblique Sun angles. Application of the technique to TM (Thematic Mapper) radiance observations from Landsat-5 over the Southern Great Plains site of the ARM (Atmospheric Radiation Measurement) program gives very similar regional statistics and histograms, but significant differences at the pixel level. NDNR can be also combined with the inverse NIPA (Nonlocal Independent Pixel Approximation) of Marshak (1998) which is applied for the first time on overcast Landsat scene subscenes. We demonstrate the sensitivity of the NIPA-retrieved cloud fields on the parameters of the method and discuss practical issues related to the optimal choice of these parameters.

Orepoulos, Lazaros↗

Quantum Interferometry

Recently, several researchers, including yours truly, have been able to demonstrate theoretically that quantum photon entanglement has the potential to also revolutionize the entire field of optical interferometry, by providing many orders of magnitude improvement in interferometer sensitivity. The quantum entangled photon interferometer approach is very general and applies to many types of interferometers. In particular, without nonlocal entanglement, a generic classical interferometer has a statistical-sampling shot-noise limited sensitivity that scales like 1/Sqrt[N], where N is the number of particles (photons, electrons, atoms, neutrons) passing through the interferometer per unit time. However, if carefully prepared quantum correlations are engineered between the particles, then the interferometer sensitivity improves by a factor of Sqrt[N] (square root of N) to scale like 1/N, which is the limit imposed by the Heisenberg Uncertainty Principle. For optical (laser) interferometers operating at milliwatts of optical power, this quantum sensitivity boost corresponds to an eight-order-of-magnitude improvement of signal to noise. Applications are to tests of General Relativity such as ground and orbiting optical interferometers for gravity wave detection, Laser Interferometer Gravity Observatory (LIGO) and the European Laser Interferometer Space Antenna (LISA), respectively.

Dowling, Jonathan P.↗

Convection in Slab and Spheroidal Geometries

Three-dimensional numerical simulations of compressible turbulent thermally driven convection, in both slab and spheroidal geometries, are reviewed and analyzed in terms of velocity spectra and mixing-length theory. The same ideal gas model is used in both geometries, and resulting flows are compared. The piecewise-parabolic method (PPM), with either thermal conductivity or photospheric boundary conditions, is used to solve the fluid equations of motion. Fluid motions in both geometries exhibit a Kolmogorov-like k(sup -5/3) range in their velocity spectra. The longest wavelength modes are energetically dominant in both geometries, typically leading to one convection cell dominating the flow. In spheroidal geometry, a dipolar flow dominates the largest scale convective motions. Downflows are intensely turbulent and up drafts are relatively laminar in both geometries. In slab geometry, correlations between temperature and velocity fluctuations, which lead to the enthalpy flux, are fairly independent of depth. In spheroidal geometry this same correlation increases linearly with radius over the inner 70 percent by radius, in which the local pressure scale heights are a sizable fraction of the radius. The effects from the impenetrable boundary conditions in the slab geometry models are confused with the effects from non-local convection. In spheroidal geometry nonlocal effects, due to coherent plumes, are seen as far as several pressure scale heights from the lower boundary and are clearly distinguishable from boundary effects.

Porter, David H.↗