Hi, As the 2D heat equation's results are quite good under the newmark system, I tried to implement the same equation in 1D to compare with an analytical solution. Not convinced by the solution values in 1D, I found that the initial condition (set to 1) was still 0 after assembling the system. Furthermore, even if the project solution tool gave satisfaction to affect initial values with 1, the final result kept going ambiguous.
I am almost sure to make a mistake but, does the newmark class is configured to work with 1D meshes ? Best regards, Maxime NB: Libmesh is configured with PETSc and the code comes from example 8 ------------------------------------------------------------------------- This SF.net email is sponsored by: Microsoft Defy all challenges. Microsoft(R) Visual Studio 2008. http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
