Hi,
For the pressure equation I have
-del.K(x)grad p + a(x)*p = f(x)
with neumann boudary conditions u.n=g at the left and right faces and
dirichlet everywhere else. Can I still use Step 43 as a starting point?
The advection equation is the same as in Step-43
If yes, are these the only two places where the code changes?
{
local_matrix(i,j) += (phi_u[i] * k_inverse_values[q] *
mobility_inverse(old_s,viscosity) * phi_u[j]
- div_phi_u[i] * phi_p[j]
- phi_p[i] * div_phi_u[j])
* darcy_fe_values.JxW(q);
}
and
{
const Tensor<1,dim>
<http://www.dealii.org/developer/doxygen/deal.II/classTensor.html>
phi_i_u = darcy_fe_face_values[velocities].value (i, q);
local_rhs(i) += -(phi_i_u *
darcy_fe_face_values.normal_vector(q) *
boundary_values[q] *
darcy_fe_face_values.JxW(q));
}
In other words, Can I use FE_Q elements for both velocity and pressure in
this case too?
Also, I don't get this part of the code in assemble_darcy_system:
for (unsigned int i=0; i<dofs_per_cell; ++i)
for (unsigned int j=i+1; j<dofs_per_cell; ++j)
local_matrix(i,j) = local_matrix(j,i);
--
Thanks,
gk
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
