Skip to content
Snippets Groups Projects
Commit 676dfba2 authored by Nicolas Marsic's avatar Nicolas Marsic
Browse files

more with markdown

parent 2225a8f8
No related branches found
No related tags found
No related merge requests found
...@@ -98,7 +98,7 @@ Python module: ...@@ -98,7 +98,7 @@ Python module:
getdp.py getdp.py
Finally, Don't forget to update your `PATH`, `PYTHONPATH` and `LD_LIBRARY_PATH`. Finally, don't forget to update your `PATH`, `PYTHONPATH` and `LD_LIBRARY_PATH`.
Enjoy! Enjoy!
......
This is a simple linear eigenvalue problem example. It consists in a rectangular
electromagnetic cavity. Simply run:
make cim
to test cim.py.
The [Makefile](Makefile) will first mesh the geometry [square.geo](square.geo)
by calling [Gmsh](http://gmsh.info).
Afterwards, cim.py is called. The [GetDP](http://getdp.info) formulation is
located in [maxwell.pro](maxwell.pro). In the [GetDP](http://getdp.info)
formulation, the variables:
angularFreqRe
angularFreqIm
x()
b()
imposeRHS
doPostpro
doApply
fileName
are reserved for cim.py.
In order to check the solution, simply run:
make ref
This will solve the linear eigenvalue problem with a classical algorithm. This
resolution is implemented in [ref.pro](ref.pro). If you use [GetDP](
http://getdp.info) with [SLEPc](http://slepc.upv.es), the default algorithm
should be [Krylov-Schur](https://dx.doi.org/10.1137/S0895479800371529).
This is a simple linear eigenvalue problem example. It consists in a rectangular
electromagnetic cavity. Simply run:
make cim
to test cim.py.
The Makefile will first mesh the geometry square.geo by calling Gmsh.
Afterwards, cim.py is called. The GetDP formulation is located in maxwell.pro.
In the GetDP formulation, the variables:
angularFreqRe
angularFreqIm
x()
b()
imposeRHS
doPostpro
doApply
fileName
are reserved for cim.py.
In order to check the solution, simply run:
make ref
This will solve the linear eigenvalue problem with a classical algorithm.
This resolution is implemented in ref.pro. If you use GetDP with SLEPc,
the default algorithm should be Krylov-Schur.
This is a non-linear eigenvalue problem example. It consists in a spherical This is a non-linear eigenvalue problem example. It consists in a spherical
electromagnetic cavity. The wall conductivity is modelled by a Leontovich electromagnetic cavity. The wall conductivity is modelled by a Leontovich
surface impedance boundary condition (SIBC). Simply run: surface impedance boundary condition (SIBC). Simply run:
make
make
to test cim.py. to test cim.py.
The Makefile will first mesh the geometry square.geo by calling Gmsh. The [Makefile](Makefile) will first mesh the geometry [sphere.geo](sphere.geo)
Afterwards, cim.py is called. The GetDP formulation is located in sphere.pro. by calling [Gmsh](http://gmsh.info). Afterwards, cim.py is called. The [GetDP](
In the GetDP formulation, the variables: http://getdp.info) formulation is located in [sphere.pro](sphere.pro). In the
angularFreqRe [GetDP](http://getdp.info) formulation, the variables:
angularFreqIm
x() angularFreqRe
b() angularFreqIm
imposeRHS x()
doPostpro b()
doApply imposeRHS
fileName doPostpro
are reserved for cim.py. The parts of the GetDP code related to cim.py are doApply
located in cimParameters.pro and cimResolution.pro. fileName
are reserved for cim.py. The parts of the [GetDP](http://getdp.info) code
related to cim.py are located in [cimParameters.pro](cimParameters.pro) and
[cimResolution.pro](cimResolution.pro).
For the default parameters: For the default parameters:
radius = 150e-6 m
conductivity = 1e15 S/m radius = 100 mm
conductivity = 1e15 S/m
the analytical resonance angular frequency for the fundamental mode should be: the analytical resonance angular frequency for the fundamental mode should be:
5.48362e12+4.24068e05j
8.22543e9+2.46361e1j
where j is the imaginary unit. This analytical result comes from the reference: where j is the imaginary unit. This analytical result comes from the reference:
S. Papantonis and S. Lucyszyn, "Lossy spherical cavity resonators for [S. Papantonis and S. Lucyszyn, "Lossy spherical cavity resonators for
stress-testing arbitrary 3D eigenmode solvers," Progress In Electromagnetics stress-testing arbitrary 3D eigenmode solvers," Progress In Electromagnetics
Research, vol. 151, pp. 151-167, 2015. Research, vol. 151, pp. 151-167, 2015.](https://dx.doi.org/10.2528/PIER15031702)
Same example as [maxwell\_sibc](../maxwell_sibc), but imposes a divergence free
source for Beyn. This is achieved by taking the curl of the random Beyn source.
This example requires a [GetDP](http://getdp.info) version compiled with [Gmsh](
http://gmsh.info) to handle Fields.
See as maxwell_sibc, but imposes a divergence free source for Beyn.
This is achieved by taking the curl of the random Beyn source.
This example requires a GetDP version compiled with Gmsh to handle Fields.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment