Hi Martin,
ok, so i corrected the condition, so now, my code is:
void boundaryTypesAtPos(BoundaryTypes &bcTypes, const GlobalPosition&
globalPos) const /*@\label{tutorial-decoupled:bctype}@*/
{
if ((globalPos[0] > 580- eps_) && (globalPos[1] > 580-
eps_) )
{
bcTypes.setDirichlet(pressEqIdx);
bcTypes.setDirichlet(satEqIdx);
//bcTypes.setAllDirichlet(); // alternative if the same
BC is used for both types of equations
}
// all other boundaries
else if ((globalPos[0] < 20- eps_) && (globalPos[1] <
20- eps_) )
{
//bcTypes.setNeumann(pressEqIdx);
//bcTypes.setDirichlet(satEqIdx);
bcTypes.setAllNeumann(); // alternative if the same BC
is used for both types of equations
}
else
bcTypes.setAllNeumann();
}
//! Value for dirichlet boundary condition at position globalPos.
/*! In case of a dirichlet BC for the pressure equation the
pressure \f$ [Pa] \f$, and for
* the transport equation the saturation [-] have to be defined
on boundaries.
*
* \param values Values of primary variables at the boundary
* \param intersection The boundary intersection
*
* Alternatively, the function dirichletAtPos(PrimaryVariables
&values, const GlobalPosition& globalPos)
* could be defined, where globalPos is the vector including the
global coordinates of the finite volume.
*/
void dirichletAtPos(PrimaryVariables &values, const
GlobalPosition& globalPos) const
{
values=0;
if ((globalPos[0] > 580- eps_) && (globalPos[1] > 580- eps_) )
{
values[pwIdx] = 3.5e7;
values[swIdx] = 0.0;
}
else if ((globalPos[0] < 20- eps_) && (globalPos[1] < 20-
eps_) )
values[swIdx] = 1.0;
}
//! Value for neumann boundary condition \f$ [\frac{kg}{m^3 \cdot
s}] \f$ at position globalPos.
/*! In case of a neumann boundary condition, the flux of matter
* is returned as a vector.
*
* \param values Boundary flux values for the different phases
* \param globalPos The position of the center of the finite volume
*
* Alternatively, the function neumann(PrimaryVariables &values,
const Intersection& intersection) could be defined,
* where intersection is the boundary intersection.
*/
void neumannAtPos(PrimaryVariables &values, const GlobalPosition&
globalPos) const /*@\label{tutorial-decoupled:neumann}@*/
{
values = 0;
if ((globalPos[0] < 20- eps_) && (globalPos[1] < 20- eps_) )
{
values[nPhaseIdx] = -1e-8;
//values[wPhaseIdx] = 0.0;
}
else
{
values[nPhaseIdx] = 0.0;
values[wPhaseIdx] = 0.0;
}
}
//! Initial condition at position globalPos.
/*! Only initial values for saturation have to be given!
*
* \param values Values of primary variables
* \param element The finite volume element
*
* Alternatively, the function initialAtPos(PrimaryVariables
&values, const GlobalPosition& globalPos)
* could be defined, where globalPos is the vector including the
global coordinates of the finite volume.
*/
void initial(PrimaryVariables &values,
const Element &element) const
/*@\label{tutorial-decoupled:initial}@*/
{
values = 0;
}
private:
const Scalar eps_;
};
} //end namespace
with eps_=1e-6
and now i have this error
Don't panic... !
Rank 0: No parameter file given. Defaulting to
'./tutorial_decoupled.input' for input file.
Initializing problem 'tutorial_decoupled'
Dune reported error: ISTLError
[apply:/home/latifa/Dumux_2.6.0/dune-istl-2.3.1/dune/istl/solvers.hh:651]:
breakdown in BiCGSTAB - rho 0 <= EPSILON 1e-80 after 8.5 iterations
i don't understand why i have this error, and what does it mean?
2015-10-21 20:06 GMT+02:00 Martin Schneider
<[email protected]
<mailto:[email protected]>>:
Hi,
the first if-query is wrong, corresponding to your boundary
conditions:
if x < 20m and y < 20 m: q_w.n = -1e-8 and Sw=1
if x > 580 and y > 580: pw= 150 bar and Sn=1
it should be
if ((globalPos[0] > 580- eps_) && (globalPos[1] >
580- eps_) )
{
bcTypes.setDirichlet(pressEqIdx);
bcTypes.setDirichlet(satEqIdx);
//bcTypes.setAllDirichlet(); // alternative if the
same BC is used for both types of equations
}
instead of:
if ((globalPos[0] > 600- eps_) && (globalPos[1] >
600- eps_) )
{
bcTypes.setDirichlet(pressEqIdx);
bcTypes.setDirichlet(satEqIdx);
//bcTypes.setAllDirichlet(); // alternative if the
same BC is used for both types of equations
}
Regards,
Martin
On 10/21/2015 07:34 PM, Ait Mahiout Latifa wrote:
Hi,
i want to change the boundary conditions in tutorial_decoupled
problem, so that in an domain 600*600, are imposed the boundary
conditions:
if x < 20m and y < 20 m: q_w.n = -1e-8 and Sw=1
if x > 580 and y > 580: pw= 150 bar and Sn=1
and no flux in the other parts of the domain.
So, i change x and y in the .input file, and i have the folowing
modifications in the tutorial_decoupled.hh:
void boundaryTypesAtPos(BoundaryTypes &bcTypes, const
GlobalPosition& globalPos) const
/*@\label{tutorial-decoupled:bctype}@*/
{
if ((globalPos[0] > 600- eps_) && (globalPos[1] >
600- eps_) )
{
bcTypes.setDirichlet(pressEqIdx);
bcTypes.setDirichlet(satEqIdx);
//bcTypes.setAllDirichlet(); // alternative if the same BC is
used for both types of equations
}
// all other boundaries
else if ((globalPos[0] < 20- eps_) &&
(globalPos[1] < 20- eps_) )
{
bcTypes.setNeumann(pressEqIdx);
bcTypes.setDirichlet(satEqIdx);
//bcTypes.setAllNeumann(); // alternative if the same BC is used
for both types of equations
}
else
bcTypes.setAllNeumann();
}
//! Value for dirichlet boundary condition at position globalPos.
/*! In case of a dirichlet BC for the pressure equation the
pressure \f$ [Pa] \f$, and for
* the transport equation the saturation [-] have to be
defined on boundaries.
*
* \param values Values of primary variables at the boundary
* \param intersection The boundary intersection
*
* Alternatively, the function
dirichletAtPos(PrimaryVariables &values, const GlobalPosition&
globalPos)
* could be defined, where globalPos is the vector including
the global coordinates of the finite volume.
*/
void dirichletAtPos(PrimaryVariables &values, const
GlobalPosition& globalPos) const
{
if ((globalPos[0] > 600- eps_) && (globalPos[1] > 600-
eps_) )
{
values[pwIdx] = 3.5e7;
values[swIdx] = 0.0;
}
else if ((globalPos[0] < 20- eps_) && (globalPos[1] <
20- eps_) )
values[swIdx] = 1.0;
}
//! Value for neumann boundary condition \f$ [\frac{kg}{m^3
\cdot s}] \f$ at position globalPos.
/*! In case of a neumann boundary condition, the flux of matter
* is returned as a vector.
*
* \param values Boundary flux values for the different phases
* \param globalPos The position of the center of the finite
volume
*
* Alternatively, the function neumann(PrimaryVariables
&values, const Intersection& intersection) could be defined,
* where intersection is the boundary intersection.
*/
void neumannAtPos(PrimaryVariables &values, const
GlobalPosition& globalPos) const
/*@\label{tutorial-decoupled:neumann}@*/
{
values = 0;
if ((globalPos[0] < 20- eps_) && (globalPos[1] < 20-
eps_) )
{
values[nPhaseIdx] = -1e-8;
values[wPhaseIdx] = 0.0;
}
else
{
values[nPhaseIdx] = 0.0;
values[wPhaseIdx] = 0.0;
}
}
//! Initial condition at position globalPos.
/*! Only initial values for saturation have to be given!
*
* \param values Values of primary variables
* \param element The finite volume element
*
* Alternatively, the function initialAtPos(PrimaryVariables
&values, const GlobalPosition& globalPos)
* could be defined, where globalPos is the vector including
the global coordinates of the finite volume.
*/
void initial(PrimaryVariables &values,
const Element &element) const
/*@\label{tutorial-decoupled:initial}@*/
{
values = 0;
}
There isn't a problem in compilation, but in execution, io have
the folowing error:
./tutorial_decoupled
Wherever he saw a hole he always wanted to know the depth of it.
To him this was important.
- Jules Verne, A journey to the center of the earth
Rank 0: No parameter file given. Defaulting to
'./tutorial_decoupled.input' for input file.
Initializing problem 'tutorial_decoupled'
Dune reported error: ISTLError
[apply:/home/latifa/Dumux_2.6.0/dune-istl-2.3.1/dune/istl/solvers.hh:679]:
h=0 in BiCGSTAB
So please, where os the problem in my definition of the boundary
conditions? An how i can arrange it?
Best regards.
_______________________________________________
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 <tel:%28%2B49%29%200711%2F%20685-69159>
Fax:(+49) 0711/ 685-60430 <tel:%28%2B49%29%200711%2F%20685-60430>
E-Mail:[email protected]
<mailto:[email protected]>
_______________________________________________
Dumux mailing list
[email protected]
<mailto:[email protected]>
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
