I am attempting to build a unsteady compressible Navier Stokes solver using Libmesh. I was planning to implement the characteristic-based split algorithm described in Volume 3 of Finite Element Methods (Zienkwicz). The algorithm essentially comprises of the following steps:
1) Compute momentum values (U*) at the new time step but without contribution of the pressure terms. (Predictor) 2) Update Pressure values using U* 3) Correct new momentum values by taking into account Pressure values. (Corrector) 4) Update Energy equation. I wanted a little advice on how I should structure my solver to ensure I make use of LibMesh's functionality. I would be solving a Linear System in all the steps. Do I declare multiple Transient Linear Systems each corresponding to the different steps in the algorithm, and solve them one after the other for a given time step? What would be the best way to implement a predictor-corrector routine with minimum hassle? Also, since the scheme is conditionally stable, I would need to evaluate h/(c+|u|) at each node, where h would be the length of the median from that node to the opposite side (unstructured). How can I evaluate this? Thank you for all your help Regards Amal ------------------------------------------------------------------------------ Open source business process management suite built on Java and Eclipse Turn processes into business applications with Bonita BPM Community Edition Quickly connect people, data, and systems into organized workflows Winner of BOSSIE, CODIE, OW2 and Gartner awards http://p.sf.net/sfu/Bonitasoft _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
