ERROR:: CRS_IncompleteLU: Memory allocation error.

Numerical methods and mathematical models of Elmer
Post Reply
studioa
Posts: 23
Joined: 23 Jul 2011, 23:39
Antispam: Yes

ERROR:: CRS_IncompleteLU: Memory allocation error.

Post by studioa »

Hello -
I'm trying Elmer on another aiflow simulation for a building interior. After I get into this issue:

ERROR:: CRS_IncompleteLU: Memory allocation error.

Running on Win 7 on a quad core machine with 6 gb of memory.

I use the version:

ELMER SOLVER (v 6.1) STARTED AT: 2011/10/11 12:09:43
MAIN:
MAIN: ==========================================
MAIN: E L M E R S O L V E R S T A R T I N G
MAIN: Library version: 6.1 (Rev: 5191)
MAIN: ==========================================


And the sif output is:

Header
CHECK KEYWORDS Warn
Mesh DB "." "."
Include Path ""
Results Directory ""
End

Simulation
Max Output Level = 4
Coordinate System = Cartesian
Coordinate Mapping(3) = 1 2 3
Simulation Type = Transient
Steady State Max Iterations = 1
Output Intervals = 10
Timestepping Method = BDF
BDF Order = 1
Timestep intervals = 5
Timestep Sizes = .5
Solver Input File = ATRIUM_0.sif
Post File = ATRIUM_0.ep
End

Constants
Gravity(4) = 0 -1 0 9.82
Stefan Boltzmann = 5.67e-08
Permittivity of Vacuum = 8.8542e-12
Boltzmann Constant = 1.3807e-23
Unit Charge = 1.602e-19
End

Body 1
Target Bodies(1) = 1
Name = "Body 1"
Equation = 1
Material = 1
End

Solver 1
Equation = Result Output
Procedure = "ResultOutputSolve" "ResultOutputSolver"
Output File Name = ATRIUM_0
Output Format = Vtu
End

Solver 2
Equation = Navier-Stokes
Procedure = "FlowSolve" "FlowSolver"
Variable = Flow Solution[Velocity:3 Pressure:1]
Exec Solver = After Timestep
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-8
Nonlinear System Max Iterations = 20
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-8
Linear System Preconditioning = ILU0
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
End

Solver 3
Equation = K-Epsilon
Procedure = "KESolver" "KESolver"
Exec Solver = After Timestep
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-8
Nonlinear System Max Iterations = 20
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-8
Linear System Preconditioning = ILU0
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
End

Equation 1
Name = "NS+KE"
Active Solvers(3) = 1 2 3
End

Material 1
Name = "AIR"
Viscosity = 1.983e-5
Heat expansion Coefficient = 3.43e-3
Heat Conductivity = 0.0257
Sound speed = 343.0
Heat Capacity = 1005.0
Density = 1.205
End

Boundary Condition 1
Target Boundaries(1) = 1
Name = "OUTLET"
End

Boundary Condition 2
Target Boundaries(1) = 2
Name = "SIDES"
Noslip wall BC = True
End

Boundary Condition 3
Target Boundaries(1) = 3
Name = "SIDES"
Noslip wall BC = True
End

Boundary Condition 4
Target Boundaries(1) = 4
Name = "SIDES"
Noslip wall BC = True
End

Boundary Condition 5
Target Boundaries(1) = 5
Name = "SIDES"
Noslip wall BC = True
End

Boundary Condition 6
Target Boundaries(1) = 6
Name = "SIDES"
Noslip wall BC = True
End

Boundary Condition 7
Target Boundaries(1) = 7
Name = "SIDES"
Noslip wall BC = True
End

Boundary Condition 8
Target Boundaries(1) = 8
Name = "SIDES"
Noslip wall BC = True
End

Boundary Condition 9
Target Boundaries(1) = 9
Name = "SIDES"
Noslip wall BC = True
End

Boundary Condition 10
Target Boundaries(1) = 10
Name = "INLET"
Velocity 1 = 10
End



I have tried incomplete lower upper 1 and 2 with BiCGStabl , but it gets to the same error.

Help & suggestions greatly appreciated.

Thanks!
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: ERROR:: CRS_IncompleteLU: Memory allocation error.

Post by raback »

Well, you might try preconditioner "none" and add "BiCGStabl Polynomial Degree = 4". So compensate the poorer preconditioner with more robust iterative method. -Peter
Post Reply