<div dir="ltr"><br clear="all"><div>Hi Folks,</div><div><br></div><div>I am thinking of trying TS module of PETSc to free myself of maintaining </div><div>time stepping part of the code and to try different </div><div>time integrators easily. The problem that I am considering is incompressible </div><div>Navier Stokes with variable coefficients solved using  projection method. </div><div>The method would involve solving equations in two parts:</div><div><br></div><div>rho u*_t - div ( sigma) = RHS( u_explicit)  ---- (1)</div><div><br></div><div>and then projecting u* onto divergence free u by solving a pressure </div><div>Poisson equation  (Eq 2).</div><div><br></div><div> Here sigma = mu ( u_ij + u_ji) is the viscous stress</div><div>tensor.</div><div><br></div><div>I have custom solvers for solving (1) and (2) using geometric MG  wrapped as</div><div>a PCSHELL to MATSHELL obtained by discretizing (1) and (2) say by </div><div>Backward Euler.</div><div><br></div><div>Is it possible to wrap Eq (1) and (2) in TS in matrix-free form and then solve for </div><div>velocity and pressure using projection method using matrix-free solvers? </div><div><br></div><div>Thanks,</div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>--Amneet <br><br></div><div><br></div><div><br></div></div></div>
</div>