CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Dirichlet–Neumann Partitioning BCs

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By Cyp

Reply
 
LinkBack Thread Tools Display Modes
Old   April 13, 2013, 08:05
Post Dirichlet–Neumann Partitioning BCs
  #1
New Member
 
Giancarlo
Join Date: Apr 2013
Location: Milan
Posts: 18
Rep Power: 4
Giancarlo_IngChimico is on a distinguished road
Hi all,
I have to implement the Dirichlet–Neumann Partitioning BCs for my solver that treats two different regions (fluid and solid).
Is there any preconfigured BC in the framework of OpenFoam?

Thanks

Regards

Giancarlo
Giancarlo_IngChimico is offline   Reply With Quote

Old   April 14, 2013, 05:14
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,592
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Giancario,

I am not familiar with the term "Dirichlet-Neumann Partitioning", however, maybe the mixedFvPatch class is what you need.

In this class you can switch (smoothly) from a dirichlet to a neumann condition depending on the instantaneous dynamics of the system. I use it in the binary fashion to switch between fixed value and zero gradient, where the choice depends on the internal value of some concentration field.

Kind regards

Niels
ngj is online now   Reply With Quote

Old   April 14, 2013, 15:45
Default
  #3
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 227
Rep Power: 9
Cyp is on a distinguished road
Hi Giancarlo,

I fear there is no Direchlet-Neumann BC coded by default in OF and that you have to develop your owns. You can find a lot of informations in http://www.personal.psu.edu/dab143/O...ven_slides.pdf and its associated tutorial http://www.personal.psu.edu/dab143/O...erTraining.tgz.

The relaxation factors are quite important in this example.

@++
Cyp
wayne14 likes this.
Cyp is offline   Reply With Quote

Old   April 14, 2013, 17:20
Default
  #4
New Member
 
Giancarlo
Join Date: Apr 2013
Location: Milan
Posts: 18
Rep Power: 4
Giancarlo_IngChimico is on a distinguished road
Thanks a lot Cyp. This materials is fantastic. It's just what I was looking for.

Regards

Giancarlo
Giancarlo_IngChimico is offline   Reply With Quote

Old   May 30, 2014, 08:48
Default
  #5
Senior Member
 
Join Date: Oct 2013
Posts: 229
Rep Power: 4
chriss85 is on a distinguished road
In case you still read here, have you got this working?

I'm currently implementing something similar for electric potential between different regions.

By now OpenFOAM has a boundary condition for the temperature that uses an approach like this, called turbulentTemperatureCoupledBaffleMixed.

I adopted it to my needs, but I don't understand the weighting factor in the mixed formulation. Why is it set to this value? Does it help convergence if this value is used instead of some arbitrary value between 0 and 1 (like 0.5)?

Furthermore, how does relaxation help here? My current solver doesn't use relaxation (it's transient), and it requires quite some iterations initially until a consistent solution is reached (up to 30 in a test case I made).

The link mentioned using Aitken's delta˛ method (see here), which is basically a quadratic extrapolation of the convergence series of the field. In its current implementation, OpenFOAM doesn't use this. Have you implemented this? Is it worth it for transient cases? I believe that the convergence should be better in transient cases because the solution changes only a bit between steps.
chriss85 is online now   Reply With Quote

Old   June 2, 2014, 09:46
Default
  #6
Senior Member
 
Join Date: Oct 2013
Posts: 229
Rep Power: 4
chriss85 is on a distinguished road
After some more investigation I figured out that relaxation didn't refer to the field relaxation, but rather to the coupling between the patches of different regions.

I also compiled the MRConjugateHeatFoam solver with OF2.3, but it doesn't calculate any sensible results and the solution of the example looks very incorrect. Has anyone tried this solver with success?
chriss85 is online now   Reply With Quote

Old   June 9, 2014, 10:24
Default MRConjugateHeatFoam
  #7
Senior Member
 
Mohammad Shakil Ahmmed
Join Date: Oct 2012
Location: AU
Posts: 122
Rep Power: 5
ahmmedshakil is on a distinguished road
Hi Chirss85,
I worked with the MRConjugateHeatFoam. In my case, I validated the solver with some experimental cases, and the results were really good. BTW, my problem was heat conduction between two solids (i.e. silicon on a glass substrates) due to laser annealing.

Cheers
shakil


Quote:
Originally Posted by chriss85 View Post
After some more investigation I figured out that relaxation didn't refer to the field relaxation, but rather to the coupling between the patches of different regions.

I also compiled the MRConjugateHeatFoam solver with OF2.3, but it doesn't calculate any sensible results and the solution of the example looks very incorrect. Has anyone tried this solver with success?
ahmmedshakil is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
ERROR #001100279 has occurred in subroutine ErrAction. smnaryal CFX 7 January 3, 2013 19:09
Dealing with BC's in OF 1.6 vkrastev OpenFOAM Running, Solving & CFD 5 September 4, 2012 11:58
Understanding Code behind BCs Linse OpenFOAM Programming & Development 8 January 9, 2012 09:58
BCs for Pressure Correction Equation (SIMPLE) Bharath Somayaji Main CFD Forum 1 March 1, 2006 07:12
Need Assistance: BCs for k-e model implementation Dean Schrage Main CFD Forum 3 August 11, 2001 01:02


All times are GMT -4. The time now is 03:37.