Search NASA⌕ Search

SEARCH · Search NASA

Results for “keyholing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Laser melting modes in metal powder bed fusion additive manufacturing

In laser powder bed fusion additive manufacturing of metals, extreme thermal conditions create many highly dynamic physical phenomena such as vaporization and recoil, Marangoni convection, and protrusion and keyhole instability. Collectively however, the full set of phenomena is too complicated for practical applications and, in reality, the melting modes are used as a guideline for printing. With increasing local material temperature beyond the boiling point, the mode can change from conduction to keyhole. These mode designations ignore laser-matter interaction details but in many cases are adequate to determine the approximate microstructures and hence the properties of the build. To date, no consistent, common, and coherent definitions have been agreed upon because of historic limitations in melt pool and vapor depression morphology measurements. Here, we distinguish process-based definitions of different melting modes from those based on postmortem evidence. The latter are mainly derived from the transverse cross-sections of the fusion zone, whereas the former come directly from time-resolved x-ray imaging of melt pool and vapor depression morphologies. These process-based definitions are more strict and physically sound, and they offer new guidelines for laser additive manufacturing practices and create new research directions. Further, we highlight the significance of the keyhole, which substantially enhances the laser energy absorption by the melt pool. Recent studies strongly suggest that stable-keyhole laser melting enables efficient, sustainable, and robust additive manufacturing. The realization of this scenario demands the development of multiphysics models, signal translations from morphology to other feasible signals, and in-process metrology across platforms and scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Physics Of Variable-Polarity Plasma Arc Welding

Report describes experimental study of some of the physical and chemical effects that occur during variable-polarity plasma arc (VPPA) keyhole welding of 2219 aluminum alloy. Comprised three major programs: (1) determination of effects of chemical additions (i.e., impurities) on structure and shape of bead and keyhole; (2) determination of flow in regions surrounding keyhole; (3) development of analog used easily to study flow in keyhole region.

Walsh, Daniel W.↗

Human-Robot Interaction

Human-robot interaction (HRI) is a discipline investigating the factors affecting the interactions between humans and robots. It is important to evaluate how the design of interfaces affect the human's ability to perform tasks effectively and efficiently when working with a robot. By understanding the effects of interface design on human performance, workload, and situation awareness, interfaces can be developed to appropriately support the human in performing tasks with minimal errors and with appropriate interaction time and effort. Thus, the results of research on human-robot interfaces have direct implications for the design of robotic systems. For efficient and effective remote navigation of a rover, a human operator needs to be aware of the robot's environment. However, during teleoperation, operators may get information about the environment only through a robot's front-mounted camera causing a keyhole effect. The keyhole effect reduces situation awareness which may manifest in navigation issues such as higher number of collisions, missing critical aspects of the environment, or reduced speed. One way to compensate for the keyhole effect and the ambiguities operators experience when they teleoperate a robot is adding multiple cameras and including the robot chassis in the camera view. Augmented reality, such as overlays, can also enhance the way a person sees objects in the environment or in camera views by making them more visible. Scenes can be augmented with integrated telemetry, procedures, or map information. Furthermore, the addition of an exocentric (i.e., third-person) field of view from a camera placed in the robot's environment may provide operators with the additional information needed to gain spatial awareness of the robot. Two research studies investigated possible mitigation approaches to address the keyhole effect: 1) combining the inclusion of the robot chassis in the camera view with augmented reality overlays, and 2) modifying the camera frame of reference. The first study investigated the effects of inclusion and exclusion of the robot chassis along with superimposing a simple arrow overlay onto the video feed of operator task performance during teleoperation of a mobile robot in a driving task. In this study, the front half of the robot chassis was made visible through the use of three cameras, two side-facing and one forward-facing. The purpose of the second study was to compare operator performance when teleoperating a robot from an egocentric-only and combined (egocentric plus exocentric camera) view. Camera view parameters that are found to be beneficial in these laboratory experiments can be implemented on NASA rovers and tested in a real-world driving and navigation scenario on-site at the Johnson Space Center.

Sandor, Aniko↗

Elucidating the effect of preheating temperature on melt pool morphology variation in Inconel 718 laser powder bed fusion via simulation and experiment

In laser powder bed fusion (L-PBF) additive manufacturing, the mechanical performance, microstructure and defects of fabricated parts are closely associated with the melt pool morphology, e.g., its dimension and shape through the building process. Past studies have largely focused on how the process parameters such as laser power and scan speed affect melt pool characteristics. In this study, the melt pool morphology variation as a function of preheating temperature in the conduction, transition, and keyhole regimes and the underlying mechanisms in each regime are investigated through ex-situ sample characterization and computation thermal fluid dynamics (CtFD) simulation. Single tracks with different combinations of laser power and scan speed are deposited on an Inconel 718 bare plate preheated to a temperature range of 100–500 °C in the experiment. Significant changes are observed in melt pool morphology as a function of preheating temperature from optical measurements of melt track cross sections. The depth of melt pool in the three regimes increases monotonically with preheating temperature, e.g., at 500 °C, the experimental melt pool depth is increased by 49% in conduction regime, 34% in transition regime and 33% in keyhole regime, respectively, while the variation of melt pool width in each regime does not all follow an increasing trend but depends on the melt pool regimes. Melt pool width variation in the conduction and transition regimes is found to depend on the enhanced heat conduction directly related to temperature dependent thermal properties. Through validated CtFD simulations, it is found that in the keyhole regime the evaporation mass, recoil pressure, and laser drilling effect is enhanced with higher preheating temperature, which gives rise to a deeper melt pool. The simulations also reveal that preheating temperature significantly elongates the melt track length due to the increased flow rate and strong recoil pressure that accelerates the backward flow.

36 MATERIALS SCIENCE↗

Real-time tracking and analysis of gas bubble dynamics in laser powder bed fusion using in-situ X-ray characterization and machine learning

Porosity defects remain a significant challenge in the laser powder bed fusion (LPBF) process, adversely affecting the mechanical properties and reliability of additively manufactured components. Here, this study investigates the real-time formation and trajectory of gas bubbles during LPBF of Al6061 alloy using advanced in-situ X-ray characterization and machine learning. The unsupervised Gaussian mixture model and particle tracking algorithm developed are able to precisely track and quantify the properties of gas bubbles and keyhole pores. Our analysis identified five distinct types of gas bubble formation and movement patterns, emphasizing the diverse origins and behaviors of these defects. It enables precise quantification of trajectories, velocities, and morphological changes of gas bubbles, offering a granular view of the subsurface dynamics within the melt pool. Additionally, we explored keyhole-induced pore dynamics, revealing the critical role of keyhole oscillation and collapse for the formation of both large and small gas pores. It defines four different regions of gas bubble movement within the melt pool, providing a clearer understanding of how local fluid dynamics affect pore behavior. The results underscore the importance of integrating in-situ experimental observation and automated machine learning to develop a more robust predictive model for defect formation in LPBF.

In-situ X-ray imaging↗

High-speed X-ray study of process dynamics caused by surface features during continuous-wave laser polishing

During high-speed X-ray imaging of laser surface polishing experiments of specimens of 316L stainless steel at Argonne National Lab's Advanced Photon Source, it was discovered that the induced keyhole changes shape and dimensions while crossing an engineered surface feature without altering process parameters. It was observed that the post-surface feature keyhole was deeper than that of the pre surface feature keyhole. Here, this work reports on the first in-situ observation of the effect of localized surface geometry on underlying melt pool behavior. This has implications for defect formation mechanisms during laser melting processes that rely on melt pool geometry.

42 ENGINEERING↗

A Gaussian Process-Based extended Goldak heat source model for finite element simulation of laser powder bed fusion additive manufacturing process

In this study, laser powder bed fusion (L-PBF) additive manufacturing (AM) is a key enabling technology to manufacture highly complex and integrated metallic structures. In L-PBF AM process, the melting of the metal powders and the layers underneath can be governed by either “conduction mode” or “keyhole mode”, with the keyhole mode reportedly leading to porosity and decreased strength and ductility by many studies. In part scale simulations, finite element (FE) model is often used to study the temperature distribution during printing and to predict the residual stress, where a volumetric heat flux with a Gaussian or a double ellipsoidal (Goldak) distribution is often applied as the laser heat source. However, the above heat source models can only capture the melt pool shape in the conduction mode, and fail to capture the transition to keyhole melting mode when the process parameters change. To overcome this inaccuracy, an extended Goldak heat source model is proposed by introducing a laser penetration term as a function of laser parameters obtained from a Gaussian-Process (GP) model. The model is validated by “2D pad” AlSi10Mg L-PBF experiments under a wide range of laser power, scan speed, and laser focus offset, and the results show the model successfully captures the measured melt pool shape in all conditions.

36 MATERIALS SCIENCE↗

Determining processing behaviour of pure Cu in laser powder bed fusion using direct micro-calorimetry

We report copper is challenging to process by laser powder bed fusion (LPBF) given its high reflectivity at common infrared laser diode wavelengths and high thermal conductivity. Successful deposition of copper in a predictable and repeatable fashion relies on understanding the development of the keyhole melting regime, as well as heating, melting, boiling and vapour formation behaviour when interacting with a laser beam within an LPBF environment. In this study, in situ optical absorptivity measurements are used to clarify the complex physics of the laser material interaction. Absorptivity of laser energy is measured using direct micro-calorimetry and compared to melt pool depth in correlation to processing parameters. The measured absorptivity for a 100 μm layer thickness of powder was found to be approximately four times higher than that of the bare polished discs. It was also shown that high laser power above 500 W and scan speed up to 150 mm/s are appropriate for effective melting of the powder layer, with these parameters overcoming the threshold required to achieve keyhole melting. This is explained by multiple reflections withing the powder particles and the lower thermal conductivity of packed powder in comparison to bare discs. Melt pool formation was found to be highly unstable and an explosive behavior was observed when in the keyhole regime, caused by high fluctuations in absorptivity values. This work demonstrates calorimetry can be used to monitor melting behaviour in a real-time fashion during processing for this challenging to process material, thereby avoiding unnecessary parametric optimisation. In addition, the parametric window for optimum processing revealed here can inform future work.

36 MATERIALS SCIENCE↗

Direct mechanistic connection between acoustic signals and melt pool morphology during laser powder bed fusion

Various nondestructive diagnostic techniques have been proposed for in situ process monitoring of laser powder bed fusion (LPBF), including melt pool pyrometry, whole-layer optical imaging, acoustic emission, atomic emission spectroscopy, high speed melt pool imaging, and thermionic emission. Correlations between these in situ monitoring signals and defect formation have been demonstrated with acoustic signals having been shown to predict pore formation with especially high confidence in recent machine learning studies. Here, in this work, time-resolved acoustic data are collected in both the conduction and keyhole welding regimes of LPBF-processed Ti-6Al-4V alloy. A non-dimensionalized Strouhal number analysis, used in whistle aeroacoustics, is applied to demonstrate that the acoustic signals recorded in the keyhole regimes can be directly associated with the vapor depression morphology. This mechanistic understanding developed from whistle aeroacoustics shows that acoustic monitoring during the LPBF process can provide a direct probe into the vapor depression dynamics and defect occurrence, especially in the keyhole regimes relevant to printing and defect formation.

36 MATERIALS SCIENCE↗

At grade optical crossover for monolithic optial circuits

Planar optical circuits may be made to cross through each other, (thus eliminating extra steps required to fabricate elevated, nonintersecting crossovers) by control of the dimensions of the crossing light conductors (10, 12) to be significantly greater than d=0.89.lambda. and the angle of crossing as nearly 90.degree. as conveniently possible. A light trap may be provided just ahead of the intersection to trap any light being reflected in the source conductor at angles greater than about 45.degree.. The light trap may take the form of triangular shaped portions (16a, 16b) on each side of the source conductor with the far side of the triangular portion receiving incident light at an angle so that incident light will be reflected to the other side, or it may take the form of windows (18a, 18b) in place of the triangular portions. Planar optical circuit boards (21-23) may be fabricated and stacked to form a keyboard (20) with intersecting conductors (26-29) and keyholes (0-9) where conductors merge at the broad side of the circuit boards. These keyholes may be prearranged to form an array or matrix of keyholes.

Jamieson, Robert S.↗

Innovative Tools Advance Revolutionary Weld Technique

The iconic, orange external tank of the space shuttle launch system not only contains the fuel used by the shuttle s main engines during liftoff but also comprises the shuttle s backbone, supporting the space shuttle orbiter and solid rocket boosters. Given the tank s structural importance and the extreme forces (7.8 million pounds of thrust load) and temperatures it encounters during launch, the welds used to construct the tank must be highly reliable. Variable polarity plasma arc welding, developed for manufacturing the external tank and later employed for building the International Space Station, was until 1994 the best process for joining the aluminum alloys used during construction. That year, Marshall Space Flight Center engineers began experimenting with a relatively new welding technique called friction stir welding (FSW), developed in 1991 by The Welding Institute, of Cambridge, England. FSW differs from traditional fusion welding in that it is a solid-state welding technique, using frictional heat and motion to join structural components without actually melting any of the material. The weld is created by a shouldered pin tool that is plunged into the seam of the materials to be joined. The tool traverses the line while rotating at high speeds, generating friction that heats and softens but does not melt the metal. (The heat produced approaches about 80 percent of the metal s melting temperature.) The pin tool s rotation crushes and stirs the plasticized metal, extruding it along the seam as the tool moves forward. The material cools and consolidates, resulting in a weld with superior mechanical properties as compared to those weld properties of fusion welds. The innovative FSW technology promises a number of attractive benefits. Because the welded materials are not melted, many of the undesirables associated with fusion welding porosity, cracking, shrinkage, and distortion of the weld are minimized or avoided. The process is more energy efficient, safe (no toxic smoke or shielding gas, liquid metal splatter, arcing, dangerous voltage, or radiation), and environmentally sound (no consumables, fumes, or noise) than fusion welding. Under computer control, an automated FSW machine can create welds with high reproducibility, improving efficiency and overall quality of manufactured materials. The process also allows for welding dissimilar metals as well as those metals considered to be "unweldable" such as the 7xxx series aluminum alloys. Its effectiveness and versatility makes FSW useful for aerospace, rail, automotive, marine, and military applications. A downside to FSW, however, is the keyhole opening left in the weld when the FSW pin tool exits the weld joint. This is a significant problem when using the FSW process to join circumferential structures such as pipes and storage containers. Furthermore, weld joints that taper in material thickness also present problems when using the conventional FSW pin tool, because the threaded pin rotating within the weld joint material is a fixed length. There must be capability for the rotating pin to both increase and decrease in length in real time while welding the tapered material. (Both circumferential and tapered thickness weldments are found in the space shuttle external tank.) Marshall engineers addressed both the keyhole and tapered material thickness problems by developing the auto-adjustable pin tool. This unique piece of equipment automatically withdraws the pin into the tool s shoulder for keyhole closeout. In addition, the auto-adjustable pin tool retracts, or shortens, the rotating pin while welding a weld joint that tapers from one thickness to a thinner thickness. This year, the impact of the Marshall innovation was recognized with an "Excellence in Technology Transfer Award" from the Federal Laboratory Consortium.

Source record↗

Uncovering acoustic signatures of pore formation in laser powder bed fusion

Abstract We present a machine learning workflow to discover signatures in acoustic measurements that can be utilized to create a low-dimensional model to accurately predict the location of keyhole pores formed during additive manufacturing processes. Acoustic measurements were sampled at 100 kHz during single-layer laser powder bed fusion (LPBF) experiments, and spatio-temporal registration of pore locations was obtained from post-build radiography. Power spectral density (PSD) estimates of the acoustic data were then decomposed using non-negative matrix factorization with custom $$\varvec{k}$$ k -means clustering (NMF $$\varvec{k}$$ k ) to learn the underlying spectral patterns associated with pore formation. NMF $$\varvec{k}$$ k returned a library of basis signals and matching coefficients to blindly construct a feature space based on the PSD estimates in an optimized fashion. Moreover, the NMF $$\varvec{k}$$ k decomposition led to the development of computationally inexpensive machine learning models which are capable of quickly and accurately identifying pore formation with classification accuracy of supervised and unsupervised label learning greater than 95% and 90%, respectively. The intrinsic data compression of NMF k , the relatively light computational cost of the machine learning workflow, and the high classification accuracy makes the proposed workflow an attractive candidate for edge computing toward in-situ keyhole pore prediction in LPBF.

36 MATERIALS SCIENCE↗

In-situ sensor monitoring of multi-class gas porosity formation in laser powder bed fusion using convolutional neural network

In-situ monitoring of defect formation remains a significant challenge in the laser powder bed fusion (LPBF) process. Recent advances have enabled real-time defect detection with machine learning and in-situ sensing technologies; however, most studies focus on binary classification of keyhole pores, limiting nuanced multi-class pore differentiation and formation mechanisms. This work introduces a multi-class pore detection framework (no pore, small pores < 15 µm, and large pores > 15 µm) by leveraging photodiode sensor data alongside high-fidelity synchrotron X-ray imaging. The 15 µm threshold is selected to distinguish between two fundamentally different defect mechanisms, following the physical size-mechanism boundary established by prior high-resolution synchrotron X-ray characterization of Al6061 LPBF. Distinguishing these classes is critical because large keyhole pores are structurally detrimental, whereas small gas pores are often benign, requiring different process control strategies. Thermal emission monitoring data collected simultaneously with high-speed X-ray imaging at the Stanford Synchrotron Radiation Lightsource (SSRL), are correlated with subsurface melt pool dynamics to establish ground truth. Continuous Wavelet Transform (CWT) with optimized parameters converts the photodiode time-series signals into time–frequency images, facilitating feature extraction. Convolutional Neural Networks (CNN) are then applied for real-time multi-class pore classification in an average inference time of 1 ms per signal window. It achieves 79% accuracy and an Area Under the Receiver Operating Characteristic curve (AUC ROC) score of 0.89 with five-fold cross-validation. The results demonstrate that coupling CWT-based feature engineering with CNN architecture enables reliable multi-class pore detection in Al6061 builds using affordable in-situ sensors. This approach advances scalable and affordable quality assurance in additive manufacturing by moving beyond binary defect detection toward more nuanced classification of porosity mechanisms with in-situ sensors and machine learning.

Laser powder bed fusion, Multi-class pores, In-sit↗

In-situ Observations of Directed Energy Deposition Additive Manufacturing Using High-Speed X-ray Imaging

In laser-based directed energy deposition (DED) additive manufacturing, interactions among the laser beam, particle flow, and melt pool influence the properties of the solidified final part. Two separate DED systems, one with high powder flow rates to represent industrial-scale DED processing and the other with low powder flow rates for individual particle tracking, were synchronized with the high-speed imaging setup at the Advanced Photon Source in Argonne National Laboratory. In-situ x-ray imaging of the DED process using both systems highlighted the influence of powder flow rates. Increased powder flow rates resulted in less laser absorption into the melt pool, leading to a transition from a keyhole mode to a melt pool without a keyhole but with surface fluctuations due to powder flow. Increased velocities of particles during powder flow resulted in a decrease in particle melting times and a greater propensity for porosity formation. Altogether, better understanding of the interactions that occur during various scales of the DED process will enable flexibility, control, and new materials development in DED-based additive manufacturing.

36 MATERIALS SCIENCE↗

Image-Based Fracture Surface Defect Characterization Methods for Additively Manufactured Ti-6Al-4V Tested in Fatigue

Abstract Fatigue initiation in additively manufactured samples/parts often occurs at processed-induced defects such as lack-of-fusion (LoF), keyhole, or other morphological/microstructural defects that have unique characteristics and measurable qualities. Attempts at identifying and minimizing such defects have utilized optimized processing conditions along with in situ and ex situ characterization that includes metallography and/or X-ray computed tomography (XCT). This paper highlights the benefits of using fracture surface analyses to detect and quantify defects that may not be detected by metallography/XCT due to sectioning and resolution limits. In addition to using manual quantification of fatigue initiating LoF and keyhole defects on fracture surfaces, image-based machine learning using convolutional neural networks such as U-Net were also used to automate the process. Statistical analyses were used to identify the extreme cases of defects that initiated and accelerated fatigue and to model the distribution of defect size and shape characteristics to distinguish the type of defect. Initial results show agreement between trained machine learning models and ground truth data in defect segmentation, and the distributions of defect characteristics are distinguishable to particular process-induced defect types.

Materials Science↗

Microstructural Development in Inconel 718 Nickel-Based Superalloy Additively Manufactured by Laser Powder Bed Fusion

Excellent weldability and high temperature stability make Inconel 718 (IN718) one of the most popular alloys to be produced by additive manufacturing. Here in this study, we investigated the effects of laser powder bed fusion (LPBF) parameters on the microstructure and relative density of IN718. The samples were fabricated with independently varied laser power (125 - 350W), laser scan speed (200 - 2200 mm/s), and laser scan rotation (0° - 90°). Archimedes’ method, optical microscopy, and scanning electron microscopy were employed to assess the influence of LPBF parameters on the relative density and microstructure. Optimal processing windows were identified for a wide range of processing parameters and relative density greater than 99.5% was achieved using volumetric energy density between 50 and 100 J/mm3. Microstructural features including melt pool geometry, lack of fusion defect, keyhole porosity, and sub-grain cellular microstructure were examined and quantified to correlate to LPBF parameters. A simple empirical model was postulated to relate relative sample density and LPBF volumetric energy density. Melt pool dimensions were quantitatively measured and compared to estimations based on Rosenthal solution, which yielded a good agreement with the width, but underestimated the depth, particularly at high energy input, due to lack of consideration for keyhole mode. In addition, the sub-grain cellular-dendritic microstructure in the as-built samples was observed to decrease with increasing laser scan speed. Quantification of the sub-micron cellular-dendritic microstructure yielded estimated cooling rate in the order of 10 5 – 10 7 K/s.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multiscale investigation of thermomechanical and compositional developments in Ni alloy 718 under laser processing

Laser processing has been widely employed in various applications due to its exceptional spatial resolution. However, the rapid temperature gradients generated in localized areas present significant challenges for experimental characterization using conventional instruments. To characterize Ni alloy 718 during laser processing, we employed in-situ synchrotron X-ray diffraction with a high-speed detector, a method particularly well-suited for probing processes with high temporal and spatial resolution. Through a series of in-situ experiments, we investigated the local variations in the evolution of microstructures and thermomechanical behaviors within a keyhole mode melt pool. The in situ macroscopic thermomechanical behaviors were quantified using an empirical model derived from diffraction patterns, with experimental results showing reasonable agreement with finite element analysis. Various laser parameters were tested to assess their influences on the residual strains in the melt pools. The results revealed that the residual strain in the keyhole mode melt pool is relatively insensitive to variations in the parameters and is smaller than that in the melt pool created under conduction mode laser scanning. Additionally, we analyzed the shapes of individual diffraction spots, providing insights into the plastic behaviors and compositional developments in the resolidified alloy. The analysis confirmed that compositional variations in a dendritic microstructure manifest as asymmetric broadening of the diffraction spots.

Ni alloy 718↗

Operando synchrotron X-ray analysis of melt pool dynamics in an Al-Sn immiscible alloy

The melt flow in an Al-50vol.% Sn immiscible alloy, produced by single-track laser melting of Al and Sn elemental powders, was studied in real time. High-speed synchrotron X-ray imaging was used to track Al and Sn liquids’ movements, and also to examine elemental distributions in the laser tracks, complimented by electron microscopy after solidification. Key aspects, including melt pool geometry, keyhole instability, and flow dynamics (flow pattern and velocity), were examined using digital image analysis. Relatively deeper melt pools formed at 400 W and 300 mm/s exhibited greater stability, with smooth surfaces, consistent outward flow, and minor vortices near the keyhole. In contrast, shallower pools produced at higher scanning speeds (>500 mm/s) demonstrated greater instability with increased surface waviness, and stronger velocity fluctuations, leading to numerous micro-vortices and increased Al-Sn heterogeneity. Velocity scale estimations, supported by experimental observations, examined the roles of vapour pressure, Marangoni effect, buoyancy, inertial, and surface tension forces in the flow. The results revealed that vapour pressure and mechanical waves dominated at high scanning speeds (shallow pools), while Marangoni forces were equally significant in deep pools at lower speeds (300 mm/s). Buoyancy was found to have minimal impact in both cases. Furthermore, the interaction between inertial and surface tension forces played a critical role in determining the degree of waviness of the pools’ surfaces. These findings offer valuable insights into melt pool dynamics during laser processing of immiscible alloys and other metallic systems using elemental powders, and provide guidance for developing high-fidelity computational fluid dynamics models.

immiscible alloys↗