Search NASA⌕ Search

SEARCH · Search NASA

Results for “Slip localization”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Thermoosmotic Transport in Nanochannels Grafted with pH-responsive Polyelectrolyte Brushes Modelled Using Augmented Strong Stretching Theory

In this paper, we develop a theory to establish that the thermoosmotic (TOS) effects, induced by the application of an axial temperature gradient, leads to a massive enhancement in liquid transport in nanochannels grafted with charged polyelectrolyte (PE) brushes. Here, we quantify the TOS transport by quantifying the induced electric field and the induced TOS flow field. The different components of the electric field, namely the ionic component, the thermal component, and the osmotic component, as well as the contributions of different ions to these components, are quantified. Furthermore, we express the TOS velocity as a combination of chemiosmotic (COS), thermal, and electroosmotic (EOS) components. The COS and the thermal components augment each other and the overall strength and direction of the TOS flow are dictated by the direction and the relative strength of the EOS component. Most importantly, we compare the cases of brush-grafted nanochannels with those of the brush-free nanochannels of identical surface charge densities: the TOS transport is massively augmented in the brush-grafted nanochannels attributed to the combination of the localization of the EDL (and hence any body force that depends on the EDL charge density) away from the nanochannel wall (i.e., the location of the maximum drag force) and the presence of a possible molecular slip (experienced by the liquid) along the brush surface.

42 ENGINEERING↗
DOE OSTI2023

Plastic homogeneity in nanoscale heterostructured binary and multicomponent metallic eutectics: An overview

Heterostructured materials comprised of relatively soft/hard disparate phases typically exhibit composite strengthening but lack plastic deformability at ambient temperatures. However, heterostructured systems comprised of nanoscale phases can simultaneously enhance yield strength and strain hardening, thereby promoting uniform distribution of plastic flow. In this review, the atomic-scale deformation mechanisms in model systems of eutectic alloys, Al-Al 2 Cu and Al-Si, refined to nanoscales via laser rapid solidification are discussed, and compared with literature on multi-component (high entropy) eutectics such as Ni-Al-Fe-based with Cr and/or Co additions. The nano-lamellar Al-Al 2 Cu structures exhibit unit defect mechanisms not reported in monolithic Al 2 Cu intermetallic: localized shear on {0 1 1} and shear-induced faults on {1 2 1} planes, constrained by closely-spaced dislocation arrays in Al confined by Al/Al 2 Cu interfaces. The unexpected plasticity mechanisms are enabled by slip continuity in nanoscale Al-Al 2 Cu eutectics associated with the orientation relationship and interface habit planes. In nano-fibrous Al-Si eutectic, tensile ductility at strength approaching 600 MPa is observed resulting from dislocation plasticity in the nano-Al channels and cracking in Si nanofibers. Molecular dynamics simulations show that Al dislocations easily cross-slip (screw) or climb (edge) along Al-Si interfaces, making slip transmission difficult. The propagation of nano-cracks is suppressed by surrounding strain hardening Al, retaining good ductility of the sample, in spite of lack of direct slip transmission. Finally, the critical unit mechanisms of slip transmission and interface-enabled plasticity observed in nanoscale eutectic binary systems can also explain the strength-ductility relationship in multi-component eutectics and homogeneously distributed plastic flow with increasing microstructural heterogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electron Nanocrystallography of Complex Materials and Processes

The proposed work is to develop Scanning CBED (SCBED) for the study of multi element crystals, including multiple principal elements alloys or high entropy alloys and relaxor-ferroelectric crystals. SCBED works like scanning transmission electron microscopy (STEM) but records the whole CBED at each probe position. Since CBED is a highly sensitive electron diffraction technique, SCBED combines the full benefits of CBED and STEM. Specifically, the SCBED based approach enables a quantitative determination of local atomic and electronic structures. The proposed work is thus to take advantage of SCBED for the study of complex multi-element crystals. The results from SCBED are expected to complement the methods of electron energy loss spectroscopy (EELS) near edge fine structure for probing local bonds, and X-ray and neutron scattering using the total scattering approach for probing local structure. Under the DOE BES support, we have developed SCBED to determine the structure of high entropy alloys, and we have coupled the work with in-situ compression study of dislocation dynamics to uncover the fundamental mechanisms of crystal slip.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Role of deformation twinning in fatigue of CrCoNi medium-entropy alloy at room temperature

Here cylindrical specimens of CrCoNi alloy with electropolished surfaces were subjected to constant total strain amplitude low cycle fatigue. The alloy exhibited an initial period of cyclic hardening followed by cyclic softening until failure occurred. At the end of hardening stage at the peak of cyclic stress, well-developed persistent slip markings (PSMs) consisting of extrusions and intrusions were associated with thin deformation twins. A sophisticated experimental workflow was designed to extract information from the surface and the bulk of tested material. A combination of SEM, EBSD, ECCI, FIB and HR-STEM was used to study the internal structure and the surface profiles around the deformation twins, which were produced during the initial period of cyclic loading. Furthermore, localized cyclic plastic strain and stress concentrations near deformation twins led not only to early, well-developed PSMs, but also to the activation of TWIP and TRIP plasticity even at low macroscopic stress amplitudes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of orientation-dependent deformation mechanisms in additively manufactured Zr using in-situ micromechanical testing: Twinning and orientation gradient

Here, orientation-dependent plasticity in zirconium (Zr)-alloy, which is produced through ultrasonic additive manufacturing (UAM) with subsequent hot-isostatic pressing (HIP), was analyzed by recording electron backscattered diffraction (EBSD) data during in-situ micromechanical testing. Three (0001) plane orientations of a hexagonal close-packed (HCP) structure were analyzed parallel to (1) the rolling direction [X||tensile direction (TD)], (2) the build direction [Z], and (3) the transverse direction [Y]. The analysis revealed that the grains with ~<0001>||TD show twin-dominant plasticity with three variants from {10 $\overline{1}$ 2}<$\overline{1}$ 011>; minor slipping with ($10\overline{1}$ 1)[$\overline{2}$ 113], and ($11\overline{2}$ 2)[$\overline{1} \overline{1}$ 23] pyramidal slip has also been observed. However, grains with orientation ~<0001>⊥TD are mainly sensitive to dislocations assisted plasticity-leading to orientation gradients formation. Furthermore, the neck formation was identified as originating from higher populated micro-crack locations and their association with localized plasticity at defect points in the UAM material. These results demonstrate that minimization of impurities to enable grain growth across prior foil interfaces makes HIP an effective methodology for producing Zr plate, with deformation characteristics expected for conventionally manufactured Zr.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Kink bands promote exceptional fracture resistance in a NbTaTiHf refractory medium-entropy alloy

Single-phase body-centered cubic (bcc) refractory medium- or high-entropy alloys can retain compressive strength at elevated temperatures but suffer from extremely low tensile ductility and fracture toughness. We examined the strength and fracture toughness of a bcc refractory alloy, NbTaTiHf, from 77 to 1473 kelvin. This alloy’s behavior differed from that of comparable systems by having fracture toughness over 253 MPa·m 1/2 , which we attribute to a dynamic competition between screw and edge dislocations in controlling the plasticity at a crack tip. Whereas the glide and intersection of screw and mixed dislocations promotes strain hardening controlling uniform deformation, the coordinated slip of <111> edge dislocations with {110} and {112} glide planes prolongs nonuniform strain through formation of kink bands. These bands suppress strain hardening by reorienting microscale bands of the crystal along directions of higher resolved shear stress and continually nucleate to accommodate localized strain and distribute damage away from a crack tip.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The effects of heat treatment on the response of WE43 Mg alloy: crystal plasticity finite element simulation and SEM-DIC experiment

The present report addressed the effects of heat treatment on the mechanical response of a WE43 Mg alloy using an integrated framework of SEM-DIC experiment and CPFE simulation. Both macroscopic responses including yield strength, ultimate strength, ductility, and microscopic responses, including local displacement and strain maps, were experimentally investigated. The focus of this work is to use the CPFE simulation as an integrated computational tool to study the effects of heat treatment. The CPFE framework was evaluated using the local fields of displacement and strain obtained from the SEM-DIC experiment rather than the conventional scheme of using macroscopic responses at different loading directions. Subsequently, the information which is available using CPFE, such as the critical resolved shear stress (CRSS) and relative slip activity, was used to study the effects of heat treatment on the response of WE43 Mg alloy. The contributions of different strengthening mechanisms on the CRSS were addressed. The results show that effects of heat treatment can be captured using the predominant mechanisms of the grain size effect and the influence of precipitates. Finally, it has been shown that classical Hall-Petch in which one constant can capture the size effects, should be modified. To do so, each deformation mode should have a unique Hall-Petch constant, which are calculated here for the WE43 Mg alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Anomalous size effect on yield strength enabled by compositional heterogeneity in high-entropy alloy nanoparticles

High-entropy alloys (HEAs), although often presumed to be random solid solutions, have recently been shown to display nanometer-scale variations in the arrangements of their multiple chemical elements. Here, we study the effects of this compositional heterogeneity in HEAs on their mechanical properties using in situ compression testing in the transmission electron microscope (TEM), combined with molecular dynamics simulations. We report an anomalous size effect on the yield strength in HEAs, arising from such compositional heterogeneity. By progressively reducing the sample size, HEAs initially display the classical “smaller-is-stronger” phenomenon, similar to pure metals and conventional alloys. However, as the sample size is decreased below a critical characteristic length (~180 nm), influenced by the size-scale of compositional heterogeneity, a transition from homogeneous deformation to a heterogeneous distribution of planar slip is observed, coupled with an anomalous “smaller-is-weaker” size effect. Atomic-scale computational modeling shows these observations arise due to compositional fluctuations over a few nanometers. These results demonstrate the efficacy of influencing mechanical properties in HEAs through control of local compositional variations at the nanoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dislocation‐Driven Formation of Oriented Macroperiodic Metastructures of Curved Single Crystal Lattices in Glass

Abstract Single crystals fabricated in glass by localized heating can develop uniquely deformed lattices stabilized by the surrounding amorphous medium. The development of lattice curvature appears to be intrinsic to the crystal growth process in some systems, while the result of the locally changing crystallography in others. In this work, a model laser‐fabricated rotating lattice Sb 2 S 3 crystal grown in stoichiometric glass is used to demonstrate fabrication of novel macroperiodic metastructures that utilize intrinsic lattice curvature superimposed with subtle crystallographic influences. The limited availability of slip systems drives the lattice curvature magnitude to vary with crystal growth direction, maximizing for lattices aligned with the predominant Burgers vector along with corresponding increases in dislocation density. Misaligned lattice orientations form smaller secondary lattice curvatures arising from misaligned Burgers vectors with further elastic contributions. Over extended crystal growth, these secondary components align the lattice to rotate about either the <001> or <010> crystal axes forming repeating metastructures of lattice orientation with periodicity 20–160 microns in length. The mechanistic approach used in this work may be expanded to other systems with known slip systems to better understand and design macroperiodic metastructures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering↗
DOE OSTI2023

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗
DOE OSTI2023

Quaternary slip rates and most recent surface rupture of the Bitterroot fault, western Montana

The Bitterroot fault is a ~100-km-long Quaternary active normal fault that bounds the eastern margin of the north–south-trending Bitterroot Mountains and accommodates extension near the Intermountain Seismic Belt. Here, new detailed mapping using high-resolution topographic data derived from light detection and ranging (lidar) along the southern Bitterroot Range documents multiple generations of fault scarps in Holocene–Pleistocene deposits with vertical off sets that increase in magnitude with age. Fault mapping indicates a complex fault geometry characterized by an en echelon pattern of discontinuous segments of 41–78° east-dipping normal faults that appear to cut the older Eocene detachment fault, and locally 70–88° west-dipping antithetic normal faults. 10 Be cosmogenic exposure dating provides in situ age control for 32 surface boulders (>1 m) sampled in glacial deposits. Near Como Dam, two Pinedale-age glacial moraine sequences yield peak age distributions of 15.0 ± 0.4 ka and 16.4 ± 0.6 ka as apparent exposure ages (e = 0), and 15.4 ± 0.4 ka and 16.8 ± 0.6 ka based on the maximum allowed boulder surface erosion rate (e = 2 mm/ka). Vertical separation of 3.5 ± 0.2 m across the ~16–17 ka glacial moraine off set by the Bitterroot fault scarp yields a fault slip rate of 0.2–0.3 mm/yr. Glacial Lake Missoula highstand shorelines, inset into the ~15 ka glacial moraine and vertically off set 4.6 ± 1.6 m by an antithetic strand of the Bitterroot fault, yield fault slip rates of 0.2–0.5 mm/yr that overlap with fault slip rates on the main strand near Lake Como (0.2–0.3 mm/yr). At the Ward Creek fan located ~15 km to the north of Lake Como, two glacial debris fan sequences yield peak age distributions of 16.6 ± 0.4 ka and 62.8 ± 1.7 ka (e = 0), and 17.0 ± 0.4 ka and 69.9 ± 2.2 ka (e = 2 mm/ka). Vertical separations of 2.4 ± 0.3 m and 4.5 ± 0.2 m on the ~17 ka and ~63–70 ka fan surfaces off set by the Bitterroot fault yield fault slip rates of 0.2–0.3 mm/yr and 0.1 mm/yr, respectively. Our results indicate broadly consistent fault slip rates for the main fault segments at Lake Como (0.2–0.3 mm/yr) and the Ward Creek fan (0.1–0.3 mm/yr) with an along-strike range of 0.1–0.3 mm/yr for the southern Bitterroot fault. Fault scaling relations and evidence of multiple late Quaternary fault surface ruptures suggest the Bitterroot fault could produce a Mw ~7.2 earthquake. Structural model constraints and our slip rate results indicate both high-angle or low-angle fault geometries are possible at depth. A seismogenic low-angle fault model could generate a larger earthquake of Mw >7.2. Earthquake history is unknown for the Bitterroot fault, but fault scarps in young glacial deposits demonstrate its seismogenic potential. Data from this study suggest seismic hazards from the Bitterroot fault may pose a high level of risk to the Missoula metropolitan area, the State’s second most populous region, and major infrastructures across the Missoula and Bitterroot Valleys.

58 GEOSCIENCES↗
DOE OSTI2022

Origin, parameters, and underlying deformation mechanisms of propagating deformation bands in irradiated 316L stainless steel

Lüders-type propagating deformation bands were observed in specimens of irradiated 316L stainless steel samples removed from Spallation Neutron Source target vessels after service. Mechanical testing with digital image correlation (DIC) and in-situ tensile testing with scanning electron microscopy electron backscatter diffraction showed that the observed Lüders-type behavior was not related to the known transformation-induced plasticity or twinning-induced plasticity behavior. Instead, the phenomenon occurs at small local strain values before a significant amount of martensite or deformation twins appear in the microstructure. Microstructural analysis and in-situ mechanical test results suggest Lüders-type band formation and propagation were related to the appearance and evolution of defect-free channels—analogous to slip bands. A modified Swift equation with a Ludwigson-like component was offered to rationalize the phenomenon and model the strain-softening processes at small strain values. Finally, the results indicate complex microstructural processes and deformation mechanisms were active at small strain values and underline the benefits of advanced mechanical test approaches such as DIC.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deformation and fracture characteristics of zirconium plate produced via ultrasonic additive manufacturing

Abstract The microstructural evolution, deformation modes, and fracture mechanisms of zirconium plate produced using ultrasonic additive manufacturing (UAM) are presented. In addition to conventional tensile testing techniques, digital image correlation captured highly variable strain accumulation in specimens loaded perpendicular or parallel to the build height (Z). When tested in parallel to Z, delamination at prior foil/foil interfaces creates strain localization noticeable in strain rate maps, whereas specimens loaded perpendicular to Z illustrate conventional strain hardening until necking accelerates delamination. Although bond strengths are statistically and spatially variable, in situ electron backscattering diffraction tests illustrate the ability for grains near interfaces to accommodate strain with twinning and slip modes consistent with conventionally produced zirconium alloys. Finally, mixtures of ductile and delamination-induced fracture highlight the interface-driven failure modes of UAM zirconium plate in the as-built condition. Graphic abstract

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fingerprinting shock-induced deformations via diffraction

Abstract During the various stages of shock loading, many transient modes of deformation can activate and deactivate to affect the final state of a material. In order to fundamentally understand and optimize a shock response, researchers seek the ability to probe these modes in real-time and measure the microstructural evolutions with nanoscale resolution. Neither post-mortem analysis on recovered samples nor continuum-based methods during shock testing meet both requirements. High-speed diffraction offers a solution, but the interpretation of diffractograms suffers numerous debates and uncertainties. By atomistically simulating the shock, X-ray diffraction, and electron diffraction of three representative BCC and FCC metallic systems, we systematically isolated the characteristic fingerprints of salient deformation modes, such as dislocation slip (stacking faults), deformation twinning, and phase transformation as observed in experimental diffractograms. This study demonstrates how to use simulated diffractograms to connect the contributions from concurrent deformation modes to the evolutions of both 1D line profiles and 2D patterns for diffractograms from single crystals. Harnessing these fingerprints alongside information on local pressures and plasticity contributions facilitate the interpretation of shock experiments with cutting-edge resolution in both space and time.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Role of interfaces on the mechanical response of accumulative roll bonded nanometallic laminates investigated via dislocation dynamics simulations

Unraveling the effects of continuous dislocation interactions with interfaces, particularly at the nanometer length scales, is key to a broader understanding of plasticity, to material design and to material certification. To this end, this work proposes a novel discrete dislocation dynamics-based model for dislocation interface interactions tracking the fate of residual dislocation on interfaces. This new approach is used to predict the impact of dislocation/interface reactions on the overall mechanical behavior of accumulative roll bonded nanometallic laminates. The framework considers the dynamic evolution of the interface concurrent with a large network of dislocations, thus, accounting for the local short and long range effects of the dislocations under the external boundary conditions. Specifically, this study focuses on two-phase Fe/Cu nanometallic laminates, and investigates the role of the underlying elastic and plastic contrast of the Fe and the Cu layers on the composite response of the material. Moreover, the role of initial microstructures, resulting from processing is also investigated. Subsequently, the model is used to examine the effect of layer thickness and interface orientation relationship on the residual stresses of the relaxed microstructure. The associated mechanical response of these laminates are compared when loaded under normal direction compression, as well as shear compression. Finally, this work predicts a dominant effect of the layer thickness, as compared to the interface orientation relationship, on the macroscopic response and on the residual stresses of these nanolaminates, while the local dislocation transmission propensity through the interface is significantly influenced by the corresponding orientation relationship.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ordering effects on deformation substructures and strain hardening behavior of a CrCoNi based medium entropy alloy

Here, a CrCoNi based medium entropy alloy with small additions of Ti, Al and Nb (denoted as (CrCoNi) 93 Al 4 Ti 2 Nb) in the as-quenched condition, exhibits tensile properties comparable to those of the equiatomic CrCoNi alloy at room temperature. Dark field transmission electron microscopy (TEM), atomic resolution high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) together with atom probe tomography (APT) show that spatially-localized long range ordering (LRO) L1 2 domains exist in this alloy. The evolution of deformation substructure with plastic deformation in this alloy was characterized using electron backscatter diffraction (EBSD), electron channeling contrast imaging (ECCI) and STEM based techniques including the recently developed weak beam dark field STEM imaging. Plastic deformation occurs by the slip of a/2<110>dislocations, which are narrowly dissociated into Shockley partial dislocations on {111} slip planes. Their dissociation distances in the (CrCoNi) 93 Al 4 Ti 2 Nb alloy are much smaller than the widths of the corresponding partials in the equiatomic CrCoNi alloy due to one or more of the minor alloying elements (Al, Ti, Nb). Dislocation slip in this alloy has a pronounced planar character. The leading dislocations in slip bands glide as pairs due to the existence of LRO domains. Multipoles were formed through the slip of dislocations with opposite signs on adjacent {111} slip planes. Those multipoles serve as building blocks for the formation of subgrain structures consisting of fine slip bands. The distances between slip bands were continuously refined during plastic deformation and dynamic refinement of slip bands plays a crucial role in strain hardening. The effects of LRO domains on planar dislocation slip, the deactivation of deformation twinning and strain hardening of this alloy are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗