We’re looking forward to start into 2023 with you!

In this CRC colloquium we're happy to welcome:

Sandra Klinge (TU Berlin)

Title: Multiscale material modeling in engineering

The presentation deals with the multiscale modeling of materials with complex microstructures and shows three illustrative examples from different fields. First, the multiscale finite element method (FEM) is used to model cancellous bone and thus to reduce the extent of necessary laboratory tests typical of other methods. The main focus here is placed on the generation of a suitable representative volume element (RVE) consisting of the solid bone phase and the fluid marrow. Macroscopic simulations deal with the ultrasonic tests and enable an analysis of the attenuation dependency on the bone density and excitation frequency.

Second example is related to the simulation of the strain induced crystallization. This model introduces a triple decomposition of the deformation gradient and two specific types of internal variables: regularity of the polymer chain network and thermal flexibility. Starting with the minimum principle of the dissipation potential, it derives evolution equations for internal variables that are able to simulate the formation and the degradation of crystalline regions and to monitor the temperature change during cyclic tensile tests.

The last example uses the multiscale finite element method to simulate the effective material behavior of calcified hydrogels. Within this framework, RVEs are generated to depict the biphasic material microstructure consisting of the organic hydrogel and anorganic calcium phosphate. The choice of the hydrogel matrix influences the type of agglomeration of the anorganic phase and thus dictate the geometry of the RVE. Most commonly, the anorganic phase appears in the form of spherical inclusions or honeycomb grids where the characteristic size of a typical unit might vary. The example also uses asymptotic homogenization technique for the study of the effective diffusivity of hydrogels.

Frederic Legoll (ENPC)

Title: Parareal algorithms for molecular dynamics simulations

In this talk, we consider parareal algorithms in the context of molecular dynamics simulations. The parareal algorithm, which has been originally proposed two decades ago, efficiently solves initial-value problems by using parallel-in-time computations. It is based on a decomposition of the time interval into subintervals and makes use of a predictor-corrector strategy. It first uses a coarse, inexpensive solver to quickly step through the whole time domain, and next refines this approximate solution by using an accurate fine solver which is applied concurrently, in parallel, over each time subinterval (assuming that many processors are available). The prediction-correction process is repeated until convergence is reached.

Although the parareal algorithm, in its original formulation, always converges, it suffers from various limitations in the context of molecular dynamics. In particular, it is observed that the algorithm does not provide any computational gain (in terms of wall-clock time compared to a standard sequential integration) in the limit of increasingly long time-horizons. This numerical observation is backed up with theoretical discussions. We then introduce a modified version of the parareal algorithm where the algorithm adaptively divides the entire time-horizon into smaller time slabs. We numerically show that the adaptive algorithm overcomes the various limitations of the standard parareal algorithm, thereby allowing for significantly improved gains. Several numerical examples (based on the Langevin equation) will be discussed, including the simulation using the LAMMPS software of defects diffusing in a tungsten lattice.

This talk is based on joint works with O. Gorynina, T. Lelievre, D. Perez and U. Sharma.