Hello,
i have following problem.
I have two sets of boundaries. The first, where is prescribed dirichlet and
second it shouldnt be prescribed. I have used this construction. It is coupled
system of NSE and convective diffusive equation for concentration. I need
no-slip condition for velocity on boundary 2. On boundary 1 is prescribed non
zero Dirichlet boundary condition for velocity and concentration. Should it be
solved only by introduction of separate system for concentration?
Thank You
Marek Capek
std::map<unsigned int,double> boundary_values;
std::vector<bool> component_mask(dim+2, true);
component_mask[dim] = false; //pressure
component_mask[dim+1]=true;
VectorTools::interpolate_boundary_values (dof_handler1,
//vtok
1,
BoundaryValues<dim>(u_d,w_d),
boundary_values,
component_mask);
MatrixTools::apply_boundary_values (boundary_values,
system_matrix1,
solution1,
system_rhs1);
component_mask[dim+1]=false;
VectorTools::interpolate_boundary_values (dof_handler1,
//gama, kde to ma reagovat
2,
BoundaryValues2<dim>(),
boundary_values,
component_mask);
MatrixTools::apply_boundary_values (boundary_values,
system_matrix1,
solution1,
system_rhs1);
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
