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.     

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.     

分子动力学模拟研究[Bmim][PF6]+水+醇三元体系的微观结构和分子间相互作用

研究离子液体体系的微观结构和分子间相互作用具有重要意义. 本文对1-丁基-3-甲基咪唑六氟磷酸盐([Bmim][PF₆])+水+乙醇和[Bmim][PF₆]+水+异丙醇三元体系进行了分子模拟研究, 计算了径向分布函数和不同组成的水-醇混合溶剂与离子液体阴阳离子间的相互作用能, 并将其分解为库仑相互作用能和Lennard-Jones(LJ)势能. 在此基础上, 研究了溶液体系的微观结构、分子间相互作用和相行为. 结果表明, 水倾向于与离子液体阴离子和阳离子极性部分作用, 醇倾向于与阴离子和阳离子非极性部分作用; 库仑力主导阴离子-溶剂相互作用, 色散力主导阳离子-溶剂相互作用, 阴阳离子的缔合状态对色散力影响较小, 对库仑力的影响非常显著.