Assessing the scalability of multiphysics tools for modeling solidification and melting processes on parallel clusters

A comprehensive solution of solidification/melting processes requires the simultaneous representation of free surface fluid flow, heat transfer, phase change, nonlinear solid mechanics and, possibly, electromagnetics together with their interactions, in what is now known as multiphysics simulation. Such simulations are computationally intensive and the implementation of solution strategies for multiphysics calculations must embed their effective parallelization. For some years, together with our collaborators, we have been involved in the development of numerical software tools for multiphysics modeling on parallel cluster systems. This research has involved a combination of algorithmic procedures, parallel strategies and tools, plus the design of a computational modeling software environment and its deployment in a range of real world applications. One output from this research is the three-dimensional parallel multiphysics code, PHYSICA. In this paper we report on an assessment of its parallel scalability on a range of increasingly complex models drawn from actual industrial problems, on three contemporary parallel cluster systems.

Computational modeling of mold filling and related free-surface flows in shape casting: an overview of the challenges involved

Accurate representation of the coupled effects between turbulent fluid flow with a free surface, heat transfer, solidification, and mold deformation has been shown to be necessary for the realistic prediction of several defects in castings and also for determining the final crystalline structure. A core component of the computational modeling of casting processes involves mold filling, which is the most computationally intensive aspect of casting simulation at the continuum level. Considering the complex geometries involved in shape casting, the evolution of the free surface, gas entrapment, and the entrainment of oxide layers into the casting make this a very challenging task in every respect. Despite well over 30 years of effort in developing algorithms, this is by no means a closed subject. In this article, we will review the full range of computational methods used, from unstructured finite-element (FE) and finite-volume (FV) methods through fully structured and block-structured approaches utilizing the cut-cell family of techniques to capture the geometric complexity inherent in shape casting. This discussion will include the challenges of generating rapid solutions on high-performance parallel cluster technology and how mold filling links in with the full spectrum of physics involved in shape casting. Finally, some indications as to novel techniques emerging now that can address genuinely arbitrarily complex geometries are briefly outlined and their advantages and disadvantages are discussed.

Multiscale modelling of dendrite growth during vacuum arc remelting

Vacuum Arc Remelting (VAR) is the accepted method for producing homogeneous, fine microstructures that are free of inclusions required for rotating grade applications. However, as ingot sizes are increasing INCONEL 718 becomes increasingly susceptible to defects such as freckles, tree rings, and white spots increases for large diameter billets. Therefore, predictive models of these defects are required to allow optimization of process parameters. In this paper, a multiscale and multi-physics model is presented to predict the development of microstructures in the VAR ingot during solidification. At the microscale, a combined stochastic nucleation approach and finite difference solution of the solute diffusion is applied in the semi-solid zone of the VAR ingot. The micromodel is coupled with a solution of the macroscale heat transfer, fluid flow and electromagnetism in the VAR process through the temperature, pressure and fluid flow fields. The main objective of this study is to achieve a better understanding of the formation of the defects in VAR by quantifying the influence of VAR processing parameters on grain nucleation and dendrite growth. In particular, the effect of different ingot growth velocities on the microstructure formation was investigated. It was found that reducing the velocity produces significantly more coarse grains.

Multiscale modeling of the onset of freckle formation during vacuum arc remelting

A multiscale model for the Vacuum Arc Remelting process (VAR) was developed to simulate dendritic microstructures during solidification and investigate the onset of freckle formation. On the macroscale, a 3D multi-physics model of VAR was used to study complex physical phenomena, including liquid metal flow with turbulence, heat transfer, and magnetohydrodynamics. The results showed that unsteady fluid flow in the liquid pool caused significant thermal perturbation at the solidification front. These results were coupled into a micromodel to simulate dendritic growth controlled by solute diffusion, including local remelting. The changes in Rayleigh number as the microstructure remelts was quantified to provide an indicator of when fluid flow channels (i.e. freckles) will initiate in the mushy zone. By examining the simulated microstructures, it was found that the Rayleigh number increased more than 300 times during remelting, which suggests that thermal perturbation could be responsible for the onset of freckle formation.

A unified 3D model of the VAR process

A 3D time-dependent model of the VAR process has been developed using CFD techniques. The model solves the coupled field equations for fluid flow, heat transfer (including phase change) and electromagnetic field, for both the electrode and the ingot. The motion of the electic arc 'preferred spot' can be specified based on observations. Correlations are sought between the local gap height, resulting from instantaneous liquid pool surface shape and electrode tip shape, and the arc motion. The detailed behaviour of the melting film on the electrode tip is studies using a spectral free surface technique, which allows investigation of the drops' detachment and drip shorts.

Effect of varying electromagnetic field on the VAR process

The three-dimensional, time-dependent electromagnetic field arising from the precession of the arc centre in a vacuum arc remelting furnace is shown (in a numerical simulation) to affect the fluid flow and heat transfer conditions near the solidification front in the upper part of the ingot.

Fundamentals of steady state, non-Newtonian rimming flow

Rimming flow on the inner surface of a horizontal rotating cylinder is investigated. Using a scale analysis, a theoretical description is obtained for steady-state non-Newtonian flow. Simple lubrication theory is applied since the Reynolds number is small and the liquid film is thin. Since the Deborah number is very small the flow is viscometric. The shear-thinning number, which characterizes the shear-thinning effect, may be small or large. A general constitutive law for this kind of flow requires only a single function relating shear stress and shear rate that corresponds to a generalized Newtonian liquid. For this case the run-off condition for rimming flow is derived. Provided the run-off condition is satisfied, the existence of a continuous steady-state solution is proved. The rheological models, which show Newtonian behavior at low shear rates with transition to power-law shear thinning at moderate shear rates, are considered. Numerical results are carried out for the Carreau and Ellis models, which exhibit Newtonian behavior near the free surface and power-law behavior near the wall of the rotating cylinder.

Gas-liquid two phase flow through a vertical 90deg elbow bend

Pressure drop data are reported for two phase air-water flow through a vertical to horizontal 90° elbow bend set in 0.026 m i.d. pipe. The pressure drop in the vertical inlet tangent showed some significant differences to that found for straight vertical pipe. This was caused by the elbow bend partially choking the inflow resulting in a build-up of pressure and liquid in the vertical inlet riser and differences in the structure of the flow regimes when compared to the straight vertical pipe. The horizontal outlet tangent by contrast gave data in general agreement with literature even to exhibiting a drag reduction region at low liquid rates and gas velocities between 1 and 2 m s -1. The elbow bend pressure drop was best correlated in terms of le/d determined using the actual pressure loss in the inlet vertical riser. The data showed a general increase with fluid rates that tapered off at high fluid rates and exhibited a negative pressure region at low rates. The latter was attributed to the flow being smoothly accommodated by the bend when it passed from slug flow in the riser to smooth stratified flow in the outlet tangent. A general correlation was presented for the elbow bend pressure drop in terms of total Reynolds numbers. A modified Lockhart-Martinelli model gave prediction of the data.

Theoretical analysis on microwave heating of oil–water emulsions supported on ceramic, metallic or composite plates

A detailed theoretical analysis has been carried out to study efficient heating due to microwaves for one-dimensional (1D) oil–water emulsion samples placed on various ceramic, metallic (reflective) and ceramic–metallic composite supports. Two typical emulsion systems are considered such as oil-in-water (o/w) and water-in-oil (w/o). A preliminary study has been carried out via average power vs emulsion thickness diagram to estimate microwave power absorption within emulsion samples for various cases. The maxima in average power, also termed as ‘resonances’, are observed for specific emulsion thicknesses and the two consecutive resonances of significant magnitudes are termed as R1 and R2 modes. For both o/w and w/o emulsions, it is observed that microwave power absorption is enhanced in presence of metallic and composite supports during both R1 and R2 modes. The efficient heating strategies characterized by ‘large heating rates’ with ‘minimal thermal runaway’ i.e. uniform temperature distributions within the sample have been assessed for each type of emulsion. Based on the detailed spatial distributions of power and temperature for various cases, SiC-metallic composite support may be recommended as an optimal heating strategy for o/w samples with higher oil fractions (0.45) whereas metallic and Alumina-metallic composite supports may be favored for samples with smaller oil fractions (=0.3) during R1 mode. For w/o samples, SiC-metallic composite support may be suitable heating strategy for all ranges of water fractions during R1 mode. During R2 mode, metallic and Alumina-metallic composite supports are favored for both o/w and w/o emulsion samples. Current study recommends the efficient way to use microwaves in a single mode waveguide and the heating strategy can be suitably extended for heating of any other emulsions for which dielectric properties are easily measurable or available in the literature.

Header design for flow equalization in microstructured reactors

To enhance the uniformity of fluid flow distribution in microreactors, a header configuration consisting of a cone diffuser connected to a thick-walled screen has been proposed. The thick-walled screen consists of two sections: the upstream section constitutes a set of elongated parallel upstream channels and the downstream section constitutes a set of elongated parallel downstream channels positioned at an angle of 90 with respect to the upstream channels. In this approach the problem of flow, equalization reduces to that of flow equalization in the first and second downstream channels of the thick-walled screen. In turn, this requires flow equalization in the corresponding cross sections of the upstream channels. The computational fluid dynamics analysis of the fluid flow maldistribution shows that eight parallel upstream channels with a width of 300-600 pm are required per 1 cm of length for flow equalization. The length to width ratio of these channels has to be > 15. The numerical results suggest that the proposed header-configuration can effectively improve the performance of the downstream microstructured devices, decreasing the ratio of the maximum flow velocity to the mean flow velocity from 2 to 1.005 for a wide range of Reynolds numbers (0.5-10). 2006 American Institute of Chemical Engineers AlChE J, 53: 28-38, 2007.

Inertia-Gravity Waves Emitted from Balanced Flow: Observations, Properties, and Consequences

This paper describes laboratory observations of inertia–gravity waves emitted from balanced fluid flow. In a rotating two-layer annulus experiment, the wavelength of the inertia–gravity waves is very close to the deformation radius. Their amplitude varies linearly with Rossby number in the range 0.05–0.14, at constant Burger number (or rotational Froude number). This linear scaling challenges the notion, suggested by several dynamical theories, that inertia–gravity waves generated by balanced motion will be exponentially small. It is estimated that the balanced flow leaks roughly 1% of its energy each rotation period into the inertia–gravity waves at the peak of their generation. The findings of this study imply an inevitable emission of inertia–gravity waves at Rossby numbers similar to those of the large-scale atmospheric and oceanic flow. Extrapolation of the results suggests that inertia–gravity waves might make a significant contribution to the energy budgets of the atmosphere and ocean. In particular, emission of inertia–gravity waves from mesoscale eddies may be an important source of energy for deep interior mixing in the ocean.