# how to solve for additional term?

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

 January 22, 2002, 11:56 how to solve for additional term? #1 Bram Guest   Posts: n/a Dear all, I have a set of transport equations for scalars with some extra terms are involved apart from convection and diffusion. One of these extra terms is something as a second order derivative for the scalars in a virual space (not (xyz)space). Is it possible to add a step such as above by solving for the n-th timestep this second order derivative and then add it to the source term (with usrsrc), so that it will be accounted for in the (n+1)time step? Has anyone of you experience with this? Looking forward to lots of replies, Bram PS I use CFX 4.3

 January 27, 2002, 17:30 Re: how to solve for additional term? #2 Astrid Guest   Posts: n/a 1) How does your virtual space depend on xyz? Not at all? Then, what is it? 2) I think it is possible. I have to check whether the second derivative is directly available. 3) I have no experience with virtual spaces, but we can always try to get you on the road. 4) Why don't you use CFX4.4?? Veel succes, Astrid

 January 28, 2002, 07:39 Re: how to solve for additional term? #3 Bram Guest   Posts: n/a Thanks a lot for your answer! It might be a bit abstract but I'll try to explain the problem as good as possible. It's a combustion problem and the 'virtual space' is mixture fraction space [0..1]. That space itself does not depend on xyz. The second order derivative of the scalar can be considered as a mixing/diffusion parameter of the scalar in that space. The mixture fraction space is discretised into a number of points. This gives intervals, and for each interval a scalar transport equation is being solved throughout the flow domain. However this transport equation contains this second order derivative in mixture fraction space as well. And that's where the bottleneck is located at the moment.... I found something in the manual which allows the user to modify the coefficients of the differencing scheme and it could be that that opens a way to solve this extra term. (The actual values of the scalar for each interval are known and the discretisation of mixture fraction space is known as well, so somehow it should be possible to calculate coeffients which could be added to the matrix..) So indeed it might be possible. What do you think? Unfortunately CFX 4.4 is not available to me. Would it make a huge difference? Bram PS Beetje veel tekst en waarschijnlijk niet echt lekker duidelijk uitgelegd, maar ik ben benieuwd!

 January 28, 2002, 17:36 Re: how to solve for additional term? #4 Astrid Guest   Posts: n/a 1) You don't need CFX 4.4 2) It sound like you have implemented the mixture fraction space succesfully. Correct? How did you implement it? As a UserScalar? 3) Indeed you can modify your discretisation scheme using USRFDS. Perhaps you need the second derivative of your mixture fraction. Therefore you might use the utility routine GRADS(.....). Then you have to use it twice to get a second order derivatie. Not sure if it will work. Therefore I need to know more of your (typisch universiteits) project. Astrid

 January 29, 2002, 05:56 Re: how to solve for additional term? #5 Bram Guest   Posts: n/a Indeed, mixture fraction space is implemented succesfully (vingers crossed). I stored it as an array in permanent workspace. The second derivate I just calculate with a finite difference calculation. The outcome of that I add to the solution by putting it into the USRSRC for each specific scalar. I skipped the idea of modifying the discretisation scheme, as the matrix being used there is one in physical space. Het is inderdaad een typisch universiteits project. Het doel is om gedeeltelijk voorgemengde methaanvlammen te kunnen modelleren. (reactants prediction, NOx prediction, temperatures etc.). Thanks for your suggestions again! Bram PS Ik heb respect voor je 'posting' tijden. Of zit je aan de andere kant van onze planeet?

 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 siw CFX 1 December 6, 2010 07:41 titio OpenFOAM Running, Solving & CFD 0 May 4, 2009 09:17 SAKTI CD-adapco 1 August 28, 2008 01:38 segersson OpenFOAM Running, Solving & CFD 5 March 3, 2006 00:06 Brian FLUENT 1 October 24, 2005 09:15

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