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]>:
> 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]>:
> >
> > > 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
> > > *********************************************
> > > _______________________________________________
> > > Dumux mailing list
> > > [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
>
>
> --
> Holger Class
>
>
> ###############################################################################
> Lehrstuhl fuer Hydromechanik und Hydrosystemmodellierung (LH2)
> Institut für Wasser- und Umweltsystemmodellierung (IWS), Universitaet
> Stuttgart
> www.iws.uni-stuttgart.de/hydrosys
> email: [email protected]
> Pfaffenwaldring 61 ** 70550 Stuttgart
> Tel.: ++49 711 / 685-64678 Fax.: ++49 711 /
> 685-60430
>
> ##################################Ecc#4,4######################################
> _______________________________________________
> Dumux mailing list
> [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
