Skip to feedback

Award Abstract # 1520862
Collaborative Research: Numerical Simulation of the Morphosynthesis of Polycrystalline Biominerals

NSF Org: DMS
Division Of Mathematical Sciences
Recipient: THE UNIVERSITY OF TEXAS AT EL PASO
Initial Amendment Date: July 24, 2015
Latest Amendment Date: August 4, 2017
Award Number: 1520862
Award Instrument: Continuing Grant
Program Manager: Leland Jameson
DMS
 Division Of Mathematical Sciences
MPS
 Directorate for Mathematical and Physical Sciences
Start Date: September 1, 2015
End Date: August 31, 2019 (Estimated)
Total Intended Award Amount: $159,999.00
Total Awarded Amount to Date: $159,999.00
Funds Obligated to Date: FY 2015 = $36,916.00
FY 2016 = $60,791.00

FY 2017 = $62,292.00
History of Investigator:
  • Natasha Sharma (Principal Investigator)
    nssharma@utep.edu
Recipient Sponsored Research Office: University of Texas at El Paso
500 W UNIVERSITY AVE
EL PASO
TX  US  79968-8900
(915)747-5680
Sponsor Congressional District: 16
Primary Place of Performance: University of Texas at El Paso
500 W University Avenue
El Paso
TX  US  79902-5816
Primary Place of Performance
Congressional District:
16
Unique Entity Identifier (UEI): C1DEGMMKC7W7
Parent UEI: C1DEGMMKC7W7
NSF Program(s): COMPUTATIONAL MATHEMATICS
Primary Program Source: 01001516DB NSF RESEARCH & RELATED ACTIVIT
01001617DB NSF RESEARCH & RELATED ACTIVIT

01001718DB NSF RESEARCH & RELATED ACTIVIT
Program Reference Code(s): 7237, 9263
Program Element Code(s): 127100
Award Agency Code: 4900
Fund Agency Code: 4900
Assistance Listing Number(s): 47.049

ABSTRACT

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.

PROJECT OUTCOMES REPORT

Disclaimer

This Project Outcomes Report for the General Public is displayed verbatim as submitted by the Principal Investigator (PI) for this award. Any opinions, findings, and conclusions or recommendations expressed in this Report are those of the PI and do not necessarily reflect the views of the National Science Foundation; NSF has not approved or endorsed its content.

Numerical Simulation of the Morphosynthesis of Polycrystalline Biominerals was focused on the modeling of the morphosynthesis of polycrystalline biominerals described by a coupled Navier-Stokes-Cahn-Hilliard (NS-CH) system and the phase field crystal (PFC) equation. State-of-the-art discretization techniques such as Interior Penalty Discontinuous Galerkin (IPDG) methods to solve the NS-CH system and the PFC equation were proposed. For the Stokes system, the adaptive mesh refinement is driven by residual-type a posteriori error estimators. In particular, the convergence of the adaptive solution process as well as its quasi-optimality in terms of the computational complexity with respect to a properly specified approximation class. In particular, a stronger property namely the contractive for this adaptive IPDG was proved. High quality software implementations will be used to test the algorithms and inform the analysis. High quality software was developed during this project. A numerical scheme for the PFC equation was also developed. This scheme was developed by employing a convex-splitting time discretization and a continuous C0 interior penalty method to spatially discretize the equation,  the numerical scheme is shown to be unconditionally energy stable and uniquely solvable. The investigator and her collaborators are establishing a profound theoretical foundation for adaptive discontinuous Galerkin methods for the numerical simulation of the polycrystalization process.


Last Modified: 01/13/2020
Modified by: Natasha Sharma

Please report errors in award information by writing to: awardsearch@nsf.gov.

Print this page

Back to Top of page