Search results

A rigorous Finite Element (FE) formulation based on an enthalpy technique is developed for solving coupled nonlinear heat conduction/mass diffusion problems with phase change. The FE formulation consists of a fully coupled heat conduction and solute diffusion formulation, with solid‐liquid phase change, where the effects of pressure and convection are neglected. A full enthalpy method is employed eliminating singularities which result from abrupt changes in heat capacity at the phase interfaces. The FE formulation is based on the fixed grid technique where the elements are two dimensional, four noded quadrilaterals with the primary variables being enthalpy and average solute concentration. Temperature and solid mass fraction are calculated on a local level at each integration point of an element. A fully consistent Newton‐Raphson method is used to solve the global coupled equations and an Euler backward difference scheme is used for the temporal discretization. The solution of the enthalpy‐temperature relationship is carried out at the integration points using a Newton‐Raphson method. A secant method employing the regula falsi technique takes into account sudden jumps or sharp changes in the enthalpy‐temperature behaviour which occur at the phase zone interfaces. The Euler backward difference integration rule is used to calculate the solid mass fraction and its derivatives. A practical example is analysed and results are presented.

The aim of the present work is to study the effect of processing conditions on solidification path and resultant microstructure and further predict the solidification behavior of gas-atomized Sn-5mass%Pb droplets.

Combined with previous models for in-flight droplet nucleation and non-equilibrium solidification, a simulation method is applied to four typical containerless solidification conditions with helium, nitrogen or argon gas at two different gas jet velocities, in the presence of 10 or 500 ppm oxygen. The simulation outputs distribution of primary dendrite composition, tip velocity and tip radius with radial distance from the nucleation point, and the fraction solid at the end of recalescence and the post-recalescence duration. Both surface and internal nucleation are considered. The possible dendritic fragmentation in the post-recalescence stage is also discussed.

Result indicates that dendritic fragmentation is not likely to occur in droplets solidifying along the paths considered in the simulation.

The simulation method applies to any droplet-based solidification process for which droplet cooling schedule is known and thus provides a scientific basis for powder quality assurance.

Proper selection of the stability parameter determines the accuracy of dendrite tip kinetics at a single crystal scale. Recently developed sophisticated phase field modelling of a single grain evolution provides evidence that this parameter is not constant during the process. Nevertheless, in the commonly used micro-macroscopic simulations of alloy solidification, it is a common practice to use a constant value of the stability parameter, resulting from the marginal stability theory. This paper aims to address the issue of how this inaccuracy in modelling crystal growth kinetics can influence numerically predicted zones of columnar and equiaxed dendrites and the macro-segregation formation.

Using the original authors’ micro-macroscopic computer simulation model of binary alloy solidification, the calculations have been performed for the Kurz-Giovanola-Trivedi (KGT) crystal growth kinetics with two different values of the stability parameter, and for two different compositions of Al-Cu alloys. The computational model is based on single domain-based formulation of transport equations, which are discretized on control-volume mesh. To identify zones of different grain structures, developing within the two-phase liquid-solid region, an envelope of columnar dendrite tips is tracked on a fixed non-orthogonal, triangular control volume grid. The models of porous and slurry media are used, along with the concept of the switching function, to account for diverse flow resistances in the columnar and equiaxed crystal zones. The numerical predictions are carefully studied to address the question of how the chosen stability parameter influences macroscopic structures of a cast, the most important issue from the engineering point of view.

The carried-out comprehensive numerical analysis shows that the value of the stability parameter of the KGT-constrained dendrite growth model does not have a direct significant impact on the macrosegregation formation. It, however, visibly influences the undercooling along the front, separating different dendritic structures and the size of the undercooled melt region where the equiaxed grains can develop. It also affects the amount of eutectic phase created.

To the best of the authors’ knowledge, this is the first attempt at estimating the influence of some inaccuracies, caused by possible ambiguities in choosing the stability constant of the KGT law, on numerically predicted macroscopic fields of solute concentration, the developing zones of columnar and equiaxed crystals and the macrosegregation patterns.

The general aim of this work is to calculate the extent of the equiaxed zone in continuously cast steel products. Free equiaxed grains can grow only in undercooled liquid regions. Undercooling of the bulk liquid occurs because the columnar dendrite tips growing from the mould reject solutes in the liquid. The specific aim of this contribution is to calculate the thermal and physical state of continuously cast steel long products assuming a columnar solidification mode, taking into account the tip undercooling at the solidification front. A 2‐D heat transfer model has been developed where the columnar solidification mode is assumed. The calculation of the undercooling at the advancing solidification front is coupled with the heat transfer equation. The comparison between the results of the present model and the classical heat transfer model indicates the importance of modelling the undercooling phenomenon. The influence of the secondary cooling has also been studied.

This study aims to simulate the dendritic growth in Stokes flow by iteratively coupling a domain and boundary type meshless method.

A preconditioned phase-field model for dendritic solidification of a pure supercooled melt is solved by the strong-form space-time adaptive approach based on dynamic quadtree domain decomposition. The domain-type space discretisation relies on monomial augmented polyharmonic splines interpolation. The forward Euler scheme is used for time evolution. The boundary-type meshless method solves the Stokes flow around the dendrite based on the collocation of the moving and fixed flow boundaries with the regularised Stokes flow fundamental solution. Both approaches are iteratively coupled at the moving solid–liquid interface. The solution procedure ensures computationally efficient and accurate calculations. The novel approach is numerically implemented for a 2D case.

The solution procedure reflects the advantages of both meshless methods. Domain one is not sensitive to the dendrite orientation and boundary one reduces the dimensionality of the flow field solution. The procedure results agree well with the reference results obtained by the classical numerical methods. Directions for selecting the appropriate free parameters which yield the highest accuracy and computational efficiency are presented.

A combination of boundary- and domain-type meshless methods is used to simulate dendritic solidification with the influence of fluid flow efficiently.

The purpose of this paper is to endorse the idea of using a special post-calculating front tracking (FT) procedure, along with the enthalpy-porosity front tracking (EP-FT) single continuum model, in order to identify zones of different dendritic microstructures developing in the mushy zone during cooling and solidification of a binary alloy.

The 2D and 3D algorithms of the FT approach along with different crystal growth laws were implemented in macroscopic calculations of binary alloy solidification with the identification of different dendrite zones developing during the process.

Direct comparison of results predicted by the FT model with that based on the concept of the critical value of the solid volume fraction shows the sensitivity of the latter on an arbitrary assumed value of the dendrite coherency point (DCP). Moreover, for a carefully chosen DCP value the second model provides results that are close to those given by the FT-based approach. It is also observed that the macro-segregation pattern obtained by the proposed method is hardly influenced by chosen dendrite tip kinetics.

To the best authors’ knowledge, for the first time the 3D FT model has been used along with the enthalpy porosity approach to simulate the development of zones of different dendrite morphology during binary alloy solidification. And, a weak influence of assumed different dendrite tip kinetics on the macro-segregation pattern has been proved, what justifies this underlying assumption of the EP-FT method.

In the present work, a numerical method, based on the well established enthalpy technique, is developed to simulate the growth of binary alloy equiaxed dendrites in presence of melt convection. The paper aims to discuss these issues.

The principle of volume-averaging is used to formulate the governing equations (mass, momentum, energy and species conservation) which are solved using a coupled explicit-implicit method. The velocity and pressure fields are obtained using a fully implicit finite volume approach whereas the energy and species conservation equations are solved explicitly to obtain the enthalpy and solute concentration fields. As a model problem, simulation of the growth of a single crystal in a two-dimensional cavity filled with an undercooled melt is performed.

Comparison of the simulation results with available solutions obtained using level set method and the phase field method shows good agreement. The effects of melt flow on dendrite growth rate and solute distribution along the solid-liquid interface are studied. A faster growth rate of the upstream dendrite arm in case of binary alloys is observed, which can be attributed to the enhanced heat transfer due to convection as well as lower solute pile-up at the solid-liquid interface. Subsequently, the influence of thermal and solutal Peclet number and undercooling on the dendrite tip velocity is investigated.

As the present enthalpy based microscopic solidification model with melt convection is based on a framework similar to popularly used enthalpy models at the macroscopic scale, it lays the foundation to develop effective multiscale solidification.

The motivation for this work is to establish a model that not only includes the main factors resulting in macrosegregation but also retains simplicity and consistency for the sake of potential application in casting practice.

A mathematical model for the numerical simulation of thermosolutal convection and macrosegregation in the solidification of multicomponent alloys is developed, in which the coupled macroscopic mass, momentum, energy and species conservation equations are solved. The conservation equations are discretized by using the control volume‐based finite difference method, in which an up‐wind scheme is adopted to deal with the convection term. The alternative direction implicit procedure and a line‐by‐line solver, based on the tri‐diagonal matrix algorithm, are employed to iteratively solve the algebraic equations. The velocity‐pressure coupling is handled by using the SIMPLE algorithm.

Based on the present study, the liquid flow near the dendritic front is believed to play an important role in large‐scale transport of the solute species. The numerical or experimental results in the literatures on the formation of channel segregation, especially those about the location of the initial flow as well as the morphology of the liquidus front, are well supported by the present investigation.

The modelling is limited to dealing with the thermosolutal convection of two‐dimensional cases. More complicated phenomena (e.g. crystal movement) and 3D geometry should be considered in future research.

The present model can be used to analyze the effects of process parameters on macrosegregation and, with further development, could be applied as a useful tool in casting practice.

The numerical simulation demonstrates the capability of the model to simulate the thermosolutal convection and macrosegregation in alloy solidification. It also shows that the present model has good application potential in the prediction and control of channel segregation.

This paper aims to tackle the problem of some ambiguity of the momentum equation formulation in the commonly used macroscopic models of two‐phase solid/liquid region, developing during alloy solidification. These different appearances of the momentum equation are compared and the issue is addressed of how the choice of the particular form affects velocity and temperature fields.

Attention is focused on the ensemble averaging method, which, owing to its stochastic nature, is a new promising tool for setting up the macroscopic transport equations in highly inhomogeneous multiphase micro‐ and macro‐structures, with morphology continuously changing in time when the solidification proceeds. The basic assumptions of the two other continuum models, i.e. based on the classical mixture theory and on the volume‐averaging technique, are also unveiled. These three different forms of the momentum equation are then compared analytically and their impact on calculated velocity and temperature distribution in the mushy zone is studied for the selected test problem of binary alloy solidification driven by diffusion and thermal natural convection in a square mould.

It is found that a chosen appearance of the momentum equation mildly affects temporal velocity/temperature, and shapes of the phase interface at longer times of the solidification.

This mainly results from small variations of the liquid fraction across the mushy zone and from a low solidification rate, and it may change drastically when anisotropic properties of the mushy zone, solutal convection, different phase densities and cooling conditions are considered. Therefore, further comprehensive study is needed.

The paper addresses how the different focus of the momentum equation for liquid flow is compared.

Additive manufacturing (AM) is revolutionizing the manufacturing industry due to several advantages and capabilities, including use of rapid prototyping, fabrication of complex geometries, reduction of product development cycles and minimization of material waste. As metal AM becomes increasingly popular for aerospace and defense original equipment manufacturers (OEMs), a major barrier that remains is rapid qualification of components. Several potential defects (such as porosity, residual stress and microstructural inhomogeneity) occur during layer-by-layer processing. Current methods to qualify AM parts heavily rely on experimental testing, which is economically inefficient and technically insufficient to comprehensively evaluate components. Approaches for high fidelity qualification of AM parts are necessary.

This review summarizes the existing powder-based fusion computational models and their feasibility in AM processes through discrete aspects, including process and microstructure modeling.

Current progresses and challenges in high fidelity modeling of AM processes are presented.

Potential opportunities are discussed toward high-level assurance of AM component quality through a comprehensive computational tool.