CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Pure conduction problem with different diffusivities

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 3, 2008, 05:09
Default Hi all, I'm trying to solve
  #1
New Member
 
Diego Angeli
Join Date: Mar 2009
Posts: 26
Rep Power: 8
diego_angeli is on a distinguished road
Hi all,

I'm trying to solve a conduction problem in materials with different diffusivities. Think about a wall made of some layers of different materials.

I started from the laplacianFoam solver, and modified it as follows:

in the createFields.H file, I changed the diffusivity DT from a constant to a field. Here's the code

volScalarField DT
(
IOobject
(
"DT",
runTime.timename(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
),
mesh
);

the solver re-compiled fine, so I suppose to have performed things in the right way

then I followed the instructions I found here:
http://www.cfd-online.com/OpenFOAM_D...tml?1214910015

I tried to use setFields to give to DT the values I wanted in the different regions of my mesh (a simple blockMesh: I have to simulate a flat wall with thermal leaks), but as I create the file "DT" in the "0" folder and I run setFields, the resulting "DT" file returns me just a uniform field with DT equal to the default value I set for it in the setFieldsDict.

Have you got any clue about that? I'm rather new to OpenFOAM, so I could be stuck on a trivial thing...

And, do you think that my way of changing the solver is the right one?

Thank you very much in advance!

Diego Angeli
Modena, Italy
diego_angeli is offline   Reply With Quote

Old   October 4, 2008, 18:44
Default Hi Diego, Sounds like the r
  #2
Member
 
Ola Widlund
Join Date: Mar 2009
Location: Sweden
Posts: 87
Rep Power: 8
olwi is on a distinguished road
Hi Diego,

Sounds like the right approach.

I think you should focus on the possibility that you're doing something wrong with setFields. Try to play around with the dambreak tutorial, for example.

Another thing you need to get good results: use harmonic interpolation of DT in the laplacian. Otherwise you get some smearing on the interface between media. It can be shown formally that harmonic rather than arithmetic (=linear) interpolation is what you need for a correct discretization.

Good luck!

/Ola
olwi is offline   Reply With Quote

Old   October 8, 2008, 06:53
Default Hi Ola Thank you for the ti
  #3
New Member
 
Diego Angeli
Join Date: Mar 2009
Posts: 26
Rep Power: 8
diego_angeli is on a distinguished road
Hi Ola

Thank you for the tip on the interpolation scheme, and for the "thumbs up".

The approach was indeed right, and I just systematically mistook the order of the two vertices defining boxes in setFields. The results I got are satisfactory!

Now I'm trying to impose a mixed condition at the boundaries and I saw that it is not so immediate... let's see

by now, many thanks!!

Diego
diego_angeli 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
Need a help concerning a pure convection problem happysmile Main CFD Forum 6 July 8, 2009 13:24
pure conduction sudhir FLUENT 4 April 17, 2008 02:44
problem pure diffusion Jorn CFX 3 July 31, 2007 12:06
conduction problem erica FLUENT 0 February 14, 2006 13:27
Convection-Conduction Problem TOM Main CFD Forum 0 February 7, 2005 10:57


All times are GMT -4. The time now is 21:53.