|
[Sponsors] |
December 2, 2011, 10:03 |
Wall function BC:s on regionCouple patches
|
#1 |
New Member
Ola Widlund
Join Date: Jun 2011
Location: Sweden
Posts: 13
Rep Power: 15 |
Hi!
I'm using the regionCouple functionality in 1.6-ext for problems in electrohydrodynamics, involving fluid flow and charge transport across fluid/solid interfaces. If my flow is turbulent, I want to use the turbulence wall function boundary conditions implemented in classes epsilonWallFunctionFvPatchScalarField, nutWallFunctionFvPatchScalarField, etc. The problem is that these BC classes require the coupled boundary to descend from the wallFvPatch class, rather than from the fvPatch class (via coupledFvPatch). In the explicit coupling which is part of the OpenCFD distribution solves this by introducing a whole line of directMappedWall classes on different levels, as opposed to the directMapped class which holds the basic inter-region coupling functionality. I think my problems is best solved in the same way. That is, by implementing a bunch of "wall" classes, bottom up, to do the same for the otherwise already excellent implicit regionCouple functionality. It's not extreme amounts of work, but anyway: Has someone already done it? Otherwise my weekend is saved... /Ola |
|
December 6, 2011, 07:10 |
|
#2 |
New Member
Ola Widlund
Join Date: Jun 2011
Location: Sweden
Posts: 13
Rep Power: 15 |
It seems I'm the only one who thinks this is an interesting topic...
Ok, I've done the new regionCouple classes needed, to make the coupled boundaries masquerade as wall patches, so that wall function BC:s work. It seems to work, no particular surprises. The new classes are: regionCoupleWallPolyPatch, regionCoupleWallFvPatch, coupledWallPolyPatch, coupledWallFvPatch and a bunch of templated regionCoupleWallFvPatchFields. It would be great if I could also make the GAMG solver work, at least for the un-coupled flow problem, when the regionCouple patches are detached. For that to work, one would presumably have to define a new GAMGinterface class to handle regionCouple patches. I see there are already GAMGinterfaces for the Ggi stuff, so it shouldn't be impossible. At a detached regionCouple patch, there is not even any action across the boundary, so it should be enough to fool GAMG to believe it's a wall. To make GAMG work for the coupled matrices as a whole is of course another story! Could it be done, or are there any fundamental issues that can't be overcome? Best regards, Ola |
|
February 1, 2013, 07:24 |
|
#3 |
Member
Elisabet Mas de les Valls
Join Date: Mar 2009
Location: Barcelona, Spain
Posts: 64
Rep Power: 17 |
Dear Ola and foamers,
I've encountered the same problem of using regionCouple and needing wall functions on the same patch. Following Ola's instructions, I've spent some hours solving the problem and creating the new b.c. that uses the wallPolyPatch class. After succeeding and testing the new b.c. my case was ready to be sent to the cluster of PC. The case is running now perfectly (as far as I know). But then, suddenly, I remembered that in the cluster I had a newer OF version (still 1.6-ext, but a recent one) where the regionCouple b.c. asks for the 'isWall' input label. I have just checked the code and found that if this 'isWall' label is set to 'yes', then wall functions can perfectly be used!!!!! So, don't repeat my errors!!! |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] mesh airfoil NACA0012 | anand_30 | OpenFOAM Meshing & Mesh Conversion | 13 | March 7, 2022 17:22 |
Need some wall function approaches! | yka8150 | Main CFD Forum | 0 | September 21, 2009 23:08 |
Wall function problem in Fluent | mefpz | FLUENT | 1 | October 10, 2007 13:43 |
Droplet Evaporation | Christian | Main CFD Forum | 2 | February 27, 2007 06:27 |
Wall function in adverse pressure gradients | stephane baralon | Main CFD Forum | 11 | September 2, 1999 04:05 |