Numerical investigation of high-speed droplet impact using a multiscale two-fluid approach

A single droplet impact onto solid surfaces remains a fundamental and challenging topic in both experimental and numerical studies with significant importance in a plethora of industrial applications, ranging from printing technologies to fuel injection in internal combustion engines. Under high-speed impact conditions, additional complexities arise as a result of the prompt droplet splashing and the subsequent violent fragmentation; thus, different flow regimes and a vast spectrum of sizes for the produced secondary flow structures coexist in the flow field. The present work introduces a numerical methodology to capture the multiscale processes involved with respect to local topological characteristics. The proposed methodology concerns a compressible Σ-Υ two-fluid model with dynamic interface sharpening based on an advanced flow topology detection algorithm. The model has been developed in OpenFOAM® and provides the flexibility of dealing with the multiscale character of droplet splashing, by switching between a sharp and a diffuse interface within the Eulerian-Eulerian framework in segregated and dispersed flow regions, respectively. An additional transport equation for the interface surface area density (Σ) introduces important information for the sub-grid scale phenomena, which is exploited in the dispersed flow regions to provide an insight into the extended cloud of secondary droplets after impact on the target. A high-speed water droplet impact case has been examined and evaluated against new experimental data; these refer to a millimetre size droplet impacting a solid dry smooth surface at velocity as high as 150m/s, which corresponds to a Weber number of ~7.6×10^5. At the investigated impact conditions compressibility effects dominate the early stages of droplet splashing. A strong shock wave forms and propagates inside the droplet, where transonic Mach numbers occur; local Mach numbers up to 2.5 are observed for the expelled surrounding gas outside the droplet. The proposed numerical approach is found to capture relatively accurately the phenomena and provide significant information regarding the produced flow structure dimensions, which is not available from the experiment

Using the generalized interpolation material point method for fluid-solid interactions induced by surface tension

This thesis is devoted to the development of new, Generalized Interpolation Material Point Method (GIMP)-based algorithms for handling surface tension and contact (wetting) in fluid-solid interaction (FSI) problems at small scales. In these problems, surface tension becomes so dominant that its influence on both fluids and solids must be considered. Since analytical solutions for most engineering problems are usually unavailable, numerical methods are needed to describe and predict complicated time-dependent states in the solid and fluid involved due to surface tension effects. Traditional computational methods for handling fluid-solid interactions may not be effective due to their weakness in solving large-deformation problems and the complicated coupling of two different types of computational frameworks: one for solid, and the other for fluid. On the contrary, GIMP, a mesh-free algorithm for solid mechanics problems, is numerically effective in handling problems involving large deformations and fracture. Here we extend the capability of GIMP to handle fluid dynamics problems with surface tension, and to develop a new contact algorithm to deal with the wetting boundary conditions that include the modeling of contact angle and slip near the triple points where the three phases -- fluid, solid, and vapor -- meet. The error of the new GIMP algorithm for FSI problems at small scales, as verified by various benchmark problems, generally falls within the 5% range. In this thesis, we have successfully extended the capability of GIMP for handling FSI problems under surface tension in a one-solver numerical framework, a unique and innovative approach.Chapter 1. Introduction -- Chapter 2. Using the generalized interpolation material point method for fluid dynamics at low reynolds numbers -- Chapter 3. On the modeling of surface tension and its applications by the generalized interpolation material point method -- Chapter 4. Using the generalized interpolation material point method for fluid-solid interactions induced by surface tension -- Chapter 5. Conclusions

Study on SPH Viscosity Term Formulations

For viscosity-dominated flows, the viscous effect plays a much more important role. Since the viscosity term in SPH-governing (Smoothed Particle Hydrodynamics) equations involves the discretization of a second-order derivative, its treatment could be much more challenging than that of a first-order derivative, such as the pressure gradient. The present paper summarizes a series of improved methods for modeling the second-order viscosity force term. By using a benchmark patch test, the numerical accuracy and efficiency of different approaches are evaluated under both uniform and non-uniform particle configurations. Then these viscosity force models are used to compute a documented lid-driven cavity flow and its interaction with a cylinder, from which the most recommended viscosity term formulation has been identified

Towards an efficient and highly accurate coupled numerical modelling approach for wave interactions with a dike on a very shallow foreshore

C

Simulation of Single Vapor Bubble Condensation with Sharp Interface Mass Transfer Model

Pure numerical simulation of phase-change phenomena such as boiling and condensation is challenging, as there is no universal model to calculate the transferred mass in all configurations. Among the existing models, the sharp interface model (Fourier model) seems to be a promising solution. In this study, we investigate the limitation of this model via a comparison of the numerical results with the analytical solution and experimental data. Our study confirms the great importance of the initial thermal boundary layer prescription for a simulation of single bubble condensation. Additionally, we derive a semi-analytical correlation based on energy conservation to estimate the condensing bubble lifetime. This correlation declares that the initial diameter, subcooled temperature, and vapor thermophysical properties determine how long a bubble lasts. The simulations are carried out within the OpenFOAM framework using the VoF method to capture the interface between phases. Our investigation demonstrates that calculation of the curvature of interface with the Contour-Based Reconstruction (CBR) method can suppress the parasitic current up to one order

Interface-Resolving Simulations of Gas-Liquid Two-Phase Flows in Solid Structures of Different Wettability

This PhD study is devoted to numerical investigations of two-phase flows on and through elementary and complex solid structures of varying wettability. The phase-field method is developed and implemented in OpenFOAM®. The numerical method/code is verified by a series of test cases of two-phase flows, and then applied to investigate: (1) droplet wetting on solid surfaces; (2) air bubble rising and interacting with cellular structures and (3) gas-liquid interfacial flows in foam structures

Large-Eddy Simulation of a Classical Hydraulic Jump: Influence of Modelling Parameters on the Predictive Accuracy

Results from large-eddy simulations of a classical hydraulic jump at inlet Froude number two are reported. The computations were performed using the general-purpose finite-volume-based code OpenFOAM\uae, and the primary goal was to evaluate the influence of the modelling parameters on the predictive accuracy, as well as establish the associated best-practice guidelines. A benchmark simulation was conducted on a grid with a 1 mm-cell-edge length to validate the solver and provide a reference solution for the parameter influence study. The remaining simulations covered different selections of the modelling parameters: geometric vs. algebraic interface capturing, three mesh resolution levels, and four choices of the convective flux interpolation scheme. Geometric interface capturing led to better accuracy, but deteriorated the numerical stability and increased the simulation times. Interestingly, numerical dissipation was shown to systematically improve the results, both in terms of accuracy and stability. Strong sensitivity to the grid resolution was observed directly downstream of the toe of the jump

Mathematical Modeling Of Pre And Post Combustion Processes In Coal Power Plant

Coal is a brownish-black sedimentary rock with organic and inorganic constituents. It has been a vital energy resource for humans for millennia. Coal accounts for approximately one quarter of the world’s energy consumption, with 65% of this is energy utilized by residential consumers, and 35% by industrial consumers. Coal operated power stations provide 42% of U.S. electricity supply. The United States hold 96% of coal reserves in North America region, out of which 26% are known for commercial usage. The coal combusted in these power generating facilities requires certain pre-combustion processing, while by-products of coal combustion go through certain post-combustion processing. The application of hydrometallurgical extraction of Rare Earth Elements (REE) from North Dakota Lignite coal feedstock can assist coal value amplification. Extraction of REE from lignite coals liberates REEs and CMs that are vital to electronics, power storage, aviation, and magnets industries. The REE extraction process also reduces the sulfur content of ND lignite coal, along with ash components that foul heat exchange surfaces and can have benefits for post-combustion scrubbing units. When coal is combusted, the exhaust gasses contain carbon dioxide (CO2), sulfur dioxide (SO2), oxides of nitrogen (NOx), water (H2O) and nitrogen (N2). Carbon dioxide comprises approximately 8-10 vol% of the flue gas and is reported to contribute to the greenhouse effect, a primary reason for climate change. Carbon Capture and Storage (CCS) involves of CO2 by use of liquid or solid absorbents to separate CO2 from combustion flue gas. Little data is available on gas-liquid interfacial area correlations in the literature for use of second generation solvents, such as MonoEthanolAmine (MEA), in structured packing absorber columns consisting of thin corrugated metal plates or gauzes, designed to force fluids on complicated paths. While mathematical model development for existing post-combustion carbon capture (PCCC) technologies, such as carbon capture simulations using computational fluid dynamics (CFD) for prediction of mass transfer coefficients is well developed, models describing the behavior of third generation solvents is lacking. Two main research opportunities exist: (i) due to the complex chemistry of coal, there is a requirement for a modeling tool that can account for the coal composition and complex hydrometallurgical extraction processes to assist in designing and sizing pre-combustion REE extraction plants; and (ii) CFD models are required that can capture the mass transfer coefficients of third generation CO2 solvents using structured packing. Two primary hypotheses have been developed to address the research opportunities: (1.) Process modeling of hydrometallurgical extraction of REE provides some theory-based understanding that is complementary to experimental validation and, with the help of chemical kinetics and percentage carboxylation existing in feedstocks, can forecast the efficiency and leachability of other feedstocks, and (2.) A detailed Volume of Fluid (VOF) simulation of coupled mass and momentum transfer problems in small intricate regions of corrugated structured and packed panels placed at 45° angle can be used to predict mass transfer coefficients for third generation solvents by using open-source numerical C/C++ based framework called Open Fields-Operations-And-Manipulations (OpenFOAM). The hydrometallurgical process modeling is developed using METSIM, a leading hydrometallurgical process modeling software tool. The steady state process model provides an overview of REE production along with equipment inventory sizing. The model also has functions to define percentage of organic carboxylic acid bonds present in coal, since, the prior research has identified that the primary association of REE in lignite coal is as weakly-bonded complexes of carboxyl groups, which are targets of the extraction technology. The CFD modeling work is expected to determine critical mass transfer coefficients for CO2 capture using structured packing columns. Further, the developed CFD model and its validity will be tested against experimental data from various industrial and literature sources