CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

Strange multicomponent source behaviour

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 12, 2007, 02:38
Default Strange multicomponent source behaviour
  #1
Zitron
Guest
 
Posts: n/a
Hi,

I am doing a multicomponent simulation of a room, in which air is slowly being converted into CO2 at a boundary wall. To simulate this, I made a wall with a mass fraction wall source that introduces CO2 at 1.1e-5 kg/s.

I have built a transient simulation at room temperature with buoyancy, which worked just fine converging to 0.0005, except:

1. For some reason CO2 comes in at a very low temperature, a couple of the cells next to the wall can be 60K below the specified wall temperature of 310K.

2. It seem like more CO2 is being added than I specified. Using the

volumeInt(CO2 Ideal Gas.Density * CO2 Ideal Gas.Mass Fraction)@Air

function, after 1 second, the total CO2 mass in the domain is ~1.5e-5 kg (should be 1.1e-5), after 600 seconds, it becomes 0.0099 kg instead of the predicted 0.0066 kg.

I've just got into multicomponent modeling, does any one know what might be going on?

Thanks,

-Z-
  Reply With Quote

Old   July 12, 2007, 09:43
Default Re: Strange multicomponent source behaviour
  #2
opaque
Guest
 
Posts: n/a
Dear Zitron,

How did you setup the boundary source? Is it a Continuity equation Mass Source? or just a Mass Fraction equation source.. They are different animals and you must understand what they do. At what temperature is the CO2 coming in?

Is it laminar or turbulent flow?

Hope it helps

Opaque

  Reply With Quote

Old   July 12, 2007, 12:19
Default Re: Strange multicomponent source behaviour
  #3
Zitron
Guest
 
Posts: n/a
Hi Opaque,

Thanks for your reply. The source on the wall is a mass fraction source, so I think there is no actual fluid introduced into the domain (no continuity source). There are also no inlets or outlets anywhere to add or remove fluid. The setup is similar to the 2D buoyant cavity tutorial, except with multicomponent fluid.

I don't know if there is a way to set the CO2 temperature, since on the mass fraction source options there is only a input for total source or flux of source. I can add an energy source, which will raise the temperature, but I can't understand why I would need to do that.

Since the flow speeds in the domain are below 0.6m/s, I chose not to use any turbulence model. I think that is a good assumption, correct? I also used 5-8 layers of inflated mesh on all walls.

By the way, I think I solved my second problem:

"volumeInt(CO2 Ideal Gas.Density * CO2 Ideal Gas.Mass Fraction)@Air" is the wrong way to do it, I should use:

"volumeInt(Density * CO2 Ideal Gas.Mass Fraction)@Air" instead. Now It gives me 0.00659 kg of CO2 after 600 sec, which is pretty close to the 0.0066 kg theoretical.

Thanks for your help,

-Z-
  Reply With Quote

Old   July 12, 2007, 14:01
Default Re: Strange multicomponent source behaviour
  #4
opaque
Guest
 
Posts: n/a
Dear Zitron,

If that is the case, then you need to account for the energy balance as well. CO2 does not show at wall for free. The control volume at the wall has a finite amount of energy before the conversion. That is, summation of Mass Fraction_i * energy _i = Energy_Value. When you change the composition distribution (mass fractions), the solver will determine a which temperature the new composition will return the same Energy_Value.. CO2 has a non-trivial enthalpy of formation; therefore, the temperature must be very low such that the energy_value constrain is satisfied.

Hopefully I explained correctly (and clear).. You need an energy source/sink.

Hope this helps,

Opaque.

  Reply With Quote

Old   July 12, 2007, 15:32
Default Re: Strange multicomponent source behaviour
  #5
Zitron
Guest
 
Posts: n/a
Ah! OK now it makes sense! So I need to add a energy source equal to the enthalpy of formation of CO2 x amount of CO2 produced, or perhaps I can simply set the enthalpy of formation of CO2 to zero, since I'm not interested in the chemistry.

Thanks a lot, that was really helpful!

-Z-
  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
swak4foam building problem GGerber OpenFOAM Installation 54 April 24, 2015 16:02
pisoFoam compiling error with OF 1.7.1 on MAC OSX Greg Givogue OpenFOAM Programming & Development 3 March 4, 2011 18:18
OpenFOAM on MinGW crosscompiler hosted on Linux allenzhao OpenFOAM Installation 127 January 30, 2009 20:08
DxFoam reader update hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51


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