CFD Online Discussion Forums

CFD Online Discussion Forums (
-   Main CFD Forum (
-   -   Traction boundary conditions (

m malik January 11, 2004 11:04

Traction boundary conditions
I have the following question on 3-dimensional FEM solution of momentum-continuity equations. Consider a simple capillary flow problem with only the pressure being specified at the two ends; velocities not specified. How does one apply the traction conditions: ,i,e, tx = (mu)(du/dz+dw/dx), ty = (mu)(dv/dz+dw/dy), and tz = -p+2(mu)(dw/dz) when velocity field is not know a-priori? (z = coordinate along the capillary axis). Thanks.

xueying January 12, 2004 19:40

Re: Traction boundary conditions
The traction equals t = -p n + mu n.(du/dx + transpose(du/dx)) where t, n, u, x are vectors

In my code, I specify the traction as follows, If the flow is fully developed, then n.(du/dx)=0, so t = -p n + mu n.(transpose(du/dx)) The traction is encoded thru the weak form integral, du/dx is computed from the finite element basis functions and coefficients.

You can also refer to papers related with open flow boundaries.

m malik January 13, 2004 09:37

Re: Traction boundary conditions
I am dealing with the polymer flows, and in my code I apply the traction boundary condition in the following manner. The nonlinear finite element equations are solved using the Picard iteration scheme. In each iteration, the traction vector is calculated using the external pressure and the mu.(dui/dxj + duj/dxi) terms from the velocity field of the previous iteration. No assumption of fully developed flow is imposed. Is this the right approach?

xueying January 13, 2004 12:56

Re: Traction boundary conditions
If it's not fully developed flow and it's a viscous flow, your imposing method should be correct, at least we impose exactly the same traction.

Since your flow is polymer flow, do you have elastic stress contribution from polymer? In our code, if it's polymer flow, then it's viscoelastic flow, the traction has contributions from pressure, viscous stress, and elastic stress.

m malik January 13, 2004 13:49

Re: Traction boundary conditions
Thanks Xueying for your response. No, I do not have elastic stress contribution. You are probably talking about the elastic dumbbell models. I am just considering steady state viscoplatic flows. However, irrespective of the model, why should one impose fully developed flow condition? For example in tube flow, the solution must approach to fully developed with increasing aspect (length-to-diameter) ratio. (This is one of the things that is bothering me.)

Rami January 14, 2004 05:15

Re: Traction boundary conditions

<head> <meta http-equiv=Content-Type content="text/html; charset=windows-1255"> <meta name=Generator content="Microsoft Word 10 (filtered)"> <title>Re: Traction boundary conditions</title>




<body lang=EN-US>

<div class=Section1>

<p class=MsoNormal style='text-align:justify;text-indent:.2in'>I think a formal and general answer would be the simplest way to answer your question. Tensor notation is used below, in which a repeated index implies summation. In addition, "," stands for covariant derivative (which, for Cartesian coordinates, degenerates to partial derivative).

<p class=MsoNormal style='text-align:justify;text-indent:.2in'>The traction vector, t<sub>k</sub>, in a given direction n<sub>i</sub> is defined as

<p class=MsoNormal style='margin-left:.5in;text-align:justify;text-indent:.2in'>t<sub>k</sub> = <span style='font-family:Symbol'>s</span><sub>ik</sub> n<sub>i</sub>

<p class=MsoNormal style='text-align:justify'>where <span style='font-family: Symbol'>s</span><sub>ik</sub> is the stress tensor.

<p class=MsoNormal style='text-align:justify;text-indent:.2in'>The constitutive relation for a fluid with no elasticity is

<p class=MsoNormal style='margin-left:.5in;text-align:justify;text-indent:.2in'><span style='font-family:Symbol'>s</span><sub>ik</sub> = -p <span style='font-family:Symbol'>d</span><sub>ik</sub> + <span style='font-family: Symbol'>m</span><sub>v</sub> <span style='font-family:Symbol'>d</span><sub>ik</sub> <span style='font-family:Symbol'>e</span><sub>mm</sub> + 2<span style='font-family:Symbol'>m</span> (<span style='font-family:Symbol'>e</span><sub>ik</sub> 1/3 <span style='font-family:Symbol'>d</span><sub>ik</sub> <span style='font-family:Symbol'>e</span><sub>mm</sub>)

<p class=MsoNormal style='text-align:justify'>where p is the pressure, <span style='font-family:Symbol'>m</span> and <span style='font-family:Symbol'>m</span><sub>v</sub> are the viscosity and the bulk viscosity, respectively, <span style='font-family:Symbol'>d</span><sub>ik</sub> is Kronecker's delta and the strain-rate tensor, <span style='font-family:Symbol'>e</span><sub>ik</sub> , is related to the velocity vector, u<sub>k</sub>, by

<p class=MsoNormal style='margin-left:.5in;text-align:justify;text-indent:.2in'><span style='font-family:Symbol'>e</span><sub>ik</sub> = 1/2 (u<sub>i,k</sub> + u<sub>k</sub><sub>,i</sub>)

<p class=MsoNormal style='text-align:justify;text-indent:.2in'>In the finite element method, the shape-functions relate the nodal values to the element distributions. Therefore, having the nodal values (at any iteration level), the traction may be calculated by the above equations.

<p class=MsoNormal style='text-align:justify;text-indent:.2in'>A developed flow assumption is not required in a pipe flow. You may definitely solve for a developing flow. For a long enough domain you will notice the flow approaches the fully developed profile asymptotically.

<p class=MsoNormal align=center style='text-align:center;text-indent:.2in'>I hope this helps,

<p class=MsoNormal align=center style='text-align:center;text-indent:.2in'>Rami




hesamking March 29, 2012 12:46

How can I apply free traction BC in CFX?
my BC is:
thank you so much

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