The improvement of numerical modeling in the solution of incompressible viscous flow problems using finite element method based on spherical Hankel shape functions

In this paper, the finite element method with new spherical Hankel shape functions is developed for simulating 2‐dimensional incompressible viscous fluid problems. In order to approximate the hydrodynamic variables, the finite element method based on new shape functions is reformulated. The governing equations are the Navier‐Stokes equations solved by the finite element method with the classic Lagrange and spherical Hankel shape functions. The new shape functions are derived using the first and second kinds of Bessel functions. In addition, these functions have properties such as piecewise continuity. For the enrichment of Hankel radial basis functions, polynomial terms are added to the functional expansion that only employs spherical Hankel radial basis functions in the approximation. In addition, the participation of spherical Bessel function fields has enhanced the robustness and efficiency of the interpolation. To demonstrate the efficiency and accuracy of these shape functions, 4 benchmark tests in fluid mechanics are considered. Then, the present model results are compared with the classic finite element results and available analytical and numerical solutions. The results show that the proposed method, even with less number of elements, is more accurate than the classic finite element method.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

A remark on similarity analysis of boundary layer equations of a class of non-Newtonian fluids

A similarity analysis of three-dimensional boundary layer equations of a class of non-Newtonian fluid in which the stress, an arbitrary function of rates of strain, is studied. It is shown that under any group of transformation, for an arbitrary stress function, not all non-Newtonian fluids possess a similarity solution for the flow past a wedge inclined at arbitrary angle except Ostwald-de-Waele power-law fluid. Further it is observed, for non-Newtonian fluids of any model only $90°$ of wedge flow leads to similarity solutions. Our results contain a correction to some flaws in Pakdemirli׳s [14] (1994) paper on similarity analysis of boundary layer equations of a class of non-Newtonian fluids.     

A modern approach to the solution and analysis of the Simha-Somcynsky model for the description of pressure-volume-temperature properties of polymer fluids

We revisit the Simha-Somcynsky model of polymer fluids with the purpose of developing novel theoretical and computational approaches to simplify and speed up its solution as well as the fitting of experimental data, and decrease its level of mathematical complexity. We report a novel method that allows us to solve one of the two equations of the model exactly, thus putting the level of mathematical difficulty on a par with the one of other models for polymer fluids. Moreover, we describe a computational algorithm capable of fitting all five parameters of the model in an unbiased way. The results obtained reproduce literature results and fit experimental pressure-volume-temperature and solubility parameter data for three polymers very accurately. Moreover, the new techniques allow for the investigation of the model at very low temperatures. Unexpectedly, the model predicts behaviors that could be interpreted as a glass transition, as routinely observed in dilatometry and differential scanning calorimetry, and a glass phase. We compared the predicted and experimental T g’s for cis poly(1,4-butadiene) and found an excellent quantitative agreement.     

Biocatalytic Feedback‐Controlled Non‐Newtonian Fluids

Non‐Newtonian fluids are ubiquitous in daily life and industrial applications. Herein, we report an intelligent fluidic system integrating two distinct non‐Newtonian rheological properties mediated by an autocatalytic enzyme reaction. Associative polyelectrolytes bearing a small amount of ionic and alkyl groups are engineered: by carefully balancing the charge density and the hydrophobic effect, the polymer solutions demonstrate a unique shear thickening property at low pH while shear thinning at high pH. The urea‐urease clock reaction is utilized to program a feedback‐induced pH change, leading to a strong upturn of the nonlinear viscoelastic properties. As long as the chemical fuel is supplied, two distinct non‐Newtonian states can be achieved with a tunable lifetime span. As a proof of concept, we demonstrate how the physical energy‐driven nonequilibrium properties can be manipulated by a chemical‐fueled process.     

A new study of associating inhomogeneous fluids with classical density functional theory

ABSTRACT

A new density functional for the study of associating inhomogeneous fluids based on Wertheim's first-order thermodynamic perturbation theory is presented and compared to the most currently used associating density functionals. This functional is developed using the weighted density approximation in the range of association of hard spheres. We implement this functional within the framework of classical density functional theory together with modified fundamental measure theory to account for volume exclusion of hard spheres. This approach is tested against molecular simulations from literature of pure associating hard spheres and mixtures of non-associationg and associating hard spheres with different number of bonding sites close to a hard uniform wall. Furthermore, we compare and review our results with the performance of associating functionals from literature, one based on fundamental measure theory and the inhomogeneous version of Wertheim's perturbation theory. Results obtained with classical DFT and the three functionals show excellent agreement with molecular simulations in systems with one hard wall. For the cases of small pores where only one or two layers of fluid are allowed discrepancies between results with classical DFT and molecular simulations were found.     

A modified SPH method to model the coalescence of colliding non-Newtonian liquid droplets

This paper develops a modified smoothed particle hydrodynamics (SPH) method to model the coalescence of colliding non-Newtonian liquid droplets. In the present SPH, a van der Waals (vdW) equation of state is particularly used to represent the gas-to-liquid phase transition similar to that of a real fluid. To remove the unphysical behavior of the particle clustering, also known as tensile instability, an optimized particle shifting technique is implemented in the simulations. To validate the numerical method, the formation of a Newtonian vdW droplet is first tested, and it clearly demonstrates that the tensile instability can be effectively removed. The method is then extended to simulate the head-on binary collision of vdW liquid droplets. Both Newtonian and non-Newtonian fluid flows are considered. The effect of Reynolds number on the coalescence process of droplets is analyzed. It is observed that the time up to the completion of the first oscillation period does not always increase as the Reynolds number increases. Results for the off-center binary collision of non-Newtonian vdW liquid droplets are lastly presented. All the results enrich the simulations of the droplet dynamics and deepen understandings of flow physics. Also, the present SPH is able to model the coalescence of colliding non-Newtonian liquid droplets without tensile instability.     

Crystalline Nitridophosphates by Ammonothermal Synthesis

Nitridophosphates are a well-studied class of compounds with high structural diversity. However, their synthesis is quite challenging, particularly due to the limited thermal stability of starting materials like P₃N₅. Typically, it requires even high-pressure techniques (e.g. multianvil) in most cases. Herein, we establish the ammonothermal method as a versatile synthetic tool to access nitridophosphates with different degrees of condensation. α-Li₁₀P₄N₁₀, β-Li₁₀P₄N₁₀, Li₁₈P₆N₁₆, Ca₂PN₃, SrP₈N₁₄, and LiPN₂ were synthesized in supercritical NH₃ at temperatures and pressures up to 1070 K and 200 MPa employing ammonobasic conditions. The products were analyzed by powder X-ray diffraction, energy dispersive X-ray spectroscopy, and FTIR spectroscopy. Moreover, we established red phosphorus as a starting material for nitridophosphate synthesis instead of commonly used and not readily available precursors, such as P₃N₅. This opens a promising preparative access to the emerging compound class of nitridophosphates.     

Understanding the Structural Differences between Spherical and Rod‐Shaped Human Insulin Nanoparticles Produced by Supercritical Fluids Precipitation

In this study, the thermal denaturation mechanism and secondary structures of two types of human insulin nanoparticles produced by a process of solution‐enhanced dispersion by supercritical fluids using dimethyl sulfoxide (DMSO) and ethanol (EtOH) solutions of insulin are investigated using spectroscopic approaches and molecular dynamics calculations. First, the temperature‐dependent IR spectra of spherical and rod‐shaped insulin nanoparticles prepared from DMSO and EtOH solution, respectively, are analyzed using principal component analysis (PCA) and 2D correlation spectroscopy to obtain a deeper understanding of the molecular structures and thermal behavior of the two insulin particle shapes. All‐atom molecular dynamics (AAMD) calculations are performed to investigate the influence of the solvent molecules on the production of the insulin nanoparticles and to elucidate the geometric differences between the two types of nanoparticles. The results of the PCA, the 2D correlation spectroscopic analysis, and the AAMD calculations clearly reveal that the thermal denaturation mechanisms and the degrees of hydrogen bonding in the spherical and rod‐shaped insulin nanoparticles are different. The polarity of the solvent might not alter the structure or function of the insulin produced, but the solvent polarity does influence the synthesis of different shapes of insulin nanoparticles.