Collaborative Research: Numerical Simulation of the Morphosynthesis of Polycrystalline Biominerals

多晶生物矿物是生物矿物的结晶稳定的晶体多晶型物，其特征在于晶体的堆叠层与相邻层之间的矿物桥。独特的晶体织构产生了特定的材料性质，例如韧性、耐腐蚀性和耐温性，这使得这些晶体对于光学纳米结构（光子带隙、衍射光栅）和特殊涂层（例如，在半导体器件技术中）。因此，目前许多材料科学家都在试图实现这类生物矿物的合成。该项目的目的是提供结晶过程的数学模型和数值模拟的算法工具，以了解该过程的机制，使实验人员能够优化他们的实验室设置，多晶生物矿物的形态合成遵循多阶段结晶过程，包括聚合物诱导的液相结晶，前体（PILP）相、由于成核而发生球晶以及镶嵌介晶薄结构的再结晶。PILP相由生物矿物和与乙醇混合的阴离子聚合物的水溶液组成，并且以液体混合物中富含聚合物的PILP液滴的液-液相分离为特征。混合是照顾由表面声波（SAW）操纵的流体流，其中SAW是由两个锥形叉指换能器在双模式下操作产生的。多晶化开始形成球晶，球晶在基材上扩散，形成均匀的球晶薄膜。连续冷却导致球晶薄膜再结晶成镶嵌多晶薄结构。描述PILP相的液-液相分离可以用不可压Navier-Stokes方程和Cahn-Hilliard方程组成的耦合系统来描述。对于数值模拟，该项目将使用基于时间上隐式离散化的分裂方案和用于空间上离散化的C 0内部惩罚间断伽辽金（C 0-IPDG）方法，相对于计算域的单纯三角剖分。研究了分裂方法的收敛性，并采用面向目标的对偶加权方法实现了分裂方法的时空适应性。作为多晶化的数学模型，该项目研究了相场模型，该模型由测量局部结晶度的动力学方程、生物矿物的浓度场、取向场和冷却过程中温度演变的热方程组成。浓度场方程是一个四阶Cahn-Hilliard型方程。同样，在时间上隐式离散化，在空间上通过C 0-IPDG方法，该项目将使用分裂方法和双加权残差进行时空自适应，其特征在于在最终时间作为目标函数的所需结晶度。模型验证将基于合作实验室提供的实验数据，并将进行系统的参数研究，以调查各种工艺参数的影响。

Polycrystalline biominerals are thermodynamically stable crystal polymorphs of biogenic minerals featuring stacked layers of crystals with mineral bridges between adjacent layers. The unique crystal texture gives rise to specific material properties such as toughness, corrosion resistance, and temperature resistance, which makes these crystals highly attractive for optical nanostructures (photonic band gaps, diffraction gratings) and for special coatings (e.g., in semiconductor device technology). Therefore, many material scientists are currently trying to realize the synthesis of such biominerals. The aim of this project is to provide both a mathematical model for the crystallization process and algorithmic tools for numerical simulations in order to understand the mechanisms of the process, to enable the experimentalists to optimize their laboratory settings, and thus to pave the way for an industrially relevant production line.The morphosynthesis of polycrystalline biominerals follows a multistage crystallization process including a polymer-induced liquid-precursor (PILP) phase, the occurrence of spherulites due to nucleation, and the recrystallization of mosaic mesocrystal thin structures. The PILP phase consists of an aqueous solution of the biomineral and an anionic polymer mixed with ethanol and features a liquid-liquid phase separation in terms of polymer-rich PILP droplets in the liquid mixture. The mixing is taken care of by a surface acoustic waves (SAWs) manipulated fluid flow where the SAWs are generated by two tapered interdigital transducers operating in dual mode. The polycrystallization sets in with the formation of spherulites that spread across the substrate to form a uniform spherulitic thin film. Continuous cooling leads to a recrystallization of the spherulitic thin film into a mosaic polycrystalline thin structure. The liquid-liquid phase separation characterizing the PILP phase can be described by a coupled system consisting of the incompressible Navier-Stokes equations and a Cahn-Hilliard equation. For the numerical simulation, the project will use a splitting scheme based on an implicit discretization in time and C0 Interior Penalty Discontinuous Galerkin (C0-IPDG) methods for discretization in space with respect to simplicial triangulations of the computational domain. The research will study the convergence of the splitting method and realize space-time adaptivity by the goal oriented dual weighted approach. As a mathematical model for the polycrystallization the project investigates a phase field model consisting of the dynamic equations for the measure of local crystallinity, the concentration field for the biomineral, the orientation field, and a heat equation for the evolution of the temperature during the cooling process. The equation for the concentration field is a fourth order Cahn-Hilliard type equation. Again, discretizing implicitly in time and by C0-IPDG methods in space, the project will use a splitting method and dual weighted residuals for space-time adaptivity featuring a desired crystallinity at final time as the objective functional. A model validation will be based on experimental data provided by cooperating laboratories and a systematic parameter study will be performed to investigate the influence of various process parameters.