A coupled smoothed particle hydrodynamics-finite volume approach for shock capturing in one-dimension
Not Available
42 ENGINEERING↗
SEARCH · Engineering Papers
Results for “smoothed particles 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.
Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design
Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.
Fahrenthold, Eric P.↗
Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science
We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Efficient smoothed particle radiation hydrodynamics II: Radiation hydrodynamics
Not Available
Bassett, Brody R. (ORCID:0000000333109451)↗
A Coupled Smoothed Particle Hydrodynamics-Finite Volume Approach for Shock Capturing in One-Dimension.
Abstract not provided.
Myers, Connor↗
qSmoothed particle interpolation
Smoothed particle hydrodynamics (SPH) discretization techniques are generalized to develop a method, smoothed particle interpolation (SPI), for solving initial value problems of systems of a nonhydrodynamical nature. Under this approach, SPH is viewed as strickly as interpolation scheme and, as such, suitable for solving general hyperbolic and parabolic equations. The SPI method is tested on (1) the wave equation with inhomogenous sound speed and (2) Burger's equation. The efficiency of SPI is studied by comparing SPI solutions to those obtained with standard finite-difference methods. It is shown that the power of SPI arises when the smoothing particles are free to move.
Laguna, Pablo↗
Co-Extrusion of Dissimilar Aluminum Alloys via Shear-Assisted Processing and Extrusion
Bimetallic tubes are used when a component requires more than one performance requirement, for instance, strength/creep resistance and oxidation/corrosion resistance. Shear-assisted processing and extrusion (ShAPE) can be used to fabricate extrudates that are comparable or superior in performance relative to conventionally manufactured extrudates. For the first time, ShAPE has been successfully employed in producing bimetallic Al tubing consisting of 6061 and 7075 Al alloys. Both light and electron microscopy techniques were used to investigate the integrity of the tubes, especially the interface between the core and cladding of the bimetallic tubes. Void-free bimetallic tubes were produced using ShAPE. Quantification of tube integrity was carried out with tensile testing under as-extruded and post-aged conditions. All the bimetallic tubes in the as-extruded samples exhibited uniform elongation above 5% with good tensile strength. Key insights such as material flow during bimetallic tube extrusion were obtained from the characterization of remnant billet and simulation results from a smoothed particle hydrodynamics model.
36 MATERIALS SCIENCE↗
Hydrodynamics of binary coalescence. 1: Polytropes with stiff equations of state
We performed a series of three-dimensional hydrodynamic calculations of binary coalescence using the smoothed particle hydrodynamics (SPH) method. The initial conditions are exact polytropic equilibrium configurations on the verge of dynamical instability. We consider synchronized equilbria only and concentrate on stiff equations of state, with adiabatic Gamma greater than 5/3. We assume that the polytropic constants (K identically equal to P/(rho(exp Gamma)) are the same for both components. These conditions apply well to models of neutron star binaries. Accordingly, we discuss our results in the context of the Laser Interferometer Gravitational-Wave Observatory (LIGO) project, and we calculate the emission of gravitational radiation in the quadruple approximation. The fully nonlinear development of the instability is followed using SPH until a new equilibrium configuration is reached by the system. We find that the properties of this final configuration depend sensitively on both the compressibility and mass ratio. An axisymmetric merged configuration is always produced when the adiabatic exponent Gamma approximately less than 2.3. As a consequence, the emission of gravitational radiation shuts off abruptly right after the onset of dynamical instability. In contrast, triaxial merged configurations are obtained when Gamma approximately greater than 2.3, and the system continues to emit gravitational waves after the final coalescence. Systems with mass ratios q not equal to 1 typically become dynamically unstable before the onset of mass transfer. Stable mass transfer from one neutron star to another in a close binary is therefore probably ruled out. For a mass ratio q approximately less than 0.5, however, dynamical mass transfer can temporarily retard the coalescence by causing a rapid reexpansion of the binary into a new, slightly eccentric but dynamically stable orbit. The maximum amplitude h(sub max) and peak luminosity L(sub max) of the gravitational waves emitted during the final coalescence are nearly independent of Gamma, but depend sensitively on the mass ratio q. The approximate scalings we find are h(sub max) varies as q(exp 2) and L(sub max) varies as q(exp 6) for q close to unity.
Rasio, Frederic A.↗
Characterizing Wheel-Soil Interaction Loads Using Meshfree Finite Element Methods: A Sensitivity Analysis for Design Trade Studies
A wheel experiencing sinkage and slippage events poses a high risk to planetary rover missions as evidenced by the mobility challenges endured by the Mars Exploration Rover (MER) project. Current wheel design practice utilizes loads derived from a series of events in the life cycle of the rover which do not include (1) failure metrics related to wheel sinkage and slippage and (2) performance trade-offs based on grouser placement/orientation. Wheel designs are rigorously tested experimentally through a variety of drive scenarios and simulated soil environments; however, a robust simulation capability is still in development due to myriad of complex interaction phenomena that contribute to wheel sinkage and slippage conditions such as soil composition, large deformation soil behavior, wheel geometry, nonlinear contact forces, terrain irregularity, etc. For the purposes of modeling wheel sinkage and slippage at an engineering scale, meshfree nite element approaches enable simulations that capture su cient detail of wheel-soil interaction while remaining computationally feasible. This study implements the JPL wheel-soil benchmark problem in the commercial code environment utilizing the large deformation modeling capability of Smooth Particle Hydrodynamics (SPH) meshfree methods. The nominal, benchmark wheel-soil interaction model that produces numerically stable and physically realistic results is presented and simulations are shown for both wheel traverse and wheel sinkage cases. A sensitivity analysis developing the capability and framework for future ight applications is conducted to illustrate the importance of perturbations to critical material properties and parameters. Implementation of the proposed soil-wheel interaction simulation capability and associated sensitivity framework has the potential to reduce experimentation cost and improve the early stage wheel design proce
terramechanics↗
Material flow behavior and microstructural refinement of AA6061 alloy during friction extrusion
In this work we used friction extrusion (FE), a solid phase processing technique, to produce dense, fully consolidated 5?mm rods of aluminum alloy 6061 (AA6061). The combination of large shear stresses and high temperatures at the tool-billet interface during extrusion produced equiaxed, dynamically recrystallized grains and finely distributed precipitates. Texture and microstructure evolution during extrusion was investigated in detail using scanning electron microscopy, electron backscatter diffraction, and transmission electron microscopy. Smoothed particle hydrodynamics simulation was performed to study complex material flow during extrusion. Simulation results suggest spiral material flow during extrusion which corroborated experimental results. Advantages of friction extruded microstructure over conventionally extruded counterparts are also explored using flash annealing for solution treatment followed by artificial aging. Mechanical properties of the as-friction extruded, and the artificially aged specimens were evaluated using tensile testing and compared with conventional extruded material properties.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Meshfree Process Modeling and Experimental Validation of Friction Riveting of Aluminum 5052 to Aluminum 6061
Friction riveting (Fric-riveting) is an innovative, fast, and energy-efficient process for spot-joining metal-metal structures. Although fric-riveting has been studied experimentally in recent years, its process modeling is rarely found in the literature primarily because of the associated large material deformation, extreme thermomechanical conditions, and complex contact conditions. In this work, a mesh-free smoothed particle hydrodynamics (SPH) framework that can well handle the abovementioned numerical challenges is used to simulate the fric-riveting of AA5052 to AA6061. Predicted material morphology, multi-point temperatures, and plunge force are thoroughly validated by experimental observations. The material severe plastic deformation zone in the vicinity of the riveting zone is further predicted by the SPH model, which indicates the material mixing and potential grain refinement zone. Based on the validated model, process parameters can be optimized which yields better performance over the baseline case.
Friction riveting (fric-riveting), Smoothed partic↗
Insights into the failure mode of the Chelyabinsk meteor from high-fidelity simulation
In 2013, an asteroid with a diameter of approximate 20 m broke up in the atmosphere above the Russian city of Chelyabinsk producing what is known as an airburst event. The event was recorded on video from multiple angles and a 500-kg fragment was recovered from Lake Chebarkul shortly thereafter. These aspects, coupled with the asteroid’s large size, distinguish the Chelyabinsk event as unique. In this study, we use the smoothed particle hydrodynamics method to examine the breakup mode of a Chelyabinsk-sized monolithic asteroid and probe whether a monolithic structure is consistent with observations from the Chelyabinsk event. We model the object as a homogeneous sphere with properties consistent with Chelyabinsk meteorites. We use planar-2D simulations to test the sensitivity of results to variations in the material model and to qualitatively inspect the cause of crack nucleation and fragmentation. We then compare results from a 3D simulation to observations from the event. Our 3D simulation predicts a burst height within 3 km of the observed peak luminosity suggesting the responsible asteroid may have been monolithic.
79 ASTRONOMY AND ASTROPHYSICS↗
An innovative and alternative approach toward gear fabrication
In this paper, we present the first demonstration of single-step gear manufacturing via an innovative in-situ friction stir forging approach. This novel process is a natural extension to the friction stir processing technique and relies on shear and normal stresses to form the complex shapes. Friction stir welding pin tool with shoulder is plunged at high revolutions per minute (rpm) into the sample to heat it to the desired temperature and enable plastic flow. Then, the tool is forged at high plunge rate to form the desired structurally-sound and defect free complex shapes. To better visualize the process the experiments are complemented by a 3D thermomechanically coupled, smoothed particle hydrodynamics (SPH) model to understand the material flow pattern both radially and through thickness in the gear teeth, plastic strain distribution, and temperature profile. We have demonstrated application of I-FSF to inherently poor formable AZ31 Mg, AA 7075-T6, and AA5083–10 vol% TiB 2 composite and fabricated 14-tooth spur gears from these materials with minimal farther machining. In conclusion, microstructure of the forged gears exhibited a directional material flow pattern due to the normal forging action and refinement of grains to 6–13 μm through dynamically recrystallization across the radial, through-thickness direction, which is in good agreement with the corresponding SPH simulation results.
42 ENGINEERING↗
Numerical simulation of asteroid geometry variance on airburst threat
For an atmospheric airburst the primary source of concern when assessing uncertainty is the size and velocity. Determining these properties provides the basis for threat assessment, as the total energy of the asteroid may then be estimated, and the threat investigated thoroughly. Even with clarity as to how much energy an asteroid may deposit, a great deal of uncertainty still exists for the actual energy deposition process. One such source of uncertainty is the geometry of the incoming asteroid. The geometry of an asteroid will alter the stress distribution during entry, which adds uncertainty to when fracture will occur. Here, in this study, we use Smoothed Particle Hydrodynamics to model the atmospheric airburst of Tunguska-scale asteroids with varying geometric profiles, including a sphere, ellipsoid, binary and superellipsoid. Each asteroid is modeled as a homogenous structure with strength. We assess uncertainty through a series of planar 2D simulation cases for each geometry, comparing the source of stochasticity across geometries. A single 3D airburst simulation for each geometry is also analyzed. Additionally, the 3D cases are compared to the highly uncertain Tunguska event, predicting variance in burst height across geometries, but all bounded by theoretical burst heights proposed for Tunguska.
Airburst↗
Modelling a Heaving Point-Absorber with a Closed-Loop Control System Using the DualSPHysics Code
The present work addresses the need for an efficient, versatile, accurate and open-source numerical tool to be used during the design stage of wave energy converters (WECs). The device considered here is the heaving point-absorber developed and tested by Sandia National Laboratories. The smoothed particle hydrodynamics (SPH) method, as implemented in DualSPHysics, is proposed since its meshless approach presents some important advantages when simulating floating devices. The dynamics of the power take-off system are also modelled by coupling DualSPHysics with the multi-physics library Project Chrono. A satisfactory matching between experimental and numerical results is obtained for: (i) the heave response of the device when forced via its actuator; (ii) the vertical forces acting on the fixed device under regular waves and; (iii) the heave response of the WEC under the action of both regular waves and the actuator force. This proves the ability of the numerical approach proposed to simulate accurately the fluid–structure interaction along with the WEC’s closed-loop control system. In addition, radiation models built from the experimental and WAMIT results are compared with DualSPHysics by plotting the intrinsic impedance in the frequency domain, showing that the SPH method can be also employed for system identification.
42 ENGINEERING↗
A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges
Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.
58 GEOSCIENCES↗
SPH modeling of biomass granular flow: Theoretical implementation and experimental validation
The commercialization of biomass-derived energy is impeded by flowability challenges arising from the feeding and handling of granular biomass materials in full-scale biorefineries. To overcome these obstacles, a robust and accurate model to simulate the flow of granular biomass is indispensable. However, conventional mesh-based numerical codes are limited by inherent mesh distortion in simulating large deformation that commonly occurs in granular biomass handling. Here, in this study, we propose a graphics processing unit (GPU)-accelerated meshless Smoothed Particle Hydrodynamics (SPH) code to model the flow of granular biomass materials. A modified void ratio-based mass conversation, a hybrid particle-to-particle/surface frictional boundary treatment, and a hypoplastic constitutive model are implemented. Four numerical examples, an elastic block sliding on inclined planes, sand column collapse, Angle of Repose, and axial compression tests for pine chips, were simulated using the developed SPH code. The results demonstrate good agreement between numerical predictions and analytical and experimental data for all four examples, validating the SPH code and increasing confidence that it can be applied to simulate more complex granular biomass handling processes, such as hopper feeding or auger conveyance.
09 BIOMASS FUELS↗
Simulating Hydrodynamics in Cosmology with CRK-HACC
Abstract We introduce CRK-HACC, an extension of the Hardware/Hybrid Accelerated Cosmology Code (HACC), to resolve gas hydrodynamics in large-scale structure formation simulations of the universe. The new framework couples the HACC gravitational N -body solver with a modern smoothed-particle hydrodynamics (SPH) approach called conservative reproducing kernel SPH (CRKSPH). CRKSPH utilizes smoothing functions that exactly interpolate linear fields while manifestly preserving conservation laws (momentum, mass, and energy). The CRKSPH method has been incorporated to accurately model baryonic effects in cosmology simulations—an important addition targeting the generation of precise synthetic sky predictions for upcoming observational surveys. CRK-HACC inherits the codesign strategies of the HACC solver and is built to run on modern GPU-accelerated supercomputers. In this work, we summarize the primary solver components and present a number of standard validation tests to demonstrate code accuracy, including idealized hydrodynamic and cosmological setups, as well as self-similarity measurements.
79 ASTRONOMY AND ASTROPHYSICS↗