[SOLVED] Comparing with Comsol: many questions for you :D

Numerical methods and mathematical models of Elmer
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by raback »

Hi

Some more information would be needed. I think that indeed Comsol uses a 2nd order elements even though when saving the .mphtxt file it looks like a linear mesh. Are we sure that it stops there. Basically when using p-elements the element order could be further increased.

If the polynomial order is the same there can be some additional numerical tricks. In Elmer we have stabilization with various methods and experimental flux corrected transport. Comsol has also something, perhaps better tricks. Maybe you can share your findings in more detail and dig into the Comsol settings what could the settings be.

-Peter
mzenker
Posts: 1999
Joined: 07 Dec 2009, 11:49
Location: Germany

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by mzenker »

Hi,

we had a COMSOL demo version to test some time ago. AFAIR Comsol can generate a report with all the settings listed. So that could be of some help here.
Just a suggestion.

Matthias
JamesMarine
Posts: 28
Joined: 10 Jul 2019, 17:12
Antispam: Yes

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by JamesMarine »

I tried to re-run this case and learn to use Elmer for acoustic cases, but when runing the solver, the below errors occurs:

SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: helmholtz equation

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0 ffffffffffffffff

Does anyone know what is the reason for this?
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by kevinarden »

It is the error that occurs when you point to the wrong body. I downloaded your mesh and sif and the mesh was body 1 not body 2

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

Changed to
Target Bodies(1) = 1

and it ran.
JamesMarine
Posts: 28
Joined: 10 Jul 2019, 17:12
Antispam: Yes

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by JamesMarine »

Thanks, kevinarden! It works.
JamesMarine
Posts: 28
Joined: 10 Jul 2019, 17:12
Antispam: Yes

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by JamesMarine »

I continue to recieve the message that "WARNING:: IterSolve: Numerical Error: Too many iterations were needed. ", is there anything to worry about this?
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by kevinarden »

It means it has not converged to the criteria specified. It may be a good solution, but the math is not converged to the tolerances requested. It could also be a solution but not a correct one.
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by kevinarden »

you specified

Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-10

which is the default. You con increase the iterations and/or decrease the convergence tolerance.

Convergence can also be improved by a finer mesh and/or higher order elements.
panosvar
Posts: 74
Joined: 29 Oct 2019, 16:02
Antispam: Yes

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by panosvar »

It takes too much time to converge or is it just me_
mika
Posts: 230
Joined: 15 Sep 2009, 07:44

Re: [SOLVED] Comparing with Comsol: many questions for you :D

Post by mika »

A key factor to obtain fast convergence is the preconditioning strategy. Do you see any improvement if you set for example

Linear System Preconditioning = ILUn

with n in the keyword value replaced by a larger integer (i.e., ILU0, ILU1, ILU2, ...)?

-- Mika
Post Reply