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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Realistic Solar Surface Convection Simulations

We perform essentially parameter free simulations with realistic physics of convection near the solar surface. We summarize the physics that is included and compare the simulation results with observations. Excellent agreement is obtained for the depth of the convection zone, the p-mode frequencies, the p-mode excitation rate, the distribution of the emergent continuum intensity, and the profiles of weak photospheric lines. We describe how solar convection is nonlocal. It is driven from a thin surface thermal boundary layer where radiative cooling produces low entropy gas which forms the cores of the downdrafts in which most of the buoyancy work occurs. We show that turbulence and vorticity are mostly confined to the intergranular lanes and underlying downdrafts. Finally, we illustrate our current work on magneto-convection.

Stein, Robert F.↗
NASA NTRS2000

Optical Control of Intersubband Absorption in a Multiple Quantum Well-Embedded Semiconductor Microcravity

Optical intersubband response of a multiple quantum well (MQW)-embedded microcavity driven by a coherent pump field is studied theoretically. The n-type doped MQW structure with three subbands in the conduction band is sandwiched between a semi-infinite medium and a distributed Bragg reflector (DBR). A strong pump field couples the two upper subbands and a weak field probes the two lower subbands. To describe the optical response of the MQW-embedded microcavity, we adopt a semi-classical nonlocal response theory. Taking into account the pump-probe interaction, we derive the probe-induced current density associated with intersubband transitions from the single-particle density-matrix formalism. By incorporating the current density into the Maxwell equation, we solve the probe local field exactly by means of Green's function technique and the transfer-matrix method. We obtain an exact expression for the probe absorption coefficient of the microcavity. For a GaAs/Al(sub x)Ga(sub 1-x)As MQW structure sandwiched between a GaAs/AlAs DBR and vacuum, we performed numerical calculations of the probe absorption spectra for different parameters such as pump intensity, pump detuning, and cavity length. We find that the probe spectrum is strongly dependent on these parameters. In particular, we find that the combination of the cavity effect and the Autler-Townes effect results in a triplet in the optical spectrum of the MQW system. The optical absorption peak value and its location can be feasibly controlled by varying the pump intensity and detuning.

Liu, Ansheng↗
NASA NTRS2000

Incremental Parallelization of Non-Data-Parallel Programs Using the Charon Message-Passing Library

Message passing is among the most popular techniques for parallelizing scientific programs on distributed-memory architectures. The reasons for its success are wide availability (MPI), efficiency, and full tuning control provided to the programmer. A major drawback, however, is that incremental parallelization, as offered by compiler directives, is not generally possible, because all data structures have to be changed throughout the program simultaneously. Charon remedies this situation through mappings between distributed and non-distributed data. It allows breaking up the parallelization into small steps, guaranteeing correctness at every stage. Several tools are available to help convert legacy codes into high-performance message-passing programs. They usually target data-parallel applications, whose loops carrying most of the work can be distributed among all processors without much dependency analysis. Others do a full dependency analysis and then convert the code virtually automatically. Even more toolkits are available that aid construction from scratch of message passing programs. None, however, allows piecemeal translation of codes with complex data dependencies (i.e. non-data-parallel programs) into message passing codes. The Charon library (available in both C and Fortran) provides incremental parallelization capabilities by linking legacy code arrays with distributed arrays. During the conversion process, non-distributed and distributed arrays exist side by side, and simple mapping functions allow the programmer to switch between the two in any location in the program. Charon also provides wrapper functions that leave the structure of the legacy code intact, but that allow execution on truly distributed data. Finally, the library provides a rich set of communication functions that support virtually all patterns of remote data demands in realistic structured grid scientific programs, including transposition, nearest-neighbor communication, pipelining, gather/scatter, and redistribution. At the end of the conversion process most intermediate Charon function calls will have been removed, the non-distributed arrays will have been deleted, and virtually the only remaining Charon functions calls are the high-level, highly optimized communications. Distribution of the data is under complete control of the programmer, although a wide range of useful distributions is easily available through predefined functions. A crucial aspect of the library is that it does not allocate space for distributed arrays, but accepts programmer-specified memory. This has two major consequences. First, codes parallelized using Charon do not suffer from encapsulation; user data is always directly accessible. This provides high efficiency, and also retains the possibility of using message passing directly for highly irregular communications. Second, non-distributed arrays can be interpreted as (trivial) distributions in the Charon sense, which allows them to be mapped to truly distributed arrays, and vice versa. This is the mechanism that enables incremental parallelization. In this paper we provide a brief introduction of the library and then focus on the actual steps in the parallelization process, using some representative examples from, among others, the NAS Parallel Benchmarks. We show how a complicated two-dimensional pipeline-the prototypical non-data-parallel algorithm- can be constructed with ease. To demonstrate the flexibility of the library, we give examples of the stepwise, efficient parallel implementation of nonlocal boundary conditions common in aircraft simulations, as well as the construction of the sequence of grids required for multigrid.

VanderWijngaart, Rob F.↗
NASA NTRS2001

Global Discrete Artificial Boundary Conditions for Time-Dependent Wave Propagation

We construct global artificial boundary conditions (ABCs) for the numerical simulation of wave processes on unbounded domains using a special non-deteriorating algorithm that has been developed previously for the long-term computation of wave-radiation solutions. The ABCs are obtained directly for the discrete formulation of the problem; in so doing, neither a rational approximation of 'non-reflecting kernels,' nor discretization of the continuous boundary conditions is required. The extent of temporal nonlocality of the new ABCs appears fixed and limited; in addition, the ABCs can handle artificial boundaries of irregular shape on regular grids with no fitting/adaptation needed and no accuracy loss induced. The non-deteriorating algorithm, which is the core of the new ABCs is inherently three-dimensional, it guarantees temporally uniform grid convergence of the solution driven by a continuously operating source on arbitrarily long time intervals, and provides unimprovable linear computational complexity with respect to the grid dimension. The algorithm is based on the presence of lacunae, i.e., aft fronts of the waves, in wave-type solutions in odd-dimension spaces, It can, in fact, be built as a modification on top of any consistent and stable finite-difference scheme, making its grid convergence uniform in time and at the same time keeping the rate of convergence the same as that of the non-modified scheme. In the paper, we delineate the construction of the global lacunae-based ABCs in the framework of a discretized wave equation. The ABCs are obtained for the most general formulation of the problem that involves radiation of waves by moving sources (e.g., radiation of acoustic waves by a maneuvering aircraft). We also present systematic numerical results that corroborate the theoretical design properties of the ABCs' algorithm.

Ryabenkii, V. S.↗
NASA NTRS1998

Modeling of Branched (L, T and Y) Carbon Nanotubes

Models for connecting two or three carbon nanotubes (CNT) using topological defects (i.e., pentagons and heptagons) are presented for the characterization of experimentally observed L, T and Y CNT junctions. The effects of the separation and orientation of the topological defects on the structures and energetics of these junctions are investigated using the nonlocal density function theory (DFT) and semi-empirical molecular orbital (AM1) calculations, and the Brenner empirical potential molecular mechanics simulations. The potential applications of L, Y and T CNT junctions in nanoelectronic devices are also discussed.

Han, Jie↗
NASA NTRS2003

Transient Interfacial Phenomena in Miscible Polymer Systems (TIPMPS)

Almost one hundred years ago Korteweg published a theory of how stresses could be induced in miscible fluids by concentration gradients, causing phenomena that would appear to be the same as with immiscible fluids. Miscible fluids could manifest a transient or effective interfacial tension (EIT). To this day, there has been no definitive experiment to confirm Korteweg's model but numerous fascinating and suggestive experiments have been reported. The goal of TIPMPS is to answer the question: Can concentration and temperature gradients in miscible materials induce stresses that cause convection? Many polymer processes involving miscible monomer and polymer systems could be affected by fluid flow and so this work could help understand miscible polymer processing, not only in microgravity, but also on earth. Demonstrating the existence of this phenomenon in miscible fluids will open up a new area of study for materials science. The science objectives of TIPMPS are: (1) Determine if convection can be induced by variation of the width of a miscible interface; (2) Determine if convection can be induced by variation of temperature along a miscible interface; (3) Determine if convection can be induced by variation of conversion along a miscible interface An interface between two miscible fluids can best be created via a spatially-selective photopolymerization of dodecyl acrylate with a photoinitiator, which allows the creation of precise and accurate concentration gradients between polymer and monomer. Optical techniques will be used to measure the refractive index variation caused by the resultant temperature and concentration fields. The viscosity of the polymer will be measured from the increase in the fluorescence of pyrene. Because the large concentration and temperature gradients cause buoyancy-driven convection that prevents the observation of the predicted flows, the experiment must be done in microgravity. In this report, we will consider our efforts to estimate the square gradient parameter, k, and our use of the estimates in modeling of the planned TIPMPS experiments. We developed a model consisting of the heat and diffusion equations with convective terms and of the Navier-Stokes equations with an additional volume force written in the form of the Korteweg stresses arising from nonlocal interaction in the fluid. The fluid's viscosity dependence on polymer conversion and temperature was taken from measurements of poly(dodecyl acrylate). Numerical modeling demonstrated that significant flows would arise for conditions corresponding to the planned experiments.

Pojman, John A.↗
NASA NTRS1996

Electrical double layers at the oil/water interface

This review presents the historical development and current status of the theory of the electrical double layer at a liquid/liquid interface. It gives rigorous thermodynamic definitions of all basic concepts related to liquid interfaces and to the electrical double layer. The difference between the surface of a solid electrode and the interface of two immiscible electrolyte solutions (ITIES) is analyzed in connection to their electrical properties. The most important classical relationships for the electrical double layer are presented and critically discussed. The generalized adsorption isotherm is derived. After a short review of the classical Gouy-Chapman and Verwey-Niessen models, more recent developments of the double layer theory are presented. These include effects of variable dielectric permittivity, nonlocal electrostatics, hydration forces, the modified Poisson-Boltzmann equation and the ion-dipole plasma. The relative merits of different theories are estimated by comparing them with computer simulation of the ITIES and electrical double layer. Special attention is given to the structure of ITIES and its variation due to adsorption of ions and amphiphilic molecules.

Non-NASA Center↗
NASA NTRS1997

An Improved Treatment of External Boundary for Three-Dimensional Flow Computations

We present an innovative numerical approach for setting highly accurate nonlocal boundary conditions at the external computational boundaries when calculating three-dimensional compressible viscous flows over finite bodies. The approach is based on application of the difference potentials method by V. S. Ryaben'kii and extends our previous technique developed for the two-dimensional case. The new boundary conditions methodology has been successfully combined with the NASA-developed code TLNS3D and used for the analysis of wing-shaped configurations in subsonic and transonic flow regimes. As demonstrated by the computational experiments, the improved external boundary conditions 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 speedup of convergence of the multigrid iterations.

Tsynkov, Semyon V.↗
NASA NTRS2002

Nonlinear Propagation of Planet-Generated Tidal Waves

The propagation and evolution of planet-generated density waves in protoplanetary disks is considered. The evolution of waves, leading to shock formation and wake dissipation, is followed in the weakly nonlinear regime. The 2001 local approach of Goodman and Rafikov is extended to include the effects of surface density and temperature variations in the disk as well as the disk cylindrical geometry and nonuniform shear. Wave damping due to shocks is demonstrated to be a nonlocal process spanning a significant fraction of the disk. Torques induced by the planet could be significant drivers of disk evolution on timescales of approx. 10(exp 6)-10(exp 7) yr, even in the absence of strong background viscosity. A global prescription for angular momentum deposition is developed that could be incorporated into the study of gap formation in a gaseous disk around the planet.

Rafikov, R. R.↗
NASA NTRS2005

A Quasi-3-D Theory for Impedance Eduction in Uniform Grazing Flow

A 2-D impedance eduction methodology is extended to quasi-3-D sound fields in uniform or shearing mean flow. We introduce a nonlocal, nonreflecting boundary condition to terminate the duct and then educe the impedance by minimizing an objective function. The introduction of a parallel, sparse, equation solver significantly reduces the wall clock time for educing the impedance when compared to that of the sequential band solver used in the 2-D methodology. The accuracy, efficiency, and robustness of the methodology is demonstrated using two examples. In the first example, we show that the method reproduces the known impedance of a ceramic tubular test liner. In the second example, we illustrate that the approach educes the impedance of a four-segment liner where the first, second, and fourth segments consist of a perforated face sheet bonded to honeycomb, and the third segment is a cut from the ceramic tubular test liner. The ability of the method to educe the impedances of multisegmented liners has the potential to significantly reduce the amount of time and cost required to determine the impedance of several uniform liners by allowing them to be placed in series in the test section and to educe the impedance of each segment using a single numerical experiment. Finally, we probe the objective function in great detail and show that it contains a single minimum. Thus, our objective function is ideal for use with local, inexpensive, gradient-based optimizers.

Watson, W. R.↗
NASA NTRS2007

Transport in the Subtropical Lowermost Stratosphere during the Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment

We use in situ measurements of water vapor (H2O), ozone (O3), carbon dioxide (CO2), carbon monoxide (CO), nitric oxide (NO), and total reactive nitrogen (NOy) obtained during the CRYSTAL-FACE campaign in July 2002 to study summertime transport in the subtropical lowermost stratosphere. We use an objective methodology to distinguish the latitudinal origin of the sampled air masses despite the influence of convection, and we calculate backward trajectories to elucidate their recent geographical history. The methodology consists of exploring the statistical behavior of the data by performing multivariate clustering and agglomerative hierarchical clustering calculations and projecting cluster groups onto principal component space to identify air masses of like composition and hence presumed origin. The statistically derived cluster groups are then examined in physical space using tracer-tracer correlation plots. Interpretation of the principal component analysis suggests that the variability in the data is accounted for primarily by the mean age of air in the stratosphere, followed by the age of the convective influence, and last by the extent of convective influence, potentially related to the latitude of convective injection (Dessler and Sherwood, 2004). We find that high-latitude stratospheric air is the dominant source region during the beginning of the campaign while tropical air is the dominant source region during the rest of the campaign. Influence of convection from both local and nonlocal events is frequently observed. The identification of air mass origin is confirmed with backward trajectories, and the behavior of the trajectories is associated with the North American monsoon circulation.

Pittman, Jasna V.↗
NASA NTRS2008

Demonstrations that the Solar Wind Is Not Accelerated by Waves

The present work uses both observations and theoretical considerations to show that hydromagnetic waves cannot produce the acceleration of the fast solar wind and the related heating of the open solar corona. Waves do exist, and can play a role in the differential heating and acceleration of minor ions, but their amplitudes are not sufficient to power the wind, as demonstrated by extrapolation of magnetic spectra from Helios and Ulysses observations. Dissipation mechanisms invoked to circumvent this conclusion cannot be effective for a variety of reasons. In particular, turbulence does not play a strong role in the corona as shown by both observations of coronal striations and theoretical considerations of line-tying to a nonturbulent photosphere, nonlocality of interactions, and the nature of the kinetic dissipation. In the absence of wave heating and acceleration, the chromosphere and transition region become the natural source of open coronal energization. We suggest a variant of the 'velocity filtration' approach in which the emergence and complex churning of the magnetic flux in the chromosphere and transition region continuously and ubiquitously produces the nonthermal distributions required. These particles are then released by magnetic carpet reconnection at a wide range of scales and produce the wind as described in kinetic approaches. Since the carpet reconnection is not the main source of the energization of the plasma, there is no expectation of an observable release of energy in nanoflares.

Roberts, Aaron↗
NASA NTRS2008

Observational and Theoretical Challenges to Wave or Turbulence Accelerations of the Fast Solar Wind

We use both observations and theoretical considerations to show that hydromagnetic waves or turbulence cannot produce the acceleration of the fast solar wind and the related heating of the open solar corona. Waves do exist as shown by Hinode and other observations, and can play a role in the differential heating and acceleration of minor ions but their amplitudes are not sufficient to power the wind, as demonstrated by extrapolation of magnetic spectra from Helios and Ulysses observations. Dissipation mechanisms invoked to circumvent this conclusion cannot be effective for a variety of reasons. In particular, turbulence does not play a strong role in the corona as shown by both eclipse observations of coronal striations and theoretical considerations of line-tying to a nonturbulent photosphere, nonlocality of interactions, and the nature of kinetic dissipation. In the absence of wave heating and acceleration, the chromosphere and transition region become the natural source of open coronal energization. We suggest a variant of the velocity filtration approach in which the emergence and complex churning of the magnetic flux in the chromosphere and transition region continuously and ubiquitously produces the nonthermal distributions required. These particles are then released by magnetic carpet reconnection at a wide range of scales and produce the wind as described in kinetic approaches. Since the carpet reconnection is not the main source of the energization of the plasma, there is no expectation of an observable release of energy in nanoflares.

Roberts, D. Aaron↗
NASA NTRS2010

Excitation of Crossflow Instabilities in a Swept Wing Boundary Layer

The problem of crossflow receptivity is considered in the context of a canonical 3D boundary layer (viz., the swept Hiemenz boundary layer) and a swept airfoil used recently in the SWIFT flight experiment performed at Texas A&M University. First, Hiemenz flow is used to analyze localized receptivity due to a spanwise periodic array of small amplitude roughness elements, with the goal of quantifying the effects of array size and location. Excitation of crossflow modes via nonlocalized but deterministic distribution of surface nonuniformity is also considered and contrasted with roughness induced acoustic excitation of Tollmien-Schlichting waves. Finally, roughness measurements on the SWIFT model are used to model the effects of random, spatially distributed roughness of sufficiently small amplitude with the eventual goal of enabling predictions of initial crossflow disturbance amplitudes as functions of surface roughness parameters.

Carpenter, Mark H.↗
NASA NTRS2010

The Writhe of Helical Structures in the Solar Corona

Context. Helicity is a fundamental property of magnetic fields, conserved in ideal MHD. In flux rope topology, it consists of twist and writhe helicity. Despite the common occurrence of helical structures in the solar atmosphere, little is known about how their shape relates to the writhe, which fraction of helicity is contained in writhe, and how much helicity is exchanged between twist and writhe when they erupt. Aims. Here we perform a quantitative investigation of these questions relevant for coronal flux ropes. Methods. The decomposition of the writhe of a curve into local and nonlocal components greatly facilitates its computation. We use it to study the relation between writhe and projected S shape of helical curves and to measure writhe and twist in numerical simulations of flux rope instabilities. The results are discussed with regard to filament eruptions and coronal mass ejections (CMEs). Results. (1) We demonstrate that the relation between writhe and projected S shape is not unique in principle, but that the ambiguity does not affect low-lying structures, thus supporting the established empirical rule which associates stable forward (reverse) S shaped structures low in the corona with positive (negative) helicity. (2) Kink-unstable erupting flux ropes are found to transform a far smaller fraction of their twist helicity into writhe helicity than often assumed. (3) Confined flux rope eruptions tend to show stronger writhe at low heights than ejective eruptions (CMEs). This argues against suggestions that the writhing facilitates the rise of the rope through the overlying field. (4) Erupting filaments which are S shaped already before the eruption and keep the sign of their axis writhe (which is expected if field of one chirality dominates the source volume of the eruption), must reverse their S shape in the course of the rise. Implications for the occurrence of the helical kink instability in such events are discussed.

WRITHE↗
NASA NTRS2010

Demonstrations that the Solar Wind is Not Accelerated by Waves or Turbulence

The present work uses observations and theoretical considerations to provide both qualitative and quantitative arguments that hydromagnetic waves, whether turbulent or not, cannot produce the acceleration of the fast solar wind and the related heating of the open solar corona. Waves do exist, and can play a role in the differential heating and acceleration of minor ions, but their amplitudes are not sufficient to power the wind, as demonstrated by extrapolation of magnetic spectra from Helios and Ulysses observations. Dissipation mechanisms invoked to circumvent this conclusion cannot be effective for a variety of reasons. In particular, turbulence does not play a strong role in the corona as shown both by observations of coronal striations and other features, and by theoretical considerations of line-tying to a nonturbulent photosphere, nonlocality of interactions, and the nature of the kinetic dissipation. We consider possible "ways out" of the arguments presented, and suggest that in the absence of wave or turbulent heating and acceleration, the chromosphere and transition region become the natural source, if yet unproven, of open coronal energization through the production of nonthermal particle distributions.

Roberts, D. Aaron↗
NASA NTRS2011

Comment on "'Elevated Heat Pump' Hypothesis for the Aerosol-Monsoon Hydroclimate Link: 'Grounded' in Observations?" by S. Nigam and M. Bollasina

In their recent paper, Nigam and Bollasina [2010] (hereinafter NB) claimed to have found observational evidences that are at variance with the elevated heat pump (EHP) hypothesis regarding the possible impacts of absorbing aerosols on the South Asian summer monsoon [Lau et al., 2006; Lau and Kim, 2006]. We found NB's arguments and inferences against the EHP hypothesis flawed, stemming from their own out of context interpretation of the hypothesis. NB argued that the simultaneous negative correlation of aerosol with rainfall, and correlations with other quantities in May, are evidence against the EHP hypothesis. Their argument cannot be justified. First, Lau and Kim [2006] (hereinafter LK06) never stated that the main rainfall response to EHP is in May. Second, the EHP is about responses of the entire Indian monsoon system that are nonlocal in space and time with respect to the aerosol forcing. As shown in Figure 4 of LK06, while the aerosol anomalies are strongest in April-May, the strongest rainfall response is in June-July, with the enhanced rainfall fed by an induced thermally driven circulation which brings additional moisture from the ocean to the Indian subcontinent. Third, the increased rainfall over the Bay of Bengal as shown in Figure 1a of NB and the increased low-level convergence in Figure 1f of NB do not necessarily reflect responses associated with EHP but rather the large ]scale circulation that provides the buildup of the aerosols before the onset of the monsoon rainfall over India. Because aerosol can only accumulate where there is little or no washout by rain, the negative correlation is a necessary condition for increased atmospheric loading of aerosols. For the same reason, the spatial distributions of rainfall and aerosol generally are offset with each other, i.e., high aerosol in regions of low rainfall. This is evident in Figure 1, which shows the climatological mean of the MODIS aerosol optical depth (AOD), and TRMM rainfall over India in May. The maximum AOD is found over the Indo-Gangetic Plain and the desert regions of northwest India and Pakistan

Lau, K. M.↗
NASA NTRS2010

A Meteorological Overview of the TC4 Mission

The TC4 mission in Central America during summer 2007 examined convective transport into the tropical Upper Troposphere/Lower Stratosphere (UTLS) and the evolution of cirrus clouds. The tropical tropopause layer (TTL) circulation is dominated by the Asian monsoon anticyclone and westward winds that stretch from the western Pacific into the Atlantic. During TC4, TTL westward flow over Central America was stronger than normal. Incidence of cold clouds over the Central American region was the third lowest out of 34 years sampled. The major factor was an incipient La Nina, specifically anomalously cold temperatures off the Pacific Coast of South America. Weakness in the low level Caribbean jet caused a shift in the coldest clouds from the Caribbean to the Pacific side of Central America. The character of tropopause temperature variability was that of upward propagating waves generated by local and nonlocal convection. These waves produced tropopause temperature variations of 3 K, with peak-to-peak variations of 8 K. At low levels in Central America, flow from the Sahara desert predominated; further south, the air came from the Amazon region. Convectively influenced air in the upper troposphere came from Central America, the northern Amazon region, the Atlantic ITCZ, and the North American monsoon. In the TTL, Asian and African convection affected the observed air masses. North of 10N in the Central American TTL, African and Asian convection may have contributed as much to the air masses as Central and South American convection. South of 8N, Asian and African convection had far less impact.

Pfister, L.↗