On 8/26/20 7:05 AM, 'simon doersam' via deal.II User Group wrote:
In Step.21 (The two phase flow problem) I replaced the domain "Hyper_cube"
with "Hyper_ball".
The numerical solution is unrealistic and appears to be numerical problems at
all points,
where the boundary intersects with one of the following planes: E1: = {x \ in
R ^ 3, x1 = x2}, E2: = {x \ in R ^ 3, x1 = x3}, E3: = {x \ in R ^ 3, x2 = x3}.
The pressure and the viscosity are partially negative at these points.
Can someone help me and explain why this is happening?
You also need to change the definition of boundary values:
template <int dim>
class SaturationBoundaryValues : public Function<dim>
{
public:
SaturationBoundaryValues()
: Function<dim>(1)
{}
virtual double value(const Point<dim> &p,
const unsigned int /*component*/ = 0) const override
{
if (p[0] == 0)
return 1;
else
return 0;
}
};
This simply doesn't make any sense with the geometry you are now using.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/72cd864c-5a53-7fff-a7c0-9d9eb6b97b66%40colostate.edu.
