Phase-field

Phase-field modelling has become a versatile method for predicting microstructural evolution and phase transformations in complex materials systems. The Pace3D simulation framework provides a high-performance, massively parallel implementation of a generalized multi-phase-field model, enabling efficient large-scale simulations on modern HPC architectures. 

Pace3D solves the evolution equations of an N-order-parameter phase-field formulation, allowing the treatment of multi-phase material systems as well as polycrystalline grain and particle structures. The order parameters may represent phases (solid, liquid, gas) or individual crystals, providing a flexible framework for a wide range of applications.

To achieve fast time-to-solution, Pace3D combines advanced model reductions, numerical optimizations, and MPI-based parallelization. The framework has been continuously developed for more than 25 years and is supported by numerous benchmark applications, validating interfacial dynamics as well as triple- and multi-junction behavior.

With its proven scalability, physical fidelity, and long-term development history, Pace3D offers a reliable platform for predictive phase-field simulations in industrial and scientific environments.

 

3D phase-field simulation results showing dendritic and liquid morphologies during growth at rotational misorientations θ_R of 45° (A–C) and 15° (D–F). The grand-chemical potential model demonstrates that secondary dendrite arms develop and persist at higher misorientation angles.

Microsegregation of alloying elements: (A) morphological evolution, and (B) concentration profiles along lines 1 and 2 (marked in A) for carbon, chromium, and nickel. Distinct microsegregation patterns emerge between open liquid channels and trapped liquid pockets. Learn more...

Schematic of the weld pool, mushy zone, and solidified weld seam during laser beam welding (LBW). The mushy zone illustrates δ-Fe dendrite orientations exhibiting cubic crystal anisotropy inherited from partially melted grains in the base material, with modeling of complex phase evolution and solute microsegregation within the dendritic microstructure.

Comparison of experimental and phase-field simulated solute segregation in a quaternary Fe-Cr-Ni-C austenitic stainless steel. Concentration distributions of alloying elements—nickel (A,B) and chromium (C,D)—across dendritic and interdendritic regions, with local distribution fields for quantitative comparison.


Diffusion & Grand Chem

Grand chemical potential-driven morphological patterns—such as those in solidification, grain coarsening, and solid-state transformations—along with diffusion processes distinguishing surface, grain boundary, and bulk pathways, are primarily investigated using multi-component, multi-phase-field simulations for industrially relevant alloy systems, including steels, Al-, and Ni-based alloys. CALPHAD-based Gibbs energy functions provide realistic microsegregation predictions for critical elements, differentiating interstitial diffusion (e.g., carbon, boron) from substitutional diffusion (e.g., manganese, nickel). These capabilities advance alloy design, defect identification, and industrial manufacturing optimization.

Heat transfer

Pace3D includes a solver module for heat transfer that can be coupled with phase-field evolution as well as with other solver modules such as solid mechanics and computational fluid dynamics. This enables the simulation of strongly coupled thermo-physical processes. The module accounts for latent heat release during phase transformations and allows the consideration of thermo-mechanical effects through coupling with solid mechanics. In addition, convective heat transport in fluid flow can be modeled through its integration with the CFD solver, enabling consistent treatment of heat transfer in multiphysics simulations involving solids, fluids, and evolving microstructures.

Heat distribution in the air–aluminum domain for the experimental (left) and the synthetic (right) open-cell foam model. Learn more...

Heat distribution including temperature isolines in air–aluminum domains of different porosity after 5 s of physical time. Learn more...

Example application of the developed method for calculating fluid flow and heat transfer in cellular solids: a quarter of an encapsulated foam structure exposed to a fluid flow and heated at the cylindrical shell; (left) temperature distribution in a cross-section and (right) temperature distribution on the surface of the foam structure as well as velocity vectors in a cross-section. Learn more...


Interaction of rigid particles with a surrounding fluid: The left image shows particles fluidizing under gravity in a fluidized bed, together with the corresponding velocity field (Reynolds number based on the velocity magnitude). The right image depicts rigid particles falling under gravity into a three-dimensional foam structure. The model captures the interaction between fluid flow and rigid body motion, allowing both freely moving particles driven by the flow and flow fields induced by forced particle motion. Learn more...

Simulation of the rise of a gas bubble in a liquid. The rise velocity over time and the bubble shape at a specific instant are compared for different viscous stress models. The simulations account for capillary effects as well as density and viscosity contrasts between the two phases. Fluid flow is modeled using a two-phase formulation coupled with phase-field evolution to capture the dynamics of the fluid interface. Learn more...

Impact of a rigid sphere onto the interface between two immiscible fluids. Two-dimensional simulations compare different surface wettabilities, represented by varying equilibrium contact angles, and their influence on interface deformation and fluid–solid interaction during impact. The simulations are performed with Pace3D using a two-phase flow formulation coupled with phase-field evolution and rigid body motion, including capillary effects and contact-angle modeling according to Young’s law. Learn more...


Fluid dynamics

Pace3D includes a computational fluid dynamics (CFD) solver module that provides both a Lattice–Boltzmann solver and a Navier–Stokes solver. These approaches are coupled with the phase-field method, enabling simulations of a wide range of multiphase flow problems. Fluid and solid regions can be treated within a single computational domain using a diffuse-domain formulation. The framework supports macroscopically immiscible multiphase flow through a continuum surface force model capturing capillary effects. It further includes a solver for rigid particulate flows based on a phase-field fictitious-domain approach. In addition, porous-media flow can be modeled using a Brinkman–Forchheimer formulation, and thermal effects are incorporated via temperature-dependent viscosities and the Boussinesq approximation.

Mechanics

Continuum-mechanical modeling is essential for predicting stress evolution, deformation, and failure in complex materials. The Pace3D simulation framework provides a high-performance, massively parallel mechanics solver embedded in a generalized multiphase-field environment. It solves the balance of linear momentum consistently coupled to evolving microstructures, enabling simulations of elastic and inelastic deformation in multiphase and polycrystalline systems. The framework supports small- and finite-strain formulations as well as thermo-mechanical and chemo-mechanical couplings, capturing diffusion-induced stresses, eigenstrains, and microcrack formation within a unified variational setting. Advanced numerical schemes, consistent interface treatments, and MPI parallelization ensure efficient simulations. The solver has been validated against analytical benchmarks and chemo-mechanical test cases, providing a scalable platform for predictive multiphysics simulations.

Simulation of the compression of a two-phase foam structure mimicking an aerogel undergoing large deformations. The model describes elastic behavior within an Eulerian framework, enabling robust simulations of porous materials subject to significant deformation. Learn more...

Schematic illustration of the simulation setup used to model the elastic response of foam structures. By combining structure generation algorithms with solid-mechanics solvers for porous media, the influence of foam morphology on the elastic behaviour of porous materials can be studied systematically. Learn more...

Simulation of domain configurations in ferroelectric thin films under varying misfit strain. The model incorporates thin-film boundary conditions and couples phase-field evolution with polarization dynamics, solid mechanics, and space-charge transport. This multiphysics approach enables the study of strain-dependent domain structure formation in ferroelectric materials. Learn more...


Concentration in uncracked (left) and cracked (right) secondary particle during PITT simulation.

Evolution of Li concentration and phases during 1C CC-CV charge: Top row shows a dense agglomerate where concentration gradients form in the radial direction and the phase transitions follow a shrinking core behaviour. The nanoporous agglomerate in the bottom row exhibits a more homogenous lithium distribution and follows a mosaic pattern during the phase transition. Patches of each phase have been labelled for reference. Learn more...

Secondary particle morphology of NMC811 with colors marking primary particles (a) and (b), crack order-parameter φc (b) and Li-ion concentration 𝑐 after electrolyte infiltration (c). The black vectors in (a) indicate the c-lattice axis of the primary particles.

Secondary particle morphology of NMC with colors marking primary particles (a), crack order-parameter φc (b)
and Li-ion concentration c after electrolyte infiltration (c)


Electrochemistry

Fast charging and long-term durability of battery electrodes are strongly governed by their microstructure. During electrochemical cycling, ion intercalation drives phase transformations, concentration gradients, and the evolution of mechanical stresses that can ultimately lead to microcrack formation.

Our multi-physics battery framework enables modelling of these processes in artificial as well as real 3D microstructures. It captures anisotropy of ion diffusion and chemo-mechanics, as well as the resulting microcrack formation. This provides detailed insight into how features such as primary and secondary particle size, texture and interfaces influence performance and degradation.