SEARCH · Engineering Papers

Results for “distributed and parallel processing”

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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Accelerated panel methods using the fast multipole method

Panel methods are commonly used in computational fluid dynamics for the solution of potential flow problems. The methods are a numerical technique based on the surface distribution of singularity elements. The solution is the process of finding the strength of the singularity elements distributed over the body's surface. This process involves the solution of the matrix problem Pq = p' for a set of unknowns q. The Fast Multipole Method is used to directly compute q without using matrix solvers. The algorithm works in O(N) time for N points, a great improvement over standard matrix solvers. In panel methods, the surface of a body is divided into a series of quadrilateral panels. The methods involve the computation of the influence of all other panels on each individual panel. The influence is based on the surface distribution, though this can be approximated by the area for distant panels. An alternative approximation, though with arbitrary accuracy, is to develop a multipole expansion about the center of the panel to describe the effect of a given panel on distant points in space. The expansion is based on the moments of the panel, thus allow the use of various surface distributions without changing the basic algorithm, just the computation of the various moments. The expansions are then manipulated in a tree walk to develop Taylor series expansions about a point in space which describe the effect of all distant panels on any point within a volume of convergence. The effect of near panels then needs to be computed directly, but the effect of all distant panels can be computed by simply evaluating the resulting expansion. The Fast Multipole Method has been applied to panel methods for the solution of source and doublet distributions. A major feature of the algorithm is that the algorithm does not change to derive the potential and velocity for sources and doublets. The same expansions can be used for both sources and doublets. Since the velocity is related to the potential, and the doublet potential is related to the z-component of the source velocity, all values can be derived from the same expansion by taking a series of partial derivatives. This requires more expansion terms to be kept since terms are lost in the process of taking partial derivatives. Thus to maintain accuracy for the doublet computation, more terms are required than if just evaluating for sources. The resulting Fast Multipole code should then parallelize better than classical panel methods due to the locality of data dependencies found in the Fast Multipole Method. Theoretically the parallelized code should execute in O(log N) time with O(N) processors, though this is not practical. Ongoing work includes implementing the parallel accelerated panel method, including methods to improve the load balancing of the problem by taking advantage of the known geometry of panels, and to encorporate sensitivity analysis into the algorithm.

Leathrum, James F., Jr.↗

NASA NTRS2006

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

NASA NTRS1960

The Origin of Tektites

Tektites are probably extraterrestrial, rather than the result of heating some terrestrial materials, because they are a chemically homogeneous group with definite peculiarities (high silica, excess of alkaline earths over alkalis, excess of potash over soda, absence of water), and because some of them (the australites) appear to have undergone ablation in flight through the atmosphere. Since comparatively slow heating is required to explain the liquefaction of the tektite material, it is suggested that the tektites arrived along orbits which were nearly parallel to the surface of the earth, and which resulted from the decay of the orbit of a natural satellite. The great meteor procession of February 9, 1913, is an example of such an object. Comparison with the reentry phenomena of the artificial satellite 1957 Beta suggests that the 1913 shower consisted of a single large stone weighing about 400 kilograms, and a few dozen smaller bodies weighing about 40 grams each, formed by ablation from the larger body. It is shown that under the observed conditions considerable liquid flow would be expected in the stone, which would be heated to about 2100 K. Objects falling from such a shower near the perigee point of the orbit would have a considerable distribution along the orbit as a result of slight variations in height or drag coefficient. The distribution in longitude would be made wider by the turning of the earth under the orbit during the time of fall. The ultimate source of the body which produces a tektite shower is probably the moon, which appears, by virtue of its polarization and the phase distribution of the returned light, to contain high-silica materials. It is suggested that the Igast object alleged to have fallen in 1855 is in fact genuine and represents an unmelted portion of the lunar crust.

OKeefe, J. A.↗

NASA NTRS1995

Runtime volume visualization for parallel CFD

This paper discusses some aspects of design of a data distributed, massively parallel volume rendering library for runtime visualization of parallel computational fluid dynamics simulations in a message-passing environment. Unlike the traditional scheme in which visualization is a postprocessing step, the rendering is done in place on each node processor. Computational scientists who run large-scale simulations on a massively parallel computer can thus perform interactive monitoring of their simulations. The current library provides an interface to handle volume data on rectilinear grids. The same design principles can be generalized to handle other types of grids. For demonstration, we run a parallel Navier-Stokes solver making use of this rendering library on the Intel Paragon XP/S. The interactive visual response achieved is found to be very useful. Performance studies show that the parallel rendering process is scalable with the size of the simulation as well as with the parallel computer.

Ma, Kwan-Liu↗

NASA NTRS1986

A comparison of the regional slope characteristics of Venus and earth - Implications for geologic processes on Venus

The range of 3 degree by 3 degree regional slopes of the Earth and Venus is similar (approximately 0.0-2.4 degrees), although the surface distribution of these values differs significantly. On earth, cratonic and abyssal plains form extensive regions of 0.0 degree slope. Within these regions a variety of features (mid-ocean ridges, volcanic island chains, subduction zones, and floded mountains) have regional slope characteristics influenced by seafloor spreading and plate recycling, as well as an active weathering regime. The plains provinces of Venus are much more rugged than earth's plains and are marked by numerous closely spaced circular and linear features (0.1-0.2 degree regional slope) concentrated into broad linear zones of global extent. Although Venus highlands are bounded by narrow zones of relatively steep slope, the margins of Aphrodite Terra and Beta Regio are not as steep as earth's continental margins and appear to be best developed parallel to the trends of major chasmata within these regions. Ishtar Terra's margins are significantly steeper and more continuous than other highland margins and are comparable to passive margins on earth. The Venus highlands do not contain appreciable smooth, flat interior regions, implying that highland topography is not significantly modified by erosion or deposition.

Sharpton, V. L.↗

NASA NTRS1994

Formation of thin sulfide rinds on the CM parent body

Many of the textures that are observed in CM chondrites have been alternately assigned nebular, parent-body, or combined nebular-parent body origins. Since it is very difficult to substantiate the production of complex textures in the nebula, an alternate approach to this problem is to rigorously determine which of the observable textures could have been produced by reasonable aqueous alteration processes on the CM parent body. Potential parent body reactions involving S deserve special attention because S-bearing phase define many important CM textures, such as rims. We have examined the possibility that the thin (about 5 microns) rinds of sulfides observed around some partially dissolved olivines within the chondrules and matrixes of CM chondrites were formed by the preferential precipitation of sulfides at or near dissolving olivine boundaries during parent-body alteration. Our model defines two infinite and parallel planes of olivine that are separated by pure water. Average separation distances between olivine grains in CM chondrites at the time of accretion (about 100-200 microns) were estimated by assuming a closed system fluid/rock ratio of 45% and varying the bimodal grain-size distribution. We restrict our calculations to the case of an isochemical system with sufficient bisulfide in solution to account for precipitation of pyrrhotite at STP. Our model examines the possibility that dissolving olivines with compositions between Fo(sub 0.9) and Fo(sub 0.4) can produce a strong gradient of Fe(2+) at pHs from 7 to 12 such that the precipitation of pyrrhotite will be initiated at the olivine-fluid boundary. Since CM phase equilibria is consistent with highly reducing conditions, Fe released by olivines would largely remain in solution as Fe(2+) until the solubility product of a Fe(2+) bearing phase was exceeded. Our calculations indicate that all examined combinations of olivine composition and pH yield a strong Fe(2+) gradient at the time the pyrrhotite solubility is exceeded near the olivine-fluid boundary. This favors the initiation of sulfide precipitation around dissolving olivines. However, sulfide rinds are not observed around all partially dissolved olivine grains in CM chondrites, which suggests that the availability of aqueous sulfide species limits the production of the thin sulfide rinds.

Browning, L. B.↗

NASA NTRS2007

Modeling of Radiative Transfer in Protostellar Disks

This program implements a spectral line, radiative transfer tool for interpreting Spitzer Space Telescope observations by matching them with models of protostellar disks for improved understanding of planet and star formation. The Spitzer Space Telescope detects gas phase molecules in the infrared spectra of protostellar disks, with spectral lines carrying information on the chemical composition of the material from which planets form. Input to the software includes chemical models developed at JPL. The products are synthetic images and spectra for comparison with Spitzer measurements. Radiative transfer in a protostellar disk is primarily affected by absorption and emission processes in the dust and in molecular gases such as H2, CO, and HCO. The magnitude of the optical absorption and emission is determined by the population of the electronic, vibrational, and rotational energy levels. The population of the molecular level is in turn determined by the intensity of the radiation field. Therefore, the intensity of the radiation field and the population of the molecular levels are inter-dependent quantities. To meet the computational challenges of solving for the coupled radiation field and electronic level populations in disks having wide ranges of optical depths and spatial scales, the tool runs in parallel on the JPL Dell Cluster supercomputer with C++ and Fortran compiler with a Message Passing Interface. Because this software has been developed on a distributed computing platform, the modeling of systems previously beyond the reach of available computational resources is possible.

VonAllmen, Paul↗

NASA NTRS2002

Nucleation and Growth According to Lysozyme

How does one take a molecule, strongly asymmetric in both shape and charge distribution, and assemble it into a crystal? We propose a model for the nucleation and crystal growth process for tetragonal lysozyme that may be very germane to other monomeric proteins. The first species formed is postulated to be a dimer. Through repeating associations involving the same intermolecular interactions this becomes the 4(sub 3) helix, that in turn serves as the basic unit for nucleation and crystal growth. High salt attenuates surface charges while promoting hydrophobic interactions. Symmetry facilitates helix self-association. Assembly stability is enhanced when a four helix structure is obtained, with each bound to two neighbors. Only two unique interactions are required. The first are those for helix formation, where the dominant interaction is the intermolecular bridging anion. The second is the anti-parallel side-by-side helix-helix interaction, guided by alternating pairs of symmetry related salt bridges along each side. At this stage all eight unique positions of the P4(sub 3)2(sub 1)2(sub 1) unit cell are filled. From the above, the process is one of a) attenuating the most strongly interacting groups, such that b) the molecules begin to self-associate in defined patterns, so that c) symmetry is obtained, which d) propagates as a growing crystal. Simple and conceptually obvious in hindsight, this tells much about what we are empirically doing when we crystallize macromolecules. By adjusting the solution parameters we are empirically balancing the intermolecular interactions, preferentially attenuating the dominant strong (for lysozyme the charged groups) while strengthening the lesser strong (hydrophobic) interactions. Lysozyme is atypical in the breadth of its crystallization conditions; many proteins only crystallize under narrowly defined conditions, pointing to the criticality of the empirical balancing process. Lack of a singularly defined association pathway leads to formation of multiple species, i.e., amorphous precipitation. Weak interactions, such as hydrogen bonds, are promiscuous, serving to strengthen rather than define specific interactions. Participation in an interaction sequesters that surface from subsequent interactions, and we expect the strongest bonds to form first. When two molecules self associate the resulting species will have an axis of symmetry. Subsequent interactions between two associated species having equivalent interactions will also have symmetry. Only a few unique sets of interactions are required to give any of the commonly found space groups for monomeric proteins. This model and what it suggests will be discussed.

Pusey, Marc L.↗

NASA NTRS1995

Monte Carlo simulations of particle acceleration at oblique shocks: Including cross-field diffusion

The Monte Carlo technique of simulating diffusive particle acceleration at shocks has made spectral predictions that compare extremely well with particle distributions observed at the quasi-parallel region of the earth's bow shock. The current extension of this work to compare simulation predictions with particle spectra at oblique interplanetary shocks has required the inclusion of significant cross-field diffusion (strong scattering) in the simulation technique, since oblique shocks are intrinsically inefficient in the limit of weak scattering. In this paper, we present results from the method we have developed for the inclusion of cross-field diffusion in our simulations, namely model predictions of particle spectra downstream of oblique subluminal shocks. While the high-energy spectral index is independent of the shock obliquity and the strength of the scattering, the latter is observed to profoundly influence the efficiency of injection of cosmic rays into the acceleration process.

Baring, M. G.↗

NASA NTRS2004

QuakeSim and the Solid Earth Research Virtual Observatory

We are developing simulation and analysis tools in order to develop a solid Earth science framework for understanding and studying active tectonic and earthquake processes. The goal of QuakeSim and its extension, the Solid Earth Research Virtual Observatory (SERVO), is to study the physics of earthquakes using state-of-the-art modeling, data manipulation, and pattern recognition technologies. We are developing clearly defined accessible data formats and code protocols as inputs to simulations, which are adapted to high-performance computers. The solid Earth system is extremely complex and nonlinear resulting in computationally intensive problems with millions of unknowns. With these tools it will be possible to construct the more complex models and simulations necessary to develop hazard assessment systems critical for reducing future losses from major earthquakes. We are using Web (Grid) service technology to demonstrate the assimilation of multiple distributed data sources (a typical data grid problem) into a major parallel high-performance computing earthquake forecasting code. Such a linkage of Geoinformatics with Geocomplexity demonstrates the value of the Solid Earth Research Virtual Observatory (SERVO) Grid concept, and advances Grid technology by building the first real-time large-scale data assimilation grid.

virtual observatory↗

NASA NTRS2017

AIRS-Only Product in Giovanni for Exploring Up-to-Date AIRS Observation and Comparing with AIRS+AMSU Product

The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) has been the home of processing, archiving, and distribution services for the Atmospheric Infrared Sounder (AIRS) mission since its launch in 2002 for global observations of the atmospheric state. Giovanni, a Web-based application developed by the GES DISC, provides a simple and intuitive way to visualize, analyze, and access vast amounts of Earth science remote sensing data without having to download the data. Most important AIRS variables, including temperature and humidity profiles, outgoing longwave radiation, cloud properties, and trace gases, are available in Giovanni. AIRS is an instrument suite comprised of a hyperspectral infrared instrument (AIRS) and two multichannel microwave instruments, the Advanced Microwave Sounding Unit (AMSU) and the Humidity Sounder for Brazil (HSB). As HSB ceased operation in the very early stages of the AIRS mission, the AIRS project operates two parallel retrieval algorithms: one using both IR and MW measurements (AIRS+AMSU) and the other using only IR measurements (AIRS-only), which covers most of the mission duration. The AIRS+AMSU product is better quality, and the variables in Giovanni are from this product. However, generation of the AIRS+AMSU product has been suspended since the AMSU instrument anomaly occurred in late September 2016. To continue exploring up-to-date AIRS observations, the same set of variables from the AIRS-only product have been added to Giovanni by the GES DSIC. This will also support comparison of AIRS-only with AIRS+AMSU retrievals. In this presentation, we demonstrate the visualization of the AIRS-only product and plots/statistics of comparison with AIRS+AMSU product using Giovanni.

AIRS↗

NASA NTRS1997

Real-Time X-Ray Transmission Microscopy of Solidifying Al-In Alloys

Real-time observations of transparent analog materials have provided insight, yet the results of these observations are not necessarily representative of opaque metallic systems. In order to study the detailed dynamics of the solidification process, we develop the technologies needed for real-time X ray microscopy of solidifying metallic systems, which has not previously been feasible with the necessary resolution, speed, and contrast. In initial studies of Al-In monotectic alloys unidirectionally solidified in an X-ray transparent furnace, in situ records of the evolution of interface morphologies, interfacial solute accumulation, and formation of the monotectic droplets were obtained for the first time: A radiomicrograph of Al-30In grown during aircraft parabolic maneuvers is presented, showing the volumetric phase distribution in this specimen. The benefits of using X-ray microscopy for postsolidification metallography include ease of specimen preparation, increased sensitivity, and three-dimensional analysis of phase distribution. Imaging of the solute boundary layer revealed that the isoconcentration lines are not parallel (as is often assumed) to the growth interface. Striations in the solidified crystal did not accurately decorate the interface position and shape. The monotectic composition alloy under some conditions grew in an uncoupled manner.

Curreri, Peter A.↗

NASA NTRS1995

Growth and microtubule orientation of Zea mays roots subjected to osmotic stress

Previous work has shown that microtubule (MT) reorientation follows the onset of growth inhibition on the lower side of graviresponding roots, indicating that growth reduction can occur independently of MT reorientation. To test this observation further, we examined whether the reduction in growth in response to osmotic stress is correlated with MT reorientation. The distribution and rate of growth in maize roots exposed to 350 mOsm sorbitol and KCl or 5 mM Mes/Tris buffer were measured with a digitizer. After various times roots were processed for indirect immunofluorescence microscopy. Application of sorbitol or KCl had no effect on the organization of MTs in the apical 2 mm of the root but resulted in striking and different effects in the basal region of the root. Sorbitol treatment caused rapid appearance of oval to circular holes in the microtubular array that persisted for at least 9 h. Between 30 min and 4 h of submersion in KCl, MTs in cortical cells 4 mm and farther from the quiescent center began to reorient oblique to the longitudinal axis. After 9 h, the alignment of MTs had shifted to parallel to the root axis but MTs of the epidermal cells remained transverse. In KCl-treated roots MT reorientation appeared to follow a pattern of development similar to that in controls but without elongation. Our data provide additional evidence that MT reorientation is not the cause but a consequence of growth inhibition.

NASA Discipline Plant Biology↗

NASA NTRS2003

Minimizing Segregation during the Controlled Directional Solidification of Dendric Alloys

Gravity-driven convection induced in the liquid by density gradients of temperature or composition disrupts uniform dendritic growth during controlled directional solidification and promotes severe macrosegregation. The solute-rich region about the dendrite tip appears to play a pivotal role in channel initiation. Allen and Hunt referred to this region as an "initial transient" or dynamic region constituting steep concentration gradients. Experimental investigation also point to the role the tip region plays in developing microstructure. Hellawell and co-workers showed that flow-through dendritic channels could be effectively disrupted, and segregation minimized, during the gradient freezing of bulk castings by rotating the melt through a slight angle with respect to Earth's gravity vector. Adapting this principle to controlled directional solidification, it has been shown" that segregation in dendritic alloys can be minimized, and properties improved, by processing the sample near horizontal in conjunction with a slow axial rotation of the crucible. It is postulated that the observed microstructural uniformity arises by maintaining the developing solute field about the dendrite tip. Solute rejected during vertical directional solidification will rise or sink parallel to the primary dendrite arms during axial rotation setting the stage for accumulation, instabilities, and segregation. In contrast, during horizontal growth, the rejected solute will sink or rise perpendicular to the primary dendrite. Now, in the presence of a slight axial rotation, solute that was initially sinking (or rising) will find itself above (or below) its parent dendrite, i.e., still about the tip region. The following is intended to experimentally demonstrate the viability of this concept in coordination with a model that gives predictive insight regarding solute distribution about growing dendrites. Alloys based on the lead-tin eutectic system were used in this study. The system is well characterized, the constituent metals are available in a very pure form, and the thermophysical properties are well known. During solidification of hypoeutectic alloys, e.g., 55 wt pct Pb, the primary dendrites reject the less dense tin, and for the hypereutectic alloys, e.g., 75 wt pct Sn, the primary dendrites reject denser lead. Alloys were prepared by melting appropriate amounts of lead and tin in a glass crucible after which the homogeneous liquid was sucked directly into 5-mm i.d. glass tubes. The sample tube, containing approximately 30 cm of alloy, was then mechanically driven into the directional solidification furnace assembly and positioned such that approx. 20 cm of the sample was remelted. Subsequently, directional solidification was initiated by withdrawing the sample through a water-cooled jacket at a constant growth velocity of 2 ,microns/s. After 5 to 6 cm of growth, the sample was quickly removed from the furnace and quenched in a water bath to preserve the solid-liquid interface. Samples were directionally solidified vertically upward, nearly horizontally, and some in conjunction with an applied axial rotation of the crucible. Temperature gradients at the solid-liquid interface were measured with an in-siru K-type thermocouple. Solidified samples were cut perpendicular and parallel to the growth direction and conventionally prepared for microscopic examination.

Grugel, Richard N.↗

NASA NTRS2010

Fatigue Resistance of the Grain Size Transition Zone in a Dual Microstructure Superalloy Disk

Mechanical property requirements vary with location in nickel-based superalloy disks. To maximize the associated mechanical properties, heat treatment methods have been developed for producing tailored microstructures. In this study, a specialized heat treatment method was applied to produce varying grain microstructures from the bore to the rim portions of a powder metallurgy processed nickel-based superalloy disk. The bore of the contoured disk consisted of fine grains to maximize strength and fatigue resistance at lower temperatures. The rim microstructure of the disk consisted of coarse grains for maximum resistance to creep and dwell crack growth at high temperatures up to 704 C. However, the fatigue resistance of the grain size transition zone was unclear, and needed to be evaluated. This zone was located as a band in the disk web between the bore and rim. Specimens were extracted parallel and transverse to the transition zone, and multiple fatigue tests were performed at 427 and 704 C. Mean fatigue lives were lower at 427 C than for 704 C. Specimen failures often initiated at relatively large grains, which failed on crystallographic facets. Grain size distributions were characterized in the specimens, and related to the grains initiating failures as well as location within the transition zone. Fatigue life decreased with increasing maximum grain size. Correspondingly, mean fatigue resistance of the transition zone was slightly higher than that of the rim, but lower than that of the bore. The scatter in limited tests of replicates was comparable for all transition zone locations examined.

Gabb, T. P.↗

NASA NTRS1984

Interactive computer modeling of combustion chemistry and coalescence-dispersion modeling of turbulent combustion

An interactive computer code for simulation of a high-intensity turbulent combustor as a single point inhomogeneous stirred reactor was developed from an existing batch processing computer code CDPSR. The interactive CDPSR code was used as a guide for interpretation and direction of DOE-sponsored companion experiments utilizing Xenon tracer with optical laser diagnostic techniques to experimentally determine the appropriate mixing frequency, and for validation of CDPSR as a mixing-chemistry model for a laboratory jet-stirred reactor. The coalescence-dispersion model for finite rate mixing was incorporated into an existing interactive code AVCO-MARK I, to enable simulation of a combustor as a modular array of stirred flow and plug flow elements, each having a prescribed finite mixing frequency, or axial distribution of mixing frequency, as appropriate. Further increase the speed and reliability of the batch kinetics integrator code CREKID was increased by rewriting in vectorized form for execution on a vector or parallel processor, and by incorporating numerical techniques which enhance execution speed by permitting specification of a very low accuracy tolerance.

Pratt, D. T.↗

NASA NTRS2012

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗

NASA NTRS1996

Efficient Process Migration for Parallel Processing on Non-Dedicated Networks of Workstations

This paper presents the design and preliminary implementation of MpPVM, a software system that supports process migration for PVM application programs in a non-dedicated heterogeneous computing environment. New concepts of migration point as well as migration point analysis and necessary data analysis are introduced. In MpPVM, process migrations occur only at previously inserted migration points. Migration point analysis determines appropriate locations to insert migration points; whereas, necessary data analysis provides a minimum set of variables to be transferred at each migration pint. A new methodology to perform reliable point-to-point data communications in a migration environment is also discussed. Finally, a preliminary implementation of MpPVM and its experimental results are presented, showing the correctness and promising performance of our process migration mechanism in a scalable non-dedicated heterogeneous computing environment. While MpPVM is developed on top of PVM, the process migration methodology introduced in this study is general and can be applied to any distributed software environment.

Chanchio, Kasidit↗