CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Interpolation across regionCouple patches in 1.6-ext

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

Like Tree2Likes
  • 1 Post By Ola Widlund
  • 1 Post By sebware

LinkBack Thread Tools Display Modes
Old   November 28, 2011, 10:58
Default Interpolation across regionCouple patches in 1.6-ext
New Member
Ola Widlund
Join Date: Jun 2011
Location: Sweden
Posts: 13
Rep Power: 8
Ola Widlund is on a distinguished road

My question concerns interpolation across regionCouple patches in 1.6-ext. There are a number of related threads in the forum, most of which you will find if you search for "regionCoupling".

I successfully use OpenFOAM for solving ion transport problems in insulation material, now also including electrohydrodynamics effects (ion winds).

When I started looking at the implicit region coupling in 1.6-ext, I thought that the coupling across regions was implemented to mimic the processor patches when the regionCoupling patch is in its "attached" state. I liked this, because it would mean that interpolation operations (fvc::interpolation) across coupled patches would use whatever interpolation scheme I select in the fvSchemes dictionary, just as for processor patches.

Now that I dig deeper into the code, and examine the results more closely, I realize that the regionCoupling patches don't work that way. The value at the boundary is hard-coded to harmonic averaging in the regionCoupling patch field class. The fvc::interpolate operator would return a "calculated" patch with this harmonic interpolate, regardless of the fvSchemes settings. This works well for most applications where you have diffusion across material interfaces, but also has some limitations. Without going into details, it is not necessarily what you want when you discretize charge transport across a material interface.

I can easily hack a fix for this (looping and setting the values I want on all boundaries), but I wonder if there is a very fundamental reason that I'm missing for choosing this implementation. Or if it could be done in the way I wanted it...

If the regionCoupled patch (in its attached state) was to work like the processor patch, it would be enough to choose "harmonic" interpolation for the laplacian to reproduce the same results for the conjugate heat transfer solver in 1.6-ext, without compromising the larger generality of the underlying coupledFvPatch class.

It would be interesting to have some comments on this.

Best regards,

mm.abdollahzadeh likes this.
Ola Widlund is offline   Reply With Quote

Old   May 8, 2012, 09:08
Default regionCoupleFvPatch Theory Questions
New Member
Join Date: Feb 2010
Posts: 9
Rep Power: 10
sebware is on a distinguished road

I plan to couple the magnetic vectorpotential with a regionCoupleFvPatch.
With the template technology it was very easy to built a vector version of the couplingRegion boundary.

But now i'll plan todo my own interpolation to calculate the patchField.
So i try to understand, on which position in the cells are the different fields are located.

I hope somebody correct me or give additional informations

I use:

1. the distance of cellCenter(owner) to interface, i think that's the patch().deltaCoeffs() isnt't it ?

2. the distance of interface to cellCenter(neigbour), i dont have a plan about this :-(

3. the field in the cellCenter(Owner) i think that's comes from this->patchInternalField()

4. the field in the cellCenter(Neigbour). That must be the this->patchNeigbourField()

But there is also a access funktion for the shadowPatchField(). Where is the location of this field?


best regards
sebware is offline   Reply With Quote

Old   May 8, 2012, 09:55
Default Additional sketches
New Member
Join Date: Feb 2010
Posts: 9
Rep Power: 10
sebware is on a distinguished road
here two little sketches to show my opinion of the field locations and my used values.

best regards
Attached Images
File Type: jpg of_coupleRegionfvPatch.jpg (78.5 KB, 129 views)
File Type: jpg of_deltaCoeffs.jpg (38.9 KB, 129 views)
mm.abdollahzadeh likes this.
sebware is offline   Reply With Quote

Old   June 18, 2013, 14:40
Senior Member
mahdi abdollahzadeh
Join Date: Mar 2011
Location: Covilha,Portugal
Posts: 151
Rep Power: 8
mm.abdollahzadeh is on a distinguished road
hi sebastian

Could please tell me how I could access in DEltaX1 and DeltaX2 in regioncoupling?

mm.abdollahzadeh is offline   Reply With Quote


coupled patches, regioncouple, regioncoupling

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

Similar Threads
Thread Thread Starter Forum Replies Last Post
Handling cyclic BC from gambit to openfoam for a cascade airfoil problem - OF 1.6 maverick OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 June 18, 2011 04:36
Test directory missing in OpenFOAM 1.6 ext andrewryan OpenFOAM 2 March 20, 2011 16:41
Surface interpolation schemes and parallelization jutta OpenFOAM Running, Solving & CFD 0 February 25, 2010 15:32
Cross-compiling OpenFOAM 1.6 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 7 January 19, 2010 16:39
OpenFOAM 1.6 CreatePatch Problem TarifaPirata OpenFOAM Bugs 1 September 10, 2009 04:35

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