Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smoothed Particle Hydrodynamics”

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.

88 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Constraining the Properties of the Eta Carinae System via 3-D SPH Models of Space-Based Observations: The Absolute Orientation of the Binary Orbit

The extremely massive (> 90 Stellar Mass) and luminous (= 5 x 10(exp 6) Stellar Luminosity) star Eta Carinae, with its spectacular bipolar "Homunculus" nebula, comprises one of the most remarkable and intensely observed stellar systems in the Galaxy. However, many of its underlying physical parameters remain unknown. Multiwavelength variations observed to occur every 5.54 years are interpreted as being due to the collision of a massive wind from the primary star with the fast, less dense wind of a hot companion star in a highly elliptical (e approx. 0.9) orbit. Using three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the binary wind-wind collision, together with radiative transfer codes, we compute synthetic spectral images of [Fe III] emission line structures and compare them to existing Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS) observations. We are thus able, for the first time, to tightly constrain the absolute orientation of the binary orbit on the sky. An orbit with an inclination of approx. 40deg, an argument of periapsis omega approx. 255deg, and a projected orbital axis with a position angle of approx. 312deg east of north provides the best fit to the observations, implying that the orbital axis is closely aligned in 3-D space with the Homunculus symmetry axis, and that the companion star orbits clockwise on the sky relative to the primary.

Madura, Thomas I.↗
NASA NTRS2011

Capabilities of the Impact Testing Facility at Marshall Space Flight Center

The test and analysis capabilities of the Impact Testing Facility at NASA's Marshall Space Flight Center are described. Nine different gun systems accommodate a wide range of projectile and target sizes and shapes at velocities from subsonic through hypersonic, to accomplish a broad range of ballistic and hypervelocity impact tests. These gun systems include ballistic and microballistic gas and powder guns, a two-stage light gas gun, and specialty guns for weather encounter studies. The ITF "rain gun" is the only hydrometeor impact gun known to be in existence in the United States that can provide single impact performance data with known raindrop sizes. Simulation of high velocity impact is available using the Smooth Particle Hydrodynamic Code. The Impact Testing Facility provides testing, custom test configuration design and fabrication, and analytical services for NASA, the Department of Defense, academic institutions, international space agencies, and private industry in a secure facility located at Marshall Space Flight Center, on the US Army's Redstone Arsenal in Huntsville, Alabama. This facility performs tests that are subject to International Traffic in Arms Regulations (ITAR) and DoD secret classified restrictions as well as proprietary and unrestricted tests for civil space agencies, academic institutions, and commercial aerospace and defense companies and their suppliers.

Finchum, Andy↗
NASA NTRS2013

PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils

PowderSim is a calculation tool that combines a discrete-element method (DEM) module, including calibrated interparticle-interaction relationships, with a mesh-free, continuum, SPH (smoothed-particle hydrodynamics) based module that utilizes enhanced, calibrated, constitutive models capable of mimicking both large deformations and the flow behavior of regolith simulants and lunar regolith under conditions anticipated during in situ resource utilization (ISRU) operations. The major innovation introduced in PowderSim is to use a mesh-free method (SPH-based) with a calibrated and slightly modified critical-state soil mechanics constitutive model to extend the ability of the simulation tool to also address full-scale engineering systems in the continuum sense. The PowderSim software maintains the ability to address particle-scale problems, like size segregation, in selected regions with a traditional DEM module, which has improved contact physics and electrostatic interaction models.

Johnson, Scott↗
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 NTRSDate not supplied

Single-Scattering Properties of Melting Precipitation for GPM Passive Microwave and Radar Remote Sensing Applications

Over the past two decades, detailed computational simulations of the intricate three dimensional structures of ice-phase crystals and aggregates of those crystals have been developed. The microwave single-scattering properties of these simulated ice particles have been computed and used to improve quantitative estimates of snow rates and to better deRne the vertical structure of snow water contents in deep convective systems, as derived from satellite-borne passive microwave and/or radar remote sensing measurements. The same icephase particles have more recently been used as the starting point for simulations of melting precipitation using computational melting methods. In the current study, a heuristic melting method, as well as a physically-based melting procedure based on smoothed-particle hydrodynamics, are applied to ice particle models to describe the full evolution of the particles from dry snow to liquid drops. The discrete dipole approximation is utilized to calculate the single-scattering properties of the mixed-phase particles throughout the melting process. Then, the properties of the particles are “mapped” into simpliRed microphysical simulations of particle size spectra in the melting layers of stratiform, precipitating clouds. The bulk single-scattering properties of the melting layers and the sensitivity of their properties to modeling assumptions are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗
NASA NTRSDate not supplied

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗
NASA NTRSDate not supplied

Case Study of Meteoroid Fluence on a Mars Sample Return Mission

This paper presents results from an analysis to assess the fluence of meteoroids on the Mars Sample Return Campaign Micrometeoroid Protection System (MMPS) and probability of damage to the Earth Entry System Thermal Protection System (TPS). Results are presented for the mission Outbound trajectory from Earth to Mars. The MEM3 meteoroid environment model was used for the assessment. Results are presented showing the reduction in fluence on cube faces and on the MMPS that can be achieved by selection of an anti-velocity orientation of the MMPS. A methodology is presented that determines the cumulative fluence on the MMPS using the MEM3 meteoroid fluxes and the spacecraft attitude profile over the trajectory. Resulting distributions of fluence over impact velocity, angle, and mass are presented. The fluence results are used to establish velocity-angle scenarios for hydrocode simulation and to assess net damage fluence to the TPS. Example tables containing Smooth Particle Hydrodynamics Code (SPHC) simulation scenarios as well as hypothetical hydrocode simulation scenarios for illustration purposes are presented. The tables present the resulting net damaging fluence to the TPS for each scenario and the total across the scenarios. Additionally, examples of CTH hydrocode simulation are included, demonstrating their use for determining critical impactor across the impactor parameter space, and damage of impactors that penetrate the MMPS.

Micrometeoroid↗
NASA NTRSDate not supplied

Extending the Applicability of Thermal Protection System Ballistic Limit Equations Beyond the Testable Regime

All spacecraft use some sort of thermal insulation in their design, which is usually placed on the spacecraft’s exterior. As such, it is highly susceptible to high-speed impacts by meteoroids and orbital debris. These high-speed impacts can damage the insulation to a point where the protection it offers is below acceptable limits. It is, therefore, important to be able to characterize expected TPS damage levels stemming from such high-speed impacts. In a recent paper, we presented the results of a study that sought to characterize the high-speed impact damage that would be sustained by two thermal protection materials that have recently gained attention for possible use in interplanetary missions. As part of the analyses performed, ballistic limit equations for these materials were developed using recent high-speed impact test data, and then assessed by comparing their predictions against test data. In this paper, we extend the applicability of the ballistic limit equation for one of those thermal materials beyond the testable regime using the results of series of Smooth Particle Hydrodynamics Code (SPHC) impact simulations at velocities in excess of 10 km/s. When the predictions of the updated and expanded BLE were compared against the results of the SPHC runs, the predictions were seen to be consistent with the data: non-perforation data points were typically below the BLE curves, while perforation data points were typically above them.

Thermal Protection System↗
NASA NTRS2026

The Role of Disruptive Impacts on Ocean Generation and Longevity in Icy Moons

Among icy moons of the outer Solar System, subsurface water oceans and large collisions both seem to be common, but the impact of the latter on the presence and persistence of the former is unclear and has rarely been investigated. Here we interface a smoothed-particle hydrodynamics model to simulate collisions with a thermal-structural evolution model to simulate the evolution of moons pre-collision, post-collision and without a collision. Overall, even such large-scale collisions affect only the ocean thickness or longevity, and the presence or absence of an ocean is affected for only a part of a moon’s history. In reaccreted moons, the ocean survives the impact and becomes much thicker inside larger moons with radius near 1,000 km, whereas an ocean that would otherwise arise inside moons with radius near 500 km is absent because the collision promotes ice–rock differentiation. Our simulations have not yielded an ocean developed post-impact—whether directly via collisional or reaccretional heating or indirectly through tidal heating due to collision-induced orbital changes—in a moon that would otherwise have remained frozen. The ocean-enhancing effect is pronounced only for late disruptive impacts onto large, 1,000-km-class targets, which are unlikely in recent solar system history.

Marc Neveu↗
NASA NTRSDate not supplied

Low-G Slosh Workshop Results From 2023: State of the Art, Gaps and Forward Work

The NASA Engineering and Safety Center (NESC) sponsored a workshop to dis-cuss low-g slosh. Over a dozen different NASA centers and industry partners participated in this workshop. The main topics were to identify the state of the art regarding test data, Computational Fluid Dynamics (CFD) analysis resources, and GN&C models that are available for this problem. The primary intent was to determine deficiencies in the collective knowledge base and to identify priorities that should be established for proactive work in the areas of model development, validation, and model integration. Recommendations from workshop participants include establishing an online slosh modeling resource, conducting ground/flight testing, developing surrogate models, and establishing co-simulation techniques. This paper summarizes the workshop findings and recommendations for the larger community.

Lattice-Boltzmann↗
NASA NTRS2009

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗
NASA NTRS1996

[Two-phase Interfaces in Weak External Fields]

Our aim has been that of understanding from first principles the behavior of two-phase interfaces in the absence of gravitational constraints. This is fundamental to our ability to deal with the fluid structures that abound in the real biological, chemical, and physical world. A substantial effort was mounted to determine how familiar hydrodynamic concepts have to be modified and interpreted to make them appropriate to the multi-level structure alluded to above. This was primarily in the context of the microscopic symmetric pressure tensor, which was, for the first time, expressed in the invaluable density functional format, and the used to follow the predictions of popular microscopic models of the energetics of interfacial systems. In the course of these investigations, the previous murky relation between pressure tensor and thermodynamics was completely clarified. The process of extending thermodynamic information to interfacial dynamics was initiated along two paths. One was from the viewpoint of an inertialess lattice gas, resulting in the surprising conclusion that at this level, all transport is governed by precisely the thermodynamic free energy, albeit with a non-trivial effective particle mobility. The other aimed at understanding the fashion in which slow macroscopic motions, accounted for by a time-varying microscopic energy, generate effective hydrodynamic parameters. By examining a solvable model system, it was found that all current procedures for doing so are deficient, and suitable alleviation suggested. The major effect of this project was to set the stage for the analysis of the substantial dynamical regimes in which extensive equilibrium information provides the dominant background. This produces a smooth junction to the models of Araki and Munakata, Giacomin and Lebowitz, and Oxtoby. It is also crucial to our understanding of the complex interfacial equilibrium configurations required for intermediate stages of two-phase separation, for which the "phase-field" techniques we have been developing are uniquely effective. And it puts us within striking range of the physical experiments that can provide incisive tests of the theoretical framework.

Percus, J. K.↗
NASA NTRS1983

Monte Carlo simulation of steady state shock structure including cosmic ray mediation and particle escape

Both hydrodynamic calculations (Drury and Volk, 1981, and Axford et al., 1982) and kinetic simulations imply the existence of thermal subshocks in high-Mach-number cosmic-ray-mediated shocks. The injection efficiency of particles from the thermal background into the diffusive shock-acceleration process is determined in part by the sharpness and compression ratio of these subshocks. Results are reported for a Monte Carlo simulation that includes both the back reaction of accelerated particles on the inflowing plasma, producing a smoothing of the shock transition, and the free escape of particles allowing arbitrarily large overall compression ratios in high-Mach-number steady-state shocks. Energy spectra and estimates of the proportion of thermal ions accelerated to high energy are obtained.

Ellison, D. C.↗
NASA NTRS1974

Quantum mechanical streamlines. I - Square potential barrier

Exact numerical calculations are made for scattering of quantum mechanical particles hitting a square two-dimensional potential barrier (an exact analog of the Goos-Haenchen optical experiments). Quantum mechanical streamlines are plotted and found to be smooth and continuous, to have continuous first derivatives even through the classical forbidden region, and to form quantized vortices around each of the nodal points. A comparison is made between the present numerical calculations and the stationary wave approximation, and good agreement is found between both the Goos-Haenchen shifts and the reflection coefficients. The time-independent Schroedinger equation for real wavefunctions is reduced to solving a nonlinear first-order partial differential equation, leading to a generalization of the Prager-Hirschfelder perturbation scheme. Implications of the hydrodynamical formulation of quantum mechanics are discussed, and cases are cited where quantum and classical mechanical motions are identical.

Hirschfelder, J. O.↗
NASA NTRS1994

The structure of mass-loading shocks. 1: Comets

A new multifluid model to describe the solar wind-cometary ion plasma in the outer cometary coma is derived. This model is distinguished from previous multifluid models in that wave-particle effects are included explicitly. By considering hydrodynamic timescales, self-consistent dissipative terms (which correspond to the spatial diffusion of cometary ions) are derived, thereby allowing the structure of a cometary bow shock to be resolved. In the hypersonic limit, our diffusive multifluid model reveals an attractive connection to the one-fluid model formulated by Biermann et al. (1967), a connection which serves as an important guide when investigating our more complicated nonhypersonic diffusive model. The detailed model consists of cometary ions diffusing in a warm solar wind which is heated by the dissipation of Alfvenic turbulence (in situ and generated by the initial ring beam distribution). The quasi-parallel cometary shock is found to be at least an order of magnitude thicker than a quasi-perpendicular cometary shock. It is furthermore found that for undisturbed solar wind (sonic) Mach numbers in excess of approximately 6, the cometary bow shock (for both perpendicular and parallel cases) is smoothed completely by the cometary ions and no 'proton' or thermal subshock is necessary. Finally, a detailed comparison of the model with the plasma observations made at Halley by Giotto is undertaken. It is found that the observed shock structure, thickness, location, and plasma parameters compare well with the predictions of the theoretical model, particularly in the case of the quasi-parallel shock.

Zank, G. P.↗
NASA NTRS2014

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗