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

Second derivative with cfx 5.7.1

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 20, 2006, 04:46
Default Second derivative with cfx 5.7.1
  #1
Alessandro
Guest
 
Posts: n/a
Hi,i want to create a source term which depends oin a second derivative of velocity. How can I get second derivative within CFX? I can get all componenets of first derivative of velocity but I don't know how to get second derivatives of these components from first derivatives? Should I compile a second fortran subroutine with User Get var? Can anyone help me?
  Reply With Quote

Old   April 20, 2006, 10:51
Default Re: Second derivative with cfx 5.7.1
  #2
Alessandro
Guest
 
Posts: n/a
I have to calculate the second derivative on velocity within CFX 5.7.1. Can anyone help me? Thank you. Alessandro
  Reply With Quote

Old   April 20, 2006, 11:43
Default Re: Second derivative with cfx 5.7.1
  #3
opaque
Guest
 
Posts: n/a
Dear Alessandro,

Would you mind sharing the origin of your source term? If it comes from a div (coef * grad Phi) form. It is better to assemble it as a summation of fluxes around a control volume than as a source term.. In that case you do not need second derivatives at all..

Good luck, Opaque

  Reply With Quote

Old   April 20, 2006, 12:21
Default Re: Second derivative with cfx 5.7.1
  #4
Alessandro
Guest
 
Posts: n/a
Hi Opaque, thank you for your answer. I'm modeling an entrained bed gasifier and I want to insert in the momentum equation a term which is neglected in CFX model. This term is the solid stress tensor
  Reply With Quote

Old   April 20, 2006, 12:31
Default Re: Second derivative with cfx 5.7.1
  #5
opaque
Guest
 
Posts: n/a
Dear Alessando,

Have you checked ANSYS CFX 10.0? I understand that the solid pressure term (you called solid stress tensor) is already included..

In any case, the solid stress tensor has the form I mentioned:

div( coef1 * (grad U + grad^T U) + coef2 * div U)

so, there is no need for second derivatives..

Good luck, Opaque

  Reply With Quote

Old   April 20, 2006, 12:39
Default Re: Second derivative with cfx 5.7.1
  #6
Alessandro
Guest
 
Posts: n/a
How can I take this term into account without second derivative evaluation? Thank you
  Reply With Quote

Old   April 20, 2006, 12:45
Default Re: Second derivative with cfx 5.7.1
  #7
opaque
Guest
 
Posts: n/a
Dear Alessandro,

ANSYS CFX uses finite volume approach; therefore, you are taking the volume integral of the term

div( coef1 * (grad U + grad^T U) + coef2 * div U)

Using Gauss theorem

volInt (div Phi) dV = areaInt (Phi dot dA)

Then, you must sum the fluxes around a control volume.. Those fluxes are

sum around CV [coef1 * (grad U + grad^T U) + coef2 * div U) dot dA ]

Of course, you can only do this via User Fortran.. It is not trivial, but that is the "correct way" to do it..

You should contact the CFX help desk for more information on how to do such work..However, I rather upgrade (besides many other benefits that 10.0 has)..

Good luck, Opaque

  Reply With Quote

Old   April 20, 2006, 13:00
Default Re: Second derivative with cfx 5.7.1
  #8
Alessandro
Guest
 
Posts: n/a
Thank you opaque. I can get solid velocity gradient through user fortran subroutines. If I refer to a summation of fluxes around a control volume then , as you already said, I don't need second derivatives at all. But a question comes, how can I insert the term in the momentum equation for the solid dispersed phase?
  Reply With Quote

Old   April 21, 2006, 10:41
Default Re: Second derivative with cfx 5.7.1
  #9
Alessandro
Guest
 
Posts: n/a
Opaque, have you any suggestion how to add the terms to the momentum equation?
  Reply With Quote

Old   April 21, 2006, 10:52
Default Re: Second derivative with cfx 5.7.1
  #10
opaque
Guest
 
Posts: n/a
Dear Alessandro,

My best advice here is to contact the CFX help desk.. This is not a trivial operation.. Here is why:

- You must create a subdomain

- Add a momemtum source term, and use User Fortran functions to specify the values

Now the User Fortran

- The source term will be called for a particular group of elements. Be careful with the ENTITY (please read theory doc + USER_CALC_INFO doc).

- Within the group of elements you must gather the Normal Area vector on integration points, velocity gradient on integration points, molecular coefficients at i.p's and multiply them as needed. For example, coef*gradU dot A_ip

- Add the fluxes to the corresponding entities as required.. For this one, you must understand how ANSYS CFX discretizes the equations.

I am still convinced that you should upgrade to ANSYS CFX 10.0.

Good luck, Opaque
  Reply With Quote

Old   April 21, 2006, 11:31
Default Re: Second derivative with cfx 5.7.1
  #11
Alessandro
Guest
 
Posts: n/a
Thank you opaque. I also thik i should upgrade to cfx 10. Unfortunately we have cfx 5.7.1 and i think we cannot upgrade.. About what you saild we have already create a subdomain, a user function and routine related to a fortran subroutine which calculates the velocity gradient.. the subdomain requires a source term to be inserted specifying the different components, x,y,z components. So, we have some chances..don't you thik so?
  Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to creat a mesh to use in CFX 5.7.1 ! Trung CFX 5 September 8, 2010 05:32
CFX 5.7.1 PRE and solver won't start daniel CFX 1 January 20, 2006 11:09
amd athlon64 cpu and cfx 5.7.1 solver cmete CFX 3 October 20, 2005 23:16
USRBCS in CFX 5 (actually, 5.7.1) ATS CFX 2 August 24, 2005 13:30
Available memory is less than 2GB in CFX 5.7.1 Korsh Mik CFX 3 June 5, 2005 20:34


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