|
[Sponsors] |
June 20, 2013, 13:26 |
RegionCoupling - openfoam 1.6ext
|
#1 |
Senior Member
mahdi abdollahzadeh
Join Date: Mar 2011
Location: Covilha,Portugal
Posts: 153
Rep Power: 15 |
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: -K1grad(T1)=Ksgrad(Ts)+sigma 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 Code:
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(patch.nf() & patch.delta()); 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 best mahdi |
|
June 20, 2013, 13:44 |
|
#2 |
Senior Member
mahdi abdollahzadeh
Join Date: Mar 2011
Location: Covilha,Portugal
Posts: 153
Rep Power: 15 |
please give me some idea
|
|
April 30, 2022, 08:08 |
|
#3 |
New Member
Zeinab Abosedaira
Join Date: Mar 2022
Posts: 24
Rep Power: 4 |
Hello Dr. Mahdi,
Have you found any solution for adding a refGradient value so we satisfy the harmonic mean in the interface? I am working on a similar case with Suzan model plasma. Thank you! |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Gmsh] 2D Mesh Generation Tutorial for GMSH | aeroslacker | OpenFOAM Meshing & Mesh Conversion | 12 | January 19, 2012 03: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 12: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 |