|March 16, 2015, 12:32||
Implementing boundary condition as additional term in differential equations
Manuel do Nascimento
Join Date: May 2012
Posts: 14Rep Power: 6
I am wondering if there is anyone who did something like this (sorry for bad english):
let's say we have a one dimensional heat convection problem dT/dt = a*d²T/dx².
Now let's say one boundary is supposed to be adiabatic -> dT/dx = 0.
The idea: I would like to implement the boundary condition as an additional term in the differential equation with the following method: dT/dt = a*d²T/dx² + A*(dT/dx)
Now the factor A has to have a value equal to zero (or extremely small) in any area, except at (or very near) the adiabatic wall where it has to have a large value. The large value of A then has the function to make sure the boundary condition is satisfied at this area to make the product in the new differential equation equal to zero.
This is just an imaginary example to get the idea. I think this method already exists but I don't know the name or how to find it in a mathematical publication.
Does anyone know a solver / a case where something like this has been done? Or an idea how to implement this in OpenFOAM? I am quite new to this software.
Best regards, and thanks
|March 25, 2015, 12:48||
Join Date: Jan 2012
Posts: 41Rep Power: 7
I think what you want to do is close to a penalization method, which is part of the family of immersed boundary methods or fictitious domain method. (You will find a lot of article for those methods in literature)
But first of all, why don't you want to use a normal boundary condition for your temperature gradient ?
|Thread||Thread Starter||Forum||Replies||Last Post|
|Question about heat transfer coefficient setting for CFX||Anna Tian||CFX||1||June 16, 2013 06:28|
|Low Mixing time Problem||Mavier||CFX||5||April 29, 2013 00:00|
|Velocity profile boundary condition||Tuca||FLOW-3D||1||April 23, 2013 12:02|
|Problems implementing new Boundary Condition||diegosene||OpenFOAM Running, Solving & CFD||5||October 31, 2011 22:32|
|asking for Boundary condition in FLUENT||Destry||FLUENT||0||July 27, 2010 00:55|