|
[Sponsors] |
Problem with implementation of SUPG in FEM for flow problems |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 1, 2015, 04:34 |
Problem with implementation of SUPG in FEM for flow problems
|
#1 |
New Member
Mahdi Salami
Join Date: Feb 2014
Posts: 6
Rep Power: 12 |
Dear Guys,
I'm trying to implement SUPG stabilization to FEM code for flow problems. I tried it in 1D and 2D for scalar variable and it is OK. Then I tried to use it for flow analysis where problems started. I tried many strategies proposed but I couldn't make it. I think I'm madding a logical mistake in implementing SUPG. I'm attaching the code to get your opinions. The code works well for low Re but when it comes to higher Re, it fails. SUPG was implemented in int_KF2.m function. The code is based on mixed formulation and 9-point quadrilateral Taylor-Hood element was used. Just rum main.m and see the results written in fem.dat file in TECPLOT. The SUPG parameter was calculated as proposed by Tezduyar. I would appreciated any comment on the code... Tnx PS: It is better to say that I implemented SU. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
FEM implementation in OpenFOAM 2.0.0 or later? | seboxx | OpenFOAM | 3 | August 15, 2011 07:24 |
non-linear k-epsilon model implementation problems | Saidi | Main CFD Forum | 2 | March 4, 2010 13:23 |
Problems with SUPG body force term | FEM question | Main CFD Forum | 0 | January 21, 2006 17:51 |
How good is CFD? | kai | Main CFD Forum | 39 | April 7, 2000 12:48 |
extremely simple problem... can you solve it properly? | Mikhail | Main CFD Forum | 40 | September 9, 1999 09:11 |