Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed parallelization”

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,675 records · Page 93

A Proposed Model for Protein Crystal Nucleation and Growth

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, based upon fluorescence, light, neutron, and X-ray scattering data, size exclusion chromatography experiments, dialysis kinetics, AFM, and modeling of growth rate data, from this and other laboratories. The first species formed is postulated to be a 'head to side' dimer. Through repeating associations involving the same intermolecular interactions this grows to a 4(sub 3) helix structure, that in turn serves as the basic unit for nucleation and subsequent crystal growth. High salt attenuates surface charges while promoting hydrophobic interactions. Symmetry facilitates subsequent helix-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(sub3)2(sub 1),2(sub 1) unit cell are filled. 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 growth 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. In the general case for proteins the lack of a singularly defined association pathway may lead 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. This model, its basis, how it fits into the currently understood osmotic second virial coefficient approach to crystallization, and what it suggests will be discussed.

Pusey, Marc↗

Particle/Continuum Hybrid Simulation in a Parallel Computing Environment

The objective of this study was to modify an existing parallel particle code based on the direct simulation Monte Carlo (DSMC) method to include a Navier-Stokes (NS) calculation so that a hybrid solution could be developed. In carrying out this work, it was determined that the following five issues had to be addressed before extensive program development of a three dimensional capability was pursued: (1) find a set of one-sided kinetic fluxes that are fully compatible with the DSMC method, (2) develop a finite volume scheme to make use of these one-sided kinetic fluxes, (3) make use of the one-sided kinetic fluxes together with DSMC type boundary conditions at a material surface so that velocity slip and temperature slip arise naturally for near-continuum conditions, (4) find a suitable sampling scheme so that the values of the one-sided fluxes predicted by the NS solution at an interface between the two domains can be converted into the correct distribution of particles to be introduced into the DSMC domain, (5) carry out a suitable number of tests to confirm that the developed concepts are valid, individually and in concert for a hybrid scheme.

Baganoff, Donald↗

Method and Apparatus of Multiplexing and Acquiring Data from Multiple Optical Fibers Using a Single Data Channel of an Optical Frequency-Domain Reflectometry (OFDR) System

A method and system for multiplexing a network of parallel fiber Bragg grating (FBG) sensor-fibers to a single acquisition channel of a closed Michelson interferometer system via a fiber splitter by distinguishing each branch of fiber sensors in the spatial domain. On each branch of the splitter, the fibers have a specific pre-determined length, effectively separating each branch of fiber sensors spatially. In the spatial domain the fiber branches are seen as part of one acquisition channel on the interrogation system. However, the FBG-reference arm beat frequency information for each fiber is retained. Since the beat frequency is generated between the reference arm, the effective fiber length of each successive branch includes the entire length of the preceding branch. The multiple branches are seen as one fiber having three segments where the segments can be resolved. This greatly simplifies optical, electronic and computational complexity, and is especially suited for use in multiplexed or branched OFS networks for SHM of large and/or distributed structures which need a lot of measurement points.

Parker, Jr., Allen R↗

Polyphony: A Workflow Orchestration Framework for Cloud Computing

Cloud Computing has delivered unprecedented compute capacity to NASA missions at affordable rates. Missions like the Mars Exploration Rovers (MER) and Mars Science Lab (MSL) are enjoying the elasticity that enables them to leverage hundreds, if not thousands, or machines for short durations without making any hardware procurements. In this paper, we describe Polyphony, a resilient, scalable, and modular framework that efficiently leverages a large set of computing resources to perform parallel computations. Polyphony can employ resources on the cloud, excess capacity on local machines, as well as spare resources on the supercomputing center, and it enables these resources to work in concert to accomplish a common goal. Polyphony is resilient to node failures, even if they occur in the middle of a transaction. We will conclude with an evaluation of a production-ready application built on top of Polyphony to perform image-processing operations of images from around the solar system, including Mars, Saturn, and Titan.

Space Exploration,↗

Digital Control Technologies for Modular DC-DC Converters

Recent trends in aerospace Power Management and Distribution (PMAD) systems focus on using commercial off-the-shelf (COTS) components as standard building blocks. This move to more modular designs has been driven by a desire to reduce costs and development times, but is also due to the impressive power density and efficiency numbers achieved by today's commercial DC-DC converters. However, the PMAD designer quickly learns of the hidden "costs" of using COTS converters. The most significant cost is the required addition of external input filters to meet strict electromagnetic interference (MIAMI) requirements for space systems. In fact, the high power density numbers achieved by the commercial manufacturers are greatly due to the lack of necessary input filters included in the COTS module. The NASA Glenn Research Center is currently pursuing a digital control technology that addresses this problem with modular DC-DC converters. This paper presents the digital control technologies that have been developed to greatly reduce the input filter requirements for paralleled, modular DC-DC converters. Initial test result show that the input filter's inductor size was reduced by 75 percent, and the capacitor size was reduced by 94 percent while maintaining the same power quality specifications.

Button, Robert M.↗

Creep and stress rupture of oxide dispersion strengthened mechanically alloyed Inconel alloy MA 754

The creep and stress rupture behavior of the mechanically alloyed oxide dispersion strengthened nickel-base alloy MA 754 was studied at 760, 982 and 1093 C. Tensile specimens with a fine, highly elongated grain structure, oriented parallel and perpendicular to the longitudinal grain direction were tested at various stresses in air under constant load. It was found that the apparent stress dependence was large, with power law exponents ranging from 19 to 33 over the temperature range studied. The creep activation energy, after correction for the temperature dependence of the elastic modulus, was close to but slightly larger than the activation energy for self diffusion. Rupture was intergranular and the rupture ductility as measured by percentage elongation was generally low, with values ranging from 0.5 to 16 pct. The creep properties are rationalized by describing the creep rates in terms of an effective stress which is the applied stress minus a resisting stress consistent with the alloy microstructure. Values of the resisting stress obtained through a curve fitting procedure are found to be close to the values of the particle by-pass stress for this oxide dispersion strengthened alloy, as calculated from the measured oxide particle distribution.

Howson, T. E.↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

Predeployment validation of fault-tolerant systems through software-implemented fault insertion

Fault injection-based automated testing (FIAT) environment, which can be used to experimentally characterize and evaluate distributed realtime systems under fault-free and faulted conditions is described. A survey is presented of validation methodologies. The need for fault insertion based on validation methodologies is demonstrated. The origins and models of faults, and motivation for the FIAT concept are reviewed. FIAT employs a validation methodology which builds confidence in the system through first providing a baseline of fault-free performance data and then characterizing the behavior of the system with faults present. Fault insertion is accomplished through software and allows faults or the manifestation of faults to be inserted by either seeding faults into memory or triggering error detection mechanisms. FIAT is capable of emulating a variety of fault-tolerant strategies and architectures, can monitor system activity, and can automatically orchestrate experiments involving insertion of faults. There is a common system interface which allows ease of use to decrease experiment development and run time. Fault models chosen for experiments on FIAT have generated system responses which parallel those observed in real systems under faulty conditions. These capabilities are shown by two example experiments each using a different fault-tolerance strategy.

Czeck, Edward W.↗

Computation of Gust-Cascade Interaction Using the CE/SE Method

The problem 2 in Category 3 of the 4th Computational Aeroacoustic(CAA) Workshop is solved using the space-time conservation element and solution element (CE/SE) method. This problem models rotor-stator interaction in a 2D cascade. It involves complex geometries and flow physics including vortex shedding and acoustic radiation. The parallel version of the 2D nonlinear Euler solver is used with an unstructured triangular mesh to solve this problem. The Giles approach is incorporated with the CE/SE method to handle non-equal pitches of the rotor and stator. Validation on the Giles approach is performed using Problem 3.1 in the 2nd CAA Workshop. The space-time CE/SE method is a finite volume method with second-order accuracy in both space and time. The flux conservation is enforced in both space and time instead of space only. It has low numerical dissipation and dispersion errors. It uses simple non-reflecting boundary conditions and is compatible with unstructured meshes. It is simple, flexible, and generate reasonably accurate solutions. The CE/SE method has been successfully applied to solve numerous practical problems, especially aeroacoustic problems. Some preliminary numerical results of the benchmark problem 3.2 of the 4th CAA Workshop are shown. The steady-state pressure contour is plotted. The mean pressure distribution on the blade surface is compared with Turbo solution showing a good agreement. The sound pressure level versus the rotor harmonic n at the six designated positions on the blade surface, three locations at inlet plane, and three locations at the outlet plane are plotted. It can be seen that the acoustic response exists only at the excitation frequencies (n = 1,2,3). On the blade surface, the acoustic wave at n = 1 is dominant, while at the inlet and outlet planes, the sound pressure level at n = 2 becomes the largest, which is similar to the results presented. The distribution of sound pressure level at different spatial modes along the z- direction is plotted for n = 1,2,3, respectively. It shows that the spatial modes m = -32 and 22 at n = 1 exponentially decay, and the spatial modes m = 10 at n = 2, m = -42 and 12 at n = 3 propagate both upstream and downstream, which agrees with the prediction based on the linearized theory. Some oscillations are observed, which needs to be investigated further. In the final paper, the numerical results will be compared with a frequency-domain solver LINFLUX solution if it is available.

Wang, X.-Y.↗

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗

The Hubble Space Telescope Medium Deep Survey with the Wide Field and Planetary Camera. 1: Methodology and results on the field near 3C 273

We present results from the Medium Deep Survey (MDS), a Key Project using the Hubble Space Telescope (HST). Wide Field Camera (WFC) images of random fields have been taken in 'parallel mode' with an effective resolution of 0.2 sec full width at half maximum (FWHM) in the V(F555W) and I(F785LP) filters. The exposures presented here were targeted on a field away from 3C 273, and resulted in approximately 5 hr integration time in each filter. Detailed morphological structure is seen in galaxy images with total integrated magnitudes down to V approximately = 22.5 and I approximately = 21.5. Parameters are estimated that best fit the observed galaxy images, and 143 objects are identified (including 23 stars) in the field to a fainter limiting magnitude of I approximately = 23.5. We outline the extragalactic goals of the HST Medium Deep Survey, summarize our basic data reduction procedures, and present number (magnitude) counts, a color-magnitude diagram for the field, surface brightness profiles for the brighter galaxies, and best-fit half-light radii for the fainter galaxies as a function of apparent magnitude. A median galaxy half-light radius of 0.4 sec is measured, and the distribution of galaxy sizes versus magnitude is presented. We observe an apparent deficit of galaxies with half-light radii between approximately 0.6 sec and 1.5 sec, with respect to standard no-evolution or mild evolution cosmological models. An apparent excess of compact objects (half-light radii approximately 0.1 sec) is also observed with respect to those models. Finally, we find a small excess in the number of faint galaxy pairs and groups with respect to a random low-redshift field sample.

Griffiths, R. E.↗

Miniature Sapphire Acoustic Resonator - MSAR

A room temperature sapphire acoustics resonator incorporated into an oscillator represents a possible opportunity to improve on quartz ultrastable oscillator (USO) performance, which has been a staple for NASA missions since the inception of spaceflight. Where quartz technology is very mature and shows a performance improvement of perhaps 1 dB/decade, these sapphire acoustic resonators when integrated with matured quartz electronics could achieve a frequency stability improvement of 10 dB or more. As quartz oscillators are an essential element of nearly all types of frequency standards and reference systems, the success of MSAR would advance the development of frequency standards and systems for both groundbased and flight-based projects. Current quartz oscillator technology is limited by quartz mechanical Q. With a possible improvement of more than x 10 Q with sapphire acoustic modes, the stability limit of current quartz oscillators may be improved tenfold, to 10(exp -14) at 1 second. The electromagnetic modes of sapphire that were previously developed at JPL require cryogenic temperatures to achieve the high Q levels needed to achieve this stability level. However sapphire fs acoustic modes, which have not been used before in a high-stability oscillator, indicate the required Q values (as high as Q = 10(exp 8)) may be achieved at room temperature in the kHz range. Even though sapphire is not piezoelectric, such a high Q should allow electrostatic excitation of the acoustic modes with a combination of DC and AC voltages across a small sapphire disk (approximately equal to l mm thick). The first evaluations under this task will test predictions of an estimated input impedance of 10 kilohms at Q = 10(exp 8), and explore the Q values that can be realized in a smaller resonator, which has not been previously tested for acoustic modes. This initial Q measurement and excitation demonstration can be viewed similar to a transducer converting electrical energy to mechanical energy and back. Such an electrostatic tweeter type excitation of a mechanical resonator will be tested at 5 MHz. Finite element calculation will be applied to resonator design for the desired resonator frequency and optimum configuration. The experiment consists of the sapphire resonator sandwiched between parallel electrodes. A DC+AC voltage can be applied to generate a force to act on a sapphire resonator. With the frequency of the AC voltage tuned to the sapphire resonator frequency, a resonant condition occurs and the sapphire Q can be measured with a high-frequency impedance analyzer. To achieve high Q values, many experimental factors such as vacuum seal, gas damping effects, charge buildup on the sapphire surface, heat dissipation, sapphire anchoring, and the sapphire mounting configuration will need attention. The effects of these parameters will be calculated and folded into the resonator design. It is envisioned that the initial test configuration would allow for movable electrodes to check gap spacing dependency and verify the input impedance prediction. Quartz oscillators are key components in nearly all ground- and space-based communication, tracking, and radio science applications. They play a key role as local oscillators for atomic frequency standards and serve as flywheel oscillators or to improve phase noise in high performance frequency and timing distribution systems. With ultra-stable performance from one to three seconds, an Earth-orbit or moon-based MSAR can enhance available performance options for spacecraft due to elimination of atmospheric path degradation.

Wang, Rabi T.↗

An experimental methodology for a fuzzy set preference model

A flexible fuzzy set preference model first requires approximate methodologies for implementation. Fuzzy sets must be defined for each individual consumer using computer software, requiring a minimum of time and expertise on the part of the consumer. The amount of information needed in defining sets must also be established. The model itself must adapt fully to the subject's choice of attributes (vague or precise), attribute levels, and importance weights. The resulting individual-level model should be fully adapted to each consumer. The methodologies needed to develop this model will be equally useful in a new generation of intelligent systems which interact with ordinary consumers, controlling electronic devices through fuzzy expert systems or making recommendations based on a variety of inputs. The power of personal computers and their acceptance by consumers has yet to be fully utilized to create interactive knowledge systems that fully adapt their function to the user. Understanding individual consumer preferences is critical to the design of new products and the estimation of demand (market share) for existing products, which in turn is an input to management systems concerned with production and distribution. The question of what to make, for whom to make it and how much to make requires an understanding of the customer's preferences and the trade-offs that exist between alternatives. Conjoint analysis is a widely used methodology which de-composes an overall preference for an object into a combination of preferences for its constituent parts (attributes such as taste and price), which are combined using an appropriate combination function. Preferences are often expressed using linguistic terms which cannot be represented in conjoint models. Current models are also not implemented an individual level, making it difficult to reach meaningful conclusions about the cause of an individual's behavior from an aggregate model. The combination of complex aggregate models and vague linguistic preferences has greatly limited the usefulness and predictive validity of existing preference models. A fuzzy set preference model that uses linguistic variables and a fully interactive implementation should be able to simultaneously address these issues and substantially improve the accuracy of demand estimates. The parallel implementation of crisp and fuzzy conjoint models using identical data not only validates the fuzzy set model but also provides an opportunity to assess the impact of fuzzy set definitions and individual attribute choices implemented in the interactive methodology developed in this research. The generalized experimental tools needed for conjoint models can also be applied to many other types of intelligent systems.

Turksen, I. B.↗

Vector constraints on western U.S. deformation from space geodesy, neotectonics, and plate motions

The rate-of-slip vector on the San Andreas fault in central California estimated from geodetic and Holocene geological data is inconsistent with the prediction of rigid plate models such as RM2. This well-known 'San Andreas discrepancy' is diagnostic of plate deformation distributed both east of the fault in the Basin and Range and west of the fault along the California continental margin. Constraints on the integrated deformation rates across these two regions consistent with: (1) the kinematical boundary conditions imposed by the rigid plate model; (2) neotectonic and paleoseismic estimates of deformation rates; (3) ground-based geodetic measurements; and (4) rates of change observed by very long baseline interferometry along seven baselines to western U.S. sites are constructed. The space-geodetic data on Basin and Range extension taken over a 4-year interval are compatible with geological observations averaged over the Holocene; the best estimate of its integrated deformation rate, provided by the joint inversion of both data types, is 9.7 + or - 2.1 mm/yr, N 56 deg W + or - 10 deg, too small and in the wrong direction to account entirely for the San Andreas discrepancy. The integral of this deformation, estimated by subtracting the Basin and Range contribution from the discrepancy vector, requires significant right-lateral shear parallel to the San Andreas (13 + or - 5 mm/yr) and some compression perpendicular to it (9 + or - 3 mm/yr).

Minster, J. Bernard↗

Relationship of Topside Ionospheric Ion Outflows to Auroral Forms and Precipitation, Plasma Waves, and Convection Observed by Polar

The POLAR satellite often observes upflowing ionospheric ions (UFIs) in and near the aurora] oval on southern perigee (approx. 5000 km altitude) passes. We present the UFI features observed by the thermal ion dynamics experiment (TIDE) and the toroidal imaging mass angle spectrograph (TIMAS) in the dusk-dawn sector under two different geomagnetic activity conditions in order to elicit their relationships with auroral forms, wave emissions, and convection pattern from additional POLAR instruments. During the active interval, the ultraviolet imager (UVI) observed a bright discrete aurora on the duskside after the substorm onset and then observed a small isolated aurora form and diffuse auroras on the dawnside during the recovery phase. The UFIs showed clear conic distributions when the plasma wave instrument (PWI) detected strong broadband wave emissions below approx. 10 kHz, while no significant auroral activities were observed by UVI. At higher latitudes, the low-energy UFI conics gradually changed to the polar wind component with decreasing intensity of the broadband emissions. V-shaped auroral kilometric radiation (AKR) signatures observed above -200 kHz by PWI coincided with the region where the discrete aurora and the UFI beams were detected. The latitude of these features was lower than that of the UFI conics. During the observations of the UFI beams and conics, the lower-frequency fluctuations observed by the electric field instrument were also enhanced, and the convection directions exhibited large fluctuations. It is evident that large electrostatic potential drops produced the precipitating electrons and discrete auroras, the UFI beams, and the AKR, which is also supported by the energetic plasma data from HYDRA. Since the intense broadband emissions were also observed with the UFIs, the ionospheric ions could be energized transversely before or during the parallel acceleration due to the potential drops.

Hirahara, M.↗

Carbon Nanotubes in Water: MD Simulations of Internal and External Flow, Self Organization

We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.

Jaffe, Richard L.↗

Relationship of Topside Ionospheric Ion Outflows to Auroral Forms and Precipitations, Plasma Waves, and Convection Observed by POLAR

The POLAR satellite often observes upflowing ionospheric ions (UFls) in and near the auroral oval on southern perigee (approximately 5000 km altitude) passes. We present the UFI features observed by the thermal ion dynamics experiment (TIDE) and the toroidal imaging mass-angle spectrograph (TIMAS) in the dusk-dawn sector under two different geomagnetic activity conditions in order to elicit their relationships with auroral forms, wave emissions, and convection pattern from additional POLAR instruments. During the active interval, the ultraviolet imager (UVI) observed a bright discrete aurora on the dusk side after the substorm onset and then observed a small isolated aurora form and diffuse auroras on the dawn side during the recovery phase. The UFls showed clear conic distributions when the plasma wave instrument (PWI) detected strong broadband wave emissions below approximately 10 kHz, while no significant auroral activities were observed by UVI. At higher latitudes, the low-energy UFI conics gradually changed to the polar wind component with decreasing intensity of the broadband emissions. V-shaped auroral kilometric radiation (AKR) signatures observed above approximately 200 kHz by PWI coincided with the region where the discrete aurora and the UFI beams were detected. The latitude of these features was lower than that of the UFI conics. During the observations of the UFI beams and conics, the lower-frequency fluctuations observed by the electric field instrument (EFI) were also enhanced, and the convection directions exhibited large fluctuations. It is evident that large electrostatic potential drops produced the precipitating electrons and discrete auroras, the UFI beams, and the AKR, which is also supported by the energetic plasma data from HYDRA. Since the intense broadband emissions were also observed with the UFIs. the ionospheric ions could be energized transversely before or during the parallel acceleration due to the potential drops.

Hirahara, M.↗