Code covered by the BSD License  

Highlights from
3D Crouzeix-Raviart mortar finite element method

3.0

3.0 | 1 rating Rate this file 16 Downloads (last 30 days) File Size: 32.8 KB File ID: #22430
image thumbnail

3D Crouzeix-Raviart mortar finite element method

by

 

13 Dec 2008 (Updated )

Implementation of 3D Crouzeix-Raviart mortar finite element

| Watch this File

File Information
Description

Software to paper (please cite when using the software)

Leszek Marcinkowski, Talal Rahman and Jan Valdman,
A 3D Crouzeix-Raviart mortar finite element.
Computing 86, No. 4, 313-330 (2009)

download: http://www.springerlink.com/content/p8nx4573480hm98h/?p=31e284c1655141f691369102e361c9a4&pi=1

This solves the Poisson problem
 in the domain Omega= (0,1) x (0,1) x (-1,1) assuming
the volume force equal to 3 pi*pi* sin(pi*x) sin(pi*y) sin(pi*z)
and the zero Dirichlet conditions.

The discrete solution is computed on two equal subdomains, Omega1= (0,1) x (0,1) x (0,1),
Omega2= (0,1) x (0,1) x (-1,0)
using 3D Crouzeix-Raviart mortar finite element method.
The subdomains are triangulated using different regular triangulation in tetrahedrons. Therefore, there are no matching grids across the intersection plane z=0.
To "connect" both subdomains solutions, a new interpolation operator is implemented and compared with a classical one.
 
To test the functionality of the algorithm, run "start.m".

It takes about 2 minutes on my PC to compute a solution on a mesh with 57216 edges. The major obstacle for the faster implementation is the use of matlab function polyxpoly works fine but is too general - it might be replaced by a faster function, which computes an intersection of two triangles in 2D.

Required Products MATLAB Compiler
MATLAB release MATLAB 7.5 (R2007b)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (1)
21 Aug 2011 Qwuasi Moore  
Updates
16 Dec 2008

improved functionality, getFaces function

04 May 2009

Description

16 Nov 2009

full article citation added, title picture changed

Contact us