Issue Information ‐ TOC

Numerical modeling of multiphase flow generally requires a special procedure at the solid wall in order to be consistent with Young's law for static contact angles. The standard approach in the lattice Boltzmann method, which consists of imposing fictive densities at the solid lattice sites, is shown to be deficient for this task. Indeed, fictive mass transfer along the boundary could happen and potentially spoil the numerical results. In particular, when the contact angle is less than 90 degrees, the deficiencies of the standard model are major. Various videos that demonstrate this behavior are provided (Supporting Information). A new approach is proposed and consists of directly imposing the contact angle at the boundaries in much the same way as Dirichlet boundary conditions are generally imposed. The proposed method is able to retrieve analytical solutions for static contact angles in the case of straight and curved boundaries even when variable density and viscosity ratios between the phases are considered. Although the proposed wetting boundary condition is shown to significantly improve the numerical results for one particular class of lattice Boltzmann model, it is believed that other lattice Boltzmann multiphase schemes could also benefit from the underlying ideas of the proposed method. The proposed algorithm is two‐dimensional, and the D2Q9 lattice is used. Copyright © 2016 John Wiley & Sons, Ltd.     

Effect of Injection Parameters and Addition of Nanoscale Materials on the Shrinkage of Polypropylene Copolymer

The effect of variation of injection conditions and addition of nano-calcium carbonate (CaCO₃), nano-silicon dioxide (SiO₂) and full-vulcanized nano-powdered styrene butadiene rubber (PSBR) on the shrinkage of injection-molded polypropylene-ethylene copolymer (90/10, co-PP) were investigated. The results showed that the shrinkage was different for different locations along the flow path. The shrinkage in the length direction of the injection-molded sample varied with the adjustment of the processing parameters, while the shrinkage in the width and thickness direction was almost unchanged. The addition of nano-CaCO₃ and PSBR decreased the shrinkage of co-PP, while the shrinkage of co-PP/ SiO₂ composite was almost unchanged.     

Optimal flow conditions of a tracheobronchial model to reengineer lung structures

The high demand for lung transplants cannot be matched by an adequate number of lungs from donors. Since fully ex-novo lungs are far from being feasible, tissue engi-neering is actively considering implantation of engineered lungs where the devitalized structure of a donor is used as scaffold to be repopulated by stem cells of the receiv-ing patient. A decellularized donated lung is treated inside a bioreactor where transport through the tracheobronchial tree (TBT) will allow for both deposition of stem cells and nour-ishment for their subsequent growth, thus developing new lung tissue. The key concern is to set optimally the boundary conditions to utilize in the bioreactor. We propose a pre-dictive model of slow liquid ventilation, which combines a one-dimensional (1-D) mathematical model of the TBT and a solute deposition model strongly dependent on fluid velocity across the tree. With it, we were able to track and drive the concentration of a generic solute across the airways, look-ing for its optimal distribution. This was given by properly adjusting the pumps' regime serving the bioreactor. A feed-back system, created by coupling the two models, allowed us to derive the optimal pattern. The TBT model can be easily invertible, thus yielding a straightforward flow/pressure law at the inlet to optimize the efficiency of the bioreactor.     

A filter‐based,mass‐conserving lattice Boltzmann method for immiscible multiphase flows

Some issues of He–Chen–Zhang lattice Boltzmann equation (LBE) method (referred as HCZ model) (J. Comput. Physics 1999; 152 :642–663) for immiscible multiphase flows with large density ratio are assessed in this paper. An extended HCZ model with a filter technique and mass correction procedure is proposed based on HCZ's LBE multiphase model. The original HCZ model is capable of maintaining a thin interface but is prone to generating unphysical oscillations in surface tension and index function at moderate values of density ratio. With a filtering technique, the monotonic variation of the index function across the interface is maintained with larger density ratio. Kim's surface tension formulation for diffuse–interface method (J. Comput. Physics 2005; 204 :784–804) is then used to remove unphysical oscillation in the surface tension. Furthermore, as the density ratio increases, the effect of velocity divergence term neglected in the original HCZ model causes significant unphysical mass sources near the interface. By keeping the velocity divergence term, the unphysical mass sources near the interface can be removed with large density ratio. The long‐time accumulation of the modeling and/or numerical errors in the HCZ model also results in the error of mass conservation of each dispersed phase. A mass correction procedure is devised to improve the performance of the method in this regard. For flows over a stationary and a rising bubble, and capillary waves with density ratio up to 100, the present approach yields solutions with interface thickness of about five to six lattices and no long‐time diffusion, significantly advancing the performance of the LBE method for multiphase flow simulations. Copyright © 2010 John Wiley & Sons, Ltd.     

A thionine-based reversible redox sensor in a sequential injection system

According to the current demands of Green Analytical Chemistry and regarding the need for lower reagent consumption with improved analytical performance, an automatic methodology with a flow-through optosensor incorporating solid-phase spectrophotometric detection was developed. The sensor used in this work was based on the redox state of thionine whose oxidized form is blue and reduced form is colorless with monitoring carried out at 621 nm. This redox indicator was immobilized on gel beads and subsequently packed into a flow-through cell. It was then assembled into a sequential injection system and was shown to be an excellent alternative to monitor enzymatic redox reactions automatically as the redox catalysis is performed by glucose dehydrogenase. This enzyme is a representative dehydrogenase enzyme and uses NAD⁺ as cofactor, promoting the oxidation of glucose to glucono-lactone and reduction of NAD⁺ to NADH. The produced NADH promotes color depletion on the surface of the sensor. The calibration graph for glucose was linear between 5.74 × 10⁻⁴ and 2.00 × 10⁻³ mol L⁻¹ and detection limit was 1.72 × 10⁻⁴ mol L⁻¹. Glucose concentration in different samples including sera, salines, perfusion solutions, powder for preparing oral solutions and solutions for hemodialysis was determined. The method proved to be reproducible with a RSD < 5% for glucose determinations.     

The effect of a market factor on information flow between stocks using the minimal spanning tree

We empirically investigated the effects of market factors on the information flow created from N(N−1)/2 linkage relationships among stocks. We also examined the possibility of employing the minimal spanning tree (MST) method, which is capable of reducing the number of links to N−1. We determined that market factors carry important information value regarding information flow among stocks. Moreover, the information flow among stocks showed time-varying properties according to the changes in market status. In particular, we noted that the information flow increased dramatically during periods of market crises. Finally, we confirmed, via the MST method, that the information flow among stocks could be assessed effectively with the reduced linkage relationships among all links among stocks from the perspective of the overall market.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.     

A study of gridless method with dynamic clouds of points for solving unsteady CFD problems in aerodynamics

When solving unsteady computational fluid dynamics problems in aerodynamics with a gridless method, a cloud of points is usually required to be regenerated due to its accommodation to moving boundaries. In order to handle this problem conveniently, a fast dynamic cloud method based on Delaunay graph mapping strategy is proposed in this paper. A dynamic cloud method makes use of algebraic mapping principles and therefore points can be accurately redistributed in the flow field without any iteration. In this way, the structure of the gridless clouds is not necessarily changed so that the clouds regeneration can be avoided successfully. The spatial derivatives of the mathematical modeling of the flow are directly determined by using weighted least‐squares method in each cloud of points, and then numerical fluxes can be obtained. A dual time‐stepping method is further implemented to advance the two‐dimensional Euler equations in arbitrary Lagarangian–Eulerian formulation in time. Finally, unsteady transonic flows over two different oscillating airfoils are simulated with the above method and results obtained are in good agreement with the experimental data. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical simulation of Rhonegletscher from 1874 to 2100

Due to climatic change, many Alpine glaciers have significantly retreated during the last century. In this study we perform the numerical simulation of the temporal and spatial change of Rhonegletscher, Swiss Alps, from 1874 to 2007, and from 2007 to 2100.     

Transient two‐layer thin‐film flow

The transient two‐layer thin‐film planar flow is investigated theoretically in this study. The interplay among inertia, viscous and surface/interfacial tension is emphasized. It is found that the film and interface profiles, as well as the flow field, are strongly influenced by the viscosity ratio, velocity and film thickness ratios at inception, and the surface‐to‐interfacial tension ratio. The nonlinear stability of the steady state reveals the formation of a solitary wave after flow inception, which propagates in the form of a convective instability, with the steady state recovered only in the tail (upstream) region of the wave. In the presence of surface/interfacial tension, surface modulation appears, which grows in wavelength and amplitude with position. The flow is found to be particularly stable for higher viscosity of the lower film layer. Copyright © 2010 John Wiley & Sons, Ltd.