I'm working with Simvascular for steady flow but after many times trying I couldn't get a good nonlinear residual number, every time I got higher than 1*10-3 as in the guide is mentioned that it's not a good sign.
I put here the details of my simulation:
Unit in MGS
Mesh
-
Number of Nodes: 182073
Number of Elems: 1069575
Number of Edges: 89568
Number of Faces: 59712
-
Inlet and Outlet BCs
-
BC Type: Prescribed Velocities
Analytic shape: parabolic
Point Number: 201
Fourier Modes: 10
Flow rate(from File):
0.0 -507163.7
1.0 -507163.7
Period: 1
Outlet, Resistance: 0.02366
-
Wall Properties
Type: Rigid
-
Solver Parameters
-
Time Step Parameters
-
Number of Timesteps: 1000
Time Step Size: 0.01
-
Output Control
-
Number of Timesteps between Restarts: 20
Output Surface Stress: True
Force Calculation Method: Velocity Based
Print Average Solution: True
Print Error Indicators: False
-
Step Construction: 15
-
Advanced Parameters
Pressure Coupling: Implicit
Backflow Stabilization Coefficient: 0.2
-
Non-linear Iteration Control
-
Residual Control: True
Residual Criteria: 0.01
Minimum Required Iterations: 3
-
Linear Solver
-
svLS Type: NS
Number of Krylov Vectors per GMRES Sweep: 100
Number of Solves per Left-hand-side Formation: 1
Tolerance on Momentum Equations: 0.005
Tolerance on Continuity Equations: 0.04
Tolerance on svLS NS Solver: 0.04
Maximum Number of Iterations for svLS NS Solver: 100
Maximum Number of Iterations for svLS Momentum Loop: 20
Maximum Number of Iterations for svLS Continuity Loop: 4000
-
Discretization Control
-
Time Integration Rule: Second order
Time Integration Rho Infinity: 0.5
Flow Advection Form: Convective
Quadrature Rule on Interior: 2
Quadrature Rule on Boundary: 3
-
Actually, I got good residual when I decreased the inflow to 10 [mm3/s] but I'd like to know that is it high inflow(507163.7[mm3/s]) for Simvascular? what is the maximum inflow that I can use in Simvascular?
I tried to change my solver parameters to have a good result but I wonder if you could possibly explain me to how can I deal with this problem?
Do you have any idea about good solver parameters?
Thank you in advance.
Kind regards,
Muhamad