Flow-induced shear stress and deformation of a core–shell-structured microcapsule in a microchannel

A numerical model was developed and validated to investigate the fluid–structure interactions between fully developed pipe flow and core–shell-structured microcapsule in a microchannel. Different flow rates and microcapsule shell thicknesses were considered. A sixth-order rotational symmetric distribution of von Mises stress over the microcapsule shell can be observed on the microcapsule with a thinner shell configuration, especially at higher flow rate conditions. It is also observed that when being carried along in a fully developed pipe flow, the microcapsule with a thinner shell tends to accumulate stress at a higher rate compared to that with a thicker shell. In general, for the same microcapsule configuration, higher flow velocity would induce a higher stress level over the microcapsule shell. The deformation gradient was used to capture the microcapsule's deformation in the present study. The effect of Young's modulus on the microcapsule shell on the microcapsule deformation was investigated as well. Our findings will shed light on the understanding of the stability of core–shell-structured microcapsule when subjected to flow-induced shear stress in a microfluidic system, enabling a more exquisite control over the breakup dynamics of drug-loaded microcapsule for biomedical applications.     

Decay rate for viscoelastic wave equation of variable coefficients with delay and dynamic boundary conditions

In this paper, we consider a viscoelastic wave equation of variable coefficients in the presence of past history with nonlinear damping and delay in the internal feedback and dynamic boundary conditions. Under suitable assumptions, we establish an explicit and general decay rate result without imposing restrictive assumption on the behavior of the relaxation function at infinity by Riemannian geometry method and Lyapunov functional method.     

Exploration and optimization of conditions for quantitative analysis of lignans in Schisandra chinensis by an online supercritical fluid extraction with supercritical fluid chromatography system

An online supercritical fluid extraction with supercritical fluid chromatography system could provide sequential extraction and quantitative analysis of lignans in Schisandra chinensis. Supercritical fluid extraction conditions were optimized at 15 MPa, 50°C, and 4 min with supercritical CO₂ adding 1% methanol; the elution volume and flow rate were set at 6 mL and 2 mL/min to blow extract out of the tank completely. The split‐flow rate was confirmed at 2.5%, which determines injection volume and accuracy of quantitative detection. The factors having negative influences on supercritical fluid chromatography retention in the online system, including sample loading forms and backpressure settings, are discussed in the paper. At last, an extraction‐quantitative method for lignans in Schisandra chinensis was developed, which could be finished within 19.5 min. The total content percentage of four lignans (Schisandrin, Schisandrin A, Schisandrin B and Schisandrol B) in four batches was respectively measured to be 1.42, 1.54, 1.62, and 1.90%.     

Scattering coefficients for a sphere in a visco-acoustic medium for arbitrary partial wave order

Analytical solutions are reported for the scattering coefficients of a solid elastic sphere suspended in a viscous fluid for arbitrary partial wave order. Expressions are derived for incident compressional and shear wave modes, taking into account the viscosity of the surrounding fluid and resultant wave mode conversion. The long compressional wavelength limit is employed to simplify the derivation, whereas no restriction is placed on the shear wavelength in the fluid compared to the particle dimension. The analytical approximations are compared with numerical results obtained from matrix inversion of the boundary equations and agree within the validity domain of the solutions.     

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.     

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.     

Effects of particle size and concentration on bubble coalescence and froth formation in a slurry bubble column

A new approach for simulating the formation of a froth layer in a slurry bubble column is proposed. Froth is considered a separate phase, comprised of a mixture of gas, liquid, and solid. The simulation was carried out using commercial flow simulation software (FIRE v2014) for particle sizes of 60–150 μm at solid concentrations of 0–40 vol%, and superficial gas velocities of 0.02–0.034 m/s in a slurry bubble column with a hydraulic diameter of 0.2 m and height of 1.2 m. Modelling calculations were conducted using a Eulerian–Eulerian multiphase approach with k–ε turbulence. The population balance equations for bubble breakup, bubble coalescence rate, and the interfacial exchange of mass and momentum were included in the computational fluid dynamics code by writing subroutines in Fortran to track the number density of different bubble sizes. Flow structure, radial gas holdup, and Sauter mean bubble diameter distributions at different column heights were predicted in the pulp zone, while froth volume fraction and density were predicted in the froth zone. The model was validated using available experimental data, and the predicted and experimental results showed reasonable agreement. To demonstrate the effect of increasing solid concentration on the coalescence rate, a solid-effect multiplier in the coalescence efficiency equation was used. The solid-effect multiplier decreased with increasing slurry concentration, causing an increase in bubble coalescence efficiency. A slight decrease in the coalescence efficiency was also observed owing to increasing particle size, which led to a decrease in Sauter mean bubble diameter. The froth volume fraction increased with solid concentration. These results provide an improved understanding of the dynamics of slurry bubble reactors in the presence of hydrophilic particles.     

喷流冷却复合陶瓷薄片激光器的热特性仿真分析及实验研究

为分析喷流冷却复合陶瓷薄片激光器的热特性,设计用于冷却复合陶瓷薄片的喷流冷却系统.利用湍流换热理论和计算流体动力学仿真方法建立喷流冷却复合陶瓷薄片激光器的流固耦合热仿真模型,定义评价其冷却能力和冷却均匀性的定量参数.根据该仿真模型得到喷流冷却系统的最优设计参数,并进行实验验证.使用163孔喷板,流量为0.2kg/s,入口温度为20℃,在1200 W泵浦时获得359 W激光输出功率,并测得复合陶瓷薄片上表面的最高温度为92℃.激光输出功率与复合陶瓷薄片上表面温度均与泵浦功率呈近似正线性关系,且温度的实验值与仿真值相符度较高.     

Sonoelectrochemical hydrogenation of safrole: A reactor design,statistical analysis and computational fluid dynamic approach

In this work, ultrasound-assisted electrocatalytic hydrogenation (US-ECHSA) of safrole was carried out in water medium, using sacrificial anode of nickel. The ultrasonic irradiation was carried out at frequency of 20 kHz ± 500 Hz with a titanium cylindrical horn (MS 73 microtip; Ti-6AI-4V alloy; 3.0 mm diameter). The optimal conditions were analyzed by statistical experimental design (fractional factorial). The influence of the sonoelectrochemical reactor design was also investigated by using computational fluid dynamics as simulation tool. Among the five parameters studied: catalyst type, use of β-cyclodextrin as inverse phase transfer catalyst, sonoelectrochemical reactor design, ultrasound mode and the temperature of the solution, only the last three were significant. The hydrogenation product, dihydrosafrole, reached 94% yield, depending on the experimental conditions applied. Data of computational fluid dynamics showed that a wing shape tube added to the sonoelectrochemical reactor can work as a cooling apparatus, during the electrochemical process. The reactional solution temperature diminishes 14 °C when compared to the four-way-type reactor. Cooper cathode, absence of β-cyclodextrin, four-way-type reactor, ultrasound continuous mode (14 W) and absence of temperature control were the most effective reaction parameters for the safrole hydrogenation using US-ECHSA method. The proposed approach represents an important contribution for understanding the hydrodynamic behavior of sonoelectrochemical reactors designs and, consequently, for the reducing of the experimental costs inherent to the sonoelectrochemical process.