BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Harald Garcke

Harald Garcke contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.APmath.NANumerical Analysismath.DGmath.OCmath.SPphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2026arXiv

A parametric finite element method for the incompressible Navier--Stokes equations on an evolving surface

In this paper we consider the numerical approximation of the incompressible surface Navier--Stokes equations on an evolving surface. For the discrete representation of the moving surface we use parametric finite elements of degree $\ell \geq 2$. In the semidiscrete continuous-in-time setting we are able to prove a stability estimate that mimics a corresponding result for the continuous problem. Some numerical results, including a convergence experiment, demonstrate the practicality and accuracy of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Convergence analysis for the Barrett--Garcke--Nurnberg method of transport type for evolving curves

In this paper, we propose a Barrett-Garcke-Nurnberg (BGN) method for evolving geometries under general flows and present the corresponding convergence analysis. Specifically, we examine the scenario where a closed curve evolves according to a prescribed background velocity field. Unlike mean curvature flow and surface diffusion, where the evolution velocities inherently exhibit parabolicity, this case is dominated by transport which poses a significant difficulty in establishing convergence proofs. To address the challenges imposed by this transport-dominant nature, we derive several discrete energy estimates of the transport type on discretized polynomial surfaces within the framework of the projection error. The use of the projection error is indispensable as it provides crucial additional stability through its orthogonality structure. We prove that the proposed method converges sub-optimally in the L2 norm, and this is the first convergence proof for a fully discrete numerical method solving the evolution of curves driven by general flows.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

On a Mullins-Sekerka model for the growth of active droplets modelling protocells: Stability analysis and numerical computations

Mullins-Sekerka models with chemical reactions can lead to scenarios where droplets grow, become unstable, split, grow and undergo further division. These grow and division cycles have been proposed as a model for protocells and are believed to play a fundamental role in living systems by providing chemical compartments which are important in the organization of living systems. This paper analyses chemically active Mullins-Sekerka models. Existence of radially symmetric solutions is shown and a detailed stability analysis in radial as well as planar situations is given. In particular, we also analyze multilayered solutions leading to shell-type situations. Finally, we introduce a numerical method based on a parametric finite element approach that explicitly accounts for topological changes, thereby allowing for droplet splitting and merging. Several numerical simulations verify the findings of the theoretical stability analysis and show complex dynamical behavior, including multiple instabilities, splittings of droplets and appearance of shell-type solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

A Cahn-Hilliard phase field model coupled to an Allen-Cahn model of viscoelasticity at large strains

We propose a new Cahn-Hilliard phase field model coupled to incompressible viscoelasticity at large strains, obtained from a diffuse interface mixture model and formulated in the Eulerian configuration. A new kind of diffusive regularization, of Allen-Cahn type, is introduced in the transport equation for the deformation gradient, together with a regularizing interface term depending on the gradient of the deformation gradient in the free energy of the system. We study the global existence of a weak solution for the model. While standard diffusive regularizations of the transport equation for the deformation gradient presented in literature allows the existence study only for simplified cases, i.e. in two space dimensions and for convex elastic free energy densities of Neo-Hookean type which are independent from the phase field variable, the present regularization allows to study more general cases. In particular, we obtain the global existence of a weak solution in three space dimensions and for generic nonlinear elastic energy densities with polynomial growth. Our analysis considers elastic free energy densities which depend on the phase field variable and which can possibly degenerate for some values of the phase field variable. By means of an iterative argument based on elliptic regularity bootstrap steps, we find the maximum allowed polynomial growths of the Cahn-Hilliard potential and the elastic energy density which guarantee the existence of a solution in three space dimensions. We propose two unconditionally energy stable finite element approximations of the model, based on convex splitting ideas and on the use of a scalar auxiliary variable, proving the existence and stability of discrete solutions. We finally report numerical results for different test cases with shape memory alloy type free energy with pure phases characterized by different elastic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Phase-Field Methods for Spectral Shape and Topology Optimization

We optimize a selection of eigenvalues of the Laplace operator with Dirichlet or Neumann boundary conditions by adjusting the shape of the domain on which the eigenvalue problem is considered. Here, a phase-field function is used to represent the shapes over which we minimize. The idea behind this method is to modify the Laplace operator by introducing phase-field dependent coefficients in order to extend the eigenvalue problem on a fixed design domain containing all admissible shapes. The resulting shape and topology optimization problem can then be formulated as an optimal control problem with PDE constraints in which the phase-field function acts as the control. For this optimal control problem, we establish first-order necessary optimality conditions and we rigorously derive its sharp interface limit. Eventually, we present and discuss several numerical simulations for our optimization problem.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Cahn-Hilliard model coupled to viscoelasticity with large deformations

We propose a new class of phase field models coupled to viscoelasticity with large deformations, obtained from a diffuse interface mixture model composed by a phase with elastic properties and a liquid phase. The model is formulated in the Eulerian configuration and it is derived by imposing the mass balance for the mixture components and the momentum balance that comes from a generalized form of the principle of virtual powers. The latter considers the presence of a system of microforces and microstresses associated to the microscopic interactions between the mixture's constituents together with a system of macroforces and macrostresses associated to their viscoelastic behavior, taking into account also the friction between the phases. The free energy density of the system is given as the sum of a Cahn-Hilliard term and an elastic polyconvex term, with a coupling between the phase field variable and the elastic deformation gradient in the elastic contribution. General constitutive assumptions complying with a mechanical version of the second law of thermodynamics in isothermal situations are taken. We study the global existence of a weak solution for a simplified and regularized version of the general model, which considers an incompressible elastic free energy of Neo-Hookean type with elastic coefficients depending on the phase field variable. The regularization is properly designed to deal with the coupling between the phase field variable and the elastic deformation gradient in the elastic energy density. The analysis is made both in two and three space dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical analysis for a Cahn-Hilliard system modelling tumour growth with chemotaxis and active transport

In this work, we consider a diffuse interface model for tumour growth in the presence of a nutrient which is consumed by the tumour. The system of equations consists of a Cahn--Hilliard equation with source terms for the tumour cells and a reaction-diffusion equation for the nutrient. We introduce a fully-discrete finite element approximation of the model and prove stability bounds for the discrete scheme. Moreover, we show that discrete solutions exist and depend continuously on the initial and boundary data. We then pass to the limit in the discretization parameters and prove convergence to a global-in-time weak solution to the model. Under additional assumptions, this weak solution is unique. Finally, we present some numerical results including numerical error investigation in one spatial dimension and some long time simulations in two and three spatial dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical analysis for the interaction of mean curvature flow and diffusion on closed surfaces

An evolving surface finite element discretisation is analysed for the evolution of a closed two-dimensional surface governed by a system coupling a generalised forced mean curvature flow and a reaction--diffusion process on the surface, inspired by a gradient flow of a coupled energy. Two algorithms are proposed, both based on a system coupling the diffusion equation to evolution equations for geometric quantities in the velocity law for the surface. One of the numerical methods is proved to be convergent in the $H^1$ norm with optimal-order for finite elements of degree at least two. We present numerical experiments illustrating the convergence behaviour and demonstrating the qualitative properties of the flow: preservation of mean convexity, loss of convexity, weak maximum principles, and the occurrence of self-intersections.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Qualitative Properties for a System Coupling Scaled Mean Curvature Flow and Diffusion

We consider a system consisting of a geometric evolution equation for a hypersurface and a parabolic equation on this evolving hypersurface. More precisely, we discuss mean curvature flow scaled with a term that depends on a quantity defined on the surface coupled to a diffusion equation for that quantity. Several properties of solutions are analyzed. Emphasis is placed on to what extent the surface in our setting qualitatively evolves similar as for the usual mean curvature flow. To this end, we show that the surface area is strictly decreasing but give an example of a surface that exists for infinite times nevertheless. Moreover, mean convexity is conserved whereas convexity is not. Finally, we construct an embedded hypersurface that develops a self-intersection in the course of time. Additionally, a formal explanation of how our equations can be interpreted as a gradient flow is included.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Short Time Existence for Coupling of Scaled Mean Curvature Flow and Diffusion

We prove a short time existence result for a system consisting of a geometric evolution equation for a hypersurface and a parabolic equation on this evolving hypersurface. More precisely, we discuss a mean curvature flow scaled with a term that depends on a quantity defined on the surface coupled to a diffusion equation for that quantity. The proof is based on a splitting ansatz, solving both equations separately using linearization and a contraction argument. Our result is formulated for the case of immersed hypersurfaces and yields a uniform lower bound on the existence time that allows for small changes in the initial value of the height function.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Volume-preserving parametric finite element methods for axisymmetric geometric evolution equations

We propose and analyze volume-preserving parametric finite element methods for surface diffusion, conserved mean curvature flow and an intermediate evolution law in an axisymmetric setting. The weak formulations are presented in terms of the generating curves of the axisymmetric surfaces. The proposed numerical methods are based on piecewise linear parametric finite elements. The constructed fully practical schemes satisfy the conservation of the enclosed volume. In addition, we prove the unconditional stability and consider the distribution of vertices for the discretized schemes. The introduced methods are implicit and the resulting nonlinear systems of equations can be solved very efficiently and accurately via the Newton's iterative method. Numerical results are presented to show the accuracy and efficiency of the introduced schemes for computing the considered axisymmetric geometric flows.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Long-time Dynamics for a Cahn-Hilliard Tumor Growth Model with Chemotaxis

Mathematical models that describe the tumor growth process have been formulated by several authors in order to understand how cancer develops and to develop new treatment approaches. In this study, it is aimed to investigate the long-time behavior of the two-phase diffuse-interface model, which was proposed in [11] to model a tumor tissue as a mixture of cancerous and healthy cells. Up to now, studies on the long-time behavior of this model have neglected chemotaxis and active transport, which have a significant effect on tumor growth. In this research, our main aim is to study this model with chemotaxis and active transport. We prove an asymptotic compactness result for the weak solutions of the problem in the whole phase-space $H^1 (Ω) \times L^2 (Ω)$. We establish the existence of the global attractor in a phase space denoted via mass conservation. We also prove that the global attractor equals to the unstable manifold emanating from the set of stationary points. Moreover, we obtain that the global attractor has a finite fractal dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On a phase field model of Cahn-Hilliard type for tumour growth with mechanical effects

Mechanical effects have mostly been neglected so far in phase field tumour models that are based on a Cahn-Hilliard approach. In this paper we study a macroscopic mechanical model for tumour growth in which cell-cell adhesion effects are taken into account with the help of a Ginzburg--Landau type energy. In the overall model an equation of Cahn-Hilliard type is coupled to the system of linear elasticity and a reaction-diffusion equation for a nutrient concentration. The highly non-linear coupling between a fourth-order Cahn-Hilliard equation and the quasi-static elasticity system lead to new challenges which cannot be dealt within a gradient flow setting which was the method of choice for other elastic Cahn-Hilliard systems. We show existence, uniqueness and regularity results. In addition, several continuous dependence results with respect to different topologies are shown. Some of these results give uniqueness for weak solutions and other results will be helpful for optimal control problems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Parametric finite element approximations of curvature driven interface evolutions

Parametric finite elements lead to very efficient numerical methods for surface evolution equations. We introduce several computational techniques for curvature driven evolution equations based on a weak formulation for the mean curvature. The approaches discussed, in contrast to many other methods, have good mesh properties that avoid mesh coalescence and very non-uniform meshes. Mean curvature flow, surface diffusion, anisotropic geometric flows, solidification, two-phase flow, Willmore and Helfrich flow as well as biomembranes are treated. We show stability results as well as results explaining the good mesh properties.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Numerical approximation of curve evolutions in Riemannian manifolds

We introduce variational approximations for curve evolutions in two-dimensional Riemannian manifolds that are conformally flat, i.e.\ conformally equivalent to the Euclidean space. Examples include the hyperbolic plane, the hyperbolic disk, the elliptic plane as well as any conformal parameterization of a two-dimensional surface in ${\mathbb R}^d$, $d\geq 3$. In these spaces we introduce stable numerical schemes for curvature flow and curve diffusion, and we also formulate a scheme for elastic flow. Variants of the schemes can also be applied to geometric evolution equations for axisymmetric hypersurfaces in ${\mathbb R}^d$. Some of the schemes have very good properties with respect to the distribution of mesh points, which is demonstrated with the help of several numerical computations.

0 citations0 reviewsNo code linkedOpen access