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

RegionCoupling - openfoam 1.6ext

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

LinkBack Thread Tools Display Modes
Old   June 20, 2013, 13:26
Default RegionCoupling - openfoam 1.6ext
Senior Member
mahdi abdollahzadeh
Join Date: Mar 2011
Location: Covilha,Portugal
Posts: 151
Rep Power: 8
mm.abdollahzadeh is on a distinguished road
Dear All

Im trying to solve my previous case which was based on chtmultiregionfoam with conjucateheatfoam.
my boundary at the interface of solid and fluid looks like:


I dont use coorectboundarycondition for K . so I dont have harmonic interpolation of K. and for T I am doing the interpolation similar to energyRegionCoupledFvPatchScalarField in openfoam 2.2

 Field<Type> fOwn = this->patchInternalField();
    Field<Type> fNei = this->patchNeighbourField();
    Field<Type> fOwnsigma = this->patchsigmaField();
    Field<Type> fNeisigma = this->patchsigmaNeighbourField();

    Field<Type> Alpha=this->patchDTField();
    scalarField fAlpha=mag(Alpha);
deltas( &;     const fvPatch& patch = regionCoupledPatch_.patch();
    const scalarField alphaDelta(fAlpha/deltasii);
    Field<Type> nbrAlpha=this->patchDTNeighbourField();
    scalarField fnbrAlpha=mag(nbrAlpha);
    const scalarField nbrAlphaDelta(fnbrAlpha/nbrDeltasii);
    scalarField weightsd(fOwn.size());  
    scalarField gweights(fOwn.size());  

    forAll (weightsd, faceI) 
        scalar di = alphaDelta[faceI];
        scalar dni = nbrAlphaDelta[faceI];
        weightsd[faceI] = di/(di + dni);
        gweights[faceI] = 1.0/(di + dni);
if (sigmaflag_==1)
    // Do interpolation
   Field<Type>::operator=(weightsd*fOwn + (1.0 - weightsd)*fNei-fOwnsigma*sigmaflag_*gweights-fNeisigma*sigmaflag_*gweights); 

which sigma is just defined on surface of the interface. in the correct result I should get some sort of generation due to sigma and the internal field should be effected by that.

the boundary values that I get after modifying regioncoupling is the same as the values at the boundary that I get by chtmultiregion .
but its not puting any effect on the internal field !!!!!.

I would appreciate any comment

Attached Images
File Type: jpg ttttt.jpg (26.9 KB, 10 views)
File Type: jpg ggg.jpg (27.0 KB, 8 views)
mm.abdollahzadeh is offline   Reply With Quote

Old   June 20, 2013, 13:44
Senior Member
mahdi abdollahzadeh
Join Date: Mar 2011
Location: Covilha,Portugal
Posts: 151
Rep Power: 8
mm.abdollahzadeh is on a distinguished road
please give me some idea
mm.abdollahzadeh is offline   Reply With Quote


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
2D Mesh Generation Tutorial for GMSH aeroslacker Open Source Meshers: Gmsh, Netgen, CGNS, ... 12 January 19, 2012 04:52
Problem installing OpenFOAM 1.5 installation on RHEL 4. vwsj84 OpenFOAM Installation 4 April 23, 2009 04:48
2009 OpenFOAM Summer School in Zagreb, Croatia hjasak OpenFOAM Announcements from Other Sources 0 March 27, 2009 13:08
64bitrhel5 OF installation instructions mirko OpenFOAM Installation 2 August 12, 2008 18:07
OpenFOAM Training and Workshop Hrvoje Jasak Main CFD Forum 0 October 7, 2005 07:14

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