Locating near wall cells for evaporation model

mehrdad_kbg · January 30, 2014, 13:34

Hi everyone,

I am trying to simulate evaporation due to high heat flux on a specific wall. I have written following lines to locate the near wall cells where the evaporation occurs.
At this moment I am only considering a constant value for evaporation (GammaGL[iCell] = 20).
I was wondering if these lines locate the near wall cells and employ the evaporation rate on the near wall cells?

Thanks for your help.

const fvPatchList& patches = mesh.boundary();

forAll(patches, patchi)
{

const fvPatch& thePatchItselfWall = patches[patchi];

if (isA<wallFvPatch>(thePatchItselfWall))
{
forAll(thePatchItselfWall, iFace)
{
label iCell = thePatchItselfWall.faceCells()[iFace];
vector iCellCentre = mesh.cellCentres()[iCell];

GammaGL[iCell] = 20;

}
}
}

Lieven · January 30, 2014, 17:09

That should work

Cheers,

Lieven

chriss85 · February 4, 2014, 02:56

I am also going to implement a similar behaviour in the future. I was planning to add a source term for rho in the first cell layer, based on the heat flux going into the wall. The wall would have a fixedValue boundary condition matching the evaporation temperature of the material. Do you think this is a possible approach or are there better methods?

January 30, 2014, 13:34	Locating near wall cells for evaporation model	#1
mehrdad_kbg New Member Join Date: Sep 2010 Location: Windsor, Ontario Posts: 17 Rep Power: 15	Hi everyone, I am trying to simulate evaporation due to high heat flux on a specific wall. I have written following lines to locate the near wall cells where the evaporation occurs. At this moment I am only considering a constant value for evaporation (GammaGL[iCell] = 20). I was wondering if these lines locate the near wall cells and employ the evaporation rate on the near wall cells? Thanks for your help. const fvPatchList& patches = mesh.boundary(); forAll(patches, patchi) { const fvPatch& thePatchItselfWall = patches[patchi]; if (isA<wallFvPatch>(thePatchItselfWall)) { forAll(thePatchItselfWall, iFace) { label iCell = thePatchItselfWall.faceCells()[iFace]; vector iCellCentre = mesh.cellCentres()[iCell]; GammaGL[iCell] = 20; } } }

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
How to handel source term in near wall cells	kaifu	OpenFOAM	1	May 15, 2011 16:59
UDF for wall slipping	HFLUENT	Fluent UDF and Scheme Programming	0	April 27, 2011 12:03
[snappyHexMesh] snappyHexMesh aborting	Tobi	OpenFOAM Meshing & Mesh Conversion	0	November 10, 2010 03:23
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug	unoder	OpenFOAM Installation	11	January 30, 2008 20:30
Wall functions	Abhijit Tilak	Main CFD Forum	6	February 5, 1999 01:16

January 30, 2014, 17:09		#2
Lieven Senior Member Lieven Join Date: Dec 2011 Location: Leuven, Belgium Posts: 299 Rep Power: 22	That should work Cheers, Lieven

February 4, 2014, 02:56		#3
chriss85 Senior Member Join Date: Oct 2013 Posts: 397 Rep Power: 18	I am also going to implement a similar behaviour in the future. I was planning to add a source term for rho in the first cell layer, based on the heat flux going into the wall. The wall would have a fixedValue boundary condition matching the evaporation temperature of the material. Do you think this is a possible approach or are there better methods?