CN112179813A - Liquid contact angle on-line measurement method based on experimental image - Google Patents

Liquid contact angle on-line measurement method based on experimental image Download PDF

Info

Publication number
CN112179813A
CN112179813A CN202010870388.6A CN202010870388A CN112179813A CN 112179813 A CN112179813 A CN 112179813A CN 202010870388 A CN202010870388 A CN 202010870388A CN 112179813 A CN112179813 A CN 112179813A
Authority
CN
China
Prior art keywords
contact angle
liquid
predicted
differential
measurement result
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010870388.6A
Other languages
Chinese (zh)
Other versions
CN112179813B (en
Inventor
李坤
吕存景
郝鹏飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN202010870388.6A priority Critical patent/CN112179813B/en
Publication of CN112179813A publication Critical patent/CN112179813A/en
Application granted granted Critical
Publication of CN112179813B publication Critical patent/CN112179813B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • G01N13/02Investigating surface tension of liquids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • G01N13/02Investigating surface tension of liquids
    • G01N2013/0208Investigating surface tension of liquids by measuring contact angle

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)
  • Measuring Arrangements Characterized By The Use Of Fluids (AREA)

Abstract

An on-line measuring method of liquid contact angle based on experimental image, the method is based on static or quasi-static real-time experimental image directly obtained in experiment, extract the interface outline with Canny method; then, solving a series of predicted contours by changing the predicted contact angle serving as a boundary condition by using a Young-Laplace equation; and comparing the predicted profile with the experimental profile to calculate a norm error so as to determine a contact angle and an error limit. The method converts the measurement of the contact angle into the error contrast of the profile, not only has higher measurement precision, but also can reflect the change of the contact angle in the experiment caused by the instability of the coating. The measuring method can be used for measuring the contact angle on a plane and spherical substrate, does not need a special measuring device, can be applied to common horizontal drop experiments, and can also be applied to the contact angle measurement in liquid bridge, capillary lifting and liquid experiments containing air holes, which cannot be realized by the conventional method.

Description

Liquid contact angle on-line measurement method based on experimental image
Technical Field
The invention relates to an on-line measuring method of a liquid contact angle, in particular to an on-line measuring method of a liquid contact angle based on an experimental image, and belongs to the technical field of physical measurement.
Background
The contact angle of solid-liquid-gas or solid-liquid is an important index for measuring the affinity of three phases of solid-liquid-gas/solid-liquid, is closely related to the wetting state of the three phases of solid-liquid-gas/solid-liquid, and is an important characterization index commonly used in the fields of industrial production and scientific research. The current method for measuring the contact angle is mainly the horizontal dropping method, namely, the contact angle is measured by using the contact profile of a liquid drop and a substrate by dropping a small liquid drop on the surface of a solid flat sample. The method needs a special measuring instrument and requires the solid sample to be in a horizontal position, and does not meet special experimental requirements. The method is used for measuring the contact angle of the surface of the material before the experiment begins, and for some surfaces with unstable coatings, the method can not accurately reflect the actual condition of the contact angle at the later stage of the experiment and can not monitor the change condition of the contact angle in real time, so that the effect of the coating in the experiment can not be fed back in time.
Disclosure of Invention
The invention aims to provide an on-line liquid contact angle measuring method based on an experimental image, aiming at the problems that a special measuring instrument is needed in the current liquid contact angle measurement, the contact angle change caused by coating instability in the experimental process cannot be reflected in real time, and a flat substrate needs to be specially processed.
The technical scheme of the invention is as follows:
the measuring process is realized by predicting the contact angle, further predicting the outline and comparing the predicted outline with the experimental outline to calculate the norm error, and the norm error between the predicted outline and the experimental outline is taken as the standard for measuring the measuring effect; finally obtaining a minimum predicted contact angle with norm error meeting the target precision, a maximum predicted contact angle and a predicted contact angle with the minimum norm error, namely an upper limit of a contact angle measurement result, a lower limit of the contact angle measurement result and an optimal value of the contact angle measurement result; the method comprises the following specific steps:
1) selecting a static or quasi-static solid-liquid-gas or solid-liquid three-phase contact area image from a real-time experiment video to obtain an experiment image;
2) selecting the position of the solid substrate: if the substrate is a flat substrate, selecting two points on the substrate, and if the substrate is a spherical substrate, selecting three points on the substrate;
3) after the position of the solid substrate is determined, calculating the gray change gradient of different positions by using a canny algorithm according to the gray distribution of an experimental image, and identifying the interface profile by using a double-threshold method;
4) performing polynomial fitting on the interface profile extracted in the step 3), and performing smoothing treatment on the interface profile;
5) by observing an experimental image, giving a value which is obviously smaller than an actual contact angle as an iteration initial value of a predicted contact angle;
6) using the predicted contact angle as a boundary condition of the predicted interface contour, and solving coordinate increment of any two adjacent interface points by using a difference format of a Young-Laplace equation so as to obtain all points on the predicted contour;
7) comparing the predicted contour obtained in the step 6) with the experimental contour processed in the step 4), calculating a norm error, if the norm error is greater than the target precision, increasing a predicted step length for the predicted contact angle, and repeating the step 6); if the norm error is smaller than the target precision, taking the predicted contact angle as the lower limit of the contact angle measurement result, and returning the lower limit of the contact angle measurement result;
8) predicting the contact angle to increase by a prediction step; using the predicted contact angle as a boundary condition of the predicted interface contour, and solving coordinate increment of any two adjacent interface points by using a difference format of a Young-Laplace equation so as to obtain all points on the predicted contour;
9) comparing the predicted contour obtained in the step 8) with the experimental contour processed in the step 4), calculating a norm error, and repeating the step 8) if the norm error is smaller than the target precision; if the norm error is larger than the target precision, subtracting a step length from the predicted contact angle to serve as the upper limit of the contact angle measurement result, stopping stepping, and returning to the upper limit of the contact angle measurement result;
10) and taking the predicted contact angle with the minimum norm error between the lower limit of the contact angle measurement result and the upper limit of the contact angle measurement result as the optimal value of the contact angle measurement result, returning the optimal value of the contact angle measurement result, and finishing the measurement.
In the above technical solution, the differential format of the Young-Laplace equation includes the following two types:
a) two-dimensional planar conditions:
Figure BDA0002650905060000021
wherein γ is the surface tension between the liquid-gas or liquid-liquid interfaces;
Figure BDA0002650905060000022
is the tangential angle differential; ds is the differential of the interface arc length; x is the abscissa; z is the ordinate; dx is the abscissa differential; dz is the ordinate differential; Δ p is the pressure difference at the reference position, which needs to satisfy the condition: Δ p ═ ρ gz0Wherein z is0Is the liquid phase depth at the reference position; rho is the density difference of liquid-gas or liquid-liquid two phases; g is the acceleration of gravity;
b) the case of axial symmetry:
Figure BDA0002650905060000031
wherein γ is the surface tension between the liquid-gas or liquid-liquid interfaces;
Figure BDA0002650905060000032
is the tangential angle differential; ds is the differential of the interface arc length; r is the abscissa; z is the ordinate; dr is the abscissa differential; dz is the ordinate differential; Δ p is the pressure difference at the reference position, which needs to satisfy the condition: Δ p ═ ρ gz0Wherein z is0Is the liquid phase depth at the reference position; rho is the density difference of liquid-gas or liquid-liquid two phases; g is the acceleration of gravity.
Compared with the prior art, the invention has the following advantages and prominent technical effects: firstly, a real-time experiment image directly obtained from an experiment is utilized, and a special measuring device is not needed; for the unstable condition of the material surface coating, the change condition of the contact angle in the experimental process can be reflected; the application field is expanded from a conventional horizontal drop experiment to a capillary lifting experiment, a liquid bridge experiment and a liquid experiment containing air holes, and the method is also suitable for measuring the contact angle on the spherical substrate with curvature; and thirdly, the predicted contour is solved by a self-adaptive method, the measurement of the contact angle is converted into the contrast of the contour, and the higher precision is achieved.
Drawings
FIG. 1 is a flow chart of the measurement method of the present invention.
FIG. 2 is a front view of the backlight of the collapse experiment of the air holes in the liquid film.
FIG. 3 is a comparison of the predicted profile and the experimental profile for the example predicted contact angle of 152 degrees.
FIG. 4 is a comparison of the predicted profile and the experimental profile for the example where the predicted contact angle is 154 degrees.
FIG. 5 is a graph comparing the predicted profile and the experimental profile when the predicted contact angle is 156 degrees in the examples.
Fig. 6 is a corresponding relationship between the predicted contact angle and the error norm in the embodiment.
Detailed Description
The present invention will be described in detail with reference to the accompanying drawings and examples.
According to the liquid contact angle on-line measuring method based on the experimental image, a special measuring device is not needed, the experimental image is directly used for realizing on-line contact angle measurement, and the method is not limited to a horizontal drop experimental image; the change of the contact angle in the experiment can be reflected in the case of unstable coating, and the method can be used for measuring the contact angle of not only a plane substrate but also a spherical substrate. The measuring process is realized by predicting the contact angle to further obtain a predicted contour and comparing the predicted contour with an experimental contour to calculate a norm error, and the norm error between the predicted contour and the experimental contour is used as a standard for measuring the measuring effect; the method finally obtains the minimum predicted contact angle, the maximum predicted contact angle and the predicted contact angle with the minimum norm error, wherein the norm error meets the target precision, namely the optimal values of the upper limit of the contact angle measurement result, the lower limit of the contact angle measurement result and the contact angle measurement result. The method specifically comprises the following steps:
1) selecting a static or quasi-static experimental image of a solid-liquid-gas or solid-liquid three-phase contact area needing to measure a contact angle from a real-time experimental video;
2) sequentially selecting points on the solid substrate, returning the enter key to the positions of the points on the substrate after the selection is finished, and treating the substrate as a flat substrate if two points are selected; if three points are selected, treating the substrate as a spherical substrate;
3) referring to the solid substrate position determined in the step 2), calculating the gray level change gradient of different positions by using a canny algorithm according to the gray level distribution of the experimental image, wherein two adjustable threshold values, namely a high threshold value and a low threshold value, are provided as parameters for controlling the boundary identification continuity: the pixel points with the gray gradient larger than the high threshold are considered as strong boundary points; pixel points with gray gradients smaller than the low threshold are not considered as boundary points; the gray gradient is between two thresholds, and the pixel points near the strong boundary point are considered as weak boundary points; the gray gradient is between two thresholds, but the pixel point which is not near the strong boundary point is not considered as the boundary point;
4) performing polynomial fitting on the interface contour extracted in the step 3) to perform smoothing processing on the boundary, wherein the polynomial fitting is used for extracting boundary points and comparing the boundary points with a predicted contour obtained by subsequent calculation, the order of the polynomial fitting can be adjusted, and the fitting condition is observed in real time;
5) according to the experimental image, giving a value which is obviously smaller than the actual contact angle as an initial predicted contact angle;
6) and (3) taking the predicted contact angle as a boundary condition of the predicted interface profile, and solving the corresponding predicted interface profile by utilizing a differential format of a Young-Laplace equation: the program can calculate the coordinate increment of any two adjacent boundary points according to the difference format, and the coordinate of each point can be calculated by increment accumulation only by knowing the predicted contact angle of the coordinate of the first point;
i. the differential format of the Young-Laplace equation comprises two conditions:
a) two-dimensional case:
Figure BDA0002650905060000041
wherein γ is the surface tension between the liquid-gas or liquid-liquid interfaces;
Figure BDA0002650905060000042
is the tangential angle differential; ds is the differential of the interface arc length; x is the abscissa; z is the ordinate; Δ p is the pressure difference at the reference position, which needs to satisfy the condition: Δ p ═ ρ gz0Wherein z is0Is the liquid phase depth at the reference position; rho is the density difference of liquid-gas or liquid-liquid two phases; g is the acceleration of gravity; dx is the abscissa differential; dz is the ordinate differential;
b) the axisymmetric case is:
Figure BDA0002650905060000051
wherein γ is the surface tension between the liquid-gas or liquid-liquid interfaces;
Figure BDA0002650905060000052
is the tangential angle differential; ds is the differential of the interface arc length; r is the abscissa; z is the ordinate; Δ p is the pressure difference at the reference position, which needs to satisfy the condition: Δ p ═ ρ gz0Wherein z is0Is the liquid phase depth at the reference position; rho is the density difference of liquid-gas or liquid-liquid two phases; g is the acceleration of gravity; dr is the abscissa differential; dz is the ordinate differential;
using the predicted contact angle as a boundary condition, and solving the predicted profile comprises the following steps:
a) the predicted contact angle is taken as the coordinate s is 0, and x is x1Z is 0 or s is 0 and r is r1The tangent angle coordinate of the point where z is 0, i.e. the tangent angle coordinate of the first point, where x1And r1The abscissa value of the outline contact part identified in the step 4;
b) solving the coordinate increment from the first point to the second point by utilizing a Young-Laplace equation difference format, and summing the coordinate increment with the coordinate of the first point to obtain the coordinate of the second point;
c) solving the coordinate increment from the second point to the third point by utilizing a Young-Laplace equation difference format, and summing the coordinate increment with the coordinate of the second point to obtain the coordinate of the third point; repeating the steps until the tangent angle is larger than 180 degrees or smaller than minus 180 degrees, and obtaining a series of points as points on the predicted contour;
when the predicted contour is solved in the step 6, the parameters of the Young-Laplace equation, namely the reference point pressure difference delta p corresponding to the predicted contact angle, need to be targeted, and the step is as follows:
a) giving an initial estimated upper limit and a lower limit of the pressure difference of the reference point;
b) respectively taking the estimated upper limit and the estimated lower limit as Young-Laplace equation parameters, calculating two groups of points on the predicted contour, wherein the tangent angle coordinates of the two groups of contour termination points meet a positive-negative condition;
c) gradually reducing the difference between the estimated upper limit and the estimated lower limit of the pressure difference by utilizing a dichotomy method, and stopping iteration until the difference between the estimated upper limit and the estimated lower limit is smaller than the target precision;
d) taking the average value of the pressure difference estimation upper limit and the estimation lower limit as the reference point pressure difference corresponding to the predicted contact angle, and taking the pressure difference as a Young-Laplace equation parameter to obtain a predicted profile, namely the predicted profile corresponding to the predicted contact angle;
7) comparing the predicted contour obtained in the step 6 with the experimental contour processed in the step 4, calculating a norm error, if the norm error is greater than the target precision, increasing a predicted step length for the predicted contact angle, and repeating the step 6; if the norm error is smaller than the target precision, taking the value as the lower limit of the contact angle measurement result, and returning the lower limit of the contact angle measurement result;
i. step 11, calculating the norm of the norm error includes the following steps:
a) 1-norm:
Figure BDA0002650905060000061
b) 2-norm:
Figure BDA0002650905060000062
c) p-norm:
Figure BDA0002650905060000063
wherein E represents a norm, x1,x2,...,xnTo be the radial coordinate of the point on the experimental profile,
Figure BDA0002650905060000064
predicting the radial coordinate of a point on the contour, wherein p is any positive integer;
8) predicting the contact angle to increase by a prediction step; solving the predicted interface contour by taking the predicted contact angle as a boundary condition of the predicted interface contour and utilizing a difference format of a Young-Laplace equation, wherein the process is similar to the step 6;
9) comparing the predicted contour obtained in the step 8 with the experimental contour processed in the step 4, calculating a norm error, and repeating the step 8 if the norm error is smaller than the target precision; if the norm error is larger than the target precision, subtracting one step length from the predicted contact angle to serve as the upper limit of the contact angle measurement result, stopping stepping, and returning to the upper limit of the contact angle measurement result;
10) and taking the predicted contact angle with the minimum norm error between the lower limit of the contact angle measurement result and the upper limit of the contact angle measurement result as the optimal value of the contact angle measurement result, returning the optimal value of the contact angle measurement result, and finishing the measurement.
Example (b):
taking the collapse experiment of the air holes in the liquid film as an example, the measurement process is as follows:
1) selecting an experimental image in the quasi-static shrinkage process of the air hole boundary before the liquid film is collapsed as an experimental image of a contact angle to be detected, as shown in fig. 2, wherein a black area is an air hole, and a white area is a liquid film;
2) selecting two points on a substrate, and determining a solid flat substrate;
3) setting the high and low threshold values of the identification image as 120 and 40 respectively;
4) the fitting order of the interface profile is set to 9
5) The contact angle in the experimental image is obviously larger than 150 degrees, 150 degrees are selected as initial prediction contact angles, and the prediction step length is 0.1 degree;
6) selecting an axisymmetric format by a Young-Laplace equation;
7) setting the estimated upper limit of the pressure difference of the reference point to be 60 Pa, the estimated lower limit to be 40 Pa and the target precision of the pressure difference to be 0.001;
8) setting a contact angle norm error to adopt a 2-norm, wherein the target precision is 10;
9) the program sequentially displays a contact angle measurement result lower limit profile comparison graph 3, a contact angle measurement result optimal value profile comparison graph 4, a contact angle measurement result upper limit profile comparison graph 5 and a norm error and predicted contact angle corresponding relation graph 6, wherein the upper limit of the contact angle measurement result, the optimal value of the contact angle measurement result and the position where the iteration of the predicted contact angle is stopped (the predicted contact angle is larger than the upper limit of the contact angle measurement result by one predicted step length) are marked in the graph 6;
10) the program returns the measurement results: the lower limit of the contact angle measurement result is 151.9 degrees, the optimal value of the contact angle measurement result is 154.4 degrees, and the upper limit of the contact angle measurement result is 156.7 degrees.

Claims (2)

1. An on-line measuring method for liquid contact angle based on experimental image is characterized in that the method comprises the following steps:
1) selecting a static or quasi-static solid-liquid-gas or solid-liquid three-phase contact area image from a real-time experiment video to obtain an experiment image;
2) selecting the position of the solid substrate: if the substrate is a flat substrate, selecting two points on the substrate, and if the substrate is a spherical substrate, selecting three points on the substrate;
3) after the position of the solid substrate is determined, calculating the gray change gradient of different positions by using a canny algorithm according to the gray distribution of an experimental image, and identifying the interface profile by using a double-threshold method;
4) performing polynomial fitting on the interface profile extracted in the step 3), and performing smoothing treatment on the interface profile;
5) by observing an experimental image, giving a value which is obviously smaller than an actual contact angle as an iteration initial value of a predicted contact angle;
6) using the predicted contact angle as a boundary condition of the predicted interface contour, and solving coordinate increment of any two adjacent interface points by using a difference format of a Young-Laplace equation so as to obtain all points on the predicted contour;
7) comparing the predicted contour obtained in the step 6) with the experimental contour processed in the step 4), calculating a norm error, if the norm error is greater than the target precision, increasing a predicted step length for the predicted contact angle, and repeating the step 6); if the norm error is smaller than the target precision, taking the predicted contact angle as the lower limit of the contact angle measurement result, and returning the lower limit of the contact angle measurement result;
8) predicting the contact angle to increase by a prediction step; using the predicted contact angle as a boundary condition of the predicted interface contour, and solving coordinate increment of any two adjacent interface points by using a difference format of a Young-Laplace equation so as to obtain all points on the predicted contour;
9) comparing the predicted contour obtained in the step 8) with the experimental contour processed in the step 4), calculating a norm error, and repeating the step 8) if the norm error is smaller than the target precision; if the norm error is larger than the target precision, subtracting a step length from the predicted contact angle to serve as the upper limit of the contact angle measurement result, stopping stepping, and returning to the upper limit of the contact angle measurement result;
10) and taking the predicted contact angle with the minimum norm error between the lower limit of the contact angle measurement result and the upper limit of the contact angle measurement result as the optimal value of the contact angle measurement result.
2. The method according to claim 1, wherein the differential format of the Young-Laplace equation comprises the following two cases:
a) two-dimensional planar conditions:
Figure FDA0002650905050000021
wherein γ is the surface tension between the liquid-gas or liquid-liquid interfaces;
Figure FDA0002650905050000022
is the tangential angle differential; ds is the differential of the interface arc length; x is the abscissa; z is the ordinate; dx is the abscissa differential; dz is the ordinate differential; Δ p is the pressure difference at the reference position, which needs to satisfy the condition: Δ p ═ ρ gz0Wherein z is0Is the liquid phase depth at the reference position; rho is the density difference of liquid-gas or liquid-liquid two phases; g is the acceleration of gravity;
b) the case of axial symmetry:
Figure FDA0002650905050000023
wherein γ is the surface tension between the liquid-gas or liquid-liquid interfaces;
Figure FDA0002650905050000024
is the tangential angle differential; ds is the differential of the interface arc length; r is the abscissa; z is the ordinate; dr is the abscissa differential; dz is the ordinate differential; Δ p is the pressure difference at the reference position, which needs to satisfy the condition: Δ p ═ ρ gz0Wherein z is0Is the liquid phase depth at the reference position; rho is the density difference of liquid-gas or liquid-liquid two phases; g is the acceleration of gravity.
CN202010870388.6A 2020-08-26 2020-08-26 Liquid contact angle on-line measurement method based on experimental image Active CN112179813B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010870388.6A CN112179813B (en) 2020-08-26 2020-08-26 Liquid contact angle on-line measurement method based on experimental image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010870388.6A CN112179813B (en) 2020-08-26 2020-08-26 Liquid contact angle on-line measurement method based on experimental image

Publications (2)

Publication Number Publication Date
CN112179813A true CN112179813A (en) 2021-01-05
CN112179813B CN112179813B (en) 2021-08-06

Family

ID=73923939

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010870388.6A Active CN112179813B (en) 2020-08-26 2020-08-26 Liquid contact angle on-line measurement method based on experimental image

Country Status (1)

Country Link
CN (1) CN112179813B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113029870A (en) * 2021-03-31 2021-06-25 华中科技大学 Three-phase contact angle measuring method of partial phase solution liquid system and application thereof

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02287139A (en) * 1989-04-28 1990-11-27 Natl Res Inst For Metals Surface tension measuring instrument
TW200728703A (en) * 2006-01-17 2007-08-01 Picvue Optoelectronics International Inc Method for the surface tension measurement of melting glasses at high temperature
CN102539289A (en) * 2010-09-30 2012-07-04 三菱综合材料株式会社 Liquid surface free energy and solid-liquid contact angle, interface free energy and prediction method for attaching energy
CN102809524A (en) * 2012-08-20 2012-12-05 华北电力大学(保定) Static contact angle calculation method based on modified ellipse fitting method
CN103247047A (en) * 2013-04-23 2013-08-14 重庆交通大学 Image edge detection method based on fractional order partial differential
CN104266940A (en) * 2014-10-21 2015-01-07 华北电力大学(保定) Method for calculating static contact angle of hydrophobic material
CN104568672A (en) * 2014-12-04 2015-04-29 上海梭伦信息科技有限公司 Device and method for testing interfacial tension and contact angle by adopting spinning drop method under ultrahigh pressure and at high temperature
CN104697900A (en) * 2013-12-05 2015-06-10 上海梭伦信息科技有限公司 Method for measuring interface tension and contact angle under vacuum controlled-temperature conditions, and apparatus thereof
CN106092833A (en) * 2016-08-28 2016-11-09 宁波新边界科学仪器有限公司 Use sessile drop method all automatic measurement liquid surface or the method for interfacial tension and device thereof
WO2018039038A1 (en) * 2016-08-26 2018-03-01 Board Of Regents, The University Of Texas System Measuring contact angles beteween a solid-fluid pair using x-ray imaging of the solid-fluid-fluid interface inside a capillary
CN107817194A (en) * 2017-10-27 2018-03-20 清华大学 Weak hydrophobic surface Static Contact angle measuring method and system
US20190347375A1 (en) * 2018-05-08 2019-11-14 Fei Company Wettability estimation by differential multi-phase simulation

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02287139A (en) * 1989-04-28 1990-11-27 Natl Res Inst For Metals Surface tension measuring instrument
TW200728703A (en) * 2006-01-17 2007-08-01 Picvue Optoelectronics International Inc Method for the surface tension measurement of melting glasses at high temperature
CN102539289A (en) * 2010-09-30 2012-07-04 三菱综合材料株式会社 Liquid surface free energy and solid-liquid contact angle, interface free energy and prediction method for attaching energy
CN102809524A (en) * 2012-08-20 2012-12-05 华北电力大学(保定) Static contact angle calculation method based on modified ellipse fitting method
CN103247047A (en) * 2013-04-23 2013-08-14 重庆交通大学 Image edge detection method based on fractional order partial differential
CN104697900A (en) * 2013-12-05 2015-06-10 上海梭伦信息科技有限公司 Method for measuring interface tension and contact angle under vacuum controlled-temperature conditions, and apparatus thereof
CN104266940A (en) * 2014-10-21 2015-01-07 华北电力大学(保定) Method for calculating static contact angle of hydrophobic material
CN104568672A (en) * 2014-12-04 2015-04-29 上海梭伦信息科技有限公司 Device and method for testing interfacial tension and contact angle by adopting spinning drop method under ultrahigh pressure and at high temperature
WO2018039038A1 (en) * 2016-08-26 2018-03-01 Board Of Regents, The University Of Texas System Measuring contact angles beteween a solid-fluid pair using x-ray imaging of the solid-fluid-fluid interface inside a capillary
CN106092833A (en) * 2016-08-28 2016-11-09 宁波新边界科学仪器有限公司 Use sessile drop method all automatic measurement liquid surface or the method for interfacial tension and device thereof
CN107817194A (en) * 2017-10-27 2018-03-20 清华大学 Weak hydrophobic surface Static Contact angle measuring method and system
US20190347375A1 (en) * 2018-05-08 2019-11-14 Fei Company Wettability estimation by differential multi-phase simulation

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
QIANG YU 等: "Droplet Image Feedback Control System in Evaportion Experiment", 《QRIGINAL ARTICLE》 *
冯艳辉: "SJ-10蒸发对流实验液滴图像处理软件的研究与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
徐志钮: "快速ADSA-P算法及用于超疏水材料接触角计算", 《仪器仪表学报》 *
王应明 等: "基于不同误差准则和范数的组合预测方法研究", 《控制与决策》 *
管宏蕊 等: "图像边缘检测经典算法研究综述", 《首都师范大学(自然科学版)》 *
高蔷 等: "利用Young-Laplace 算法程序计算液体表面张力", 《科学技术与工程》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113029870A (en) * 2021-03-31 2021-06-25 华中科技大学 Three-phase contact angle measuring method of partial phase solution liquid system and application thereof

Also Published As

Publication number Publication date
CN112179813B (en) 2021-08-06

Similar Documents

Publication Publication Date Title
CN117173189B (en) Visual inspection system for polishing effect of aluminum alloy surface
Cabezas et al. A new method of image processing in the analysis of axisymmetric drop shapes
CN105404861A (en) Training and detecting methods and systems for key human facial feature point detection model
CN108447070B (en) Industrial part defect detection algorithm based on pixel vector invariant relation features
CN112179813B (en) Liquid contact angle on-line measurement method based on experimental image
CN108629312A (en) A kind of dynamic hand tracking method and system
CN116758077B (en) Online detection method and system for surface flatness of surfboard
CN106289678B (en) A kind of bent detection device of cylindrical body abnormity drum and recognition methods
CN115096206A (en) Part size high-precision measurement method based on machine vision
CN107610174B (en) Robust depth information-based plane detection method and system
CN115239661A (en) Mechanical part burr detection method and system based on image processing
CN116883412B (en) Graphene far infrared electric heating equipment fault detection method
CN116664584B (en) Intelligent feedback regulating system for production of thin-wall zinc alloy die casting die
CN113095323A (en) SIFT-improvement-based digital image correlation method real-time detection method
CN110514110B (en) Platform leveling control method
KR20030034274A (en) Generalized photometric stereo for the hybrid reflectance model using an image sequences
CN109727258B (en) Image segmentation method based on regional gray heterogeneous energy
CN113327232B (en) Relay contact set contact surface alignment method, apparatus and computer storage medium
Barone et al. Grid pattern for in-plane strain measurements by digital image processing
CN115761017B (en) Visual determination method for alkyd resin paint color
CN112734768B (en) Image vanishing point searching method and device and electronic equipment
CN116681697B (en) Cobalt removal measuring method and device for diamond compact and electronic equipment
CN117058130B (en) Visual inspection method for coating quality of optical fiber drawing surface
CN117929226A (en) Pore structure analysis method, device, equipment and medium of scanning electron microscope image
CN117788423A (en) Fitting contour degree determining method, device, equipment and storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant