
[Sponsors] 
Diffusive transport equation for additional variable 

LinkBack  Thread Tools  Search this Thread  Display Modes 
April 21, 2018, 14:45 
Diffusive transport equation for additional variable

#1 
New Member
Swapnil Rane
Join Date: Sep 2017
Posts: 23
Rep Power: 8 
I'm trying to solve a porous media problem with the 'Diffusive transport equation for additional variable' equation.
But I'm facing trouble assigning diffusion value and source value for the model Can anyone help me? 

April 22, 2018, 06:06 

#2 
Senior Member
GertJan
Join Date: Oct 2012
Location: Europe
Posts: 1,827
Rep Power: 27 
What troubles? Need more info.
Don't you know how to do it? You tried it but CFX gives an error? You get unrealistic results? Divergence? What is your problem? 

April 22, 2018, 06:10 
Boundary conditions

#3 
New Member
Swapnil Rane
Join Date: Sep 2017
Posts: 23
Rep Power: 8 
How can I set up the source term to be zero value according to the equation?
What’s the value of additional variable i.e diffusion at the outlets? 

April 22, 2018, 07:21 

#4 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,695
Rep Power: 143 
If you do not add a source term then the source term is zero. Is that what you mean?
When you define the diffusion coefficient that diffusion value is used everywhere. There is nothing special at outlets (except that values are convected from the upstream control volume, as normal for outlets).
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

April 22, 2018, 07:29 
Change the diffusion value

#5 
New Member
Swapnil Rane
Join Date: Sep 2017
Posts: 23
Rep Power: 8 
Yes. The source term to be zero.
Also, even though when I change the diffusion coefficient by any factor, the mass flow at outlet remains same. How can I fix this? 

April 22, 2018, 07:41 

#6 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,695
Rep Power: 143 
If you want the source term to be zero then don't add any source terms. Simple.
Mass flow at the outlet is a function of the pressure difference over your modelled domain. The additional variable usually does not affect this. So you can change the diffusion constant all day and the mass flow rate will remain unchanged.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

April 22, 2018, 07:47 

#7  
Senior Member
GertJan
Join Date: Oct 2012
Location: Europe
Posts: 1,827
Rep Power: 27 
Quote:
Apparently, you don't have convective flow, since you are only solving the diffusion equation for your scalar, not? Then the massflow does not change indeed since this is determined by convection of your fluid, which is absent. But does the scalar variable at your outlet change? If your setup is correct, I think it should change........ Last edited by GertJan; April 22, 2018 at 09:37. 

April 22, 2018, 10:34 
Convective flow

#8 
New Member
Swapnil Rane
Join Date: Sep 2017
Posts: 23
Rep Power: 8 
I’m working on wine cork material which mostly has diffusive flow of air through it. How can make the diffusive term have more impact on the mass flow at outlets?
Should I add convective flow ? And how can I do it? 

April 22, 2018, 15:56 

#9 
Senior Member
GertJan
Join Date: Oct 2012
Location: Europe
Posts: 1,827
Rep Power: 27 
Diffusion is about Concentration differences. So, as long as you define mass fractions (or Mole fraction), you have to count the molecules that travel through the cork and convert that into an amount of mass.
I don't expect CFX to be the most appropriate software to solve this. It sounds like a quasi 1D/2D diffusion problem. You'd either solve it analytical or numerical using Matlab or similar.......... 

April 22, 2018, 20:49 
Model

#10 
New Member
Swapnil Rane
Join Date: Sep 2017
Posts: 23
Rep Power: 8 
I've the solution through MATLAB.
But I'm trying to model the flow pattern using CFD software. How can I add concentration of oxygen to the body? i.e concentration of oxygen is 1.487 times that of outlets. Can't I use diffusive transport equation for my problem? If Yes, Can anyone set up a simple problem and upload it? 

April 23, 2018, 08:33 

#11 
Senior Member
Join Date: Jun 2009
Posts: 1,800
Rep Power: 32 
If you already solved in MATLAB, what mathematical problem did you solve ?
If you state the mathematical problem, i.e. the equations and the boundary conditions, those in the forum may be able to help. Otherwise, the description is getting lost in translation. 

April 23, 2018, 12:46 
model

#12 
New Member
Swapnil Rane
Join Date: Sep 2017
Posts: 23
Rep Power: 8 
Cork is a material containing oxygen gas in its pores. At the time of production of wine bottles, the cork is compressed and placed inside the bottle.
I'm using the diffusive transport equation to simulate the flow pattern. Boundary conditions are 1) The cork is compressed which leads to a pressure (1.487 atm/1.487 oxygen concentration) in the cork body. 2) The two outlets, one in the bottle and other to atmosphere are at atmospheric pressure (1 atm/zero oxygen concentration) 3) The diffusion coefficient is 2E10 [m^2 s^1] 4) Total time 6 months, timestep is one month To get 1) The total mass flow graph at the outlets to be same. 2) The start point and end point of the total mass flow graph to be same, even after changing the diffusion coefficient by any factor. (just the curve should change) Do you need any more information? 

April 23, 2018, 15:01 

#13 
Senior Member
Join Date: Jun 2009
Posts: 1,800
Rep Power: 32 
Thank you for the physical description of your model, interesting.
However, you have not provided the mathematical model you are using. What equation are you trying to solve? In your MATLAB model, which PDE's are you solving? Are you solving for the mass fraction/concentration of O2 in the cork only? Which PDE equation, which terms are active? If your setup is an elliptic PDE for the mass conservation of O2 in the cork only, i.e. d (density * Y_O2)/dt  div (gamma * grad (Y_O2)) = 0. NOTICE, there is no advection. Your boundary conditions must be setup as either mass fraction for O2, or mass fluxes for O2. The total mass flow of O2 will be a function of those mass fractions. If your setup includes the NavierStokes equations in the porous media, are you including the bottle in the calculation, or just the cork? 

April 23, 2018, 17:00 
Equation

#14 
New Member
Swapnil Rane
Join Date: Sep 2017
Posts: 23
Rep Power: 8 
I'm just trying to solve the diffusive transport equation (ANSYS CFX solver theory guide 11.0 ) Eq 137 page 41 (I've attached a screenshot of the equation)
It's NOT a porous media problem. I just want to solve the equation using diffusion coefficient of oxygen through wine cork. 

April 23, 2018, 18:18 

#15 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,695
Rep Power: 143 
This does not sound like a CFD question, it looks like just a diffusion equation. The Navier Stokes equations are for when fluid inertia is significant and that does not appear to be the case for you.
If this is correct and you want to do this in CFX I would turn off the fluid solvers using expert parameters and just leave the additional variables on. Set up an additional variable to do the diffusion equation and this should be a simple and quick simulation to perform. And by the way  do not say the time step is 1 month unless you have verified that 1 month is acceptable. A very common mistake is to use too large a time step. You should determine the time step size based on a sensitivity study. Finally: Please do not post multiple identical posts.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

April 23, 2018, 21:17 
expert parameter

#16 
New Member
Swapnil Rane
Join Date: Sep 2017
Posts: 23
Rep Power: 8 
Which parameters should I turn off? How can I do it?
In comsol, this problem is mainly solved under the section of transport of concentrated species with the above equation. Can I use any other ansys product for this problem? I know the timestep and total time sounds absurd, but the flow of oxygen through wine cork is important for a year of storage. 

April 23, 2018, 21:41 

#17 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,695
Rep Power: 143 
You can turn the solvers off with expert parameters. Look in the documentation under expert parameters for how to do this. Then you will just be left with the additional variable solver operating. But note you will still have to set up some form of flow condition so CFXPre accepts it  it just won't do anything with it if you turn the solver off.
You missed my point with my comment about the time step. The total simulation time is set by the physics of what you are looking at. If you want to see what happens after a year then the total time will be a year. But the time step is set by what the numerical model needs to run accurately. You need to determine what time step will give an accurate simulation. It is a very common beginner mistake to guess what the user thinks is a good time step and it turns out to be far too large  don't fall into this trap.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

April 25, 2018, 08:00 

#18 
Senior Member
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,167
Rep Power: 23 
Sounds like a 1 dimensional diffusion problem? Center of the cork is a symmetry plane? Since oxygen goes out both directions to 0%?
I solved this years ago, (very early in my career) and made an equation based off the results of numerical models. There already existed an equation for diffusion of semiinfinite thickness, but I made an equation for finite thickness material, which to my knowledge didn't exist before. My analysis uses water concentration and outgassing rate from teflon, but this can really be any scalar: heat, mass, etc. My paper: http://microboonedocdb.fnal.gov/cgi....pdf&version=1 Basic form of Equation: Your "thickness" would be 1/2 of the cork thickness, since it is a symmetry plane at the center (if that is correct) 

April 25, 2018, 08:07 
Model

#19 
New Member
Swapnil Rane
Join Date: Sep 2017
Posts: 23
Rep Power: 8 
Can you upload a archieve file for reference? I really got confused with all discussions and solutions.


April 25, 2018, 08:42 

#20 
Senior Member
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,167
Rep Power: 23 
Here is your solution:
Note the volume porosity of 3% that I'm just guessing on, but should be part of your equation they way I'm understanding it. You don't care about all the stuff with the diffusion coefficient changing with temperature. You just need the basic form of the equation which I posted earlier. 

Tags 
additional variable, diffusive term 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Adding diffusion term to interFoam transport equation  Gearb0x  OpenFOAM Programming & Development  3  February 14, 2023 04:16 
Solving Transport equation for additional variable using CFX  creddy_trddc  CFX  11  April 2, 2015 21:07 
solve only transport equation and not NavierStokes  cosine  CFX  2  November 19, 2011 02:58 
UDF Transport equation for transition momentum thickness Reynolds number  Wantami  FLUENT  0  April 9, 2011 02:50 
error in COMSOL:'ERROR:6164 Duplicate Variable'  bhushas  COMSOL  1  May 30, 2008 04:35 