Search NASA⌕ Search

SEARCH · Search NASA

Results for “SPH”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

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.↗
NASA NTRS1994

Numerical simulations of catastrophic disruption: Recent results

Numerical simulations have been used to study high velocity two-body impacts. In this paper, a two-dimensional Largrangian finite difference hydro-code and a three-dimensional smooth particle hydro-code (SPH) are described and initial results reported. These codes can be, and have been, used to make specific predictions about particular objects in our solar system. But more significantly, they allow us to explore a broad range of collisional events. Certain parameters (size, time) can be studied only over a very restricted range within the laboratory; other parameters (initial spin, low gravity, exotic structure or composition) are difficult to study at all experimentally. The outcomes of numerical simulations lead to a more general and accurate understanding of impacts in their many forms.

Benz, W.↗
NASA NTRS1994

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↗
NASA NTRS2000

Modeling the Compression of Merged Compact Toroids by Multiple Plasma Jets

A fusion propulsion scheme has been proposed that makes use of the merging of a spherical distribution of plasma jets to dynamically form a gaseous liner. The gaseous liner is used to implode a magnetized target to produce the fusion reaction in a standoff manner. In this paper, the merging of the plasma jets to form the gaseous liner is investigated numerically. The Los Alamos SPHINX code, based on the smoothed particle hydrodynamics method is used to model the interaction of the jets. 2-D and 3-D simulations have been performed to study the characteristics of the resulting flow when these jets collide. The results show that the jets merge to form a plasma liner that converge radially which may be used to compress the central plasma to fusion conditions. Details of the computational model and the SPH numerical methods will be presented together with the numerical results.

Thio, Y. C. Francis↗
NASA NTRS2003

Giant Molecular Cloud Structure and Evolution

Bodenheimer and Burkert extended earlier calculations of cloud core models to study collapse and fragmentation. The initial condition for an SPH collapse calculation is the density distribution of a Bonnor-Ebert sphere, with near balance between turbulent plus thermal energy and gravitational energy. The main parameter is the turbulent Mach number. For each Mach number several runs are made, each with a different random realization of the initial turbulent velocity field. The turbulence decays on a dynamical time scale, leading the cloud into collapse. The collapse proceeds isothermally until the density has increased to about 10(exp 13) g cm(exp -3). Then heating is included in the dense regions. The nature of the fragmentation is investigated. About 15 different runs have been performed with Mach numbers ranging from 0.3 to 3.5 (the typical value observed in molecular cloud cores is 0.7). The results show a definite trend of increasing multiplicity with increasing Mach number (M), with the number of fragments approximately proportional to (1 + M). In general, this result agrees with that of Fisher, Klein, and McKee who published three cases with an AMR grid code. However our results show that there is a large spread about this curve. For example, for M=0.3 one case resulted in no fragmentation while a second produced three fragments. Thus it is not only the value of M but also the details of the superposition of the various velocity modes that play a critical role in the formation of binaries. Also, the simulations produce a wide range of separations (10-1000 AU) for the multiple systems, in rough agreement with observations. These results are discussed in two conference proceedings.

Hollenbach, David↗
NASA NTRS2003

Water Impact Test and Simulation of a Composite Energy Absorbing Fuselage Section

In March 2002, a 25-ft/s vertical drop test of a composite fuselage section was conducted onto water. The purpose of the test was to obtain experimental data characterizing the structural response of the fuselage section during water impact for comparison with two previous drop tests that were performed onto a rigid surface and soft soil. For the drop test, the fuselage section was configured with ten 100-lb. lead masses, five per side, that were attached to seat rails mounted to the floor. The fuselage section was raised to a height of 10-ft. and dropped vertically into a 15-ft. diameter pool filled to a depth of 3.5-ft. with water. Approximately 70 channels of data were collected during the drop test at a 10-kHz sampling rate. The test data were used to validate crash simulations of the water impact that were developed using the nonlinear, explicit transient dynamic codes, MSC.Dytran and LS-DYNA. The fuselage structure was modeled using shell and solid elements with a Lagrangian mesh, and the water was modeled with both Eulerian and Lagrangian techniques. The fluid-structure interactions were executed using the fast general coupling in MSC.Dytran and the Arbitrary Lagrange-Euler (ALE) coupling in LS-DYNA. Additionally, the smooth particle hydrodynamics (SPH) meshless Lagrangian technique was used in LS-DYNA to represent the fluid. The simulation results were correlated with the test data to validate the modeling approach. Additional simulation studies were performed to determine how changes in mesh density, mesh uniformity, fluid viscosity, and failure strain influence the test-analysis correlation.

Fasanella, Edwin L.↗
NASA NTRS2005

Gas Density Discontinuities in Merging Clusters

Chandra has discovered a new phenomenon in galaxy clusters, the sharp gas density edges. Depending on the sign of the temperature jump across the edge, these features may either be bow shocks or cold fronts. While bow shocks obviously are driven by merging sub-clusters, what causes cold fronts is not entirely clear, as they are observed both in mergers and in relaxed clusters. The purpose of the XMM study of A3376, an interesting cluster with density edges, is to understand the origin of cold fronts and to look for possible shocks. The XMM data for A3376 have been mostly analyzed (the X-ray edge turned out to be a cold front). Preliminary results have been shown at a conference and a paper is in preparation. We also have Chandra data for this cluster, and are comparing and combining the two datasets. In the course of analyzing the X-ray data for this cluster as well as several others, it has become apparent that we need the help of hydrodynamic simulations to study the precise mechanism by which cold fronts are formed, the main goal of the present project. A postdoc (Yago Ascasibar) is currently running SPH simulations of an idealized sub- cluster merger. These advanced simulations are nearing completion and two papers with their results are in preparation.

Mushotzky, Richard↗
NASA NTRS2010

Mapping and Modeling the Extended Winds of the Massive Interacting Binary, Eta Carinae

The combination HST/STIS high spatial and moderate spectral resolutions have revealed the massive interacting wind structure of Eta Carinae by forbidden lines of singly and doubly ionized elements. Throughout the 5.54-year period, lines of Fe++, Ne++, Ar++, S++ and N+ reveal the interacting wind structures, near critical electron densities of 10(exp 5) to 3 x 10(exp 7)cu cm, photoionized by the hot secondary, Eta Car B, Lines of Fe+ and Ni+ trace the denser (>10(exp 7)cu cm. less-ionized (< 8 eV) primary wind of Eta Car A as it wraps around the interacting binary stars. For 5 years of the 5.54 year period, the FUV radiation from Eta Car B escapes the orbital region, ionizing the boundaries of the expanding wind structures. But for three to six months, Eta Car B plunges into the primary wind approaching to within 1 to 2 AU, leading to cutoff of FUV and X-ray fluxes. The interacting wind structure, resolved out to 0.8", drops io ionization and then rebuilds as Eta Car B emerges from the primary wind envelope. Solid Particle Hydrodynamical(SPH) models have been developed extending out to 2000 AU and adapted to include FUV radiation effects of the winds. In turn, synthetic spectroimages of selected forbidden lines have been constructed and compared to the spectroimages recorded by the HST/STIS throughout 1998.0 to 2004.3, extending across the 1998 and 2003.5 minima. By this method, we show that the orbital axis of the binary system must bc within 15 degrees of the Homunculus axis of symmetry and that periastron occurs with Eta Car B passing on the far side of Eta Car B. This result ties the current binary orbit with the bipolar ejection with intervening skirt and leads to implications that the binary system influenced the mass ejection of the l840s and the lesser ejection of the 1890s.

Gull, Ted↗
NASA NTRS2011

Constraining the Absolute Orientation of eta Carinae's Binary Orbit: A 3-D Dynamical Model for the Broad [Fe III] Emission

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA(theta) that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38deg, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i approx. = 130deg to 145deg, Theta approx. = -15deg to +30deg, and an orbital axis projected on the sky at a P A approx. = 302deg to 327deg east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta(sub B), thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.

Madura, T. I.↗
NASA NTRS2012

Multi-Wavelength Implications of the Companion Star in eta Carinae

Eta-Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e approximately 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of eta-Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of eta-Car's observables at multiple wavelengths.

Madura, Thomas I.↗
NASA NTRS2011

Water Ice Clouds and Dust in the Martian Atmosphere Observed by Mars Climate Sounder

The water ice clouds are primarily controlled by the temperature structure and form at the water condensation level. Clouds in all regions presented show day/night differences. Cloud altitude varies between night and day in the SPH and tropics: (1) NPH water ice opacity is greater at night than day at some seasons (2) The diurnal thermal tide controls the daily variability. (3) Strong day/night changes indicate that the amount of gas in the atmosphere varies significantly. See significant mixtures of dust and ice at the same altitude planet-wide (1) Points to a complex radiative and thermal balance between dust heating (in the visible) and ice heating or cooling in the infrared. Aerosol layering: (1) Early seasons reveal a zonally banded spatial distribution (2) Some localized longitudinal structure of aerosol layers (3) Later seasons show no consistent large scale organization

Mars Climate Sounder (MCS) instrument↗
NASA NTRS2013

Effects of the LBV Primary's Mass-loss Rate on the 3D Hydrodynamics of eta Carinae's Colliding Winds

At the heart of eta Carinae's spectacular "Homunculus" nebula lies an extremely luminous (L(sub Total) greater than approximately 5 × 10(exp 6) solar luminosity) colliding wind binary with a highly eccentric (e approximately 0.9), 5.54-year orbit (Figure 1). The primary of the system, a Luminous Blue Variable (LBV), is our closest (D approximately 2.3 kpc) and best example of a pre-hypernova or pre-gamma ray burst environment. The remarkably consistent and periodic RXTE X-ray light curve surprisingly showed a major change during the system's last periastron in 2009, with the X-ray minimum being approximately 50% shorter than the minima of the previous two cycles1. Between 1998 and 2011, the strengths of various broad stellar wind emission lines (e.g. Halpha, Fe II) in line-of-sight (l.o.s.) also decreased by factors of 1.5 - 3 relative to the continuum2. The current interpretation for these changes is that they are due to a gradual factor of 2 - 4 drop in the primary's mass-loss rate over the last approximately 15 years1, 2. However, while a secular change is seen for a direct view of the central source, little to no change is seen in profiles at high stellar latitudes or reflected off of the dense, circumbinary material known as the "Weigelt blobs"2, 3. Moreover, model spectra generated with CMFGEN predict that a factor of 2 - 4 drop in the primary's mass-loss rate should lead to huge changes in the observed spectrum, which thus far have not been seen. Here we present results from large- (plus or minus 1620 AU) and small- (plus or minus 162 AU) domain, full 3D smoothed particle hydrodynamics (SPH) simulations of eta Car's massive binary colliding winds for three different primary-star mass-loss rates (2.4, 4.8, and 8.5 × 10(exp -4) solar mass/yr). The goal is to investigate how the mass-loss rate affects the 3D geometry and dynamics of eta Car's optically-thick wind and spatially-extended wind-wind collision (WWC) regions, both of which are known sources of observed X-ray, optical, UV, and near-IR emission and absorption. We use two domain sizes in order to better understand how the primary's mass-loss rate influences the various observables that form at different length scales. The 3D simulations provide information important for helping constrain η Car's recent mass-loss history and future state.

lavestigate↗
NASA NTRS2015

Launch Environment Water Flow Simulations Using Smoothed Particle Hydrodynamics

This paper describes the use of Smoothed Particle Hydrodynamics (SPH) to simulate the water flow from the rainbird nozzle system used in the sound suppression system during pad abort and nominal launch. The simulations help determine if water from rainbird nozzles will impinge on the rocket nozzles and other sensitive ground support elements.

Smoothed Particle Hyrodynamics↗
NASA NTRS2015

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.

Liou, Frank↗
NASA NTRS2015

A Coordinated X-Ray and Optical Campaign of the Nearby Massive Binary Sigma Orionis Aa: X-Ray Variability - II

We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution gratings spectral dataset of the Sigma Ori Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of approximately 479 kiloseconds and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range 5-25 angstroms is confirmed, with maximum amplitude of about plus or minus 15 percent within a single approximately 125 kiloseconds observation. Periods of 4.76 days and 2.04 days are found in the total X-ray flux, as well as an apparent overall increase in flux level throughout the 9-day observational campaign. Using 40 kiloseconds contiguous spectra derived from the original observations, we investigate variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phi equals 0.0 when the secondary Aa2 is at inferior conjunction. We use the results of an SPH radiative transfer code model, customized for this project, to relate the presence of a low density cavity in the primary stellar wind embedded shock that is associated with the secondary star to the emission line width variability.

binaries: closeaEuro"binaries: eclipsingaEuro"star↗
NASA NTRS2023

A Recent Impact Origin of Saturn’s Rings and Mid-Sized Moons

Saturn’s rings appear to be remarkably young, ∼100 Myr old, ruling out most previous hypotheses for their origin and incentivising new explanations. This could include a recent formation also for at least some of Saturn’s mid-sized icy moons. We present high-resolution SPH simulations that reveal how a range of collision scenarios driven by recent resonant instabilities in a precursor satellite system can scatter large masses of fragments and debris throughout the system that could evolve to both form the rings and reaccrete into present-day moons.

J Kegerreis↗
NASA NTRS2023

A Recent Impact Origin of Saturn’s Rings and Mid-Sized Moons

We simulate the collision of precursor icy moons analogous to Dione and Rhea as a possible origin forSaturn’s remarkably young rings. Such an event could have been triggered a few hundred million years ago by resonant instabilities in a previous satellite system. Using high-resolution SPH simulations, we find that this kind of impact can produce a wide distribution of massive objects and scatter material throughout the system. This includes the direct placement of pure-ice ejecta onto orbits that enterSaturn’s Roche limit, which could form or rejuvenate rings. In addition, fragments and debris of rock and ice totalling more than the mass of Enceladus can be placed onto highly eccentric orbits that would intersect with any precursor moons orbiting in the vicinity of Mimas, Enceladus, or Tethys. This could prompt further disruption and facilitate a collisional cascade to distribute more debris for potential ring formation, the re-formation of the present-day moons, and to evolve into an eventual cratering population of planeto-centric impactors.

Saturnian satellites↗