Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase flow”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Uncertainty Quantification Analysis of Both Experimental and CFD Simulation Data of a Bench-scale Fluidized Bed Gasifier

Adequate assessment of the uncertainties in modeling and simulation is becoming an integral part of the simulation based engineering design. The goal of this study is to demonstrate the application of non-intrusive Bayesian uncertainty quantification (UQ) methodology in multiphase (gas-solid) flows with experimental and simulation data, as part of our research efforts to determine the most suited approach for UQ of a bench scale fluidized bed gasifier. UQ analysis was first performed on the available experimental data. Global sensitivity analysis performed as part of the UQ analysis shows that among the three operating factors, steam to oxygen ratio has the most influence on syngas composition in the bench-scale gasifier experiments. An analysis for forward propagation of uncertainties was performed and results show that an increase in steam to oxygen ratio leads to an increase in H2 mole fraction and a decrease in CO mole fraction. These findings are in agreement with the ANOVA analysis performed in the reference experimental study. Another contribution in addition to the UQ analysis is the optimization-based approach to guide to identify next best set of additional experimental samples, should the possibility arise for additional experiments. Hence, the surrogate models constructed as part of the UQ analysis is employed to improve the information gain and make incremental recommendation, should the possibility to add more experiments arise. In the second step, series of simulations were carried out with the open-source computational fluid dynamics software MFiX to reproduce the experimental conditions, where three operating factors, i.e., coal flow rate, coal particle diameter, and steam-to-oxygen ratio, were systematically varied to understand their effect on the syngas composition. Bayesian UQ analysis was performed on the numerical results. As part of Bayesian UQ analysis, a global sensitivity analysis was performed based on the simulation results, which shows that the predicted syngas composition is strongly affected not only by the steam-to-oxygen ratio (which was observed in experiments as well) but also by variation in the coal flow rate and particle diameter (which was not observed in experiments). The carbon monoxide mole fraction is underpredicted at lower steam-to-oxygen ratios and overpredicted at higher steam-to-oxygen ratios. The opposite trend is observed for the carbon dioxide mole fraction. These discrepancies are attributed to either excessive segregation of the phases that leads to the fuel-rich or -lean regions or alternatively the selection of reaction models, where different reaction models and kinetics can lead to different syngas compositions throughout the gasifier. To improve quality of numerical models used, the effect that uncertainties in reaction models for gasification, char oxidation, carbon monoxide oxidation, and water gas shift will have on the syngas composition at different grid resolution, along with bed temperature were investigated. The global sensitivity analysis showed that among various reaction models employed for water gas shift, gasification, char oxidation, the choice of reaction model for water gas shift has the greatest influence on syngas composition, with gasification reaction model being second. Syngas composition also shows a small sensitivity to temperature of the bed. The hydrodynamic behavior of the bed did not change beyond grid spacing of 18 times the particle diameter. However, the syngas concentration continued to be affected by the grid resolution as low as 9 times the particle diameter. This is due to a better resolution of the phasic interface between the gases and solid that leads to stronger heterogeneous reactions. This report is a compilation of three manuscripts published in peer-reviewed journals for the series of studies mentioned above.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2018

Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation

Tight oil formations such as those in the Bakken petroleum system are known to hold hundreds of billions of barrels of oil in place; however, the primary recovery factor for these plays is typically less than 10%. Tight oil formations, including the Bakken Formation, therefore, may be attractive candidates for enhanced oil recovery (EOR) using CO 2 . Multiphase fluid behavior and flow in fluid-rich shales can vary substantially depending on the size of pore throats, and properties such as fluid viscosity and density are much different in nanoscale pores than in macroscale pores. Thus it is critical to understand the nature and distribution of nano-, micro-, and macroscale pores and fracture networks. To address these issues, the Energy & Environmental Research Center (EERC) has been conducting a research program entitled “Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation.” The objectives of the project are 1) the use of advanced characterization methods to better understand and quantify the petrophysical and geomechanical factors that control CO 2 and oil mobility within tight oil formation samples, 2) the determination of CO 2 permeation and oil extraction rates in tight reservoir rocks and organic-rich shales of the Bakken, and 3) the integration of the laboratory-based CO 2 permeation and oil extraction data and the characterization data into geologic models and dynamic simulations to develop predictions of CO 2 storage resource and EOR in the Bakken tight oil formation. A combination of standard and advanced petrophysical characterization techniques were applied to characterize samples of Bakken Formation tight reservoir rock and shales from multiple wells. Techniques included advanced computer tomography (CT) imaging, scanning electron microscopy (SEM) techniques, whole-core and micro x-ray CT imaging, field emission (FE) SEM, and focused ion beam (FIB) SEM. Selected samples were also analyzed for geomechanical properties. X-ray CT imaging yielded information on the occurrence of fractures, bedding planes, fossils, and bioturbation in core, as well as data on bulk density and photoelectric factor logs, which were used to interpret porosity, organic content, and mineralogy. FESEM was used for characterization of nano- and microscale features, including nanoscale pore visualization and micropore and pore throat mineralogy. FIBSEM yielded micro- to nanoscale visualization of fracture networks, porosity and pore-size distribution, connected versus isolated porosity, and distribution of organics. Results from the characterization activities provide insight on nanoscale fracture properties, pore throat mineralogy and connectivity, rock matrix characteristics, mineralogy, and organic content. Laboratory experiments demonstrated that CO 2 can permeate the tight matrix of Bakken shale and nonshale reservoir samples and mobilize oil from those samples. Geologic models were created at scales ranging from the core plug to the reservoir, and dynamic simulations were conducted. The data from the characterization and laboratory-based activities were integrated into modeling research activities to determine the fundamental mechanisms controlling fluid transport in the Bakken, which support EOR scheme design and estimation of CO 2 storage potential in tight oil formations. Simulation results suggest a CO 2 storage resource estimate range of 169 million to 1.5 billion tonnes for the Bakken in North Dakota, possibly resulting in 1.8 billion to 16 billion barrels of incremental oil.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

A CONDUCTION-BASED HEAT PIPE MODEL FOR ANALYZING THE ENTIRE PROCESS OF LIQUID-METAL HEAT PIPE STARTUP

The Heat Pipe-cooled Microreactor (HPM) is one of the micro nuclear reactor designs under active study at the U.S. Idaho National Laboratory. Among the major concerns of HPM research is to understand the startup behavior of the heat pipe-cooled system associated with the startup of the liquid-metal heat pipes initialing from frozen state. The startup of liquid-metal heat pipes typically involves a number of nonlinear mass and heat transport processes including the phase change from solid to liquid and vapor. Hence, it is still a huge challenge to simulate the liquid-metal heat pipe startup using conventional CFD methods and software. The major difficulties of numerical CFD modeling come from the phase-change process, multiphase interaction, microporous wick flow, and compressible gas dynamics that occur during startup of the liquid-metal heat pipes. This paper proposes a simplified conduction-based method to provide practical insights into the entire startup process of the liquid-metal heat pipes while mitigating the challenges of addressing all the complex physics. We discuss the theoretical basis and modeling assumptions to analyze the liquid-metal heat pipe startup from frozen state based solely on heat-conduction equations. Then, the proposed model is implemented into the commercial CFD software to verify the model performance. The model prediction results are discussed via the comparison with the experimental data obtained from sodium heat-pipe startup experiments.

42 ENGINEERING↗
DOE OSTI2018

Multiphase Large Eddy Simulation of a Pulsed Sieve-Plate Extraction Column - 18125

Presented here is research resulting from a framework of single-phase and multiphase modelling trials. These are intended to accurately simulate the internal flow phenomena arising in pulsed sieve-plate extraction columns representative of those found in industry. Three-dimensional unsteady flow calculations using large eddy simulation coupled with dynamic sub-grid scale modelling have been used to capture the transient effects of sinusoidal flow conditions using the open source CFD code OpenFOAM{sup R}. Multiphase interface interactions are modelled using the volume of fluid method with appropriate heavy and light phase fluidic properties capturing surface tension, density and viscous effects. The system studied consists of a 150 mm diameter column with two perforated plates with a fractional free area of 25 %, operational under 1 Hz pulsing frequency at a 50 mm amplitude and plate spacing of 300 mm. The materials simulated include water (single-phase) and 3 M nitric acid and dodecane / tributyl phosphate mixture at 30 vol% (multiphase). Results are compared with appropriate correlations from literature. The dispersive mixing efficiency of the system was also evaluated via a mixing index value which allows quantification of the ratio of rotational and irrotational flow components. This study has found that the predictive methods employed are effective at capturing the different scales of turbulence present within pulsed sieve-plate extraction columns, and its influence on their operation. As such, the work highlights the ability of large eddy simulation-based modelling to predict the complex flows in such column geometries. Furthermore, the techniques employed provide valuable insight into the hydrodynamic flow characteristics and mixing mechanisms arising within pulsed sieve-plate extraction columns. (authors)

42 ENGINEERING↗
DOE OSTI2015

On the use of Darcy's law and invasion-percolation approaches for modeling large-scale geologic carbon sequestration

Most large-scale flow and transport simulations for geologic carbon sequestration (GCS) applications are carried out using simulators that solve flow equations arising from Darcy's law. Recently, the computational advantages of invasion-percolation (IP) modeling approaches have been presented. We show that both the Darcy's-law- and the gravity-capillary balance solved by IP approaches can be derived from the same multiphase continuum momentum equation. More specifically, Darcy's law arises from assuming creeping flow with no viscous momentum transfer to stationary solid grains, while it is assumed in the IP approach that gravity and capillarity are the dominant driving forces in a quasi-static two-phase (or more) system. There is a long history of use of Darcy's law for large-scale GCS simulation. However, simulations based on Darcy's law commonly include significant numerical dispersion as users employ large grid blocks to keep run times practical. In contrast, the computational simplicity of IP approaches allows large-scale models to honor fine-scale hydrostratigraphic details of the storage formation which makes these IP models suitable for analyzing the impact of small-scale heterogeneities on flow. However, the lack of time-dependence in the IP models is a significant disadvantage, while the ability of Darcy's law to simulate a range of flows from single-phase- and pressure-gradient-driven flows to buoyant multiphase gravity-capillary flow is a significant advantage. We believe on balance that Darcy's law simulations should be the preferred approach to large-scale GCS simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

An extended trajectory-mechanics approach for calculating two-phase flow paths

A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

The Plant Water Management Experiments: Soil

A simple means of watering plants in the low-g environment aboard orbiting spacecraft is not obvious. Since the beginning of spaceflight, numerous approaches have been pursued to water plants that seek to maximize plant viability and system reliability, while minimizing crew time and system complexity. We are not there yet. The Plant Water Management (PWM) Soil experiments seek to apply recent advances in low-g capillary fluidics phenomena to the challenges faced by plant growth operations aboard spacecraft. The primary challenge is to establish earth-like flows minimizing low-g specific adaptations required of the plants. This is difficult due to the ever-present fluid physics challenges of poorly-wetting multiphase inertial-visco-capillary flows in geometrically complex conduits and containers. In this paper, we present recent flight results for the PWM Soil experiments where arcillite ‘soil reservoirs’ are arranged in a non-wetting host soil that serves as an O2-breathing wetting barrier. In this way, a largely terrestrial water-soil environment is mimicked where, as liquid is evapo-transpired through the growing plant foliage, the effective water table passively ‘falls’ reducing viscous lengths and increasing water uptake for the plant. We present data from 6 days of 24-7 experiments on the ISS testing 3 different plant root models. We also present and correlate a capillary flow model which captures the primary features of the flow. Our summary is valued for the assessment of current and future low-g plant watering systems employing soil media.

microgravity↗
DOE OSTI2017

The Fascinating and Complex Dynamics of Geyser Eruptions

Geysers episodically erupt liquid and vapor. Despite two centuries of scientific study, basic questions persist-why do geysers exist? What determines eruption intervals, durations, and heights? What initiates eruptions? Through monitoring eruption intervals, analyzing geophysical data, taking measurements within geyser conduits, performing numerical simulations, and constructing laboratory models, some of these questions have been addressed. Geysers are uncommon because they require a combination of abundant water recharge, magmatism, and rhyolite flows to supply heat and silica, and large fractures and cavities overlain by low-permeability materials to trap rising multiphase and multicomponent fluids. Eruptions are driven by the conversion of thermal to kinetic energy during decompression. Larger and deeper cavities permit larger eruptions and promote regularity by isolating water from weather variations. Furthermore, the ejection velocity may be limited by the speed of sound of the liquid + vapor mixture.

58 GEOSCIENCES↗
DOE OSTI2017

Development of a Multiphase Adjoint Capability in OpenFOAM

Thermal hydraulic predictions of flow regimes within nuclear reactor cores require significant computational resources and accuracy to ensure that the core design does not violate the thermal limits of materials. Multiphase computational fluid dynamics (CFD) has become a valuable computational capability for predicting flow characteristics as they relate to applications such as nuclear reactor safety and mixing phenomena around spacer grids within fuel assemblies. Multiphase CFD requires closure relationships for both the effects of turbulence as well as the transfer of mass, momentum, and energy between phases. Understanding the sensitivity of the parameters used in closure relationships is important for improving accuracy, but performing rigorous sensitivity analysis is computationally expensive. Adjoint mathematics have been utilized for reactor physics and fluid dynamics applications with respect to verification, sensitivity analysis, adaptive mesh refinement, and optimization. The usage of adjoint methods to predict the sensitivities of a variety of closure relationship parameters within a given problem has been shown to return sensitivity responses with significantly fewer computations than normal 'large-scale' methods. As the complexity of multiphase CFD methods increases, the demand for accurate verification tools and improvements in computational efficiency also increases, making the development of adjoint capabilities an attractive option. Closure relationships that describe interphase momentum physics in multiphase CFD modelling rely heavily on prescribed constants determined from experimental data. Due to the computational efficiency of an adjoint approach, the time to derive these closure relationships could be dramatically reduced, and data assimilation techniques could be used to return probability density functions of closure models derived using experimental data. This work develops a multiphase adjoint capability for the Eulerian-Eulerian two field, two phase equations. Specifically, the adjoint capability has been developed for the sub-cooled nucleate boiling solver boilEulerFoam developed by Dr. Nam Dinh of NC State University within the open-source CFD library OpenFOAM.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1992

Numerical studies of dilute and dense spray characteristics

Several issues involving the improvement of physical submodels and the computational efficiency in modeling dilute and dense spray combustion are discussed. First, the implementations of a dispersion width approach accounting for turbulent dispersion within each computational parcel is discussed. This is essentially a statistical transport model and the testings of this model confirm the capability of accurately representing dispersion in nearly-homogeneous and inhomogeneous turbulent flows with improved efficiency over the delta function stochastic separated flow model. To account for the dense spray effects, an existing drop collision and coalescence model and a Taylor analogy breakup (TAB) model were employed. These models were incorporated into a state-of-the-art multiphase all-speed transient flow solution procedure. Several examples including nonevaporating, evaporating, and burning dense spray cases were studied. The numerical results show reasonably good comparisons with available experimental data in terms of spray penetration, drop sizes, and overall configuration of a spray flame.

Kim, Y. M.↗
NASA NTRS1992

Studies on fuel spray characteristics in high-pressure environment

The present study deals with several issues involving the improvement of physical submodels and the computational efficiency in modeling dense fuel sprays. To improve the computational efficiency, a parcel PDF approach is implemented which can account for turbulent dispersion within each computational parcel. The advantage of a parcel PDF tracking method is to reduce the number of computational parcels representing the spray dynamics as well as to obtain grid-independent solutions for two-phase flows. To account for the dense spray effects, an existing drop collision and coalescence model, two breakup models, and a Reitz's wave instability model were used. These models were incorporated into a state-of-the-art multiphase all-speed transient flow solution procedure. Comparative performance for two breakup models as well as the turbulence modulation effects are also studied. Validation cases include the nonevaporating and evaporating solid-cone dense sprays. The predictions show a reasonably good agreement with available experimental results in terms of spray penetration, drop sizes, gas and drop mean velocities, and gas and drop rms velocities. The numerical results indicate that the present parcel PDF model has the capability of accurately representing drop dispersion in dense sprays with manageable number of computational parcels.

Shang, H. M.↗
DOE OSTI2019

FEARCE: Fast, Easy, Accurate, and Robust Continuum Engineering Improving fuel efficiency and reducing emissions in combustion engines 2019 R&D 100 Awards

FEARCE investigates and simulates turbulent chemically reactive flows in fuel combustion, such as fuel injection and fuel-air mixing, thermodynamic combustion losses, and combustion–emission formation processes in engines. Our objective is to understand turbulent flows and how they react and interact with other complicating variables, including highly nonlinear heat and mass transfer, small-scale velocity, and chemical kinetics or reactions. Our emphasis is on the fluid dynamics and chemical reactions that take place in all kinds of internal combustion engines, including the reactions between multiple species (such as air, gasoline, carbon dioxide, water, and hundreds of other species) and multiphase (liquid and gas) fluid flow.

42 ENGINEERING↗
DOE OSTI2019

Modeling Three-Dimensional Fluid-Driven Propagation of Multiple Fractures using TOUGH-FEMM

Here, a novel numerical simulation tool named TOUGH-FEMM is presented and applied to model three-dimensional (3D) hydraulic fracturing in porous rock. The fluid flow in both fractures and porous rock is modeled using TOUGH2, which is a well-established code for analysis of multiphase and multi-component fluid flow. Rock deformations associated with fracture propagation are modeled using finite element-meshfree method (FEMM). FEMM is an approach to simulate fracture propagation without remeshing, in which the fracture path does not need to be predetermined. Fracture mechanics with mixed-mode stress intensity factors are employed to detect fracture instability and determine the direction of fracture propagation. TOUGH-FEMM is confirmed for modeling fluid-driven fracture propagation in 3D through a number of simulation examples, including modeling of hydraulic fracturing laboratory experiments and by comparison to independent numerical simulation results for multiple interacting hydraulic fractures at ten to hundred meter scale.

58 GEOSCIENCES↗
DOE OSTI2013

Influence of Bubbles on the Turbulence Anisotropy

Direct numerical simulation (DNS) with interface tracking of turbulent bubbly flows is becoming a major tool in advancing our knowledge in the area of multiphase modeling research. A comprehensive analysis of the turbulent flow structure allows us to evaluate the state-of-the-art computational multiphase fluid dynamics (CMFD) models and to propose new closure laws. The presented research will demonstrate how the multiphase DNS data can inform the development of computational fluid dynamics (CFD) models. In particular, the Reynolds stress distribution will be evaluated for single- and two-phase bubbly flows and the level of turbulence anisotropy will be measured in several scenarios. The results will help determine if the isotropic turbulent models are suitable for bubbly flow applications or if there is a strong need to apply and develop Reynolds-stress turbulent models for two-phase flow CFD modeling.

Engineering↗
DOE OSTI2021

Robust hierarchical dispatch for residential distribution network management considering home thermal flexibility and model predictive control

In the transactive energy (TE) paradigm, the devices of participative consumers, or prosumers, may be aggregated and employed to drive operational objectives at the network level. Home heating, ventilation and air-conditioning (HVAC) systems in particular are well-suited to modulate their behaviours based on both home thermal flexibility and requests from the utility grid. This paper develops a robust, hierarchical power dispatch scheme in the context of a residential distribution network. The formulation couples a unique, multiphase linear distribution optimal power flow (OPF) at the upper level with model predictive control (MPC)-based HVAC fleet controllers at the lower level. The proposed approach is tested on nearly 2000 homes with a three-phase distribution network in an intraday market setting, where two major applications are explored and analysed.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Novel Membrane-Based Corrosion Sensors For Operating In Natural Gas Pipelines

Membrane-based electrochemical sensors (MBES) provide an attractive option for real-time continuous humidity and corrosion rate measurements inside of natural gas pipelines. The MBES, can use electrochemical techniques such as linear polarization resistance and electrochemical impedance spectroscopy even when a conductive aqueous phase is intermittent or unavailable. We now present field-test performance data of our MBES operated inside a multiphase water and natural gas flow facility. The MBES measured internal corrosion and conductivity at the top, bottom, and side of the pipe spool. Tests were carried out with flow rates of 4.6, 9.1, and 18.3 m/s, temperatures of 21, 29, and 35 °C, and a pressure of 1000 psia. All sensors followed similar trends in open circuit potential, corrosion rate, and conductivity over time. Neither flow rate, temperature, nor build-up of metal-based solids on the sensor surfaces impacted the quality of the electrochemical measurements. The field-test results confirm that the membrane-based electrochemical sensors are suitable for monitoring internal corrosion of operational natural-gas pipelines.

Duffy, Timothy↗
DOE OSTI2015

DNS of turbulent flow with hemispherical wall roughness

The present study of the effect of roughness density on the mean flow turbulence parameters is motivated by the need for new generation of boundary conditions for multiphase computational multiphase fluid dynamics (CMFD) models applied to boiling flows. Effect of roughness element density on the turbulent flow in a channel is quantified through direct numerical simulations (DNSs). The Navier--Stokes equations are solved using finite element method and bubbles are approximated as rigid near-hemispherical obstacles at the wall. Six different cases were analysed including channel flow with smooth wall and channel flow with rough wall for five different bubble nucleation site densities. Friction factor and the law of the wall was calculated and compared with the previously published results. Existing correlations for nucleating bubble site density dependency on a wall heat flux were used to obtain a relation between the heat flux and the friction factor, leading to the law of the wall dependency on the heat flux. This separate effect study provides new guidelines on how the heat flux in subcooled boiling regime affects the turbulence behaviour near the wall and guides the computational fluid dynamics model development for boiling two-phase flows.

Mechanics↗
DOE OSTI2018

Entropic multiple-relaxation-time multirange pseudopotential lattice Boltzmann model for two-phase flow

Here, an entropic multiple-relaxation-time lattice Boltzmann approach is coupled to a multirange Shan-Chen pseudopotential model to study the two-phase flow. Compared with previous multiple-relaxation-time multiphase models, this model is stable and accurate for the simulation of a two-phase flow in a much wider range of viscosity and surface tension at a high liquid-vapor density ratio. A stationary droplet surrounded by equilibrium vapor is first simulated to validate this model using the coexistence curve and Laplace’s law. Then, two series of droplet impact behavior, on a liquid film and a flat surface, are simulated in comparison with theoretical or experimental results. Droplet impact on a liquid film is simulated for different Reynolds numbers at high Weber numbers. With the increase of the Sommerfeld parameter, onset of splashing is observed and multiple secondary droplets occur. The droplet spreading ratio agrees well with the square root of time law and is found to be independent of Reynolds number. Moreover, shapes of simulated droplets impacting hydrophilic and superhydrophobic flat surfaces show good agreement with experimental observations through the entire dynamic process. The maximum spreading ratio of a droplet impacting the superhydrophobic flat surface is studied for a large range of Weber numbers. Results show that the rescaled maximum spreading ratios are in good agreement with a universal scaling law. This series of simulations demonstrates that the proposed model accurately captures the complex fluid-fluid and fluid-solid interfacial physical processes for a wide range of Reynolds and Weber numbers at high density ratios.

42 ENGINEERING↗