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

Periodicity boundary implementation

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 8, 2017, 02:38
Default Periodicity boundary implementation
  #1
New Member
 
Join Date: Sep 2016
Posts: 3
Rep Power: 6
Jbsl is on a distinguished road
Hello friends,

I am working on software for solving Eulers equations in MATLAB, on structured grid, using Finite volume method with MacCormack Scheme and Jameson's artifical viscosity. I successfully ran program for GAMM channel test case. Now im trying to apply my software on different geometry, namely turbine cascade. For that i need to used periodicity boundary and i feel like I am implementing it wrong. So far structure of my program looks something like this:


-set inlet/outlet boundary
-set periodicity so that for grid enlarged by two ghost cells i have 1:ny+2 and i set W(:,1) = W(:,ny+1) and W(:,ny+2) = W(:,2) where W is vector of conservative values
-calculate flux
-calculate predictor
-set periodicity again for predicted values same way
-recalculate flux using values from predictor
-calculate corrector
-set ghost cells on inlet/outlet
-set periodicity for enlarged matrice with 4 ghost cells (for Jameson) same way as above.
-calculate Jameson
-repeat

I also apply wall boundary, but since it worked for me before i don't think there should be problem.
At this point program does about 2000 iterations, residual is decresing, then it drop to complex values.
My question is if my implementation of boundary is right. If yes, i would be glad for any advice on where can be mistake in my software, since i'm really stuck.

Thanks for all the help and i'm sorry for my English.
Jbsl is offline   Reply With Quote

Old   May 8, 2017, 09:47
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,416
Rep Power: 58
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I suggest to search for similar posts already answered
FMDenaro is offline   Reply With Quote

Old   May 8, 2017, 15:09
Default
  #3
New Member
 
Join Date: Sep 2016
Posts: 3
Rep Power: 6
Jbsl is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I suggest to search for similar posts already answered
Hello,
that's the point, i searched the forum and thought i understand how to implement the boundary. It even seems to be working, however i still get complex values. Thats why i posted here I don't get why residual is decresing and then suddenly i get complex values. All i wanted was to check, if my boundaries were implemented correctly.
Jbsl is offline   Reply With Quote

Reply

Tags
boundary condition, matlab, periodicity

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
My radial inflow turbine Abo Anas CFX 27 May 11, 2018 02:44
Multiphase flow - incorrect velocity on inlet Mike_Tom CFX 6 September 29, 2016 02:27
sliding mesh problem in CFX Saima CFX 45 September 22, 2015 11:53
CFX13 Post Periodic interface EtaEta CFX 7 December 8, 2011 18:15
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Meshing & Mesh Conversion 2 July 15, 2005 05:15


All times are GMT -4. The time now is 13:05.