Next Article in Journal
Dual-Parameter Prediction of Downhole Supercritical CO2 with Associated Gas Using Levenberg–Marquardt (LM) Neural Network
Next Article in Special Issue
Computational Fluid Dynamics Study on Bottom-Hole Multiphase Flow Fields Formed by Polycrystalline Diamond Compact Drill Bits in Foam Drilling
Previous Article in Journal
Numerical Study of Laminar Flow and Vortex-Induced Vibration on Cylinder Subjects to Free and Forced Oscillation at Low Reynolds Numbers
Previous Article in Special Issue
Slow Translation of a Composite Sphere in an Eccentric Spherical Cavity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability or Instability of a Static Liquid Bridge Appearing in Shaped Crystal Growth from Melt via the Pulling-Down Method

by
Andreea V. Cojocaru
1,† and
Stefan Balint
2,*,†
1
Independent Researcher, 300223 Timisoara, Romania
2
Department of Computer Science, West University of Timisoara, 300223 Timisoara, Romania
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fluids 2024, 9(8), 176; https://doi.org/10.3390/fluids9080176
Submission received: 19 June 2024 / Revised: 26 July 2024 / Accepted: 29 July 2024 / Published: 31 July 2024
(This article belongs to the Special Issue Multiphase Flow and Granular Mechanics)

Abstract

:
This study presents sufficient conditions for the stability or instability of the static liquid bridge appearing in crystal growth from the melt of micro-fibers, thin plates, and hollow micro-tubes of predetermined sizes using the pulling-down method. The case in which the contact angle and the growth angle verify the inequality α c > π / 2 α g is considered. Experimentally, only stable static liquid bridges can be created; unstable static liquid bridges exist just in theory, because in reality they collapse. The results of this study are significant for shaped crystal growth from melted materials, with given macroscopic dimensions, and using specific equipment. This is because the obtained inequalities represent limits for what can and cannot be achieved experimentally.

1. Introduction

Small liquid amounts connecting two solid surfaces, separated by a short distance, are defined as liquid bridges (see Figure 1). These solid surfaces can be surrounded by a gas phase or can be immersed in another liquid immiscible with that of the liquid bridge.
Theoretical and experimental study of these liquid formations has been very extensive during the last century. The shape of liquid bridges of pure liquids is determined by intermolecular forces (liquid–liquid molecular interactions determine surface tension whereas the three-phase contact line is additionally influenced by liquid–solid molecular forces), by gravity, and by external fields such as pressure or electrostatic fields. Although liquid bridges refer to small liquid volumes, they can be found in a range of size scales. Regarding the upper size scale, there is a limit at the physical dimension of about 1 cm. At this size, gravity causes spontaneous collapse of liquid bridges. Even in the absence of gravity (space station experiments or by density matching with the surrounding fluid), there is the so-called Rayleigh stability limit at which some oscillation modes of the bridge become unstable, leading to its disintegration into droplets. For bridges smaller than about 1 mm, the effect of gravity is negligible, and intermolecular forces completely dominate their shape. There is no lower limit to the bridge’s physical size beyond that imposed by the molecular nature of the liquid. The situation in the presence of surfactants is much more complex. These molecules not only affect the equilibrium surface tension of the liquid but also bring about a dynamic interfacial behavior of the bridge as they diffuse slowly from the bulk of the bridge towards the gas–liquid interface; see [1].
There are many technological applications where liquid bridges are encountered; see [1,2,3,4,5,6,7,8]. A non-exhaustive list includes processes such as granulation, flotation, and coating. In mechanical engineering, liquid bridges allow control of forces in micro-gripping processes. In geology, liquid bridges determine the properties of moist soil. They are also related to tertiary oil recovery from porous media. A particular recent area of interest refers to self-assembly of particles for which surface tension is the main driving force. New technologies for fabricating micro-electromechanical systems based on liquid bridges have been proposed. Liquid bridges have been employed as dispensing devices in medical applications. Liquid bridges formed between the tip of the atomic force microscopy probe and the analyzed sample (in humid environments) significantly affect the measurements. Despite the plethora of experimental and theoretical studies on liquid bridges, it is surprising that very few refer to their stability and instability, which factor in shaped single crystal growth from melt material via the edge-defined-film-fed-growth (EFG) and micro-pulling-down (μ-PD) methods. The liquid bridge is a small amount of melt connecting the shaper and the crystal being grown, and in crystal growth this is called the meniscus.
The basic growth methods available for crystal growth are, broadly, growth from the melt, growth from vapor, and growth from solution. The Springer Handbook of Crystal Growth [9] presents in 1816 pages the state of the art in crystal growth until the year 2010. In Chapter 40 pg. 1379–1402 of this book, the authors Th. George, St. Balint, and L. Braescu present several mathematical models describing processes that take place in the case of crystal growth from the melt by Bridgman-Stockbarger (BC) and by the edge-defined-film-fed-growth (EFG) method. For the BC method see Figure 2, and for the EFG method see Figure 3, Figure 4, Figure 5 and Figure 6.
First, a short description is provided of the real processes, which are modeled, along with the equations, boundary conditions, and initial values defining the model. After that, numerical results are provided, obtained from computations in the framework of the model, making a comparison between the computed results and those obtained from other models and with the experimental data. Predictions are made concerning the impurity repartition in the crystal in the framework of each model. The philosophy behind building the mathematical model is to provide a quantitative description of the mechanism by balancing what is wanted with what can be achieved. The quantitative description is made in terms of a certain number of variables (called the model variables) such that the mathematical model is a set of equations concerning these variables. Analysis of the mathematical model leads to results or predictions that can be tested against observations. It is important to realize that all models are idealizations and, hence, are limited in their applicability; experimental verifications can establish the authenticity of the model. In [9] Chapter 40, six models for BS growth (Figure 2) and four models for EFG growth are presented regarding the impurity distribution in the crystal. This is because the quality of the crystal being grown depends to a considerable extent on the homogeneity and distribution of both detrimental and specially added impurities. Variations in the concentrations of impurities along the crystal length and cross section lead to variations in the mechanical, electrical, and optical properties in the mass of the crystal. The impurity repartition in the crystal is determined by the processes that take place in the melt during the growth and mass transfer. Quantitative descriptions of these processes permit the prediction of such repartitioning in the crystal. Consequently, computer modeling of crystal growth processes has often been described as an art as well as a science, to acknowledge the seemingly endless difficulties that arise in the application of numerical methods.
Three core competencies are required for the effective use of modeling in crystal growth. One is the thorough grasp of the fundamentals of continuum transport phenomena. The second is the building of an appropriate mathematical model. The third is a general understanding of the numerical methods necessary for solving the governing equations of transport phenomena. For readers interested in more details regarding the comparison of different crystal growth methods, their advantages and disadvantages, we recommend book [9]. A significant advantage of the EFG method, in comparison with the BS method, is that the crystal is grown without interaction with the crucible, which considerably improves the structural quality of the material: there are fewer residual stresses, dislocations, spurious nucleation, or twins. In the case of the BS method, the melt is encapsulated in a crucible, and the crystallization of the melt takes place in conditions of permanent contact between the melt and the crystal with the inner wall of the crucible (Figure 2). In the case of the EFG method, there is a liquid bridge between the crystal and the shaper (die) called the meniscus. The melt is in a crucible, from which it flows through a capillary tube onto the surface of the shaper (Figure 3, Figure 4, Figure 5 and Figure 6). Here, a liquid bridge is formed between the shaper surface and the crystal. The crystallization takes place on the so-called “crystallization front”, which is the border line between the upper part of the liquid bridge and the bottom of the crystal.
Therefore, the crystal is grown without interaction with the crucible, which considerably improves the structural quality of the material. For more detail, see [10]. “Dewetting” refers to a phenomenon that has occurred spontaneously during many experiments of Bridgman solidification of semiconductors in space (see reviews [11,12,13]). It also refers to a process developed for crystal growth on Earth (see review in [14]). In both cases, the crystal is grown without interaction with the crucible, which considerably improves the structural quality of the material: there are fewer residual stresses, dislocations, spurious nucleation, or twins. The origin of the gap between the crystal and the crucible comes from a small liquid meniscus at the level of the solid–liquid interface [4]. While this phenomenon is spontaneous under microgravity conditions, because of the lack of hydrostatic pressure, it has been adapted on Earth by applying on the liquid a gas pressure difference, of the order of the hydrostatic pressure, in order to create and maintain the meniscus. For the geometry of the growth system and the main dimensions, angles, temperatures, and pressures of interest in the process, see Figure 7.
Many experiments under microgravity have shown that the gap, which is typically smaller than 100 μm, is remarkably constant for several hours of growth. Similarly, dewetted crystals obtained on Earth demonstrate that, under given conditions (essentially a bad wetting of the liquid on the crucible), the crystal radius stays spontaneously constant, while it is almost impossible to obtain a dewetted crystal for other configurations [15,16].
It appears that the process is extremely stable (the grown crystal does not reattach to the crucible wall) in certain cases, while it shows high instability under other conditions. The thin gap thickness is directly linked to the meniscus shape and position, which depend on capillary forces, hydrostatic and hydrodynamic pressures, and heat transfer, all of which likely fluctuate with time. Therefore, in order to master the growth process, it is necessary to perform dynamic stability analysis. Dynamic stability analyses were reported in [17,18].
The micro-pulling-down (μ-PD) technique is a kind of “inverse edge-defined film-fed growth” technique. This means that the melt is pulled downward from a vertical shaper, which is on the bottom part of the crucible. This technique was developed by Fukuda’s laboratory in Japan [19,20,21,22,23,24]. The advantage in the case of the μ-PD technique is that the thermal convection in the liquid bridge (meniscus melt) is zero. The thermal convection can be significant in the case of the edge-defined-film-fed growth technique, even in the case of a static liquid bridge (pulling rate equal to zero). Inequalities presented here are sufficient conditions for the stability or instability of the static liquid bridge appearing in the case of micro-fibers, thin plates, and hollow micro-tube growth using the μ-PD technique. The case in which the contact angle and the growth angle verify the inequality α c > π / 2 α g is considered. The novelty of this article consists of the obtained inequalities that represent limits for what can and cannot be experimentally achieved with a predetermined material, equipment, and single crystal size. These inequalities can be obtained using second-order conditions for the minimum value of the free energy function. To the best of our knowledge, these types of results have never been published in the literature for materials for which the contact angle and the growth angle verify the inequality α c > π / 2 α g .
In Section 2, the stability or instability of the liquid bridge in the case of micro-fiber growth is analyzed. Section 3 is dedicated to a similar analysis in the case of thin plate growth. Section 4 presents the same kind of analysis in the case of hollow micro-tube growth. In Section 5 we present results, and in Section 6 we offer short comments and conclusions.
We emphasize that this study only investigates the stability and instability of the liquid bridge; our research excludes the study of the growth of shaped single crystals. For this reason, the reference list is short in comparison with the existing literature, and the thermal aspect is completely disregarded.

2. Liquid Bridge in Micro-Fiber Growth

The prototype of an axis-symmetric fiber growth system using the μ-PD technique is represented in Figure 8.
In the hydrostatic approximation, the equation describing the equilibrium capillary surface (shape of the liquid bridge surface) is the Laplace capillary equation (see [25] p. 132). However, to calculate the meniscus shape for the capillary-shaping techniques, it is convenient to employ the Laplace equation in its explicit differential form. To obtain this form of the equation following [25] (p. 132), the variation problem of the total free energy minimum of a liquid column will be solved. For our system, the total free energy is composed of the surface free energy and the gravity field energy. The minimum should be calculated with the additional condition of liquid volume constancy. Solving this problem presents some difficulties. Therefore, our study will be restricted to considering a meniscus (liquid bridge) possessing axial symmetry (see Figure 8 and Formula (1)). Such menisci are obtained during melt pulling of straight circular cylinder-shaped crystals (see [10] p. 24). It should be mentioned that the Laplace capillary equation follows from the necessary condition of the minimum of Function (1). Not all the liquid bridge that represents the solution of the Laplace capillary equation will be stable, i.e., will be physically realized. The stability of this liquid bridge (the static one) should be distinguished from the dynamic stability of the crystallization process. For a statically stable liquid bridge, not only necessary (first order) but also sufficient (second order) conditions of Function (1) should be established. The first-order conditions are the Euler equations, and the second-order conditions are the Jacobi differential equation with the Legendre condition.
The energy function of the melt column behind the free surface in terms of the meridian curve profile z ( r ) of the static axis-symmetric liquid bridge is as follows:
I ( z ) = R c R d γ 1 + z 2 1 / 2 1 2 ρ g z 2 p z r d r
where γ is the melt surface tension, ρ is the melt density, g is the gravity, and p = P a + ρ · g · h m is the pressure difference across the free surface (see [26]).
The stability condition of the axis-symmetric static liquid bridge is that z ( r ) minimizes (1).
The first-order condition of the minimum value is the Euler equation
d d r ( F z ) F z = 0
for the following function:
F ( z , z , r ) = γ · 1 + z 2 1 2 1 2 · ρ · g · z 2 p · z · r
In computing the Euler equation for the meridian curve, the following equation is found:
z = ρ g z + p γ 1 + z 2 3 2 1 r 1 + z 2 z R c r R d
The function z ( r ) describes the meridian curve profile and verifies the following conditions:
z ( R c ) = tan π 2 α g ;   z ( R c ) = h c > 0
z ( R d ) = tan α c ;   z ( R d ) = 0
z ( r )   is   strictly   decreasing   on   R c , R d
The first equality in (3) expresses that at the three-phase point R c , h c , where the thermal condition for solidification is realized ( h c is the meniscus height), the angle between the tangent line to the meridian curve of the free surface and the vertical is equal to the growth angle α g .
The second equality in (3) expresses that the meniscus meridian curve of the meniscus is fixed to the outer edge of the crystal at the three-phase point R c , h c .
The first equality in (4) expresses that at the outer edge of the shaper, point ( R d , 0 ) , the angle between the tangent line to the meridian curve of the free surface and the horizontal is equal to the contact angle α c , which is a material “constant” depending on the shaper and the melt.
The second equality in (4) expresses that the meridian curve is fixed at the outer edge of the shaper, i.e., at the point ( R d , 0 ) , where the meridian curve touches the outer edge of the shaper.
Condition (5) concerns the shape of the meridian curve. This formally means z ( r ) < 0 , or concavity. The inequality α c > π / 2 α g is a consequence of (5).
The analysis of stability of the static liquid bridge in the case α c > π / 2 α g is justified by the fact that there exist shapes and crystals for which this inequality works and the stability analysis does not exist.
In the case of microfiber growth by the micro-pulling -down method, our understanding of the concavity condition of the meridian curve z ( r ) is z ( r ) < 0 , i.e., the second derivative is strictly negative. According to this condition, the first-order derivative z r is strictly decreasing. Taking into account Formulas (3) and (4), it follows that α c > π / 2 α g . It can be shown that inequality α c > π / 2 α g implies the first-order derivative of z ( r ) is strictly decreasing from where z ( r ) < 0 in a neighborhood of the point Rd. We point out two things: (a) In Figure 3, a concave meridian curve is represented; in Figure 4, Figure 5, Figure 6, convex meridian curves are represented that correspond to the inequality αc ≤ π/2 − αg. (b) In the case of Figure 1, this kind of understanding of the concepts “concave” and convex is not applicable; that is, because the curve represented in Figure 1 is not the graphic of the function z ( r ) . In the case of Figure 1, to a value of r corresponds two values for z; one is situated in the upper half part of the curve, which is concave, and the other value is situated in the lower half part of the curve, which is convex.
In case of convex liquid bridges, αcπ/2 − αg the second-order derivative of the function z ( r ) is strictly positive, and the first-order derivative of z ( r ) is strictly increasing. The stability of the static liquid bridge in the αcπ/2 − αg case was analyzed in [26,27]. The stability of the liquid bridge in the α c > π / 2 α g case is an open problem.
Using (2)–(5) and the concavity condition for z ( r ) ( z ( r ) < 0 ), the following inequalities can be determined:
γ · α c + α g π 2 R d · n n 1 · cos α c ρ · g R d · n 1 n · tan α c + γ R d · cos α g p γ · α c + α g π 2 R d · n n 1 · sin α g + γ R d · n · sin α c
where n = R d R c > 1 .
Inequalities (6) are necessary conditions for the existence of z ( r ) . This has to be verified by pressure difference p = P a + ρ · g · h m . For stability or instability of the liquid bridge, the second-order conditions of the minimum of (1), i.e., the Legendre condition 2 F z z 0 and the Jacobi equation
d d r [ 2 F z z · z ] [ 2 F z z d d r ( 2 F z z ) ] · η = 0
have to be analyzed.
The Legendre condition in this case is the inequality:
r · γ [ 1 + ( z ) 2 ] 3 2 0
which is obviously verified. The Jacobi equation in this case is the following:
( r · γ [ 1 + ( z ) 2 ] 3 2 · η ) + ρ · g · r · η = 0
To obtain a stability result for the liquid bridge, it is sufficient to identify a “Sturm-type upper bound” ([27] Chapter 11) for the Jacobi Equation (8). For this purpose, we note that the following inequalities concerning the coefficients of Equation (8) hold:
r · γ [ 1 + ( z ) 2 ] 3 2 R d · γ n · cos 3 α c , ρ · g · r ρ · g · R d
It follows that the equation as follows:
[ R d · γ n · cos 3 α c ] · ξ + [ ρ · g · R d ] · ξ = 0
Is a “Sturm-type upper bound“ for Equation (8). An arbitrary solution of Equation (10) is of the form ξ ( r ) = A · sin ( ω t + ϕ ) , where A and ϕ are arbitrary constants and ω 2 = ρ · g · n γ · cos 3 α c . The half period of any non-zero solution of (10) is π ω = π · γ 1 / 2 · cos 3 / 2 α c ρ 1 / 2 · n 1 / 2 · g 1 / 2 . If the half period is more than R d R c , then any non-zero solution of (10) vanishes at most once on the interval [ R c , R d ] .
In other words, if the following inequality holds n 1 n 1 / 2 < π · 1 R d · γ 1 / 2 · cos 3 / 2 α c ρ 1 / 2 · g 1 / 2 , then any non-zero solution ξ ( r ) of (10) vanishes at most once on the interval [ R c , R d ] . Hence, according to [27] Chapter 11, the solution η ( r ) of the Jacobi equation that verifies η ( R d ) = 0 and η ( R d ) = 1 has only one zero on the interval [ R c , R d ] . This means that the Jacobi condition for the weak minimum value is verified (see [27]).
Theorem 1.
If for  n = R d R c > 1  the following inequality holds:
n 1 n 1 / 2 < π · 1 R d · γ 1 / 2 · cos 3 / 2 α c ρ 1 / 2 · g 1 / 2
 then the liquid bridge is stable.
To obtain instability results for the Jacobi equation, a “Sturm-type lower bound equation” ([27] Chapter 11) has to be found. For this purpose, we note that for the coefficients of the Jacobi equation, the following inequalities hold:
r · γ [ 1 + ( z ) 2 ] 3 2 R d γ · sin 3 α g ,   ρ · g · r ρ · g · R d / n
It follows that the following equation:
[ R d · γ · sin 3 α g ] · ς + [ ρ · g · R d / n ] · ζ = 0
is a “Sturm-type lower bound equation” for the Jacobi Equation (8). An arbitrary solution of Equation (13) is of the form ζ ( r ) = A · sin ( ω t + ϕ ) , where A and ϕ are arbitrary constants and ω 2 = ρ · g n · γ · sin 3 α g . The period of any non-zero solution of (13) is 2 · π ω = 2 · π · γ 1 / 2 · n 1 / 2 · sin 3 / 2 α g ρ 1 / 2 · g 1 / 2 . If the period is less than R d R c , then any non-zero solution of (13) vanishes at least twice on the interval [ R c , R d ] . In other words, if the inequality n 1 n 3 / 2 > 2 · π · 1 R d · γ 1 / 2 · sin 3 / 2 α g ρ 1 / 2 · g 1 / 2 holds, then any non-zero solution ζ ( r ) of (13) vanishes at least twice on the interval [ R c , R d ] .
Hence, according to [27] Chapter 11, the solution η ( r ) of the Jacobi equation that verifies η ( R d ) = 0 and η ( R d ) = 1 vanishes at least twice on the interval [ R c , R d ] . This means that the Jacobi condition for the weak minimum value is not satisfied (see [27]).
Theorem 2.
If for  n = R d R c > 1  the following inequality holds:
n 1 n 3 / 2 > 2 · π · 1 R d · γ 1 / 2 · sin 3 / 2 α g ρ 1 / 2 · g 1 / 2
 then the liquid bridge is unstable.

3. Liquid Bridge in Thin Plate Growth

The prototype of a YOZ plane-symmetric plate growth system using the μ-PD technique is represented in Figure 9.
In this case, the energy function of the melt column behind the liquid bridge-free surface in terms of the liquid bridge profile curve z x is the following:
I z = x c x d γ · 1 + z 2 1 2 1 2 · ρ · g · z 2 p · z d x
with p = P a + ρ · g · h m (see [28]).
The stability condition of the YOZ plane-symmetric static liquid bridge is that z ( x ) minimizes (15). The first-order condition of the minimum is the Euler equation
d d x ( F z ) F z = 0
for the following function:
F ( z , z , x ) = γ · 1 + z 2 1 2 1 2 · ρ · g · z 2 p · z
Computing the Euler equation in this case for the liquid bridge profile curve, the following Equation is found:
z = ρ · g · z + p γ 1 + z 2 3 / 2 ; x c x x d
The function z ( x ) describing the profile curve has to verify the following boundary conditions:
z x c = tan π / 2 α g ;   z ( x c ) = h c > 0
z x d = tan α c ;   z x d = 0
z ( r )   is   strictly   decreasing   on   x c , x d
Using (16)–(19) and the concavity condition for z ( x ) ( z ( x ) < 0 ), the following inequalities can be found:
γ · α c + α g π 2 x d · n n 1 · cos α c ρ · g x d · n 1 n · tan α c p γ · α c + α g π 2 x d · n n 1 · sin α g
where n = x d x c > 1 . Inequalities (20) represent necessary conditions for the existence of z ( x ) . The pressure difference p = P a + ρ · g · h m has to verify (20). For stability or instability of the liquid bridge, the second-order conditions of the minimum of (16), i.e., the Legendre condition
2 F z z 0
and the Jacobi equation
d d x [ 2 F z z · η ] [ 2 F z z d d x ( 2 F z z ) ] · η = 0
have to be analyzed.
The Legendre condition in this case is the following inequality:
γ [ 1 + ( z ) 2 ] 3 2 0
which is obviously verified.
The Jacobi equation in this case is the following equation:
( γ [ 1 + ( z ) 2 ] 3 2 · η ) + ρ · g · η = 0
We note that the following inequalities concerning the coefficients of Equation (22) hold:
γ [ 1 + ( z ) 2 ] 3 2 γ · cos 3 α c ;   ρ · g ρ · g
Therefore, the following equation:
[ γ · cos 3 α c ] · ξ + [ ρ · g ] · ξ = 0
is a “Sturm-type upper bound” for Equation (22). Hence, the conclusion follows:
Theorem 3.
If for  n = x d x c > 1  the following inequality holds:
n 1 n < π · 1 x d · γ 1 / 2 · cos 3 / 2 α c ρ 1 / 2 · g 1 / 2
 then the liquid bridge is stable.
Analogously, the following equation:
[ γ · sin 3 α g ] · ζ + [ ρ · g ] · ζ = 0
is a “Sturm-type lower bound” for Equation (23). Hence, the conclusion follows:
Theorem 4.
If for  n = x d x c > 1  the following, inequality holds:
n 1 n > 2 · π · 1 x d · γ 1 / 2 · sin 3 / 2 α g ρ 1 / 2 · g 1 / 2
 then the liquid bridge is unstable.

4. Liquid Bridge in Thin Hollow Tube Growth

The prototype of an axis-symmetric hollow tube growth system using the μ-PD technique is represented in Figure 10.
In the case of the hollow tube growth method, there exist two independent liquid bridges: one in the interior of the tube and one on the exterior of the tube. The meridian curve of the bridge located in the interior of the tube is denoted by z i r , and the meridian curve of the bridge located on the exterior of the tube is denoted by z e ( r ) . In terms of these meridian curves, the functional energies of the melt columns behind the liquid bridge-free surface area are as follows (see [29]):
I i z i = R g i R c i γ · 1 + z i 2 1 2 1 2 · ρ · g · z i 2 p i · z i · r · d r
I e z e = R c e R g e γ · 1 + z e 2 1 2 1 2 · ρ · g · z e 2 p e · z e · r · d r
The stability condition of the axis-symmetric static liquid bridges is that z i ( x ) minimizes (28) and z e ( x ) minimizes (29).
The first-order necessary condition of the minimum is the Euler equation
d d x ( F z ) F z = 0
for the following functions:
F i ( z i , z i , r ) = γ · 1 + z i 2 1 2 1 2 · ρ · g · z i 2 p i · z i · r
F e ( z e , z e , r ) = γ · 1 + z e 2 1 2 1 2 · ρ · g · z e 2 p e · z e · r
with p i = P i a + ρ · g · h and p e = P e a + ρ · g · h .
Computing the Euler equation in this case for the interior and exterior liquid bridge profile curves, the following equations are found:
z i = ρ · g · z i + p i γ 1 + z i 2 3 2 1 r · 1 + z i 2 · z i ; R g i r R c i
z e = ρ · g · z e + p e γ 1 + z e 2 3 2 1 r · 1 + z e 2 · z e ;   R c e r R g e
Here, R g i , R g e represent the inner and outer radii of the shaper, respectively; R c i ,   R c e represent the inner and outer radii of the tube being grown, respectively, and satisfy the following:
0 < R g i < R c i < ( R g i + R g e ) / 2 < R c e < R g e
The function z i r describes the interior bridge meridian curve profile and verifies the following boundary conditions:
z i R c i = tan π / 2 α g ;   z i ( R c i ) = h i c > 0 .
z i R g i = tan α c ;   z i R g i = 0 .
z   is   strictly   increasing   on   [ R g i ,   R c i ]
The function z e r describes the exterior bridge meridian curve profile and verifies the following boundary conditions:
z e R c e = tan π / 2 α g ;   z e ( R c e ) = h e c > 0 .
z e R g e = tan α c ;   z e R g e = 0
z   is   strictly   decreasing   on   R c e , R g e .
Using (30), (32)–(34), and the concavity condition for z i ( r ) ( z i ( r ) < 0 ), the following inequalities can be found:
γ · α c + α g π / 2 R g i · 1 m 1 · cos α c ρ · g · R g i · ( m 1 ) · tan α c γ R g i · sin α c p i γ · α c + α g π / 2 R g i · 1 m 1 · sin α g γ R g i · 1 m · cos α g
where R c i / R g i = m > 1 .
Using (30), (35)–(37), and the concavity condition for z e ( r ) ( z e ( r ) < 0 ), the following inequalities can be found:
γ · α c + α g π 2 R g e · n n 1 · cos α c ρ · g R g e · n 1 n · tan α c + γ R g e · cos α g p e γ · α c + α g π 2 R g e · n n 1 · sin α g + γ R g e · n · sin α c
where n = R g e / R c e > 1 .
Inequalities (38) and (39) represent necessary conditions for the existence of z i ( r ) and z e ( r ) . The pressure differences p i = P i a + ρ · g · h and p e = P e a + ρ · g · h have to verify these inequalities. For stability or instability of the interior and exterior liquid bridge, the second-order conditions of the minima of (28) and (29), i.e., the Legendre condition
2 F z z 0
and the Jacobi equation
d d r [ 2 F z z · η ] [ 2 F z z d d r ( 2 F z z ) ] · η = 0
have to be analyzed for the following functions:
F i ( z i , z i , r ) = γ · 1 + z i 2 1 2 1 2 · ρ · g · z i 2 p i · z i · r
F e ( z e , z e , r ) = γ · 1 + z e 2 1 2 1 2 · ρ · g · z e 2 p e · z e · r
The Legendre conditions for the interior and exterior bridges are the following:
2 F i z i   z i = γ · [ 1 + ( z i ) 2 ] 3 / 2 · r 0 , r [ R g i ,   R c i ]
2 F e z e z e = γ · [ 1 + ( z e i ) 2 ] 3 / 2 · r 0 , r [ R c e ,   R g e ]
We note that these conditions are obviously verified.
The Jacobi equations for the interior and exterior liquid bridges are as follows:
γ · d d r [ 1 + ( z i ) 2 ] 3 / 2 · r · ξ ] + ρ · g · r · ξ = 0   r [ R g i ,   R c i ]
γ · d d r [ 1 + ( z e ) 2 ] 3 / 2 · r · ζ ] + ρ · g · r · ζ = 0   r R c e ,   R g e
For the Jacobi equation in the case of the interior liquid bridge (40), the following Equation (42):
γ · d d r [ R g i · cos 3 α c · ξ ] + ρ · g · R c i · ξ = 0
is a Sturm-type upper bound, and the Equation (43)
γ · d d r [ R c i · sin 3 α g · ζ ] + ρ · g · R g i · ζ = 0
is a Sturm-type lower bound. Therefore, for the inner liquid bridge, the following conclusion holds:
Theorem 5.
If for  m = R c i / R g i  the following, inequality holds:
( m 1 ) · m 1 / 2 < π R g i · γ 1 / 2 cos α c   3 / 2 ρ 1 / 2 · g 1 / 2
 then the internal liquid bridge is stable.
If for m = R c i / R g i the following inequality holds:
( m 1 ) · m 1 / 2 > 2 π R g i · γ 1 / 2 sin 3 / 2 α g ρ 1 / 2 · g 1 / 2
then the internal liquid bridge is unstable.
In the case of the exterior liquid bridge for the Jacobi Equation (41), the following Equation (46):
γ · d d r [ R c e · cos 3 α c · ξ ] + ρ · g · R g e · ξ = 0
is a Sturm-type upper bound, and the following Equation (47):
γ · d d r [ R g e · sin 3 α g · ζ ] + ρ · g · R c e · ζ = 0
is a Sturm-type lower bound. Therefore, for the exterior liquid bridge, the following conclusion holds:
Theorem 6.
If for  n = R g e / R c e  the following, inequality holds:
( n 1 ) · n 1 / 2 < π R g e · γ 1 / 2 cos α c   3 / 2 ρ 1 / 2 · g 1 / 2
 then the external liquid bridge is stable.
If for n = R g e / R c e the following inequality holds:
( n 1 ) · n 3 / 2 > 2 π R g e · γ 1 / 2 sin α g   3 / 2 ρ 1 / 2 · g 1 / 2
then the external liquid bridge is unstable.

5. Results

In the case of thin fibers for liquid bridges, the following results hold:
  • If, for n = R d R c > 1 the inequality n 1 n 1 / 2 < π · 1 R d · γ 1 / 2 · cos 3 / 2 α c ρ 1 / 2 · g 1 / 2 holds, then the liquid bridge is stable.
  • If, for n = R d R c > 1 the inequality n 1 n 3 / 2 > 2 · π · 1 R d · γ 1 / 2 · sin 3 / 2 α g ρ 1 / 2 · g 1 / 2 holds, then the liquid bridge is unstable.
In the case of thin plates for liquid bridges, the following results hold:
  • If for n = x d x c > 1 the inequality n 1 n < π · 1 x d · γ 1 / 2 · cos 3 / 2 α c ρ 1 / 2 · g 1 / 2 holds, then the liquid bridge is stable.
  • If for n = x d x c > 1 the inequality n 1 n > 2 · π · 1 x d · γ 1 / 2 · sin 3 / 2 α g ρ 1 / 2 · g 1 / 2 holds, then the liquid bridge is unstable.
In the case of thin hollow tubes for liquid bridges, the following results hold:
  • If for m = R c i / R g i the inequality ( m 1 ) · m 1 / 2 < π R g i · γ 1 / 2 cos α c   3 / 2 ρ 1 / 2 · g 1 / 2 holds, then the internal liquid bridge is stable.
  • If for m = R c i / R g i the inequality ( m 1 ) · m 1 / 2 > 2 π R g i · γ 1 / 2 sin 3 / 2 α g ρ 1 / 2 · g 1 / 2 holds, then the internal liquid bridge is unstable.
  • If for n = R g e / R c e the inequality ( n 1 ) · n 1 / 2 < π R g e · γ 1 / 2 cos α c   3 / 2 ρ 1 / 2 · g 1 / 2 holds, then the external liquid bridge is stable.
  • If for n = R g e / R c e the inequality ( n 1 ) · n 3 / 2 > 2 π R g e · γ 1 / 2 sin α g   3 / 2 ρ 1 / 2 · g 1 / 2 holds, then the external liquid bridge is unstable.

6. Comments and Conclusions

(i)
In the case of materials for which αcπ/2 − αg in papers [26,28,29,30,31], sufficient conditions for the theoretical existence and static stability of liquid bridges in thin fiber, thin plate, and thin hollow tube growth by the PD method were reported. If the sufficient conditions for stability of liquid bridges reported in [26,28,29] are not verified, this does not imply that the liquid bridge is unstable. As far as we know, there are no instability results reported in the literature for the αcπ/2 − αg case.
In the case of materials for which αcπ/2 − αg in [32,33,34], the dynamics of the growth process of thin fiber, thin plate, and thin hollow tube are described when the liquid bridge is stable. In these papers, the flow field, the thermal field, and the impurity distributions are also computed and numerical simulations are provided.
(ii)
In the case of materials for which α c > π / 2 α g , our findings reveal sufficient conditions for the stability or instability of a static liquid bridge appearing in the case of thin fiber, thin plate, and thin hollow tube growth by the pulling-down method. The imposed boundary conditions originate in crystal growth. If the reported sufficient condition of stability of the liquid bridge is not verified, this does not imply that the liquid bridge is unstable. The main novelty in this article consists in the obtained inequalities. These represent limits for what can and cannot be achieved. Experimentally, only stable static liquid bridges can be created. Unstable static liquid bridges exist just in theory; in reality, they collapse; therefore, they are not appropriate for crystal growth.

Author Contributions

Conceptualization, A.V.C. and S.B.; Formal analysis, S.B.; Writing—original draft, A.V.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kostoglou, M.; Karapantsios, T.D. On the identification of liquid surface properties using liquid bridges. Adv. Colloid Interface Sci. 2015, 222, 436–445. [Google Scholar] [CrossRef] [PubMed]
  2. Roisman, V.; Abboud, M.; Brockmann, P.; Berner, F.; Berger, R.; Brumm, P.R.; Sauer, H.M.; Dörsam, E.; Hussong, J. Forced flows in liquid bridges. Curr. Opin. Colloid Interface Sci. 2023, 67, 101738. [Google Scholar] [CrossRef]
  3. Gómez, F.G.; Flórez, V.M.; Morales, M.; Blanco, A.; López, C. Colloidal crystals and water: Perspectives on liquid–solid nanoscale phenomena in wet particulate media. Adv. Colloid Interface Sci. 2016, 234, 142–160. [Google Scholar] [CrossRef] [PubMed]
  4. Morawetz, K. Theory of water and charged liquid bridges. Phys. Rev. A At. Mol. Opt. Phys. 2012, 86, 026302. [Google Scholar] [CrossRef] [PubMed]
  5. Peregrine, D.H.; Shoker, G.; Symon, A. The Bifurcation of Liquid Bridges; Cambridge University Press: Cambridge, UK, 2006. [Google Scholar]
  6. Huang, C.H.; Carvalho, M.S.; Kumar, S. Stretching liquid bridges with moving contact lines: Comparison of liquid-transfer predictions and experiments. Soft Matter 2016, 12, 7457–7469. [Google Scholar] [CrossRef] [PubMed]
  7. Mastrangeli, M. The Fluid Joint: The Soft Spot of Micro- and Nanosystems; WILEY-VCH: Hoboken, NJ, USA, 2015. [Google Scholar] [CrossRef]
  8. Meseguer, J.; Slobozhanin, L.A.; Perales, J.M. A Review on the Stability of Liquid Bridges. Adv. Space Res. 1995, 16, 5–14. [Google Scholar] [CrossRef]
  9. Dhanaraj, G.; Byrappa, K.; Prasad, V.; Dudley, M. (Eds.) Springer Handbook of Crystal Growth; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  10. Tatarchenko, V.A. Shaped Crystal Growth; Kluwer Academic Publishers: Philip Drive Norwell, MA, USA, 1993. [Google Scholar]
  11. Wilcox, W.R.; Regel, L.L. Detached solidification in microgravity—A review. Microgravity Sci. Technol. 1998, 11, 152–166. [Google Scholar]
  12. Duffar, T. Bulk crystal growth under microgravity conditions. In Bulk Crystal Growth of Electronic, Optical and Optoelectronics Materials; Capper, P., Ed.; Wiley: Chichester, UK, 2005; pp. 477–524. [Google Scholar]
  13. Duffar, T.; Boiton, P.; Dusserre, P.; Abadie, J. Crucible de-wetting during Bridgman growth in microgravity. II. Smooth crucibles. J. Cryst. Growth 1997, 179, 397–409. [Google Scholar] [CrossRef]
  14. Duffar, T.; Sylla, L. Vertical Bridgman and dewetting. In Crystal Growth Processes Based on Capillarity; Duffar, T., Ed.; Wiley-Blackwell: Hoboken, NJ, USA, 2010; pp. 355–408. ISBN 978-0-470-71244-3. [Google Scholar]
  15. Sylla, L.; Paulin, J.P.; Vian, G.; Garnier, C.; Duffar, T. Effect of residual impurities on the dewetting of antimonide melts in fused silica crucibles in the case of bulk crystal growth. Mater. Sci. Eng. A 2008, 495, 208–214. [Google Scholar] [CrossRef]
  16. Sylla, L. Etude Expérimentale et Thermodynamique du Pocédé de Démouillage Appliqué aux Semiconducteurs. Ph.D. Thesis, Grenoble Institute of Technology, Saint-Martin-d’Hères, France, 2008. (In French). [Google Scholar]
  17. Balint, S.; Epure, S.; Duffar, T.; Braescu, L. Dewetted Bridgman crystal growth: Practical stability over a bounded time period in a forced regime. J. Eng. Math. 2012, 75, 191–208. [Google Scholar] [CrossRef]
  18. Balint, S.; Epure, S. Non-Lyapunov Type Stability in a Model of the Dewetted Bridgman Crystal Growth under Zero Gravity Conditions. Microgravity Sci. Technol. 2011, 23 (Suppl. S1), S49–S58. [Google Scholar] [CrossRef]
  19. Fukuda, T.; Chani, V.I. (Eds.) Shaped Crystals: Growth by Micro-Pulling-down Technique. Springer: New York, NY, USA, 2007. [Google Scholar]
  20. Yoshikawa, A.; Chani, V.I. Growth of Optical Crystals by the Micro-Pulling-down Method. MRS Bull. 2009, 34, 266–270. [Google Scholar] [CrossRef]
  21. Fukuda, T.; Rudolf, R.; Uda, S. (Eds.) Fiber Crystal Growth from the Melt; Springer: New York, NY, USA, 2004. [Google Scholar] [CrossRef]
  22. Zhuravleva, M.; Chani, V.I.; Yanagida, T.; Yoshikawa, A. The micro-pulling-down growth of Bi4Si3O12 (BSO) and Bi4Ge3O12 (BGO) fiber crystals and their scintillation efficiency. J. Crystal Growth 2008, 310, 2152–2156. [Google Scholar] [CrossRef]
  23. Yokota, Y.; Chani, V.I.; Sato, M.; Tota, K.; Onodera, K.; Yanagida, T.; Yoshikawa, A. Growth and crystallinity of shaped and multiple sapphire crystals by a micro-pulling-down method. J. Crystal Growth 2011, 318, 983–986. [Google Scholar] [CrossRef]
  24. Yokota, Y.; Tanaka, H.; Sato, M.; Chani, V.; Tota, K.; Onodera, K.; Yanagida, T.; Yoshikawa, A. Growth and Characterization of Shape-Controlled Single Crystals by a Micro-pulling-down Method. MRS Proc. 2011, 1309, 33–38. [Google Scholar] [CrossRef]
  25. Landau, L.D.; Lifchits, E.M. Mecanique de Fluides; Mir: Moscow, Russia, 1971; p. 132. [Google Scholar]
  26. Balint, A.M.; Balint, S.; Tanasie, L. An axis symmetric 2D description of a single crystal micro fiber growth process from the melt by micro-pulling-down method. Part 1. Math. Eng. Sci. Aerosp. 2019, 10, 759–770. [Google Scholar]
  27. Hartman, P. Ordinary Differential Equations; John Wiley & Sons: New York, NY, USA, 1964. [Google Scholar]
  28. Balint, A.M.; Stefan, B. A 2D description of the single crystal thin plate growth from the melt by micro-Pulling-down method. Part 1. INCAS Bull. 2018, 10, 37–52. [Google Scholar] [CrossRef]
  29. Balint, A.M.; Balint, S.; Tanasie, L. An axis symmetric 2D description of the growth process of a single crystal hollow tube, grown from the melt by pulling down method. Part. 1. Math. Eng. Sci. Aerosp. 2019, 10, 143–162. [Google Scholar]
  30. Balint, A.M.; Balint, S. Existence and Stability of the Solution of a Nonlinear Boundary Value Problem. Abstr. Appl. Anal. 2012, 2012, 582746. [Google Scholar] [CrossRef]
  31. Balint, A.M.; Stefan, B. Dynamics and control of the shape and size of a sitting drop like meniscus, occurring in single crystal growth. J. Frankl. Inst. 2010, 347, 30–52. [Google Scholar] [CrossRef]
  32. Balint, A.M.; Balint, S.; Tanasie, L. An axis symmetric 2D description of a single crystal micro fiber growth process from the melt by micro-pulling-down method. Part 2. Math. Eng. Sci. Aerosp. 2019, 10, 771–780. [Google Scholar]
  33. Balint, A.M.; Stefan, B. A 2D description of the single crystal thin plate growth from the melt by micro-Pulling-down method. Part 2. INCAS Bull. 2018, 10, 53–64. [Google Scholar] [CrossRef]
  34. Balint, A.M.; Balint, S.; Tanasie, L. An axis symmetric 2D description of the growth process of a single crystal hollow tube, grown from the melt by pulling down method. Part. 2. Math. Eng. Sci. Aerosp. 2019, 10, 163–174. [Google Scholar]
Figure 1. Concave capillary bridge between two planes (schematic representation).
Figure 1. Concave capillary bridge between two planes (schematic representation).
Fluids 09 00176 g001
Figure 2. Cylindrical bar BC growth.
Figure 2. Cylindrical bar BC growth.
Fluids 09 00176 g002
Figure 3. Ribbon EFG growth, case α c > π / 2 α g .
Figure 3. Ribbon EFG growth, case α c > π / 2 α g .
Fluids 09 00176 g003
Figure 4. Ribbon EFG growth, case αcπ/2 − αg.
Figure 4. Ribbon EFG growth, case αcπ/2 − αg.
Fluids 09 00176 g004
Figure 5. Cylindrical bar EFG growth, case αcπ/2 − αg.
Figure 5. Cylindrical bar EFG growth, case αcπ/2 − αg.
Fluids 09 00176 g005
Figure 6. Hollow tube EFG growth, case αcπ/2 − αg.
Figure 6. Hollow tube EFG growth, case αcπ/2 − αg.
Fluids 09 00176 g006
Figure 7. Schematic of the dewetted Bridgman technique.
Figure 7. Schematic of the dewetted Bridgman technique.
Fluids 09 00176 g007
Figure 8. Schematic presentation of an axis-symmetric fiber growth system using the μ-PD technique.
Figure 8. Schematic presentation of an axis-symmetric fiber growth system using the μ-PD technique.
Fluids 09 00176 g008
Figure 9. Schematic presentation of a YOZ-plane symmetric plate growth system using the μ-PD method.
Figure 9. Schematic presentation of a YOZ-plane symmetric plate growth system using the μ-PD method.
Fluids 09 00176 g009
Figure 10. Schematic presentation of an axis-symmetric hollow tube growth system using the μ-PD method.
Figure 10. Schematic presentation of an axis-symmetric hollow tube growth system using the μ-PD method.
Fluids 09 00176 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cojocaru, A.V.; Balint, S. Stability or Instability of a Static Liquid Bridge Appearing in Shaped Crystal Growth from Melt via the Pulling-Down Method. Fluids 2024, 9, 176. https://doi.org/10.3390/fluids9080176

AMA Style

Cojocaru AV, Balint S. Stability or Instability of a Static Liquid Bridge Appearing in Shaped Crystal Growth from Melt via the Pulling-Down Method. Fluids. 2024; 9(8):176. https://doi.org/10.3390/fluids9080176

Chicago/Turabian Style

Cojocaru, Andreea V., and Stefan Balint. 2024. "Stability or Instability of a Static Liquid Bridge Appearing in Shaped Crystal Growth from Melt via the Pulling-Down Method" Fluids 9, no. 8: 176. https://doi.org/10.3390/fluids9080176

APA Style

Cojocaru, A. V., & Balint, S. (2024). Stability or Instability of a Static Liquid Bridge Appearing in Shaped Crystal Growth from Melt via the Pulling-Down Method. Fluids, 9(8), 176. https://doi.org/10.3390/fluids9080176

Article Metrics

Back to TopTop