A flexible multiscale algorithm based on an improved smoothed particle hydrodynamics method for complex viscoelastic flows

Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.

NUMERICAL SIMULATION OF MORPHOLOGY OF POLYMER CHAIN COILS IN COMPLEX FLOWS

A new coupled finite element formulation is proposed to calculate a conformation tensor model in two complex flows： a planar contraction flow and a planar flow around a symmetrically placed cylinder. The components of conformation tensor are first computed together with the velocity and pressure to describe the change of morphology of polymer chain coils in flow fields. Macroscopic quantities of viscoelastic flow are then calculated based on the conformation tensor. Comparisons between the numerical simulations and experiments for stress patterns and velocity profiles are carried out to prove the validity of the method.

Hybrid Methods for Computing the Streamfunction and Velocity Potential for Complex Flow Fields over Mesoscale Domains

Three types of previously used numerical methods are revisited for computing the streamfunctionψand velocity potentialχfrom the horizontal velocity v in limited domains.The first type,called the SOR-based method,uses a classical successive over-relaxation(SOR)scheme to computeψ(orχ)first with an arbitrary boundary condition(BC)and thenχ(orψ)with the BC derived from v.The second type,called the spectral method,uses spectral formulations to construct the inner part of(ψ,χ)-the inversion of(vorticity,divergence)with a homogeneous BC,and then the remaining harmonic part of(ψ,χ)with BCs from v.The third type,called the integral method,uses integral formulas to compute the internally induced(ψ,χ)-the inversion of domain-internal(vorticity,divergence)using the free-space Greenꞌs function without BCs and then the remaining harmonicψ(orχ)with BCs from v minus the internally-induced part.Although these methods have previously been successfully applied to flows in large-scale and synoptic-scale domains,their accuracy is compromised when applied to complex flows over mesoscale domains,as shown in this paper.To resolve this problem,two hybrid approaches,the integral-SOR method and the integral-spectral method,are developed by combining the first step of the integral method with the second step adopted from the SOR-based and spectral methods,respectively.Upon testing these methods on real-case complex flows,the integral-SOR method is significantly more accurate than the integral-spectral method,noting that the latter is still generally more accurate than the three previously-used methods.The integral-SOR method is recommended for future applications and diagnostic studies of complex flows.

Experimental Study of Sediment Incipience Under Complex Flows

Sediment incipience under flows passing a backward-facing step was studied. A series of experiments were conducted to measure scouring depth, probability of sediment incipience, and instantaneous flow velocity field downstream of a backward-facing step. Instantaneous flow velocity fields were measured by using Particle Image Velocimetry (PIV), and an image processing method for determining probability of sediment incipience was employed to analyze the experimental data. The experimental results showed that the probability of sediment incipience was the highest near the reattachment point, even though the near-wall instantaneous flow velocity and the Reynolds stress were both much higher further downstream of the backward-facing step. The possible me- chanisms are discussed for the sediment incipience near the reattachment point.

Advances in LES of Two-phase Combustion(Ⅱ) LES of Complex Gas-Particle Flows and Coal Combustion

Large-eddy simulation(LES) is under its rapid development and is recognized as a possible second generation of CFD methods used in engineering.Large-eddy simulation of two-phase flows and combustion is particularly important for engineering applications.Some investigators,including the present authors,give their review on LES of spray combustion in gas-turbine combustors and internal combustion engines.However,up to now only a few papers are related to the state-of-the-art on LES of gas-particle flows and combustion.In this paper a review of the advances in LES of complex gas-particle flows and coal combustion is presented.Different sub-grid scale(SGS) stress models and combustion models are described,some of the main results are summarized,and some research needs are discussed.

Application Progress of Computational Fluid Dynamic Techniques for Complex Viscous Flows in Ship and Ocean Engineering

Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.

Investigation on Liquid Holdup in Vertical Zero Net-Liquid Flow

Zero net-liquid flow (ZNLF) is a special case of upward gas-liquid two-phase flow. It is a phenomenon observed as a gas-liquid mixture flows in a conduit but the net liquid flow rate is zero. Investigation on the liquid holdup of ZNLF is conducted in a vertical ten-meter tube with diameter of 76 mm, both for Newtonian and nonNewtonian fluids. The gas phase is air. The Newtonian fluid is water and the non-Newtonian fluids are water-based guar gel solutions. The correlations developed for predicting liquid holdup on the basis of Lockhart-Martinelli parameter are not suitable to ZNLF. A constitutive correlation for the liquid holdup of vertical ZNLF was put forward by using the mass balance. It is found that the liquid holdup in ZNLF is dependent on both the gas flow rate and the flow distribution coefficient.

DIGITAL ANALYSIS TECHNOLOGY FOR MORPHOLOGY OF POLYMER CHAIN COILS IN FLOW FIELDS

Polymer chain coils with entanglement is a crucial scale of structures in polymer materials since their relaxationtimes are matching practical processing times.Based on the phenomenological model of polymer chain coils and a new finiteelement approach,we have designed a computer software including solver,pre-and post-processing modules,and developeda digital analysis technology for the morphology of polymer chain coils in flow fields(DAMPC).Using this technology wemay simulate the morphology development of chain coils in various flow fields,such as simple shear flow,elongational flow,and any complex flow at transient or steady state.The applications made up to now show that the software predictions arecomparable with experimental results.

Laser Doppler Measurements of Twin Impinging Jets Aligned with a Crossflow

This paper presents a detailed analysis of the complex flow beneath two impinging jets aligned with a low-velocity crossflow which is relevant for the future F-35 VSTOL configuration, and provides a quantitative picture of the main features of interest for impingement type of flows. The experiments were carried out for a Reynolds number based on the jet exit conditions of Rej = 4.3 × 10^4, an impingement height of 20.1 jet diameters and for a velocity ratio between the jet exit and the crossflow VR = V/Uo of 22.5. The rear jet is located at S = 6 D downstream of the first jet. The results show a large penetration of the first （upstream）jet that is deflected by the crossflow and impinges on the ground, giving rise to a ground vortex due to the collision of the radial wall and the crossflow that wraps around the impinging point like a scarf. The rear jet （located downstream） it is not so affected by the crossflow in terms of deflection, but due to the downstream wall jet that flows radially from the impinging point of the first jet it does not reach the ground. The results indicate a new flow pattern not yet reported so far, that for a VSTOL aircraft operating in ground vicinity with front wind or small forward movement may result in enhanced under pressures in the aft part of the aircraft causing a suction down force and a change of the pitching moment towards the ground.

Turbulent Energy Budgets of a Ground Vortex Flow

Turbulent kinetic energy budgets are presented for a highly curved flow generated by the collision of plane wall turbulent jet with a low-velocity boundary layer. The different terms are obtained in the vertical plane of symmetry by quadratic interpolation of the LDV （Laser Doppler Velocimetry） measurements, for a wall jet-to-boundary layer velocity ratio of 2. The results, which have relevance to flows encountered in powered-lift aircraft operating in ground effect, quantify the structure of the complex ground vortex flow. The analysis of turbulent energy equation terms using the measured data revealed that production by normal and shear stresses are both very important to the turbulent structure of the impact zone of the ground vortex. This is an indication that the modeling of turbulence of a ground vortex requires a good representation of the production by normal stresses which is most important in the collision zone.

AN EXPLORATION TO THE MODELLING OF TWO-DIMENSIONAL COMPLEX TURBULENT SHEAR FLOWS

The regions with shear stress and mean velocity gradient of opposite sign often exist in complex turbulent shear flows.In these cases,the eddy viscosity hypothesis breaks down.Hinze regards the,departure from eddy viscosity hypothesis as a result from transportation of mean momentum over distance by the large structures and arrives at a shear stress expression including the second order derivatives of the mean velocity.However,his expression greatly overestimates the shear stress.This implies that the flow particles are unlikely to have enough memory of the mean momentum over distance.By assuming the departure from eddy viscosity hypothesis as a result from transportation of the shear stress contained in smaller eddies over distance by the large structures,the present author has arrived at a new shear stress expression.The shear stress estimated so far is in good agreement with the experiments.

Complex granular flows of sticky-wet material on flip-flow screens:Calibration of discrete element simulations

Discrete element method(DEM)is an effective approach for studying the screening process of flip-flow screens.However,there have been few studies focusing on the thick layer of sticky-wet particles on flip-flow screens.To achieve accurate simulations of the thick layer of sticky-wet particles on a flip-flow screen,firstly,the movement law of particle flow was studied,and a multi-regime combination cali-bration method based on characteristics of particle flow regimes was proposed.Based on the Plackett-Burman experiment,the curse of dimensionality caused by multi-state and multi-contact parameters was overcome.Subsequently,the lifting cylinder,rotating drum,and trampoline tests were carried out to obtain macroscopic reference values under various granular flow regimes.The calibration results were then determined using the response surface method and climbing algorithm.Finally,the calibration results were tested at both macroscopic and mesoscopic scales and compared with a commonly used calibration method.The results demonstrate that the calibration method,which considers the multi-state characteristics,improves simulation accuracy by 2%-10%and reduces the simulation error to less than 10%,thus meeting the requirements for engineering optimization of flip-flow screens.

带内部支撑结构的低压排汽缸气动性能分析

The static pressure recovery of low-pressure exhaust hood is important for the overall effectiveness of steam turbines.The tubular and plate stiffeners inside the exhaust contribute to the structural safety of exhaust,which affect the aerodynamic performance.Given the complicated exhaust model coupled with the last stage of turbine,this paper intends to investigate the aerodynamic performance of exhaust hood with individual stiffeners using highfidelity numerical simulations in order to figure out the corresponding effects.The results show that(1)the types of stiffeners have different effects on the aerodynamic performance;and(2)different installation positions and types of plate stiffeners have different effects on aerodynamic performance.The above investigations highlight the future demand regarding reasonable layout and quantity of stiffeners to improve the aerodynamic performance of exhaust as well as maintaining the structural safety.

Molecular Functions of Oxygen-Evolving Complex Family Proteins in Photosynthetic Electron Flow

Oxygen-evolving complex （OEC） protein is the original name for membrane-peripheral subunits of photosystem （PS） II. Recently, multiple isoforms and homologs for OEC proteins have been iden- tified in the chloroplast thylakoid lumen, indicating that functional diversification has occurred in the OEC family. Gene expression profiles suggest that the Arabidopsis OEC proteins are roughly categorized into three groups： the authentic OEC group, the stressresponsive group, and the group including proteins related to the chloroplast NAD（P）H dehydrogenase （NDH） complex involved in cyclic electron transport around PSI. Based on the above gene expression profiles, molecular functions of the OEC family proteins are discussed together with our current knowledge about their functions.

SIMULATIONS OF FLOW INDUCED CORROSION IN API DRILLPIPE CONNECTOR

Drillpipe failure is an outstanding issue in drilling engineering, often involving great financial losses. In view of the special features of the flow channel in the high failure zone, this article analyzes the drillpipe failure mechanism from the point of view of flow induced corrosion. Based on the Eulerian-Langrangian method and the discrete phase model, a numerical simulation method is used to investigate the flows of the drilling fluid in the drillpipe connector during the operation of three typical drilling methods （mud drilling, air drilling and foam drilling）. From the flow field in the drillpipe connector, especially, the velocity and pressure distributions in the threaded nipple and the thickened intermediate belt, one may detect the existence of the flow induced corrosion. Then, some structural optimization measures for the drillpipe connector are proposed, and the optimization effects are compared.

ENERGY FLOWS IN COMPLEX ECOLOGICAL SYSTEMS:A REVIEW

Energy flow drives the complex systems to evolve.The allometric scaling as the universalenergy flow pattern has been found in different scales of ecological systems.It reflects the generalpower law relationship between flow and store.The underlying mechanisms of energy flow patterns areexplained as the branching transportation networks which can be regarded as the result of systematicoptimization of a biological target under constraints.Energy flows in the ecological system may bemodelled by the food web model and population dynamics on the network.This paper reviews thelatest progress on the energy flow patterns,explanatory models for the allometric scaling and modellingapproach of flow and network evolution dynamics in ecology.Furthermore,the possibility of generalizingthese flow patterns,modelling approaches to other complex systems is discussed.

A novel multi-fidelity coupled simulation method for flow systems

For the numerical simulation of flow systems with various complex components, the traditional one-dimensional （1D） network method has its comparative advantage in time consuming and the CFD method has its absolute advantage in the detailed flow capturing. The proper coupling of the advantages of different dimensional methods can strike balance well between time cost and accuracy and then significantly decrease the whole design cycle for the flow systems in modern machines. A novel multi-fidelity coupled simulation method with numerical zooming is developed for flow systems. This method focuses on the integration of one-, two-and three-dimensional codes for various components. Coupled iterative process for the different dimensional simulation cycles of sub-systems is performed until the concerned flow variables of the whole system achieve convergence. Numerical zooming is employed to update boundary data of components with different dimen-sionalities. Based on this method, a highly automatic, multi-discipline computing environment with integrated zooming is developed. The numerical results of Y-Junction and the air system of a jet engine are presented to verify the solution method. They indicate that this type of multi-fidelity simulationmethod can greatly improve the prediction capability for the flow systems.

A hybrid multidimensional Riemann solver to couple self-similar method with MULTV method for complex flows

Since proposed,the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows.However,it needs numerous logical operations in supersonic cases,which limit the method’s applicability in engineering problems greatly.In order to overcome this defect,a hybrid multidimensional Riemann solver,called HMTHS(Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures),is proposed.It simulates the strongly interacting zone by adopting the MHLLES(Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures)scheme at subsonic speeds,which is with a high resolution by considering the second moment in the similarity variables.Also,it adopts the MULTV(Multidimensional Toro and Vasquez)scheme,which is with a high resolution in capturing discontinuities,to simulate the flux at supersonic speeds.Systematic numerical experiments,including both one-dimensional cases and twodimensional cases,are conducted.One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves,shock waves,and linear contact discontinuities.Two-dimensional cases,such as the double Mach reflection case,the supersonic shock/boundary layer interaction case,the hypersonic flow over the cylinder case,and the hypersonic viscous flow over the double-ellipsoid case,indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows.Therefore,it is promising to be widely applied in both scholar and engineering areas.

Lattice Boltzmann Simulation for Complex Flow in a Solar Wall

In this letter,we present a lattice Boltzmann simulation for complex flow in a solar wall system which includes porous media flow and heat transfer,specifically for solar energy utilization through an unglazed transpired solar air collector(UTC).Besides the lattice Boltzmann equation(LBE) for time evolution of particle distribution function for fluid field,we introduce an analogy,LBE for time evolution of distribution function for temperature.Both temperature fields of fluid(air) and solid(porous media) are modeled.We study the effects of fan velocity,solar radiation intensity,porosity,etc.on the thermal performance of the UTC.In general,our simulation results are in good agreement with what in literature.With the current system setting,both fan velocity and solar radiation intensity have significant effect on the thermal performance of the UTC.However,it is shown that the porosity has negligible effect on the heat collector indicating the current system setting might not be realistic.Further examinations of thermal performance in different UTC systems are ongoing.The results are expected to present in near future.

BLOCKAGE COEFFICIENT TECHNIQUE FOR SQUARE ENCLOSURE FLOW WITH NO THICKNESS SPACERS

In this paper,the face blockage coefficient technique and finite volume approach are used to solve numerically the complete N. -S. (Navier-Stokes) equations in a three-dimensional square enclosure with no thickness spacers. The pressure-velocity correction algorithm (SIM PLE -C, Semi-Implicit Method for Pressure Linked Equation Consistent) is adopted in a direct calculation of all physical variables. The computations are performed in three grids with different density. The applicability, feasibility and efficiency of the face blockage coefficient technique applied to simulate three-dimensional viscous flows are investigated and discussed. The computational results of the square enclosure with no thickness spacers have been compared with the numerical results of the enclosure with thickness spacers,which were obtained by using both the face and volume blockage coefficients. The present research shows that the blockage coefficient technique is a simple and convenient numerical technique,with encouraging prospects, for computing the velocity,pressure and other physical quantities of the complex flows in a channel with thin internal spacers and barricades.