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

: This study presents sufﬁcient conditions for the stability or instability of the static liquid bridge appearing in crystal growth from the melt of micro-ﬁbers, 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 signiﬁcant for shaped crystal growth from melted materials, with given macroscopic dimensions, and using speciﬁc equipment. This is because the obtained inequalities represent limits for what can and cannot be achieved experimentally.


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.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.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-filmfed-growth (EFG) method.For the BC method see Figure 2, and for the EFG method see Figures 3-6.
Fluids 2024, 9, x FOR PEER REVIEW Despite the plethora of experimental and theoretical studies on liquid bridges, it prising that very few refer to their stability and instability, which factor in shaped crystal growth from melt material via the edge-defined-film-fed-growth (EFG) and pulling-down (μ-PD) methods.The liquid bridge is a small amount of melt connec shaper and the crystal being grown, and in crystal growth this is called the menis The basic growth methods available for crystal growth are, broadly, growth f melt, growth from vapor, and growth from solution.The Springer Handbook of Growth [9] presents in 1816 pages the state of the art in crystal growth until the ye In Chapter 40 pg.1379-1402 of this book, the authors Th.George, St. Balint, and L. B present several mathematical models describing processes that take place in the crystal growth from the melt by Bridgman-Stockbarger (BC) and by the edge-d film-fed-growth (EFG) method.For the BC method see Figure 2, and for the EFG see Figures 3-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, Fluids 2024, 9, 176 5 of 17 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 (Figures 3-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. 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 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.

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.
uid bridge; our research excludes the study of the growth of shaped single crystals.this reason, the reference list is short in comparison with the existing literature, and thermal aspect is completely disregarded.

Liquid Bridge in Micro-Fiber Growth
The prototype of an axis-symmetric fiber growth system using the μ-PD techniqu represented in Figure 8.In the hydrostatic approximation, the equation describing the equilibrium capill surface (shape of the liquid bridge surface) is the Laplace capillary equation (see [25 132).However, to calculate the meniscus shape for the capillary-shaping techniques, convenient to employ the Laplace equation in its explicit differential form.To obtain t form of the equation following [25] (p.132), the variation problem of the total free ene minimum of a liquid column will be solved.For our system, the total free energy is co posed of the surface free energy and the gravity field energy.The minimum should calculated with the additional condition of liquid volume constancy.Solving this probl presents some difficulties.Therefore, our study will be restricted to considering a men 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: 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 for the following function: In computing the Euler equation for the meridian curve, the following equation is found: The function z(r) describes the meridian curve profile and verifies the following conditions: Fluids 2024, 9, 176 8 of 17 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.
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 R d .We point out two things: (a) In Figure 3, a concave meridian curve is represented; in Figures 4-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: 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 have to be analyzed.The Legendre condition in this case is the inequality: Fluids 2024, 9, 176 9 of 17 which is obviously verified.The Jacobi equation in this case is the following: 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: It follows that the equation as follows: 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 The half period of any non-zero solution of ( 10) 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 In other words, if the following inequality holds , 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]).
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: It follows that the following equation: 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 . The period of any non-zero solution of ( 13) is . 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 holds, then any non-zero solution ζ(r) of ( 13) vanishes at Fluids 2024, 9, 176 10 of 17 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]).
then the liquid bridge is unstable.

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 ( ) x z is the following: (see [28]).The stability condition of the YOZ plane-symmetric static liquid bridge is that ) (x z minimizes (15).The first-order condition of the minimum is the Euler equation Computing the Euler equation in this case for the liquid bridge profile curve, the following Equation is found: 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: 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 for the following function: Computing the Euler equation in this case for the liquid bridge profile curve, the following Equation is found: The function z(x) describing the profile curve has to verify the following boundary conditions: Using ( 16)-( 19) and the concavity condition for z(x) (z (x) < 0), the following inequalities can be found: 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 have to be analyzed.The Legendre condition in this case is the following inequality: which is obviously verified.The Jacobi equation in this case is the following equation: We note that the following inequalities concerning the coefficients of Equation ( 22) hold: Therefore, the following equation: 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: Fluids 2024, 9, 176 12 of 17 then the liquid bridge is stable.
Analogously, the following equation: 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: then the liquid bridge is unstable.

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.Analogously, the following equation: is a "Sturm-type lower bound" for Equation (23).Hence, the conclusion follows: the following, inequality holds: then the liquid bridge is unstable.

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 .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]): 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]): Fluids 2024, 9, 176 13 of 17 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 for the following functions: Computing the Euler equation in this case for the interior and exterior liquid bridge profile curves, the following equations are found: Here, R gi , R ge represent the inner and outer radii of the shaper, respectively; R ci , R ce represent the inner and outer radii of the tube being grown, respectively, and satisfy the following: The function z i (r) describes the interior bridge meridian curve profile and verifies the following boundary conditions: The function z e (r) describes the exterior bridge meridian curve profile and verifies the following boundary conditions: Using ( 30), ( 32)- (34), and the concavity condition for z i (r) (z i (r) < 0), the following inequalities can be found: where R ci /R gi = m > 1.Using ( 30), ( 35)-(37), and the concavity condition for z e (r) (z e (r) < 0), the following inequalities can be found: where n = R ge /R ce > 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 a i + ρ • g • h and p e = −P a e + ρ • 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 and the Jacobi equation have to be analyzed for the following functions: The Legendre conditions for the interior and exterior bridges are the following: We note that these conditions are obviously verified.The Jacobi equations for the interior and exterior liquid bridges are as follows: For the Jacobi equation in the case of the interior liquid bridge (40), the following Equation (42): is a Sturm-type lower bound.Therefore, for the inner liquid bridge, the following conclusion holds: Theorem 5.If for m = R ci /R gi the following, inequality holds: then the internal liquid bridge is stable.
If for m = R ci /R gi the following inequality holds: then the internal liquid bridge is unstable.
In the case of the exterior liquid bridge for the Jacobi Equation (41), the following Equation (46): is a Sturm-type upper bound, and the following Equation (47): is a Sturm-type lower bound.Therefore, for the exterior liquid bridge, the following conclusion holds: Theorem 6.If for n = R ge /R ce the following, inequality holds: then the external liquid bridge is stable.
If for n = R ge /R ce the following inequality holds: then the external liquid bridge is unstable.

Results
In the case of thin fibers for liquid bridges, the following results hold: holds, then the liquid bridge is stable.
holds, then the liquid bridge is unstable.
In the case of thin plates for liquid bridges, the following results hold:

Figure 8 .
Figure 8. Schematic presentation of an axis-symmetric fiber growth system using the μ-PD te nique.

Figure 8 .
Figure 8. Schematic presentation of an axis-symmetric fiber growth system using the µ-PD technique.

Figure 9 .
Figure 9. Schematic presentation of a YOZ-plane symmetric plate growth system using the μ-PD method.

Figure 9 .
Figure 9. Schematic presentation of a YOZ-plane symmetric plate growth system using the µ-PD method.

Figure 10 .
Figure 10.Schematic presentation of an axis-symmetric hollow tube growth system using the μ-PD method.
meridian curve of the bridge located on the exterior of the tube is denoted by ) (r z e

Figure 10 .
Figure 10.Schematic presentation of an axis-symmetric hollow tube growth system using the µ-PD method.
Sturm-type upper bound, and the Equation (43)