-
Isotropic finite-difference approximations for phase-field simulations of polycrystalline alloy solidification
Authors:
Kaihua Ji,
Amirhossein Molavi Tabrizi,
Alain Karma
Abstract:
Phase-field models of microstructural pattern formation during alloy solidification are commonly solved numerically using the finite-difference method, which is ideally suited to carry out computationally efficient simulations on massively parallel computer architectures such as Graphic Processing Units. However, one known drawback of this method is that the discretization of differential terms in…
▽ More
Phase-field models of microstructural pattern formation during alloy solidification are commonly solved numerically using the finite-difference method, which is ideally suited to carry out computationally efficient simulations on massively parallel computer architectures such as Graphic Processing Units. However, one known drawback of this method is that the discretization of differential terms involving spatial derivatives introduces a spurious lattice anisotropy that can influence the solid-liquid interface dynamics. We find that this influence is significant for the case of polycrystalline dendritic solidification, where the crystal axes of different grains do not generally coincide with the reference axes of the finite-difference lattice. In particular, we find that with the commonly used finite-difference implementation of the quantitative phase-field model of binary alloy solidification, both the operating state of the dendrite tip and the dendrite growth orientation are strongly affected by the lattice anisotropy. To circumvent this problem, we use known methods in both real and Fourier space to derive finite-difference approximations of leading differential terms in 2D and 3D that are isotropic at order $h^2$ of the lattice spacing $h$. Importantly, those terms include the divergence of the anti-trapping current that is found to have a critical influence on pattern selection. The 2D and 3D discretizations use an approximated form of the anti-trapping current that facilitates the Fourier-space derivation of the associated isotropic differential operator at $O(h^2)$, but we also derive a 2D discretization of the standard form of this current. Finally, we present 2D and 3D phase-field simulations of alloy solidification, showing that the isotropic finite-difference implementations can dramatically reduce spurious lattice anisotropy effects.
△ Less
Submitted 24 October, 2021;
originally announced October 2021.
-
Spatiotemporal Organization of Electromechanical Phase Singularities during High-Frequency Cardiac Arrhythmias
Authors:
A. Molavi Tabrizi,
A. Mesgarnejad,
M. Bazzi,
S. Luther,
J. Christoph,
A. Karma
Abstract:
Ventricular fibrillation (VF) is a lifethreatening electromechanical dysfunction of the heart associated with complex spatiotemporal dynamics of electrical excitation and mechanical contraction of the heart muscle. It has been hypothesized that VF is driven by three-dimensional (3D) rotating electrical scroll waves, which can be characterized by filament-like electrical phase singularities (EPS).…
▽ More
Ventricular fibrillation (VF) is a lifethreatening electromechanical dysfunction of the heart associated with complex spatiotemporal dynamics of electrical excitation and mechanical contraction of the heart muscle. It has been hypothesized that VF is driven by three-dimensional (3D) rotating electrical scroll waves, which can be characterized by filament-like electrical phase singularities (EPS). Recently, it was shown that rotating excitation waves during VF are associated with rotating waves of mechanical deformation. 3D mechanical scroll waves and mechanical filaments describing their rotational core were observed in the ventricles by using high-resolution ultrasound. The findings suggest that the spatiotemporal organization of cardiac fibrillation may be assessed from waves of mechanical deformation. However, the complex relationship between excitation and mechanical waves during VF is currently not understood. Here, we study the fundamental nature of mechanical phase singularities (MPS), their spatiotemporal organization and relation with EPS. We demonstrate the existence of two fundamental types of MPS: "paired singularities", which are co-localized with EPS, and "unpaired singularities", which can form independently. We show that the unpaired singularities emerge due to the anisotropy of the active force field, generated by fiber anisotropy in cardiac tissue, and the non-locality of elastic interactions, which jointly induce strong spatiotemporal inhomogeneities in the strain fields. The inhomogeneities lead to the breakup of deformation waves and create MPS, even in the absence of EPS, which are typically associated with excitation wave break. We exploit these insights to develop an approach to discriminate paired and unpaired MPS. Our findings provide a fundamental understanding of the complex spatiotemporal organization of electromechanical waves in the heart.
△ Less
Submitted 14 April, 2022; v1 submitted 19 July, 2021;
originally announced July 2021.
-
Predicting Solvation Free Energies and Thermodynamics in Polar Solvents and Mixtures Using a Solvation-Layer Interface Condition
Authors:
Amirhossein Molavi Tabrizi,
Spencer Goossens,
Ali Mehdizadeh Rahimi,
Matthew G. Knepley,
Jaydeep P. Bardhan
Abstract:
We demonstrate that with two small modifications, the popular dielectric continuum model is capable of predicting, with high accuracy, ion solvation thermodynamics in numerous polar solvents, and ion solvation free energies in water--co-solvent mixtures. The first modification involves perturbing the macroscopic dielectric-flux interface condition at the solute--solvent interface with a nonlinear…
▽ More
We demonstrate that with two small modifications, the popular dielectric continuum model is capable of predicting, with high accuracy, ion solvation thermodynamics in numerous polar solvents, and ion solvation free energies in water--co-solvent mixtures. The first modification involves perturbing the macroscopic dielectric-flux interface condition at the solute--solvent interface with a nonlinear function of the local electric field, giving what we have called a solvation-layer interface condition (SLIC). The second modification is a simple treatment of the microscopic interface potential (static potential). We show that the resulting model exhibits high accuracy without the need for fitting solute atom radii in a state-dependent fashion. Compared to experimental results in nine water--co-solvent mixtures, SLIC predicts transfer free energies to within 2.5 kJ/mol. The co-solvents include both protic and aprotic species, as well as biologically relevant denaturants such as urea and dimethylformamide. Furthermore, our results indicate that the interface potential is essential to reproduce entropies and heat capacities. The present work, together with previous studies of SLIC illustrating its accuracy for biomolecules in water, indicates it as a promising dielectric continuum model for accurate predictions of molecular solvation in a wide range of conditions.
△ Less
Submitted 14 November, 2016; v1 submitted 3 November, 2016;
originally announced November 2016.
-
Generalizing The Mean Spherical Approximation as a Multiscale, Nonlinear Boundary Condition at the Solute--Solvent Interface
Authors:
Amirhossein Molavi Tabrizi,
Matthew G. Knepley,
Jaydeep P. Bardhan
Abstract:
In this paper we extend the familiar continuum electrostatic model with a perturbation to the usual macroscopic boundary condition. The perturbation is based on the mean spherical approximation (MSA), to derive a multiscale hydration-shell boundary condition (HSBC). We show that the HSBC/MSA model reproduces MSA predictions for Born ions in a variety of polar solvents, including both protic and ap…
▽ More
In this paper we extend the familiar continuum electrostatic model with a perturbation to the usual macroscopic boundary condition. The perturbation is based on the mean spherical approximation (MSA), to derive a multiscale hydration-shell boundary condition (HSBC). We show that the HSBC/MSA model reproduces MSA predictions for Born ions in a variety of polar solvents, including both protic and aprotic solvents. Importantly, the HSBC/MSA model predicts not only solvation free energies accurately but also solvation entropies, which standard continuum electrostatic models fail to predict. The HSBC/MSA model depends only on the normal electric field at the dielectric boundary, similar to our recent development of an HSBC model for charge-sign hydration asymmetry, and the reformulation of the MSA as a boundary condition enables its straightforward application to complex molecules such as proteins.
△ Less
Submitted 14 July, 2016;
originally announced July 2016.