A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes. Marco Verani
|
|
|
- Bryce Bailey
- 5 years ago
- Views:
Transcription
1 A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes Marco Verani MOX, Department of Mathematics, Politecnico di Milano Joint work with: P. F. Antonietti (MOX Politecnico di Milano, Italy) L. Beirão da Veiga (Università di Milano, Italy) S. Scacchi (Università di Milano, Italy) Polytopal Element Methods, Georgia Tech, October 26th, 2015 M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
2 Outline 1 Cahn-Hilliard equation 2 C 1 -VEM for Cahn-Hilliard 3 Numerical results 4 Conclusions M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
3 Cahn-Hilliard (CH) problem Let Ω R 2 be an open bounded domain. Let ψ(x) = (1 x 2 ) 2 /4 and φ(x) = ψ (x). t u ( φ(u) γ 2 u(t) ) = 0 in Ω [0, T ], u(, 0) = u 0 ( ) in Ω, ( n u = n φ(u) γ 2 u(t) ) = 0 on Ω [0, T ], where n denotes the (outward) normal derivative and γ R +, 0 < γ 1, represents the interface parameter. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
4 Cahn-Hilliard (CH) problem Let Ω R 2 be an open bounded domain. Let ψ(x) = (1 x 2 ) 2 /4 and φ(x) = ψ (x). t u ( φ(u) γ 2 u(t) ) = 0 in Ω [0, T ], u(, 0) = u 0 ( ) in Ω, ( n u = n φ(u) γ 2 u(t) ) = 0 on Ω [0, T ], where n denotes the (outward) normal derivative and γ R +, 0 < γ 1, represents the interface parameter. CH models the evolution of interfaces specified as the level set of a smooth continuos function u exhibiting large gradients across the interface [Korteweg 1901, Cahn and Hilliard 1958, Ginzburg and Landau 1965, van der Waals 1979] Many applications: phase separation in binary alloys, tumor growth, origin of Saturn s rings, separation of di-block copolymers, population dynamics, image processing, clustering of mussels... M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
5 Huge literature on numerical methods for CH... M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
6 Huge literature on numerical methods for CH... Here we limit ourvselves to FEM as the main properties (and limitations) of these schemes are instrumental to motivate the introduction of our new approach. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
7 Huge literature on numerical methods for CH... Here we limit ourvselves to FEM as the main properties (and limitations) of these schemes are instrumental to motivate the introduction of our new approach. CH is a fourth order nonlinear problem C 1 FEM (see [Elliott Zheng, 1986], [Elliott French, 1987]). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
8 Huge literature on numerical methods for CH... Here we limit ourvselves to FEM as the main properties (and limitations) of these schemes are instrumental to motivate the introduction of our new approach. CH is a fourth order nonlinear problem C 1 FEM (see [Elliott Zheng, 1986], [Elliott French, 1987]). To avoid the well known difficulty met in the implementation of C 1 FEM, another possibility is the use of non-conforming (see [Elliott French, 1989] or discontinuous (see [Wells Kuhl Garikipati, 2006]) methods. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
9 Huge literature on numerical methods for CH... Here we limit ourvselves to FEM as the main properties (and limitations) of these schemes are instrumental to motivate the introduction of our new approach. CH is a fourth order nonlinear problem C 1 FEM (see [Elliott Zheng, 1986], [Elliott French, 1987]). To avoid the well known difficulty met in the implementation of C 1 FEM, another possibility is the use of non-conforming (see [Elliott French, 1989] or discontinuous (see [Wells Kuhl Garikipati, 2006]) methods. Alternatively, use mixed methods (see e.g. [Elliot French Milner, 1989], [Elliott Larsson, 1992] and [Kay Styles Suli, 2009] for the continuous and discontinuous setting, respectively). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
10 Huge literature on numerical methods for CH... Here we limit ourvselves to FEM as the main properties (and limitations) of these schemes are instrumental to motivate the introduction of our new approach. CH is a fourth order nonlinear problem C 1 FEM (see [Elliott Zheng, 1986], [Elliott French, 1987]). To avoid the well known difficulty met in the implementation of C 1 FEM, another possibility is the use of non-conforming (see [Elliott French, 1989] or discontinuous (see [Wells Kuhl Garikipati, 2006]) methods. Alternatively, use mixed methods (see e.g. [Elliot French Milner, 1989], [Elliott Larsson, 1992] and [Kay Styles Suli, 2009] for the continuous and discontinuous setting, respectively). Recently, isogeometric analysis [Gómez Calo Bazilevs Hughes, 2008]. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
11 Huge literature on numerical methods for CH... Here we limit ourvselves to FEM as the main properties (and limitations) of these schemes are instrumental to motivate the introduction of our new approach. CH is a fourth order nonlinear problem C 1 FEM (see [Elliott Zheng, 1986], [Elliott French, 1987]). To avoid the well known difficulty met in the implementation of C 1 FEM, another possibility is the use of non-conforming (see [Elliott French, 1989] or discontinuous (see [Wells Kuhl Garikipati, 2006]) methods. Alternatively, use mixed methods (see e.g. [Elliot French Milner, 1989], [Elliott Larsson, 1992] and [Kay Styles Suli, 2009] for the continuous and discontinuous setting, respectively). Recently, isogeometric analysis [Gómez Calo Bazilevs Hughes, 2008]. Here we propose a C 1 -VEM scheme (see, e.g., [Beirão da Veiga Brezzi Cangiani Manzini Marini Russo, 2013] for an intro to VEM for second order pbs and [Brezzi Marini, 2013] for plate bending pb.) M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
12 Weak formulation of CH a (v, w) = ( 2 v) : ( 2 w) dx v, w H 2 (Ω), Ω a (v, w) = ( v) ( w) dx v, w H 1 (Ω), Ω a 0 (v, w) = v w dx v, w L 2 (Ω), Ω M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
13 Weak formulation of CH a (v, w) = ( 2 v) : ( 2 w) dx v, w H 2 (Ω), Ω a (v, w) = ( v) ( w) dx v, w H 1 (Ω), a 0 (v, w) = Ω v w dx v, w L 2 (Ω), Ω and the semi-linear form r(z; v, w) = φ (z) v w dx z, v, w H 2 (Ω). Ω M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
14 Weak formulation of CH a (v, w) = ( 2 v) : ( 2 w) dx v, w H 2 (Ω), Ω a (v, w) = ( v) ( w) dx v, w H 1 (Ω), Ω a 0 (v, w) = v w dx v, w L 2 (Ω), and the semi-linear form r(z; v, w) = φ (z) v w dx Let Ω Ω z, v, w H 2 (Ω). V = { v H 2 (Ω) : n u = 0 on Ω } The weak formulation reads as: find u(, t) V such that { a 0 ( t u, v) + γ 2 a (u, v) + r(u; u, v) = 0 v V, u(, 0) = u 0 ( ). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
15 Towards C 1 -VEM discretization of CH Ω h decomposition of Ω into polygons E M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
16 Towards C 1 -VEM discretization of CH Ω h decomposition of Ω into polygons E Intermediate local space (see [Brezzi Marini, 2013]): { Ṽ h E = v H 2 (E) : 2 v P 2 (E), v E C 0 ( E), v e P 3 (e) e E, } v E [C 0 ( E)] 2, n v e P 1 (e) e E, M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
17 Towards C 1 -VEM discretization of CH Ω h decomposition of Ω into polygons E Intermediate local space (see [Brezzi Marini, 2013]): { Ṽ h E = v H 2 (E) : 2 v P 2 (E), v E C 0 ( E), v e P 3 (e) e E, } v E [C 0 ( E)] 2, n v e P 1 (e) e E, Linear Operators: D1 evaluation of v at the vertexes of E; D2 evaluation of v at the vertexes of E. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
18 The operator Π E The projection operator Π E : Ṽh E P 2 (E) is defined by { a E (Π E v h, q) = a E (v h, q) q P 2 (E) for all v h Ṽh E where ((v h, w h )) E = ((Π E v h, q)) E = ((v h, q)) E ν vertexes of E v h (ν) w h (ν) q P 1 (E), v h, w h C 0 (E). The operator Π E is uniquely determined by D1 and D2, i.e. Π E v h depends only on the values of v h and v h at the vertices of E. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
19 Virtual local space W h E = { v Ṽh E : Π E (v h) q dx = E E v h q dx q P 2 (E) }. [Ahmad Alsaedi Brezzi Marini Russo, 2013] M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
20 Virtual local space W h E = { v Ṽh E : Π E (v h) q dx = E E v h q dx q P 2 (E) }. Properties: The set of operators D1 and D2 constitutes a set of degrees of freedom for the space W h E. [Ahmad Alsaedi Brezzi Marini Russo, 2013] M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
21 Virtual local space W h E = { v Ṽh E : Π E (v h) q dx = E E v h q dx q P 2 (E) }. Properties: The set of operators D1 and D2 constitutes a set of degrees of freedom for the space W h E. P 2 (E) W h E. [Ahmad Alsaedi Brezzi Marini Russo, 2013] M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
22 Virtual local space W h E = { v Ṽh E : Π E (v h) q dx = E E v h q dx q P 2 (E) }. Properties: The set of operators D1 and D2 constitutes a set of degrees of freedom for the space W h E. P 2 (E) W h E. The L 2 -projection Π 0 E : W h E P 2 (E) defined by a 0 E (Π0 E v h, q) = a 0 E (v h, q) q P 2 (E), is computable (only) on the basis of the values of the degrees of freedom D1 and D2. [Ahmad Alsaedi Brezzi Marini Russo, 2013] M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
23 Global VEM space Properties: W h = { v V : v E W h E E Ω h } W h H 2 (Ω) conforming VEM solution. The global degrees of freedom are: values of v h at the vertexes of the mesh Ω h ; values of v h at the vertexes of the mesh Ω h. Dimension of W h is three times the number of vertexes in the mesh. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
24 The operator Π E We define Π E : W h E P 2 (E) by ae (Π E v h, q) = ae (v h, q) Π E v h dx = v h dx. E E q P 2 (E) Since the bilinear form ae (, ) has a non trivial kernel (given by the constant functions) we added a second condition in order to keep the operator Π E well defined. The operator Π E is uniquely determined by D1 and D2 M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
25 Local discrete bilinear forms We propose the following discrete (and symmetric) local forms ah,e (v h, w h ) = ae (Π E v h, Π E w h) + h 2 E s E (v h Π E v h, w h Π E w h), ah,e (v h, w h ) = ae (Π E v h, Π E w h) + s E (v h Π E v h, w h Π E w h), a 0 h,e (v h, w h ) = a 0 E (Π0 E v h, Π 0 E w h) + h 2 E s E (v h Π 0 E v h, w h Π 0 E w h), for all v h, w h W h E, where s E (v h, w h ) = ν vertexes of E ( ) v h (ν) w h (ν) + (h ν ) 2 v h (ν) w h (ν) being h ν some characteristic mesh size lenght associated to the node ν (for instance the maximum diameter among the elements having ν as a vertex). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
26 Properties of the discrete bilinear forms The symbol stands for the symbol, or 0. (Consistency) a h,e (p, v h) = a E (p, v h) p P 2 (E), v h W h E, (Stability) 1 There exist two positive constants c, c independent of the element E {Ω h } h such that c a E (v h, v h ) a h,e (v h, v h ) c a E (v h, v h ) v h W h E, The bilinear forms a h,e (, ) are continuous with respect to the relevant norm: H 2 (for = ), H 1 (for = ) and L 2 (for = 0). 1 It holds under regularity assumption on the polygonal partition M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
27 Global discrete bilinear forms The global discrete bilinear form is assembled as usual a h (v h, w h ) = E Ω h a h,e (v h, w h ) v h, w h W h, with the usual multiple meaning of the symbol. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
28 Discrete semilinear form Continuous semilinear form r(z; v, w) = r E (z; v, w) z, v, w H 2 (Ω), E Ω h r E (z; v, w) = (3z(x) 2 1) v(x) w(x) dx E Ω h. E M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
29 Discrete semilinear form Continuous semilinear form r(z; v, w) = r E (z; v, w) z, v, w H 2 (Ω), E Ω h r E (z; v, w) = (3z(x) 2 1) v(x) w(x) dx E Ω h. E Idea: on each element E, we approximate the term z(x) 2 with its average M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
30 Discrete semilinear form Continuous semilinear form r(z; v, w) = r E (z; v, w) z, v, w H 2 (Ω), E Ω h r E (z; v, w) = (3z(x) 2 1) v(x) w(x) dx E Ω h. E Idea: on each element E, we approximate the term z(x) 2 with its average Discrete local semilinear form r h,e (z h ; v h, w h ) = φ (z h ) E a h,e (v h, w h ) z h, v h, w h W h E where φ (z h ) E =3 E 1 a 0 h,e (z h, z h ) 1 M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
31 Discrete semilinear form Continuous semilinear form r(z; v, w) = r E (z; v, w) z, v, w H 2 (Ω), E Ω h r E (z; v, w) = (3z(x) 2 1) v(x) w(x) dx E Ω h. E Idea: on each element E, we approximate the term z(x) 2 with its average Discrete local semilinear form r h,e (z h ; v h, w h ) = φ (z h ) E a h,e (v h, w h ) z h, v h, w h W h E where φ (z h ) E =3 E 1 a 0 h,e (z h, z h ) 1 The global form is then assembled as usual r h (z h ; v h, w h ) = E Ω h r h,e (z h ; v h, w h ) w h, r h, v h W h. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
32 C 1 -VEM semi-discrete scheme Discrete global space with boundary conditions: W 0 h = W h V = { v W h : n u = 0 on Ω }. The semi-discrete problem is: Find u h (, t) in Wh 0 such that ah 0 ( tu h, v h ) + γ 2 ah (u h, v h ) + r h (u h, u h ; v h ) = 0, u h (0, ) = u 0,h ( ), for all v h W 0 h being u 0,h W 0 h a suitable approximation of u 0. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
33 Convergence result for the semi-discrete scheme Theorem (Antonietti, Beirao da Veiga, Sacchi, V., 2015) Let u be the (sufficiently regular) exact solution of the Cahn-Hilliard problem and u h the solution of the semi-discrete VEM problem. Then for all t [0, T ] it holds u u h L 2 (Ω) h 2. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
34 Numerical results Exact solution: u(x, y, t) = t cos(2πx) cos(2πy) on Ω = (0, 1) 2 Table: Test 1: H 2, H 1 and L 2 errors and convergence rates α computed on four quadrilateral meshes discretizing the unit square. h u h u H 2 (Ω) α u h u H 1 (Ω) α u h u L 2 (Ω) α 1/ e e e-2 1/ e e e / e e e / e e e Parameters: γ = 1/10 and t = 1e 7. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
35 Evolution of an ellipse (I) (a) Quadrilateral mesh of = elements (49923 dofs). Figure: Snapshots at three temporal frames (t = 0, 0.5, 1). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
36 Evolution of an ellipse (II) (a) Triangular mesh of 8576 elements (13167 dofs). Figure: Snapshots at three temporal frames (t = 0, 0.5, 1). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
37 Voronoi polygonal meshes Figure: Examples of Voronoi polygonal meshes (quadrilaterals, pentagons, hexagons) with 10 (left) and 100 (right) elements. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
38 Evolution of a cross (I) (a) Quadrilateral mesh of = elements (49923 dofs). Figure: Snapshots at three temporal frames (t = 0, 0.05, 1). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
39 Evolution of a cross (II) (a) Triangular mesh of 8576 elements (13167 dofs). Figure: Snapshots at three temporal frames (t = 0, 0.05, 1). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
40 Evolution of a cross (III) (a) Voronoi polygonal mesh of elements (59490 dofs). Figure: Snapshots at three temporal frames (t = 0, 0.05, 1). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
41 Spinoidal decomposition (a) Voronoi polygonal mesh of elements (59490 dofs). Figure: Snapshots at three temporal frames (t = 0.01, 0.05, 5). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
42 Conclusions C 1 -VEM of minimal degree for the approximation of CH equation. Theoretical L 2 -convergence of the semi-discrete scheme Numerical tests Reference P. F. Antonietti, L. Beirao da Veiga, S. Scacchi and M. Verani, A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes, accepted on SINUM, M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
43 Spinoidal decomposition (I) (a) Quadrilateral mesh of = elements (49923 dofs). Figure: Snapshots at three temporal frames (t = 0.01, 0.05, 5). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
44 Spinoidal decomposition (II) (a) Triangular mesh of 8576 elements (13167 dofs). Figure: Snapshots at three temporal frames (t = 0.075, 0.25, 1.25). M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
45 VEM vs mixed-dg DG VEM h dof u h u H 1 (Ω) u h u L 2 (Ω) dof u h u H 1 (Ω) u h u L 2 (Ω) 1/ e e e e-4 1/ e e e e-4 1/ e e e e-5 1/ e e e e-5 CH with exact solution, γ = 1/10 and t = 1e 7. The simulation is run for 10 time steps. The H 1 and L 2 errors are computed at the final time step on four triangular (for the DG method) and quadrilateral (for the VEM) meshes discretizing the unit square. M. Verani (MOX - PoliMi) VEM Cahn-Hilliard POEMs / 28
The Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
arxiv: v1 [math.na] 11 Feb 2015
![arxiv: v1 [math.na] 11 Feb 2015](/thumbs/74/70570447.jpg "arxiv: v1 [math.na] 11 Feb 2015")
A C 1 VIRTUAL LMNT MTHOD FOR TH CAHN-HILLIARD QUATION WITH POLYGONAL MSHS PAOLA F. ANTONITTI, LOURNCO BIRÃO DA VIGA, SIMON SCACCHI, AND MARCO VRANI arxiv:1502.03259v1 [math.na] 11 Feb 2015 Abstract. In
Virtual Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
arxiv: v1 [math.na] 27 Jan 2016
![arxiv: v1 [math.na] 27 Jan 2016](/thumbs/92/109630157.jpg "arxiv: v1 [math.na] 27 Jan 2016")
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
A $C^1$ virtual element method for the Cahn-Hilliard equation with polygonal meshes
MOX-Report No. 12/2015 A $C^1$ virtual element method for the Cahn-Hilliard equation with polygonal meshes Antonietti, P. F.; Beirao da Veiga, L.; Scacchi, S.; Verani, M. MOX, Dipartimento di Matematica
arxiv: v1 [math.na] 27 Jan 2016
![arxiv: v1 [math.na] 27 Jan 2016](/thumbs/95/124616312.jpg "arxiv: v1 [math.na] 27 Jan 2016")
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES On the approximation of a virtual coarse space for domain decomposition methods in two dimensions Juan Gabriel Calvo Preprint No. 31-2017 PRAHA 2017
Divergence free Virtual Elements for the Stokes problem on polygonal meshes
Divergence free Virtual Elements for the Stokes problem on polygonal meshes L. Beirão da Veiga, C. Lovadina and G. Vacca Abstract arxiv:1510.01655v1 [math.na] 6 Oct 2015 In the present paper we develop
Finite Element Multigrid Framework for Mimetic Finite Difference Discretizations
Finite Element Multigrid Framework for Mimetic Finite ifference iscretizations Xiaozhe Hu Tufts University Polytopal Element Methods in Mathematics and Engineering, October 26-28, 2015 Joint work with:
arxiv: v1 [math.na] 29 Feb 2016
![arxiv: v1 [math.na] 29 Feb 2016](/thumbs/92/110614295.jpg "arxiv: v1 [math.na] 29 Feb 2016")
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
Basic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
H(div) and H(curl)-conforming Virtual Element Methods
Noname manuscript No. (will be inserted by the editor) H(div) and H(curl)-conforming Virtual Element Methods L. Beirão da Veiga F. Brezzi L. D. Marini A. Russo the date of receipt and acceptance should
Numerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions Juan Gabriel Calvo Preprint No. 25-2017 PRAHA 2017 A VIRTUAL
VEM s for the numerical solution of PDE s. F. Brezzi
VEM s for the numerical solution of PDE s F. Brezzi IUSS-Istituto Universitario di Studi Superiori, Pavia, Italy IMATI-C.N.R., Pavia, Italy From Joint works with L. Beirão da Veiga, L.D. Marini, and A.
THE NONCONFORMING VIRTUAL ELEMENT METHOD
THE NONCONFORMING VIRTUAL ELEMENT METHOD BLANCA AYUSO DE DIOS, ONSTANTIN LIPNIOV, AND GIANMARCO MANZINI Abstract. We introduce the nonconforming Virtual Element Method (VEM) for the approximation of second
Serendipity Finite Element Methods for Cardiac Electrophysiology Models
Serendipity Finite Element Methods for Cardiac Electrophysiology Models Andrew Gillette Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/ Andrew Gillette -
Dual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations
Dual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations Christoph Hofer and Ulrich Langer Doctoral Program Computational Mathematics Numerical
Discontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
arxiv: v1 [math.na] 3 Jun 2018
![arxiv: v1 [math.na] 3 Jun 2018](/thumbs/94/120852596.jpg "arxiv: v1 [math.na] 3 Jun 2018")
SUPG stabilization for the nonconforming virtual element method for advection-diffusion-reaction equations S Berrone a, A Borio a, G Manzini b,c arxiv:186879v1 [mathna] 3 Jun 18 Abstract a Dipartimento
Virtual Elements for plates and shells
Virtual Elements for plates and shells Claudia Chinosi Dipartimento di Scienze e Innovazione Tecnologica, Università del Piemonte Orientale, Alessandria, Italia E-mail: claudia.chinosi@unipmn.it eywords:
arxiv: v1 [math.na] 11 Sep 2017
![arxiv: v1 [math.na] 11 Sep 2017](/thumbs/73/68884687.jpg "arxiv: v1 [math.na] 11 Sep 2017")
Dual virtual element method in presence of an inclusion Alessio Fumagalli arxiv:1709.03519v1 [math.na] 11 Sep 2017 Abstract We consider a Darcy problem for saturated porous media written in dual formulation
Isogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
Discontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
On Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
Weak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
Mimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS
THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS Alain Miranville Université de Poitiers, France Collaborators : L. Cherfils, G. Gilardi, G.R. Goldstein, G. Schimperna,
Yongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
Geometric Decomposition of Serendipity Finite Element Spaces
Geometric Decomposition of Serendipity Finite Element Spaces Andrew Gillette Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/ Joint Mathematics Meetings 2013
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
February 18, :41 WSPC/INSTRUCTION FILE Brezzi-Corrections VIRTUAL ELEMENT METHOD FOR GENERAL SECOND ORDER ELLIPTIC PROBLEMS ON POLYGONAL MESHES
February 8, 206 7:4 WSPC/INSTRUCTION FIL Brezzi-Corrections Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company VIRTUAL LMNT MTHOD FOR GNRAL SCOND ORDR LLIPTIC PROBLMS
Hybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
On atomistic-to-continuum couplings without ghost forces
On atomistic-to-continuum couplings without ghost forces Dimitrios Mitsoudis ACMAC Archimedes Center for Modeling, Analysis & Computation Department of Applied Mathematics, University of Crete & Institute
AMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
Applied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
A Virtual Element Method for the Transmission Eigenvalue Problem
A Virtual Element Method for the Transmission Eigenvalue Problem David Mora a,b, Iván Velásquez b,c a GIMNAP, Departamento de Matemática, Universidad del Bío-Bío, Casilla 5-C, Concepción, Chile. arxiv:183.1979v1
Mesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain
Mesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain Stephen Edward Moore Johann Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences,
arxiv: v1 [math.na] 8 Jan 2019
![arxiv: v1 [math.na] 8 Jan 2019](/thumbs/88/117564040.jpg "arxiv: v1 [math.na] 8 Jan 2019")
arallel solvers for virtual element discretizations of elliptic equations in mixed form F. Dassi a, S. Scacchi b a Dipartimento di Matematica, Università degli Studi di Milano-Bicocca, Via Cozzi, 2133
Benchmarking high order finite element approximations for one-dimensional boundary layer problems
Benchmarking high order finite element approximations for one-dimensional boundary layer problems Marcello Malagù,, Elena Benvenuti, Angelo Simone Department of Engineering, University of Ferrara, Italy
An introduction to Isogeometric Analysis
An introduction to Isogeometric Analysis A. Buffa 1, G. Sangalli 1,2, R. Vázquez 1 1 Istituto di Matematica Applicata e Tecnologie Informatiche - Pavia Consiglio Nazionale delle Ricerche 2 Dipartimento
ABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
The Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
Some applications of discontinuous Galerkin methods in solid mechanics
Some applications of discontinuous Galerkin methods in solid mechanics Adrian Lew, Alex Ten Eyck, and Ramsharan Rangarajan Mechanical Engineering, Stanford University, USA. Email: {lewa,alabute,rram}@stanford.edu
On a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
Preconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
Numerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
![[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,](/thumbs/88/117157910.jpg "[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,")
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
Weak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
Well-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
Isogeometric Analysis of the Navier Stokes Cahn Hilliard equations with application to incompressible two-phase flows
Isogeometric Analysis of the Navier Stokes Cahn Hilliard equations with application to incompressible two-phase flows Babak S. Hosseini 1, Stefan Turek 1, Matthias Möller 2 1 TU Dortmund University, Institute
Relaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
Juan Vicente Gutiérrez Santacreu Rafael Rodríguez Galván. Departamento de Matemática Aplicada I Universidad de Sevilla
Doc-Course: Partial Differential Equations: Analysis, Numerics and Control Research Unit 3: Numerical Methods for PDEs Part I: Finite Element Method: Elliptic and Parabolic Equations Juan Vicente Gutiérrez
Some New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
Small energy regularity for a fractional Ginzburg-Landau system
Small energy regularity for a fractional Ginzburg-Landau system Yannick Sire University Aix-Marseille Work in progress with Vincent Millot (Univ. Paris 7) The fractional Ginzburg-Landau system We are interest
The Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
UC Irvine UC Irvine Previously Published Works
UC Irvine UC Irvine Previously Published Works Title Some error analysis on virtual element methods Permalink https://escholarship.org/uc/item/1ph610x7 Journal CALCOLO, 55(1) ISSN 0008-0624 Authors Chen,
with deterministic and noise terms for a general non-homogeneous Cahn-Hilliard equation Modeling and Asymptotics
12-3-2009 Modeling and Asymptotics for a general non-homogeneous Cahn-Hilliard equation with deterministic and noise terms D.C. Antonopoulou (Joint with G. Karali and G. Kossioris) Department of Applied
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART I: ERROR ANALYSIS UNDER MINIMUM REGULARITIES
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM PART I: ERROR ANALYSIS UNDER MINIMUM REGULARITIES XIAOBING FENG AND ANDREAS PROHL Abstract. In this first part
WEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED
C 0 IPG Method for Biharmonic Eigenvalue Problems
C IPG Method for Biharmonic Eigenvalue Problems Susanne C. Brenner, Peter Monk 2, and Jiguang Sun 3 Department of Mathematics and Center for Computation & Technology, Louisiana State University, Baton
c 2013 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 51, No. 1, pp. 654 675 c 2013 Society for Industrial and Applied Mathematics HIERARCHICAL A POSTERIORI ERROR ESTIMATORS FOR THE MIMETIC DISCRETIZATION OF ELLIPTIC PROBLEMS PAOLA
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
Mixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
An Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
Computer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
Introduction to Finite Element computations
Non Linear Computational Mechanics Athens MP06/2012 Introduction to Finite Element computations Vincent Chiaruttini, Georges Cailletaud vincent.chiaruttini@onera.fr Outline Continuous to discrete problem
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
Local discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
New class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
Well-posedness and asymptotic analysis for a Penrose-Fife type phase-field system
0-0 Well-posedness and asymptotic analysis for a Penrose-Fife type phase-field system Buona positura e analisi asintotica per un sistema di phase field di tipo Penrose-Fife Salò, 3-5 Luglio 2003 Riccarda
The CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
On some nonlinear parabolic equation involving variable exponents
On some nonlinear parabolic equation involving variable exponents Goro Akagi (Kobe University, Japan) Based on a joint work with Giulio Schimperna (Pavia Univ., Italy) Workshop DIMO-2013 Diffuse Interface
On second order sufficient optimality conditions for quasilinear elliptic boundary control problems
On second order sufficient optimality conditions for quasilinear elliptic boundary control problems Vili Dhamo Technische Universität Berlin Joint work with Eduardo Casas Workshop on PDE Constrained Optimization
arxiv: v1 [math.na] 8 Feb 2018
![arxiv: v1 [math.na] 8 Feb 2018](/thumbs/94/120327685.jpg "arxiv: v1 [math.na] 8 Feb 2018")
arxiv:180.094v1 [math.na] 8 Feb 018 The nonconforming virtual element method for eigenvalue problems F. Gardini a, G. Manzini b, and G. Vacca c a Dipartimento di Matematica F. Casorati, Università di Pavia,
WEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
The space of C 1 -smooth isogeometric spline functions on trilinearly parameterized volumetric two-patch domains. Katharina Birner and Mario Kapl
The space of C 1 -smooth isogeometric spline functions on trilinearly parameterized volumetric two-patch domains Katharina Birner and Mario Kapl G+S Report No. 69 March 2018 The space of C 1 -smooth isogeometric
CRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
From nonlocal to local Cahn-Hilliard equation. Stefano Melchionna Helene Ranetbauer Lara Trussardi. Uni Wien (Austria) September 18, 2018
From nonlocal to local Cahn-Hilliard equation Stefano Melchionna Helene Ranetbauer Lara Trussardi Uni Wien (Austria) September 18, 2018 SFB P D ME S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
Numerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
A Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
Discontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
Space-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
BACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
arxiv: v2 [math.na] 29 Nov 2016
![arxiv: v2 [math.na] 29 Nov 2016](/thumbs/91/106508701.jpg "arxiv: v2 [math.na] 29 Nov 2016")
Noname manuscript No. (will be inserted by the editor) Multigrid algorithms for hp-version Interior Penalty Discontinuous Galerkin methods on polygonal and polyhedral meshes P. F. Antonietti P. Houston
Paola F. Antonietti 1, Luca Formaggia 2, Anna Scotti 3, Marco Verani 4 and Nicola Verzotti 5. Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher MIMETIC FINITE DIFFERENCE APPROXIMATION OF FLOWS IN FRACTURED POROUS MEDIA Paola
Isogeometric modeling of Lorentz detuning in linear particle accelerator cavities
Isogeometric modeling of Lorentz detuning in linear particle accelerator cavities Mauro Bonafini 1,2, Marcella Bonazzoli 3, Elena Gaburro 1,2, Chiara Venturini 1,3, instructor: Carlo de Falco 4 1 Department
Optimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
Domain Decomposition Preconditioners for Isogeometric Discretizations
Domain Decomposition Preconditioners for Isogeometric Discretizations Luca F Pavarino, Università di Milano, Italy Lorenzo Beirao da Veiga, Università di Milano, Italy Durkbin Cho, Dongguk University,
A-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements
A-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements Joint work with Lourenco Beirão da Veiga (Milano), Claudia Chinosi (Alessandria) and Carlo Lovadina (Pavia) Previous
Finite Element Methods for Fourth Order Variational Inequalities
Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2013 Finite Element Methods for Fourth Order Variational Inequalities Yi Zhang Louisiana State University and Agricultural
Maximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
Algebraic Coarse Spaces for Overlapping Schwarz Preconditioners
Algebraic Coarse Spaces for Overlapping Schwarz Preconditioners 17 th International Conference on Domain Decomposition Methods St. Wolfgang/Strobl, Austria July 3-7, 2006 Clark R. Dohrmann Sandia National
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
Well-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
A New time Integration Scheme for Cahn-hilliard Equations
A New time Integration Scheme for Cahn-hilliard Equations Item Type Conference Paper Authors Schaefer, R.; Smol-ka, M.; Dalcin, L; Paszyn'ski, M. Citation A New time Integration Scheme for Cahn-hilliard