# Time step vs Negative concentration

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

 April 17, 2006, 10:38 Time step vs Negative concentration #1 James Guest   Posts: n/a Hi,all For problems including both fast and slow reactions, taking a big timestep may results in negative concentrations. This always happens in my study. Some people said that using a dimensionless concentration(mass fraction) can solve the above problem. Is this correct? Can someone help clarify this and give theoritical explaination? Thanks! Best regards! James

 April 17, 2006, 18:47 Re: Time step vs Negative concentration #2 Glenn Horrocks Guest   Posts: n/a Hi, It is a boundedness issue. This is spurious numerical results, often caused by high order differencing schemes, but highly non-linear processes like chemistry can cause it too. How to fix it depends on the exact problem. An alternative approach might work, also simply putting bounds on the variable values (ie fixing it to greater than or equal to zero) can also work. Glenn Horrocks

 April 18, 2006, 15:28 Re: Time step vs Negative concentration #3 James Guest   Posts: n/a Thanks, Glenn I am actually using step fuctions in cfx to limit the concentration above zero, simialr as you suggested. However I am interested in if using dimensionless concentrations could be another option. Best regards! James

 April 18, 2006, 17:45 Re: Time step vs Negative concentration #4 Glenn Horrocks Guest   Posts: n/a Hi, How are you using the step functions? Is the concentration an algebraic function or is a convected variable? If it is a convected variable (which is probably what it is) then the step functions will not limit it. You will have to use the variable bonudednesses limiters which I think are available under the solver tab or expert parameters. Glenn Horrocks

 April 18, 2006, 18:29 Re: Time step vs Negative concentration #5 James Guest   Posts: n/a Thanks again, Glenn I define the concentrations as volumetric additional variables(scalars)(I guess it's convected variable in your term?). I am also using 'step function' to terminate the source term when the concentration becomes less than zero. What is the 'variable bonudednesses limiters'? Do you mean 'scalar diffusion scheme' or 'bounded bnd diffusion tets', or others? I found following information in help document, but haven't yet understand how to use them. 'scalar diffusion scheme' Default Value 2 (tetrahedral) 1 (hexagonal) Description Specify whether the standard central scheme default) or the positive definite scheme is applied to the scalar equations. This parameter may be of use in obtaining convergence with poor quality meshes. A value of 2 sets the scalar diffusion scheme to be positive definite. bounded bnd diffusion tets Description: Controls whether or not a bounded diffusion scheme is used at boundaries for tetrahedral elements. The diffusion scheme for other element types is bounded by default through the default setting of 5 for the expert parameters "cht diffusion scheme" and "scalar diffusion scheme". Tetrahedral elements are controlled separately because the bounded scheme may be lead to less accurate results. Parameter Type: String Allowed Values: t,f Default Value: f Do you have more suggestions? Best regards! James

 April 19, 2006, 18:40 Re: Time step vs Negative concentration #6 Glenn Horrocks Guest   Posts: n/a Hi, Yes, I mean volumetric additional variables. I don't think just limiting the source term will stop negative values. It is not a true boundedness limit. I have had a look around and cannot find a way of setting the boundedness of a variable. The 'scalar diffusion scheme' or 'bounded bnd diffusion tets' you mention are different things and not what I mean. Have a look in the book "Computational Fluid Dynamics" by Roache as he has a good discussion of this issue. Glenn Horrocks

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post CKH OpenFOAM Running, Solving & CFD 12 March 21, 2016 14:05 pUl| FLUENT 31 August 21, 2015 04:46 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03 nileshjrane OpenFOAM Running, Solving & CFD 8 August 26, 2010 12:50 panda60 OpenFOAM Bugs 4 July 8, 2010 10:23

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