[Solved] Electrostatic: Charge spere in vacuum

Numerical methods and mathematical models of Elmer
Wil
Posts: 72
Joined: 10 Jan 2014, 13:26
Antispam: Yes
Location: Austria

Re: [Solved] Electrostatic: Charge spere in vacuum

Post by Wil »

Dear manu1905,
I guess you fixed potential is coming from a conductive material below. Do you also consider to stimulate the conductive domain and couple the electrostatic solver with the static current conduction as discussed here viewtopic.php?f=3&t=4362&start=0
.
BR,
Wil
mickyjohn
Posts: 1
Joined: 26 Mar 2018, 12:40
Antispam: Yes

Re: [Solved] Electrostatic: Charge spere in vacuum

Post by mickyjohn »

thanks for this sharing.
sabrina
Posts: 27
Joined: 16 Sep 2019, 14:26
Antispam: Yes

Re: [Solved] Electrostatic: Charge spere in vacuum

Post by sabrina »

Hello everyuone, I know this is an old post, but I was trying to reproduce your example using your sif and your gmsh (I am new to Elmer), but I run into this error when running the solver:

ELMER SOLVER (v 8.4) STARTED AT: 2019/09/16 13:21:46
ParCommInit: Initialize #PEs: 1
MAIN:
MAIN: =============================================================
MAIN: ElmerSolver finite element software, Welcome!
MAIN: This program is free software licensed under (L)GPL
MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd.
MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi
MAIN: Version: 8.4 (Rev: 141b308, Compiled: 2019-09-13)
MAIN: Running one task without MPI parallelization.
MAIN: Running with just one thread per task.
MAIN: HYPRE library linked in.
MAIN: MUMPS library linked in.
MAIN: =============================================================
MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: case.sif
LoadInputFile: Scanning input file: case.sif
LoadInputFile: Loading input file: case.sif
WARNING:: LoadInputFile: > Material 2 < not used in any Body!
Loading user function library: [StatElecSolve]...[StatElecSolver_Init0]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init0]
LoadMesh: Base mesh name: ./point_charge

Program received signal SIGBUS: Access to an undefined portion of a memory object.

Backtrace for this error:
#0 0x10ce4df3d
#1 0x10ce4d34d
#2 0x7fff74d52b5c
#3 0x108177a04
#4 0x1081bb2ff
#5 0x1081c6936
#6 0x1081ba57e
#7 0x1084150c6
#8 0x10843a4f5
#9 0x108200f39
#10 0x1085c24e9
#11 0x1080caba7
#12 0x1080cad2e
Bus error: 10

Can anyone advise me on what is wrong?

Thank you very much in advance!

Sabrina
annier
Posts: 1168
Joined: 27 Aug 2013, 13:51
Antispam: Yes

Re: [Solved] Electrostatic: Charge spere in vacuum

Post by annier »

Dear Sabrina,
Sabrina wrote:LoadMesh: Base mesh name: ./point_charge
Is your meshfile name point_charge?

Code: Select all

Header
  CHECK KEYWORDS Warn
  Mesh DB "." "."  !The mesh files are in the same place as the sif file
  !Mesh DB "." "Mesh_1"  !The directory Mesh_1( containing mesh files) and sif file are inside the same folder
  Include Path ""
  Results Directory ""
End
If point_charge.msh is the gmsh mesh format, did you run ElmerGrid command to convert this mesh format into ElmerSolver readable mesh format (http://www.nic.funet.fi/pub/sci/physics ... Manual.pdf).

Code: Select all

 $ ElmerGrid 14 2 point_charge.msh -autoclean 
numbers :
14 - gmsh file format of input mesh
2- ElmerSolver readable mesh file format

The summary is:
ElmerSolver finds problem when loading the mesh in your case.
You can create your own geometry/ mesh , and run the test case.



Yours Sincerely,
Anil Kunwar
Anil Kunwar
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
sabrina
Posts: 27
Joined: 16 Sep 2019, 14:26
Antispam: Yes

Re: [Solved] Electrostatic: Charge spere in vacuum

Post by sabrina »

Hi Anil,

Thank you so much for your reply.

Just to make myself more clear, these are the steps I did:

1. I opened with gmsh the geo file posted here and then created the mesh and saved it in a file called point_charge.msh.

2. I run elmergrid to convert it to a format that Elmer understands. I did exactly what you wrote: Elmergrid 14 2 point_charge.msh -autoclean.

3. The new folder called point_charge was created with the corresponding mesh.[stuff] files.

4. I copied the sif file and put it in the same directory as the folder containing the mesh files (and so I changed the line in the sif so that it finds the mesh folder). Since that did not work, I also tried putting the sif file without editing it inside the folder with the mesh files. That did not work either...

I am trying other tests from elmer (like the capacitance of two balls) and the same thing happens. However, I noticed that when I run 2D examples, it works.

Any hint is very much appreciated!

Best regards,

Sabrina
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: [Solved] Electrostatic: Charge spere in vacuum

Post by kevinarden »

Posting your files would help to debug.

If the sif file is in the same file as the mesh files then the line is

Mesh DB "." "."

it means look in the current folder for the mesh files

If the mesh is in a sub folder to the sif file then it is

Mesh DB "." "point_charge"

it means look in a subfolder called point_charge to the current folder for the mesh files.

ElmerGrid 14 2 point_charge.msh -autoclean put the mesh in a subfolder to the current folder call point_charge
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: [Solved] Electrostatic: Charge spere in vacuum

Post by kevinarden »

Mesh DB "." "."

more specifically the first "." is the path, the . just means current directory
the second "." means the directory in the path . just means the current directory
if you wanted to fully specify it it could be

Mesh DB "/home/titan/fea" "point_charge"

then it would look in /home/titan/fea/point_charge for the mesh files
sabrina
Posts: 27
Joined: 16 Sep 2019, 14:26
Antispam: Yes

Re: [Solved] Electrostatic: Charge spere in vacuum

Post by sabrina »

Hi Kevin, thank you for your reply!

The paths are not the issue here, I am using them correctly according to what you explained. Also because, as I mentioned in the post, the program seems to work fine with 2D problems. It is only when I try any 3D case that I run into this issue. I did not post the sif and mesh files because they are an exact copy of the OP's.

Best regards,

Sabrina
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: [Solved] Electrostatic: Charge spere in vacuum

Post by kevinarden »

I downloaded the files an ran them with no issues.
Post Reply