skip to main content


Title: Improved time integration for phase‐field crystal models of solidification
Abstract

We optimize a numerical time‐stabilization routine for a class of phase‐field crystal (PFC) models of solidification. By numerical experiments, we demonstrate that our simple approach can improve the accuracy of underlying time integration schemes by a few orders of magnitude. We investigate different time integration schemes. Moreover, as a prototypical example for applications, we extend our numerical approach to a PFC model of solidification with an explicit temperature coupling.

 
more » « less
Award ID(s):
2012634
PAR ID:
10482807
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
PAMM
Date Published:
Journal Name:
PAMM
Volume:
23
Issue:
1
ISSN:
2603-5243
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Summary

    While system dynamics are usually derived in continuous time, respective model‐based optimal control problems can only be solved numerically, ie, as discrete‐time approximations. Thus, the performance of control methods depends on the choice of numerical integration scheme. In this paper, we present a first‐order discretization of linear quadratic optimal control problems for mechanical systems that is structure preserving and hence preferable to standard methods. Our approach is based on symplectic integration schemes and thereby inherits structure from the original continuous‐time problem. Starting from a symplectic discretization of the system dynamics, modified discrete‐time Riccati equations are derived, which preserve the Hamiltonian structure of optimal control problems in addition to the mechanical structure of the control system. The method is extended to optimal tracking problems for nonlinear mechanical systems and evaluated in several numerical examples. Compared to standard discretization, it improves the approximation quality by orders of magnitude. This enables low‐bandwidth control and sensing in real‐time autonomous control applications.

     
    more » « less
  2. Neural ordinary differential equations (NODEs) -- parametrizations of differential equations using neural networks -- have shown tremendous promise in learning models of unknown continuous-time dynamical systems from data. However, every forward evaluation of a NODE requires numerical integration of the neural network used to capture the system dynamics, making their training prohibitively expensive. Existing works rely on off-the-shelf adaptive step-size numerical integration schemes, which often require an excessive number of evaluations of the underlying dynamics network to obtain sufficient accuracy for training. By contrast, we accelerate the evaluation and the training of NODEs by proposing a data-driven approach to their numerical integration. The proposed Taylor-Lagrange NODEs (TL-NODEs) use a fixed-order Taylor expansion for numerical integration, while also learning to estimate the expansion's approximation error. As a result, the proposed approach achieves the same accuracy as adaptive step-size schemes while employing only low-order Taylor expansions, thus greatly reducing the computational cost necessary to integrate the NODE. A suite of numerical experiments, including modeling dynamical systems, image classification, and density estimation, demonstrate that TL-NODEs can be trained more than an order of magnitude faster than state-of-the-art approaches, without any loss in performance. 
    more » « less
  3. Abstract

    This work focuses on modeling the interaction between an incompressible, viscous fluid and a poroviscoelastic material. The fluid flow is described using the time‐dependent Stokes equations, and the poroelastic material using the Biot model. The viscoelasticity is incorporated in the equations using a linear Kelvin–Voigt model. We introduce two novel, noniterative, partitioned numerical schemes for the coupled problem. The first method uses the second‐order backward differentiation formula (BDF2) for implicit integration, while treating the interface terms explicitly using a second‐order extrapolation formula. The second method is the Crank–Nicolson and Leap‐Frog (CNLF) method, where the Crank–Nicolson method is used to implicitly advance the solution in time, while the coupling terms are explicitly approximated by the Leap‐Frog integration. We show that the BDF2 method is unconditionally stable and uniformly stable in time, while the CNLF method is stable under a CFL condition. Both schemes are validated using numerical simulations. Second‐order convergence in time is observed for both methods. Simulations over a longer period of time show that the errors in the solution remain bounded. Cases when the structure is poroviscoelastic and poroelastic are included in numerical examples.

     
    more » « less
  4. Abstract

    In this paper, we introduce a new framework for deriving partitioned implicit-exponential integrators for stiff systems of ordinary differential equations and construct several time integrators of this type. The new approach is suited for solving systems of equations where the forcing term is comprised of several additive nonlinear terms. We analyze the stability, convergence, and efficiency of the new integrators and compare their performance with existing schemes for such systems using several numerical examples. We also propose a novel approach to visualizing the linear stability of the partitioned schemes, which provides a more intuitive way to understand and compare the stability properties of various schemes. Our new integrators are A-stable, second-order methods that require only one call to the linear system solver and one exponential-like matrix function evaluation per time step.

     
    more » « less
  5. Abstract

    We introduce a non-isothermal phase-field crystal model including heat flux and thermal expansion of the crystal lattice. The fundamental thermodynamic relation between internal energy and entropy, as well as entropy production, is derived analytically and further verified by numerical benchmark simulations. Furthermore, we examine how the different model parameters control density and temperature evolution during dendritic solidification through extensive parameter studies. Finally, we extend our framework to the modeling of open systems considering external mass and heat fluxes. This work sets the ground for a comprehensive mesoscale model of non-isothermal solidification including thermal expansion within an entropy-producing framework, and provides a benchmark for further meso- to macroscopic modeling of solidification.

     
    more » « less