|
[Sponsors] |
November 4, 2013, 10:08 |
Electrostatic Multiregion Foam
|
#1 |
New Member
Join Date: Apr 2013
Posts: 24
Rep Power: 13 |
Hi,
I am trying to develop a solver to calculate electrostatic fields in a simulation volume containing multiple dielelectric materials with different properties, e.g. conductivity and permittivity. In this case (let's assume absence of charges) I need to solve the Poisson Equation: div (epsilon grad(phi)) = 0 (with phi being the electric potential).. As far as I have understood my problem so far (), I need to define boundary conditions in which I put a constant potential on the (lets assume two) electrodes in my geometry. So far so good (nothing which couldn't be solved with electrostaticFoam). If I want to introduce several regions with different dielectric properties my field lines are not distributed equally in the volume, as they are "bend" at the interface between the two materials according to the following rules:
What I do not understand so far is how they implement the coupling between the two regions. From looking at the tutorial case, the coupling is technically implemented as a boundary condition at the region interface (which makes sense). However, from looking at the source file of the boundary condition (two be honest: I have never looked at the source file of a BC before so I am note quite sure what to look for... actually I have the impression the major work is done by a function called "updateCoeffs()"—am I right?) I cannot figure out how they actually achieve the desired coupling. I looked at turbulentTemperatureCoupledBaffleMixedFvPatchScala rField (quite a mouthful), as this is the one used in the tutorial cases of the solver. However: This seems to be a rather complicated BC as it is neither a pure FixedValue nor a FixedGradient BC, but something in between (I don't actually care, because I need FixedGradient anyway )... An intensive search revealed that "in the old days" this was easier, because people achieved this before by applying FixedGradient Boundary Conditions which do not exist anymore. (Or which I could not find).. Coming to the question I want to ask: Can someone give me a hint for a coupling condition between two regions which is more or less achieving what I want—so that I could have a look inside and understand how they do it? Or—if this is impossible, because such BCs do not exist anymore—I would be very grateful if you could give me some advice of how turbulentTemperatureCoupledBaffleMixedFvPatchScala rField actually works.. Thank you so much for reading thru all of this and for answering in advance! Best wishes Georg |
|
May 24, 2016, 11:52 |
|
#2 |
New Member
dklogo
Join Date: May 2015
Posts: 11
Rep Power: 11 |
Hi GPesch,
I am facing a similar problem with you. Actually I want to implement a coupling condition for electric field betwen the fluid and solid region. If you solve your problem; can you suggest me a solution ? How I can create a coupling condition ? |
|
May 25, 2016, 01:17 |
|
#3 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 |
As the original poster said, have a look at chtMultiRegionFoam. You will find all the building blocks you need.
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. |
|
November 14, 2021, 07:36 |
|
#4 |
New Member
Jelmer
Join Date: Sep 2012
Posts: 3
Rep Power: 13 |
Dusting off this old topic here...
I'm facing a similar problem. I want to calculate a steady electric field in the absence of any charges in multiple materials. This is needed for the calculation of stray currents in seawater and the seabed from a voltage source. What I'm currently doing, is using chtMultiRegionSimpleFoam, and just pretend I'm doing a heat conduction problem instead of an electrical conduction problem. Instead of fixed voltage boundary conditions, I now have fixed temperture boundary conditions. Both problems are governed by the Laplace equation, so results should be identical, right? Curious to hear whether this approach is fully valid...! |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] mesh airfoil NACA0012 | anand_30 | OpenFOAM Meshing & Mesh Conversion | 13 | March 7, 2022 17:22 |
[blockMesh] error message with modeling a cube with a hold at the center | hsingtzu | OpenFOAM Meshing & Mesh Conversion | 2 | March 14, 2012 09:56 |
[blockMesh] BlockMesh FOAM warning | gaottino | OpenFOAM Meshing & Mesh Conversion | 7 | July 19, 2010 14:11 |
[blockMesh] Axisymmetrical mesh | Rasmus Gjesing (Gjesing) | OpenFOAM Meshing & Mesh Conversion | 10 | April 2, 2007 14:00 |
[Gmsh] Import gmsh msh to Foam | adorean | OpenFOAM Meshing & Mesh Conversion | 24 | April 27, 2005 08:19 |