Discovered in 1960 by Duwartz et al.  at Caltech, Metallic Glasses  may be defined as disordered atomic-scale structural arrangement of atoms formed as a result of rapid cooling of complex alloy systems directly from their melt state to below room temperature with a large undercooling and a suppressed kinetics in such a way that the supercooled state is retained/frozen    . This results in the formation of “glassy structure”. The process is very much similar to inorganic/oxide glass formation in which large oxide molecules (silicates/borides/aluminates/and sulphides) form a regular network retained in its frozen/supercooled liquid state . The only difference is that metallic glasses are comprised of metallic atoms rather than inorganic metallic compounds. Their atomic arrangement is based on mismatch of atomic size and quantity (minimally three) , is based on short    to medium range order    or long range disorder  (unlike metals—well defined long range order) and can be explained by other advanced theories or mechanisms such as liquid-liquid transition   , phase separation , confusion , order in disorder     , frustration , and vitrification  . Important features characterizing them are their amorphous microstructure and unique mechanical properties. Owing to absence of mechanisms promoting development and movement of dislocations, no plasticity is exhibited by bulk metallic glass matrix composites. This results in very high yield strength and elastic strain limit as there is no easy plane of flow. From fundamental definition standpoint, metallic glasses are typically different from bulk metallic glasses, in that the former is glassy structure less than 1 mm thick, while later is greater than 1 mm  . Till date, the largest bulk metallic glass produced in “as cast” condition is 80 mm diameter and 85 mm in length . There are reports of making large thin castings as casing of smart phones but they are typically less than 1 mm . Furthermore, they are characterised by special properties such as glass forming ability (GFA), and metastability.
2. Formation and Stability—Three Laws
The formation and stability of bulk metallic glasses is described by their ability to retain glassy state at room temperature. In initial days of metallic glass and bulk metallic glass development, it was very difficult to retain glassy structure at room temperature. Only alloys falling into window of very narrow compositional contrast cooled at extremely high cooling rate can form glassy structure     . A deviation from any of these parameters severely hampers their glassy state and crystallisation occurs   . This property is known as glass forming ability (GFA) . This is single most important property in metallic glass family of alloys which governs their formation and evolution. GFA has been increasingly studied and considerable progress has been made in its improvement     by alterations and improvement in both composition and processing condition’s window   . Now, alloys having multicomponent composition can be cast in glassy state even at slow cooling rates owing to their superior glass forming ability   -  which in turn is governed by various theories    -  and analytical models  .
Fundamentally, research over the period of time has yielded three basic laws which are now considered universal for forming any bulk metallic glass system . These are described as follows. Any glass forming system consist of elements which must:
1) Be three in number (at minimum). (Elements greater than 3 are considered beneficial).
2) Differ in their atomic size by 12% among three elements. (Atoms of elements with large size are considered to exhibit superior glass forming ability).
3) Have negative heat of mixing among three elements. (This ensures tendency to demix or confuse  ensuring retention of glassy structure at room temperature).
This result in new structure with high degree of densely packed atomic configurations, which in turn, results in completely new atomic configuration at local level with long range homogeneity and attractive interaction. In general bulk metallic glasses or bulk glassy alloys are typically designed around alloy systems that exhibit: 1) a deep eutectic, which decreases the amount of undercooling needed to vitrify the liquid, and 2) alloys that exhibit a large atomic size mismatch, which creates lattice stresses that frustrate crystallisation . An important way to arrive at optimum glass forming composition and then selecting alloying elements is based on proper choice of eutectic or off-eutectic composition, diameter and heat of mixing . These laws were first proposed by Abdullah S. Jabri and Douglas C. Hoffmann at Jet Propulsion Laboratory, Caltech  and Prof. Akisha Inoue at WPI-IMR, Tohoku University, Japan  independently.
3. Ductile Bulk Metallic Glass and Bulk Metallic Glass Matrix Composites
Inherent brittleness, which is a big disadvantage in manufacturing and application of bulk metallic glasses at large scale can be overpowered by inducing plasticity in glassy structure while retaining its high strength at the same time    . This can be done by various mechanisms including exploitation of intrinsic ability of glass to exhibit plasticity at very small (nano) length scale  , introduction of external impulse (obstacles) to shear band formation and propagation (ex-situ composites)  , self or externally assisted multiplication of shear bands  , formation of ductile phases in brittle glassy matrix during solidification (in-situ composites)     and transformation inside a ductile crystalline phase e.g. B2 - B19’ transformation in Zr based systems (stress/transformation induced plasticity (TRIP)    . It was known thermodynamically and observed after simulation  and experimentation     since early days that structurally constrained glass either relaxes (losses free volume) or underwent crystal structure change (devitrification) during heating. The drive for devitrification   did not came out of ingenious effort but was an outcome of natural impulse of stressed structure to attain a state of lower stress by undergoing structural change and evolution of new structures (phases) when subjected to temperature effect similar to heat treatment for crystalline metallic alloys. This result in new class of bulk metallic glass called ductile bulk metallic glass  - . Other approach which encompasses nucleation and growth of primary ductile phase in brittle glassy matrix during solidification gave rise to second type of ductile material known as ductile bulk metallic glass matrix composites. These were reported first time in 2001 by Prof. William L. Johnson’s group at Caltech  when they successfully incorporated primary ductile phase reinforcements in glassy matrix in the form of precipitates in-situ nucleating during solidification thus giving birth to “so called” family of in-situ dendrite/metallic glass matrix composites  -  Later reinforcements were observed to exhibit different morphologies other than dendrites (e.g. spheroid, whiskers and plates) but fundamental mechanism of reinforcements remained same (i.e. metal matrix composites). These materials are formed as a result of conventional solute partitioning mechanisms as observed in other metallurgical alloys resulting in copious formation of ductile phase (Ti-Zr-Nb β in case of Ti based composites , B2 CuZr in case of Zr based composites  -  or transformed B2 (B19’ martensite) in case of shape memory alloys (a special class of bulk metallic glass matrix composites)       ). This ductile phase can appear in the form of three dimensional dendrites or spheroids emerging directly from liquid during solidification. Structural changes in these alloys can be explained by “phase separation”  -  before liquid-solid or solid-solid transformation or formation and retention of “quenched in” nuclei during solidification. This is another very important route for the manufacture of these alloys. Very recently, another important phenomenon known as liquid-liquid transition (LLT)  has been observed in these alloys. This has become observable due to ability to characterize these alloys under Synchrotron radiation     and in microgravity conditions   employing container less levitation techniques .
These impart their special properties (enhanced plasticity and compressive strength) not otherwise attainable by conventional processing routes or in simple binary and ternary compositions. This however, is seldom the case and is not readily observed as compared to solid state phase separation  which is dominant mechanism in these alloys. More advanced mechanisms of forming these materials are by local microstructural evolution by phase separation right at shear bands . It narrates that solid-solid phase separation (spinodal decomposition) occurs at the onset of shear band which is the cause of microstructural evolution. Few notable class of alloys in these types of ductile composites are Ti based BMGMCs  - , Ti based shape memory BMGMC , Zr-Cu-Al-Ti  , Zr-Cu-Al-Ni , and Zr-Cu-Al-Co shape memory BMG composites . Each has its own mechanisms of formation and individual phases are formed by liquid-solid or solid-solid phase transformations.
4. Additive Manufacturing of Bulk Metallic Glasses and their Composites
Although, considerable progress has been made in advancing “as cast” sizes in bulk metallic glass and their composites, still, maximum possible diameter and length which has been produced by conventional means till date , is far from limit of satisfaction to be used in any structural engineering application. This primarily is associated with mechanical cooling rate achievable as a result of quenching effect from water cooled walls of Cu container which in itself in not enough to overcome critical cooling rate (Rc) of alloy (~0.067 K/s ) to produce a uniform large bulk glassy structure. In addition to this, occurrence of this bulk glassy structure is limited to compositions with excellent inherent glass forming ability (GFA)  . This is not observed in compositions which are strong candidates to be exploited for making large scale industrial structural components           with higher critical cooling rates (Rc) (10 K/s ). This poses a limitation to this conventional technique and urges the need of advanced manufacturing method which does not encompass these shortcomings.
Additive Manufacturing has emerged as potential technique   to fulfil this gap and produce bulk metallic glass matrix composites   in single step across a range of compositions virtually covering all spectrums    . It achieves this by exploiting very high cooling rates available in very short period transient liquid melt pool    in a small region where laser strikes the sample (LSM/LSF (solid), SLM/LENSÒ (powder)). This, when coupled with superior glass forming ability (GFA) of bulk metallic glass matrix composites (BMGMC), efficiently overcome dimensional limitation as virtually any part carrying glassy structure can be fabricated. In addition, incipient pool formation  and its rapid cooling results in extremely versatile and beneficial properties in final fabricated part such as high strength, hardness, toughness, controlled microstructure, dimensional accuracy, consolidation and integrity. The mechanism underlying this is layer-by-layer (LBL) formation, which ensures glass formation in each layer during solidification before proceeding to next layer. That’s how; a large monolithic glassy structure can be produced. This layer-by-layer (LBL) formation also helps in development of secondary phases in a multicomponent alloy    as layer preceding fusion layer (which is solidified) undergoes another heating cycle (heat treatment) below melting temperature (Tm) somewhat in the nose region of TTT diagram  which not only assist in phase transformation   but also helps in increase of toughness, homogenisation and compaction of part. This is a new, promising and growing technique of rapidly forming metal , plastic , ceramic or composite  parts by fabricating a near-net shape out of raw materials either by powder method or wire method (classified on the basis of additives used). The movement of energy source (laser or electron beam) is dictated by a CAD geometry which is fed to a computer at the back end and manoeuvred by CNC   system. Process has wide range of applicability across various industrial sectors ranging from welding  - , repair  , and cladding  -  to full scale part development.
However, there is dearth of knowledge about exact mechanisms of formation (nucleation and growth and/or liquid-liquid transition   ) of ductile phase dendrites in-situ during solidification of BMGMC happening inside liquid melt pool of additive manufacturing which is essential to further advance improvement in the process and suggest its optimisation. Modelling and simulation techniques especially those employing finite element methods (FEM) (phase field  - ), and CAFE  - ) at part scale are very helpful in explaining the evolution of microstructure and grain size development in metals and alloys. They have been extensively used in predicting solidification behaviour of various types of melts during conventional production methods    . However, their use in additive manufacturing applications     specially related to BMGMC is still in its infancy. Virtually no effort has been made to understand nucleation and growth of ductile phase crystalline dendrites in-situ during solidification in BMGMC by modelling and simulation. A step forward is taken in present study to address these gaps and bring together the strengths of different techniques and methodologies at one platform. A comprehensive review of existing literature is presented and a tale of manufacturing of these alloys is described with step-by-step arrival at set of equations which are proposition of solution methodology. Ductile bulk metallic glass metal matrix composites are proposed to be formed by taking advantage of
1) Materials chemistry: A multicomponent alloy. Its glass forming ability (GFA) is used as a measure to manipulate composition and vice-versa.
2) Solidification processing: Liquid melt pool formation, its size, shape and geometry, role of number density, size and distribution of ductile phase in resultant glassy alloy matrix. It is taken as a function of type, size and amount of nucleates (inoculant).
3) Additive manufacturing: Use of very high cooling rate inherently available in the process to (a) not only form glassy matrix but use liquid melt pool formed at very high temperature to trigger nucleation (liquid-solid transformation) of ductile phase in the form of dendrites from within the pool “in-situ” (This is done by controlling machine parameters in such a way that optimised cooling rate satisfying narrow window of “quenching” bulk metallic glasses is achieved) (b) take advantage of heating (heat treatment) of preceding layer to trigger solid-solid transformation (devitrification) again to form ductile phase and achieve homogeneity, consolidation and part integrity eliminating the need of post processing or after treatment and
4) Modelling and Simulation: Strong and powerful mathematical modelling techniques based on
a) Transient heat transfer for “liquid melt pool formation as a result of laser-matter interaction” and
b) Its “evolution-solidification” by
i) Deterministic (modified classical nucleation theory, JMAK Correction and Rappaz Modification) or
ii) Stochastic/probabilistic (3D CAFE model for nucleation and growth (solute diffusion and capillary action driven))
Modelling of microstructure evolution comprising of ductile phase (in some case equiaxed dendrites) evolution in glassy melt will be used to simulate the conditions in liquid melt pool of BMGMC during additive manufacturing. It is envisaged that effect of number density, size and distribution of ductile phase dendrites will be evaluated/verified by simulation of melt pool developed using different values of aforementioned parameters.
Note: Additive Manufacturing (AM) methods can also be classified on the basis of energy source used (i.e. Laser based or Electron Beam based).
5. A Brief Introduction to Modelling and Simulation
Although, in use since ancient Roman times , modelling and simulation picked up interest and achieved pinnacle in modern day scientific and engineering sectors with the advent of computer technology which came not more than two decades ago. Now, it has proved itself to be important integral part of product and part design, product development as well as prediction, utilisation and enhancement of properties. Various branches of modelling and simulation, ranging from part scale modelling which involves development of codes of theorems in advanced computing platform such as FORTRANÒ, JavaÒ, C++, MATLABÒ and PythonÒ to their stand alone or parallel simulations in customised simulation languages/packages such as ABAQUSÒ, SolidworksÒ, AnsysÒ and CatiaÒ to performing complex atomistic simulations in dedicated proprietary software, have now become integral part of design procedure in major industrial clusters. Its use in research and development is also becoming an important part of whole process to eliminate so called “trial and error” methods which are not only time consuming but are energy, materials and resources extensive. In materials science and engineering mainly two of its branches are routinely used. These are “part scale modelling and simulation” and “atomistic modelling and simulation”. The former is used for the complete design of complex machinery segments, equipment, assemblies, sub-assemblies, their materials of fabrication and property prediction in different regions as a function of extrinsic parameters such as heat, velocity, pressure, time whilst the latter is used for the prediction, estimation and improvement in atomic-scale properties using theories of atomic configuration and arrangement mainly relying on intrinsic parameters (such as specific heat/latent heat, heat capacity, density, heat of fusion and atomic fit or misfit). The unique ability of atomistic modelling and simulation is that it uses atomic functions and their variables to generate knowledge about their behaviour under various impulses. In both cases, the use of these methods are big help and support in saving time, materials, resources as well as improve functionality and in-service property development and behaviour prediction.
The exponential rise in the use of modelling and simulation with the advent and progress of computer technology and increase of computing power of machines gave rise to greater flexibility in the design and development process. Many difficult, or in some cases, impossible to envisage problems can now be simulated using these computing platforms. These include simulation of water flow and its patterns in rivers and channels, simulation of interior of sun, stars and other heavenly bodies, cosmic events and nuclear engineering problems. However, despite of these advantages, there are still situations and applications which limits the use of modelling and simulation techniques. These include, unavailability of strong efficient computing algorithms (with lesser approximations) needed for the replication of actual real-world situations, unavailability of real world experimental data (physical constants and thermo-physical properties) needed to simulate a particular problem, unavailability of more accurate deterministic or accurate (non-probabilistic) models using actual situations rather than basing their outcome on statistics. Owing to these reasons, there is still need for further investigation and removal of bottlenecks from modelling and simulation techniques and it is envisaged that their popularity is still at arm’s length. However, progress in these is notable and healthy.
6. Current State of art in Modelling and Simulation of Additive Manufacturing of Bulk Metallic Glass and their Composites
Nucleation and growth phenomena in single component (pure metals), binary and multicomponent alloys is rather well understood. Classical nucleation theory (CNT)  provides many answers to the behaviour of these melts. Traditionally, bulk metallic glass and their composites were produced using conventional methods (Cu mould suction casting    and twin roll casting ) in which their metastable phase (glass) and any in-situ ductile precipitates (stable phase) are nucleated based on their ability to surpass activation energy barrier. In addition, these processes, impart very high cooling rate to castings which is essential for retention of supercooled liquid (glass) at room temperature. However, they have their limitations which pose limits to their applicability and use for making complex and intricate parts and components at large scale.
Very recently, with the advent and popularity of additive manufacturing (AM), interest has sparked to exploit the inherent and fundamental advantages present in this unique process to produce BMG and BMGMCs. Additive manufacturing techniques are useful in achieving this objective as very high cooling rate in fused liquid melt pool is already present inherently to assist the formation of glassy structure which is suppression of “kinetics” and prolonging of undercooling (“thermodynamics”)—two main phenomena responsible for any phase transformation. However, the in-situ nucleation of primary phase equiaxed dendrites or spheroids during solidification and then microstructural evolution (solute diffusion and capillary action assisted) is not satisfactorily explained by classical nucleation theory alone.
Either some modification is needed in classical nucleation theory (CNT) or more reliable probabilistic microstructure evolution models (e.g. J-M-A-K Correction  or Rappaz modification) are needed to explain nucleation and growth (and other phenomena e-g liquid-liquid transition (LLT)    and phase separations ) in BMGMCs. In this work, an effort has been made to meet both requirements. Following are propositions and state of art;
1) At present scenario, there is no single hybrid/combined model which explain phenomena of heat transfer (liquid melt pool formation as a result of laser-matter interaction and its evolution-solidification) and coupled this with nucleation and growth (NG) (solute diffusion  and capillary action driven) at microscale to predict microstructure and grain size in BMGMC as melt cools in additive manufacturing melt pool.
a) Only one study has been conducted to model the same phenomena (solidification only) during Cu mold suction casting which will serve as base  in addition to very recent attempts  in which emphasis is laid on development of generalised theory rather than solving a problem.
b) Only one study has been reported on microstructure formation during twin roll casting using cellular automaton-finite element (CAFÉ)  but that is not aimed at BMGMCs, is carried out using commercial software package and does not involve any mathematical modelling at the back end. Software embedded (nucleation and growth and heat transfer) models are used only.
c) Four prominent studies namely by Zhou et al. , Zhang et al. , Zinoviev et al.  and a group at Shenyang, China  , have been reported very recently using CAFE but these are based on modelling microstructure evolution in modified additive manufacturing (HDMR , LAMP  on 316L SS , 2D CAFE , Cladding  ) processes.
d) Few studies in the past have been conducted employing selective laser melting (SLM) using CAFE    , CAPF, CAFVM , modified CAFE  approaches but none have been carried out on BMGMC.
2) No effort has been made to correlate the effect of edge-to-edge matching (E2EM)   with direction of easy heat flow and crystallographic orientation (an effect that can be inherently used in additive manufacturing).
3) No substantial study has been reported about evolution of microstructure in three dimensions in BMGMCs in additive manufacturing.
4) No effort has been made to combine the effect of changing properties with decrease of temperature (transient conditions). Most of models till now predict solutions in terms of steady state phenomena.
5) Very few studies have been carried out to combine cellular automaton (CA) with finite element (FE) in case of additive manufacturing (AM) while it is routine approach to predict grain size in case of other processes (casting, welding).
In light of this research gap, present review is compiled. An effort has been put together to overcome these shortcomings and propose a methodology for the modelling and simulation of solidification phenomena during additive manufacturing of BMGMCs. A model system Zr47.5Cu45.5Al5Co2 has been proposed owing to evolution of distinct ductile phase (CuZr-B2) in it during solidification and its tendency to show shape memory effect (exhibiting two types of martensitic phase evolution from B2 phase). Further, the method is proposed to be applied to conventional wedge shape casting geometry along with its final application to melt pool in additive manufacturing making use of powers of deterministic, probabilistic and their coupled modelling approaches. This route is proposed to get maximum benefit from application of modelling and simulation to understand nucleation and growth phenomena during solidification both in conventional as well as modern processing technologies (AM). It is envisaged that application of hybrid CAFE model by programming on FORTRANÒ/MATLABÒ and parallel simulation on ABAQUSÒ will help understand solidification in BMGMC in much better way not done elsewhere previously.
This analysis is divided into two sections. The first section deals with the evolution of the melt pool as a result of the interaction of highly localised, focused laser light with matter (metal powder). This results in the formation of a melt pool whose shape, size, geometry and transient behaviour is very much a function of the heat transfer coefficients (HTC) evolving at every step of its formation (melting and homogenisation) and dissipation (solidification). Solidification in this section is considered by a modified general (classical) nucleation theory (CNT) . Once formed, this pool travels as the laser traverses its path all along the powder bed dictated by CAD geometry at the back end. The second section deals with the microstructural evolution during solidification which is primarily a solute diffusion and capillary action dominated phenomena. This is dealt with by microscopic 2D and 3D probabilistic cellular automaton models which model nucleation and equiaxed dendritic growth resulting in the formation of the microstructure within the liquid melt pool as it solidifies. Advantages and disadvantages of cellular automaton method are also described. (Note: only “Vitrification (glass formation)” effects are considered and devitrification (glass to crystal precipitation) is not considered). The evolution of microstructure is checked against the variation of number density (nucleation density), volume fraction, size and distribution of ductile phase in the glassy matrix. Inoculants for ductile phase formation were selected previously by edge to edge matching (E2EM)  .
8. Modelling and Simulation of Heat Transfer in Liquid Melt Pool-Solidification
As the microstructure formed during selective laser melting (SLM) is mostly columnar , it is a good indicator that heat flux transfer from melt is highly unidirectional thus heat transfer from bottom is a transient 1D process. Although, heat is lost from the material in x-y plane i.e. perpendicular to the z-direction (perpendicular to build direction), its contribution is so low that it can be safely ignored. However, this was an old concept. New experimental observations have proposed a new concept according to which during SLM, a melt pool is formed, where the shape of this pool is a function of;
1) Laser power (laser beam intensity).
2) Presence of thermocapillary convection (marangoni convection).
In even more advanced and recent models,   the transfer of heat after its generation is considered by three main parameters:
1) Heat transfer due to convection.
2) Evaporation (i.e. formation of plasma) (this results in re-radiation (inverse radiation)).
3) Conduction from the bottom and the side walls.
This is very recent and advanced approach which, however, ignores marangoni convection effects. Overall, the heat transfer phenomena associated with the solidification of metal in a liquid melt pool in AM is associated with three processes:
1) Generation of heat (laser matter interaction).
2) Assimilation of heat (melting and stages of solidification).
3) Extraction of heat.
8.1. Generation of Heat (Laser-Matter Interaction)
This is the first stage of additive manufacturing in which heat is generated. The problem in this stage is related with impingement of light of certain intensity (I) on a solid surface for a certain amount of time which may results in production of heat. This interaction can be explained in terms of law known as the “Beer lambert law”.
Beer Lambert Law for Additive Manufacturing
Consider a thin layer of powder with thickness d1, on a flat disk substrate of refractory metal with thickness d2 and radius r uniformly illuminated by light of intensity I.
For absorptivity of powder (or melt) assuming uniform temperature throughout the disk, the temperature evolution is
where A(T) = Absorptivity; Q(T) = Thermal loss (convective and radiative); I = Intensity; ρ1 = Density of powder; ρ2 = Density of substrate; c1 = Specific heat of powder; c2 = Specific heat of substrate; d1 = Thickness of powder; d2 = Thickness of substrate.
Heat generated by this process is used for melt pool generation (its morphology, homogenisation, and holding (generation of supercooled liquid (SCL) region and its progression)).
8.2. Assimilation of Heat (Melting and Stages of Solidification)
As the heat generated above interacts with metal powder, it causes its melting and generation of liquid melt pool. The behaviour of a certain metal/alloy in the melt pool can be explained by its cooling curve which is briefly described below.
8.2.1. General form of Cooling Curve
A cooling curve of a metal/alloy is a plot of the variation of temperature with time. It has different regions which embodied various types of information. Cooling curves can have different shapes depending on the metal or alloy type. A schematic cooling curve is shown in Figure 1 for a single component pure metal (without any inoculants).
Its distinct regions are explained as follows;
Region above A1: This is the region in which metal is in its complete liquid-state and can be described by only melting and liquid-state homogenisation. Heat carried by metal in this region is “super heat” only and lost in the form of specific heat (mcpΔT). This homogenisation in turn depends on type of melting (gas/solid (coal)/liquid (oil) fired crucible furnace melting, electric (resistance/induction/arc) melting) and subsequent melt treatment. (Note: Homogenisation is required by some external means in case of all modes of melting. Only induction furnace is manifested by self-homogenisation due to phenomena of induction currents).
Region A1 - A: This is a region which is characterised by the loss of super heat until the first arrest point A. (Point at which the first nucleas form—explained in detail in later sections). This is also called the start of solidification. In pure metals it is a sharp point (melting point) while in alloys, it can be a range
Figure 1. Cooling curve for a single component pure metal (without any inoculants).
(melting range). In BMGMCs/multicomponent alloys, it is also called start of the super cooled region (SCL). This region is followed by undercooling (ΔTn) region which is described below
Region A - D: This is the most important region of cooling curve (present case) for pure metals. In this region, metal cools down to a specific temperature characterised by a certain minimum amount of energy (activation energy for nucleation) needed to overcome a barrier of energy (energy barrier to nucleation) to create a liquid-solid (L-S) interface eventually leading to formation of a stable nuclei out of the melt. This region is further divided into two regions A - C and C - D.
Region A - C: This is region in which undercooling occurs, heat is extracted, the temperature drops and shape of cooling curve goes down. This is characterised by two energies described in the above paragraph.
Region C - D: This is the region in which heat energy is absorbed, temperature is gained and shape of curve goes up. This is called recalescence.
1) Recalescence is gain in temperature as a result of thermal fluctuations caused by phase transformations occurring within solidifying melt/alloy. In present case, phase change is solid formation within undercooled liquid while thermal fluctuations are described by release of heat in the form of heat of fusion.
2) Region A - C is characterised by another point. Point “B” occurring in the middle of cooling curve. This is specifically shown in Figure 1 as intermediate point of Supercooled liquid region (SCL). For the present case model (transient heat transfer conditions will be modelled at this point as well to get better understanding of phenomena occurring in SCL in BMGMCs).
Region D - E: This is the region at which (after arrest point D), metal losses all its heat of fusion (mHf). In this region transformation occurs at constant temperature in such a way that all liquid gets transformed into complete solid (all fine equiaxed grains formation at mould wall (Cu mould casting)/at surface of inoculant (heterogeneous nucleation, not present case), “equiaxed-columnar” transition, growth of columnar dendrites, CET and growth of all equiaxed dendrites accomplishes). This is also called the solidification time.
Region E - F: This is the region in which solid cools. That is, after all liquid gets transformed into solid, the solid casting cools down to room temperature. This again occurs after a sharp invariant point (point F) in case of pure metals and after a range in case of multicomponent alloys.
8.2.2. Cooling Curve for Well Inoculated Zr-Based in-situ Dendrite BMGMCs
Shape of cooling curve changes its form as melt is changed from single component to binary to multicomponent alloys. This can be explained in the form of various cases.
Case 1: Well inoculated single component melt:
In these types of alloys, undercooling/undercooled region (ΔTn) diminishes and is almost absent. Inoculation with potent nuclei serves as active nucleation sites and triggers heterogeneous nucleation as the alloy reaches its first invariant point. Thus, no undercooling happens, and solid alloy directly starts cooling as all liquid gets transformed to solid at constant temperature.
Case II: Binary alloys without inoculants (slowly cooled)
In these types of alloys cooling occurs in following steps
1) Distinct undercooling occurs (characterised by drop and gain (recalescence) of temperature)
2) It is followed by region of constant temperature cooling which is called local solidification. This is only visible in case of very fluid alloys in which mushy region is very fluid/less viscous (not BMGMCs). This region is absent in most multicomponent (industrial) alloys as their solidification is dominated by mushy zone. (Note: BMGMCs are special case of alloys in which mushy region is extensively dominated but another phenomenon known as “sluggishness” governs the solidification. In these alloys, three laws  which describe BMGMC formation and evolution make sure that not only sluggishness dominates kinetics but it also ensures “glass formation” (i.e. retaining supercooled liquid at room temperature).
3) Alloy solidification range (it depends on alloy. In slowly cooled binary alloys (most laboratory conditions), this is very clearly marked (usually bears an intermediate shape)).
4) At the end of this range, alloy becomes stable momentarily at constant temperature (usually negligible in most industrial castings) at which nuclei (dendrite arm branches) grow and fills interdendritic arm spacing and other small liquid pockets. This is marked by end of solidification (in some cases, it is also characterised by start of CET and then growth of equiaxed grains).
5) Following this point, solid alloy cools to room temperature or below room temperature (in case of cryogenic cooling).
Note: For theoretical analysis, cooling curve can be of any type of combination between type of alloy (single component, binary and multicomponent), method of cooling (slow or fast) and inoculation (zero inoculation and well inoculated). All these can be drawn following rules of thermal transitions and kinetics. For simplicity and sufficiency, we will jump to cooling curve of multicomponent alloy (BMGMCs) fast cooled and well inoculated (present case).
Case III: Multicomponent alloys with inoculants (fast cooled) (present case BMGMCs)
In these types of alloys, cooling can occur following below steps (Figure 2).
1) No undercooling occurs (as there is sufficient amount (number) of potent nuclei which serve as sites for active nucleation triggering heterogeneous nucleation prior to loss of temperature (drop of cooling curve), and gain of temperature (recalescence—rise of cooling curve).
Figure 2. Cooling curve for a multicomponent alloy with inoculants (fast cooled) (BMGMC).
2) This is followed by region of constant temperature at which all liquid get transformed into solid. However, in these alloys, this region is very small (because of presence of marked mushy zone).
3) Instantly, after this region, alloy enters in “solidification range”. As the alloy is very fast cooled, this region is again not very clearly identified which is typical behaviour in case of fast cooled castings.
4) Following this, again alloy momentarily enters in brief constant temperature zone which marks starts of CET and growth of equiaxed grains (B2 CuZr phase equiaxed dendrites) until all liquid gets transformed into solid (end of solidification). This again is not very distinct as other phenomena (suppressing kinetics) dominate.
5) Finally, after this, BMGMC solidifies to room temperature.
Note: Shape of cooling curve in case of slowly cooled and fast cooled alloys is the slope of curve towards the end of cooling which is very steep in case of very fast cooled alloys (liquid melt pools (present case)).
8.3. Extraction of Heat-Determination of Heat Transfer Coefficients (HTCs)
In the development of model, heat transfer coefficients will be determined at every point of cooling curve following earlier defined one dimensional (1D) schemes . These will ensure, time of solidification calculation during cooling following above cooling curve and helps in determining shape of melt pool and its transient behaviour during cooling.
Final Time of Solidification
Final time of solidification is sum of time in each region/section of cooling curve of a particular alloy/melt. It will be determined using standard transport equations and will be used empirically to assess the conformability of additive manufacturing process. Time of solidification gives other parameters as well such as fraction of mass solidified after a time, t, which is direct measure of microstructure evolved during that time. It can be qualitatively (extrapolation) used to predict further (type (equiaxed, columnar, mix, CET) and amount) evolution of microstructure with time.
8.4. Modelling and Simulation of Nucleation (Heterogeneous) in Liquid Melt Pool-Microstructural Development
Modelling and simulation of microstructural development in liquid melt pool can be described by macroscopic and microscopic models of heat and mass transfer depending on type of alloy, its nature, number of elements, cooling curve, undercoolings (constitutional (solute/particulate), thermal, curvature, interfacial), thermal and kinetic limitations, behaviour of mushy zone, presence or absence of inoculants. These can be broadly divided into macroscopic and microscopic models  which are explained as follows;
8.4.1. Macroscopic Models
By following the regimes of macroscopic models, finite element method (FEM) and finite difference methods (FDM) can be used to explain microstructural development both during steady and transient state transport processes.
Both FEM and FDM based models cannot fully describe mushy region, its behaviour and evolution during solidification as they do not account for microscopic
1) solute diffusion and
2) capillary effects
which are primarily responsible for scale at which microstructure forms (which is very small as compared to macroscopic methods based on average continuity equations     in which it is assumed that solidification starts at liquidus and finishes at solidus/eutectic temperatures (A case of BMGMCs having good match of GFA and eutectic temperature  ). In order to overcome these limitations, microscopic models were proposed.
8.5. Microscopic Models of Microstructure Evolution/Formation during Solidification
Stage 1 Model: These models take into account the mechanism of 1) grain nucleation and 2) grain growth in alloys which are solidifying with equiaxed dendrite or eutectic microstructures . These do not account for alloys which are solidifying with columnar dendritic and planar interfaces. A modification of these accounts for equiaxed-columnar (at mould wall) and columnar to equiaxed transition (CET) in bulk of liquid. These can be used to “describe microstructures” and “prediction of grain size” in case of eutectic compositions of BMGMC. Majority of these is based on “analytical/deterministic approaches” which can be described as follows;
1) Choose a time “t” (initially non-zero value)
2) At this time t, density of grains (which have nucleated in bulk) is a function of undercooling
f(ΔTn) is difficult to be found from theoretical considerations alone. It needs to be found experimentally i.e. form a set of experiments e.g.
Method 1 Measurement of cooling curve.
This has been explained in detail in Section 7.2.1. and 7.2.2.
Method 2 Measurement of grain density (optical micrograph of cross section (using Image JÒ/ASTM 562 - 11 Method)) for specimens solidified at various cooling rates .
As soon as grain has nucleated, and its growth can be explained by special modified case of classical nucleation theory (CNT) for BMGMC (A detailed treatment of modified CNT for BMGMC is given in Appendix A) and its distribution can be explained by constitutional supercooling zone/Interdependence theory (propagation of liquid-solid interface/liquid-solid spherical front)—a possibility which is still under investigation by author for suitability for additive manufacturing processes), it grows with an interface velocity which is also a function of undercooling.
8.5.3. Velocity of Growth
Velocity of growth may be written as
In this case, there is no need to determine solidification kinetics of dendrite tip/eutectic (spherical front) interface by cooling curve or grain size but it can be determined by theoretical models developed (by using basic laws of physics)   as applied to BMGMC only under transient condition.
8.5.4. Impingement and Columnar to Equiaxed Transition (CET)
Impingement of grains as they grow is another important phenomenon which for all practical reasons governs the shape of grain after columnar to equiaxed transition (CET) (CET in additive manufacturing is recently explained by Amrita Basak et al.  which is combined with present model and is explained in detail in Appendix B). This phenomenon is not remarkably present in bulk metallic glass and their composites due to their sluggish nature and slightly formed crystal grains as compared to huge glassy matrix. However, despite these drawbacks, this is mainly responsible for equiaxed dendritic grain formation even in glassy alloys, especially in eutectic compositions which is assumed to be the case for present research.
This has been typically treated by
1) Standard J M A K   correction or by
2) Geometrical   or
3) Random grain arrangement models 
These “microscopic” solidification models have been coupled with “macroscopic” transient one-dimensional (1D) heat flow calculations to successfully predict “microstructural features” specially “grain size” at the scale of whole process (part scale)  .
These deterministic models have their following limitations
1) Grain selection
They cannot account for the “grain selection” which occurs
a) Close to mould region/surface giving rise to columnar dendritic microstructure (in case of conventional Cu mould suction casting/twin roll casting (TRC)) or
b) At surface of external inoculant particles (precursors of heterogeneous nucleation) in case of well inoculated melts (present case)) giving rise to onset of columnar dendritic microstructure (at a very small length scale) since they almost neglect any aspect which is related to crystallographic effects.
2) “Equiaxed-Columnar” Transition
They cannot predict the co-called “equiaxed-columnar” transition which occurs very near to mould wall  or variation of transverse size of columnar grains  (also known as columnar dendritic arm branching). This is explained in detail in individual cases for each type of metal (crystal structure).
Case 1: Cubic Metals
It is well established facts that for cubic metals, this “grain selection” is based upon a criterion of best alignment” of the < 100 > crystallographic axes of grain with heat flow direction   . Thus, this method cannot account for this anisotropic behaviour of heat flow. A solution to this problem could be proposed by determining best fit direction by use of recent developments in crystallography and their application to solidification. Edge to edge matching (E2EM): One way is to use Edge-to-Edge Matching (E2EM) technique at inoculant-ductile phase level (in case of Zr based BMGMC) (present research). This gives rise to selection of suitable potent nuclei of certain size and specific preferred orientation (i.e. along a defined easy crystallographic plane (e.g. (001)). If this crystallographic plane direction could be used in conjunction with macroscopic heat flow models, it can give rise to “prediction or selection of grain”. In other words, if matching crystallographic axes (suitable for a potent inoculant selection for CuZr-B2 ductile phase’s preferred precipitation (in case of BMGMC)) could be best aligned with heat flow direction (or heat flow direction could be assigned to this preferred matching crystallographic axes) a best “grain selection: could be determined (one of aim of present research—not done previously elsewhere). This type of phenomena is particularly important in
a) Directional solidification (DS) or
b) Production of single crystal dendritic alloys for aerospace applications or
c) Production of BMGMC by Bridgeman solidification
Note: This is in addition to use of E2EM for selection of potent nuclei
Case 2: BCC Metals
These methods are also ineffective in predicting “equiaxed-columnar” and then “branching of dendrite arms” in bcc metals (i.e. grain selection) as best alignment between heat flow and crystallographic direction is not well known. Only assumptions are possible (i.e. in case of bcc best heat flow direction could be assigned to close packed direction).
Case 3: FCC Metals
These methods are again ineffective in predicting the “equiaxed - columnar”, “CET” and then branching of dendrite arms in fcc metals (i.e. grain selection) as best alignment between heat flow and close packed direction (111) could only be assumed (to a satisfactory qualitative level). More quantitative experimentation is needed to determine best directions along which heat flow occurred or revert to more advanced models.
3) Extension of a grain into an open region of liquid.
They cannot account for extension of a grain into an open region of liquid.
4) Columnar-to-Equiaxed Transition (CET)
Finally, when very fine equiaxed grains at a region very close to mould wall/right at the interface of inoculant and melt are converted to columnar grains, which when grow, there comes a point/plane at which columnar grains gets converted to not so fine equiaxed grains. This point is known as columnar to equiaxed transition (CET). These equiaxed grains finally extend towards centre of casting (wedge shape/melt pool centreline in case of additive manufacturing). CET primarily happens as a result of thermal fluctuations which happen at melt (liquid) and solid (solidified melt) interface which are triggered by solutal effects as well as heat extraction or absorption due to phase changes occurring at a micro-scale (explained in subsequent sections). CET is dominant when thermal gradient is small.
8.6. Evolution of Probabilistic Models
The solution to above four problems is presented first by Brown and Spittle  . They developed probabilistic models. They used monte carlo (MC) procedure for explaining solidification phenomena developed in earlier research . MC method is based upon minimising of interfacial energy (which is practically calculated by using physical properties of material (Zr- and Fe-based BMGMC)) from literature and earlier published data or inference from extrapolation or interpolation of data as needed). Procedurally, this minimisation is achieved by
1) Considering the energy of “unlike sites” (e.g. (a) “liquid/solid sites” or (b) “sites belonging to different grains” and
2) By allowing transition between these states to occur according to randomly generated numbers
By using this method, Brown and Spittle merely able to produce computed 2 D microstructures which resembled very closely to those observed in real micrographic cross section. In particular
a) The selection of grains in the columnar zone and
b) Columnar to equiaxed transition (CET)
were nicely reproduced using this technique also
a) the effect of solute concentration or
b) melt superheat upon the resultant microstructure
was determined “qualitatively” in a nice way. Their quantitative representation was not achieved.
These methods suffer consistently from lack of physical basis and thus cannot be used to analyse quantitatively the effect of various physical phenomena (happening within the phase transformations). For example, to illustrate this, consider the following example.
1) During one monte carlo (MC) time step, Consider N sites where N is number of sites whose evolution is calculated and is chosen from another N (total number) sites. Therefore, not all sites of interest (i.e. those located near to solid-liquid interface) are investigated. This in turn, leads to algorithm predicted grain competition in columnar region, which does not at all reflect the physical mechanisms observed in organic alloys.
2) Furthermore, the results are sensitive to type of Monte Carlo (MC) network itself which is used for computations. Thus, a single powerful model is presented in present work which combines “advantages of probabilistic methods with those of deterministic approaches” to predict more accurately the grain structure in a casting.
8.7. Two-Dimensional Cellular Automaton (CA) Method
For this purpose, for now, a 2D Cellular Automaton model is developed which is based upon physical mechanisms of nucleation and growth of dendritic grains. Its salient features are as follows;
1) Heterogeneous nucleation; which was modelled by means of a nucleation site distribution in deterministic solidification models, is treated in a similar way in present probabilistic approach.
2) If total density of grains which nucleate at a given undercooling is obtained from an average distribution (dc = average (distribution)), the location of these sites is chosen randomly
where are distributions of grains 1, 2, 3 to n.
where n = R (R = Real numbers).
3) Crystallographic orientation of a newly nucleated grain is also taken into account at random.
4) The growth kinetics of (a) dendrite tip and (b) of side branches are also incorporated into the model in such a way that final simulated microstructure is independent of the “cellular automaton network” which is used for computations.
Although, it produces micrographic cross sections very much similar to those already obtained by Brown and Spittle, present model has a “sound physical basis” and can thus reflect effects of (a) cooling rate” or (b) “solute concentration” quantitatively.
8.7.1. Detailed Description (Phase 1—Application of CAFE to Conventional Casting)
Physical background: Consider a BMGMC wedge shape casting as shown in Figure 3 below.
Figure 3(c) is typical 2D cross section of cast eutectic Zr-based BMGMC
(a) (b) (c) (d)
Figure 3. (a) Schematic 3D (b) Optical Micrograph of cross section (etched) (c) 2D Schematic showing regions (d) A specific region (from B2 dendrites) showing B19’ twins (B2 - B19’ TRIP).
solidified in water cooled wedge shape Cu-mould  . Their dendritic grains which have various crystallographic orientations appear as zones of different colours (Figure 3(b)). Most common regions encountered in any casting appear here   and are marked all along cross-section. On the top end of wedge shape ingot coarse grains are present as this region was exposed to air. Its more detailed explanation will follow after characterizing region chronologically from bottom to top.
Bottom region glass: The tip of casting is 100% glass (monolithic BMG). This region is classified as glass and no crystal structure is observed here because cooling rate is maximum here which results in extraction of heat at a very high rate resulting in retaining supercooled liquid state at room temperature.
Bottom region: columnar dendrites: This region marks the beginning of “equiaxed columnar” first transition. This consists of very fine layer in which this transition happens, and then columnar grains grow (primarily) in random 3D orientation) because of still rapid rate of heat transfer which is complemented by sluggish nucleation on growth mechanisms of BMGMC. These grains are not very long as heat flow pattern is somewhat exponential because of wedge shape casting which triggers next transition too quickly before extension of growth as predicted by kinetics. This helps in retaining glassy matrix all throughout the casting. Otherwise 100% crystallisation would have happened.
Bottom region: columnar to equiaxed transition (CET): This is the region in which columnar dendritic grains which have developed/grown to a satisfactory level, transit to equiaxed grains, known as CET. This is triggered by various phenomena such as solute diffusion, solute-solvent partitioning, shape of liquid-solid propagation front, thermal fluctuations happening at the tip of liquid-solid propagating interface.
Fine equiaxed dendrites (CuZr-B2): Once CET happens, equiaxed dendrites are formed all throughout the casting. Only their shape differs. In this region, they are fine sized while in Top Region, their size is even more reduced due to presence of IMCs. Casting scum and other impurities coupled with faster cooling rate from open top (convection and radiation) and side walls (conduction).
NOTE: In case of bulk metallic glass matrix composites (BMGMC) not only inoculant particles serve as sites for heterogeneous nucleation, but grain boundaries also serve this purpose  . Other defects and solidification microstructure also serve as sites for heterogeneous nucleation (their effects in total solidification (nucleation and growth model) are to be considered in final model).
Despite their recent popularity and emblem to be exploited as potential structural engineering material for extreme applications, still, an in-depth understanding of underlaying mechanisms responsible for formation of glass and nucleation of ductile phase during solidification in bulk metallic glass and their composites is at an arm’s length from the limit of satisfaction. Various theories and mechanisms such as phase separation, liquid-liquid transition, frustration, ordering and confusion have been proposed to account for development of ductility in these diverse classes of materials, but none have proved out to be fully satisfactory. These studies have been backed by strong and powerful experimental techniques such as observation under synchrotron light employing container less levitated conditions and observation of nucleation and growth under micro and zero gravity conditions on board international space station. An innovative technique known as additive manufacturing, has recently emerged as competitive method for the large-scale production of metals and alloys including bulk metallic glass matrix composites. It intrinsically, engulfs applicability and implementation of all important phenomena responsible for development of in-situ dendrite reinforced microstructure in these alloys. However, like conventional methods, it also poses challenges to produce high quality components quickly. Evolution of melt pool and microstructure in it are two main important parameters to control. Part and material loss is inevitable to investigate the root cause of these problems. This study is an attempt to bridge both these gaps. A concise review of existing literature on bulk metallic glass and their composites is followed by introduction of coupled microscale simultaneous heat and mass model accounting for development of liquid melt pool in additive manufacturing and microstructure in it. Two-dimensional probabilistic cellular automaton method is combined with one dimensional non-linear heat transfer method based on finite element method to predict heat transfer pattern and microstructure in rapidly solidifying liquid melt pool in additive manufacturing. It is shown systematically, that step-by-step application of model with the use of actual heat and mass transfer (diffusion) coefficients at every step of solidification is the solution of problem. Transient conditions are adopted by implementation of iterative technique using FORTRANÒ and MATLAB languages for programming and commercial simulation package ABAQUSÒ for parallel visual simulation. This is the first study of its kind emphasizing the need and importance of modelling and simulation in predicting solidification phenomena (heat transfer pattern and microstructure based on solute diffusion and capillary action) in additive manufacturing of bulk metallic glass and their composites.
Author greatly acknowledges the time spent in reading and correcting the manuscript by Dr. Mark Gibson and for his helpful suggestions and positive constructive criticism. He also thanks and acknowledges the moral support, always welcoming attitude and time for technical discussions provided by Prof. Milan Brandt.
Rafique, M.M.A. (2018) Modelling and Simulation of Solidification Phenomena during Additive Manufacturing of Bulk Metallic Glass Matrix Composites (BMGMC)—A Brief Review and Introduction of Technique. Journal of Encapsulation and Adsorption Sciences, 8, 67-116. https://doi.org/10.4236/jeas.2018.82005
Heterogeneous nucleation and growth in very fluid alloys (as per CNT) .
Heterogeneous Nucleation rate per unit volume is defined as
where Ns = No. of atoms in contact with substrate; υ = vibrational frequency; ΔGc = Activation energy for nucleation (Critical energy of nucleus formation (i.e. creation of liquid-solid interface); ΔGd = Activation energy of diffusion (Diffusional activation energy).
Rearranging Equation (15) using definition of υ vibrational frequency
N can be neglected as during initial stages there is no nucleation event.
According to CNT, a minimum energy value is needed to create a solid-liquid interface eventually leading to stable nuclei out of melt. This is known as “activation energy”. This activation energy is the energy to overcome ΔG*—the energy barrier to nucleation. Now, as solid-liquid interface grows to form stable nuclei, atoms must be transported through liquid thus another temperature dependent activation energy must be overcome known as ΔGd (activation energy for diffusion).
The net effect is that CNT predicts a nucleation rate (I) given by
It is the nature of difference between ΔG* and ΔGd that dictates whether solidification will be crystalline or glassy. For crystalline solids, ΔGd has a significant value while for glassy solids there is no diffusion thus ΔGd can be neglected. Thus
where kβ is constant dictated by nature and type of liquid composition and measured experimentally. ΔGd is also zero in case of small undercooling (i.e. well inoculated liquids/multicomponent alloys (Metallic Glasses inoculated with/without potent nuclei (present research) .
1) Vibrational frequency
Definition used in Equation (7).
2) The difference between frequency and rate is that frequency is “occurrence of an event per unit time” while rate is total number of that event (in terms of numerical value) per unit time.
Thus, from Equation (6)
No = Total number of heterogeneous substrate particles originally available per unit volume; N = Number that have already nucleated; Io = constant.
Value of Io can be calculated from Equation (5) using another term known as “liquid diffusion coefficient”.
where a = atomic diameter = 0.4 nm; υ = frequency
which gives for small values of undercooling (well inoculated melts/multicomponent alloys)
where ΔT = undercooling.
Thus, Nucleation rate is Equation (6)
where u is a constant
The value of u can be measured from
Method 1: T (heterogeneous nucleation temperature). This is defined as temperature, where there is an initial nucleation rate of one nucleus/cm3/sec.
Method 2: Second method to calculate u is
Taking natural log of equation (15) both sides
because lnI and lnN1020 can be neglected
where ΔTN = undercooling at heterogeneous nucleation temperature.
Time is user defined input and temperature comes from user defined value initially as well. Then its every new value is assigned back to Equation 5. With temperature and time, k changes and assigned back to Equation (15). Also, with time, υ (vibrational frequency) changes and assigned back to original equation (5). Similarly, the value of u also changes with time and temperature. Below table (Table 1) summarises the values which are user defined and which change as a function of transience as programs runs.
1) In BMGMCs, in some cases due to slow motion of large atoms, only nucleation happens, and growth never happens. In these cases, a new phenomenon known as soft impingement effects of crystals must be considered. These could be solutal/thermal. However, this is quite rare.
2) In general, in case of BMGMC, classical nucleation theory (CNT) cannot be applied alone to describe nucleation and growth.
3) Constitutional supercooling zone (CSZ) and Interdependence models can
Table 1. Summary of user defined, and program determined functions used in CNT modified for BMGMC.
not be applied because of very high (η) viscosity of BMGMCs (and their sluggish nature). CSZ and Interdependence theories are for less viscous/more fluid alloys.
4) A new concept, known as complex inter diffusion tensor  is much more helpful to describe nucleation and growth in BMGMCs.
5) Fick’s law in its native form (i.e. linear form) is not sufficient.
1) Special case of growth of “Columnar microstructures”
The growth of columnar dendrites, which is initiated by nuclei that form at the mould interface (Cu mould casting/Twin roll casting of BMGMC) (only if constitutional supercooling zone (CSZ) is suppressed—not present case) is usually simulated in a much simpler way. Again, in this case, there is no need to use cooling curve measurements or grain size measurement but same growth kinetics models   can be used to determine.
Undercooling of eutectic front (ΔTn eutectic) or
Undercooling of dendrite tips (ΔTn dendrite tip) as well as
Undercooling of lamellae or dendrite trunk spacing (ΔTn lamellae/trunk spacing)
This undercooling is determined by direct measurement of
a) Thermal gradient and
b) Speed of corresponding isotherm (eutectic or liquidus respectively i.e. speed of eutectic isotherm and speed of liquidus isotherm)
The later values are obtained from a macroscopic (part scale) heat flow calculations  . The secondary arm spacing of both equiaxed and columnar dendritic microstructures are deduced from a local solidification time.
2) Columnar structure growth in well inoculated BMGMC:
Growth of columnar dendrites can also occur at surface of external inoculants (well inoculated deeply undercooled melts—present case of BMGMC development. However, it should also be noted that another condition for growth of columnar dendrite to occur is suppression of CSZ which is clash with aforementioned condition for onset of this phenomena at external potent nuclei of inoculant. That’s why; still there is dispute about application of this concept to deeply undercooled well inoculated melts (BMGMC) whose solution is under investigation.
3) Columnar to equiaxed transition (CET) 
Growth rate of solid-liquid interface
where S = Scan speed
Temperature gradient parallel to dendrite growth direction can be calculated using
where ψ = Angle between “normal vector” and “possible dendrite growth orientation” at the solid-liquid interface. This is evaluated by CFX-Post in AnsysÒ.
A modification known as Rappaz modification is applied to predict CET. This is as follows
where Vhkl = dendrite growth velocity = ; n = material constant determined from literature  ; φ = equiaxed fraction (critical value = 0.066%); No = nucleation density; ΔTtip = tip undercooling; ΔTn = nucleation undercooling.
This will be incorporated in present model at point where CET is determined. However, this model does not give true 3D representation output.
NOTE: In general Phase Field (PF) Methods are for microstructure evolution (its type (planer front, spherical front), morphology (precipitates and dendrites)) while Cellular Automaton (CA) Methods are for grain size determination (equiaxed/columnar dendritic) and its prediction. If both are combined   , it is possible to get full map of microstructure evolution and grain size.
Below a comparison of “strengths and capabilities” and “evolution of different theories over time” which have enabled a better understanding of nucleation and growth phenomena in bulk metallic glass matrix composites, is tabulated. The aim is to present reader with a concise smart workable data for first hand use and reference for solving nucleation and growth problems in bulk metallic glass matrix composites by modelling and simulation. This will help professional programmer, working engineer and a researcher to effectively find previously done research till now with its strengths and capabilities at one platform.
1) Strengths and Capabilities
Below (Table 2) a comparison of strength, capabilities and shortcomings of both deterministic and probabilistic methods are described. It highlights and chalks out parameters and certain segments of each technique which could possibility advantageously used over others for modelling and simulation of bulk metallic glass matrix composites.
Note: N/A is an abbreviation to “Not Applicable”.
2) Evolution of theories (Table 3)
Table 2. Comparison of strength and capabilities of modelling and simulation techniques as applied to nucleation and growth problem of bulk metallic glass matrix composites (BMGMCS).
Table 3. Evolution of theories of modelling and simulation as applied to nucleation and growth problem of bulk metallic glass matrix composites (BMGMCS).
 Güntherodt, H.J. (1977) Metallic Glasses. In: Treusch, J., Ed., Festkörperprobleme 17: Plenary Lectures of the Divisions “Semiconductor Physics” “Metal Physics” “Low Temperature Physics” “Thermodynamics and Statistical Physics” “Crystallography” “Magnetism” “Surface Physics” of the German Physical Society Münster, Springer, Berlin, Heidelberg, 25-53.
 Shi, Y. and Falk, M.L. (2006) Does Metallic Glass Have a Backbone? The Role of Percolating Short Range Order in Strength and Failure. Scripta Materialia, 54, 381-386.
 Zhang, F., et al. (2014) Composition-Dependent Stability of the Medium-Range Order Responsible for Metallic Glass Formation. Acta Materialia, 81, 337-344.
 Cheng, Y.Q., Ma, E. and Sheng, H.W. (2009) Atomic Level Structure in Multicomponent Bulk Metallic Glass. Physical Review Letters, 102, Article ID: 245501.
 Lan, S., et al. (2016) Structural Crossover in a Supercooled Metallic Liquid and the Link to a Liquid-to-Liquid Phase Transition. Applied Physics Letters, 108, Article ID: 211907.
 Akhtar, D., Cantor, B. and Cahn, R.W. (1982) Measurements of Diffusion Rates of Au in Metal-Metal and Metal-Metalloid Glasses. Acta Metallurgica, 30, 1571-1577.
 Lu, Z.P., Liu, Y. and Liu, C.T. (2008) Evaluation of Glass-Forming Ability. In: Miller, M. and Liaw, P., Eds., Bulk Metallic Glasses, Springer, Boston, MA, 87-115.
 Park, E.S. and Kim, D.H. (2005) Design of Bulk Metallic Glasses with High Glass Forming Ability and Enhancement of Plasticity in Metallic Glass Matrix Composites: A Review. Metals and Materials International, 11, 19-27.
 Park, E.S., Chang, H.J. and Kim, D.H. (2008) Effect of Addition of Be on Glass-Forming Ability, Plasticity and Structural Change in Cu-Zr Bulk Metallic Glasses. Acta Materialia, 56, 3120-3131.
 Guo, G.-Q., Wu, S.-Y. and Yang, L. (2016) Structural Origin of the Enhanced Glass-Forming Ability Induced by Microalloying Y in the ZrCuAl Alloy. Metals, 6, 67.
 Cheng, Y.Q., Ma, E. and Sheng, H.W. (2008) Alloying Strongly Influences the Structure, Dynamics, and Glass Forming Ability of Metallic Supercooled Liquids. Applied Physics Letters, 93, Article ID: 111913.
 Jia, P., et al. (2006) A New Cu-Hf-Al Ternary Bulk Metallic Glass with High Glass Forming Ability and Ductility. Scripta Materialia, 54, 2165-2168.
 Miracle, D.B., et al. (2010) An Assessment of Binary Metallic Glasses: Correlations between Structure, Glass Forming Ability and Stability. International Materials Reviews, 55, 218-256.
 Lu, Z. and Liu, C.T. (2004) A New Approach to Understanding and Measuring Glass Formation in Bulk Amorphous Materials. Intermetallics, 12, 1035-1043.
 Kim, Y.C., et al. (2003) Glass Forming Ability and Crystallization Behavior of Ti-Based Amorphous Alloys with High Specific Strength. Journal of Non-Crystalline Solids, 325, 242-250.
 Wu, J., et al. (2014) New Insight on Glass-Forming Ability and Designing Cu-Based Bulk Metallic Glasses: The Solidification Range Perspective. Materials & Design, 61, 199-202.
 Shen, T.D., Sun, B.R. and Xin, S.W. (2015) Effects of Metalloids on the Thermal Stability and Glass Forming Ability of Bulk Ferromagnetic Metallic Glasses. Journal of Alloys and Compounds, 631, 60-66.
 Li, P., et al. (2014) Glass Forming Ability, Thermodynamics and Mechanical Properties of Novel Ti-Cu-Ni-Zr-Hf Bulk Metallic Glasses. Materials & Design, 53, 145-151.
 Li, P., et al. (2012) Glass Forming Ability and Thermodynamics of New Ti-Cu-Ni-Zr Bulk Metallic Glasses. Journal of Non-Crystalline Solids, 358, 3200-3204.
 Fan, C., et al. (2001) Effects of Nb Addition on Icosahedral Quasicrystalline Phase Formation and Glass-Forming Ability of Zr-Ni-Cu-Al Metallic Glasses. Applied Physics Letters, 79, 1024-1026.
 Xu, D., Duan, G. and Johnson, W.L. (2004) Unusual Glass-Forming Ability of Bulk Amorphous Alloys Based on Ordinary Metal Copper. Physical Review Letters, 92, Article ID: 245504.
 Zhang, K., et al. (2013) Computational Studies of the Glass-Forming Ability of Model Bulk Metallic Glasses. The Journal of Chemical Physics, 139, Article ID: 124503.
 Amokrane, S., Ayadim, A. and Levrel, L. (2015) Structure of the Glass-Forming Metallic Liquids by Ab-Initio and Classical Molecular Dynamics, a Case Study: Quenching the Cu60Ti20Zr20 Alloy. Journal of Applied Physics, 118, Article ID: 194903.
 Hays, C.C., Kim, C.P. and Johnson, W.L. (2000) Microstructure Controlled Shear Band Pattern Formation and Enhanced Plasticity of Bulk Metallic Glasses Containing in Situ Formed Ductile Phase Dendrite Dispersions. Physical Review Letters, 84, 2901-2904.
 Wu, Y., et al. (2014) Designing Bulk Metallic Glass Composites with Enhanced Formability and Plasticity. Journal of Materials Science & Technology, 30, 566-575.
 Jang, D. and Greer, J.R. (2010) Transition from a Strong-Yet-Brittle to a Stronger-and-Ductile State by Size Reduction of Metallic Glasses. Nature Materials, 9, 215-219.
 Choi-Yim, H., et al. (2002) Processing, Microstructure and Properties of Ductile Metal Particulate Reinforced Zr57Nb5Al10Cu15.4Ni12.6 Bulk Metallic Glass Composites. Acta Materialia, 50, 2737-2745.
 Lee, M.L., Li, Y. and Schuh, C.A. (2004) Effect of a Controlled Volume Fraction of Dendritic Phases on Tensile and Compressive Ductility in La-Based Metallic Glass Matrix Composites. Acta Materialia, 52, 4121-4131.
 Wu, D.Y., et al. (2016) Glass-Forming Ability, Thermal Stability of B2 CuZr Phase, and Crystallization Kinetics for Rapidly Solidified Cu-Zr-Zn Alloys. Journal of Alloys and Compounds, 664, 99-108.
 Kim, C., et al. (2011) Realization of High Tensile Ductility in a Bulk Metallic Glass Composite by the Utilization of Deformation-Induced Martensitic Transformation. Scripta Materialia, 65, 304-307.
 Gao, W.-H., et al. (2015) Effects of Co and Al Addition on Martensitic Transformation and Microstructure in ZrCu-Based Shape Memory Alloys. Transactions of Nonferrous Metals Society of China, 25, 850-855.
 Zhai, H., Wang, H. and Liu, F. (2016) A Strategy for Designing Bulk Metallic Glass Composites with Excellent Work-Hardening and Large Tensile Ductility. Journal of Alloys and Compounds, 685, 322-330.
 Tian, F., Li, Z. and Song, J. (2016) Solidification of Laser Deposition Shaping for TC4 Alloy Based on Cellular Automation. Journal of Alloys and Compounds, 676, 542-550.
 Zhou, M., Rosakis, A.J. and Ravichandran, G. (1998) On the Growth of Shear Bands and Failure-Mode Transition in Prenotched Plates: A Comparison of Singly and Doubly Notched Specimens. International Journal of Plasticity, 14, 435-451.
 Pan, D., et al. (2008) Experimental Characterization of Shear Transformation Zones for Plastic Flow of Bulk Metallic Glasses. Proceedings of the National Academy of Sciences, 105, 14769-14772.
 Abdeljawad, F., Fontus, M. and Haataja, M. (2011) Ductility of Bulk Metallic Glass Composites: Microstructural Effects. Applied Physics Letters, 98, Article ID: 031909.
 Pekarskaya, E., Kim, C.P. and Johnson, W.L. (2001) In Situ Transmission Electron Microscopy Studies of Shear Bands in a Bulk Metallic Glass Based Composite. Journal of Materials Research, 16, 2513-2518.
 Zhang, Q., Zhang, H.F., Zhu, Z.W. and Hu, Z.Q. (2005) Formation of High Strength In-Situ Bulk Metallic Glass Composite with Enhanced Plasticity in Cu50Zr47.5Ti2.5 Alloy. Materials Transactions, 46, 730-733.
 Zhu, Z., et al. (2010) Ta-Particulate Reinforced Zr-Based Bulk Metallic Glass Matrix Composite with Tensile Plasticity. Scripta Materialia, 62, 278-281.
 Wu, F.F., et al. (2007) Effect of Annealing on the Mechanical Properties and Fracture Mechanisms of a Zr56.2Ti13.8Nb5.0Cu6.9Ni5.6Be12.5 Bulk-Metallic-Glass Composite. Physical Review B, 75, Article ID: 134201.
 Basu, J., et al. (2003) Microstructure and Mechanical Properties of a Partially Crystallized La-Based Bulk Metallic Glass. Philosophical Magazine, 83, 1747-1760.
 Fan, C., et al. (2006) Properties of As-Cast and Structurally Relaxed Zr-Based Bulk Metallic Glasses. Journal of Non-Crystalline Solids, 352, 174-179.
 Gu, J., et al. (2013) Effects of Annealing on the Hardness and Elastic Modulus of a Cu36Zr48Al8Ag8 Bulk Metallic Glass. Materials & Design, 47, 706-710.
 Tan, J., et al. (2013) Correlation between Internal States and Strength in Bulk Metallic Glass, in PRICM. John Wiley & Sons, Inc., New York, 3199-3206.
 Song, K.K., et al. (2011) Strategy for Pinpointing the Formation of B2 CuZr in Metastable CuZr-Based Shape Memory Alloys. Acta Materialia, 59, 6620-6630.
 Ding, J., et al. (2014) Large-Sized CuZr-Based Bulk Metallic Glass Composite with Enhanced Mechanical Properties. Journal of Materials Science & Technology, 30, 590-594.
 Jiang, F., et al. (2007) Microstructure Evolution and Mechanical Properties of Cu46Zr47Al7 Bulk Metallic Glass Composite Containing CuZr Crystallizing Phases. Materials Science and Engineering: A, 467, 139-145.
 Liu, J., et al. (2010) Microstructure and Compressive Properties of In-Situ Martensite CuZr Phase Reinforced ZrCuNiAl Metallic Glass Matrix Composite. Materials Transactions, 51, 1033-1037.
 Liu, Z., et al. (2012) Microstructural Tailoring and Improvement of Mechanical Properties in CuZr-Based Bulk Metallic Glass Composites. Acta Materialia, 60, 3128-3139.
 Liu, Z.Q., et al. (2014) Microstructural Percolation Assisted Breakthrough of Trade-Off between Strength and Ductility in CuZr-Based Metallic Glass Composites. Scientific Reports, 4, 4167.
 Schryvers, D., et al. (1997) Unit Cell Determination in CuZr Martensite by Electron Microscopy and X-Ray Diffraction. Scripta Materialia, 36, 1119-1125.
 Seo, J.W. and Schryvers, D. (1998) TEM Investigation of the Microstructure and Defects of CuZr Martensite. Part I: Morphology and Twin Systems. Acta Materialia, 46, 1165-1175.
 Seo, J.W. and Schryvers, D. (1998) TEM Investigation of the Microstructure and Defects of CuZr Martensite. Part II: Planar Defects. Acta Materialia, 46, 1177-1183.
 Song, K.K., et al. (2013) Correlation between the Microstructures and the Deformation Mechanisms of CuZr-Based Bulk Metallic Glass Composites. AIP Advances, 3, Article ID: 012116.
 Sun, L., et al. (2016) Phase Separation and Microstructure Evolution of Zr48Cu36Ag8Al8 Bulk Metallic Glass in the Supercooled Liquid Region. Rare Metal Materials and Engineering, 45, 567-570.
 Antonowicz, J., et al. (2008) Early Stages of Phase Separation and Nanocrystallization in Al-Rare Earth Metallic Glasses Studied Using SAXS/WAXS and HRTEM Methods. Reviews on Advanced Materials Science, 18, 454-458.
 Michalik, S., et al. (2014) Structural Modifications of Swift-Ion-Bombarded Metallic Glasses Studied by High-Energy X-Ray Synchrotron Radiation. Acta Materialia, 80, 309-316.
 Mu, J., et al. (2013) In Situ High-Energy X-Ray Diffraction Studies of Deformation-Induced Phase Transformation in Ti-Based Amorphous Alloy Composites Containing Ductile Dendrites. Acta Materialia, 61, 5008-5017.
 Oh, Y.S., et al. (2011) Microstructure and Tensile Properties of High-Strength High-Ductility Ti-Based Amorphous Matrix Composites Containing Ductile Dendrites. Acta Materialia, 59, 7277-7286.
 Wang, Y., et al. (2014) Investigation of the Microcrack Evolution in a Ti-Based Bulk Metallic Glass Matrix Composite. Progress in Natural Science: Materials International, 24, 121-127.
 Wang, Y.S., et al. (2014) The Role of the Interface in a Ti-Based Metallic Glass Matrix Composite with in Situ Dendrite Reinforcement. Surface and Interface Analysis, 46, 293-296.
 Zhang, T., et al. (2014) Dendrite Size Dependence of Tensile Plasticity of in Situ Ti-Based Metallic Glass Matrix Composites. Journal of Alloys and Compounds, 583, 593-597.
 Chu, M.Y., et al. (2015) Quasi-Static and Dynamic Deformation Behaviors of an in-Situ Ti-Based Metallic Glass Matrix Composite. Journal of Alloys and Compounds, 640, 305-310.
 Hofmann, D.C., et al. (2008) New Processing Possibilities for Highly Toughened Metallic Glass Matrix Composites with Tensile Ductility. Scripta Materialia, 59, 684-687.
 Cheng, J.L. and Chen, G. (2013) Glass Formation of Zr-Cu-Ni-Al Bulk Metallic Glasses Correlated with L → Zr2Cu + ZrCu Pseudo Binary Eutectic Reaction. Journal of Alloys and Compounds, 577, 451-455.
 Biffi, C.A., Figini Albisetti, A. and Tuissi, A. (2013) CuZr Based Shape Memory Alloys: Effect of Cr and Co on the Martensitic Transformation. Materials Science Forum, Trans Tech Publications, Zurich.
 Inoue, A., Zhang, T. and Masumoto, T. (1990) Zr-Al-Ni Amorphous Alloys with High Glass Transition Temperature and Significant Supercooled Liquid Region. Materials Transactions, JIM, 31, 177-183.
 Chen, G., et al. (2009) Enhanced Plasticity in a Zr-Based Bulk Metallic Glass Composite with in Situ Formed Intermetallic Phases. Applied Physics Letters, 95, Article ID: 081908.
 Jeon, C., et al. (2015) Effects of Effective Dendrite Size on Tensile Deformation Behavior in Ti-Based Dendrite-Containing Amorphous Matrix Composites Modified from Ti-6Al-4V Alloy. Metallurgical and Materials Transactions A, 46, 235-250.
 Gibson, I., Rosen, W.D. and Stucker, B. (2010) Development of Additive Manufacturing Technology, in Additive Manufacturing Technologies: Rapid Prototyping to Direct Digital Manufacturing. Springer, Boston, MA, 36-58.
 Spears, T.G. and Gold, S.A. (2016) In-Process Sensing in Selective Laser Melting (SLM) Additive Manufacturing. Integrating Materials and Manufacturing Innovation, 5, 1-25.
 Li, X.P., et al. (2016) Selective Laser Melting of Zr-Based Bulk Metallic Glasses: Processing, Microstructure and Mechanical Properties. Materials & Design, 112, 217-226.
 Zheng, B., et al. (2009) Processing and Behavior of Fe-Based Metallic Glass Components via Laser-Engineered Net Shaping. Metallurgical and Materials Transactions A, 40, 1235-1245.
 Olakanmi, E.O., Cochrane, R.F. and Dalgarno, K.W. (2015) A Review on Selective Laser Sintering/Melting (SLS/SLM) of Aluminium Alloy Powders: Processing, Microstructure, and Properties. Progress in Materials Science, 74, 401-477.
 Li, Y. and Gu, D. (2014) Thermal Behavior during Selective Laser Melting of Commercially Pure Titanium Powder: Numerical Simulation and Experimental Study. Additive Manufacturing, 1-4, 99-109.
 Romano, J., et al. (2015) Temperature Distribution and Melt Geometry in Laser and Electron-Beam Melting Processes—A Comparison among Common Materials. Additive Manufacturing, 8, 1-11.
 Sun, H. and Flores, K.M. (2010) Microstructural Analysis of a Laser-Processed Zr-Based Bulk Metallic Glass. Metallurgical and Materials Transactions A, 41, 1752-1757.
 Zhang, Y., et al. (2015) Microstructural Analysis of Zr55Cu30Al10Ni5 Bulk Metallic Glasses by Laser Surface Remelting and Laser Solid Forming. Intermetallics, 66, 22-30.
 Baufeld, B., Brandl, E. and van der Biest, O. (2011) Wire Based Additive Layer Manufacturing: Comparison of Microstructure and Mechanical Properties of Ti-6Al-4V Components Fabricated by Laser-Beam Deposition and Shaped Metal Deposition. Journal of Materials Processing Technology, 211, 1146-1158.
 Kawahito, Y., et al. (2008) High-Power Fiber Laser Welding and Its Application to Metallic Glass Zr55Al10Ni5Cu30. Materials Science and Engineering: B, 148, 105-109.
 Wang, H.S., et al. (2010) Combination of a Nd:YAG Laser and a Liquid Cooling Device to (Zr53Cu30Ni9Al8)Si0.5 Bulk Metallic Glass Welding. Materials Science and Engineering: A, 528, 338-341.
 Wang, H.-S., et al. (2011) The Effects of Initial Welding Temperature and Welding Parameters on the Crystallization Behaviors of Laser Spot Welded Zr-Based Bulk Metallic Glass. Materials Chemistry and Physics, 129, 547-552.
 Acharya, R. and Das, S. (2015) Additive Manufacturing of IN100 Superalloy Through Scanning Laser Epitaxy for Turbine Engine Hot-Section Component Repair: Process Development, Modeling, Microstructural Characterization, and Process Control. Metallurgical and Materials Transactions A, 46, 3864-3875.
 Harooni, A., et al. (2016) Processing Window Development for Laser Cladding of Zirconium on Zirconium Alloy. Journal of Materials Processing Technology, 230, 263-271.
 Wu, X., Xu, B. and Hong, Y. (2002) Synthesis of Thick Ni66Cr5Mo4Zr6P15B4 Amorphous Alloy Coating and Large Glass-Forming Ability by Laser Cladding. Materials Letters, 56, 838-841.
 Yue, T.M. and Su, Y.P. (2008) Laser Cladding of SiC Reinforced Zr65Al7.5Ni10Cu17.5 Amorphous Coating on Magnesium Substrate. Applied Surface Science, 255, 1692-1698.
 Zhang, P., et al. (2011) Synthesis of Fe-Ni-B-Si-Nb Amorphous and Crystalline Composite Coatings by Laser Cladding and Remelting. Surface and Coatings Technology, 206, 1229-1236.
 Zhu, Q., et al. (2007) Synthesis of Fe-Based Amorphous Composite Coatings with Low Purity Materials by Laser Cladding. Applied Surface Science, 253, 7060-7064.
 Emmerich, H. (2009) Phase-Field Modelling for Metals and Colloids and Nucleation Therein—An Overview. Journal of Physics: Condensed Matter, 21, Article ID: 464103.
 Emmerich, H., et al. (2012) Phase-Field-Crystal Models for Condensed Matter Dynamics on Atomic Length and Diffusive Time Scales: An Overview. Advances in Physics, 61, 665-743.
 Gránásy, L., et al. (2014) Phase-Field Modeling of Polycrystalline Solidification: From Needle Crystals to Spherulites—A Review. Metallurgical and Materials Transactions A, 45, 1694-1719.
 Wang, T. and Napolitano, R.E. (2012) A Phase-Field Model for Phase Transformations in Glass-Forming Alloys. Metallurgical and Materials Transactions A, 43, 2662-2668.
 Rappaz, M. and Gandin, C.A. (1993) Probabilistic Modelling of Microstructure Formation in Solidification Processes. Acta Metallurgica et Materialia, 41, 345-360.
 Charbon, C. and Rappaz, M. (1993) 3D Probabilistic Modelling of Equiaxed Eutectic Solidification. Modelling and Simulation in Materials Science and Engineering, 1, 455.
 Gandin, C.A., Schaefer, R.J. and Rappax, M. (1996) Analytical and Numerical Predictions of Dendritic Grain Envelopes. Acta Materialia, 44, 3339-3347.
 Gandin, C.A. and Rappaz, M. (1994) A Coupled Finite Element-Cellular Automaton Model for the Prediction of Dendritic Grain Structures in Solidification Processes. Acta Metallurgica et Materialia, 42, 2233-2246.
 Chen, S., Guillemot, G. and Gandin, C.-A. (2014) 3D Coupled Cellular Automaton (CA)-Finite Element (FE) Modeling for Solidification Grain Structures in Gas Tungsten Arc Welding (GTAW). ISIJ International, 54, 401-407.
 Chen, S. (2014) Three Dimensional Cellular Automaton-Finite Element (CAFE) Modeling for the Grain Structures Development in Gas Tungsten/Metal Arc Welding Processes. Ecole Nationale Supérieure des Mines de Paris.
 Tsai, D.-C. and Hwang, W.-S. (2011) A Three Dimensional Cellular Automaton Model for the Prediction of Solidification Morphologies of Brass Alloy by Horizontal Continuous Casting and Its Experimental Verification. Materials Transactions, 52, 787-794.
 Zinoviev, A., et al. (2016) Evolution of Grain Structure during Laser Additive Manufacturing. Simulation by a Cellular Automata Method. Materials & Design, 106, 321-329.
 Zhou, X., et al. (2016) Simulation of Microstructure Evolution during Hybrid Deposition and Micro-Rolling Process. Journal of Materials Science, 51, 6735-6749.
 Zhang, J., et al. (2013) Probabilistic Simulation of Solidification Microstructure Evolution during Laser-Based Metal Deposition. Proceedings of 2013 Annual International Solid Freeform Fabrication Symposium: An Additive Manufacturing Conference, 739-748.
 Christian, J.W. (2002) Chapter 10—The Classical Theory of Nucleation. In: The Theory of Transformations in Metals and Alloys, Pergamon, Oxford, 422-479.
 Inoue, A. and Zhang, T. (1995) Fabrication of Bulky Zr-Based Glassy Alloys by Suction Casting into Copper Mold. Materials Transactions, JIM, 36, 1184-1187.
 Browne, D.J., Kovacs, Z. and Mirihanage, W.U. (2009) Comparison of Nucleation and Growth Mechanisms in Alloy Solidification to Those in Metallic Glass Crystallisation—Relevance to Modeling. Transactions of the Indian Institute of Metals, 62, 409-412.
 Li, H.Q., Yan, J.H. and Wu, H.J. (2009) Modelling and Simulation of Bulk Metallic Glass Production Process with Suction Casting. Materials Science and Technology, 25, 425-431.
 Markl, M. and Körner, C. (2016) Multiscale Modeling of Powder Bed-Based Additive Manufacturing. Annual Review of Materials Research, 46, 93-123.
 Fan, Z., et al. (2007) Numerical Simulation of the Evolution of Solidification Microstructure in Laser Deposition. Proceedings of the 18th Annual Solid Freeform Fabrication Symposium, Austin, TX, 6-8 August 2007, 256-265.
 Zhu, M.F., Lee, S.Y. and Hong, C.P. (2004) Modified Cellular Automaton Model for the Prediction of Dendritic Growth with Melt Convection. Physical Review E, 69, Article ID: 061610.
 Zheng, B., et al. (2008) Thermal Behavior and Microstructural Evolution during Laser Deposition with Laser-Engineered Net Shaping: Part I. Numerical Calculations. Metallurgical and Materials Transactions A, 39, 2228-2236.
 King, W., et al. (2015) Overview of Modelling and Simulation of Metal Powder Bed Fusion Process at Lawrence Livermore National Laboratory. Materials Science and Technology, 31, 957-968.
 Khairallah, S.A., et al. (2016) Laser Powder-Bed Fusion Additive Manufacturing: Physics of Complex Melt Flow and Formation Mechanisms of Pores, Spatter, and Denudation Zones. Acta Materialia, 108, 36-45.
 Rafique, M.M.A. and Iqbal, J. (2009) Modeling and Simulation of Heat Transfer Phenomena during Investment Casting. International Journal of Heat and Mass Transfer, 52, 2132-2139.
 Bennon, W.D. and Incropera, F.P. (1987) A Continuum Model for Momentum, Heat and Species Transport in Binary Solid-Liquid Phase Change Systems—I. Model Formulation. International Journal of Heat and Mass Transfer, 30, 2161-2170.
 Voller, V.R., Brent, A.D. and Prakash, C. (1989) The Modelling of Heat, Mass and Solute Transport in Solidification Systems. International Journal of Heat and Mass Transfer, 32, 1719-1731.
 Ganesan, S. and Poirier, D.R. (1990) Conservation of Mass and Momentum for the Flow of Interdendritic Liquid during Solidification. Metallurgical Transactions B, 21, 173-181.
 Lee, D.M., et al. (2012) A Deep Eutectic Point in Quaternary Zr-Ti-Ni-Cu System and Bulk Metallic Glass Formation near the Eutectic Point. Intermetallics, 21, 67-74.
 Basak, A., Acharya, R. and Das, S. (2016) Additive Manufacturing of Single-Crystal Superalloy CMSX-4 Through Scanning Laser Epitaxy: Computational Modeling, Experimental Process Development, and Process Parameter Optimization. Metallurgical and Materials Transactions A, 47, 3845-3859.
 Rappaz, M. and Blank, E. (1986) Simulation of Oriented Dendritic Microstructures Using the Concept of Dendritic Lattice. Journal of Crystal Growth, 74, 67-76.
 Spittle, J.A. and Brown, S.G.R. (1989) Computer Simulation of the Effects of Alloy Variables on the Grain Structures of Castings. Acta Metallurgica, 37, 1803-1810.
 Brown, S.G.R. and Spittle, J.A. (1989) Computer Simulation of Grain Growth and Macrostructure Development during Solidification. Materials Science and Technology, 5, 362-368.
 Song, H., et al. (2016) Simulation Study of Heterogeneous Nucleation at Grain Boundaries during the Austenite-Ferrite Phase Transformation: Comparing the Classical Model with the Multi-Phase Field Nudged Elastic Band Method. Metallurgical and Materials Transactions A, 1-9.
 Acharya, R. and Das, S. (2015) Additive Manufacturing of IN100 Superalloy Through Scanning Laser Epitaxy for Turbine Engine Hot-Section Component Repair: Process Development, Modeling, Microstructural Characterization, and Process Control. Metallurgical and Materials Transactions A—Physical Metallurgy and Materials Science, 46a, 3864-3875.
 Choudhury, A., et al. (2012) Comparison of Phase-Field and Cellular Automaton Models for Dendritic Solidification in Al-Cu Alloy. Computational Materials Science, 55, 263-268.
 Zaeem, M.A., Yin, H. and Felicelli, S.D. (2012) Comparison of Cellular Automaton and Phase Field Models to Simulate Dendrite Growth in Hexagonal Crystals. Journal of Materials Science & Technology, 28, 137-146.
 Tan, W., Bailey, N.S. and Shin, Y.C. (2011) A Novel Integrated Model Combining Cellular Automata and Phase Field Methods for Microstructure Evolution during Solidification of Multi-Component and Multi-Phase Alloys. Computational Materials Science, 50, 2573-2585.
 St John, D.H., et al. (2011) The Interdependence Theory: The Relationship between Grain Formation and Nucleant Selection. Acta Materialia, 59, 4907-4921.
 Acharya, R., et al. (2012) Computational Modeling and Experimental Validation of Microstructural Development in Superalloy Cmsx-4 Processed through Scanning Laser Epitaxy. Proceedings of the 8th Pacific Rim International Congress on Advanced Materials and Processing, 2711-2722.
 Gandin, C.-A., et al. (1999) A Three-Dimensional Cellular Automation-Finite Element Model for the Prediction of Solidification Grain Structures. Metallurgical and Materials Transactions A, 30, 3153-3165.
 Blázquez, J.S., et al. (2008) Instantaneous Growth Approximation Describing the Nanocrystallization Process of Amorphous Alloys: A Cellular Automata Model. Journal of Non-Crystalline Solids, 354, 3597-3605.
 Gránásy, L. (1993) Quantitative Analysis of the Classical Nucleation Theory on Glass-Forming Alloys. Journal of Non-Crystalline Solids, 156, 514-518.
 Nandi, U.K., et al. (2016) Composition Dependence of the Glass Forming Ability in Binary Mixtures: The Role of Demixing Entropy. The Journal of Chemical Physics, 145, Article ID: 034503.
 Raabe, D. (2004) Overview of the Lattice Boltzmann Method for Nano- and Microscale Fluid Dynamics in Materials Science and Engineering. Modelling and Simulation in Materials Science and Engineering, 12, R13.
 Sun, D.K., et al. (2011) Modelling of Dendritic Growth in Ternary Alloy Solidification with Melt Convection. International Journal of Cast Metals Research, 24, 177-183.
 Sun, D.K., et al. (2011) Lattice Boltzmann Modeling of Dendritic Growth in Forced and Natural Convection. Computers & Mathematics with Applications, 61, 3585-3592.
 Eshraghi, M., Jelinek, B. and Felicelli, S.D. (2015) Large-Scale Three-Dimensional Simulation of Dendritic Solidification Using Lattice Boltzmann Method. JOM, 67, 1786-1792.
 Wang, W., Luo, S. and Zhu, M. (2016) Numerical Simulation of Three-Dimensional Dendritic Growth of Alloy: Part II—Model Application to Fe-0.82WtPctC Alloy. Metallurgical and Materials Transactions A, 47, 1355-1366.
 Wang, W., Luo, S. and Zhu, M. (2016) Numerical Simulation of Three-Dimensional Dendritic Growth of Alloy: Part I—Model Development and Test. Metallurgical and Materials Transactions A, 47, 1339-1354.
 Zhang, X., et al. (2012) A Three-Dimensional Cellular Automaton Model for Dendritic Growth in Multi-Component Alloys. Acta Materialia, 60, 2249-2257.
 Zhu, M.F. and Stefanescu, D.M. (2007) Virtual Front Tracking Model for the Quantitative Modeling of Dendritic Growth in Solidification of Alloys. Acta Materialia, 55, 1741-1755.
 Vermolen, F.J. (2006) Zener Solutions for Particle Growth in Multi-Component Alloys. Delft University of Technology, Faculty of Electrical Engineering, Mathematics and Computer Science, Delft Institute of Applied Mathematics.
 Nestler, B. and Choudhury, A. (2011) Phase-Field Modeling of Multi-Component Systems. Current Opinion in Solid State and Materials Science, 15, 93-105.
 Wang, X. and Xu, X. (2001) Molecular Dynamics Simulation of Heat Transfer and Phase Change during Laser Material Interaction. Journal of Heat Transfer, 124, 265-274.
 Wood, W.W. and Jacobson, J.D. (1957) Preliminary Results from a Recalculation of the Monte Carlo Equation of State of Hard Spheres. The Journal of Chemical Physics, 27, 1207-1208.
 Pusztai, L. and Sváb, E. (1993) Structure Study of Ni62Nb38 Metallic Glass Using Reverse Monte Carlo Simulation. Journal of Non-Crystalline Solids, 156, 973-977.
 Hwang, J., et al. (2012) Nanoscale Structure and Structural Relaxation in Zr50Cu45Al5 Bulk Metallic Glass. Physical Review Letters, 108, Article ID: 195505.
 Parr, R.G., Craig, D.P. and Ross, I.G. (1950) Molecular Orbital Calculations of the Lower Excited Electronic Levels of Benzene, Configuration Interaction Included. The Journal of Chemical Physics, 18, 1561-1563.
 Perim, E., et al. (2016) Spectral Descriptors for Bulk Metallic Glasses Based on the Thermodynamics of Competing Crystalline Phases. Nature Communications, 7, Article ID: 12315.
 Zheng, G.-P. (2012) A Density Functional Theory Study on the Deformation Behaviors of Fe-Si-B Metallic Glasses. International Journal of Molecular Sciences, 13, 10401-10409.
 Daw, M.S. and Baskes, M.I. (1984) Embedded-Atom Method: Derivation and Application to Impurities, Surfaces, and Other Defects in Metals. Physical Review B, 29, 6443.