Example slides from a few Prosper-based seminar presentations · 2001. 11. 13. · Results(displayedasapseudo-animation) Surface-tension-driven collapse of a liquid-lined elastic

Example slides from a few Prosper-basedseminar presentations

...with lots of graphics

Matthias Heil & Andrew L. Hazel.

M.Heil@maths.man.ac.uk & ahazel@maths.man.ac.uk

http://www.maths.man.ac.uk/

�

mheil

Department of Mathematics

University of Manchester

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.1/4

Lagrangian wall mechanics

R

(ζ)

ζ

nt ww

w

w(ζ)

r (ζ)v

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.2/4

Lagrangian wall mechanics

R

(ζ)

ζ

nt ww

w

w(ζ)

r (ζ)v

�

Principle of virtual displacements for a linearly elastic ring of undeformed radius��� and thickness

�

, subject to load

�

:

����

� �

�

� ����

� �� �� �

�

� ����

� � � � ����

�� � � ���

� � ����� �

�� � ���

� ! " #%$

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.2/4

Lagrangian wall mechanics

R

(ζ)

ζ

nt ww

w

w(ζ)

r (ζ)v

�

Principle of virtual displacements for a linearly elastic ring of undeformed radius��� and thickness

�

, subject to load

�

:

����

� �

�

� ����

� �� �� �

�

� ����

� � � � ����

�� � � ���

� � ����� �

�� � ���

� ! " #%$

� and� represent the ring’s mid-plane strain and change of curvature, respectively.

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.2/4

Lagrangian wall mechanics

R

(ζ)

ζ

nt ww

w

w(ζ)

r (ζ)v

�

Principle of virtual displacements for a linearly elastic ring of undeformed radius��� and thickness

�

, subject to load

�

:

����

� �

�

� ����

� �� �� �

�

� ����

� � � � ����

�� � � ���

� � ����� �

�� � ���

� ! " #%$

�

Load non-dimensionalised by bending stiffness

�

, i.e.

� � " � �

where

� "

�

� � � � ��� �

���� �

$

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.2/4

Fluids: Weak form of equations

Momentum: ��

� ���� �

��

��� �� ��� �

��� ��

��� ���� �

��� �

d

� " #

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.3/4

Fluids: Weak form of equations

Momentum: ��

� ���� �

��

��� �� ��� �

��� ��

��� ���� �

��� �

d

� " #[...integrate by parts.]

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.3/4

Fluids: Weak form of equations

Momentum: �

� � ���� �

�� ��� �

��� ��

��� ���� �

� � � ���� �

�

d

�

�

�� � � � �

� ��� ���� �

���� �

��� ��

� ��

� �

d

� " #

[...split the surface integral into d

� " d

��� �

d

��� ��� ]

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.3/4

Fluids: Weak form of equations

Momentum: �

� � ���� �

�� ��� �

��� ��

��� ���� �

� � � ���� �

�

d

�

�

�� � � � �

� ��� ���� �

���� �

��� ��

� ��

� �

d

��� ���

�

�� � � � �

� ��� ���� �

���� �

��� ��

� ��

� �d

�� " #

[Note:

� � " #

on

��� � � ]

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.3/4

Fluids: Weak form of equations

Momentum: �

� � ���� �

�� ��� �

��� ��

��� ���� �

� � � ���� �

�

d

�

�

�� � � � �

� ��� ���� �

���� �

��� ��

� ��

� �

d

�� " #

[...use the traction boundary condition on the free surface� � .]

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.3/4

Fluids: Weak form of equations

Momentum: �

� � ���� �

�� ��� �

��� ��

��� ���� �

� � � ���� �

�

d

�

�

�� �

�

Ca

��

� �� � d

�� " #[...apply Weatherburn’s surface divergence theorem to the surface integral.]

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.3/4

Fluids: Weak form of equations

Momentum: �

� � ���� �

�� ��� �

��� ��

��� ���� �

� � � ���� �

�

d

�

�

Ca

�

��� � � � � � �� � �

� ! � � � � �� � �

� ! �

�

���� � � � � � �

� � �� ! � � � � �

� � �� ! � d

��

� � �

� �� � d

�� �

Ca

� � � � d � " #

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.3/4

Fluids: Weak form of equations

Momentum: �

� � ���� �

�� ��� �

��� ��

��� ���� �

� � � ���� �

�

d

�

�

Ca

�

��� � � � � � �� � �

� ! � � � � �� � �

� ! �

�

���� � � � � � �

� � �� ! � � � � �

� � �� ! � d

��

� � �

� �� � d

�� �

Ca

� � � � d � " #

�

�

�

�

� �

� �

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.3/4

Fluids: Weak form of equations

Momentum: �

� � ���� �

�� ��� �

��� ��

��� ���� �

� � � ���� �

�

d

�

�

Ca

�

��� � � � � � �� � �

� ! � � � � �� � �

� ! �

�

���� � � � � � �

� � �� ! � � � � �

� � �� ! � d

��

� � �

� �� � d

�� �

Ca

� � � � d � " #

Conservation of mass: ��� ���� �

� � �

d� " #

Non-penetration on free surface:

� � � � � � �

d

�� " #

�

Discretise using Taylor–Hood elements

�

Solve matrix equations directly by frontal method (HSL 2000)

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.3/4

Results (displayed as a pseudo-animation)

Surface-tension-driven collapse of a liquid-lined elastic ring:

x1

x 2

0 0.2 0.4 0.6 0.8 1 1.2 1.40

0.2

0.4

0.6

0.8

1P

-1.093

-1.159

-1.226

-1.292

-1.359

= 5.65 x 10-2 t / TL = 17.55* ^

�

Fully coupled discretisation of the free-surface Navier-Stokes equations and theequations of large-displacement shell theory.

�

Solution by the Newton-Raphson method.Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.4/4

Results (displayed as a pseudo-animation)

Surface-tension-driven collapse of a liquid-lined elastic ring:

x1

x 2

0 0.2 0.4 0.6 0.8 1 1.2 1.40

0.2

0.4

0.6

0.8

1P

-0.400

-0.656

-0.912

-1.169

-1.425

= 6.96 x 10-3 t / TL = 19.40* ^

�

Fully coupled discretisation of the free-surface Navier-Stokes equations and theequations of large-displacement shell theory.

�

Solution by the Newton-Raphson method.Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.4/4

Results (displayed as a pseudo-animation)

Surface-tension-driven collapse of a liquid-lined elastic ring:

x1

x 2

0 0.2 0.4 0.6 0.8 1 1.2 1.40

0.2

0.4

0.6

0.8

1P

-0.283

-0.643

-1.003

-1.362

-1.722

=1.78 t / TL = 21.19* ^

�

Fully coupled discretisation of the free-surface Navier-Stokes equations and theequations of large-displacement shell theory.

�

Solution by the Newton-Raphson method.Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.4/4

Results (displayed as a pseudo-animation)

Surface-tension-driven collapse of a liquid-lined elastic ring:

x1

x 2

0 0.2 0.4 0.6 0.8 1 1.2 1.40

0.2

0.4

0.6

0.8

1P

2.000

1.000

0.000

-7.000

-14.000

=39.7 t / TL = 21.19* ^

�

Fully coupled discretisation of the free-surface Navier-Stokes equations and theequations of large-displacement shell theory.

�

Solution by the Newton-Raphson method.Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.4/4

Results (displayed as real animation)

Self-excited oscillations during finite Reynolds number flow in a collapsible channel.

Start external animation

If this doesn’t work, check that (i) the fileStream.avi is located in the same directory asyour pdf file when you display these slides and(ii) that acroread knows how to display avi files(maybe download a later version...)

Back to the previous page?

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.5/4

Results (displayed as real animation)

Self-excited oscillations during finite Reynolds number flow in a collapsible channel.

Start external animation

If this doesn’t work, check that (i) the fileStream.avi is located in the same directory asyour pdf file when you display these slides and(ii) that acroread knows how to display avi files(maybe download a later version...)

Back to the previous page?

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.5/4

Results (displayed as real animation)

Self-excited oscillations during finite Reynolds number flow in a collapsible channel.

Start external animation

If this doesn’t work, check that (i) the fileStream.avi is located in the same directory asyour pdf file when you display these slides and(ii) that acroread knows how to display avi files(maybe download a later version...)

Back to the previous page?

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.5/4

Results (displayed as real animation)

Self-excited oscillations during finite Reynolds number flow in a collapsible channel.

Start external animation

If this doesn’t work, check that (i) the fileStream.avi is located in the same directory asyour pdf file when you display these slides and(ii) that acroread knows how to display avi files(maybe download a later version...)

Back to the previous page?

Matthias Heil & Andrew L. Hazel, Department of Mathematics, University of Manchester, UK http://www.maths.man.ac.uk/

�

mheil – p.5/4