******************* This email originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list https://spam.ic.ac.uk/SpamConsole/Senders.aspx to disable email stamping for this address. ******************* Hi all, Sorry to disturb you. I have a full 3D simulation with large number of DOF, so I have to do something to reduce the computation cost. Could you please answer my questions on it? 1) The default number of iterative solver tolerance is 1e-10 for velocity components and 1e-9 for pressure components. Under the circumstance where the number of DOF is large, the tolerance is too small and the computation will be costly. So I alternatively use a larger tolerance with 1e-5 for velocity components and 1e-4 for pressure components. Is it accurate enough for simulations? Larger tolerance means less accuracy, but I have to take the computation time into consideration. So Could you please tell me is my set reasonable? 2) When the order expansion for pressure "p" is larger than 1, the CPU time will rise suddenly. I used to set the pressure space order one lower than velocity, but when the order I need is larger than 3 or 4, the cost is still unbearable. So could I set the pressure space order 3 or 4(maybe larger) lower than the velocity? Will it cause some numerical oscillations and incorrect solutions? 3) In user guide, I see that the expansion basis can be specified in detail as a combination of one-dimensional bases. So the heterogeneous polynomial order can be used to increase the resolution in a specified direction like spanwise. Similarly, could I use more quadrature points in a direction without increasing the polynomial order to make the result more accurate? I also encounter a problem when I set the expansions as follow: <EXPANSIONS> <E COMPOSITE="C[6]" BASISTYPE="Modified_A,Modified_A,Modified_A" NUMMODES="5,5,6" POINTSTYPE="GaussLobattoLegendre,GaussLobattoLegendre,GaussLobattoLegendre" NUMPOINTS="6,6,7" FIELDS="u,v,w" /> <E COMPOSITE="C[6]" BASISTYPE="Modified_A,Modified_A,Modified_A" NUMMODES="2,2,2" POINTSTYPE="GaussLobattoLegendre,GaussLobattoLegendre,GaussLobattoLegendre" NUMPOINTS="6,6,7"FIELDS="p" /> </EXPANSIONS> My mesh type is hexahedron and the simulation can't start with a error "NaN found during time integration". I don't konw how to solve it. 4) When I apply variable polynomial order in my field, the use of preconditioner may cause CFL number rise rapidly and the simulation will soon break down. Since the preconditioner is a good choice to reduce the CPU time, could you please tell me how to solve it? I will be very grateful to your reply. Many thanks, Zhaoyu Wang Shanghai Jiao Tong University