Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2009 / Article

Research Article | Open Access

Volume 2009 |Article ID 348538 | 22 pages | https://doi.org/10.1155/2009/348538

On the Creation of a Stable Drop-Like Static Meniscus, Appropriate for the Growth of a Single Crystal Tube with Prior Specified Inner and Outer Radii

Academic Editor: Ben T. Nohara
Received22 Feb 2009
Revised06 May 2009
Accepted24 May 2009
Published29 Jul 2009

Abstract

A theoretical procedure for the creation of a stable drop-like static meniscus, appropriate for the growth of a single crystal tube, with a priori specified inner and outer radius, is presented. The method locates the controllable part 𝑝 of the pressure difference across the free surface. It consists in a set of calculus, which leads to the determination of the melt column height (between the horizontal crucible melt level and the shaper top level) in function of the pressure of the gas flow (introduced in the furnace for release the heat) in order to obtain the desired meniscus. The procedure is presented in general and is numerically illustrated for InSb tubes. The novelty is the algorithm for the exact determination of 𝑝, which has to be used, the determination of the melt column height, and the evaluation of the effect of shaper radii. The setting of the thermal conditions, which assure that for the obtained static meniscus the solidification conditions are satisfied at the β€œright” places, is not considered here.

1. Introduction

The conventional melt growth techniques, as Bridgman growth [1–3] or Czochralski pulling [4–6] of single crystals, typically produce ingots of circular or square cross-sections which need to be cut in hundreds of slices to produce wafers. Using these processes, it is difficult to produce thin wafers from an ingot without wasting 40%–50% of material as kerfs during the cutting process. For this reason the E.F.G. technology can be more appropriate to produce single crystals with prescribed shapes and sizes which can be used without additional machining.

The growth of silicon tubes by E.F.G. process was first reported by Erris et al. [7]. In [7] a theory of tube growth by E.F.G. process is developed to show the dependence of tube wall thickness on the growth variables. The theory uses approximation reported in [8, 9], and it has been shown to be a useful tool understanding the feasible limits of the wall thickness control. A more accurate predictive model would require an increase of the acceptable tolerance range introduced by approximation.

Later, the heat flow in a tube growth system was analyzed in [10–19].

The state of the arts at the time 1993-1994, concerning the calculation of the meniscus shape in general in the case of the growth by E.F.G. method is summarized in [20]. According to [20], for the general differential equation describing the free surface of a liquid meniscus, possessing axial symmetry, there are no complete analysis and solution. For the general equation only numerical integrations were carried out for a number of process parameter values that were of practical interest at the moment. The authors of [21, 22] consider automated crystal growth processes based on weight sensors and computers. They give an expression for the weight of the meniscus, contacted with crystal and shaper of arbitrary shape, in which there are two terms related to the hydrodynamic factor.

In [23] it is shown that the hydrodynamic factor is too small to be considered in the automated crystal growth. In [24] a theoretical and numerical study of meniscus dynamics, under symmetric and asymmetric configurations, is presented. A meniscus dynamics model is developed to consider meniscus shape and its dynamics, heat and mass transfer around the die top and meniscus. Analysis reveals the correlations among tube thickness, effective melt height, pull rate, die top temperature, and crystal environmental temperature.

In [25] the effect of the controllable part of the pressure difference on the free surface shape of the static meniscus is analyzed for the tube growth by E.F.G. method for materials for which 0<𝛼𝑐<πœ‹/2; 0<𝛼𝑔<πœ‹/2; 𝛼𝑐>πœ‹/2βˆ’π›Όπ‘”.

The present paper concerns also the shape and the stability of the free surface of a static meniscus (pulling rate equal to zero). More precisely, it is shown in which kind the explicit formulas reported in [25] can be combined in order to create a stable static meniscus having a free surface with prescribed size and shape, which is appropriate for the growth of a single crystal tube having a priori specified inner and outer radii. The free surface of a static meniscus is appropriate for the growth of a single crystal tube of constant inner radius π‘Ÿπ‘– and constant outer radius π‘Ÿπ‘’ if the angle between the tangent lines to the free surface at the points (π‘Ÿπ‘–,𝑧𝑖(π‘Ÿπ‘–)),(π‘Ÿπ‘’,𝑧𝑒(π‘Ÿπ‘’)) (Figure 1 where the solidification conditions have to be assured) and the vertical is equal to the growth angle 𝛼𝑔. Moreover, the function describing the free surface has to minimize the energy functional of the melt column (i.e., the meniscus has to be stable). In this paper we give a procedure for the choice of the melt column height, between the horizontal crucible melt level and shaper top level and of the pressure of the gas flow introduced in the furnace (for release the heat), in order to create a static meniscus of which free surface is appropriate for the growth of a single crystal tube of constant inner radius π‘Ÿπ‘– and outer radius π‘Ÿπ‘’.The thermal problem concerning the setting of the thermal conditions, which assure that for the obtained static meniscus at the level 𝑧𝑖(π‘Ÿπ‘–),𝑧𝑒(π‘Ÿπ‘’) the solidification conditions are satisfied is not considered in this paper. The novelty consists in the fact that the free surface is not approximated by an arc with constant curvature, the computation takes into account the pressure of the gas flow, and the stability of the free surface is assured.

2. The Free Surfaces Equations and the Pressure Difference Limits

For a single crystal tube growth by E.F.G. technique, in hydrostatic approximation, the outer free surface equation of the static meniscus is

π‘§π‘’ξ…žξ…ž=πœŒβ‹…π‘”β‹…π‘§π‘’βˆ’π‘π‘’π›Ύξ‚ƒξ€·π‘§1+ξ…žπ‘’ξ€Έ2ξ‚„3/2βˆ’1π‘Ÿβ‹…ξ‚ƒξ€·π‘§1+ξ…žπ‘’ξ€Έ2ξ‚„β‹…π‘§ξ…žπ‘’ξ‚Έπ‘…;π‘Ÿβˆˆπ‘”π‘–+𝑅𝑔𝑒2,𝑅𝑔𝑒,(2.1) and the inner free surface equation is π‘§π‘–ξ…žξ…ž=πœŒβ‹…π‘”β‹…π‘§π‘–βˆ’π‘π‘–π›Ύξ‚ƒξ€·π‘§1+ξ…žπ‘–ξ€Έ2ξ‚„3/2βˆ’1π‘Ÿβ‹…ξ‚ƒξ€·π‘§1+ξ…žπ‘–ξ€Έ2ξ‚„β‹…π‘§ξ…žπ‘–ξ‚Έπ‘…;π‘Ÿβˆˆπ‘”π‘–,𝑅𝑔𝑖+𝑅𝑔𝑒2ξ‚Ή.(2.2) Here, Ξ³ is the surface tension of the melt; ρ is the melt density; 𝑔 is the gravitational acceleration; 𝑧𝑒,𝑧𝑖 are the coordinates with respect to the Oz axis, directed vertically upwards; π‘Ÿ is the radial coordinate with respect to the π‘‚π‘Ÿ axis, oriented horizontal; 𝑅𝑔𝑒,𝑅𝑔𝑖 are the outer and inner radius of the shaper, respectively; 𝑝𝑒,𝑝𝑖 are the pressure difference across the outer and inner free surface, respectively:

𝑝𝑒=π‘π‘šβˆ’π‘π‘’π‘”βˆ’πœŒβ‹…π‘”β‹…π»,𝑝𝑖=π‘π‘šβˆ’π‘π‘–π‘”βˆ’πœŒβ‹…π‘”β‹…π».(2.3) In (2.3), π‘π‘š denotes the hydrodynamic pressure in the meniscus melt due to the thermal and Marangoni convection; 𝑝𝑒𝑔,𝑝𝑖𝑔 are the pressure of the gas flow introduced in the exterior and in the interior of the tube, respectively, for releasing the heat from the inner and outer side of the tube wall; 𝐻 denotes the melt column β€œheight” between the horizontal crucible melt level, and the shaper top level (Figure 1). 𝐻 is positive when the crucible melt level is under the shaper top level and it is negative when the shaper top level is under the crucible melt level.

The solution 𝑧𝑒=𝑧𝑒(π‘Ÿ) of (2.1) has to satisfy the following conditions:

π‘§ξ…žπ‘’ξ€·π‘Ÿπ‘’ξ€Έξ‚€πœ‹=βˆ’tan2βˆ’π›Όπ‘”ξ‚π‘§,(2.4a)ξ…žπ‘’ξ€·π‘…π‘”π‘’ξ€Έ=βˆ’tan𝛼𝑐𝑧,(2.4b)𝑒𝑅𝑔𝑒=0andπ‘§π‘’ξ€Ίπ‘Ÿ(π‘Ÿ)isstrictlydecreasingon𝑒,𝑅𝑔𝑒,(2.4c) where π‘Ÿπ‘’βˆˆ((𝑅𝑔𝑖+𝑅𝑔𝑒)/2,𝑅𝑔𝑒) is the tube outer radius; 𝛼𝑔 is the growth angle; 𝛼𝑐is the contact angle between the outer free surface and the outer edge of the shaper top and 0<𝛼𝑐<πœ‹/2;0<𝛼𝑔<πœ‹/2;πœ‹/2βˆ’π›Όπ‘”<𝛼𝑐 (Figure 1).

Condition (2.4a) expresses that at the point (π‘Ÿπ‘’,𝑧𝑒(π‘Ÿπ‘’)) (the left end of the outer free surface), where the solidification has to be realized, the angle between the tangent line to the free surface and the vertical is equal to the growth angle 𝛼𝑔 (i.e., the tangent to the tube outer wall is vertical).

Condition (2.4b) expresses that at the point (𝑅𝑔𝑒,0) (the right end of the outer free surface, where the free surface is attached to the shaper edge), the angle between the tangent line to the free surface and horizontal (i.e., the contact angle) is equal to 𝛼𝑐.

Condition (2.4c) expresses that at the point (𝑅𝑔𝑒,0) the free surface is attached to the shaper edge.

Moreover, the solution 𝑧𝑒=𝑧𝑒(π‘Ÿ) has to minimize the energy functional of the melt column: πΌπ‘’ξ€œ(𝑧)=π‘…π‘”π‘’π‘Ÿπ‘’ξ‚»ξ‚ƒξ€·π‘§π›Ύβ‹…1+ξ…žξ€Έ2ξ‚„1/2+12β‹…πœŒβ‹…π‘”β‹…π‘§2βˆ’π‘π‘’ξ‚Όπ‘§ξ€·π‘Ÿβ‹…π‘§β‹…π‘Ÿβ‹…π‘‘π‘Ÿ,𝑒=β„Žπ‘’ξ€·π‘…>0,𝑧𝑔𝑒=0.(2.5) The solution 𝑧𝑖=𝑧𝑖(π‘Ÿ) of (2.2) has to satisfy the following conditions: π‘§ξ…žπ‘–ξ€·π‘…π‘”π‘–ξ€Έ=tan𝛼𝑐𝑧,(2.6a)ξ…žπ‘–ξ€·π‘Ÿπ‘–ξ€Έξ‚€πœ‹=tan2βˆ’π›Όπ‘”ξ‚π‘§,(2.6b)𝑖𝑅𝑔𝑖=0and𝑧𝑖𝑅(π‘Ÿ)isstrictlyincreasingon𝑔𝑖,π‘Ÿπ‘–ξ€»,(2.6c) where: π‘Ÿπ‘–βˆˆ(𝑅𝑔𝑖,(𝑅𝑔𝑖+𝑅𝑔𝑒)/2) is the tube inner radius; 𝛼𝑔 is the growth angle; 𝛼𝑐is the contact angle between the inner free surface and the inner edge of the shaper top and 0<𝛼𝑐<πœ‹/2;0<𝛼𝑔<πœ‹/2;πœ‹/2βˆ’π›Όπ‘”<𝛼𝑐 (Figure 1).

Condition (2.6a) expresses that at the point (𝑅𝑔𝑖,0) (the left end of the inner free surface, where the free surface is attached to the shaper edge), the angle between the tangent line to the free surface and horizontal (i.e., the contact angle) is equal to 𝛼𝑐.

Condition (2.6b) expresses that at the point (π‘Ÿπ‘–,𝑧𝑖(π‘Ÿπ‘–)) (the right end of the inner free surface), where the solidification has to be realized, the angle between the tangent line to the free surface and the vertical is equal to the growth angle 𝛼𝑔 (i.e., the tangent to the tube inner wall is vertical).

Condition (2.6c) expresses that at the point (𝑅𝑔𝑖,0) the free surface is attached to the shaper edge.

Moreover, the solution 𝑧𝑖=𝑧𝑖(π‘Ÿ) has to minimize the energy functional of the melt column: πΌπ‘–ξ€œ(𝑧)=π‘Ÿπ‘–π‘…π‘”π‘–ξ‚»ξ‚ƒξ€·π‘§π›Ύβ‹…1+ξ…žξ€Έ2ξ‚„1/2+12β‹…πœŒβ‹…π‘”β‹…π‘§2βˆ’π‘π‘–ξ‚Όπ‘§ξ€·π‘…β‹…π‘§β‹…π‘Ÿβ‹…π‘‘π‘Ÿ,π‘”π‘–ξ€Έξ€·π‘Ÿ=0,𝑧𝑖=β„Žπ‘–>0.(2.7)

Based on the mathematical theorems, rigorously proven in [25], the following statements, regarding the creation of an appropriate meniscus, can be formulated.

Statement 1
If the solution of the initial value problem (IVP) π‘§π‘’ξ…žξ…ž=πœŒβ‹…π‘”β‹…π‘§π‘’βˆ’π‘π‘’π›Ύξ‚ƒξ€·π‘§1+ξ…žπ‘’ξ€Έ2ξ‚„3/2βˆ’1π‘Ÿβ‹…ξ‚ƒξ€·π‘§1+ξ…žπ‘’ξ€Έ2ξ‚„β‹…π‘§ξ…žπ‘’,𝑧𝑒𝑅𝑔𝑒=0,π‘§ξ…žπ‘’ξ€·π‘…π‘”π‘’ξ€Έ=βˆ’tan𝛼𝑐(2.8) is convex, then it does not represent the outer free surface of an appropriate drop-like static meniscus.

Comment 1. The above statement shows that the solutions of IVP (2.8), which represent the outer free surface of an appropriate drop-like static meniscus, can be obtained only for those values of 𝑝𝑒 for which the solution is not globally convex.

Statement 2
If 𝑝𝑒<(𝛾/𝑅𝑔𝑒)β‹…sin𝛼𝑐, then the solution of the IVP (2.8) is globally convex, and it does not represent the outer free surface of an appropriate drop-like static meniscus.

Comment 2. The above statement locates a set of 𝑝𝑒 values which are not appropriate for the creation of a drop-like meniscus.

Statement 3
If the solution 𝑧𝑒=𝑧𝑒(π‘Ÿ) of the nonlinear boundary value problem (NLBVP) (2.1), (2.4a), (2.4b), and (2.4c) represents the outer free surface of an appropriate concave static meniscus on the closed interval [𝑅𝑔𝑒/𝑛,𝑅𝑔𝑒], with 1<𝑛<(𝑅𝑔𝑖+𝑅𝑔𝑒)/2, then the following inequalities hold: π‘›π›Όπ‘›βˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑒⋅cos𝛼𝑐+𝛾𝑅𝑔𝑒⋅cosπ›Όπ‘”β‰€π‘π‘’β‰€π‘›π›Όπ‘›βˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑒⋅sin𝛼𝑔+π‘›βˆ’1π‘›β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘’β‹…tan𝛼𝑐𝛾+𝑛⋅𝑅𝑔𝑒⋅sin𝛼𝑐.(2.9)

Comment 3. The above statement shows that in order to obtain the outer free surface of a concave static meniscus, appropriate for the growth of a tube of outer radius π‘Ÿπ‘’=𝑅𝑔𝑒/𝑛, 𝑝𝑒 has to be searched in the range defined by the inequalities (2.9).

Consequence 1
If 𝑛=2⋅𝑅𝑔𝑒/(𝑅𝑔𝑖+𝑅𝑔𝑒) (π‘Ÿπ‘’ is the middle point of the interval [𝑅𝑔𝑖,𝑅𝑔𝑒]), then π‘Ÿπ‘’=𝑅𝑔𝑒/𝑛=(𝑅𝑔𝑖+𝑅𝑔𝑒)/2 and inequalities (2.9) become 𝛼2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…cos𝛼𝑐+𝛾𝑅𝑔𝑒⋅cos𝛼𝑔≀𝑝𝑒𝛼≀2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…sin𝛼𝑔+ξ€·π‘…πœŒβ‹…π‘”β‹…π‘”π‘’βˆ’π‘…π‘”π‘–ξ€Έ2β‹…tan𝛼𝑐+2⋅𝛾𝑅𝑔𝑒+𝑅𝑔𝑖⋅sin𝛼𝑐.(2.10)

Comment 4. The above consequence shows that in order to obtain the outer free surface of a concave static meniscus, appropriate for the growth of a tube of outer radius π‘Ÿπ‘’=(𝑅𝑔𝑖+𝑅𝑔𝑒)/2, 𝑝𝑒 has to be searched in the range defined by the inequalities (2.10).

Consequence 2
If 𝑝𝑒 verifies the inequality 𝑝𝑒𝛼<2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…cos𝛼𝑐+𝛾𝑅𝑔𝑒⋅cos𝛼𝑔,(2.11) then there is no π‘Ÿπ‘’ in the closed interval [(𝑅𝑔𝑖+𝑅𝑔𝑒)/2,𝑅𝑔𝑒] for which the NLBVP (2.1), (2.4a), (2.4b), and (2.4c) possesses a concave solution.

Comment 5. The above consequence shows that if 𝑝𝑒 is in the range defined by the inequality (2.11), then it is impossible to obtain a static meniscus having concave outer free surface which is appropriate for the growth of a tube of outer radius π‘Ÿπ‘’ situated in the range [(𝑅𝑔𝑖+𝑅𝑔𝑒)/2,𝑅𝑔𝑒].

Consequence 3
If 𝑛→1, then π‘Ÿπ‘’=𝑅𝑔𝑒/𝑛→𝑅𝑔𝑒 and, according to (2.9), 𝑝𝑒→+∞.

Comment 6. The above consequence shows that in order to obtain the outer free surface of a concave static meniscus, appropriate for the growth of a tube of outer radius π‘Ÿπ‘’β‰ˆπ‘…π‘”π‘’, 𝑝𝑒 has to be very high.

Statement 4
If 𝑛 and 𝑝𝑒 verify the inequalities 1<𝑛<2⋅𝑅𝑔𝑒𝑅𝑔𝑒+𝑅𝑔𝑖,𝑝𝑒>π‘›π›Όπ‘›βˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑒⋅sin𝛼𝑔+π‘›βˆ’1π‘›β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘’β‹…tan𝛼𝑐𝛾+𝑛⋅𝑅𝑔𝑒⋅sin𝛼𝑐,(2.12) then there exist π‘Ÿπ‘’βˆˆ[𝑅𝑔𝑒/𝑛,𝑅𝑔𝑒] and a concave solution of the NLBVP (2.1), (2.4a), (2.4b), and (2.4c) on the interval [π‘Ÿπ‘’,𝑅𝑔𝑒].

Comment 7. The above consequence shows that if 𝑛 and 𝑝𝑒 verify (2.12), then a static meniscus having concave outer free surface, appropriate for the growth of a tube of outer radius π‘Ÿπ‘’ situated in the range [𝑅𝑔𝑒/𝑛,𝑅𝑔𝑒], is obtained.

Consequence 4
If for 1<π‘›ξ…ž<𝑛<2⋅𝑅𝑔𝑒/(𝑅𝑔𝑖+𝑅𝑔𝑒) and 𝑝𝑒 the following inequalities hold: π‘›π›Όπ‘›βˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑒⋅sin𝛼𝑔+π‘›βˆ’1π‘›β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘’β‹…tan𝛼𝑐𝛾+𝑛⋅𝑅𝑔𝑒⋅sin𝛼𝑐<𝑝𝑒<π‘›ξ…žπ‘›ξ…žπ›Όβˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑒⋅cos𝛼𝑐+𝛾𝑅𝑔𝑒⋅cos𝛼𝑔,(2.13) then there exist π‘Ÿπ‘’ in the interval [𝑅𝑔𝑒/𝑛,𝑅𝑔𝑒/π‘›ξ…ž] and a concave solution of the NLBVP (2.1), (2.4a), (2.4b), and (2.4c) on the interval [π‘Ÿπ‘’,𝑅𝑔𝑒].

Comment 8. The above consequence shows that if 𝑝𝑒 is in the range defined by the inequalities (2.13), then a static meniscus having concave outer free surface, appropriate for the growth of a tube of outer radius π‘Ÿπ‘’ situated in the range [𝑅𝑔𝑒/𝑛,𝑅𝑔𝑒/π‘›ξ…ž], is obtained.

Statement 5
A concave solution 𝑧𝑒=𝑧𝑒(π‘Ÿ) of the NLBVP (2.1), (2.4a), (2.4b), and (2.4c) is a weak minimum of the energy functional of the melt column.

Comment 9. The above consequence shows that a static meniscus having a concave outer free surface, appropriate for the growth of a tube of outer radius π‘Ÿπ‘’ situated in the range [(𝑅𝑔𝑖+𝑅𝑔𝑒)/2,𝑅𝑔𝑒], is stable.

Statement 6
If 𝑛 and 𝑝𝑒 satisfy the inequalities 1<𝑛<2⋅𝑅𝑔𝑒𝑅𝑔𝑒+𝑅𝑔𝑖,𝑝𝑒>π‘›βˆ’1π‘›β‹…π‘”β‹…πœŒβ‹…π‘…π‘”π‘’β‹…tan𝛼𝑐𝛾+𝑛⋅𝑅𝑔𝑒,(2.14) then the solution 𝑧𝑒=𝑧𝑒(π‘Ÿ) of the IVP (2.8) is concave on the interval 𝐼∩[𝑅𝑔𝑒/𝑛,𝑅𝑔𝑒] where 𝐼 is the maximal interval of the existence of 𝑧𝑒(π‘Ÿ).

Comment 10. The above consequence shows that if 𝑝𝑒 is in the range defined by the inequalities (2.14), then eventually it can be used for the creation of a drop-like static meniscus, appropriate for the growth of a tube of outer radius π‘Ÿπ‘’ in the range [𝑅𝑔𝑒/𝑛,𝑅𝑔𝑒].

Statement 7
If 𝑝𝑒>(𝛾/𝑅𝑔𝑒)β‹…sin𝛼𝑐 and for a value π‘Ÿπ‘’=𝑅𝑔𝑒/𝑛, which satisfies π‘Ÿπ‘’βˆˆ((𝑅𝑔𝑖+𝑅𝑔𝑒)/2,𝑅𝑔𝑒), a static meniscus appropriate for the growth of a tube of outer radius π‘Ÿπ‘’ exists, then for 𝑝𝑒 the following inequalities hold: 𝛾𝑅𝑔𝑒⋅sin𝛼𝑐<𝑝𝑒<πœŒβ‹…π‘”β‹…π‘›βˆ’1𝑛⋅𝑅𝑔𝑒⋅tan𝛼𝑐𝛾+𝑛⋅𝑅𝑔𝑒⋅cos𝛼𝑐.(2.15)

Comment 11. The above statement locates those values of 𝑝𝑒 for which eventually nonglobally concave (drop-like) static meniscus, appropriate for the growth of a tube of outer radius π‘Ÿπ‘’=𝑅𝑔𝑒/𝑛, exists.

Statement 8
If the solution of the IVP π‘§π‘–ξ…žξ…ž=πœŒβ‹…π‘”β‹…π‘§π‘–βˆ’π‘π‘–π›Ύξ‚ƒξ€·π‘§1+ξ…žπ‘–ξ€Έ2ξ‚„3/2βˆ’1π‘Ÿβ‹…ξ‚ƒξ€·π‘§1+ξ…žπ‘–ξ€Έ2ξ‚„β‹…π‘§ξ…žπ‘–,𝑧𝑖𝑅𝑔𝑖=0,π‘§ξ…žπ‘–ξ€·π‘…π‘”π‘–ξ€Έ=tan𝛼𝑐(2.16) is convex, then it does not represent the inner free surface of an appropriate drop-like static meniscus.

Comment 12. The above statement shows that the inner free surface of an appropriate drop-like static meniscus can be obtained only for those values of 𝑝𝑖 for which the IVP (2.16) is not globally convex.

Statement 9
If 𝑝𝑖 satisfies the inequality 𝑝𝑖𝛾<βˆ’π‘…π‘”π‘–,(2.17) then the solution 𝑧𝑖=𝑧𝑖(π‘Ÿ) of the IVP (2.15) is convex on the interval 𝐼∩[𝑅𝑔𝑖,(𝑅𝑔𝑖+𝑅𝑔𝑒)/2], where 𝐼 is the maximal interval of the existence of 𝑧𝑖(π‘Ÿ).

Comment 13. The above statement shows that if 𝑝𝑖 is in the range defined by the inequality (2.17), then it is impossible to obtain a static meniscus, appropriate for the growth of a tube of inner radius π‘Ÿπ‘– situated in the range [𝑅𝑔𝑖,(𝑅𝑔𝑖+𝑅𝑔𝑒)/2].

Statement 10
If the solution 𝑧𝑖=𝑧𝑖(π‘Ÿ) of the NLBVP (2.2) and (2.5) represents the inner free surface of an appropriate concave static meniscus on the closed interval [𝑅𝑔𝑖,π‘šβ‹…π‘…π‘”π‘–] with 1<π‘š<(𝑅𝑔𝑖+𝑅𝑔𝑒)/2⋅𝑅𝑔𝑖, then the following inequalities hold: 1π›Όπ‘šβˆ’1𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑖⋅cosπ›Όπ‘βˆ’π›Ύπ‘…π‘”π‘–β‹…sin𝛼𝑐≀𝑝𝑖≀1(π›Όπ‘šβˆ’1)𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑖⋅sin𝛼𝑔+(π‘šβˆ’1)β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘–β‹…tanπ›Όπ‘βˆ’π›Ύπ‘šβ‹…π‘…π‘”π‘–β‹…cos𝛼𝑔.(2.18)

Comment 14. The above statement shows that in order to obtain the inner free surface of a concave static meniscus, appropriate for the growth of a tube of inner radius π‘Ÿπ‘–=π‘šβ‹…π‘…π‘”π‘–, 𝑝𝑖 has to be searched in the range defined by the inequalities (2.18).

Consequence 5
If π‘š=(𝑅𝑔𝑖+𝑅𝑔𝑒)/2⋅𝑅𝑔𝑖, then π‘Ÿπ‘–=(𝑅𝑔𝑖+𝑅𝑔𝑒)/2 (i.e., π‘Ÿπ‘– is the middle point of the interval βŒŠπ‘…π‘”π‘–,π‘…π‘”π‘’βŒ‹) and the inequalities (2.18) become: 𝛼2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…cosπ›Όπ‘βˆ’π›Ύπ‘…π‘”π‘–β‹…sin𝛼𝑐≀𝑝𝑖𝛼≀2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…sin𝛼𝑔+ξ€Ίπ‘…πœŒβ‹…π‘”β‹…π‘”π‘’βˆ’π‘…π‘”π‘–ξ€»2β‹…tanπ›Όπ‘βˆ’2⋅𝛾𝑅𝑔𝑒+𝑅𝑔𝑖⋅cos𝛼𝑔.(2.19)

Comment 15. The above consequence shows that in order to obtain the inner free surface of a concave static meniscus, appropriate for the growth of a tube of inner radius π‘Ÿπ‘–=(𝑅𝑔𝑖+𝑅𝑔𝑒)/2, 𝑝𝑖 has to be searched in the range defined by the inequalities (2.19).

Consequence 6
If 𝑝𝑖 verifies: 𝑝𝑖𝛼<2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…cosπ›Όπ‘βˆ’π›Ύπ‘…π‘”π‘–β‹…sin𝛼𝑐,(2.20) then there is no π‘Ÿπ‘– in the closed interval [𝑅𝑔𝑖,(𝑅𝑔𝑖+𝑅𝑔𝑒)/2] for which the NLBVP (2.2) and (2.5) possesses a concave solution.

Comment 16. The above consequence shows that if 𝑝𝑖 is in the range defined by the inequalities (2.19), then it is impossible to obtain a static meniscus, having concave inner free surface and appropriate for the growth of a tube of inner radius π‘Ÿπ‘– situated in the range [𝑅𝑔𝑖,((𝑅𝑔𝑖+𝑅𝑔𝑒)/2)].

Consequence 7
If π‘šβ†’1, then π‘Ÿπ‘–=π‘šβ‹…π‘…π‘”π‘–β†’π‘…π‘”π‘– and 𝑝𝑖→+∞.

Comment 17. The above consequence shows that in order to obtain the inner free surface of a concave static meniscus appropriate for the growth of a tube of inner radius π‘Ÿπ‘–β‰ˆπ‘…π‘”π‘–, 𝑝𝑖 has to be very high.

Statement 11
If π‘š and 𝑝𝑖 verify the inequalities 𝑅1<π‘š<𝑔𝑖+𝑅𝑔𝑒2⋅𝑅𝑔𝑖,𝑝𝑖>1π›Όπ‘šβˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑖⋅sin𝛼𝑔+(π‘šβˆ’1)β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘–β‹…tanπ›Όπ‘βˆ’π›Ύπ‘šβ‹…π‘…π‘”π‘–β‹…cos𝛼𝑔,(2.21) then there exist π‘Ÿπ‘– in the closed interval [𝑅𝑔𝑖,π‘šβ‹…π‘…π‘”π‘–] and a concave solution of the NLBVP (2.2) and (2.5) on the interval [𝑅𝑔𝑖,π‘Ÿπ‘–].

Comment 18. The above statement shows that if π‘š and 𝑝𝑖 verify (2.21), then a static meniscus having concave inner free surface and appropriate for the growth of a tube of inner radius π‘Ÿπ‘– in the range (𝑅𝑔𝑖,π‘šβ‹…π‘…π‘”π‘–), is obtained.

Consequence 8
If for 1<π‘šξ…ž<π‘š<(𝑅𝑔𝑖+𝑅𝑔𝑒)/2⋅𝑅𝑔𝑖 and 𝑝𝑖 the following inequalities hold: 1π›Όπ‘šβˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑖⋅sin𝛼𝑔+(π‘šβˆ’1)β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘–β‹…tanπ›Όπ‘βˆ’π›Ύπ‘šβ‹…π‘…π‘”π‘–β‹…cos𝛼𝑔<𝑝𝑖<1π‘šξ…žπ›Όβˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑖⋅cos𝛼𝑐+𝛾𝑅𝑔𝑖⋅sin𝛼𝑐,(2.22) then there exist π‘Ÿπ‘– in the interval [π‘šξ…žβ‹…π‘…π‘”π‘–,π‘šβ‹…π‘…π‘”π‘–] and a concave solution of the NLBVP (2.2) and (2.5) on the interval [𝑅𝑔𝑖,π‘Ÿπ‘–].

Comment 19. The above consequence shows that if 𝑝𝑖 is in the range defined by the inequalities (2.22), then a static meniscus having a concave inner free surface and appropriate for the growth of a tube of inner radius π‘Ÿπ‘–, situated in the range [π‘šξ…žβ‹…π‘…π‘”π‘–,π‘šβ‹…π‘…π‘”π‘–], is obtained.

Statement 12
A concave solution 𝑧𝑖=𝑧𝑖(π‘Ÿ) of the NLBVP (2.2) and (2.5) is a weak minimum of the energy functional of the melt column.

Comment 20. The above statement shows that a static meniscus having a concave inner free surface appropriate for the growth of a tube of inner radius π‘Ÿπ‘–, situated in the range [𝑅𝑔𝑖,(𝑅𝑔𝑖+𝑅𝑔𝑒)/2], is stable.

Statement 13
If 𝑝𝑖 satisfies 𝑝𝑖<(βˆ’π›Ύ/𝑅𝑔𝑖)β‹…sin𝛼𝑐 and there exists π‘Ÿπ‘–βˆˆ(𝑅𝑔𝑖,(𝑅𝑔𝑖+𝑅𝑔𝑒)/2) such that on the interval [𝑅𝑔𝑖,π‘Ÿπ‘–] an appropriate static meniscus exists, then 𝑝𝑖 verifies: βˆ’π›Ύπ‘…π‘”π‘–<𝑝𝑖𝛾<βˆ’π‘…π‘”π‘–β‹…sin𝛼𝑐.(2.23)

Comment 21. Inequality (2.23) locates the values of 𝑝𝑖 for which eventually nonglobally concave (convex-concave) static meniscus can be obtained.

Statement 14
If (βˆ’π›Ύ/𝑅𝑔𝑖)β‹…sin𝛼𝑐<𝑝𝑖 and there exists π‘Ÿπ‘–=π‘šβ‹…π‘…π‘”π‘–;π‘šβˆˆ(1,(𝑅𝑔𝑖+𝑅𝑔𝑒)/2⋅𝑅𝑔𝑖) such that on [𝑅𝑔𝑖,π‘Ÿπ‘–] an appropriate non-globally concave static meniscus exists, then 𝑝𝑖 satisfies the inequalities βˆ’π›Ύπ‘…π‘”π‘–β‹…sin𝛼𝑐<𝑝𝑖<πœŒβ‹…π‘”β‹…(π‘šβˆ’1)⋅𝑅𝑔𝑖⋅tan𝛼𝑐+𝛾𝑅𝑔𝑖.(2.24)

Comment 22. Inequality (2.24) locates the values of 𝑝𝑖 for which eventually concave-convex static menisci can be obtained.

Statement 15
If π‘š and 𝑝𝑖 satisfy the inequalities 𝑅1<π‘š<𝑔𝑖+𝑅𝑔𝑒2⋅𝑅𝑔𝑖,𝑝𝑖>𝛾𝑅𝑔𝑖+πœŒβ‹…π‘”β‹…(π‘šβˆ’1)⋅𝑅𝑔𝑖⋅tan𝛼𝑐,(2.25) then the solution 𝑧𝑖(π‘Ÿ) of the IVP (2.16) is concave on the interval 𝐼∩[𝑅𝑔𝑖,π‘šβ‹…π‘…π‘”π‘–] where 𝐼 is the maximal interval of the existence of 𝑧𝑖(π‘Ÿ).

Comment 23. Inequality (2.25) locates the values of 𝑝𝑖 for which we can obtain eventually appropriate static menisci.

3. Creation of an Appropriate Drop-Like Static Meniscus for the Growth of a Tube of Inner Radius π‘Ÿπ‘–and Outer Radius π‘Ÿπ‘’

In this section it will be shown in which kind the explicit formulas presented in the above section can be used for the creation of an appropriate drop-like meniscus when 𝛼𝑐,𝛼𝑔,𝜌,𝛾,𝑅𝑔𝑖,𝑅𝑔𝑒,π‘Ÿπ‘–andπ‘Ÿπ‘’ are given a priori. In the same time the melt column β€œheight” 𝐻, which has to be used, is found in function of the pressure of the gas flow introduced in the furnace for release the heat.

3.1. Creation of the Outer Free Surface

For the creation of the appropriate outer free surface, the following limits were considered: 𝐿𝑒1𝑛(𝑛)=π›Όπ‘›βˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑒⋅cos𝛼𝑐+𝛾𝑅𝑔𝑒⋅cos𝛼𝑔,𝐿𝑒2𝑛(𝑛)=π›Όπ‘›βˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑒⋅sin𝛼𝑔+π‘›βˆ’1π‘›β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘’β‹…tan𝛼𝑐𝛾+𝑛⋅𝑅𝑔𝑒⋅sin𝛼𝑐,𝐿𝑒3(𝑛)=π‘›βˆ’1π‘›β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘’β‹…tan𝛼𝑐𝛾+𝑛⋅𝑅𝑔𝑒⋅cos𝛼𝑐,𝐿𝑒4(𝑛)=π‘›βˆ’1π‘›β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘’β‹…tan𝛼𝑐𝛾+𝑛⋅𝑅𝑔𝑒,𝑙𝑒1𝛼=2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…cos𝛼𝑐+𝛾𝑅𝑔𝑒⋅cos𝛼𝑔,𝑙𝑒2𝛼=2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…sin𝛼𝑔+ξ€·π‘…πœŒβ‹…π‘”β‹…π‘”π‘’βˆ’π‘…π‘”π‘–ξ€Έ2β‹…tan𝛼𝑐+2β‹…π›Ύπ‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…sin𝛼𝑐,𝑙𝑒3=𝛾𝑅𝑔𝑒⋅sin𝛼𝑐(3.1) for 𝑛>1.

These limits are represented for π‘›βˆˆ(1;1.1) in Figure 2. The computations were performed in MathCAD V.13 using the following numerical data 𝑅𝑔𝑒=4.8β‹…10βˆ’3[m];𝑅𝑔𝑖=4.2β‹…10βˆ’3[m];𝛼𝑐=63.80[];𝛼=1.1135rad𝑔=28.90[]ξ€Ί=0.5044rad;𝜌=6582kg/m3ξ€»;𝛾=4.2β‹…10βˆ’1[]ξ€ΊN/m;𝑔=9.81m/s2ξ€».(3.2) When 𝑛1=1.03226, that is, π‘Ÿ1𝑒=4.65β‹…10βˆ’3[m], we have the following.

(i)If there exists a concave outer free surface, appropriate for the growth of a tube of outer radius π‘Ÿ1𝑒=4.65β‹…10βˆ’3[m], then according to Statement 3 this can be obtained for a value of 𝑝𝑒 which is in the range (𝐿𝑒1(𝑛),𝐿𝑒2(𝑛))=(134.85;164.49)[Pa].(ii)Taking into account the above fact, in order to create a concave outer free surface, appropriate for the growth of a tube of which outer radius is equal to π‘Ÿ1𝑒=4.65β‹…10βˆ’3[m], we have solved the IVP (2.8) for different values of 𝑝𝑒in the range (134.85;164.49)[Pa].

More precisely, we have integrated the following system: π‘‘π‘§π‘’π‘‘π‘Ÿ=βˆ’tan𝛼𝑒𝑑𝛼𝑒1π‘‘π‘Ÿ=βˆ’cosπ›Όπ‘’β‹…ξ‚Έπ‘”β‹…πœŒβ‹…π‘§π‘’βˆ’π‘π‘’π›Ύ+1π‘Ÿβ‹…sin𝛼𝑒(3.3) for 𝑧𝑒(𝑅𝑔𝑒)=0, π‘§ξ…žπ‘’(𝑅𝑔𝑒)=βˆ’tan𝛼𝑐 and different 𝑝𝑒. The obtained outer radii π‘Ÿπ‘’ versus 𝑝𝑒 are represented in Figure 3, which shows that the desired outer radius π‘Ÿ1𝑒=4.65β‹…10βˆ’3[m] is obtained for π‘ξ…žπ‘’=149.7[Pa].

Actually, as it can be seen on the same figure, for π‘ξ…žπ‘’=149.7[Pa], we can obtain also a second outer radius π‘Ÿ2𝑒=3.8β‹…10βˆ’3[m], which is not anymore in the desired range ((𝑅𝑔𝑖+𝑅𝑔𝑒)/2,𝑅𝑔𝑒). Moreover, the outer free surface of this meniscus is not globally concave; it is a convex-concave meniscus (Figure 4).

Taking into account π‘π‘šβ‰ˆ0 [7, 23, 24], the melt column height in this case is π»ξ…žπ‘’=βˆ’[1/πœŒβ‹…π‘”]β‹…[π‘ξ…žπ‘’+𝑝𝑒𝑔], where 𝑝𝑒𝑔β‰₯0 is the pressure of the gas flow (introduced in the furnace for release the heat from the outer side of the tube wall). When 𝑝𝑒𝑔=0 [7, 24], then π»ξ…žπ‘’ is negative, π»ξ…žπ‘’=βˆ’2.31β‹…10βˆ’3[m]; that is, the crucible melt level has to be with βˆ’π»ξ…žπ‘’=2.31β‹…10βˆ’3[m] above the shaper top level. When 𝑝𝑒𝑔=800[Pa], then π»ξ…žπ‘’=βˆ’15β‹…10βˆ’3[m]; that is, the crucible melt level has to be with βˆ’π»ξ…žπ‘’=15β‹…10βˆ’3[m] above the shaper top level.

Additional remarks are as follows.

(i)For 𝑝𝑒 in the range (βˆ’βˆž,𝑙𝑒3)=(βˆ’βˆž,78.51)[Pa], according to Statement 2, the outer free surface of the meniscus is globally convex, and, according to Statement 1, such a meniscus is not appropriate for the growth of a tube of outer radius π‘Ÿ1𝑒=4.65β‹…10βˆ’3[m]. (ii)For 𝑝𝑒 in the range (𝑙𝑒3,𝑙𝑒1)=(78.51;105.73)[Pa], according to Consequence 2, it is impossible to obtain concave outer free surface, which is appropriate for the growth of a tube of outer radius π‘Ÿ1𝑒=4.65β‹…10βˆ’3[m].(iii)For 𝑝𝑒 in the range (𝑙𝑒1,𝐿𝑒1(𝑛))=(105.73;134.85)[Pa], according to the Statement 3, it is impossible to obtain concave outer free surface, which is appropriate for the growth of a tube of outer radius π‘Ÿ1𝑒=4.65β‹…10βˆ’3[m].(iv)For 𝑝𝑒 in the range (𝐿𝑒2(𝑛);+∞)=(164.49;+∞)[Pa], according to the Statement 4, the outer free surface is concave and is appropriate for the growth of a tube which outer radius π‘Ÿπ‘’ is higher than π‘Ÿ1𝑒=4.65β‹…10βˆ’3[m].
3.2. Creation of an Appropriate Inner Free Surface

For the creation of the appropriate inner free surface, the following limits were considered 𝐿𝑖11(π‘š)=π›Όπ‘šβˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑖⋅cosπ›Όπ‘βˆ’π›Ύπ‘…π‘”π‘–β‹…sin𝛼𝑐,𝐿𝑖21(π‘š)=π›Όπ‘šβˆ’1⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2𝑅𝑔𝑖⋅sin𝛼𝑔+(π‘šβˆ’1)β‹…πœŒβ‹…π‘”β‹…π‘…π‘”π‘–β‹…tanπ›Όπ‘βˆ’1π‘šβ‹…π›Ύπ‘…π‘”π‘–β‹…cos𝛼𝑔,𝐿𝑖3(π‘š)=πœŒβ‹…π‘”β‹…(π‘šβˆ’1)⋅𝑅𝑔𝑖⋅tan𝛼𝑐+𝛾𝑅𝑔𝑖,𝑙𝑖1𝛼=2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…cosπ›Όπ‘βˆ’π›Ύπ‘…π‘”π‘–β‹…sin𝛼𝑐,𝑙𝑖2𝛼=2⋅𝛾⋅𝑐+π›Όπ‘”βˆ’πœ‹/2π‘…π‘”π‘’βˆ’π‘…π‘”π‘–β‹…sin𝛼𝑔+ξ€·π‘…πœŒβ‹…π‘”β‹…π‘”π‘’βˆ’π‘…π‘”π‘–ξ€Έ2β‹…tanπ›Όπ‘βˆ’2⋅𝛾𝑅𝑔𝑒+𝑅𝑔𝑖⋅cos𝛼𝑔,𝑙𝑖3𝛾=βˆ’π‘…π‘”π‘–β‹…sin𝛼𝑐,𝑙𝑖4𝛾=βˆ’π‘…π‘”π‘–,(3.4) for π‘š>1.

These limits are represented for π‘šβˆˆ(1;1.1) in Figure 5. The computations were performed in MathCAD V.13 using the same numerical data as for the creation of the appropriate outer free surface.

When π‘š1=1.03571, that is, π‘Ÿξ…žπ‘–=0.00435[m], we obtain the following.

(i)If there exists a concave inner free surface for the growth of a tube of inner radius π‘Ÿξ…žπ‘–=0.00435[m], then, according to the Statement 10, this can be obtained for a value of 𝑝𝑖 which is in the range (𝐿𝑖1(π‘š1),𝐿𝑖2(π‘š1))=(βˆ’31.46,βˆ’1.07)[Pa].

Taking into account the above fact, in order to create a concave inner free surface, appropriate for the growth of a tube which inner radius is equal to π‘Ÿξ…žπ‘–=0.00435[m], we have solved the IVP (2.16) for different values of 𝑝𝑖 in the range (𝐿𝑖1(π‘š1),𝐿𝑖2(π‘š1))=(βˆ’31.46,βˆ’1.07)[Pa]. More precisely, we have integrated the following system: π‘‘π‘§π‘–π‘‘π‘Ÿ=tan𝛼𝑖,𝑑𝛼𝑖=1π‘‘π‘Ÿcosπ›Όπ‘–β‹…ξ‚Έπ‘”β‹…πœŒβ‹…π‘§π‘–βˆ’π‘π‘–π›Ύβˆ’1π‘Ÿβ‹…sin𝛼𝑖(3.5) for 𝑧𝑖(𝑅𝑔𝑖)=0, π‘§ξ…žπ‘–(𝑅𝑔𝑖)=tan𝛼𝑐 and different 𝑝𝑖. The obtained inner radii π‘Ÿπ‘– versus 𝑝𝑖 are represented in Figure 6, which shows that the desired inner radius π‘Ÿ1𝑖=4.35β‹…10βˆ’3[m] is obtained for π‘ξ…žπ‘–=βˆ’16.2[Pa].

Taking π‘π‘šβ‰ˆ0 [7, 23, 24], the melt column height in this case is π»ξ…žπ‘–=βˆ’(1/πœŒβ‹…π‘”)β‹…[π‘ξ…žπ‘–+𝑝𝑖𝑔], where 𝑝𝑖𝑔β‰₯0 is the pressure of the gas flow (introduced in the furnace for release the heat from the inner side of the tube wall). When 𝑝𝑖𝑔=0 [7, 24], then π»ξ…žπ‘– is positive, π»ξ…žπ‘–=0.25β‹…10βˆ’3[m]; that is, the crucible melt level has to be with π»ξ…žπ‘–=0.25β‹…10βˆ’3[m] under the shaper top level. When 𝑝𝑖𝑔=800[Pa], then π»ξ…žπ‘–=βˆ’12.1β‹…10βˆ’3[m]; that is, the crucible melt level has to be with βˆ’π»ξ…žπ‘’=12.1β‹…10βˆ’3[m] above the shaper top level.

Additional remarks are as follows.

(i)For 𝑝𝑖 in the range (βˆ’βˆž,𝑙𝑖4)=(βˆ’βˆž,βˆ’100)[Pa], according to Statement 9, the meniscus inner free surface is convex on the interval 𝐼∩[𝑅𝑔𝑖,((𝑅𝑔𝑖+𝑅𝑔𝑒)/2)], and according to Statement 8, such a meniscus is not appropriate for the growth of a tube of inner radius π‘Ÿξ…žπ‘–=0.00435[m].(ii)For 𝑝𝑖 in the range (𝑙𝑖4,𝑙𝑖3)=(βˆ’100,89.72)[Pa] if an appropriate static meniscus exists, then the inner free surface of the meniscus is not globally concave.(iii)If 𝑝𝑖 is in the range (𝑙𝑖3,+∞)=(89.72,+∞)[Pa] and an appropriate nonglobally concave inner free surface exists, then 𝑝𝑖 is less than 𝐿𝑖3(π‘š1)=119.68[Pa].(iv)For 𝑝𝑖 in the range (𝐿𝑖2(π‘š),+∞)=(βˆ’1.07,+∞)[Pa] the obtained inner free surface is appropriate for the growth of a tube of inner radius π‘Ÿπ‘– which is less than the desired radius π‘Ÿξ…žπ‘–=0.00435[m].
3.3. Creation of a Concave Meniscus, Appropriate for the Growth of a Tube of Inner Radius π‘Ÿπ‘– and Outer Radius π‘Ÿπ‘’

For creating a concave meniscus, appropriate for the growth of a tube with outer radius π‘Ÿξ…žπ‘’=0.00465 [m] (𝑛1=𝑅𝑔𝑒/π‘Ÿπ‘’=1.03226) and inner radius π‘Ÿξ…žπ‘–=0.00435 [m] (π‘š1=𝑅𝑔𝑖/π‘Ÿπ‘–=1.03571) the melt column heights (with respect to the crucible melt level) have to be 𝐻1𝑖1=βˆ’β‹…ξ€Ίπ‘πœŒβ‹…π‘”1𝑖+𝑝𝑖𝑔1=βˆ’β‹…ξ€ΊπœŒβ‹…π‘”βˆ’16.2+𝑝𝑖𝑔,𝐻1𝑒1=βˆ’β‹…ξ€Ίπ‘πœŒβ‹…π‘”1𝑒+𝑝𝑒𝑔1=βˆ’β‹…ξ€ΊπœŒβ‹…π‘”149.7+𝑝𝑒𝑔.(3.6)

In the following we will present three different procedures to create such a meniscus.

Procedure A [𝐻1𝑒=𝐻1𝑖; 𝑝𝑖𝑔≠𝑝𝑒𝑔]
When the shaper outer top is at the same level as the shaper inner top, with respect to the crucible melt level, then the relation 𝐻1𝑒=𝐻1𝑖 holds and hence: (1/πœŒβ‹…π‘”)β‹…[βˆ’16.2+𝑝𝑖𝑔]=(1/πœŒβ‹…π‘”)β‹…[149.7+𝑝𝑒𝑔]. It follows that the pressure of the gas flow, introduced in the furnace for releasing the heat from the inner wall of the tube, 𝑝𝑖𝑔, has to be higher than the pressure of the gas flow introduced in the furnace for release the heat from the outer wall of the tube, 𝑝𝑒𝑔; π‘π‘–π‘”βˆ’π‘π‘’π‘”=149.7+16.2=165.9[Pa].The results of the integration of the system (3.3) and (3.5) for 𝑝𝑒=𝑝1𝑒=149.7[Pa] and 𝑝𝑖=𝑝1𝑖=βˆ’16.2[Pa] when the shaper outer top is at the same level as the shaper inner top, with respect to the crucible melt level, are represented in Figure 7.Figure 7 shows also that the inner free surface is higher than the outer free surface: 𝑧𝑖(π‘Ÿ1𝑖)=β„Žπ‘Žπ‘–=2.865β‹…10βˆ’4[m];𝑧𝑒(π‘Ÿ1𝑒)=β„Žπ‘Žπ‘’=2.860β‹…10βˆ’4[m];β„Žπ‘Žπ‘–βˆ’β„Žπ‘Žπ‘’=0.005β‹…10βˆ’4[m].

Procedure B [𝑝𝑒𝑔=𝑝𝑖𝑔; 𝐻1𝑒≠𝐻1𝑖]
When the pressure of the gas flow introduced in the furnace for releasing the heat from the outer wall of the tube, 𝑝𝑒𝑔, is equal to the pressure of the gas flow introduced in the furnace for release the heat from the inner wall of the tube 𝑝𝑖𝑔, then the equalities 𝑝𝑒𝑔=βˆ’π‘1π‘’βˆ’πœŒβ‹…π‘”β‹…π»1𝑒=βˆ’149.7βˆ’πœŒβ‹…π‘”β‹…π»1𝑒 and 𝑝𝑖𝑔=βˆ’π‘1π‘–βˆ’πœŒβ‹…π‘”β‹…π»1𝑖=16.2βˆ’πœŒβ‹…π‘”β‹…π»1𝑖 imply the equality 𝐻1π‘–βˆ’π»1𝑒=(1/πœŒβ‹…π‘”)β‹…165.9=2.569β‹…10βˆ’3[m].(3.7)

In other words, when 𝑝𝑒𝑔=𝑝𝑖𝑔, then the level difference between the shaper inner top and outer top has to be 𝐻1π‘–βˆ’π»1𝑒=2.569β‹…10βˆ’3[m]. In this condition the inner surface top level β„Žπ‘π‘– with respect to the shaper outer top level is β„Žπ‘π‘–=𝐻1π‘–βˆ’π»1𝑒+β„Žπ‘Žπ‘–=(1/πœŒβ‹…π‘”)β‹…(βˆ’π‘1𝑒+𝑝1𝑖)+β„Žπ‘Žπ‘–=2.855β‹…10βˆ’3[m] (see Figure 8).

Hence the level difference between the inner free surface top β„Žπ‘π‘– and outer free surface top β„Žπ‘π‘’=β„Žπ‘Žπ‘’ in this case is β„Žπ‘π‘–βˆ’β„Žπ‘π‘’=(1/πœŒβ‹…π‘”)β‹…(βˆ’π‘1𝑒+𝑝1𝑖)+(β„Žπ‘Žπ‘–βˆ’β„Žπ‘Žπ‘’)=2.5695β‹…10βˆ’3[m].

Due to the difficulties which can appear in this case in the creation of an appropriate thermal field, which assures the solidification in the right places, it is more appropriate to use 𝐻1𝑒=𝐻1𝑖 and gas flows having the property (Procedure A) 𝑝𝑖𝑔=𝑝𝑒𝑔+(𝑝1π‘–βˆ’π‘1𝑒)=𝑝𝑒𝑔+165.9[Pa]. For this choice the difference between the inner free surface height and outer free surface height with respect to the shaper top level is only β„Žπ‘–βˆ’β„Žπ‘’=0.0005β‹…10βˆ’3[m], and the creation of the thermal field, which assures the solidification in the right places, can be easier.

Procedure C [𝑝𝑒𝑔≠𝑝𝑖𝑔; 𝐻1𝑒≠𝐻1𝑖]
Finally, to create an appropriate concave meniscus for which the inner meniscus is equal to the outer meniscus height, we will proceed by changing both the pressures 𝑝𝑖𝑔, 𝑝𝑒𝑔 of the gas flows and the inner and outer top levels 𝐻𝑖,𝐻𝑒 of the shaper.
For example, if we take 𝑝𝑒𝑔=100[Pa] for creating an appropriate outer free surface, then the crucible melt level has to be with βˆ’π»1𝑒=(1/πœŒβ‹…π‘”)β‹…[149.7+100][m]=38.67β‹…10βˆ’4[m] above the shaper outer top level.
Now we take 𝐻1𝑖=𝐻1π‘’βˆ’(β„Žπ‘Žπ‘–βˆ’β„Žπ‘Žπ‘’)=βˆ’38.67β‹…10βˆ’4βˆ’0.005β‹…10βˆ’4=βˆ’38.675β‹…10βˆ’4[m]. After that, in order to satisfy 𝑝𝑖=βˆ’π‘π‘–π‘”βˆ’πœŒβ‹…π‘”β‹…π»1𝑖=βˆ’16.2[Pa], we take the pressure 𝑝𝑖𝑔 equal to 𝑝𝑖𝑔=16.2βˆ’πœŒβ‹…π‘”β‹…(βˆ’38.675β‹…10βˆ’4)=265.92[Pa]. For this choice, we will have an appropriate drop-like meniscus for which β„Žπ‘π‘–=β„Žπ‘π‘’, as it is shown in Figure 9, obtained by integration of systems (3.3) and (3.5) for the given numerical data.

4. The Effect of the Choice of the Shaper Radii

In this sequence we intend to evaluate numerically the effect of the choice of the shaper radii for the growth of a tube which outer radius is equal to π‘Ÿπ‘’, and inner radius is equal to π‘Ÿπ‘–.

For this purpose, the procedure described in sequence 3 will be applied for the growth of a tube having the same outer radius π‘Ÿξ…žπ‘’=0.00465 [m] and inner radius π‘Ÿξ…žπ‘–=0.00435 [m], but the shaper radii changed: the inner radius of the shaper now is 𝑅𝑔𝑖=4β‹…10βˆ’3[m] (instead of 4.2β‹…10βˆ’3[m]), and the outer radius of the shaper now is 𝑅𝑔𝑒=5β‹…10βˆ’3[m] (instead of 4.8β‹…10βˆ’3[m]).

Applying the procedure in this case we find that to create a static meniscus which outer free surface is appropriate for the growth, 𝑝𝑒 has to be 𝑝𝑒=124.38[Pa], and to create a static meniscus which inner free surface is appropriate for the growth, 𝑝𝑖 has to be 𝑝𝑖=βˆ’41.35[Pa]. The so-obtained outer and inner free surfaces are stable. The inner free surface height is β„Žπ‘–=6.57β‹…10βˆ’4[m], and the outer free surface height is β„Žπ‘’=6.603β‹…10βˆ’4[m].

When we take the shaper inner top at the same level as the shaper outer top, then the pressure difference between the pressure of the inner gas flow 𝑝𝑖𝑔 and the outer gas flow 𝑝𝑒𝑔 has to be π‘π‘–π‘”βˆ’π‘π‘’π‘”=π‘π‘’βˆ’π‘π‘–=165.73[Pa], and the level difference between the inner and outer free surface top is equal to β„Žπ‘–βˆ’β„Žπ‘’=βˆ’0.033β‹…10βˆ’4[m].

When we take the pressure of the inner gas flow 𝑝𝑖𝑔 equal to the pressure of the outer gas flow 𝑝𝑒𝑔, then the level difference between the shaper inner top 𝐻𝑖 and the shaper outer top 𝐻𝑒 has to be π»π‘–βˆ’π»π‘’=(1/πœŒβ‹…π‘”)β‹…[π‘π‘–βˆ’π‘π‘’]=2.567β‹…10βˆ’3[m], and the level difference between the inner and outer free surface top is equal to the difference 𝐻𝑖+β„Žπ‘–βˆ’(𝐻𝑒+β„Žπ‘’)=(1/πœŒβ‹…π‘”)β‹…[π‘π‘–βˆ’π‘π‘”]+β„Žπ‘–βˆ’β„Žπ‘’=2.5637β‹…10βˆ’3[m].

5. Conclusions

(i)Stable and drop-like static meniscus, appropriate for the growth of a single-crystal cylindrical tube with a priori specified inner and outer radii, can be created by the choice of the difference between the inner top level and outer top level of the shaper, or by the choice of the pressure of the gas flows, introduced inside and outside of the tube for release the heat from the tube walls.(ii)The values of the differences between the inner and outer top level of the shaper, or the pressure in the inner and outer gas flows, depend on the shaper radii. This dependence can be relevant.(iii)The realization of the thermal conditions, which assure that for the obtained static meniscus the solidification conditions are satisfied at the right places, is easier when we take the shaper inner top at the same level as the shaper outer top, and we use inner and outer gas flows having different pressures (𝑝𝑖𝑔≠𝑝𝑒𝑔). The setting of the thermal conditions is not discussed in this paper.

Acknowledgments

The authors thank Professor Koichi Kakimoto for the useful discussions and suggestions concerning the tube growth process. They are grateful to the Romanian National Authority for Research for supporting the research under the Grant ID 354 no. 7/2007.

References

  1. T. A. Lograsso, D. L. Schlagel, and A. O. Pecharsky, β€œSynthesis and characterization of single crystalline Gd5(SixGe1-x)4 by the Bridgman method,” Journal of Alloys and Compounds, vol. 393, no. 1-2, pp. 141–146, 2005. View at: Publisher Site | Google Scholar
  2. C. MarΓ­n and A. G. Ostrogorsky, β€œGrowth of Ga-doped Ge0.98Si0.02 by vertical Bridgman with a baffle,” Journal of Crystal Growth, vol. 211, no. 1, pp. 378–383, 2000. View at: Publisher Site | Google Scholar
  3. K. Kadokura and Y. Takano, β€œGermanium-silicon single crystal growth using an encapsulant in a silica ampoule,” Journal of Crystal Growth, vol. 171, no. 1-2, pp. 56–60, 1997. View at: Publisher Site | Google Scholar
  4. O. V. Smirnova, V. V. Kalaev, Yu. N. Makarov, N. V. Abrosimov, H. Riemann, and V. N. Kurlov, β€œThree-dimensional unsteady modeling analysis of silicon transport in melt during Cz growth of Ge1-xSix bulk crystals,” Journal of Crystal Growth, vol. 303, no. 1, pp. 141–145, 2007. View at: Publisher Site | Google Scholar
  5. B. Depuydt, A. Theuwis, and I. Romandic, β€œGermanium: from the first application of Czochralski crystal growth to large diameter dislocation-free wafers,” Materials Science in Semiconductor Processing, vol. 9, no. 4-5, pp. 437–443, 2006. View at: Publisher Site | Google Scholar
  6. D. Yang, J. Chen, H. Li et al., β€œMicro-defects in Ge doped Czochralski grown Si crystals,” Journal of Crystal Growth, vol. 292, no. 2, pp. 266–271, 2006. View at: Publisher Site | Google Scholar
  7. L. Erris, R. W. Stormont, T. Surek, and A. S. Taylor, β€œThe growth of silicon tubes by the EFG process,” Journal of Crystal Growth, vol. 50, pp. 200–211, 1980. View at: Publisher Site | Google Scholar
  8. J. C. Swartz, T. Surek, and B. Chalmerst, β€œThe EFG process applied to the growth of silicon ribbons,” Journal of Electronic Materials, vol. 4, no. 2, pp. 255–279, 1975. View at: Publisher Site | Google Scholar
  9. T. Surek, B. Chalmers, and A. I. Mlavsky, β€œThe edge-defined film-fed growth of controlled shape crystals,” Journal of Crystal Growth, vol. 42, pp. 453–465, 1977. View at: Publisher Site | Google Scholar
  10. J. P. Kalejs, A. A. Menna, R. W. Stormont, and J. W. Hutchinson, β€œStress in thin hollow silicon cylinders grown by the edge-defined film-fed growth technique,” Journal of Crystal Growth, vol. 104, no. 1, pp. 14–19, 1990. View at: Publisher Site | Google Scholar
  11. S. Rajendran, C. C. Chao, D. P. Hill, J. P. Kalejs, and V. Overbye, β€œMagnetic and thermal field model of EFG system,” Journal of Crystal Growth, vol. 109, no. 1-4, pp. 82–87, 1991. View at: Publisher Site | Google Scholar
  12. S. Rajendran, M. Larrousse, B. R. Bathey, and J. P. Kalejs, β€œSilicon carbide control in the EFG system,” Journal of Crystal Growth, vol. 128, no. 1–4, pp. 338–342, 1993. View at: Publisher Site | Google Scholar
  13. S. Rajendran, K. Holmes, and A. Menna, β€œThree-dimensional magnetic induction model of an octagonal edge-defined film-fed growth system,” Journal of Crystal Growth, vol. 137, no. 1-2, pp. 77–81, 1994. View at: Publisher Site | Google Scholar
  14. A. Roy, B. MacKintosh, J. P. Kalejs, Q.-S. Chen, H. Zhang, and V. Prasad, β€œNumerical model for inductively heated cylindrical silicon tube growth system,” Journal of Crystal Growth, vol. 211, no. 1, pp. 365–371, 2000. View at: Publisher Site | Google Scholar
  15. A. Roy, H. Zhang, V. Prasad, B. Mackintosh, M. Ouellette, and J. P. Kalejs, β€œGrowth of large diameter silicon tube by EFG technique: modeling and experiment,” Journal of Crystal Growth, vol. 230, no. 1-2, pp. 224–231, 2001. View at: Publisher Site | Google Scholar
  16. D. Sun, C. Wang, H. Zhang, B. Mackintosh, D. Yates, and J. Kalejs, β€œA multi-block method and multi-grid technique for large diameter EFG silicon tube growth,” Journal of Crystal Growth, vol. 266, no. 1–3, pp. 167–174, 2004. View at: Publisher Site | Google Scholar
  17. H. Behnken, A. Seidl, and D. Franke, β€œA 3D dynamic stress model for the growth of hollow silicon polygons,” Journal of Crystal Growth, vol. 275, no. 1-2, pp. e375–e380, 2005. View at: Publisher Site | Google Scholar
  18. B. MacKintosh, A. Seidl, M. Ouellette, B. Bathey, D. Yates, and J. Kalejs, β€œLarge silicon crystal hollow-tube growth by the edge-defined film-fed growth (EFG) method,” Journal of Crystal Growth, vol. 287, no. 2, pp. 428–432, 2006. View at: Publisher Site | Google Scholar
  19. H. Kasjanow, A. Nikanorov, B. Nacke, H. Behnken, D. Franke, and A. Seidl, β€œ3D coupled electromagnetic and thermal modelling of EFG silicon tube growth,” Journal of Crystal Growth, vol. 303, no. 1, pp. 175–179, 2007. View at: Publisher Site | Google Scholar
  20. V. A. Tatarchenko, Shaped Crystal Growth, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1993.
  21. A. V. Borodin, V. A. Borodin, V. V. Sidorov, and I. S. Pet'kov, β€œInfluence of growth process parameters on weight sensor readings in the Stepanov (EFG) technique,” Journal of Crystal Growth, vol. 198-199, part 1, pp. 215–219, 1999. View at: Publisher Site | Google Scholar
  22. A. V. Borodin, V. A. Borodin, and A. V. Zhdanov, β€œSimulation of the pressure distribution in the melt for sapphire ribbon growth by the Stepanov (EFG) technique,” Journal of Crystal Growth, vol. 198-199, part 1, pp. 220–226, 1999. View at: Publisher Site | Google Scholar
  23. S. N. Rossolenko, β€œMenisci masses and weights in Stepanov (EFG) technique: ribbon, rod, tube,” Journal of Crystal Growth, vol. 231, no. 1-2, pp. 306–315, 2001. View at: Publisher Site | Google Scholar
  24. B. Yang, L. L. Zheng, B. Mackintosh, D. Yates, and J. Kalejs, β€œMeniscus dynamics and melt solidification in the EFG silicon tube growth process,” Journal of Crystal Growth, vol. 293, no. 2, pp. 509–516, 2006. View at: Publisher Site | Google Scholar
  25. St. Balint and L. Tanasie, β€œNonlinear boundary value problems for second order differential equations describing concave equilibrium capillary surfaces,” Nonlinear Studies, vol. 15, no. 4, pp. 277–296, 2008. View at: Google Scholar | MathSciNet

Copyright © 2009 Stefan Balint and Agneta M. Balint. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

554Β Views | 479Β Downloads | 0Β Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.