An official website of the United States government
Abstract We propose two mass and heat energy conservative, unconditionally stable, decoupled numerical algorithms for solving the Cahn–Hilliard–Navier–Stokes–Darcy–Boussinesq system that models thermal convection of two‐phase flows in superposed free flow and porous media. The schemes totally decouple the computation of the Cahn–Hilliard equation, the Darcy equations, the heat equation, the Navier–Stokes equations at each time step, and thus significantly reducing the computational cost. We rigorously show that the schemes are conservative and energy‐law preserving. Numerical results are presented to demonstrate the accuracy and stability of the algorithms.
more »
« less
A computational study of lateral phase separation in biological membranes
Abstract Conservative and non‐conservative phase‐field models are considered for the numerical simulation of lateral phase separation and coarsening in biological membranes. An unfitted finite element method is proposed to allow for a flexible treatment of complex shapes in the absence of an explicit surface parametrization. For a set of biologically relevant shapes and parameter values, the paper compares the dynamic coarsening produced by conservative and non‐conservative numerical models, its dependence on certain geometric characteristics and convergence to the final equilibrium.
more »
« less
- PAR ID:
- 10453972
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical Methods in Biomedical Engineering
- Volume:
- 35
- Issue:
- 3
- ISSN:
- 2040-7939
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract In a recent article, one of the authors developed a multigrid technique for coarse‐graining dynamic powergrid models. A key component in this technique is a relaxation‐based coarsening of the graph Laplacian given by the powergrid network and its weighted graph, which is represented by the admittance matrix. In this article, we use this coarsening strategy to develop a multigrid method for solving a static system of nonlinear equations that arises through Ohm's law, the so‐called powerflow equations. These static equations are tightly knitted to the dynamic model in that the full powergrid model is an algebraic‐differential system with the powerflow equations describing the algebraic constraints. We assume that the dynamic model corresponds to a stable operating powergrid, and thus, the powerflow equations are associated with a physically stable system. This stability permits the coarsening of the powerflow equations to be based on an approximate graph Laplacian, which is embedded in the powerflow system. By algebraically constructing a hierarchy of approximate weighted graph Laplacians, a hierarchy of nonlinear powerflow equations immediately becomes apparent. This latter hierarchy can then be used in a full approximation scheme (FAS) framework that leads to a nonlinear solver with generally a larger basin of attraction than Newton's method. Given the algebraic multigrid (AMG) coarsening of the approximate Laplacians, the solver is an AMG‐FAS scheme. Alternatively, using the coarse‐grid nodes and interpolation operators generated for the hierarchy of approximate graph Laplacians, a multiplicative‐correction scheme can be derived. The derivation of both schemes will be presented and analyzed, and numerical examples to demonstrate the performance of these schemes will be given.more » « less
-
Ward, C (Ed.)Computational microstructure design aims to fully exploit the precipitate strengthening potential of an alloy system. The development of accurate models to describe the temporal evolution of precipitate shapes and sizes is of great technological relevance. The experimental investigation of the precipitate microstructure is mostly based on two-dimensional micrographic images. Quantitative modeling of the temporal evolution of these microstructures needs to be discussed in three-dimensional simulation setups. To consistently bridge the gap between 2D images and 3D simulation data, we employ the method of central moments. Based on this, the aspect ratio of plate-like particles is consistently defined in two and three dimensions. The accuracy and interoperability of the method is demonstrated through representative 2D and 3D pixel-based sample data containing particles with a predefined aspect ratio. The applicability of the presented approach in integrated computational materials engineering (ICME) is demonstrated by the example of γ″ microstructure coarsening in Ni-based superalloys at 730 °C. For the first time, γ″ precipitate shape information from experimental 2D images and 3D phase-field simulation data is directly compared. This coarsening data indicates deviations from the classical ripening behavior and reveals periods of increased precipitate coagulation.more » « less
-
In a companion paper, we put forth a thermodynamic model for complex formation via a chemical reaction involving multiple macromolecular species, which may subsequently undergo liquid–liquid phase separation and a further transition into a gel-like state. In the present work, we formulate a thermodynamically consistent kinetic framework to study the interplay between phase separation, chemical reaction, and aging in spatially inhomogeneous macromolecular mixtures. A numerical algorithm is also proposed to simulate domain growth from collisions of liquid and gel domains via passive Brownian motion in both two and three spatial dimensions. Our results show that the coarsening behavior is significantly influenced by the degree of gelation and Brownian motion. The presence of a gel phase inside condensates strongly limits the diffusive transport processes, and Brownian motion coalescence controls the coarsening process in systems with high area/volume fractions of gel-like condensates, leading to the formation of interconnected domains with atypical domain growth rates controlled by size-dependent translational and rotational diffusivities.more » « less
-
Abstract The dynamics of active smectic liquid crystals confined on a spherical surface is explored through an active phase field crystal model. Starting from an initially randomly perturbed isotropic phase, several types of topological defects are spontaneously formed, and then annihilate during a coarsening process until a steady state is achieved. The coarsening process is highly complex involving several scaling laws of defect densities as a function of time where different dynamical exponents can be identified. In general the exponent for the final stage towards the steady state is significantly larger than that in the passive and in the planar case, i.e. the coarsening is getting accelerated both by activity and by the topological and geometrical properties of the sphere. A defect type characteristic for this active system is a rotating spiral of evolving smectic layering lines. On a sphere this defect type also determines the steady state. Our results can in principle be confirmed by dense systems of synthetic or biological active particles.more » « less
