Hi Beatrix,
yes, the amount of a produced fluid phase depends on its mobility within
the cell where you have your well. Thus, the amount of a produced
component is the sum over all (produced) phases multiplied by the
concentration. For a two phase model it looks like this:
for each component in your model:
values[conti0EqIdx+compIdx]=(fs.moleFraction(nPhaseIdx,compIdx)*mobilityGas_*WellIndex_*(productionPressure_-fs.pressure(nPhaseIdx))
+
fs.moleFraction(wPhaseIdx,compIdx)*mobilityWater_*WellIndex_*(productionPressure_-fs.pressure(wPhaseIdx)))*(1/fvGeometry.subContVol[scvIdx].volume);
Regards Birger
Am 30.10.2015 um 21:10 schrieb Beatrix:
Hi Birger,
out of curiosity and because I might need it in the future too, what
changes if you want to produce? I guess you need to take the presence
of another phase into account somehow because it might get produced as
well.
Cheers
Beatrix
On 30.10.2015 12:50, Birger Hagemann wrote:
Hi Tri Dat,
exactly like this we have implemented a well model (Peaceman) in
solDependentSource. Here you finde the code for an H2 injection well:
Scalar
mobilityGas_=elemVolVars[scvIdx].relativePermeability(nPhaseIdx)/fs.viscosity(nPhaseIdx)*fs.molarDensity(nPhaseIdx);
Scalar
WellIndex_=((2*M_PI*elemVolVars[scvIdx].permX()*cellHeight_)/(log(effectiveWellRadius_/(WellDiameter))));
values[contiH2EqIdx]=mobilityGas_*WellIndex_*(InjectionPressure_-fs.pressure(nPhaseIdx))*(1/fvGeometry.subContVol[scvIdx].volume);
Or do you want to produce?
Regards Birger
On Fri, 30 Oct 2015 16:54:47 +0100
Kilian Weishaupt <[email protected]> wrote:
Hello Tri Dat,
I think you can try the solDependentSource for this purpose. If you
specify an injection rate like q= factor*(p - p_well), where p is
your actual pressure of your control volume and p_well is the
desired pressure, your pressure should go towards your desired
value. The factor might be a constant but it can also depend on
things like saturation, ect.
Best regards
Kilian
On 10/30/2015 04:33 PM, Tri Dat NGO wrote:
Hello Dumuxers,
I want to implement a simplified well model in which a constant
pressure is set for some given cells within a 3D domain (not at
boundaries). So I would like to know if it is possible to do this
in 2p (or 2p2c)implicit CC model.
I really appreciate any help you can provide.
Kind regards,
Tri Dat
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Kilian Weishaupt M.Sc.
Institut für Wasser- und Umweltsystemmodellierung (IWS)
Lehrstuhl für Hydromechanik und Hydrosystemmodellierung
Universität Stuttgart, Pfaffenwaldring 61, 70569 Stuttgart
Email: [email protected]
Telefon: 0049 711 685-60461 ** fax: 0049-711-685-60430
http://www.hydrosys.uni-stuttgart.de
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
_______________________________________________
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
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
