CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

buoyantPimpleDyMFoam dimension mismatch between Uf and phi

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 16, 2015, 04:30
Default buoyantPimpleDyMFoam dimension mismatch between Uf and phi
  #1
New Member
 
Sumit Singh
Join Date: Nov 2014
Posts: 5
Rep Power: 9
Sumit07 is on a distinguished road
I have combined buoyantPimpleFoam with pimpleDyMFoam to create new solver buoyantPimpleDyMFoam. But when running case using this solver then it gives an error that
Different dimensions for =
dimensions : [1 0 -1 0 0 0 0] = [0 3 -1 0 0 0 0]
at line where phi is calculating
phi = mesh.Sf() & Uf;
Now when I look in this equation Uf have dimension of [0 1 -1 0 0 0 0] and mesh.Sf() have dimension of [0 2 1 0 0 0 0].
But as this equation is a dot product of mesh.Sf() & Uf , so phi should have dimension of [0 1 -1 0 0 0 0] but instead it has dimension [1 0 -1 0 0 0 0].
Can any body knows why this is coming?
Sumit07 is offline   Reply With Quote

Old   June 16, 2015, 16:02
Default
  #2
Senior Member
 
Przemek
Join Date: Jun 2011
Posts: 241
Rep Power: 13
gaza is on a distinguished road
Hi Sumit
Dimensions [1 0 -1 0 0 0 0] are kg/s. It seems that phi is calculated as for compressible fluid. Phi in incompressible flow has dimensions [0 3 -1 0 0 0 0]. BuoyantPimpleFoam is for compressible fluid, while pimpleDyMFoam is for incompressible. You should merge solvers for the same type of fluid. So it looks like LHS is calculated for compressible fluid and RHS for incompressible.
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   June 17, 2015, 04:55
Default
  #3
New Member
 
Sumit Singh
Join Date: Nov 2014
Posts: 5
Rep Power: 9
Sumit07 is on a distinguished road
Thanks Przemek
I have found my error.
The error is present in createUf.H and fvc::makeRelative(phi, U) which are incopressible library.
So I created my own compressible library compressibleCreateUf.H and compressible version of makeRelative(phi, U) i.e. fvc::makeRelative(phi, rho, U).
Sumit07 is offline   Reply With Quote

Reply

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



All times are GMT -4. The time now is 20:34.