Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solidification”

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 1,099 records · Page 61

The microstructure of a selective laser melting (SLM)-fabricated NiTi shape memory alloy with superior tensile property and shape memory recoverability

A selective laser melting (SLM)-fabricated NiTi with superior tensile property and shape memory recoverability was obtained by using a unique stripe rotation scanning strategy. Here, the alloy characteristics, formation mechanisms and evolution in terms of twins, dislocations and precipitations of the alloy were systematically studied. Compared to the conventional smelting-followed-by-machining, the SLM fabrication process involves rapid solidification and repeated heating, which confer distinctive characteristics to the microstructures of SLM-fabricated NiTi alloys. Rapid solidification promotes the formation of a supersaturation solid solution matrix containing a high concentration vacancies, which in turn aggregate to generate a high density of dislocations. During subsequent repeated heating stages, these dislocations occur thermal motion along three directions of <001 > , <111> and <110 >, leading to the formation of thermal kinks, helical dislocations and wave morphology. Simultaneously, precipitated particles Ti 3 Ni 4 repeatedly nucleate and heterogeneous grow with the movement of dislocation. Such precipitation behavior, termed repeated precipitation, has not been previously reported in the conventional NiTi alloys, suggesting that it could be a unique characteristic of such alloy. After martensitic transformation, only two twins, {1 1¯ 1} type I twin and compound twin, are detected. The twinning lamellae of these two twins, where precipitation and dislocation pile-ups exist, often have uneven thickness and chaotic arrangement. Besides, the unique self-accommodated microstructures, such as secondary {1 1¯ 1} type I twin and compound twin with “herring-bone’’ lamellae, which often appear in the deformed or nanocrystalline NiTi, can also be observed. These unique microstructures may confer the distinctive properties to the SLM-fabricated NiTi.

36 MATERIALS SCIENCE↗

A sensitivity analysis of the columnar-to-equiaxed transition for Ni-based superalloys in electron beam additive manufacturing

Solidification conditions experienced during Selective-Electron Beam Melting (S-EBM) of Nickel-Base superalloys generally lead to epitaxial grain growth from previous layers, resulting in highly textured columnar grain structures, and consequently, highly anisotropic mechanical properties. However, in recent work, researchers have been able to produce equiaxed grain structures by changing the scan strategy to influence the solidification conditions, thereby obtaining isotropic properties. The competition between columnar and equiaxed grain formation is governed by both the local thermal conditions and alloy thermodynamics. The above examples demonstrate successful manipulation of the thermal conditions through changes in process parameters. However, there is no guarantee that process conditions that encourage the formation of equiaxed grains for a given alloy will coincide with the production of defect free material. The goal of this work is therefore to understand the influence of alloy composition on the columnar-to-equiaxed transition (CET) so that alloys may be designed such that grain structure control is accessible to a wider range of S-EBM process conditions. A CET model originally developed by Hunt and extended by Gäumann was used to relate processing conditions to resulting microstructure. A sensitivity analysis of the CET model was performed using the Element Effects Methods. Furthermore, the results lead to a comparison of the relative effects of various alloying additions on the CET, as well as an assessment of other model parameters, including, most significantly, the nucleation density of equiaxed grains.

36 MATERIALS SCIENCE↗

Toucan: A performance portable, scalable implementation of the DECA algorithm

In the field of additive manufacturing (AM), cellular automata (CA) is extensively used to simulate microstructural evolution during solidification. However, while traditional CA approaches are relatively fast, they still require a substantial number of time steps, are limited to moderate volumes, and are relatively difficult to improve through parallelism due to the highly localized nature of the solidification front. Here, to address these issues of time to solution and load balancing, we introduce Toucan, a parallel, performance-portable, and scalable code written in C++ with the Kokkos library that leverages the discrete event inspired cellular automata (DECA) algorithm to perform parallel-in-time (PinT) grain growth simulations. Toucan effectively mitigates load balancing issues by distributing the computational workload more evenly across processors, enhancing scalability and efficiency. We conduct both strong and weak scaling studies on up to 64 GPUs on the Frontier supercomputer, demonstrating that Toucan significantly outperforms the current state-of-the-art, time-stepped CA code, ExaCA, on both single and multi-GPU simulations. Even in AM-specific weak scaling scenarios, Toucan maintains near-ideal scaling, in contrast to the linear increase observed with ExaCA due to the moving laser raster pattern. This study highlights Toucan’s potential to transform microstructural simulations in AM by radically improving both efficiency and scalability over existing methods.

36 MATERIALS SCIENCE↗

Isotopic evidence for a young lunar magma ocean

Herein, mare basalt sources and ferroan anorthosite suite cumulates define a linear array on a 146 Sm/ 144 Nd versus 142 Nd/ 144 Nd isochron plot demonstrating these materials were derived from a common reservoir at 4336 +31 / –32 Ma. The minimum proportion of the Moon that was in isotopic equilibrium at this time is estimated to be 1-3% of its entire volume based on the geographic extent from which the analyzed samples were collected and the calculated depths from which the samples were derived. Scenarios in which large portions of the Moon were molten to depths of many hundreds of kilometers are required to produce the observed Sm-Nd isotopic equilibrium between the mantle and crustal rocks at 4.34 Ga. This is a consequence of the fact that limited heating of a solid Moon above the blocking temperature of the Sm-Nd isotopic system is insufficient to diffusively homogenize radiogenic Nd throughout the mantle and crust. There are three scenarios that might account for global-scale isotopic equilibrium on the Moon relatively late in Solar System history including: (1) Sm-Nd re-equilibration of a solid Moon resulting from widespread melting in response to mantle overturn or a very large impact, (2) early accretion of the Moon followed by delayed cooling due to the presence of an additional heat source that kept a large portion of the Moon molten until 4.34 Ga, or (3) late accretion of the Moon followed by rapid cooling of the magma ocean late in Solar System history. Neither density-driven overturn of the mantle, nor a large impact, are likely to homogenize the mantle and crust to the extent required by the Sm-Nd isochron. Likewise, secondary heating mechanisms, such as tidal heating or radioactive decay, are not efficient enough to keep the Moon molten to the depth of the mare basalt source regions for many tens to hundreds of millions of years. Instead, the age of equilibrium between such a compositionally diverse set of rocks, produced on a global scale, likely records the time of primordial solidification of the Moon from a magma ocean. This scenario accounts for both the petrogenetic characteristics of lunar rock suites, as well as their Sm-Nd isotopic systematics. It is supported by the preponderance of ~4.35 Ga ages obtained for other hypothetical magma ocean crystallization products, such as ferroan anorthosite suite rocks and K, REE, and P enriched cumulates that are thought to represent flotation cumulates of the magma ocean and the last vestiges of magma ocean solidification, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermodynamic and kinetic analysis of the melt spinning process of Fe-6.5 wt.% Si alloy

Here, the microstructural evolution of Fe-6.5 wt.% Si alloy during rapid solidification was studied over a quenching rate of 4 × 10 4 K/s to 8 ×10 5 K/s. The solidification and solid-state diffusional transformation processes during rapid cooling were analyzed via thermodynamic and kinetic calculations. The Allen-Cahn theory was adapted to model the experimentally measured bcc_B2 antiphase domain sizes under different cooling rates. The model was calibrated based on the experimentally determined bcc_B2 antiphase domain sizes for different wheel speeds and the resulting cooling rates. Good correspondence of the theoretical and experimental data was obtained over the entire experimental range of cooling rates. Along with the asymptotic domain size value at the infinite cooling rates, the developed model represents a reliable extrapolation for the cooling rate > 10 6 K/s and allows one to optimize the quenching process.

36 MATERIALS SCIENCE↗

Surface oxidation and crystallization of FeNi-Based soft magnetic nanocrystalline and amorphous nanocomposite alloys

In this work, we present characterization of surface oxidation and crystallization behavior of Metal Amorphous Nanocomposite (MANC) (Fe 70 Ni 30 ) 80 Nb 4 Si 2 B 14 alloys via complementary x-ray diffraction, electron microscopy, and x-ray photoelectron spectroscopy techniques. An adherent amorphous oxide of roughly 15–20 nm in thickness is observed on the wheel side of the ribbon, enriched in Fe relative to Ni and with significant content of Si and Nb as would be expected based upon a comparison of the Gibbs free energy of formation of the stable oxide phases. The amorphous character of the oxide is likely to result from the significant Si-content coupled with the amorphous and nanocrystalline character of the underlying metallic layer. Higher density of FeNi-enriched nanocrystals is also observed within ~100 nm of the wheel-side surface, presumably due to reduction in local glass former concentration. A rough, crystalline copper oxide layer is also present and forms through contact between the copper wheel and the solidifying melt during the rapid solidification process. On the air side of the ribbon, a crystalline oxide layer of ~15–20 nm is formed as well as intergranular oxidation that penetrates along the grain boundaries between highly textured surface crystallites of ~50–100 nm in size. The crystalline nature of the oxide on these surfaces is asserted to result from the lack of significant Si-content and the pre-existing surface crystallites that form during the rapid solidification process which promote a crystalline oxide surface layer. Furthermore, these results shed new insights into the details of surface oxidation and can help to inform processing and alloy composition adjustments which seek to tailor surface properties.

36 MATERIALS SCIENCE↗

Spatial mapping of the localized corrosion behavior of a magnesium alloy AZ31B tungsten inert gas weld

Sections of a magnesium alloy, AZ31B, joined with tungsten inert gas (TIG) welding, were examined with scanning electrochemical microscopy (SECM) and scanning Kelvin probe force microscopy (SKPFM) to investigate corrosion mechanisms by correlating observed corrosion behavior with weld-affected microstructural variations. Insight into the changing nature of the galvanic couples between weld zones and at localized microgalvanic sites were investigated using SECM and SKPFM to map both electrochemically active regions and Volta potential differences across the weld-affected zones. The formation of an Al-Zn solidification network in the fusion zone (FZ) at and near the TIG weld epicenter differs from the outer heat-affected zone (HAZ), where intermetallic particles (IMPs) are the notable secondary phase from the magnesium matrix. These microstructures were mapped with SKPFM before and after brief exposure to a salt solution, revealing micro-galvanic couples as the main driving force to corrosion initiation and propagation within each zone. The IMPs and Al-Zn solidification network act as strong cathodes and govern the corrosion processes. The galvanic coupling and evolution of the intrinsic corrosion behavior between the weld zones is explained by monitoring the hydrogen evolution reaction (HER) with SECM over time. Anodically induced cathodic activation is confirmed for this welded material, as micro-galvanic couples between microstructural features are found to transition over time to broad electrochemically active areas within the weld-affected zones, resulting in polarity reversal as time of exposure proceeds.

Anodically induced cathodic activation↗

Additive manufacturing process parameter determination for a new Fe-C-Cu alloy

As a potential replacement for stainless steel alloys commonly used to print parts with laser powder bed fusion, a new Cu precipitation strengthened ferrous alloy, with composition Fe-0.2C-6Cu (wt%), was recently developed. This material is Co- and Ni-free, printable, has mechanical properties comparable to that of high strength stainless steels (approximately 1300 MPa UTS), and is cost-advantaged relative to existing low-alloy steels for additive manufacturing of parts with complex shapes. To gain traction for broader application, optimal laser powder bed fusion (LPBF) processing conditions to produce nearly defect-free parts without compromising mechanical strength are needed. Here, an optimal processing window based upon laser speeds that achieve minimal porosity was quantified via comparisons of measured melt pool penetration, scan speeds, printed material density and laser beam energy density on a commercial LPBF unit operating at 350 W, 80 μm hatch spacing, and 50 μm layer thickness based upon optimal parameters for 17-4PH. Furthermore, the window is 300 to 500 mm/s to achieve >99.8 % dense Fe-0.2C-6Cu (wt%) specimens. Other defects such as burning, spatter, and lack of fusion were also avoided in this window. The window was validated with in-situ synchrotron X-ray imaging that enabled visualization of the vapor cavity, melting, and solidification during a single laser track scan on a miniature powder bed sample. In addition, in-situ infrared imaging provided temperature fields, cooling rate and solidification range, and confirmed minimal printing defects and spatter within the processing window.

36 MATERIALS SCIENCE↗

Localized orientation gradients in additively manufactured stainless steel 316H structures

The high solidification rates during additive manufacturing cause highly localized thermal and strain gradients. The effect of these gradients on the evolution of local orientation misorientations within a grain is not well understood. In this study, stainless steel 316H parts were fabricated via laser powder bed fusion using three different energy densities: 43, 71, and 135 J/mm 3 . Electron backscatter diffraction showed that the maximum misorientations of the grains can be up to 25° along the build direction. Misorientation gradients (RM g ) within grains are process-dependent and can change from 0.036°/μm to 0.015°/μm with increased volumetric energy densities. In conclusion, the characterized misorientation gradients are an indication of the level of dislocations and, to an extent, the plastic deformation resulting from the rapid solidification during laser powder bed fusion.

316H↗

Multimodal and multiscale strengthening mechanisms in Al-Ni-Zr-Ti-Mn alloy processed by laser powder bed fusion additive manufacturing

The unique thermokinetics of laser-powder bed fusion additive manufacturing (L-PBFAM) has been exploited for development of a novel high-strength Al-Ni-Ti-Zr-Mn alloy. The addition of 0.5 wt% Mn leads to extraordinary improvement in ultimate tensile strength (502 MPa) and work hardening due to the activation of two Mn-induced strengthening mechanisms. First, by a bimodal particle strengthening effect due to Al 31 Mn 6 Ni 12 nano-quasi-crystal particles rejected in inter-dendritic spaces and fibrous Al 3 Ni eutectic dendritic channels, which predominately contributes to the strength improvement, and second by solid solution strengthening from remaining Mn entrapped in Al. These mechanisms supplement the particle strengthening effect imparted by coherent and incoherent Al 3 (Ti,Zr) co-precipitates present at melt pool boundaries, dislocation strengthening due to solidification induced strain, and Hall-Petch and backstress strengthening effect due to heterogenous grain size distribution occurring at various length scales. In conclusion, the solidification dynamics and hierarchical heat distribution that are associated with L-PBFAM resulted in complex spatial variations in these strengthening phenomena and were investigated via a high-throughput multiscale structure–property correlation that involved thermokinetic simulation, X-ray diffraction, high-resolution nanoindentation mapping, and site-specific transmission electron microscopy of the alloy.

36 MATERIALS SCIENCE↗

On the evolution of microstructure and defect control in 316L SS components fabricated via directed energy deposition

To identify the critical issues that affect the evolution of microstructure during additive manufacturing, we investigated the influence of process parameters on the evolution of the dimensional and surface quality, microstructure, internal defects, and mechanical properties in 316L stainless steel (SS) components fabricated using laser engineered net shaping (LENS®), a directed energy deposition (DED) additive manufacturing (AM) technique. The results show that the accumulation of un-melted powder particles on the side walls of deposited sections can be avoided by selecting a laser under-focused condition. Moreover, we report that the variation of melt pool width is more sensitive to laser power than to the depth of the melt pool. The formation of a so-called “hierarchical” microstructure with cellular morphology is attributable to a combination of layer deposition and rapid solidification, which are characteristics of AM. Finally, we discuss microstructure evolution and defect formation, particularly the formation of multiple interfaces and the presence of un-melted powder particles and pores, in light of the dynamic convective fluid flow and rapid solidification that occur in the melt pool. X-ray computed tomography (X-CT) was used to precisely map the spatial distribution of pores in the DED components. The evolution of microstructure during DED is discussed in the context of related thermal phenomena in an effort to provide fundamental insight into the mechanisms that govern defect formation.

36 MATERIALS SCIENCE↗

Microstructure and creep properties of cast near-eutectic Al–Ce–Ni alloys

This study investigates the as-cast and aged microstructures, thermal stability, ambient temperature strengthening, and creep resistance of three ternary Al–Ce–Ni alloys (wt%): near-eutectic Al–10Ce–5Ni (with both eutectic and hypoeutectic regions), hypoeutectic Al-7.5Ce-3.75Ni (with numerous primary Al dendrites), and hypereutectic Al-12.5Ce-6.25Ni (with coarse, blocky primary Al3Ni and Al11Ce3 precipitates and some primary Al dendrites). Depending on the alloy composition and local solidification conditions, the following eutectic morphologies are found: (i) coarse Al–Ce eutectic colonies where Al 11 Ce 3 is in the form of “Chinese script”, (ii) intermingled regions of binary Al–Ce and Al–Ni eutectic colonies, with finer Al 3 Ni and Al 11 Ce 3 fibers, (iii) large ternary eutectic colonies, where the binary Al 3 Ni and Al 11 Ce 3 phases are alternating or intertwining within the individual, fine fibers (diameters of ~60–170 nm, depending on solidification rates), and (iv) ternary eutectic zones (between primary Al dendrites), where fine Al 3 Ni and Al 11 Ce 3 build up a 3D-interconnected network. The high volume fraction of intermetallic phases and extremely fine eutectic spacing/fiber diameter both contribute to high ambient strengthening (higher as-cast microhardness than binary Al–Ce or Al–Ni), and also provide enhanced creep resistance at 300 and 350 °C. Additionally, the alloys are coarsening-resistant up to 425 °C for extended periods, with a gradual decrease in microhardness. The alloys aged at 400 °C to 1050 h show fiber fragmentation and coarsening of the resulting particles, with the faster-diffusing Ni driving more rapid coarsening of the Al 3 Ni particles which engulf finer, more stable Al 11 Ce 3 particles. Severe overaging (performed at 590 °C for 24 h) leads to Al 3 Ni and Al 11 Ce 3 spheroids which remain submicron-sized in eutectic colonies, but micron-sized at colony boundary and at Al dendrite-eutectic interface. Creep resistance at 300 °C of overaged Al–10Ce–5Ni remains substantial, consistent with load-transfer based composite strengthening being an important strengthening mechanism in these alloys, making them excellent candidates for replacement of heavier steel or titanium parts operating under stress up to 300 °C.

36 MATERIALS SCIENCE↗

Microstructural heterogeneities in additively manufactured refractory alloy C103 and their implications for room and elevated temperature mechanical behavior

Considering the vast component design space enabled by fusion-based additive manufacturing (F-BAM) processes, e.g., directed energy deposition (DED), the scale-up manufacturing of Nb-alloys with F-BAM is advantageous for structural applications. However, varying thermokinetic parameters-induced microstructural heterogeneities are prevalent within the F-BAM processed alloys. Such microstructural heterogeneities can have significant implications for the room and elevated temperature mechanical behavior. While a few studies investigating F-BAM processed alloy C103 are available, none of these studies investigate the microstructural heterogeneities – including those associated with solidification growth modes and second phase particles – and the effect thereof on the mechanical behavior. To this end, we investigate the microstructurally heterogeneous regions with varying solidification growth morphologies, segregation behavior, and second phase particle attributes within the laser-DED processed alloy C103. The implications of such heterogeneous regions for room- and elevated-temperature tensile behavior and damage mechanisms are revealed. Particularly, the interface between the cellular and planar region is identified as susceptible to deformation localization. The implications of hot isostatic pressing (HIP) for the consolidation behavior, microstructural evolution, and resulting mechanical behavior are also discussed. Although the recrystallization and grain growth led to a reduced yield strength in the HIPed condition, the homogenization of microstructure alleviated the deformation localization sites, such as the planar/cellular interface within the melt pool. The homogenized microstructure alongside the enhanced consolidation upon HIP led to an enhanced elongation to failure. Findings establish microstructural design considerations in F-BAM processed Nb alloys and also facilitate design of post-processing heat treatments for achieving improved mechanical properties.

36 MATERIALS SCIENCE↗

Microstructural evolution, defect mitigation, and precipitation behavior in AA6061 via laser powder bed fusion with high-temperature substrate heating

Defects, particularly solidification cracking, remain persistent challenges in the laser powder bed fusion (LPBF) processing of AA6061 aluminum alloy. This study systematically investigates defect mitigation, microstructural evolution, and mechanical properties associated with high-temperature substrate preheating at 500 °C. Comprehensive microstructural analyses, including characterization of defects, grain structures, and precipitation behavior, were performed on samples in both as-built and T6 heat-treated states. Elevated preheating substantially reduced solidification cracking across a wide processing window, while demonstrating decreased crack sensitivity to laser parameters. Columnar cracks along the build direction were observed despite substrate preheating. Lack-of-fusion and keyhole porosity were effectively eliminated, though gas-induced microporosity persisted at higher powers. In-depth characterization of two distinct laser power and speed conditions confirmed the formation of micron-sized, non-coherent Mg 2 Si precipitates under heated substrate conditions, alongside α-AlFeCrMnSi intermetallic phases indicative of in-situ thermal effects during fabrication. Subsequent T6 heat treatment revealed the formation of fine, coherent needle-shaped β″ strengthening precipitates. Despite substantial differences in processing parameters, comparable mechanical properties were measured in the as-built samples (∼52 MPa yield strength, ∼130 MPa tensile strength), primarily due to reduced strain hardening effects and consistent precipitation characteristics. Meanwhile, the T6 heat treatment led to significant improvement in properties, enhancing yield strength by over 400%, aligning closely with the performance of conventional wrought AA6061-T6. These findings underscore that high-temperature substrate preheating offers an effective means to suppress cracks, control precipitation, and enhance mechanical performance in LPBF-processed AA6061.

36 MATERIALS SCIENCE↗

Model for the Solid–Liquid Interfacial Free Energy at High Pressures

The free energy involved in the formation of an interface between two phases (e.g., a solid–liquid interface) is referred to as the interfacial free energy. For the case of solidification, the interfacial free energy dictates the height of the energy barrier required to nucleate stable clusters of the newly forming solid phase and is essential for producing an accurate solidification kinetics model using classical nucleation theory (CNT)-based methods. While various methods have been proposed for modeling the interfacial free energy for solid–liquid interfaces in prior literature, many of these formulations involve making restrictive assumptions or approximations, such as the system being at or near equilibrium (i.e., the system temperature is approximately equal to the melt temperature) or that the system is at pressures close to atmospheric. However, these approximations and assumptions may break down in highly non-equilibrium situations, such as in dynamic-compression experiments where metastable liquids that are undercooled by hundreds of kelvin or over-pressurized by several gigapascals or more are formed before eventually solidifying. We derive a solid–liquid interfacial free-energy model for such high-pressure conditions by considering the enthalpies of interactions between pairs of atoms or molecules. We also consider the contribution of interface roughness (disordering) by incorporating a multilayer interface model known as the Temkin n-layer model. Our formulation is applicable to a diverse variety of materials, and we demonstrate it by developing models specifically for two different materials: water and gallium. We apply our interfacial free-energy formulation to CNT-based kinetics simulations of several suites of dynamic-compression experiments that cause liquid water to solidify to the high-pressure solid polymorph ice VII and have found good agreement to the observed kinetics with only minor empirical fitting.

36 MATERIALS SCIENCE↗

Thermodynamics-guided alloy and process design for additive manufacturing

In conventional processing, metals go through multiple manufacturing steps including casting, plastic deformation, and heat treatment to achieve the desired property. In additive manufacturing (AM) the same target must be reached in one fabrication process, involving solidification and cyclic remelting. The thermodynamic and kinetic differences between the solid and liquid phases lead to constitutional undercooling, local variations in the solidification interval, and unexpected precipitation of secondary phases. These features may cause many undesired defects, one of which is the so-called hot cracking. The response of the thermodynamic and kinetic nature of these phenomena to high cooling rates provides access to the knowledge-based and tailored design of alloys for AM. Here, we illustrate such an approach by solving the hot cracking problem, using the commercially important IN738LC superalloy as a model material. The same approach could also be applied to adapt other hot-cracking susceptible alloy systems for AM.

36 MATERIALS SCIENCE↗

Harnessing metastability for grain size control in multiprincipal element alloys during additive manufacturing

Abstract Controlling microstructure in fusion-based metal additive manufacturing (AM) remains a significant challenge due to the many parameters that directly impact solidification condition. Multiprincipal element alloys (MPEAs), also known as high entropy alloys, offer a vast compositional space to design for microstructural engineering due to their chemical complexity and exceptional properties. Here, we use the FeMnCoCr system as a model platform for exploring alloy design in MPEAs for AM. By exploiting the decreasing stability of the face-centered cubic phase with increasing Mn content, we achieve notable grain refinement and breakdown of epitaxial columnar grain growth. We employ a multifaceted approach encompassing thermodynamic modeling, operando synchrotron X-ray diffraction, multiscale microstructural characterization, and mechanical testing to gain insight into the solidification physics and its ramifications on the resulting microstructure of FeMnCoCr MPEAs. This work aims toward tailoring desirable grain sizes and morphology through targeted manipulation of phase stability, thereby advancing microstructure control in AM applications.

Wakai, Akane↗

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗