CFD Online Logo CFD Online URL
Home > Forums > FLUENT

Simple gas dispersion help

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

LinkBack Thread Tools Display Modes
Old   June 23, 2013, 10:56
Exclamation Simple gas dispersion help
New Member
Join Date: May 2013
Posts: 5
Rep Power: 6
masands is on a distinguished road
Hi all,

I am relatively new to CFD. I'm trying to validate a simple gas dispersion model in the presence of an obstacle. It is Thorney Island trial number 26.

It is a transient based simulation. The gas is freon 12 (31.6% w/w) and nitrogen (68.4% w/w) and is released instantaneously 50m upwind from a small building. The gas is denser than air with relative density of 2.0.

So I have enabled the species model and K-e standard model. I defined a log profile for the wind (Velocity Inlet) and ran a steady state with air only and converged. The wind speed is 1.9 m/s.

The side and up boundaries are Symmetry and the outlet is set as Outflow.

Then I patched the species to a cylindrical region, 50m upwind of the obstacle. I enabled gravity (-9.81 m/s^2) and enabled operating density (default). I ran the transient for 250 seconds with 1 second time step. I measured the concentration at the building (mol fraction of freon-12).

The shape of the graph looks exactly like the experimental which is good but unfortunately the %v/v of Freon-12 is exactly 10 times less than the experimental.

For example, the experimental peak is 2% at 50 s. Mine is 0.18% at 50 seconds.

I can't understand what's wrong.

Do I need to define any more boundary conditions?


It appears that I have now discovered the source of my problem. The mass fraction of Freon I was patching was too low. I know that the total volume of gas released is 2000 m^3. Therefore based on the fact:

Density Freon = 1305.8 kg/m3
Density N2 = 1.135 kg/m3
Total volume = 2000 m3
Required relative density = 2


1.135 (1-x) + 1305.8x = 2*1.225

This gives vol fraction. From there I can calculate the total mass and hence mass fraction.

Although before I start the simulation, the highest density in the domain is 2 kg/m3 which is still less than the relative. How is fluent calculating the density?

Last edited by masands; June 23, 2013 at 21:48.
masands is offline   Reply With Quote


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
Gas dispersion in air error message redhoax OpenFOAM Running, Solving & CFD 27 October 25, 2013 05:31
Inlet boundary condition for gas dispersion in low Re channel flow vainilreb OpenFOAM Pre-Processing 0 February 22, 2013 04:03
Gas Dispersion Simulation Fluid Engineer FLUENT 2 May 15, 2010 06:18
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
Two-Phase Buoyant Flow Issue Miguel Baritto CFX 4 August 31, 2006 12:02

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