Home  |  About Us  |  Link To Us  |  FAQ  |  Contact

# 3D Crouzeix-Raviart mortar finite element method 1.0

Date Added: July 05, 2013  |  Visits: 206

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=1This solves the Poisson problem in the domain Omega= (0,1) x (0,1) x (-1,1) assumingthe 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.

 Requirements: No special requirements Platforms: Matlab Keyword: Algorithm,  Classical,  Compared,  Functionality,  Implemented,  Interpolation,  Operator,  Quotconnectquot,  Quotstartmquot,  Solutions Users rating: 0/10