Hi Soon, Try VectorTools:interpolate - I found that when I was applying discontinuous initial conditions this implemented them better than the "project" function.
Thanks, Katie Katie Leonard DPhil student in Computational Biology, The University of Oxford. ________________________________________ From: [email protected] [[email protected]] on behalf of Soon Hoe Lim [[email protected]] Sent: 23 June 2012 13:52 To: [email protected] Subject: [deal.II] VectorTools::project for discontinuous initial condition Hi all. For solving time dependent problems, at t=0 we need to project the initial conditions into the finite element space. But when we have a discontinuos initial condition, say u(x,y,0)=1200 on [0.4,0.6]x[0.4,0.6], u(x,y,0)=0 elsewhere, where the domain is [0,1]x[0,1], I observed that the VectorTools::project function tries to smooth the discontinuity at t=0. This causes the solution to violate the initial condition because say at (0.4, 0.4) the solution at t=0 is no longer 1200 as desired but some smaller value. In this case, how should I deal with the initial condition? Thank you for any help. Sincerely, Soon Hoe _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
