Thank-you, the way you get the outflow value is smarter than the mine.
The problem I had is that I wanted to know if the scheme preserved the
mass, so I calculated the inflow and the outflow and compared whether
both of them had the same values.
I didn't know whether I calculated the inflow and outflow in the correct
way since I didn't obtain the same value. Because of that, I put my
question on the list.
The problem was that I calculated the flows just after one time
iteration and I didn't wait for the convergence to be obtained.
Thanks again.
Best regards.
Wolfgang Bangerth escribió:
So, in "valueOUT" variable, I have the outflow (m³/s), don't I?
I suppose that you ask this question because this doesn't work. Or does it?
Why don't you describe what is happening?
BTW, The way I would write this would be as follows:
FEFaceValues fe_values_total;
FEValuesViews::Extractors velocity(0);
std::vector<Tensor<1,dim> > velocity_values (n_q_points);
for (cell=...)
for (face=0...)
{
fe_values_total[velocity].get_function_values (total_solution,
velocity_values);
for (q_point=0...)
outflow += velocity_values[q_point] *
fe_values_total.normal_vector(q_point);
}
W.
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
