Dear Gerado,
I dont know very well how to use de return values of this function.
Depending on whether useMoles is set on true or false, the flux has to
be given either in * kg/(m^3*s) or mole/(m^3*s) in the input file. So
you have to divide your extraction rate by the volume of the element(s)
you are injecting into.
If you are not sure whether you have the right position for injection,
just put an output into your if statement:
if (globalPos[1] < 21 && globalPos[1] > 20 && globalPos[0] > 9 &&
globalPos[0] < 10)
std::cout<<"Im in the injection spot"<<std::endl;
values= -1;
else
values= 0;}
best regards,
Alex
On 02/13/2014 04:01 PM, Gerardo Zegers R wrote:
Hi.
I want to put a extraction well of water in a model. To do this i use
function sourceAtPos ( below is the function), the model is running
but it´s not extracting only water. How can I put that extraction of
water qw=-1 and extraction of air qa=0 . I dont know very well how to
use de return values of this function.
Thanks
Gerardo
void sourceAtPos(PrimaryVariables &values,
const GlobalPosition &globalPos) const
{
const Scalar time = this->timeManager().time() +
this->timeManager().timeStepSize();
if(time>5000)
{
if (globalPos[1] < 21 && globalPos[1] > 20 && globalPos[0] > 9 &&
globalPos[0] < 10)
values= -1;
else
values= 0;}
else
{values= 0;}
}
2014-02-12 12:24 GMT+01:00 Holger Class
<[email protected]
<mailto:[email protected]>>:
Dear Gerardo,
problems with convergence can have thousands of reasons, most
common is
some problem with non-physical boundary and initial values or
other model
parameters.
Try to get a simple problem running, increase step by step the
complexity
of your physics, e.g. taking a very small value of the source term
first,
and so on.
Without knowing this, we cannot help in this situation. At least I
cannot.
Best regards,
Holger
On Wed, Feb 12, 2014 at 12:15:07PM +0100, Gerardo Zegers R wrote:
> Hi thanks for the help and sorry for not give more details.
>
> When i used the function SourceAtPos , the problem can compile
without
> problem but when I run the file It cant converge.
>
> The describing area is inside the grid and they are nodes inside
because
> dx=0.4 m and dy=0.4 m and iam using now a subarea of 1x1m to put
the source
> term.
>
> What can i do to get convergence? i have tried with differents
timesteps.
>
> Thanks
>
> Gerardo
>
>
>
> 2014-02-12 9:19 GMT+01:00 Christoph Grüninger <
> [email protected]
<mailto:[email protected]>>:
>
> > Hi Gerardo,
> > "does not work" is kind of an useless problem description.
What do you
> > expect? What happens actually?
> >
> > Having a look into my crystal ball unveils the following hint:
> >
> > - Adding a std::cout to your function to proof that it is actually
> > evaluated.
> > - Your grid must match the conditions. It must be fine enough
that some
> > nodes are inside the described area and that this area is
inside the
> > grid. Please double check. Maybe a std::cout inside the first
case would
> > be helpful, too.
> > - Have a look at the resulting linear system, especially the
right hand
> > side. Does your source term have an effect?
> >
> > Bye
> > Christoph
> >
> > --
> > ... as Sir Cyril Hinshelwood has observed ...fluid dynamicists
> > were divided into hydraulic engineers who observed things that
> > could not be explained and mathematicians who explained things
> > that could not be observed. -- James Lighthill
> > *********************************************
> > CMWR 2014: 10th - 13th June 2014 in Stuttgart
> > Please visit www.cmwr14.de <http://www.cmwr14.de>
> > *********************************************
> > _______________________________________________
> > Dumux mailing list
> > [email protected]
<mailto:[email protected]>
> > https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
> >
>
>
>
> --
> Gerardo Zegers R
> 89850305
> _______________________________________________
> Dumux mailing list
> [email protected]
<mailto:[email protected]>
> https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
Holger Class
###############################################################################
Lehrstuhl fuer Hydromechanik und Hydrosystemmodellierung (LH2)
Institut für Wasser- und Umweltsystemmodellierung (IWS),
Universitaet Stuttgart
www.iws.uni-stuttgart.de/hydrosys
<http://www.iws.uni-stuttgart.de/hydrosys>
email: [email protected] <mailto:[email protected]>
Pfaffenwaldring 61 ** 70550 Stuttgart
Tel.: ++49 711 / 685-64678
<tel:%2B%2B49%20%20711%20%2F%20685-64678> Fax.: ++49
711 / 685-60430 <tel:%2B%2B49%20%20711%20%2F%20685-60430>
##################################Ecc#4,4######################################
_______________________________________________
Dumux mailing list
[email protected]
<mailto:[email protected]>
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
Gerardo Zegers R
89850305
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
*******************************************************
!!!! CMWR 2014: 10th - 13th June 2014 in Stuttgart !!!!
Please visit www.cmwr14.de
*******************************************************
Alexander Kissinger
Institut für Wasser- und Umweltsystemmodellierung
Lehrstuhl für Hydromechanik und Hydrosystemmodellierung
Pfaffenwaldring 61
D-70569 Stuttgart
Telefon: +49 (0) 711 685-64729
E-Mail: [email protected]
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
