Hi Birger, I think you are right. In Dumux 3.0 the volume is scv.volume() though. fvGeometry.subContVol[scvIdx].volume would be correct in Dumux 2.12 and earlier.
Timo Viele Grüße, Timo > Am 22.08.2018 um 13:08 schrieb Birger Hagemann > <[email protected]>: > > Dear Martin, > > The unit for the well index is m³. Consequently the unit for your expression > for values[contiWEEqIdx] is kg/s. However, the source term has the unit > kg/m³/s. I think you have to divide by the cell volume: > values[contiWEqIdx] = densityW/viscosityW * WI *(pbh - > pw)/fvGeometry.subContVol[scvIdx].volume; > > Am I wrong? > > Regards Birger > > Von: Dumux [mailto:[email protected]] Im Auftrag von > Martin Schneider > Gesendet: Mittwoch, 22. August 2018 10:03 > An: DuMuX User Forum <[email protected]>; Ranjeet Singh > <[email protected]> > Betreff: Re: [DuMuX] How to implement Peaceman well model in DuMuX? > > Dear Ranjeet, > > you can implement the Peaceman well model as a solution-dependent source, > using the following function in your problem file > > NumEqVector source(const Element &element, > const FVElementGeometry& fvGeometry, > const ElementVolumeVariables& > elemVolVars, > const SubControlVolume &scv) const > > (see for example porousmediumflow/2pnc/implicit/fuelcellproblem.hh) > > With the elemVolVars you have access to the pressure, density, etc.: > const auto& volVars = elemVolVars[scv]; > Scalar pw = volVars.pressure(wPhaseIdx); > Scalar densityW = volVars.density(wPhaseIdx); > Scalar viscosityW = volVars.viscosity(wPhaseIdx); > > Such that the source term looks as follows (if you inject water and neglect > gravity) > NumEqVector values(0.0); > values[contiWEqIdx] = densityW/viscosityW * WI *(pbh - pw); > > where WI is the well index (see the Peaceman well model); pbh the bottom hole > pressure, etc. > If gravity is included you have to use the phase potentials. > > Kind regards, > Martin > > > On 08/09/2018 04:03 AM, Ranjeet Singh wrote: > Hi, > I want use Peaceman well model in DuMuX?. Could you please some > information (or any example if already implemented) so that I can implement > the well model? > > > Thank you! > > > Regards, > Ranjeet > > > > > > _______________________________________________ > Dumux mailing list > [email protected] > https://listserv.uni-stuttgart.de/mailman/listinfo/dumux > > > -- > M.Sc. Martin Schneider > University of Stuttgart > Institute for Modelling Hydraulic and Environmental Systems > Department of Hydromechanics and Modelling of Hydrosystems > Pfaffenwaldring 61 > D-70569 Stuttgart > Tel: (+49) 0711/ 685-69159 > Fax: (+49) 0711/ 685-60430 > E-Mail: [email protected] > _______________________________________________ > Dumux mailing list > [email protected] > https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
_______________________________________________ Dumux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
