The molten pool shape and thermocapillary convection during melting or welding of metals or alloys are self-consistently predicted from scale analysis. Determination of the molten pool shape and transport variables is crucial due to its close relationship with the strength and properties of the fusion zone. In this work, surface tension coefficient is considered to be negative, indicating an outward surface flow, whereas high Prandtl number represents a thinner thickness of the thermal boundary layer than that of momentum boundary layer. Since Marangoni number is usually very high, the domain of scaling is divided into the hot, intermediate and cold corner regions, boundary layers on the solid-liquid interface and ahead of the melting front. The results find that the width and depth of the pool, peak and secondary surface velocity, and maximum temperatures in the hot and cold corner regions can be explicitly and separately determined as functions of working variables or Marangoni, Prandtl, Peclet, Stefan, and beam power numbers. The scaled results agree with numerical data, different combinations among scaled equations, and available experimental data.

The molten pool shape and thermocapillary convection in melting or welding of metals or alloys having negative surface tension coefficients and Prandtl number greater than unity are determined from a scale analysis. Negative surface tension coefficient indicates that the surface flow is in outward direction, while Prandtl number greater than unity represents that boundary layer thickness of conduction is less than that of momentum. Determination of the molten pool shape is crucial due to its close relationship with the strength, microstructure and properties of the fusion zone. Since Marangoni and Reynolds number are usually greater than ten thousands, transport processes can be determined by scale analysis. In this work, the molten pool is divided into the hot, intermediate and cold corner regions on the flat free surface, boundary layers on the solid-liquid interface and ahead of the melting front for analysis. The results find that the pool shape, surface speed and temperature profiles can be self-consistently evaluated as functions of Marangoni, Prandtl, Peclet, Stefan, and beam power numbers. The predictions agree with numerical computations and experimental data in the literature.

The lack of significant buoyancy effects in zero-gravity conditions poses an issue with fluid transfer in a stagnant liquid. In this thesis, the movement of a bubble or droplet in both stagnant and rotating liquids is analysed and presented numerically using computational fluid dynamics (CFD). The governing continuum conservation equations for two-phase flow are solved using the commercial software package (2011). The Volume of Fluid (VOF) method is used to track the liquid/gas interface in 2D and 3D domains. User-Defined Functions (UDFs) are employed in order to include the effect of surface tension gradient and fluid properties as a function of temperature, with a view to efficiently investigating temperature effects on the properties of the two phases. The flow is driven via Marangoni influence induced by the surface tension gradient, which in turn drives the bubble/droplet from the cold to the hot region. For stationary liquid, the results indicate that the scaled velocity of the bubble decreases with an increase in the Marangoni number, which agrees with the results of previous space experiments. An expression for predicting the scaled velocity of a bubble has been regressed based on the obtained data from the present numerical study for thermal Marangoni numbers up to 10,721. An expression for predicting the scaled velocity of a Fluorinert droplet migrating in oil has also been presented for an MaT range from 24.05 to 2771. The interactions of two droplets in thermocapillary motion have also been studied and compared with the results obtained for the isolated droplet. The results have shown that the leading droplet will not move faster than if it were isolated, as the trailing droplet has no influence on the velocity of the leading droplet. Three-dimensional results show that no bubbles broke in any of the cases observed and agglomeration could occur during thermocapillary migration for bubbles placed side by side. The results of the motion of a singular and multiple bubbles incorporating thermocapillary forces in a rotating liquid in a zero-gravity environment have been presented for the first time. When the Rossby number is 1, the effects of rotation are important. Furthermore, the deflection of the gas bubble motion increases towards the axis of rotation with a decrease in the Rossby number (Ro). Bubble population balance modelling has been investigated in normal gravity using Luo kernels for breakage and agglomeration and two different laminar kernels for zero-gravity conditions. The simulations covered a wide range of scenarios and results are presented as a bell and histogram shapes for number density and particle percentage distribution, respectively.

No description available.

This study focuses ontransient heat flow behavior in which centralizing energy on themetal makes metal surface come to aheat molten state with centralized heat source . This flow field is two-dimensional transient model, using Phase-field method and Two-phase flow to simulatemetal surface. In this study is under considerations of the mass conservation equation, momentum equation, energy equation and the level-set equation, regardless of the impact due to the concentration diffusion. At last it will show the flow of the molten zone caused by temperature, and the flows in molten zone forced by thermocapillary which is caused byvariation of temperature.

The research presented in this thesis was motivated by the desire to understand the flow field within a new digital microfluidic device currently under development. This required an investigation of the dynamics of a droplet migrating along the surface of another fluid due to interfacial surface tension gradients. The quantitative analysis of the flow field presented in this thesis provides the first known solution for the velocity field in a migrating droplet confined to an interface. The first step towards gaining insight into the flow field was accomplished by using the method of reflections to obtain an analytical model for a submerged droplet migrating near a free surface. The submerged droplet model enabled the analysis of the velocity field and droplet migration speed and their dependence on the fluid properties. In general, the migration velocity of a submerged droplet was found to differ dramatically from the classic problem of thermocapillary migration in an unbounded substrate. A boundary-collocation scheme was developed to determine the flow field and migration velocity of a droplet floating trapped at the air-substrate interface. The numerical method was found to produce accurate solutions for the velocity and temperature fields for nearly all parameters. This numerical scheme was used to judge the accuracy of the flow field obtained by the submerged droplet model. In particular, the model was tested using parameter values taken from a digital microfluidic device. It was determined that the submerged droplet model captured most of the flow structure within the microfluidic droplet. However, for a slightly different choice of parameters, agreement between the two methods was lost. In this case, the numerical scheme was used to uncover novel flow structures.

Ce travail de thèse est consacré à l’étude des écoulements thermogravitaires et thermocapillaires induits par chauffage laser dans des couches liquides. Le chauffage d’un système à deux liquides donne naissance à deux effets thermiques. Le premier est dû à la variation de la masse volumique avec la température et le second à la variation de la tension interfaciale avec la température. Nous avons dans un premier temps étudié ces deux contributions séparément. En confrontant expériences, théorie et simulations numériques nous avons démontré que la morphologie des écoulements thermogravitaires dépend de l’épaisseur de la couche liquide ainsi que de la largeur du champ de température. Ensuite nous nous sommes intéressés à l’étude théorique et numérique de l’effet thermocapillaire. Cette étude a révélé qu’il est possible d’étudier les écoulements à partir des déformations d’interface induites par ces mêmes écoulements. Nous avons montré qu’il existe deux régimes de déformations en fonction du rapport de hauteurs et du signe de la variation de la tension interfaciale avec la température. Enfin, nous nous sommes intéressés à la compétition entre ces deux mécanismes. L’analyse des déformations d’interface et la comparaison avec un modèle à une dimension a montré qu’en fonction du rapport de hauteurs des couches liquides nous avons une transition d’un régime d’écoulements thermocapillaires vers une régime d’écoulements thermogravitaires. De plus, nous avons montré expérimentalement et numériquement qu’il est possible d’induire une instabilité d’origine thermogravitaire conduisant à la formation d’un pont liquide. / This thesis work is dedicated to thermocapillary and thermogravitary flows induced bylaser-heating in liquid layers. The laser-heating of a two fluids system induces two thermaleffects. The first effect come from the variation of the density with temperature andthe second one is due to the variation of the interfacial tension with the temperature.In a first part, we study separately these two mechanism. With a comparison betweenexperiments, theory and numerical simulations we demonstrated that the morphology ofthe flows depends on the thickness of the layer and the width of the temperature field.Then, we studied numerically and theoretically the thermocapillary effect. That revealedit is possible to understand the flows from the interface deformations. We highligth theseare two deformation regimes depending the thickness ratio and the sign of the interfacialtension gradient. Finally, we studied the competition between the two mechanism and theexperiments revealed a transition between two flow regimes. The comparison of these resultsand a theoritical model showed there is a transition from a thermocapillary regime to athermogravitary regime. In addition, we showed the possibility to induce an instibility bythermogravitary effect which can lead to the formation of a liquid bridge.

Cette étude se concentre sur la compréhension de la physique des instabilités dans différents écoulements de cisaillement, particulièrement la cavité entraînée et la cavité thermocapillaire, où l'écoulement d'un fluide incompressible est assuré soit par le mouvement d’une ou plusieurs parois, soit par des contraintes d’origine thermique.Un code spectral a été validé sur le cas très étudié de la cavité entrainée par une paroi mobile. Il est démontré dans ce cas que l'écoulement transit d'un régime stationnaire à un instationnaire au-delà d'une valeur critique du nombre de Reynolds. Ce travail est le premier à donner une interprétation physique de l'évolution non monotonique du nombre de Reynolds critique en fonction du facteur d'aspect. Lorsque le fluide est entraîné par deux parois mobiles, la cavité entraînée possède un plan de symétrie particulièrement sensible. Des solutions asymétriques peuvent être observés en plus de la solution symétrique au-dessus d'une certaine valeur du nombre de Reynolds. La transition oscillatoire entre la solution symétrique et les solutions asymétriques est expliquée physiquement par les forces en compétition. Dans le cas asymétrique, l'évolution de la topologie permet à l'écoulement de rester stationnaire avec l'augmentation du nombre de Reynolds. Lorsque l'équilibre est perdu une instabilité se manifeste par l'apparition d'un régime oscillatoire dans l'écoulement asymétrique.Dans une cavité thermocapillaire rectangulaire avec une surface libre, Smith et Davis prévoient deux types d'instabilités convectives thermiques: des rouleaux longitudinaux stationnaires et des ondes hydrothermales instationnaires. L'apparition de ses instabilités a été mis en évidence à plusieurs reprises expérimentalement et numériquement. Alors que les applications impliquent souvent plus d'une surface libre, il semble qu'il y ait peu de connaissances sur l'écoulement thermocapillaire entraînée avec deux surfaces libres. Un film liquide libre soumis à des contraintes thermocapillaires possède un plan de symétrie particulier comme dans le cas de la cavité entrainée par deux parois mobiles. Une étude de stabilité linéaire avec deux profils de vitesse pour le film liquide libre est présentée avec différents nombres de Prandtl. Au-delà d'un nombre de Marangoni critique, il est découvert que ces états de base sont sensibles à quatre types d'instabilités convectives thermiques qui peuvent conserver ou briser la symétrie du système. Les mécanismes qui permettent de prédire ces instabilités sont également découverts et interpréter en fonction de la valeur du nombre de Prandtl du fluide. La comparaison avec les travaux de Smith et Davis est faite. Une simulation numérique directe permet de valider les résultats obtenus avec l'étude de stabilité de linéaire. / This study focuses on the understanding of the physics of different instabilities in driven cavities, specifically the lid-driven cavity and the thermocapillarity driven cavity where flow in an incompressible fluid is driven either due to one or many moving walls or due to surface stresses that appear from surface tension gradients caused by thermal gradients. A spectral code is benchmarked on the well-studied case of the lid-cavity driven by one moving wall. In this case, It is shown that the flow transit form a steady regime to unsteady regime beyond a critical value of the Reynolds number. This work is the first to give a physical interpretation of the non-monotonic evolution of the critical Reynolds number versus the size of the cavity. When the fluid is driven by two facing walls moving in the same direction, the cavity possesses a plane of symmetry particularly sensitive. Thus, asymmetrical solutions can be observed in addition to the symmetrical solution above a certain value of the Reynolds number. The oscillatory transition between the symmetric solution and asymmetric solutions is explained physically by the forces in competition. In the asymmetric case, the change of the topology allows the flow to remain steady with increasing the Reynolds number. When the equilibrium is lost, an instability manifests by the appearance of an oscillatory regime in the asymmetric flow. In a rectangular cavity thermocapillary with a free surface, Smith and Davis found two types of thermal convective instabilities: steady longitudinal rolls and unsteady hydrothermal waves. The appearance of its instability has been highlighted repeatedly experimentally and numerically. While applications often involve more than a free surface, it seems that there is little knowledge about the thermocapillary driven flow with two free surfaces. A free liquid film possesses a particular plane of symmetry as in the case of the two-sided lid-driven cavity. A linear stability analysis for the free liquid film with two velocity profiles is presented with various Prandtl numbers. Beyond a critical Marangoni number, it is observed that these basic states are sensitive to four types of thermal convective instabilities, which can keep or break the symmetry of the system. Mechanisms that predict these instabilities are discovered and interpreted according to the value of the Prandtl number of the fluid. Comparison with the work of Smith and Davis is made. A direct numerical simulation is done to validate the results obtained with the linear stability analysis.

A parallel finite element code with adaptive meshing was developed and used to study three dimensional, time-dependent fluid flows caused by thermocapillary convection as well as temperature and dopant distribution in fusion welding and floating zone crystal growth. A comprehensive numerical model of the three dimensional time-dependent fluid flows in a weld pool had been developed. This model considered most of the physical mechanisms involved in gas tungsten arc welding. The model helped obtaining the actual chaotic time-dependent melt flow. It was found that the fluid flow in the weld pool was highly complex and influenced the weld pool’s depth and width. The physicochemical model had also been studied and applied numerically in order to simulate the surfactant adsorption onto the surface effect to the surface tension of the metal liquid in a weld pool. Another model, a three dimensional time-dependent, with adaptive mesh refinement and coarsening was applied for simulating the effect of weak flow on the radial segregation in floating zone crystal growth. The phase change equation was also included in this model in order to simulate the real interface shape of floating zone. In the new parallel code, a scheme that keeps the level of node and face instead of the complete history of refinements was utilized to facilitate derefinement. The information was now local and the exchange of information between each and every processor during the derefinement process was minimized. This scheme helped to improve the efficiency of the parallel adaptive solver. / QC 20100527

Recently the computational speed and the functions of the numerical methods are advancing rapidly. It is the future trend that using the computational fluid dynamics (CFD) to perform simulation for making up the experimental deficiency, reducing the risk, improving the quality of the product, and saving the cost of research and development. A two-dimensional simulation was employed to study the melt/air and melt/solid interface shapes of the miniature molten zone formed in the laser-heated pedestal growth (LHPG) system. Using non-orthogonal body-fitting grid system with control-volume finite difference method, the interface shape can be determined both efficiently and accurately. During stable growth, the dependence of the molten-zone length and shape on the heating CO2 laser is examined in detail under both the maximum and the minimum allowed powers with various growth speeds. The effect of gravity for the miniature molten zone is also simulated, which reveals the possibility for a horizontally oriented LHPG system. Such a horizontal system is good for the growth of long crystal fibers. After comparing with the shape of the molten zone in terms of the experiment and the analysis of the simulation shown as above. Heat transfer and fluid flow in the LHPG system are analyzed near the deformed interfaces. The global thermal distributions of the crystal fiber, the melt, and the source rod are described by temperature and its axial gradient within length of ~10 mm. As compared with the growth of bulk crystal of several centimeters in dimension, natural convection drops six orders in magnitude due to smaller melt volume; therefore, conduction rather than convection determines the temperature distribution in the molten zone. Moreover, thermocapillary convection rather than mass-transfer convection becomes dominant. The symmetry and mass flow rate of double eddy pattern are significantly influenced by the molten-zone shape due to the diameter reduction and the large surface-tension-temperature coefficient in the order of 10-4~10-3. According to the analysis shown as above, the results could be further extended for the analysis of the concentration profile and study of horizontal growth.