-
BiLO: Bilevel Local Operator Learning for PDE inverse problems
Authors:
Ray Zirui Zhang,
Xiaohui Xie,
John S. Lowengrub
Abstract:
We propose a new neural network based method for solving inverse problems for partial differential equations (PDEs) by formulating the PDE inverse problem as a bilevel optimization problem. At the upper level, we minimize the data loss with respect to the PDE parameters. At the lower level, we train a neural network to locally approximate the PDE solution operator in the neighborhood of a given se…
▽ More
We propose a new neural network based method for solving inverse problems for partial differential equations (PDEs) by formulating the PDE inverse problem as a bilevel optimization problem. At the upper level, we minimize the data loss with respect to the PDE parameters. At the lower level, we train a neural network to locally approximate the PDE solution operator in the neighborhood of a given set of PDE parameters, which enables an accurate approximation of the descent direction for the upper level optimization problem. The lower level loss function includes the L2 norms of both the residual and its derivative with respect to the PDE parameters. We apply gradient descent simultaneously on both the upper and lower level optimization problems, leading to an effective and fast algorithm. The method, which we refer to as BiLO (Bilevel Local Operator learning), is also able to efficiently infer unknown functions in the PDEs through the introduction of an auxiliary variable. Through extensive experiments over multiple PDE systems, we demonstrate that our method enforces strong PDE constraints, is robust to sparse and noisy data, and eliminates the need to balance the residual and the data loss, which is inherent to the soft PDE constraints in many existing methods.
△ Less
Submitted 15 August, 2024; v1 submitted 27 April, 2024;
originally announced April 2024.
-
Individualizing Glioma Radiotherapy Planning by Optimization of Data and Physics-Informed Discrete Loss
Authors:
Michal Balcerak,
Jonas Weidner,
Petr Karnakov,
Ivan Ezhov,
Sergey Litvinov,
Petros Koumoutsakos,
Ray Zirui Zhang,
John S. Lowengrub,
Bene Wiestler,
Bjoern Menze
Abstract:
Brain tumor growth is unique to each glioma patient and extends beyond what is visible in imaging scans, infiltrating surrounding brain tissue. Understanding these hidden patient-specific progressions is essential for effective therapies. Current treatment plans for brain tumors, such as radiotherapy, typically involve delineating a uniform margin around the visible tumor on pre-treatment scans to…
▽ More
Brain tumor growth is unique to each glioma patient and extends beyond what is visible in imaging scans, infiltrating surrounding brain tissue. Understanding these hidden patient-specific progressions is essential for effective therapies. Current treatment plans for brain tumors, such as radiotherapy, typically involve delineating a uniform margin around the visible tumor on pre-treatment scans to target this invisible tumor growth. This "one size fits all" approach is derived from population studies and often fails to account for the nuances of individual patient conditions. We present the GliODIL framework, which infers the full spatial distribution of tumor cell concentration from available multi-modal imaging, leveraging a Fisher-Kolmogorov type physics model to describe tumor growth. This is achieved through the newly introduced method of Optimizing the Discrete Loss (ODIL), where both data and physics-based constraints are softly assimilated into the solution. Our test dataset comprises 152 glioblastoma patients with pre-treatment imaging and post-treatment follow-ups for tumor recurrence monitoring. By blending data-driven techniques with physics-based constraints, GliODIL enhances recurrence prediction in radiotherapy planning, challenging traditional uniform margins and strict adherence to the Fisher-Kolmogorov partial differential equation (PDE) model, which is adapted for complex cases.
△ Less
Submitted 14 April, 2024; v1 submitted 8 December, 2023;
originally announced December 2023.
-
Phase field modeling and computation of vesicle growth or shrinkage
Authors:
Xiaoxia Tang,
Shuwang Li,
John S. Lowengrub,
Steven M. Wise
Abstract:
We present a phase field model for vesicle growth or shrinkage induced by an osmotic pressure due to a chemical potential gradient. The model consists of an Allen-Cahn equation describing the evolution of phase field and a Cahn-Hilliard equation describing the evolution of concentration field. We establish control conditions for vesicle growth or shrinkage via a common tangent construction. During…
▽ More
We present a phase field model for vesicle growth or shrinkage induced by an osmotic pressure due to a chemical potential gradient. The model consists of an Allen-Cahn equation describing the evolution of phase field and a Cahn-Hilliard equation describing the evolution of concentration field. We establish control conditions for vesicle growth or shrinkage via a common tangent construction. During the membrane deformation, the model ensures total mass conservation and satisfies surface area constraint. We develop a nonlinear numerical scheme, a combination of nonlinear Gauss-Seidel relaxation operator and a V-cycles multigrid solver, for computing equilibrium shapes of a 2D vesicle. Convergence tests confirm an $\mathcal{O}(t+h^2)$ accuracy. Numerical results reveal that the diffuse interface model captures the main feature of dynamics: for a growing vesicle, there exist circle-like equilibrium shapes if the concentration difference across the membrane and the initial osmotic pressure are large enough; while for a shrinking vesicle, there exists a rich collection of finger-like equilibrium morphologies.
△ Less
Submitted 27 September, 2022;
originally announced September 2022.
-
Numerical study on viscous fingering using electric fields in a Hele-Shaw cell
Authors:
Meng Zhao,
Pedro Anjos,
John Lowengrub,
Wenjun Ying,
Shuwang Li
Abstract:
We investigate the nonlinear dynamics of a moving interface in a Hele-Shaw cell subject to an in-plane applied electric field. We develop a spectrally accurate boundary integral method where a coupled integral equation system is formulated. Although the stiffness due to the high order spatial derivatives can be removed, the long-time simulation is still expensive since the evolving velocity of the…
▽ More
We investigate the nonlinear dynamics of a moving interface in a Hele-Shaw cell subject to an in-plane applied electric field. We develop a spectrally accurate boundary integral method where a coupled integral equation system is formulated. Although the stiffness due to the high order spatial derivatives can be removed, the long-time simulation is still expensive since the evolving velocity of the interface drops dramatically as the interface expands. We remove this physically imposed stiffness by employing a rescaling scheme, which accelerates the slow dynamics and reduces the computational cost. Our nonlinear results reveal that positive currents restrain finger ramification and promote overall stabilization of patterns. On the other hand, negative currents make the interface more unstable and lead to the formation of thin tail structures connecting the fingers and a small inner region. When no flux is injected, and a negative current is utilized, the interface tends to approach the origin and break up into several drops. We investigate the temporal evolution of the smallest distance between the interface and the origin and find that it obeys an algebraic law $\displaystyle (t_*-t)^b$, where $t_*$ is the estimated pinch-off time.
△ Less
Submitted 15 January, 2022;
originally announced January 2022.
-
Nonlinear simulation of vascular tumor growth with chemotaxis and the control of necrosis
Authors:
Min-Jhe Lu,
Wenrui Hao,
Chun Liu,
John Lowengrub,
Shuwang Li
Abstract:
In this paper, we develop a sharp interface tumor growth model in two dimensions to study the effect of both the intratumoral structure using a controlled necrotic core and the extratumoral nutrient supply from vasculature on tumor morphology. We first show that our model extends the benchmark results in the literature using linear stability analysis. Then we solve this generalized model numerical…
▽ More
In this paper, we develop a sharp interface tumor growth model in two dimensions to study the effect of both the intratumoral structure using a controlled necrotic core and the extratumoral nutrient supply from vasculature on tumor morphology. We first show that our model extends the benchmark results in the literature using linear stability analysis. Then we solve this generalized model numerically using a spectrally accurate boundary integral method in an evolving annular domain, not only with a Robin boundary condition on the outer boundary for the nutrient field which models tumor vasculature, but also with a static boundary condition on the inner boundary for pressure field which models the control of tumor necrosis. The discretized linear systems for both pressure and nutrient fields are shown to be well-conditioned through tracing GMRES iteration numbers. Our nonlinear simulations reveal the stabilizing effects of angiogenesis and the destabilizing ones of chemotaxis and necrosis in the development of tumor morphological instabilities if the necrotic core is fixed in a circular shape. When the necrotic core is controlled in a non-circular shape, the stabilizing effects of proliferation and the destabilizing ones of apoptosis are observed. Finally, the values of the nutrient concentration with its fluxes and the pressure level with its normal derivatives, which are solved accurately at the boundaries, help us to characterize the corresponding tumor morphology and the level of the biophysical quantities on interfaces required in keeping various shapes of the necrotic region of the tumor. Interestingly, we notice that when the necrotic region is fixed in a 3-fold non-circular shape, even if the initial shape of the tumor is circular, the tumor will evolve into a shape corresponding to the 3-fold symmetry of the shape of the fixed necrotic region.
△ Less
Submitted 22 March, 2022; v1 submitted 24 July, 2021;
originally announced July 2021.
-
Sharp-interface problem of the Ohta-Kawasaki model for symmetric diblock copolymers
Authors:
Amlan K. Barua,
Ray Chew,
Shuwang Li,
John Lowengrub,
Andreas Münch,
Barbara Wagner
Abstract:
The Ohta-Kawasaki model for diblock-copolymers is well known to the scientific community of diffuse-interface methods. To accurately capture the long-time evolution of the moving interfaces, we present a derivation of the corresponding sharp-interface limit using matched asymptotic expansions, and show that the limiting process leads to a Hele-Shaw type moving interface problem. The numerical trea…
▽ More
The Ohta-Kawasaki model for diblock-copolymers is well known to the scientific community of diffuse-interface methods. To accurately capture the long-time evolution of the moving interfaces, we present a derivation of the corresponding sharp-interface limit using matched asymptotic expansions, and show that the limiting process leads to a Hele-Shaw type moving interface problem. The numerical treatment of the sharp-interface limit is more complicated due to the stiffness of the equations. To address this problem, we present a boundary integral formulation corresponding to a sharp interface limit of the Ohta-Kawasaki model. Starting with the governing equations defined on separate phase domains, we develop boundary integral equations valid for multi-connected domains in a 2D plane. For numerical simplicity we assume our problem is driven by a uniform Dirichlet condition on a circular far-field boundary. The integral formulation of the problem involves both double- and single-layer potentials due to the modified boundary condition. In particular, our formulation allows one to compute the nonlinear dynamics of a non-equilibrium system and pattern formation of an equilibrating system. Numerical tests on an evolving slightly perturbed circular interface (separating the two phases) are in excellent agreement with the linear analysis, demonstrating that the method is stable, efficient and spectrally accurate in space.
△ Less
Submitted 9 July, 2021;
originally announced July 2021.
-
Complex Far-Field Geometries Determine the Stability of Solid Tumor Growth with Chemotaxis
Authors:
Min-Jhe Lu,
Chun Liu,
John Lowengrub,
Shuwang Li
Abstract:
In this paper, we develop a sharp interface tumor growth model to study the effect of the tumor microenvironment using a complex far-field geometry that mimics a heterogeneous distribution of vasculature. Together with different nutrient uptake rates inside and outside the tumor, this introduces variability in spatial diffusion gradients. Linear stability analysis suggests that the uptake rate in…
▽ More
In this paper, we develop a sharp interface tumor growth model to study the effect of the tumor microenvironment using a complex far-field geometry that mimics a heterogeneous distribution of vasculature. Together with different nutrient uptake rates inside and outside the tumor, this introduces variability in spatial diffusion gradients. Linear stability analysis suggests that the uptake rate in the tumor microenvironment, together with chemotaxis, may induce unstable growth, especially when the nutrient gradients are large. We investigate the fully nonlinear dynamics using a spectrally accurate boundary integral method. Our nonlinear simulations reveal that vascular heterogeneity plays an important role in the development of morphological instabilities that range from fingering and chain-like morphologies to compact, plate-like shapes in two-dimensions.
△ Less
Submitted 11 March, 2020; v1 submitted 10 December, 2019;
originally announced December 2019.
-
Higher-order accurate diffuse-domain methods for partial differential equations with Dirichlet boundary conditions in complex, evolving geometries
Authors:
Fei Yu,
Zhenlin Guo,
John Lowengrub
Abstract:
The diffuse-domain, or smoothed boundary, method is an attractive approach for solving partial differential equations in complex geometries because of its simplicity and flexibility. In this method the complex geometry is embedded into a larger, regular domain. The original PDE is reformulated using a smoothed characteristic function of the complex domain and source terms are introduced to approxi…
▽ More
The diffuse-domain, or smoothed boundary, method is an attractive approach for solving partial differential equations in complex geometries because of its simplicity and flexibility. In this method the complex geometry is embedded into a larger, regular domain. The original PDE is reformulated using a smoothed characteristic function of the complex domain and source terms are introduced to approximate the boundary conditions. The reformulated equation, which is independent of the dimension and domain geometry, can be solved by standard numerical methods and the same solver can be used for any domain geometry. A challenge is making the method higher-order accurate. For Dirichlet boundary conditions, which we focus on here, current implementations demonstrate a wide range in their accuracy but generally the methods yield at best first order accuracy in $ε$, the parameter that characterizes the width of the region over which the characteristic function is smoothed. Typically, $ε\propto h$, the grid size. Here, we analyze the diffuse-domain PDEs using matched asymptotic expansions and explain the observed behaviors. Our analysis also identifies simple modifications to the diffuse-domain PDEs that yield higher-order accuracy in $ε$, e.g., $O(ε^2)$ in the $L^2$ norm and $O(ε^p)$ with $1.5\le p\le 2$ in the $L^{\infty}$ norm. Our analytic results are confirmed numerically in stationary and moving domains where the level set method is used to capture the dynamics of the domain boundary and to construct the smoothed characteristic function.
△ Less
Submitted 27 November, 2019; v1 submitted 9 December, 2018;
originally announced December 2018.
-
Non-stiff methods for Airy flow and the modified Korteweg-de Vries equation
Authors:
Mariano Franco-de-León,
John Lowengrub
Abstract:
In this paper, we implement non-stiff interface tracking methods for the evolution of 2-D curves that follow Airy flow, a curvature-dependent dispersive geometric evolution law. The curvature of the curve satisfies the modified Korteweg-de Vries equation, a dispersive non-linear soliton equation. We present a fully discrete space-time analysis of the equations (proof of convergence) and numerical…
▽ More
In this paper, we implement non-stiff interface tracking methods for the evolution of 2-D curves that follow Airy flow, a curvature-dependent dispersive geometric evolution law. The curvature of the curve satisfies the modified Korteweg-de Vries equation, a dispersive non-linear soliton equation. We present a fully discrete space-time analysis of the equations (proof of convergence) and numerical evidence that confirms the accuracy, convergence, efficiency, and stability of the methods.
△ Less
Submitted 30 August, 2017;
originally announced August 2017.
-
Mass Conservative and Energy Stable Finite Difference Methods for the Quasi-incompressible Navier-Stokes-Cahn-Hilliard system: Primitive Variable and Projection-Type Schemes
Authors:
Zhenlin Guo,
Ping Lin,
Steven Wise,
John Lowengrub
Abstract:
In this paper we describe two fully mass conservative, energy stable, finite difference methods on a staggered grid for the quasi-incompressible Navier-Stokes-Cahn-Hilliard (q-NSCH) system governing a binary incompressible fluid flow with variable density and viscosity. Both methods, namely the primitive method (finite difference method in the primitive variable formulation) and the projection met…
▽ More
In this paper we describe two fully mass conservative, energy stable, finite difference methods on a staggered grid for the quasi-incompressible Navier-Stokes-Cahn-Hilliard (q-NSCH) system governing a binary incompressible fluid flow with variable density and viscosity. Both methods, namely the primitive method (finite difference method in the primitive variable formulation) and the projection method (finite difference method in a projection-type formulation), are so designed that the mass of the binary fluid is preserved, and the energy of the system equations is always non-increasing in time at the fully discrete level. We also present an efficient, practical nonlinear multigrid method - comprised of a standard FAS method for the Cahn-Hilliard equation, and a method based on the Vanka-type smoothing strategy for the Navier-Stokes equation - for solving these equations. We test the scheme in the context of Capillary Waves, rising droplets and Rayleigh-Taylor instability. Quantitative comparisons are made with existing analytical solutions or previous numerical results that validate the accuracy of our numerical schemes. Moreover, in all cases, mass of the single component and the binary fluid was conserved up to 10 to -8 and energy decreases in time.
△ Less
Submitted 12 August, 2017; v1 submitted 20 March, 2017;
originally announced March 2017.
-
Convergence Analysis for Second Order Accurate Convex Splitting Schemes for the Periodic Nonlocal Allen-Cahn and Cahn-Hilliard Equations
Authors:
Zhen Guan,
John Lowengrub,
Cheng Wang
Abstract:
In this paper we provide a detailed convergence analysis for fully discrete second order (in both time and space) numerical schemes for nonlocal Allen-Cahn (nAC) and nonlocal Cahn-Hilliard (nCH) equations. The unconditional unique solvability and energy stability ensures $\ell^4$ stability. The convergence analysis for the nAC equation follows the standard procedure of consistency and stability es…
▽ More
In this paper we provide a detailed convergence analysis for fully discrete second order (in both time and space) numerical schemes for nonlocal Allen-Cahn (nAC) and nonlocal Cahn-Hilliard (nCH) equations. The unconditional unique solvability and energy stability ensures $\ell^4$ stability. The convergence analysis for the nAC equation follows the standard procedure of consistency and stability estimate for the numerical error function. For the nCH equation, due to the complicated form of the nonlinear term, a careful expansion of its discrete gradient is undertaken and an $H^{-1}$ inner product estimate of this nonlinear numerical error is derived to establish convergence. In addition, an a-priori $W^{1,\infty}$ bound of the numerical solution at the discrete level is needed in the error estimate. Such a bound can be obtained by performing a higher order consistency analysis by using asymptotic expansions for the numerical solution. Following the technique originally proposed by Strang (e.g., 1964), instead of the standard comparison between the exact and numerical solutions, an error estimate between the numerical solution and the constructed approximate solution yields an $O( s^3 + h^4)$ convergence in $\ell^\infty (0, T; \ell^2)$ norm, which leads to the necessary bound under a standard constraint $s \le C h$. Here, we also prove convergence of the scheme in the maximum norm under the same constraint.
△ Less
Submitted 21 October, 2016;
originally announced October 2016.
-
An Energy Stable Finite-Difference Scheme for Functionalized Cahn-Hilliard Equation and its Convergence Analysis
Authors:
Wenqiang Feng,
Zhen Guan,
John Lowengrub,
Cheng Wang,
Steven M. Wise
Abstract:
We present and analyze an unconditionally energy stable and convergent finite difference scheme for the Functionalized Cahn-Hilliard equation. One key difficulty associated with the energy stability is based on the fact that one nonlinear energy functional term in the expansion appears as non-convex, non-concave. To overcome this subtle difficulty, we add two auxiliary terms to make the combined t…
▽ More
We present and analyze an unconditionally energy stable and convergent finite difference scheme for the Functionalized Cahn-Hilliard equation. One key difficulty associated with the energy stability is based on the fact that one nonlinear energy functional term in the expansion appears as non-convex, non-concave. To overcome this subtle difficulty, we add two auxiliary terms to make the combined term convex, which in turns yields a convex-concave decomposition of the physical energy. As a result, an application of the convex splitting methodology assures both the unique solvability and the unconditional energy stability of the proposed numerical scheme. To deal with a 4-Laplacian solver in an $H^{-1}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm to solve the corresponding nonlinear systems. In addition, a global in time $H_{\rm per}^2$ stability of the numerical scheme is established at a theoretical level, which in turn ensures the full order convergence analysis of the scheme. A few numerical results are presented, which confirm the stability and accuracy of the proposed numerical scheme.
△ Less
Submitted 7 October, 2016;
originally announced October 2016.
-
Energy Stable Multigrid Method for Local and Non-local Hydrodynamic Models for Freezing
Authors:
Arvind Baskaran,
Zhen Guan,
John Lowengrub
Abstract:
In this paper we present a numerical method for hydrodynamic models that arise from time dependent density functional theories of freezing. The models take the form of compressible Navier-Stokes equations whose pressure is determined by the variational derivative of a free energy, which is a functional of the density field. We present unconditionally energy stable and mass conserving implicit fini…
▽ More
In this paper we present a numerical method for hydrodynamic models that arise from time dependent density functional theories of freezing. The models take the form of compressible Navier-Stokes equations whose pressure is determined by the variational derivative of a free energy, which is a functional of the density field. We present unconditionally energy stable and mass conserving implicit finite difference methods for the models. The methods are based on a convex splitting of the free energy and that ensures that a discrete energy is non-increasing for any choice of time and space step. The methods are applicable to a large class of models, including both local and non-local free energy functionals. The theoretical basis for the numerical method is presented in a general context. The method is applied to problems using two specific free energy functionals: one local and one non-local functional. A nonlinear multigrid method is used to solve the numerical method, which is nonlinear at the implicit time step. The non-local functional, which is a convolution operator, is approximated using the Discrete Fourier Transform. Numerical simulations that confirm the stability and accuracy of the numerical method are presented.
△ Less
Submitted 22 October, 2015; v1 submitted 25 March, 2015;
originally announced March 2015.
-
Analysis of the diffuse-domain method for solving PDEs in complex geometries
Authors:
Karl Yngve Lervåg,
John Lowengrub
Abstract:
In recent work, Li et al.\ (Comm.\ Math.\ Sci., 7:81-107, 2009) developed a diffuse-domain method (DDM) for solving partial differential equations in complex, dynamic geometries with Dirichlet, Neumann, and Robin boundary conditions. The diffuse-domain method uses an implicit representation of the geometry where the sharp boundary is replaced by a diffuse layer with thickness $ε$ that is typically…
▽ More
In recent work, Li et al.\ (Comm.\ Math.\ Sci., 7:81-107, 2009) developed a diffuse-domain method (DDM) for solving partial differential equations in complex, dynamic geometries with Dirichlet, Neumann, and Robin boundary conditions. The diffuse-domain method uses an implicit representation of the geometry where the sharp boundary is replaced by a diffuse layer with thickness $ε$ that is typically proportional to the minimum grid size. The original equations are reformulated on a larger regular domain and the boundary conditions are incorporated via singular source terms. The resulting equations can be solved with standard finite difference and finite element software packages. Here, we present a matched asymptotic analysis of general diffuse-domain methods for Neumann and Robin boundary conditions. Our analysis shows that for certain choices of the boundary condition approximations, the DDM is second-order accurate in $ε$. However, for other choices the DDM is only first-order accurate. This helps to explain why the choice of boundary-condition approximation is important for rapid global convergence and high accuracy. Our analysis also suggests correction terms that may be added to yield more accurate diffuse-domain methods. Simple modifications of first-order boundary condition approximations are proposed to achieve asymptotically second-order accurate schemes. Our analytic results are confirmed numerically in the $L^2$ and $L^\infty$ norms for selected test problems.
△ Less
Submitted 15 May, 2015; v1 submitted 25 July, 2014;
originally announced July 2014.
-
A numerical method for the quasi-incompressible Cahn-Hilliard-Navier-Stokes equations for variable density flows with a discrete energy law
Authors:
Zhenlin Guo,
Ping Lin,
John S. Lowengrub
Abstract:
In this paper, we investigate numerically a diffuse interface model for the Navier-Stokes equation with fluid-fluid interface when the fluids have different densities \cite{Lowengrub1998}. Under minor reformulation of the system, we show that there is a continuous energy law underlying the system, assuming that all variables have reasonable regularities. It is shown in the literature that an energ…
▽ More
In this paper, we investigate numerically a diffuse interface model for the Navier-Stokes equation with fluid-fluid interface when the fluids have different densities \cite{Lowengrub1998}. Under minor reformulation of the system, we show that there is a continuous energy law underlying the system, assuming that all variables have reasonable regularities. It is shown in the literature that an energy law preserving method will perform better for multiphase problems. Thus for the reformulated system, we design a $C^0$ finite element method and a special temporal scheme where the energy law is preserved at the discrete level. Such a discrete energy law (almost the same as the continuous energy law) for this variable density two-phase flow model has never been established before with $C^0$ finite element. A Newton's method is introduced to linearise the highly non-linear system of our discretization scheme. Some numerical experiments are carried out using the adaptive mesh to investigate the scenario of coalescing and rising drops with differing density ratio. The snapshots for the evolution of the interface together with the adaptive mesh at different times are presented to show that the evolution, including the break-up/pinch-off of the drop, can be handled smoothly by our numerical scheme. The discrete energy functional for the system is examined to show that the energy law at the discrete level is preserved by our scheme.
△ Less
Submitted 19 February, 2014; v1 submitted 6 February, 2014;
originally announced February 2014.
-
Energy Stable and Efficient Finite-Difference Nonlinear Multigrid Schemes for the Modified Phase Field Crystal Equation
Authors:
Arvind Baskaran,
Peng Zhou,
Zhengzheng Hu,
Cheng Wang,
Steven M. Wise,
John S. Lowengrub
Abstract:
In this paper we present two unconditionally energy stable finite difference schemes for the Modified Phase Field Crystal (MPFC) equation, a sixth-order nonlinear damped wave equation, of which the purely parabolic Phase Field Crystal (PFC) model can be viewed as a special case. The first is a convex splitting scheme based on an appropriate decomposition of the discrete energy and is first order a…
▽ More
In this paper we present two unconditionally energy stable finite difference schemes for the Modified Phase Field Crystal (MPFC) equation, a sixth-order nonlinear damped wave equation, of which the purely parabolic Phase Field Crystal (PFC) model can be viewed as a special case. The first is a convex splitting scheme based on an appropriate decomposition of the discrete energy and is first order accurate in time and second order accurate in space. The second is a new, fully second-order scheme that also respects the convex splitting of the energy. Both schemes are nonlinear but may be formulated from the gradients of strictly convex, coercive functionals. Thus, both are uniquely solvable regardless of the time and space step sizes. The schemes are solved by efficient nonlinear multigrid methods. Numerical results are presented demonstrating the accuracy, energy stability, efficiency, and practical utility of the schemes. In particular, we show that our multigrid solvers enjoy optimal, or nearly optimal complexity in the solution of the nonlinear schemes.
△ Less
Submitted 8 August, 2012;
originally announced August 2012.
-
Convergence Analysis of a Second Order Convex Splitting Scheme for the Modified Phase Field Crystal Equation
Authors:
Arvind Baskaran,
John Lowengrub,
Cheng Wang,
Steve Wise
Abstract:
In this paper we provide a detailed convergence analysis for an unconditionally energy stable, second-order accurate convex splitting scheme for the Modified Phase Field Crystal equation, a generalized damped wave equation for which the usual Phase Field Crystal equation is a special degenerate case. The fully discrete, fully second-order finite difference scheme in question was derived in a recen…
▽ More
In this paper we provide a detailed convergence analysis for an unconditionally energy stable, second-order accurate convex splitting scheme for the Modified Phase Field Crystal equation, a generalized damped wave equation for which the usual Phase Field Crystal equation is a special degenerate case. The fully discrete, fully second-order finite difference scheme in question was derived in a recent work [2]. An introduction of a new variable ψ, corresponding to the temporal derivative of the phase variable φ, could bring an accuracy reduction in the formal consistency estimate, because of the hyperbolic nature of the equation. A higher order consistency analysis by an asymptotic expansion is performed to overcome this difficulty. In turn, second order convergence in both time and space is established in a discrete L^\infty (0,T; H^3) norm.
△ Less
Submitted 7 August, 2012;
originally announced August 2012.
-
Analysis of a Mixture Model of Tumor Growth
Authors:
John Lowengrub,
Edriss S. Titi,
Kun Zhao
Abstract:
We study an initial-boundary value problem (IBVP) for a coupled Cahn-Hilliard-Hele-Shaw system that models tumor growth. For large initial data with finite energy, we prove global (local resp.) existence, uniqueness, higher order spatial regularity and Gevrey spatial regularity of strong solutions to the IBVP in 2D (3D resp.). Asymptotically in time, we show that the solution converges to a consta…
▽ More
We study an initial-boundary value problem (IBVP) for a coupled Cahn-Hilliard-Hele-Shaw system that models tumor growth. For large initial data with finite energy, we prove global (local resp.) existence, uniqueness, higher order spatial regularity and Gevrey spatial regularity of strong solutions to the IBVP in 2D (3D resp.). Asymptotically in time, we show that the solution converges to a constant state exponentially fast as time tends to infinity under certain assumptions.
△ Less
Submitted 30 May, 2012;
originally announced May 2012.