22.11.2014 Views

Download full text - ELSA - Europa

Download full text - ELSA - Europa

Download full text - ELSA - Europa

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Permanent FSI treatment (2)<br />

• Upon discretization, contact pressure is replaced by interaction force r<br />

• Interaction force<br />

is the resultant of<br />

the contact pressure<br />

at each node of the<br />

interface<br />

• For the moment, assume nodal conformity at the F-S interface<br />

• Velocity compatibility condition vFin = vSin<br />

is of the form Cv = b<br />

therefore one may use Lagrange multipliers to find r (see Part 1)<br />

r<br />

• However, the following apparently simple question arises:<br />

How does one define<br />

“the” normal to a<br />

discrete F-S interface?<br />

11<br />

The 2-D plane case<br />

Pioneering work (1980s) by Donea, Giuliani, Halleux:<br />

n 1<br />

n<br />

n 2<br />

L 2<br />

• Physical intuition: the<br />

discrete normal is some<br />

average of the adjacent<br />

sides normals: n = ( n1+ n2)/<br />

n1+<br />

n2<br />

y<br />

x<br />

L 1<br />

• Works well only in 2D<br />

plane cases, and for uniform<br />

mesh: L1 = L2<br />

• In more general cases, the mass balance is incorrect: some fluid is<br />

“gained” or “lost” at the interface corners<br />

12<br />

6

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!