galina lyutskanova kiril mihaylov vasil kolev instructor: tihomir ivanov

24
Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Upload: buddy-cunningham

Post on 05-Jan-2016

228 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Galina LyutskanovaKiril Mihaylov

Vasil Kolev

Instructor: Tihomir Ivanov

Page 2: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Surface tensionSurface tension is a property of fluids that makes them take a shape that minimizes their surface area.

Page 3: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

ExperimentA drop of a given fluid is subjected to rotation or to gravitational forces. Then, a photo of the drop is taken. This is the experimental profile of the drop.

Page 4: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Problem

Software that is used to control the measuring devices finds the value of the surface tension for which the theoretical and experimental profiles coincide.

Page 5: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Drop shape analysis methods

The drop shape analysis are :easy to do;can be used in difficult experimental conditions;require only a small amount of the liquid material;can be used for real-time estimations.

Page 6: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Axisymmetric drop shape analysis (ADSA)

ADSA is a powerful drop shape method. It is fast, easy to handle and produce accurate results. The corner stone of this method is the fact that any given drop is axisymmetric. Using this assumption with the help of the Young-Laplace equation we can efficiently analyze the shape of any given drop.

ADSA is used in various systems such as tensiometers.

Page 7: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Implementation of ADSATwo possible settings are pendant drop and rotating drop. In both the settings we approach the problem in the same way.We create a program that processes the image of the drop that we are going to examine in order to obtain the cloud of points. Then, we acquire differential equations that describe the influence of the interfacial tension and the gravity on the shape of the drop. These equations are dependent on parameters with the help of which we can identify the interfacial tension.

Page 8: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Give approximate discrete solutions of the differential equations using the Euler and RK methods.

By optimizing the error of our approximations we identify the parameters on which our equations depend and find the interfacial tension.

The only difference is in the differential equations but it has no substantial effects on any of the steps.

Page 9: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov
Page 10: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Steps for Contour Extraction

Page 11: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov
Page 12: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

The Young – Laplace Equation of Capillarity

Page 13: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

The model for the pendant drop

Page 14: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov
Page 15: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Substituting (2) and (3) into (1) and parameterizing the curve with the arc length s, we obtain

z

With initial conditions

b

Page 16: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Euler method Runge-Kutta method

Page 17: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

ODE45 – One-step solver, based on a Runge-Kutta methodODE113 – Multistep solver, based on the Adams-Bashforth Methods

big numerical

errors

Page 18: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

b = 1 c = [ 0.05 – 2.5]

Page 19: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov
Page 20: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov
Page 21: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

OptimizationGradient descent Gauss-Newton

Page 22: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Gradient descent Gauss-Newton

Page 23: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Gradient descent Gauss-Newton

Page 24: Galina Lyutskanova Kiril Mihaylov Vasil Kolev Instructor: Tihomir Ivanov

Thank you for your attention!