Marek, > 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.
Your construction looks basically alright (though you don't have to call apply_boundary_values in the middle: it is enough if you simply let the two calls to interpolate_boundary_values add new entries to the boundary_values map). As you mention in your second mail, since the concentration should be left untouched on boundary 2, you need to set the boundary_mask to false for this component -- but you already do that. In your mail, you don't ask a question: you don't say whether this code works as you'd expect or not. What do you need to know? Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
