[aspect-devel] Q1 P0 finite element

[Thieulot, C.A.P. (Cedric)]

Hi Harsha, 

it is, and if I am not mistaken here’s what you need in the input file:

set Stokes velocity polynomial degree = 1 

and 

set Use locally conservative discretization = true

Cheers, 

Ce/





On 31 Aug 2016, at 00:58, Harsha Lokavarapu <lokavarapuh at gmail.com> wrote:

> Hi,
> 
> In the paper, Discretization errors in the hybrid finite element particle-in-cell by M. Thielmann, D.A. May, and B.J. P. Kaus, they analyze both cell averaging and a bilinear least squares as an interpolationscheme for the benchmarks SolCx, SolKz, and SolVI. They further analyze the algorithms by measuring the convergence rates for different finite elements including Q1xP0 and Q2xP{-1}.
> 
> Is it possible to run ASPECT using Q1xP0 finite element? 
> 
> Thank You,
> -Harsha
> 
> 
> 
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel