CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Why half angle in strain rate tensor derivation?

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 15, 2017, 13:18
Default Why half angle in strain rate tensor derivation?
  #1
Member
 
Jonny
Join Date: Aug 2009
Posts: 72
Rep Power: 17
Jonny6001 is on a distinguished road
I'm working through derivation of the fundamental equations and need to clarify one point regarding elements of the strain rate tensor.

Why do we use half the combination of (dU/dy+dV/dx) etc for the off diagonal terms, rather than d/dy(mu.dU/dy)? I'm finding it difficult to justify this since it's essentially saying that a change in the 'V' velocity is imparting a force on the fluid inside a control volume in the 'X' direction.

I've looked through many text books and websites but have not come across any resource that goes through the derivation entirely. I see how the maths works but not the overall concept of why we need to combine and then split the terms in half.

I look forward to any replies.
Jonny6001 is offline   Reply With Quote

Old   October 15, 2017, 13:56
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Jonny6001 View Post
I'm working through derivation of the fundamental equations and need to clarify one point regarding elements of the strain rate tensor.

Why do we use half the combination of (dU/dy+dV/dx) etc for the off diagonal terms, rather than d/dy(mu.dU/dy)? I'm finding it difficult to justify this since it's essentially saying that a change in the 'V' velocity is imparting a force on the fluid inside a control volume in the 'X' direction.

I've looked through many text books and websites but have not come across any resource that goes through the derivation entirely. I see how the maths works but not the overall concept of why we need to combine and then split the terms in half.

I look forward to any replies.

I am not sure about your question, are you asking why the Newtonian model defines the stress tensor in terms of the symmetric gradient velocity?
FMDenaro is offline   Reply With Quote

Old   October 15, 2017, 15:47
Default
  #3
Member
 
Jonny
Join Date: Aug 2009
Posts: 72
Rep Power: 17
Jonny6001 is on a distinguished road
Hi FMDenaro,

If we look at the 'X' momentum equation, and the shear stress part:

d(tau_xx)/dx + d(tau_yx)/dy + d(tau_zx)/dz

where d(tau_yx)/dy represent shear stress acting on face normal to 'Y' axis in 'X' direction.

Using the Newton's constitutive law: tau_yx = mu.du/dy

We have d(tau_yx)/dy = d/dy(mu.du/dy) = mu d2y/dy2 (assuming constant viscosity)

This is the formulation that's obtained for incompressible flow. But for compressible flow we end up with:

(tau_yx) = mu.(du/dy + dv/dx)

And (tau_xy) = mu.(du/dy + dv/dx)

Why have we included the dv/dx in the to tau_yx term?

I appreciate that a fluid element deforms in both direction (in 2D) and that the angle that the parallelogram forms with the original vertical side is given by:

(dU/dy).dy.dt

But in my mind, only the vertical sides of the fluid element will be rotated when a shear stress is applied on the top face. I get that the horizontal sides are also rotated, further reducing the angle of the lower left corner of the fluid element. But this requires a shear stress in the 'Y' direction, why does it appear in the 'X' momentum equation?

Apologies if my question is still not clear.
Jonny6001 is offline   Reply With Quote

Old   October 15, 2017, 15:58
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Because in the incompressible flow the divergence-free constrain applies on the velocity field and some terms elide. Therefore you get (mu constant)

Div (2mu S) = mu Lap v
FMDenaro is offline   Reply With Quote

Old   October 15, 2017, 16:30
Default
  #5
Member
 
Jonny
Join Date: Aug 2009
Posts: 72
Rep Power: 17
Jonny6001 is on a distinguished road
Yes, I appreciate all that. The issue is not with the incompressible version.

Perhaps I'm not explaining my train of though very well. In summary, why does dv/dx appear in the 'X' momentum equation?

I'm looking for a more detailed answer than "because the shear stress is given by (tau_yx) = mu.(du/dy + dv/dx)"
Which is what all the textbooks that I've picked up say.

My previous reply has more detail about the problem that I have with the derivation.

I think the answer is along the lines that the angular deformation of the element is averaged over the X and Y directions (in 2D) to maintain zero net moment on the fluid element. But that's not exactly correct.

Thank you!
Jonny6001 is offline   Reply With Quote

Old   October 15, 2017, 16:51
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
The presence of the two derivatives is in the nature of the symmetry of S so that the key of the question is why the Newton law link the stress tensor to the simmetric part of the gradient. To answer, one needs to recall the conservation of the angular momentum for a fluids that leads to the symmetry condition of the stress tensor. That should answer why we have the mutual action of the variation of u and v, along y and x, respectively.
FMDenaro is offline   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
[openSmoke] libOpenSMOKE Tobi OpenFOAM Community Contributions 562 January 25, 2023 10:21
How to define strain rate and shear rate magnitude in CFD-Post? khhan0543 Fluent UDF and Scheme Programming 3 January 13, 2017 15:38
How to determine the mean filtered strain rate for the Smagorinsky model levinperson Main CFD Forum 5 November 30, 2015 14:00
flow rate for half geometry alecon88 FLUENT 0 February 16, 2012 07:17
[Netgen] Import netgen mesh to OpenFOAM hsieh OpenFOAM Meshing & Mesh Conversion 32 September 13, 2011 06:50


All times are GMT -4. The time now is 17:46.