BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ricardo H. Nochetto

Ricardo H. Nochetto appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.OCNumerical Analysismath.APgr-qc

Catalog footprint

What is connected

31works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

31 published item(s)

preprint2022arXiv

Fractional Elliptic Problems on Lipschitz Domains: Regularity and Approximation

This survey hinges on the interplay between regularity and approximation for linear and quasi-linear fractional elliptic problems on Lipschitz domains. For the linear Dirichlet integral Laplacian, after briefly recalling Hölder regularity and applications, we discuss novel optimal shift theorems in Besov spaces and their Sobolev counterparts. These results extend to problems with finite horizon and are instrumental for the subsequent error analysis. Moreover, we dwell on extensions of Besov regularity to the fractional $p$-Laplacian, and review the regularity of fractional minimal graphs and stickiness. We discretize these problems using continuous piecewise linear finite elements and derive global and local error estimates for linear problems, thereby improving some existing error estimates for both quasi-uniform and graded meshes. We also present a BPX preconditioner which turns out to be robust with respect to both the fractional order and the number of levels. We conclude with the discretization of fractional quasi-linear problems and their error analysis. We illustrate the theory with several illuminating numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finite element discretizations of nonlocal minimal graphs: convergence

In this paper, we propose and analyze a finite element discretization for the computation of fractional minimal graphs of order~$s \in (0,1/2)$ on a bounded domain $Ω$. Such a Plateau problem of order $s$ can be reinterpreted as a Dirichlet problem for a nonlocal, nonlinear, degenerate operator of order $s + 1/2$. We prove that our numerical scheme converges in $W^{2r}_1(Ω)$ for all $r<s$, where $W^{2s}_1(Ω)$ is closely related to the natural energy space. Moreover, we introduce a geometric notion of error that, for any pair of $H^1$ functions, in the limit $s \to 1/2$ recovers a weighted $L^2$-discrepancy between the normal vectors to their graphs. We derive error bounds with respect to this novel geometric quantity as well. In spite of performing approximations with continuous, piecewise linear, Lagrangian finite elements, the so-called {\em stickiness} phenomenon becomes apparent in the numerical experiments we present.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantitative Stability and Error Estimates for Optimal Transport Plans

Optimal transport maps and plans between two absolutely continuous measures $μ$ and $ν$ can be approximated by solving semi-discrete or fully-discrete optimal transport problems. These two problems ensue from approximating $μ$ or both $μ$ and $ν$ by Dirac measures. Extending an idea from [Gigli, On Hölder continuity-in-time of the optimal transport map towards measures along a curve], we characterize how transport plans change under perturbation of both $μ$ and $ν$. We apply this insight to prove error estimates for semi-discrete and fully-discrete algorithms in terms of errors solely arising from approximating measures. We obtain weighted $L^2$ error estimates for both types of algorithms with a convergence rate $O(h^{1/2})$. This coincides with the rate in [Berman, Convergence rates for discretized Monge--Ampère equations and quantitative stability of Optimal Transport, Theorem 5.4] for semi-discrete methods, but the error notion is different.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A diffuse interface model for two-phase ferrofluid flows

We develop a model describing the behavior of two-phase ferrofluid flows using phase field-techniques and present an energy-stable numerical scheme for it. For a simplified, yet physically realistic, version of this model and the corresponding numerical scheme we prove, in addition to stability, convergence and as by-product existence of solutions. With a series of numerical experiments we illustrate the potential of these simple models and their ability to capture basic phenomenological features of ferrofluids such as the Rosensweig instability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fast multilevel solvers for a class of discrete fourth order parabolic problems

In this paper, we study fast iterative solvers for the solution of fourth order parabolic equations discretized by mixed finite element methods. We propose to use consistent mass matrix in the discretization and use lumped mass matrix to construct efficient preconditioners. We provide eigenvalue analysis for the preconditioned system and estimate the convergence rate of the preconditioned GMRes method. Furthermore, we show that these preconditioners only need to be solved inexactly by optimal multigrid algorithms. Our numerical examples indicate that the proposed preconditioners are very efficient and robust with respect to both discretization parameters and diffusion coefficients. We also investigate the performance of multigrid algorithms with either collective smoothers or distributive smoothers when solving the preconditioner systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High-Order AFEM for the Laplace-Beltrami Operator: Convergence Rates

We present a new AFEM for the Laplace-Beltrami operator with arbitrary polynomial degree on parametric surfaces, which are globally $W^1_\infty$ and piecewise in a suitable Besov class embedded in $C^{1,α}$ with $α\in (0,1]$. The idea is to have the surface sufficiently well resolved in $W^1_\infty$ relative to the current resolution of the PDE in $H^1$. This gives rise to a conditional contraction property of the PDE module. We present a suitable approximation class and discuss its relation to Besov regularity of the surface, solution, and forcing. We prove optimal convergence rates for AFEM which are dictated by the worst decay rate of the surface error in $W^1_\infty$ and PDE error in $H^1$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On p-Robust Saturation for hp-AFEM

We consider the standard adaptive finite element loop SOLVE, ESTIMATE, MARK, REFINE, with ESTIMATE being implemented using the $p$-robust equilibrated flux estimator, and MARK being Dörfler marking. As a refinement strategy we employ $p$-refinement. We investigate the question by which amount the local polynomial degree on any marked patch has to be increase in order to achieve a $p$-independent error reduction. The resulting adaptive method can be turned into an instance optimal $hp$-adaptive method by the addition of a coarsening routine.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimizing the Kelvin force in a moving target subdomain

In order to generate a desired Kelvin (magnetic) force in a target subdomain moving along a prescribed trajectory, we propose a minimization problem with a tracking type cost functional. We use the so-called dipole approximation to realize the magnetic field, where the location and the direction of the magnetic sources are assumed to be fixed. The magnetic field intensity acts as the control and exhibits limiting pointwise constraints. We address two specific problems: the first one corresponds to a fixed final time whereas the second one deals with an unknown force to minimize the final time. We prove existence of solutions and deduce local uniqueness provided that a second order sufficient condition is valid. We use the classical backward Euler scheme for time discretization. For both problems we prove the $H^1$-weak convergence of this semi-discrete numerical scheme. This result is motivated by $Γ$-convergence and does not require second order sufficient condition. If the latter holds then we prove $H^1$-strong local convergence. We report computational results to assess the performance of the numerical methods. As an application, we study the control of magnetic nanoparticles as those used in magnetic drug delivery, where the optimized Kelvin force is used to transport the drug to a desired location.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A PDE Approach to Numerical Fractional Diffusion

Fractional diffusion has become a fundamental tool for the modeling of multiscale and heterogeneous phenomena. However, due to its nonlocal nature, its accurate numerical approximation is delicate. We survey our research program on the design and analysis of efficient solution techniques for problems involving fractional powers of elliptic operators. Starting from a localization PDE result for these operators, we develop local techniques for their solution: a priori and a posteriori error analyses, adaptivity and multilevel methods. We show the flexibility of our approach by proposing and analyzing local solution techniques for a space-time fractional parabolic equation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A PDE approach to space-time fractional parabolic problems

We study solution techniques for parabolic equations with fractional diffusion and Caputo fractional time derivative, the latter being discretized and analyzed in a general Hilbert space setting. The spatial fractional diffusion is realized as the Dirichlet-to-Neumann map for a nonuniformly elliptic problem posed on a semi-infinite cylinder in one more spatial dimension. We write our evolution problem as a quasi-stationary elliptic problem with a dynamic boundary condition. We propose and analyze an implicit fully-discrete scheme: first-degree tensor product finite elements in space and an implicit finite difference discretization in time. We prove stability and error estimates for this scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Adaptive Spectral Galerkin Methods with Dynamic Marking

The convergence and optimality theory of adaptive Galerkin methods is almost exclusively based on the Dörfler marking. This entails a fixed parameter and leads to a contraction constant bounded below away from zero. For spectral Galerkin methods this is a severe limitation which affects performance. We present a dynamic marking strategy that allows for a super-linear relation between consecutive discretization errors, and show exponential convergence with linear computational complexity whenever the solution belongs to a Gevrey approximation class.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bilayer Plates: Model Reduction, $Γ$-Convergent Finite Element Approximation and Discrete Gradient Flow

The bending of bilayer plates is a mechanism which allows for large deformations via small externally induced lattice mismatches of the underlying materials. Its mathematical modeling, discussed herein, consists of a nonlinear fourth order problem with a pointwise isometry constraint. A discretization based on Kirchhoff quadrilaterals is devised and its $Γ$-convergence is proved. An iterative method that decreases the energy is proposed and its convergence to stationary configurations is investigated. Its performance, as well as reduced model capabilities, are explored via several insightful numerical experiments involving large (geometrically nonlinear) deformations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convergence and Optimality of hp-AFEM

We design and analyze an adaptive $hp$-finite element method (hp-AFEM) in dimensions $n=1,2$. The algorithm consists of iterating two routines: hp-NEARBEST finds a near-best $hp$-approximation of the current discrete solution and data to a desired accuracy, and REDUCE improves the discrete solution to a finer but comparable accuracy. The former hinges on a recent algorithm by Binev for adaptive $hp$-approximation, and acts as a coarsening step. We prove convergence and instance optimality.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal Control of a Free Boundary Problem with Surface Tension Effects: A Priori Error Analysis

We present a finite element method along with its analysis for the optimal control of a model free boundary problem with surface tension effects, formulated and studied in \cite{HAntil_RHNochetto_PSodre_2014a}. The state system couples the Laplace equation in the bulk with the Young-Laplace equation on the free boundary to account for surface tension. We first prove that the state and adjoint system have the requisite regularity for the error analysis (strong solutions). We discretize the state, adjoint and control variables via piecewise linear finite elements and show optimal $O(h)$ error estimates for all variables, including the control. This entails using the second order sufficient optimality conditions of \cite{HAntil_RHNochetto_PSodre_2014a}, and the first order necessary optimality conditions for both the continuous and discrete systems. We conclude with two numerical examples which examine the various error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The equations of ferrohydrodynamics: modeling and numerical methods

We discuss the equations describing the motion of ferrofluids subject to an external magnetic field. We concentrate on the model proposed by R. Rosensweig, provide an appropriate definition for the effective magnetizing field, and explain the simplifications behind this definition. We show that this system is formally energy stable, and devise a numerical scheme that mimics the same stability estimate. We prove that solutions of the numerical scheme always exist and, under further simplifying assumptions, that the discrete solutions converge. We also discuss alternative formulations proposed in pre-existing work, primarily involving a regularization of the magnetization equation and supply boundary conditions which lead to an energy stable system. We present a series of numerical experiments which illustrate the potential of the scheme in the context of real applications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Stokes problem with Navier slip boundary condition: Minimal fractional Sobolev regularity of the domain

We prove well-posedness in reflexive Sobolev spaces of weak solutions to the stationary Stokes problem with Navier slip boundary condition over bounded domains $Ω$ of $\mathbb{R}^n$ of class $W^{2-1/s}_s$, $s>n$. Since such domains are of class $C^{1,1-n/s}$, our result improves upon the recent one by Amrouche-Seloula, who assume $Ω$ to be of class $C^{1,1}$. We deal with the slip boundary condition directly via a new localization technique, which features domain, space and operator decompositions. To flatten the boundary of $Ω$ locally, we construct a novel $W^2_s$ diffeomorphism for $Ω$ of class $W^{2-1/s}_s$. The fractional regularity gain, from $2-1/s$ to $2$, guarantees that the Piola transform is of class $W^1_s$. This allows us to transform $W^1_r$ vector fields without changing their regularity, provided $r\le s$, and preserve the unit normal which is Hölder. It is in this sense that the boundary regularity $W^{2-1/s}_s$ seems to be minimal.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A PDE approach to fractional diffusion: a posteriori error analysis

We derive a computable a posteriori error estimator for the $α$-harmonic extension problem, which localizes the fractional powers of elliptic operators supplemented with Dirichlet boundary conditions. Our a posteriori error estimator relies on the solution of small discrete problems on anisotropic cylindrical stars. It exhibits built-in flux equilibration and is equivalent to the energy error up to data oscillation, under suitable assumptions. We design a simple adaptive algorithm and present numerical experiments which reveal a competitive performance.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence rates for the classical, thin and fractional elliptic obstacle problems

We review the finite element approximation of the classical obstacle problem in energy and max-norms and derive error estimates for both the solution and the free boundary. On the basis of recent regularity results we present an optimal error analysis for the thin obstacle problem. Finally, we discuss the localization of the obstacle problem for the fractional Laplacian and prove quasi-optimal convergence rates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multilevel methods for nonuniformly elliptic operators

We develop and analyze multilevel methods for nonuniformly elliptic operators whose ellipticity holds in a weighted Sobolev space with an $A_2$--Muckenhoupt weight. Using the so-called Xu-Zikatanov (XZ) identity, we derive a nearly uniform convergence result, under the assumption that the underlying mesh is quasi-uniform. We also consider the so-called $α$-harmonic extension to localize fractional powers of elliptic operators. Motivated by the scheme proposed in [R.H. Nochetto, E. Otarola, and A.J. Salgado. A PDE approach to fractional diffusion in general domains: a priori error analysis. arXiv:1302.0698, 2013] we present a multilevel method with line smoothers and obtain a nearly uniform convergence result on anisotropic meshes. Numerical experiments reveal a competitive performance of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimal Control of a Free Boundary Problem: Analysis with Second Order Sufficient Conditions

We consider a PDE-constrained optimization problem governed by a free boundary problem. The state system is based on coupling the Laplace equation in the bulk with a Young-Laplace equation on the free boundary to account for surface tension, as proposed by P.\ Saavedra and L.\ R.\ Scott \cite{PSaavedra_RScott_1991}. This amounts to solving a second order system both in the bulk and on the interface. Our analysis hinges on a convex control constraint such that the state constraints are always satisfied. Using only first order regularity we show that the control to state operator is twice continuously Fréchet differentiable. We improve slightly the regularity of the state variables and exploit it to show existence of a control together with second order sufficient optimality conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Piecewise polynomial interpolation in Muckenhoupt weighted Sobolev spaces and applications

We develop a constructive piecewise polynomial approximation theory in weighted Sobolev spaces with Muckenhoupt weights for any polynomial degree. The main ingredients to derive optimal error estimates for an averaged Taylor polynomial are a suitable weighted Poincare inequality, a cancellation property and a simple induction argument. We also construct a quasi-interpolation operator, built on local averages over stars, which is well defined for functions in $L^1$. We derive optimal error estimates for any polynomial degree on simplicial shape regular meshes. On rectangular meshes, these estimates are valid under the condition that neighboring elements have comparable size, which yields optimal anisotropic error estimates over $n$-rectangular domains. The interpolation theory extends to cases when the error and function regularity require different weights. We conclude with three applications: nonuniform elliptic boundary value problems, elliptic problems with singular sources, and fractional powers of elliptic operators.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A PDE approach to fractional diffusion in general domains: a priori error analysis

The purpose of this work is the study of solution techniques for problems involving fractional powers of symmetric coercive elliptic operators in a bounded domain with Dirichlet boundary conditions. These operators can be realized as the Dirichlet to Neumann map for a degenerate/singular elliptic problem posed on a semi-infinite cylinder, which we analyze in the framework of weighted Sobolev spaces. Motivated by the rapid decay of the solution of this problem, we propose a truncation that is suitable for numerical approximation. We discretize this truncation using first degree tensor product finite elements. We derive a priori error estimates in weighted Sobolev spaces. The estimates exhibit optimal regularity but suboptimal order for quasi-uniform meshes. For anisotropic meshes, instead, they are quasi-optimal in both order and regularity. We present numerical experiments to illustrate the method's performance.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Adaptive Finite Element Methods for Elliptic Problems with Discontinuous Coefficients

Elliptic partial differential equations (PDEs) with discontinuous diffusion coefficients occur in application domains such as diffusions through porous media, electro-magnetic field propagation on heterogeneous media, and diffusion processes on rough surfaces. The standard approach to numerically treating such problems using finite element methods is to assume that the discontinuities lie on the boundaries of the cells in the initial triangulation. However, this does not match applications where discontinuities occur on curves, surfaces, or manifolds, and could even be unknown beforehand. One of the obstacles to treating such discontinuity problems is that the usual perturbation theory for elliptic PDEs assumes bounds for the distortion of the coefficients in the $L_\infty$ norm and this in turn requires that the discontinuities are matched exactly when the coefficients are approximated. We present a new approach based on distortion of the coefficients in an $L_q$ norm with $q<\infty$ which therefore does not require the exact matching of the discontinuities. We then use this new distortion theory to formulate new adaptive finite element methods (AFEMs) for such discontinuity problems. We show that such AFEMs are optimal in the sense of distortion versus number of computations, and report insightful numerical results supporting our analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The micropolar Navier-Stokes equations: A priori error analysis

The unsteady Micropolar Navier-Stokes Equations (MNSE) are a system of parabolic partial differential equations coupling linear velocity and pressure with angular velocity: material particles have both translational and rotational degrees of freedom. We propose and analyze a first order semi-implicit fully-discrete scheme for the MNSE, which decouples the computation of the linear and angular velocities, is unconditionally stable and delivers optimal convergence rates under assumptions analogous to those used for the Navier-Stokes equations. With the help of our scheme we explore some qualitative properties of the MNSE related to ferrofluid manipulation and pumping. Finally, we propose a second order scheme and show that it is almost unconditionally stable.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Two-step greedy algorithm for reduced order quadratures

We present an algorithm to generate application-specific, global reduced order quadratures (ROQ) for multiple fast evaluations of weighted inner products between parameterized functions. If a reduced basis (RB) or any other projection-based model reduction technique is applied, the dimensionality of integrands is reduced dramatically; however, the cost of approximating the integrands by projection still scales as the size of the original problem. In contrast, using discrete empirical interpolation (DEIM) points as ROQ nodes leads to a computational cost which depends linearly on the dimension of the reduced space. Generation of a reduced basis via a greedy procedure requires a training set, which for products of functions can be very large. Since this direct approach can be impractical in many applications, we propose instead a two-step greedy targeted towards approximation of such products. We present numerical experiments demonstrating the accuracy and the efficiency of the two-step approach. The presented ROQ are expected to display very fast convergence whenever there is regularity with respect to parameter variation. We find that for the particular application here considered, one driven by gravitational wave physics, the two-step approach speeds up the offline computations to build the ROQ by more than two orders of magnitude. Furthermore, the resulting ROQ rule is found to converge exponentially with the number of nodes, and a factor of ~50 savings, without loss of accuracy, is observed in evaluations of inner products when ROQ are used as a downsampling strategy for equidistant samples using the trapezoidal rule. While the primary focus of this paper is on quadrature rules for inner products of parameterized functions, our method can be easily adapted to integrations of single parameterized functions, and some examples of this type are considered.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A posteriori error estimates for the Electric Field Integral Equation on polyhedra

We present a residual-based a posteriori error estimate for the Electric Field Integral Equation (EFIE) on a bounded polyhedron. The EFIE is a variational equation formulated in a negative order Sobolev space on the surface of the polyhedron. We express the estimate in terms of square-integrable and thus computable quantities and derive global lower and upper bounds (up to oscillation terms).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Total Variation Diminishing Interpolation Operator and Applications

We construct, on continuous $Q_1$ finite elements over Cartesian meshes, an interpolation operator that does not increase the total variation. The operator is stable in $L^1$ and exhibits second order approximation properties. With the help of it we provide improved error estimates for discrete minimizers of the total variation denoising problem and for total variation flows.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Adaptive Fourier-Galerkin Methods

We study the performance of adaptive Fourier-Galerkin methods in a periodic box in $\mathbb{R}^d$ with dimension $d\ge 1$. These methods offer unlimited approximation power only restricted by solution and data regularity. They are of intrinsic interest but are also a first step towards understanding adaptivity for the $hp$-FEM. We examine two nonlinear approximation classes, one classical corresponding to algebraic decay of Fourier coefficients and another associated with exponential decay. We study the sparsity classes of the residual and show that they are the same as the solution for the algebraic class but not for the exponential one. This possible sparsity degradation for the exponential class can be compensated with coarsening, which we discuss in detail. We present several adaptive Fourier algorithms, and prove their contraction and optimal cardinality properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Contraction and optimality properties of adaptive Legendre-Galerkin methods: the 1-dimensional case

As a first step towards a mathematically rigorous understanding of adaptive spectral/$hp$ discretizations of elliptic boundary-value problems, we study the performance of adaptive Legendre-Galerkin methods in one space dimension. These methods offer unlimited approximation power only restricted by solution and data regularity. Our investigation is inspired by a similar study that we recently carried out for Fourier-Galerkin methods in a periodic box. We first consider an "ideal" algorithm, which we prove to be convergent at a fixed rate. Next we enhance its performance, consistently with the expected fast error decay of high-order methods, by activating a larger set of degrees of freedom at each iteration. We guarantee optimality (in the non-linear approximation sense) by incorporating a coarsening step. Optimality is measured in terms of certain sparsity classes of the Gevrey type, which describe a (sub-)exponential decay of the best approximation error.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Discrete Total Variation Flows Without Regularization

We propose and analyze an algorithm for the solution of the $L^2$-subgradient flow of the total variation functional. The algorithm involves no regularization, thus the numerical solution preserves the main features that motivate practitioners to consider this type of energy. We propose an iterative scheme for the solution of the arising problems, show that the iterations converge, and develop a stopping criterion for them. We present numerical experiments which illustrate the power of the method, explore the solution behavior, and compare with regularized flows.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Diffuse Interface Model for Electrowetting with Moving Contact Lines

We introduce a diffuse interface model for the phenomenon of electrowetting on dielectric and present an analysis of the arising system of equations. Moreover, we study discretization techniques for the problem. The model takes into account different material parameters on each phase and incorporates the most important physical processes, such as incompressibility, electrostatics and dynamic contact lines; necessary to properly reflect the relevant phenomena. The arising nonlinear system couples the variable density incompressible Navier-Stokes equations for velocity and pressure with a Cahn-Hilliard type equation for the phase variable and chemical potential, a convection diffusion equation for the electric charges and a Poisson equation for the electric potential. Numerical experiments are presented, which illustrate the wide range of effects the model is able to capture, such as splitting and coalescence of droplets.

0 citations0 reviewsNo code linkedOpen access