Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermomechanical effects”

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 451 records · Page 25

Incorporation of metal ions into polyimides

The effects of the incorporation of metal ions into various polyimides on polyimide properties are investigated. Polyimide films derived from 3,3',4,4'-benzophenone tetracarboxylic acid dianhydride (BDTA) 3,3'-diaminobenzophenone (m,m'-DABP), 4,4'-diaminobenzophenone (p,p'-DABP) or 4,4'-oxydianiline were prepared with the concurrent addition of approximately 20 metals in a variety of forms. In general, it is found that the films derived from BDTA + p,p'-DABP were brittle and of poor quality, with brittle films also produced in most of the BDTA + m, m'-DABP polyimides regardless of whether the added metal was hydrate or anhydrous. Thermomechanical analysis, torsional braid analysis, thermal gravimetric analysis, infrared spectral analysis and isothermal studies on many of the polyimide films produced indicate that the softening temperature is generally increased upon the addition of metal ions, at the expense of thermal stability, while no changes in chemical functionality are observed. The best system studied in regard to polymer property enhancement appears to be tri(acetylacetonato)aluminum(III) added to the m, m'-DABP polyamide, which has been found to exhibit four times the lap shear strength of the polyimide alone.

Taylor, L. T.↗

Thermochemical/Thermomechanical Synergies in High Temperature Solid Particle Erosion of CMAS-Exposed EBCs

Environmental barrier coatings (EBCs) are an enabling technology for the use of SiC-based ceramic matrix composites in next generation gas turbine engines. In the extreme engine environment, EBCs must be able to withstand a variety of individual damage mechanisms and their interactions with each other. Ingested particulates/debris can cause both thermochemical and thermomechanical degradation of EBCs. Siliceous debris primarily based on calcium magnesium aluminosilicates (CMAS) can melt and infiltrate and/or react with EBCs above >1200°C. Similarly, ingested debris can lead to mechanical damage and recession of coatings due to particulate erosion. Both modes of degradation can occur simultaneously during engine operation, and it is crucial to comprehensively understand the mechanisms of coating failure due to high-temperature particulate interactions. This study assesses the erosion durability of Yb2Si2O7-based EBCs exposed to CMAS of various loads in NASA Glenn’s Erosion Burner Rig Facility. CMAS exposures and erosion testing were carried out at 1316°C. The effects of CMAS loading and exposure time on EBC erosion durability were evaluated using Al2O3 as an erodent material.

EBC↗

Active interlocking metasurfaces enabled by shape memory alloys

Interlocking metasurfaces (ILMs) are a newly developed joining technology that relies on arrays of interlocking features that transmit force and constrain motion between adjoining bodies in one or more directions. This study explores harnessing the shape memory effect (SME) in Nickel-Titanium shape memory alloys (NiTi SMAs) in structures fabricated using additive manufacturing (AM) to advance the development of active ILMs by creating unit cells that open or close at specific temperatures. The study encompasses designing and fabricating two distinct interlocking array configurations using near-equiatomic NiTi powder and the laser powder bed fusion (L-PBF) AM technique, following a previously developed AM process optimization framework to manufacture defect-free parts. To guide the design process, finite element analysis (FEA) was employed to predict strain values during engage-disengage cycles. The martensitic transformation characteristics of the ILMs were characterized. Thermomechanical testing revealed that the ILMs demonstrate high locking force once engaged, coupled with complete shape recovery and good cyclic stability. Digital image correlation (DIC) was also employed to validate the FEA predictions during the engage-disengage cycles. The results indicate that NiTi SMA-based ILMs can be designed and fabricated into complex shapes using L-PBF. By leveraging the SME, the functionality of an ILM can be improved upon. The combination of computational modeling, additive manufacturing, and thermomechanical and physical property characterization provides a framework for designing future ILMs out of active materials.

Additive manufacturing↗

Spatiotemporal Learning in Power Modules: Wavelet-Enhanced Forecasting of Thermomechanical Degradation

Detecting internal defects in power electronics packages is critical for their performance and reliability, especially under extreme operating conditions, as these defects can lead to catastrophic failure if not properly addressed. Confocal scanning acoustic microscopy (C-SAM) plays a key role in the nondestructive evaluation of bond layer degradation within a power electronics package by detecting defects such as delamination, voids, and cracks. However, accurately quantifying and predicting these defects from C-SAM images remains a significant challenge due to the low noise-to-signal ratio, which typically arises from both imaging process and bond patterns itself. In this paper, we explore machine learning strategies for processing C-SAM images and providing predictive models of defect growth. We use C-SAM images of sintered copper and sintered silver samples, which are obtained under accelerated thermal experiments, as the representative dataset for our study. We investigate the effect of Fourier transforms and wavelet transforms on these datasets to remove high-frequency noise and address noise across multiple scales with histogram equalization to enhance the contrast and improve the visibility of defects. As a result, defect boundaries can be clearly distinguished, enabling more accurate tracking of their growth over time. We then employ different time-series forecasting algorithms on the denoised images to formulate an image-based lifetime prediction model. Statistical models and deep-learning techniques are trained on images obtained in the early stages of thermal shock, and defect growth in the later stages is predicted. Our work serves as a preliminary attempt to improve the accuracy of lifetime prediction models of power electronics packages, which is critical under extreme operating environments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Quasi-elastic neutron scattering study on dynamically asymmetric polymer blends

Compatibility between polymers with different glass transition temperatures controls the thermomechanical properties of blends. This work explores the segmental dynamics of poly(methyl acrylate) (PMA) chains when they are blended with polymers of different rigidities and miscibility. Inspired by the intriguing dynamic asymmetry of chains within the interfacial layer of nanoparticles, this study aims to understand dynamic heterogeneity in dynamically asymmetric blends of PMA/poly(methyl methacrylate) (PMMA), PMA/polystyrene (PS), and PMA/poly(ethylene oxide) (PEO) using the differential scanning calorimetry (DSC) and quasi-elastic neutron scattering (QENS) measurements below and above the glass transition temperature of PMMA or PS. Further, results revealed that the segmental jump distance of PMA increased on blending due to volume enhancement. The reduced effective diffusivity of PMA in PMMA is attributed to PMMA's enhanced flexibility (lower characteristic ratio, C ∞ ) and superior interaction (lower ) when compared with the PS environment. The results demonstrate that the chain rigidity and miscibility affect the free volume, interchain cooperativity, and segmental dynamics in dynamically asymmetric blends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reduced‐Order Modeling of Energetic Materials Using Physics‐Aware Recurrent Convolutional Neural Networks in a Latent Space (LatentPARC)

Physics-aware deep learning (PADL) has gained popularity for use in spatiotemporal dynamics simulations, such as those in computational modeling of energetic materials (EM). We show that the challenge PADL methods face while learning complex field evolution problems can be simplified and accelerated by decoupling it into two tasks: learning complex geometric features in evolving fields and modeling dynamics over these features in a lower-dimensional feature space. We build upon our previous work on physics-aware recurrent convolutional neural networks (PARC). PARC embeds knowledge of underlying physics into its neural network architecture for more robust and accurate prediction of evolving physical fields. PARC was shown to effectively learn complex nonlinear features such as the formation of hotspots and coupled shock fronts in various initiation scenarios of EMs, as a function of microstructures, serving effectively as a microstructure-aware burn model. Here, we further accelerate PARC and reduce its computational cost by projecting the original dynamics onto a lower-dimensional invariant manifold, or “latent space.” The projected latent representation encodes the complex geometry of evolving fields (e.g., temperature and pressure) in a set of data-driven features. The reduced dimension of this latent space allows us to learn the dynamics during the initiation of EM with a lighter and more efficient model. We observe a significant decrease in training and inference time while maintaining results comparable to PARC at inference. This work takes steps towards enabling rapid prediction of EM thermomechanics at larger scales and characterization of EM structure–property–performance linkages at a full application scale.

Mathematics and Computing↗

Thermochemical/Thermomechanical Synergies in High-Temperature Solid Particle Erosion of CMAS-Exposed EBCs

Environmental barrier coatings (EBCs) are an enabling technology for the use of SiC-based ceramic matrix composites in next generation gas turbine engines. In the extreme engine environment, EBCs must be able to withstand a variety of individual damage mechanisms and their interactions with each other. Ingested particulates/debris can cause both thermochemical and thermomechanical degradation of EBCs. Siliceous debris primarily based on calcium magnesium aluminosilicates (CMAS) can melt and infiltrate and/or react with EBCs above 1200°C. Similarly, ingested debris can lead to mechanical damage and recession of coatings due to particulate erosion. Both modes of degradation can occur simultaneously during engine operation, and it is crucial to comprehensively understand the mechanisms of coating failure due to high-temperature particulate interactions. This study assesses the erosion durability of Yb 2 Si 2 O 7 -based EBCs exposed to CMAS of various loads in NASA Glenn’s Erosion Burner Rig Facility. CMAS exposures and erosion testing were carried out at 1316°C. The effects of CMAS loading and exposure time on EBC erosion durability were evaluated using Al 2 O 3 as an erodent material.

CMAS↗

Thermochemical/Thermomechanical Synergies in High Temperature Solid Particle Erosion of CMAS Exposed EBCs

Environmental barrier coatings (EBCs) are an enabling technology for the use of SiC-based ceramic matrix composites in next generation gas turbine engines. In the extreme engine environment, EBCs must be able to withstand a variety of individual damage mechanisms and their interactions with each other. Ingested particulates/debris can cause both thermochemical and thermomechanical degradation of EBCs. Siliceous debris primarily based on calcium magnesium aluminosilicates (CMAS) can melt and infiltrate and/or react with EBCs above >1200°C. Similarly, ingested debris can lead to mechanical damage and recession of coatings due to particulate erosion. Both modes of degradation can occur simultaneously during engine operation, and it is crucial to comprehensively understand the mechanisms of coating failure due to high-temperature particulate interactions. This study assesses the erosion durability of Yb 2 Si 2 O 7 -based EBCs exposed to CMAS of various loads in NASA Glenn’s Erosion Burner Rig Facility. CMAS exposures and erosion testing were carried out at 1316°C. The effects of CMAS loading and exposure time on EBC erosion durability were evaluated using Al2O3 as an erodent material.

EBC↗

International Research Project on the Effects of Chemical Ageing of Polymers on Performance Properties: Chemical and Thermal Analysis

Work during the past six months has included significant research in several areas aimed at further clarification of the aging and chemical failure mechanism of thermoplastics (PVDF or Tefzel) pipes. Among the areas investigated were the crystallinity changes associated with both the Coflon and Tefzel after various simulated environmental exposures using X-ray diffraction analysis. We have found that significant changes in polymer crystallinity levels occur as a function of the exposures. These crystallinity changes may have important consequences on the fracture, fatigue, tensile, and chemical resistance of the materials. We have also noted small changes in the molecular weight distribution. Again these changes may result in variations in the mechanical and chemical properties in the material. We conducted numerous analytical studies with methods including X-ray Diffraction, Gel Permeation Chromatography, Fourier Transform Infrared Spectroscopy, Ultra- Violet Scanning Analysis, GC/Mass Spectrometry, Differential Scanning Calorimetry and Thermomechanical Analysis. In the ultra-violet analysis we noted the presence of an absorption band indicative of triene formation. We investigated a number of aged samples of both Tefzel and Coflon that were forwarded from MERL. We also cast films at SWT and subjected these films to a refluxing methanol 1% ethylene diamine solution. An updated literature search was conducted using Dialog and DROLLS to identify any new papers that may have been published in the open literature since the start of this project. The updated literature search and abstracts are contained in the Appendix section of this report.

Bulluck, J. W.↗

Quantitative measurements of dislocations in metals for advancing predictive simulations

LLNL applications require scientists to predict how materials evolve under various thermomechanical conditions. While this is achieved through physics-based simulations, uncertainty in the predictions of mechanical properties remains a serious challenge that limits the predictive capabilities of models because we lack methods to compare predictions of atomic-scale defects (dislocations) with experimental measurements. High energy X-ray diffraction (HEXRD) is the most relevant technique that can provide the necessary statistical information on dislocations. However, this technique is not yet quantitative because we lack a precise understanding of the relationship between X-ray diffraction patterns and the underlying material dislocation content and arrangements. To address this need, we used our novel computational X-ray diffraction method to simulate the effect of dislocations on the diffraction patterns. We compared virtual and experimental diffraction patterns. Results allowed us to clearly establish the relationship between X-ray diffraction patterns and the underlying dislocation structures, proving that it is feasible to quantitatively measure dislocation statistics with HEXRD. This project delivered a method that can provide the missing piece to LLNL’s mechanical property simulations in advanced metals by obtaining experimentally long-needed quantitative dislocation data, which could fully enable predictive capabilities.

36 MATERIALS SCIENCE↗

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗

Software For Advanced Large-scale Analysis Of Magnetic Confinement For Numerical Design, Engineering & Research (salamander)

As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems

Simon, Pierre-Clement [Idaho National Laboratory (↗

Review on Preprocessing Strategies, Deactivation, Thermal Safety, and Future Perspectives in Lithium-Ion Battery Recycling

The rapid growth in the use of lithium-ion batteries (LIBs) in electric vehicles, consumer electronics, and renewable energy storage has made effective end-of-life management essential. Recycling LIBs is critical not only for resource recovery and environmental protection but also for ensuring safety and economic viability. This review focuses on the preprocessing technologies that precede typical recovery processes, including disassembly, sorting, discharging, electrolyte removal, dismantling, thermal treatment, separation, and flotation. These steps play a foundational role in determining the efficiency, safety, and environmental impact of LIB recycling. LIBs pose substantial fire and explosion risks due to residual charge, flammable electrolytes, and reactive materials. The conditions and successive progression of the exothermic reactions which lead to thermal runaway has been discussed. It also explores secure deactivation techniques such as external circuit discharge, saline immersion, and thermomechanical methods, alongside fire prevention strategies including the use of flame retardants, elimination of oxidants, and reduction of heat generation and accumulation. Challenges and future directions are outlined, highlighting the need for standardized designs, automation, and safer, more sustainable recycling infrastructure. Furthermore, this review is distinguished by its focused analysis of preprocessing and deactivation steps, with particular attention to the thermal safety engineering aspects of LIB recycling.

Battery deactivation↗

Prototype Modeling for a Light-Trapping Planar-Cavity Enclosed Particle Solar Receiver

Concentrating solar thermal (CST) systems present a promising avenue for affordable and reliable energy production. Solar receivers are key components that determine the efficiency and longevity of these systems. Particle-based solar receivers have emerged as a compelling alternative to traditional technologies, offering several advantages that address limitations in current CST systems. This is especially true as next-generation CST technologies target applications including electricity generation, thermochemical processes, and industrial process heat, many of which necessitate higher operating temperatures than current commercial molten salt systems. Molten-salt thermal energy storage (TES) systems, commonly used in CSP, face challenges related to freezing and corrosion. Particle-based TES systems, in contrast, do not experience these issues, as particles are stable at high temperatures, exceeding 1000 degrees Celsius. This capability allows for a wider range of applications, including those requiring higher temperatures for industrial processes and efficient electricity generation. A novel innovation in particle-based solar receiver technology is the light-trapping planar cavity receiver (LTPCR) configuration developed by NREL. The LTPCR design consists of small cavity-like structures using opaque planar surfaces, enabling efficient capture and absorption of solar energy. A high incident flux concentration at the cavity aperture is absorbed on the receiver walls, and subsequently transferred to particles on the inside of cavities. The particles flow through the system, forming a fluidized bed inside of the receiver panels, effectively capturing the absorbed solar heat. Air is used as a fluidizing medium in this process to enhance particle heat transfer and mixing. The effectiveness of this design lies in its ability to manage solar flux conditions and ensure high solar-to-thermal receiver efficiency. A 100-kW prototype is currently being tested at the King Saud University in Saudi Arabia to assess the receiver performance. A range of modeling analyses for the optical, thermal, and mechanical effects were conducted to assess the performance of the receiver under on-sun conditions. The solar flux resulting from the KSU heliostat field was modeled using NREL SolTrace software and produced up to 600 kW/m2 at the receiver aperture. The solar flux absorbed on the receiver walls was then used within a computational fluid dynamics (CFD) model to predict wall temperature distributions along with radiation and convection loss. A two-phase CFD model was developed for the fluidized bed of silica sand inside the receiver panels to predict local wall-to-particle heat transfer coefficients, particle temperature distributions, and outlet temperature of the particles. We have also conducted analyses to understand the thermomechanical behavior of these innovative enclosed light-trapping solar receivers optimized for particle heating. We used finite element analysis (FEA) to predict the receiver's performance using temperature distributions obtained from CFD and based on the resulting stress profiles, evaluated creep-fatigue damage with a goal of achieving a 30-year service life. Analysis showed a significant impact of the particle-to-wall heat transfer coefficients (HTCs) on receiver performance, with higher HTCs resulting in reduced stress and increased lifespan. For instance, when using Inconel 740H, increasing the HTC from 800 W/m2 K to 1400 W/m2 K increased the creep life from 4,000 hours to over 100,000 hours. This highlights the importance of understanding and optimizing heat transfer in the design of high-efficiency receivers.

14 SOLAR ENERGY↗

Nonlinear structural analysis for fiber-reinforced superalloy turbine blades

A computational capability for predicting the nonlinear thermomechanical structural response of fiber-reinforced superalloy (FRS) turbine blades is described. This capability is embedded in a special purpose computer code (COBSTRAN) developed at the NASA Lewis Research Center. Special features of this computational capability include accounting for: fiber/matrix reaction, nonlinear and anisotropic material behavior, complex stress distribution due to local and global heterogeneity, and residual stresses due to initial fabrication and/or inelastic behavior during subsequent missions. Numerical results are presented from analyses of a hypothetical FRS turbine blade subjected to a fabrication process and subsequent mission cycle. The results demonstrate the capabilities of this computational tool to; predict local stress/strain response and capture trends of local nonlinear and anisotropic material behavior, relate the effects of this local behavior to the global response of a multilayered fiber-composite turbine blade, and trace material history from fabrication through successive missions.

Hopkins, D. A.↗

Stress-rupture behavior of small diameter polycrystalline alumina fibers

Continuous length polycrystalline alumina fibers are candidates as reinforcement in high temperature composite materials. Interest therefore exists in characterizing the thermomechanical behavior of these materials, obtaining possible insights into underlying mechanisms, and understanding fiber performance under long term use. Results are reported on the time-temperature dependent strength behavior of Nextel 610 and Fiber FP alumina fibers with grain sizes of 100 and 300 nm, respectively. Below 1000 C and 100 hours, Nextel 610 with the smaller grain size had a greater fast fracture and rupture strength than Fiber FP. The time exponents for stress-rupture of these fibers were found to decrease from approximately 13 at 900 C to below 3 near 1050 C, suggesting a transition from slow crack growth to creep rupture as the controlling fracture mechanism. For both fiber types, an effective activation energy of 690 kJ/mol was measured for rupture. This allowed stress-rupture predictions to be made for extended times at use temperatures below 1000 C.

Yun, Hee Mann↗

Topology-Informed Design Rules for Deconstructable Thermoset Copolymer Networks

Existing models of thermoset deconstruction facilitated by incorporating cleavable comonomers rely on a mean-field reverse gel point paradigm, which predicts network dissolution once cleavable bonds reach a critical stoichiometric threshold, but does not account for where those bonds reside within the network architecture. Using reactive coarse-grained molecular dynamics simulations coupled with graph-theoretic analysis, we extend this stoichiometric picture to show that deconstructability is governed by the curing-imprinted network topology rather than stoichiometry alone. This topological organization is hierarchical: at the local scale, the elastic effectiveness of cross-link junctions determines which cross-links constitute the load-bearing scaffold; at the mesoscale, the cross-linking rate kinetically templates that scaffold into topologically modular communities─densely cross-linked clusters connected by sparse bridging strands that sustain network connectivity. Using betweenness centrality to identify nodes that disproportionately lie on intercommunity shortest paths, we demonstrate that effective deconstruction of the network into macromolecular fragments requires cleavable comonomers to intercept these high-centrality bridging strands. We further find that under uniform, disassortative comonomer incorporation, this topological requirement provides a mechanistic basis for extending the reverse gel point to incorporate network topology. We also show that modularity imposes a fundamental limit on fragment uniformity that persists even when the centrality requirement is met. Finally, we demonstrate that chain stiffness provides a nearly independent lever to suppress mechanically redundant cross-links and raise the glass transition temperature without significantly altering the deconstruction outcome. Together, these findings reframe the thermoset design space around network topology and provide actionable guidelines for engineering thermoset copolymers with predictable deconstructability and targeted thermomechanical performance.

coarse-grained molecular dynamics↗

An implementation of the phase-field model based on coupled thermomechanical finite element solvers for large-strain twinning, explicit dynamic fracture and the classical Stefan problem

The implementation of a phase-field model in finite elements usually requires significant expertise and involves the development of a user element with additional degrees of freedom. An alternative implementation of the phase-field model within a thermo-mechanical finite element simulation package was presented in (Cho et al 2012 Int. J. Solids Struct. 49 1973–1992), where the phase-field variable is treated as the temperature degree of freedom. However, this approach has only been used for small strain phase-field modelling of martensitic transformations and quasistatic phase-field modelling of fracture. Here, we present a phase-field finite element implementation via the temperature degree of freedom for several additional cases from the literature: (i) the large-strain phase-field description of deformation twinning presented in (Clayton and Knap 2011 Physica D 240 841–858), (ii) phase-field description of brittle fracture with inertial effects based on the theory from (Molnár and Gravouil 2017, Finite Elem. Anal. Des. 130 27–38) and (Miehe et al 2010 Int. J. Numer. Methods Eng. 83 1273–1311) and (iii) the classical Stefan problem of solidification presented in (Mackenzie and Robertson 2002 J. Comput. Phys. 181 526–544). The last problem involves the temperature and phase-field variables as unknowns.

36 MATERIALS SCIENCE↗