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
