SCIRun User Mailing List

Text archives Help

[SCIRUN-USERS] mixed formulation

Chronological Thread 
  • From: Konstantin Skaburskas <>
  • To:
  • Subject: [SCIRUN-USERS] mixed formulation
  • Date: Sun, 1 Feb 2004 16:00:37 +0200 (EET)

Dear All,

I have a question to ask.

The topic I'm working on is the calculation of the power absorption due to
the induction of eddy currents in the human thorax.

[brief description:]
I reformulated my problem a bit, and re-wrote governing equations I want
to solve. Before I did it, the problem was to find certain potential field
distribution (homog. Neumann problem for Poisson's eq.), then take
gradient of it and multiply the result by conductivity of the medium to
get the flux (representing irrotational part of total currents). But now I
want to find the total flux at once. This is my mixed formulation:

 |(1/sigma)*J + grad(Phi) = f
 |div(J)                  = 0
  b.c.s. : homogeneous Neumann

to find J - flux (current density) [and additionally as a "garbage" I'll
get Phi - potential field, in which I'm not interested]. Where "sigma" is
the conductivity of the medium and "f" is a forcing function.

So, to be able to solve that, I think I'll have to use whether Mixed FEs -
the Raviart-Thomas-Ne'de'lec (3D) or divergence-free (3D) or some
LBB-stable (e.g. triquadratic Lagrangian 27-noded hexahedral) elements. Do
you have such elements in SCIRun package and would it be possible to solve
my mixed problem using SCIRun?

Thank you in advance!

Best regards,
Konstantin Skaburskas

- Konstantin Skaburskas, Ph.D. student         -
- Inst. of Experimental Physics & Technology      -
- University of Tartu                           phone:  +372 7 374843 -
- Tahe 4-234, Tartu, 51010, Estonia             fax:    +372 7 375858 -

== The SCIRun Users mailing list: send email to   ==
== for more details.                                                     ==
== Please acknowledge use of SCIRun in your papers and reports:          ==
==   see         ==

  • [SCIRUN-USERS] mixed formulation, Konstantin Skaburskas, 02/01/2004

Archive powered by MHonArc 2.6.16.

Top of page