CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

icoBlockedCellFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 20, 2011, 13:40
Default icoBlockedCellFoam
  #1
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Here you can find an example how to manipulate matrices in OpenFOAM.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   April 21, 2011, 03:53
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hi,

isn't this similar to what setValues can do? http://foam.sourceforge.net/docs/cpp...rce.php#l00386

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   April 21, 2011, 04:00
Default
  #3
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
At this point it is the same, but will have further extension. Thank you for comment)
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   April 21, 2011, 04:51
Default
  #4
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Interesting. I have been hacking on that a bit myself in the past, since setValues does not take care of .A(), for some unknown reason :-)
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   May 1, 2011, 17:53
Default
  #5
New Member
 
lazio
Join Date: Apr 2010
Posts: 22
Rep Power: 7
lrunber is on a distinguished road
Thanks Dr. Alexander for your useful program, it helped me understand the method of blocking cells from the solver. I am eagerly awaiting the extension! Perhaps multiple region based blocking...

Quote:
Originally Posted by alberto View Post
Interesting. ...since setValues does not take care of .A(), for some unknown reason :-)
Alberto, could you please elaborate on "does not take care of .A()"? I am currently using setValues on internal velocity fields and would like to know of these implications I was under the impression that setValues does work on the coefficient matrix A, for the cells under question being set to the identity matrix.
lrunber is offline   Reply With Quote

Old   May 1, 2011, 18:23
Default
  #6
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by lrunber View Post
Alberto, could you please elaborate on "does not take care of .A()"? I am currently using setValues on internal velocity fields and would like to know of these implications I was under the impression that setValues does work on the coefficient matrix A, for the cells under question being set to the identity matrix.
setValues blocks the solution in the cell, however the central coefficient A() is not defined. Try for example to consider an equation where A can go to zero, remove the cells where this happens with setValues, and then compute 1/A: this will lead to a segmentation fault because A is not properly set in those cells.

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   May 2, 2011, 04:09
Default
  #7
New Member
 
lazio
Join Date: Apr 2010
Posts: 22
Rep Power: 7
lrunber is on a distinguished road
Quote:
Originally Posted by alberto View Post
setValues blocks the solution in the cell, however the central coefficient A() is not defined. Try for example to consider an equation where A can go to zero, remove the cells where this happens with setValues, and then compute 1/A: this will lead to a segmentation fault because A is not properly set in those cells.
Thanks for explaining
lrunber is offline   Reply With Quote

Reply

Tags
matrix operations

Thread Tools
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 On
Pingbacks are On
Refbacks are On



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