Dear all,
I am trying to replicate the tutorial for vectorial helmholtz, but with no luck. I am sure that I am missing something along the way, but couldn't figure out what. Really appreciate any inputs.
Attached is my .sif file and the log file. The Energy Functional value that I obtain is: -3037342.7152417749 504503.76872787910
I really don't know what I have done wrong.... I have already redone this tutorial a couple of times, but no luck in replicating the results....
Thank you very much for any suggestions!
Vectorial Helmholtz - Bent waveguide
-
- Posts: 2312
- Joined: 25 Jan 2019, 01:28
- Antispam: Yes
Re: Vectorial Helmholtz - Bent waveguide
It is not converging.
The tutorial has BiCGStabl your sif has BiCGStab
The tutorial has BiCGStabl your sif has BiCGStab
Re: Vectorial Helmholtz - Bent waveguide
Thanks!!!!
Re: Vectorial Helmholtz - Bent waveguide
Sorry for bringing up this question again. I am still struggling to get the same results as in the tutorial.... I am sure it is a silly mistake that I cannot figure out. May someone please help me?
I am currently using the BiCGStabl with vanka preconditioner to solve the Vector Helmholtz Equation. For the Vector Helmholtz Post Process, I am using BiCGStab (as the tutorial don't mention anything about changing it). Although, I have to say, I have also tried to change it, with no luck in replicating the results of Energy Functional value.
The tutorial says: Energy Functional value: -11284.937620324963 453999.53923919413
I keep getting: Energy Functional value: -2241.5418742151874 447371.37659609719
I am sure it is something really minor, but I cannot figure out. I would very much appreciate any comments on that, please. The sif file I am currently running is attached here.
Thanks again and sorry for the inconvenience.
Regards
I am currently using the BiCGStabl with vanka preconditioner to solve the Vector Helmholtz Equation. For the Vector Helmholtz Post Process, I am using BiCGStab (as the tutorial don't mention anything about changing it). Although, I have to say, I have also tried to change it, with no luck in replicating the results of Energy Functional value.
The tutorial says: Energy Functional value: -11284.937620324963 453999.53923919413
I keep getting: Energy Functional value: -2241.5418742151874 447371.37659609719
I am sure it is something really minor, but I cannot figure out. I would very much appreciate any comments on that, please. The sif file I am currently running is attached here.
Thanks again and sorry for the inconvenience.
Regards
- Attachments
-
- case.sif
- with BiBGStabl solver
- (3.58 KiB) Downloaded 192 times
-
- Posts: 2312
- Joined: 25 Jan 2019, 01:28
- Antispam: Yes
Re: Vectorial Helmholtz - Bent waveguide
Not sure what to do with this in the tutorial since it is not an option in the GUI. meaning the inverse entry option.
- Attachments
-
- gui.png
- (135.51 KiB) Not downloaded yet
Re: Vectorial Helmholtz - Bent waveguide
Interesting. As a test, I manually added that entry (Inverse Relative Permeability = 1) to the waveguide sif and this is the error message:
Also, that entry doesn't exist in Solver.Keywords, and doesn't show up in Models Manual.
Rich.
Code: Select all
MAIN: Reading Model: case.sif
LoadInputFile: Scanning input file: case.sif
LoadInputFile: Scanning only size info
LoadInputFile: First time visiting
LoadInputFile: Reading base load of sif file
ERROR:: LoadInputFile:
ERROR:: LoadInputFile: Unknown specifier: [1]
ERROR:: LoadInputFile: In section: [material 1]
ERROR:: LoadInputFile: For property name:[inverse relative permeability]
Rich.
Re: Vectorial Helmholtz - Bent waveguide
Hi,
I think the vectorial Helmholtz solver doesn't read any keyword starting as "Inverse ...". Some earlier version perhaps did which may be the source of confusion. I believe the keywords described in Models Manual work.
-- Mika
I think the vectorial Helmholtz solver doesn't read any keyword starting as "Inverse ...". Some earlier version perhaps did which may be the source of confusion. I believe the keywords described in Models Manual work.
-- Mika
-
- Posts: 2312
- Joined: 25 Jan 2019, 01:28
- Antispam: Yes
Re: Vectorial Helmholtz - Bent waveguide
I went through the tutorial and my answer was the same as yours.
There may be some help in the test cases
https://github.com/ElmerCSC/elmerfem/tr ... zWaveguide
There may be some help in the test cases
https://github.com/ElmerCSC/elmerfem/tr ... zWaveguide
Re: Vectorial Helmholtz - Bent waveguide
To get back to the original question, here are a few comments on the sif file.
Either add 'Post File = case.vtu' in Simulation, or add the Result Output solver, but don't add both. In this case, I'd remove the Result Output solver.
To be sure that the vector helmholtz post process runs after vector helmholtz, change from
Note that these two points are minor and probably don't have an impact on your issue.
Rich
Either add 'Post File = case.vtu' in Simulation, or add the Result Output solver, but don't add both. In this case, I'd remove the Result Output solver.
To be sure that the vector helmholtz post process runs after vector helmholtz, change from
toExec Solver = Always
in the vector helmholtz post process solver section.Exec Solver = After Simulation
Note that these two points are minor and probably don't have an impact on your issue.
Rich
Re: Vectorial Helmholtz - Bent waveguide
So I downloaded your second case.sif and added the elmer mesh.* files from the waveguide tutorial folder.
I commented out the 'Post File = case.vtu' line in Simulation and ran elmersolver.
Rich.
I commented out the 'Post File = case.vtu' line in Simulation and ran elmersolver.
The log.txt file is attached.VectorHelmholtzSolver: Energy Functional value: -11285.445827401698 453999.47783037461
Rich.
- Attachments
-
- log.txt
- (6.6 KiB) Downloaded 177 times