icoBlockedCellFoam

makaveli_lcf · April 20, 2011, 14:40

Here you can find an example how to manipulate matrices in OpenFOAM.

alberto · April 21, 2011, 04:53

Hi,

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

Best,

makaveli_lcf · April 21, 2011, 05:00

At this point it is the same, but will have further extension. Thank you for comment)

alberto · April 21, 2011, 05:51

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 :-)

lrunber · May 1, 2011, 18:53

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

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.

alberto · May 1, 2011, 19:23

Quote:

Originally Posted by lrunber

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,

lrunber · May 2, 2011, 05:09

Quote:

Originally Posted by alberto

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

April 20, 2011, 14:40	icoBlockedCellFoam	#1
makaveli_lcf Senior Member Dr. Alexander Vakhrushev Join Date: Mar 2009 Posts: 256 Blog Entries: 1 Rep Power: 19	Here you can find an example how to manipulate matrices in OpenFOAM. __________________ Best regards, Dr. Alexander VAKHRUSHEV Christian Doppler Laboratory for "Metallurgical Applications of Magnetohydrodynamics" Simulation and Modelling of Metallurgical Processes Department of Metallurgy University of Leoben http://smmp.unileoben.ac.at

April 21, 2011, 04:53		#2
alberto Senior Member Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,912 Rep Power: 36	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 in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.

April 21, 2011, 05:00		#3
makaveli_lcf Senior Member Dr. Alexander Vakhrushev Join Date: Mar 2009 Posts: 256 Blog Entries: 1 Rep Power: 19	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 "Metallurgical Applications of Magnetohydrodynamics" Simulation and Modelling of Metallurgical Processes Department of Metallurgy University of Leoben http://smmp.unileoben.ac.at

April 21, 2011, 05:51		#4
alberto Senior Member Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,912 Rep Power: 36	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 in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.