|
[Sponsors] |
need help about ConjugateHeatSimpleFoam solver source term at the interface |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Gaoqiang Yang
Join Date: Apr 2018
Posts: 14
Rep Power: 8 ![]() |
Hi, Foamers
In the class chtRcTemperatureFvPatchScalarField'' there is a member function: Source(), defined as follows. { const fvPatch& p = patch(); const magLongDelta& mld = magLongDelta::New(p.boundaryMesh().mesh()); const scalarField TcOwn = patchInternalField(); const scalarField TcNei = patchNeighbourField(); const scalarField& Tw = *this; const chtRegionCoupleBase& K = dynamic_cast<const chtRegionCoupleBase&> ( p.lookupPatchField<volScalarField, scalar>(kName_) ); const scalarField k = K*p.deltaCoeffs(); const scalarField kOwn = K.originalPatchField()/(1 - p.weights())/mld.magDelta(p.index()); return kOwn*(Tw - TcOwn) - k*(TcNei - TcOwn); } I am confused about the last line of code. As far as I understand, kOwn = thermalConductivity at cellCenter divided by the distance from the inrerface to faceCellCenter. k = thermalConductivity at the interface divided by the distance from Owner to Neighbour cell center. Then I should return only - k*(TcNei - TcOwn); Can anybody help me? Thanks a lot. |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] Installation Problem with OF 6 version | Aurel | OpenFOAM Community Contributions | 14 | November 18, 2020 16:18 |
Question about adaptive timestepping | Guille1811 | CFX | 25 | November 12, 2017 17:38 |
[swak4Foam] build problem swak4Foam OF 2.2.0 | mcathela | OpenFOAM Community Contributions | 14 | April 23, 2013 13:59 |
[swak4Foam] funkySetFields compilation error | tayo | OpenFOAM Community Contributions | 39 | December 3, 2012 05:18 |
DxFoam reader update | hjasak | OpenFOAM Post-Processing | 69 | April 24, 2008 01:24 |