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]<mailto:[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]<mailto:[email protected]>
_______________________________________________ Dumux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
