Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dissipative particle dynamics”

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 55 records · Page 3

Orbital resonance in a dissipative medium

Orbital characteristics of solar system bodies are reviewed with reference to their degree of circularity or eccentricity. Attention is given to orbital perturbations in the vicinity of Jupiter, Saturn, and the Asteroid Belt. The trapping of orbital periods is explained on the basis of relative velocity and gravitational interaction between particles. Enhancement of orbital eccentricity through resonance with other bodies is dampened by the drag-effect of small particles in the orbital path. This mechanism is seen operating in interactions between Titan and Hyperion, and between individual asteroids in the vicinity of Jupiter; with implications for the accretion dynamics of a proto-planetary nebula.

Greenberg, R.↗

Time dependent simulation of cosmic-ray shocks including Alfven transport

Time evolution of plane, cosmic-ray modified shocks was simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven wave energy transport equations. Nonlinear feedback from the cosmic-rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the time scale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic-rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗

Effects of gas on the global stability of galactic disks - Radial flows

We study numerically the effect of gas on the global stability of a two-component self-gravitating galactic disk embedded in a live halo. The stars are evolved by using a 3D collisionless N-body code, and the gas is represented by an ensemble of finite size inelastic particles. The gravitational interaction of stars and gas is calculated using a TREE method. We find that the evolution of the gaseous distribution in the globally unstable disks can be described by two different regimes. When the gas mass fraction is less than about 10 percent, the gas is channeled toward the galactic center by a growing stellar bar. For higher gas fractions, the gas becomes highly inhomogeneous, and the bar instability in the disk is heavily damped. The gas falls toward the inner kpc due to dynamical friction. Domains of both regimes depend on the efficiency of dissipation in the gas. We also discuss the relevance of the Jeans instability and give an empirical criterion for the global bar instability in a two-component self-gravitating disk.

Shlosman, Isaac↗

A Particle-in-Cell Simulation for the Traveling Wave Direct Energy Converter (TWDEC) for Fusion Propulsion

A Particle-in-cell simulation model has been developed to study the physics of the Traveling Wave Direct Energy Converter (TWDEC) applied to the conversion of charged fusion products into electricity. In this model the availability of a beam of collimated fusion products is assumed; the simulation is focused on the conversion of the beam kinetic energy into alternating current (AC) electric power. The model is electrostatic, as the electro-dynamics of the relatively slow ions can be treated in the quasistatic approximation. A two-dimensional, axisymmetric (radial-axial coordinates) geometry is considered. Ion beam particles are injected on one end and travel along the axis through ring-shaped electrodes with externally applied time-varying voltages, thus modulating the beam by forming a sinusoidal pattern in the beam density. Further downstream, the modulated beam passes through another set of ring electrodes, now electrically oating. The modulated beam induces a time alternating potential di erence between adjacent electrodes. Power can be drawn from the electrodes by connecting a resistive load. As energy is dissipated in the load, a corresponding drop in beam energy is measured. The simulation encapsulates the TWDEC process by reproducing the time-dependent transfer of energy and the particle deceleration due to the electric eld phase time variations.

Chap, Andrew↗

Carbon Dioxide Clouds at High Altitude in the Tropics and in an Early Dense Martian Atmosphere

We use a time dependent, microphysical cloud model to study the formation of carbon dioxide clouds in the Martian atmosphere. Laboratory studies by Glandor et al. show that high critical supersaturations are required for cloud particle nucleation and that surface kinetic growth is not limited. These conditions, which are similar to those for cirrus clouds on Earth, lead to the formation of carbon dioxide ice particles with radii greater than 500 micrometers and concentrations of less than 0.1 cm(exp -3) for typical atmospheric conditions. Within the current Martian atmosphere, CO2 cloud formation is possible at the poles during winter and at high altitudes in the tropics during periods of increased atmospheric dust loading. In both cases, temperature perturbations of several degrees below the CO2 saturation temperature are required to nucleate new cloud particles suggesting that dynamical processes are the most common initiators of carbon dioxide clouds rather than diabatic cooling. The microphysical cloud model, coupled to a two-stream radiative transfer model, is used to reexamine the impact of CO2 clouds on the surface temperature within a dense CO2 atmosphere. The formation of carbon dioxide clouds leads to a warmer surface than what would be expected for clear sky conditions. The amount of warming is sensitive to the presence of dust and water vapor in the atmosphere, both of which act to dampen cloud effects. The radiative warming associated with cloud formation, as well as latent heating, work to dissipate the clouds when present. Thus, clouds never last for periods much longer than several days, limiting their overall effectiveness for warming the surface. The time average cloud optical depth is approximately unity leading to a 5-10 K warming, depending on the surface pressure. However, the surface temperature does not rise about the freezing point of liquid water even for pressures as high as 5 bars, at a solar luminosity of 75% the current value.

Colaprete, Anthony↗

Application of Three-Component PIV to a Hovering Rotor Wake

The key to accurate predictions of rotorcraft aerodynamics, acoustics, and dynamics lies in the accurate representation of the rotor wake. The vortical wake computed by rotorcraft CFD analyses typically suffer from numerical dissipation before the first blade passage. With some a priori knowledge of the wake trajectory, grid points can be concentrated along the trajectory to minimize the dissipation. Comprehensive rotorcraft analyses based on lifting-line theory rely on classical vortex models and/or semi-empirical information about the tip vortex structure. Until the location, size, and strength of the trailed tip vortex can be measured over a range of wake ages, the analyses will continue to be adjusted on a trial and error basis in order to correctly predict blade airloads, acoustics, dynamics, and performance. Using the laser light sheet technique, tip vortex location can be acquired in a straightforward manner. Measuring wake velocities and vortex core size, however, has been difficult and tedious using point-measurement techniques such as laser velocimetry. Recently, the Particle Image Velocimetry (PIV) technique has proven to be an efficient method for acquiring velocity measurements over relatively large areas and volumes of a rotor wake. The work reported to date, however, has been restricted to 2-component velocity measurements of the rotor wake. Three-component velocity measurements of a hovering rotor wake were acquired at NASA Ames Research Center in May 1999. This experiment represents a major step toward understanding the detailed structure of a rotor wake. This paper will focus primarily on the experimental technique used in acquiring this data. The accuracy and limitations of the current technique will also be discussed. Representative velocity field measurements will be included.

Yamauchi, Gloria K.↗

Diurnal Forcing of Planetary Atmospheres

This project investigated and developed models which will assist in developing an understanding of the diurnal forcings in the atmospheres of Mars, Venus, Earth and Jupiter. The differences in the temperatures during the day, and night, generate very deep diurnal boundary layers and strong slope-related winds. A general circulation model (GCM) was adapted for computation of Martian atmospheric dynamics, to simulate the global transport of these atmospheric tracers. The model was applied to the study of the Martian seasonal water cycle. The principles have been applied to venus. The region of study is the upper cloud level, and above rather than the area close to the surface. The model was also applied in efforts to study the evolution of volcanic clouds in the Earth's stratosphere. In particular this was used to model the dispersion of sulfuric acid particles in the atmosphere. A brief study of the application of the tidal theory to the planet Jupiter has been undertaken. The goal of the study is to deduce the source of the tidal dissipation which is thought to have led to the orbital resonances amongst the Galilean satellites and considerable heating of the interiors of Io, Europa, and Ganymede.

Houben, Howard C.↗

Overview of the relevant CFD work at Thiokol Corporation

An in-house developed proprietary advanced computational fluid dynamics code called SHARP (Trademark) is a primary tool for many flow simulations and design analyses. The SHARP code is a time dependent, two dimensional (2-D) axisymmetric numerical solution technique for the compressible Navier-Stokes equations. The solution technique in SHARP uses a vectorizable implicit, second order accurate in time and space, finite volume scheme based on an upwind flux-difference splitting of a Roe-type approximated Riemann solver, Van Leer's flux vector splitting, and a fourth order artificial dissipation scheme with a preconditioning to accelerate the flow solution. Turbulence is simulated by an algebraic model, and ultimately the kappa-epsilon model. Some other capabilities of the code are 2-D two-phase Lagrangian particle tracking and cell blockages. Extensive development and testing has been conducted on the 3-D version of the code with flow, combustion, and turbulence interactions. The emphasis here is on the specific applications of SHARP in Solid Rocket Motor design. Information is given in viewgraph form.

Chwalowski, Pawel↗

Dynamics of the earth's ring current - Theory and observation

The development of currents within an arbitrary distribution of particles trapped in the geomagnetic field is described. These currents combine to form the earth's ring current and thus are responsible for the worldwide depressions of surface magnetic field strength during periods of magnetic activity known as magnetic storms. Following a brief review of trapped particle motion in magnetic fields, ring current development is described and presented in terms of basic field and particle distribution parameters. Experimental observations then are presented and discussed within the theoretical framework developed earlier. New results are presented which, in the area of composition and charge state observations, hold high promise in solving many long standing ring current problems. Finally, available experimental results will be used to assess the present understanding as to ring current sources, generation, and dissipation.

Williams, D. J.↗

Effects of low frequency waves and spiky electric fields in the magnetotail

Intense low frequency waves and large amplitude spiky electric fields are commonly observed in the active magnetotail. These fields are of dynamic significance for substorms because the amplitude of the waves is large enough to provide the dissipation necessary for reconnection. The waves modify the trajectories of ions in the magnetotail, resulting in enhanced energization and pitch angle scattering in comparison with the trajectories obtained without waves. The spikes may represent the nonlinear evolution of the waves. Observations are presented of low frequency waves and spiky electric fields measured in the plasma sheet. The relationship between the waves and the spikes is compared to wave evolution theories. The Lundquist number is evaluated for the magnetotail based on the wave observations and previously published studies of the current sheet scale sizes. The effects of the waves on the ion trajectories are discussed using the results from a particle tracing code.

Cattell, C.↗

CDAW 8 observations of plasmoid signatures in the geomagnetic tail - An assessment

Magnetotail observations from the ISEE 3 distant (1983) tail mission taken during the Coordinated Data Analysis Workshop 8 (CDAW 8) A and G events are investigated. The ISEE 3 magnetic field, plasma, and energetic particle measurements taken in these two plasmoids have been analyzed and compared with various equilibrium structures and propagating waves/tail oscillation modes. Results indicate general agreement with either the closed-loop (Hones, 1977) or very small pitch angle flux rope (Hughes and Sibeck, 1987; Birn et al., 1989) models of plasmoid structure and poorer agreement with other hypotheses. Calculations based upon typical plasmoid and tail parameters are presented, indicating that the J and B force associated with the disconnected lobe field lines may be sufficient to accelerate plasmoids up to the speeds observed by ISEE 3. Overall, the energy expended in accelerating the plasmoids down the tail appears comparable to that dissipated in the inner magnetosphere and ionosphere. The study produces strong evidence in favor of the plasmoid model of substorm tail dynamics.

Slavin, J. A.↗

Radiative energy receiver for high performance energy conversion cycles

An analysis of gas dynamic processes pertinent to the functioning of earth-based and space-based solar electric power plants is presented, with attention given to potassium vapor as the working fluid. A device is described which features focused photon absorption by a nontransparent flowing gas. The feed flow is effected around the outside walls of a cavity receiver to raise efficiencies by trapping reemitted energy. A theoretical study of the interaction of a photon flux with a coaxial particle flux was performed, with the receiver flow treated as a Graetz flow. The critical parameters were defined, including a figure of merit as the gas enthalpy increase to absorbable energy ratio. A small-scale laboratory model was tested in comparison with the theoretically obtained values. Less than 15% of the absorbed energy was lost through dissipation while an 80% conversion efficiency was attained.

Rault, D.↗

What Might We Learn About Magnetospheric Substorms at the Earth from the MESSENGER Measurements at Mercury?

Satellite observations at the Earth, supported by theory and modeling, have established a close connection between the episodes of intense magnetospheric convection termed substorms and the occurrence of magnetic reconnection. Magnetic reconnection at the dayside magnetopause results in strong energy input to the magnetosphere. This energy can either be stored or used immediately to power the magnetospheric convection that produces the phenomena that collectively define the 'substorm.' However, many aspects of magnetic reconnection and the dynamic response of the coupled solar wind - magnetosphere - ionosphere system at the Earth during substorms remain poorly understood. For example, the rate of magnetic reconnection is thought to be proportional to the local Alfven speed, but the limited range of changes in this solar wind parameter at 1 AU have made it difficult to detect its influence over energy input to the Earth's magnetosphere. In addition, the electrical conductance of the ionosphere and how it changes in response to auroral charged particle precipitation are hypothesized to play a critical role in the development of substorms, but the nature of this electrodynamic interaction remain difficult to deduce from Earth observations alone. The amount of energy the terrestrial magnetosphere can store in its tail, the duration of the storage, and the trigger(s) for its dissipation are all thought to be determined by not only the microphysics of the cross-tail current layer, but also the properties of the coupled magnetosphere - ionosphere system. Again, the separation of microphysics effects from system response has proved very difficult using measurements taken only at the Earth. If MESSENGER'S charged particle and magnetic field measurements confirm the occurrence of terrestrial-style substorms in Mercury's miniature magnetosphere, then it may be possible to determine how magnetospheric convection, field-aligned currents, charged particle acceleration, reconnection, and tail energy storage are influenced by the intense magnetic reconnection expected to be associated with solar wind conditions at 0.3 - 0.5 AU from the Sun and the simplified electrodynamic feed-back anticipated for a planet lacking an ionosphere. MESSENGER observations from its 14 January 2008 Mercury flyby relevant to the occurrence of magnetic reconnection and substorms at Mercury will be discussed.

Slavin, James A.↗

Auroral ovals versus flare ribbons: Optical footprints of explosive magnetic field reconfiguration in Earth’s magnetosphere and the Sun’s corona

Despite drastically different plasma conditions in the two systems, eruptive energy release events in solar active regions and in Earth’s magnetosphere have a lot in common. Both solar active region eruptions and the substorm activity in the nightside magnetosphere are known to be associated with an abrupt large-scale reorganization of the supporting magnetic systems enabling the conversion and dissipation of the stored free magnetic energy. The rate of energy accumulation in either system tends to be by orders of magnitude slower than the energy release rate, leading to the formation of an impulsive loading-unloading cycle involving a relatively long pre-reconnection phase (substorm growth phase, pre-flare stressing of an active region magnetic field), a short and explosive reconnection phase (substorm breakup and expansion, eruptive flares and CMEs), and a prolonged post-reconnection phase (substorm recovery phase and the post-flare active region reconfiguration). In Earth’s magnetosphere, this loading-unloading cycle is diagnosed by monitoring characteristic changes in the particle precipitation pattern in the northern and southern auroral ovals marking the boundary between the open and close magnetic flux. The position, shape, structure, and dynamics of the auroral borealis were shown to communicate critically important information about an upcoming substorm, the location and timing of the initial reconnection onset, and the subsequent return of the magnetosphere to a lower-energy state. On the Sun, a similar role is played by the chromospheric flare ribbons which enable a detailed diagnostics of the overlying reconnecting magnetic field. In this talk, we’ll discuss similarities and differences in using auroral and flare ribbon observations for monitoring eruptive energy release in respectively Earth’s magnetosphere and solar corona. Special emphasis will be placed on the quantitative interpretation of the two types of optical signatures for empirical diagnostics of the underlying reconnection activity and for constraining its theoretical models. We’ll also compare forecasting capabilities of the two types of optical signatures in the context of space weather prediction objectives.

V. M. Uritskiy↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗

Overview of Results from the Cassini Magnetospheric Imaging Instrument (MIMI) During the First Year of Operations

The MIMI investigation comprises three sensors covering the indicated energy ranges: the Ion and Neutral Camera (INCA) -- 7 keV/nuc 3 Mev) and protons (1.6 < E < 160 Mev) from the back end of the dual field-of-view telescope. The Saturn observation sequences began in January, 2004 and culminated in Saturn Orbit Insertion on July 1, 2004. The MIMI sensors observed substantial activity in interplanetary space for several months prior to SOI, including several interplanetary shocks associated with corotating interaction regions, numerous increases most likely originating from particle streams in the vicinity of the Saturnian bow shock and energetic neutral atoms (ENA) emanating from Saturn s magnetosphere. Results following SOI revealed: a dynamical magnetosphere with a day-night asymmetry and an 11-hour periodicity; several water-product ions (O+, OH+, H2O+), but little N+; inferred quantities of neutral gas sufficient to cause major losses in the trapped ions and electrons in the middle and inner magnetosphere; a Titan exosphere that is a copious source of ENA; INCA imaging through ENA has also revealed a previously unknown radiation belt residing inward of the D-ring that is most likely the result of double charge-exchange between the main radiation belt and the upper layers of Saturn s exosphere. Finally, there is ample evidence for the presence of substorm-like injections of plasma that subsequently corotates for a number of days before dissipating on the night-side magnetotail. The observations will be presented and discussed in the context of current theoretical models.

Krimigis, S. M.↗