Saumil,
If you are using PETSc, you can tack on the command line
option -ksp_monitor_singular_value and check whether the linear solver
is running into trouble.
If that is indeed the case, you can try increasing the number of GMRES
restart steps or use modified Gram Schmidt.
Another ap
On Tue, Jun 28, 2011 at 8:03 AM, Saumil Patel wrote:
> Thank You, John. This makes sense. I was able to get results for Re 100.
>
> I have been trying to get results for Re = 1000. I've tried two approaches:
> 1) to multiply the diffusion term by 1/Re and 2) Change the boundary
> condition to U=
Thank You, John. This makes sense. I was able to get results for Re 100.
I have been trying to get results for Re = 1000. I've tried two approaches:
1) to multiply the diffusion term by 1/Re and 2) Change the boundary
condition to U=1000. BC's are still based on the Penalty formulation.
In bot
On Sun, Jun 26, 2011 at 1:37 PM, Saumil Patel wrote:
> Hello,
>
> I am getting familiar with the LIBMESH code and I was wondering if the FEM
> formulation for the Unsteady NSE in example 13 are in a non-dimensional
> form? If so, I would just need to multiply the diffusion terms in the Kuu
> and
Hello,
I am getting familiar with the LIBMESH code and I was wondering if the FEM
formulation for the Unsteady NSE in example 13 are in a non-dimensional
form? If so, I would just need to multiply the diffusion terms in the Kuu
and Kvv submatrices by the inverse Reynolds number factor, correct?
F