CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Parallel plate flow - SIMPLE Algorithm - stuck residuals

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 2, 2020, 13:21
Post Parallel plate flow - SIMPLE Algorithm - stuck residuals
  #1
New Member
 
Vaibhav
Join Date: Feb 2020
Posts: 5
Rep Power: 2
vaibhavkhanna is on a distinguished road
Hey, I am quite new to CFD programming and was writing my own Matlab code for 2-D flow between parallel plates to get more familiar with the simple algorithm which I will be using for my thesis. I am trying to simulate the entrance development region which results in the parabolic velocity profile.
My mesh size is 250x10, Length of the plate is 0.5m, and distance between plates is 0.01m.
Pressure under relaxation is 0.001-0.1 for different simulations shown in the per attached.
The issue I am facing with my simulations is that the continuity residuals (L2 norm of mass imbalance) for different Re and different fluids all seem to oscillate about a fixed value instead of decreasing as the iterations proceed. The termination criteria for my simulation is never met.
Here are some of the residual plots from my simulation. The velocity values are oscillating about the right solution (the centreline velocity is 1.5 times the inlet velocity) but the residuals do not decrease with iterations.

https://drive.google.com/file/d/16fb...IUvoA7YRG/view

Last edited by vaibhavkhanna; June 2, 2020 at 21:05.
vaibhavkhanna is offline   Reply With Quote

Old   June 4, 2020, 06:48
Default
  #2
Senior Member
 
Troy Snyder
Join Date: Jul 2009
Location: Akron, OH
Posts: 205
Rep Power: 15
tas38 is on a distinguished road
Quote:
Originally Posted by vaibhavkhanna View Post
Hey, I am quite new to CFD programming and was writing my own Matlab code for 2-D flow between parallel plates to get more familiar with the simple algorithm which I will be using for my thesis. I am trying to simulate the entrance development region which results in the parabolic velocity profile.
My mesh size is 250x10, Length of the plate is 0.5m, and distance between plates is 0.01m.
Pressure under relaxation is 0.001-0.1 for different simulations shown in the per attached.
The issue I am facing with my simulations is that the continuity residuals (L2 norm of mass imbalance) for different Re and different fluids all seem to oscillate about a fixed value instead of decreasing as the iterations proceed. The termination criteria for my simulation is never met.
Here are some of the residual plots from my simulation. The velocity values are oscillating about the right solution (the centreline velocity is 1.5 times the inlet velocity) but the residuals do not decrease with iterations.

https://drive.google.com/file/d/16fb...IUvoA7YRG/view

This issue in monitoring residual convergence may be a result of monitoring "absolute" rather than "scaled" residuals. Your plot of the residuals for three different fluids suggests as much to me.


I suggest you modify your code to monitor scaled/normalized residuals of the form:


R^{\phi} = \frac{\sum_{\mbox{cells}} \left| \sum_{nb}a_{nb} \phi_nb + b - a_p \phi_p \right| }{\sum_{\mbox{cells}} \left| a_p \phi_p \right| }


Also, using pressure relaxation on the order of 0.01 or 0.001 is overly small. If this is necessary, it may be a result of some bug/error within the code for which you are compensating.
tas38 is offline   Reply With Quote

Reply

Tags
convergence., parallel plates, residuals., simple algorithm, stuck residuals.

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Flow over a flat plate as an immersed solid hamed.majeed CFX 4 September 8, 2016 15:40
Nusselt number for flow through parallel plate under periodic boundary condition Vipul_Patel FLUENT 0 October 18, 2014 15:24
Triangular mesh - fully developed flow between parallel plate Amir_Ghasemi Main CFD Forum 3 November 28, 2010 05:52
flow thru parallel plate with solid heat generater cindy FLUENT 1 October 10, 2005 23:54
Question on 3D potential flow Adrin Gharakhani Main CFD Forum 13 June 21, 1999 06:18


All times are GMT -4. The time now is 09:50.