Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydrodynamical simulations”

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 433 records · Page 24

IUTAM Symposium on Hydrodynamic Diffusion of Suspended Particles

The focus of the symposium was on multiparticle hydrodynamic interactions which lead to fluctuating motion of the particles and resulting particle migration and dispersion or diffusion. Implications of these phenomena were described for sedimentation, fluidization, suspension flows, granular flows, and fiber suspensions. Computer simulation techniques as well as experimental techniques were described. Each session had an invited leadoff talk which overviewed the session topic as well as described the speaker's own related research. Ample time for discussion was included after each talk as well as at the end of each session. The symposium started with a keynote talk on the first evening on What is so puzzling about hydrodynamic diffusion?, which set the tone for the rest of the meeting by emphasizing both recent advances and unanswered issues.

Davis, R. H.↗

Improved Calibration of Modeled Discharge and Storage Change in the Atchafalaya Floodplain Using SAR Interferometry

This study focuses on the feasibility of using SAR interferometry to support 2D hydrodynamic model calibration and provide water storage change in the floodplain. Two-dimensional (2D) flood inundation modeling has been widely studied using storage cell approaches with the availability of high resolution, remotely sensed floodplain topography. The development of coupled 1D/2D flood modeling has shown improved calculation of 2D floodplain inundation as well as channel water elevation. Most floodplain model results have been validated using remote sensing methods for inundation extent. However, few studies show the quantitative validation of spatial variations in floodplain water elevations in the 2D modeling since most of the gauges are located along main river channels and traditional single track satellite altimetry over the floodplain are limited. Synthetic Aperture Radar (SAR) interferometry recently has been proven to be useful for measuring centimeter-scale water elevation changes over the floodplain. In the current study, we apply the LISFLOOD hydrodynamic model to the central Atchafalaya River Basin, Louisiana, during a 62 day period from 1 April to 1 June 2008 using two different calibration schemes for Manning's n. First, the model is calibrated in terms of water elevations from a single in situ gauge that represents a more traditional approach. Due to the gauge location in the channel, the calibration shows more sensitivity to channel roughness relative to floodplain roughness. Second, the model is calibrated in terms of water elevation changes calculated from ALOS PALSAR interferometry during 46 days of the image acquisition interval from 16 April 2008 to 1 June 2009. Since SAR interferometry receives strongly scatters in floodplain due to double bounce effect as compared to specular scattering of open water, the calibration shows more dependency to floodplain roughness. An iterative approach is used to determine the best-fit Manning's n for the two different calibration approaches. Results suggest similar floodplain roughness but slightly different channel roughness. However, application of SAR interferometry provides a unique view of the floodplain flow gradients, not possible with a single gauge calibration. These gradients, allow improved computation of water storage change over the 46-day simulation period. Overall, the results suggest that the use of 2D SAR water elevation changes in the Atchafalaya basin offers improved understanding and modeling of floodplain hydrodynamics.

Jung, Hahn Chul↗

Impact on comet Shoemaker-Levy 9 on Jupiter

Three-dimensional numerical simulations of the impact of Comet Shoemaker - Levy 9 on Jupiter and the resulting vapor plume expansion were conducted using the Smoothed Particle Hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the atmosphere between -100 km to -250 km. This energy is converted to potential energy of the resulting gas plume. The unconfined plume expands vertically and has a peak radiative power approximately equal to the total radiation from Jupiter's disk. The plume rises a few tens of atmospheric scale heights in approximately 10(exp 2) seconds. The rising plume reaches the altitude of approximately 3000 km, but no atmospheric gas is accelerated to the escape velocity (approximately 60 km/s).

Ahrens, Thomas J.↗

Impact of comet Shoemaker-Levy 9 on Jupiter

We have employed three-dimensional numerical simulations of the impact of Comet Shoemaker-Levy 9 (SL9) on Jupiter and the resulting vapor plume expansion using the smoothed particle hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the atmosphere between -100 to -250 km. This energy is converted to potential energy of the resulting gas plume. The unconfined plume expands vertically and has a peak radiative power approximately equal to the total radiation from Jupiter's disc. The plume rises a few tens of atmospheric scale heights in approximately 10(exp 2) seconds. The rising plume reaches the altitude of approximately 3000 km; however, no atmospheric gas is accelerated to the escape velocity (approximately 60 km/s).

Takata, Toshiko↗

Thermonuclear processes and accretion onto neutron star envelopes - X-ray burst and transient sources

A Lagrangian, fully implicit, one-dimensional hydrodynamic computer code is used to investigate the evolution of thermonuclear runaways in the thick accreted hydrogen-rich envelopes of 1.0-solar-mass neutron stars having radii of 10 km and 20 km. The simulations produce outbursts that range in time scale from about 2000 seconds to more than a day. For the 10-km study, the peak effective temperature is 3.3 x 10 to the 7th K, and the peak luminosity is 2 x 10 to the 5th solar luminosities. The 20-km neutron star produces a peak effective temperature and luminosity of 5.3 x 10 to the 6th K and 5.9 x 10 to the 2nd solar luminosities. Also investigated are the effects of changes in the rates of the O-14(alpha, proton) and O-15(alpha, gamma) reactions on the evolution. Hydrodynamic expansion on the 10-km neutron star produces a precursor lasting about 10 to the -6th sec. The evolution of a gas cloud impacting the surface of a 20-km, 1-solar-mass neutron star is studied in an attempt to simulate the magnetospheric gate model of the X-ray burst sources.

Starrfield, S.↗

Simulation of Relativistic Shocks and Associated Radiation from Turbulent Magnetic Fields

Using our new 3-D relativistic particle-in-cell (PIC) code, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value as predicted by hydrodynamic compression. Behind the bow shock, in the jet shock, strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. In order to go beyond the standard synchrotron model used in astrophysical objects we have used PIC simulations and calculated radiation based on first principles. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the technique to calculate emission from electrons based on simulations with a small system. We obtain spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system may generate a jitter/synchrotron spectrum.

Nishikawa, K.-I.↗

Numerical simulations of the convective flame in white dwarfs

A first step toward better understanding of the mechanism driving convective flames in exploding white dwarfs is presented. The propagation of the convective flame is examined using a two-dimensional implicit hydrodynamical code. The large scales of the instability are captured by the grid while the scales that are smaller than the grid resolution are approximated by a mixing-length approximation. It is found that largescale perturbations (of order of the pressure scale height) do grow significantly during the expansion, leading to a very nonspherical burning front. The combustion rate is strongly enhanced (compared to the unperturbed case) during the first second, but later the expansion of the star suppresses the flame speed, leading to only partial incineration of the nuclear fuel. Our results imply that large-scale perturbations by themselves are not enough to explain the mechanism by which convective flames are driven, and a study of the whole spectrum of relevant perturbations is needed. The implications of these preliminary results on future simulations, in the context of current models for Type Ia supernovae, are discussed.

Livne, Eli↗

Dynamic simulation of coronal mass ejections

A model is developed for the formation and propagation through the lower corona of the loop-like coronal transients in which mass is ejected from near the solar surface to the outer corona. It is assumed that the initial state for the transient is a coronal streamer. The initial state for the streamer is a polytropic, hydrodynamic solution to the steady-state radial equation of motion coupled with a force-free dipole magnetic field. The numerical solution of the complete time-dependent equations then gradually approaches a stationary coronal streamer configuration. The streamer configuration becomes the initial state for the coronal transient. The streamer and transient simulations are performed completely independent of each other. The transient is created by a sudden increase in the pressure at the base of the closed-field region in the streamer configuration. Both coronal streamers and coronal transients are calculated for values of the plasma beta (the ratio of thermal to magnetic pressure) varying from 0.1 to 100.

Steinolfson, R. S.↗

Hydrodynamic water impact

The hydrodynamic impact of a falling body upon a viscous incompressible fluid was investigated by numerically solving the equations of motion. Initially the mathematical model simulated the axisymmetric impact of a rigid right circular cylinder upon the initially quiescent free surface of a fluid. A compressible air layer exists between the falling cylinder and the liquid free surface. The mathematical model was developed by applying the Navier-Stokes equations to the incompressible air layer and the incompressible fluid. Assuming the flow to be one dimensional within the air layer, the average velocity, pressure and density distributions were calculated. The liquid free surface was allowed to deform as the air pressure acting on it increases. For the liquid the normalized equations were expressed in two-dimensional cylindrical coordinates. The governing equations for the air layer and the liquid were expressed in finite difference form and solved numerically. For the liquid a modified version of the Marker-and-Cell method was used. The mathematical model has been reexamined and a new approach has recently been initiated. Essentially this consists of examining the impact of an inclined plate onto a quiesent water surface with the equations now formulated in cartesian coordinates.

Kettleborough, C. F.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative, manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. These results are highlighted. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable. Another one of our goals is to examine the stability of our predicted flows. Because time-dependent three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, Robert W.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program was to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, with the funding from this NASA program, we concentrated on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we developed unique global models that allowed us to study the coupling between the different regions. These results are highlighted in the next section. Another important aspect of our NASA Theory Program concerned the effect that localized 'structure' had on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkland current patterns) or time variations in these input due to storms and substorms. Also, some of the plasma flows that we predicted with our macroscopic models could be unstable, and another one of our goals was to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another goal of our NASA Theory Program was to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and hydrodynamic predictions for a given polar wind scenario or it could involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations provides insight into when the various models can be used with confidence.

Schunk, Robert W.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗

Energy transport and dynamics

We report findings concerning energy transport and dynamics in flares during the impulsive and gradual phases based on new ground-based and space observations (notably from Yohkoh). A preheating sometimes occurs during the impulsive phase. Ca XIX line shifts are confirmed to be good tracers of bulk plasma motions, although strong blue shifts are not as frequent as previously claimed. They often appear correlated with hard X-rays but, for some events, the concept that electron beams provide the whole energy input to the thermal component seems not to apply. Theory now yields: new diagnostics of low-energy proton and electric beams; accurate hydrodynamical modeling of pulse beam heating of the atmosphere; possible diagnostics of microflares (based on X-ray line ratio or on loop variability); and simulated images of chromospheric evaporation fronts. For the gradual phase, the continual reorganization of magnetic field lines over active regions determines where and when magnetic reconnection, the mechanism favored for energy release, will occur. Spatial and temporal fragmentation of the energy release, observed at different wavelengths, is considered to be a factor as well in energy transport and plasma dynamics.

Schmieder, Brigitte↗

Hematopoietic Stem Cell Therapy to Countermeasure Cancer in Astronauts during Exploration of Deep Space

Exposure to cosmic radiation can cause chromosomal mutations, which may lead to cancer in astronauts engaged in space exploration. Therefore, our goals are to develop countermeasures to prevent space-induced cancer using hematopoietic stem cell therapy (HSCT) and gene therapy. This presentation focuses on HSCT for cancer. Our previous experiments on a simulated, space-induced immuno-deficiency model (mouse hind limb unloading ) indicated that transplanted hematopoietic stem cells (HSCs) could enhance the host's immunity by effectively eliminating bacterial infection (Ohi S, et. al. J Grav Physiol 10, P63-64, 2003; Ohi S, et. al. Proceedings of the Space Technology and Applications International Forum (STAIF) . American Institute of Physics, New York, pp. 938-950, 2004). Hence, we hypothesized that the HSCs might be effective in combating cancer as well. Studies of cocultured mouse HSCs with beta-galactosidase marked rat gliosarcoma spheroids (9L/lacZ), a cancer model, indicated antagonistic interactions , resulting in destruction of the spheroids by HSCs. Trypan Blue dye-exclusion assays were consistent with the conclusion. These results show potential usehlness of HSCT for cancer. Currently, the NASA Hydrodynamic Focusing Bioreactor (HFB), a space analog tissue/cell culture system, is being used to study invasion of the gliosarcoma (GS) spheroids into mouse brain with or without co-cultured HSCs. This may simulate the metastasis of gliosarcoma to brain. There is a tendency for the HSCs to inhibit invasion of GS spheroids into brain, as evidenced by the X-gal staining.

Ohi, S.↗

Collapse of primordial gas clouds and the formation of quasar black holes

The formation of quasar black holes during the hydrodynamic collapse of protogalactic gas clouds is discussed. The dissipational collapse and long-term dynamical evolution of these systems is analyzed using three-dimensional numerical simulations. The calculations focus on the final collapse stages of the inner baryonic component and therefore ignore the presence of dark matter. Two types of initial conditions are considered: uniformly rotating spherical clouds, and iirotational ellipsoidal clouds. In both cases the clouds are initially cold, homogeneous, and not far from rotational support (T/(absolute value of W) approximately equals 0.1). Although the details of the dynamical evolution depend sensitively on the initial conditions, the qualitative features of the final configurations do not. Most of the gas is found to fragment into small dense clumps, that eventually make up a spheroidal component resembling a galactic bulge. About 5% of the initial mass remains in the form of a smooth disk of gas supported by rotation in the gravitational potential potential well of the outer spheroid. If a central seed black hole of mass approximately greater than 10(exp 6) solar mass forms, it can grow by steady accretion from the disk and reach a typical quasar black hole mass approximately 10(exp 8) solar mass in less than 5 x 10(exp 8) yr. In the absence of a sufficiently massive seed, dynamical instabilities in a strongly self-gravitating inner region of the disk will inhibit steady accretion of gas and may prevent the immediate formation of quasar.

Loeb, Abraham↗

Orbital Debris Shape and Orientation Effects on Ballistic Limits

The SPHC hydrodynamic code was used to evaluate the effects of orbital debris particle shape and orientation on penetration of a typical spacecraft dual-wall shield. Impacts were simulated at near-normal obliquity at 12 km/sec. Debris cloud characteristics and damage potential are compared with those from impacts by spherical projectiles. Results of these simulations indicate the uncertainties in the predicted ballistic limits due to modeling uncertainty and to uncertainty in the impactor orientation.

Evans, Steven W.↗

Convective effects in float-zone and Czochralski melts

The hydrodynamics of crystal-growth melts is a relatively new research area. Numerical modeling of these processes is necessary. The work discussed herein is in two parts: numerical simulations of the flow in a Czochralski melt, and also of that in a float zone. In addition, for the float-zone case, energy stability theory will be used to determine stability bounds for the onset of oscillatory thermo-capillary flow. Convective effects in crystal-growth melts arise from a variety of mechanisms. Temperature gradients both in the direction of gravity and normal to it give rise to convection due to buoyancy effects. Rotation of the crucible and/or crystal causes a forced convection which may augment or oppose the buoyancy-driven flow. Finally, thermo-capillary forces (due to the variation of surface tension with temperature) drive surface motions which in turn generate convection in the bulk fluid. All of these mechanisms are present in either Czochralski or float-zone growth. The objective of the Czochralski modeling is to develop an accurate numerical simulation of the flow in a Czochralski silicon melt and to investigate the effects of various parameters on the flow properties. Like some earlier investigations, the intent is to simulate the effects of buoyancy, forced and thermo-capillary convection, including unsteady effects. Unlike earlier work, the aim is to include the effects of a variable free surface and freezing interface and, possibly incorporate nonaxisymmetric effects.

Neitzel, G. P.↗

Three-dimensional compressible hydrodynamic convection in the sun and stars

A pseudospectral code is used to obtain numerical solutions of 3D compressible hydrodynamic convection in a stratified medium where the boundaries are open and where radiative losses are strongest at the top. Results are presented for a series of three simulations in progressively deeper boxes. The code does not make the anelastic approximation, but follows sound waves explicitly. Spatial organization is found in the form of vertical 'stacking' of smaller structures on top of larger ones, with each structure being roughly Hp to a vertical extent. Compressibility effects are apparent in the density: the snapshots indicate that, on a horizontal plane, the density at certain points may be of an order of 10 percent larger or smaller than elsewhere on that plane.

Hossain, Murshed↗