A High Order Accurate Bound-Preserving Compact Finite Difference Scheme for Two-Dimensional Incompressible Flow

For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.

Numerical Study of a Three-Dimensional Laminar Flow in a Rectangular Channel with a 180-Degree Sharp Turn

This study presents a numerical analysis of three-dimensional steady laminar flow in a rectangular channel with a 180-degree sharp turn. The Navier-Stokes equations are solved by using finite difference method for Re = 900. Three-dimensional streamlines and limiting streamlines on wall surface are used to analyze the three-dimensional flow characteristics. Topological theory is applied to limiting streamlines on inner walls of the channel and two-dimensional streamlines at several cross sections. It is also shown that the flow impinges on the end wall of turn and the secondary flow is induced by the curvature in the sharp turn.

Identifying the enhancement mechanism of Al/MoO_(3) reactive multilayered films on the ignition ability of semiconductor bridge using a one-dimensional gas-solid two-phase flow model

Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.

Three-dimensional simulation of pore scale fluid flow in granular ore media with realistic geometry

The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.

Three-dimensional numerical simulation of the flow past a circular cylinder based on LES method

The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.

One-dimensional large-strain consolidation of soft clay with non-Darcian flow and nonlinear compression and permeability of soil

Geometrical nonlinearity of the soft soil and the deviation of water flow in the soft clay from Darcy's law have been well recognized in practice. However, the theory of consolidation, which can account for both the geometrical nonlinearity and the non-Darcian flow, has not been reported so far. In this contribution, a model for the consolidation of soft clay which can allow for these two factors simultaneously is proposed. Utilizing the finite difference method, the numerical model for this problem is developed. With the numerical model, the effects of the geometrical nonlinearity and the non-Darcian flow on the consolidation of the soft soil are investigated. The results show that when the self-weight stress is calculated by the same method, the rate of the non-Darcian consolidation for the large-strain case is larger than that for the small-strain case, but the difference between them is limited. However, the difference between the consolidation rates caused by the non-Darcian and Darcian flows is significant. Therefore, when the geometrical nonlinearity of the soft clay is considered in calculating the consolidation settlement, due to the complexity of the large-strain assumption, the small-strain assumption can be used to replace it if the self-weight stress for the small-strain assumption is calculated by considering its sedimentation. However, due to the aforementioned large difference between the consolidation rates with consideration of the non-Darcian flow in soft clay or not, it is better to consider the non-Darcian flow law for both the small and large stain assumptions.

A fast multipole boundary element method for three-dimensional potential flow problems

A fast multipole methodology (FMM) is developed as a numerical approach to reduce the computational cost andmemory requirements in solving large-scale problems. It is applied to the boundary element method (BEM) for three-dimensional potential flow problems. The algorithm based on mixed multipole expansion and numerical integration isimplemented in combination with an iterative solver. Numerical examinations, on Dirichlet and Neumann problems,are carried out to demonstrate the capability and accuracy of the present method. It has been shown that the methodhas evident advantages in saving memory and computing time when used to solve huge-scale problems which may beprohibitive for the traditional BEM implementation.

ON MULTI-DIMENSIONAL SONIC-SUBSONIC FLOW

In this paper, a compensated compactness framework is established for sonicsubsonic approximate solutions to the n-dimensional （n ≥ 2） Euler equations for steady irrotational flow that may contain stagnation points. This compactness framework holds provided that the approximate solutions are uniformly bounded and satisfy Hloc^-1（Ω） compactness conditions. As illustration, we show the existence of sonic-subsonic weak solution to n-dimensional （n ≥ 2） Euler equations for steady irrotational flow past obstacles or through an infinitely long nozzle. This is the first result concerning the sonic-subsonic limit for n-dimension （n ≥ 3）.

Analytic homotopy solution of generalized three-dimensional channel flow due to uniform stretching of the plate

In this communication a generalized three- dimensional steady flow of a viscous fluid between two infinite parallel plates is considered. The flow is generated due to uniform stretching of the lower plate in x- and y-directions. It is assumed that the upper plate is uniformly porous and is subjected to constant injection. The governing system is fully coupled and nonlinear in nature. A complete analytic solution which is uniformly valid for all values of the dimensionless parameters β Re and λ is obtained by using a purely analytic technique, namely the homotopy analysis method. Also the effects of the parameters β Re and λ on the velocity field are discussed through graphs.

Four-dimensional flow magnetic resonance imaging incirrhosis

Since its introduction in the 1970’s,magnetic resonance imaging(MRI)has become a standard imaging modality.With its broad and standardized application,it is firmly established in the clinical routine and an essential element in cardiovascular and abdominal imaging.In addition to sonography and computer tomography,MRI is a valuable tool for diagnosing cardiovascular and abdominal diseases,for determining disease severity,and for assessing therapeutic success.MRI techniques have improved over the last few decades,revealing not just morphologic information,but functional information about perfusion,diffusion and hemodynamics as well.Four-dimensional(4D)flow MRI,a time-resolved phase contrast-MRI with three-dimensional(3D)anatomic coverage and velocity encoding along all three flow directions has been used to comprehensively assess complex cardiovascular hemodynamics in multiple regions of the body.The technique enables visualization of 3D blood flow patterns and retrospective quantification of blood flow parameters in a region of interest.Over the last few years,4D flow MRI has been increasingly performed in the abdominal region.By applying different acceleration techniques,taking 4D flow MRI measurements has dropped to a reasonable scanning time of 8 to 12 min.These new developments have encouraged a growing number of patient studies in the literature validating the technique’s potential for enhanced evaluation of blood flow parameters within the liver’s complex vascular system.The purpose of this review article is to broaden our understanding of 4D flow MRI for the assessment of liver hemodynamics by providing insights into acquisition,data analysis,visualization and quantification.Furthermore,in this article we highlight its development,focussing on the clinical application of the technique.

THE OPTIMAL TRUNCATED LOW-DIMENSIONAL DYNAMICAL SYSTEMS BASED ON FLOW DATABASES

A new theory on the construction of optimal truncated Low-Dimensional Dynamical Systems (LDDSs) with different physical meanings has been developed, The physical properties of the optimal bases are reflected in the user-defined optimal conditions, Through the analysis of linear and nonlinear examples, it is shown that the LDDSs constructed by using the Proper Orthogonal Decomposition (POD) method are not the optimum. After comparing the errors of LDDSs based on the new theory POD and Fourier methods, it is concluded that the LDDSs based on the new theory are optimally truncated and catch the desired physical properties of the systems.

THREE-DIMENSIONAL COUPLED IMPELLER-VOLUTE SIMULATION OF FLOW IN CENTRIFUGAL PUMP AND PERFORMANCE PREDICTION

A three-dimensional turbulent flow through an entire centrifugal pump is simulated using k-ε turbulence model modified by rotation and curvature, SIMPLEC method and body-fitted coordinate. The velocity and pressure fields are obtained for the pump under various working conditions, which is used to predict the head and hydraulic efficiency of the pump, and the results correspond well with the measured values. The calculation results indicate that the pressure is higher on the pressure side than that on the suction side of the blade; The relative velocity on the suction side gradually decreases from the impeller inlet to the outlet, while increases on the pressure side, it finally results in the lower relative velocity on the suction side and the higher one on the pressure side at the impeller outlet; The impeller flow field is asymmetric, i.e. the velocity and pressure fields arc totally different among all channels in the impeller; In the volute, the static pressure gradually increases with the flow route, and a large pressure gratitude occurs in the tongue; Secondary flow exists in the rear part of the spiral.

CHAOTIC MOTIONS AND LIMIT CYCLE FLUTTER OF TWO-DIMENSIONAL WING IN SUPERSONIC FLOW

Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contains both structural and aerodynamic nonlinearities. Hopf bifurcation theory is used to analyze the flutter speed of the system. The effects of system parameters on the flutter speed are studied. The 4th order Runge-Kutta method is used to calculate the stable limit cycle responses and chaotic motions of the aeroelastic system. Results show that the number and the stability of equilibrium points of the system vary with the increase of flow speed. Besides the simple limit cycle response of period 1, there are also period-doubling responses and chaotic motions in the flutter system. The route leading to chaos in the aeroelastic model used here is the period-doubling bifurcation. The chaotic motions in the system occur only when the flow speed is higher than the linear divergent speed and the initial condition is very small. Moreover, the flow speed regions in which the system behaves chaos axe very narrow.

A Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows

The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller （HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257-311 （2004））. By analyzing the distribution of the finite-time Lyapunov exponent （FTLE） along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE. This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.

NUMERICAL ANALYSIS OF THREE-DIMENSIONAL FLUID FLOW AND HEAT TRANSFER IN TIG WELD POOL WITH FULL PENETRATION

A model is established to analyze three-dimensional fluid flow and heat transfer in TICweld pools with full penetration.It considers the deformation of the molten pool surfaceat the condition of full penetrated workpieees,takes the are pressure as the drivingforce of the pool surface deformation,and determines the surface configuration of weldpool based on the dynamic balance of arc pressure,pool gravity and surface tension atdeformed weld pool surface. The SIMPLER algorithm is used to calculate the fluid flowfield and temperature distribution in TIG weld pools of stainless steel workpieces.TIGwelding experiments are made to verify the validity of the model.It shows the calculatedresults by the model are in good agreement with experimental measurements. professor,Dept of Welding Engineering,Harbin Institute of Technology,Harbin 150006,China

Numerical simulation of two-dimensional granular shearing flows and the friction force of a moving slab on the granular media

This paper studies some interesting features of two-dimensional granular shearing flow by using molecular dynamic approach for a specific granular system. The obtained results show that the probability distribution function of velocities of particles is Gaussian at the central part, but diverts from Gaussian distribution nearby the wall. The macroscopic stress along the vertical direction has large fluctuation around a constant value, the non-zero average velocity occurs mainly near the moving wall, which forms a shearing zone.. In the shearing movement, the volume of the granular material behaves in a random manner. The equivalent fl＇iction coefficient between moving slab and granular material correlates with the moving speed at low velocity, and approaches constant as the velocity is large enough.

Three-dimensional Simulation of Gas/Solid Flow in Spout-fluid Beds with Kinetic Theory of Granular Flow

A three-dimensional Eulerian multiphase model, with closure law according to the kinetic theory of granular flow, was used to study the gas/solid flow behaviors in spout-fluid beds. The influences of the coefficient of restitution due to non-ideal particle collisions on the simulated results were tested. It is demonstrated that the simulated result is strongly affected by the coefficient of restitution. Comparison of simulations with experiments in a small spout-fluid bed showed that an appropriate coefficient of restitution of 0.93 was necessary to simulate the flow characteristics in an underdesigned large size of spout-fluid bed coal gasifier with diameter of lm and height of 6m. The internal jet and gas/solid flow patterns at different operating conditions were obtained. The simulations show that an optimal gas/solid flow pattern for coal gasification is found when the spouting gas flow rate is equal to the fluidizing gas flow rate and the total of them is two and a half times the minimum fluidizing gas flow rate. Besides, the radial distributions of particle velocity and gas velocity show similar tendencies; the radial distributions of particle phase pressure due to particle collisions and the particle pseudo-temperature corresponding to the macroscopic kinetic energy of the random particle motion also show similar tendencies. These indicate that both gas drag force and particle collisions dominate the movement of particles.

NUMERICAL ANALYSIS ON THREE-DIMENSIONAL FLOW FIELD OF TURBINE IN TORQUE CONVERTER

Three-dimensional flow field of turbine in torque converter is simulated by numerical calculation in order to improve the performance of torque converter. Calculation model of a torque converter is presented based on the mixing-plane technology. In the calculation of flow field, the 3D N-S equations are separated by finite-volume method and solved by semi-implicit method for pressure-linked equations（SIMPLE）. Based on flow field calculation, the flow field of turbine is simulated. The velocity and pressure in the flow field of turbine are analyzed. The external performance of the torque converter is also calculated. Results of flow simulation show that there are secondary flow, off flow and velocity gradient in turbine passage. The validity of numerical simulation is verified by comparing the results of numerical simulation with experiment data.

THE STUDY ON THREE-DIMENSIONAL MATHEMATICAL MODEL OF RIVER BED EROSION FOR WATER-SEDIMENT TWO-PHASE FLOW

Based on the tensor analysis of water-sediment two-phase how, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turbulent how. Furthermore, corresponding boundary conditions are also presented in connection with the composition and movement of non-uniform bed material. The theoretical results are applied to the calculation of the float open caisson in the construction period and good results are obtained.

Effects of Street-Bottom and Building-Roof Heating on Flow in Three-Dimensional Street Canyons

Using a computational fluid dynamics （CFD） model, the effects of street-bottom and building-roof heating on flow in three-dimensional street canyons are investigated. The building and street-canyon aspect ratios are one. In the presence of street-bottom heating, as the street-bottom heating intensity increases, the mean kinetic energy increases in the spanwise street canyon formed by the upwind and downwind buildings but decreases in the lower region of the streamwise street canyon. The increase in momentum due to buoyancy force intensifies mechanically induced flow in the spanwise street canyon. The vorticity in the spanwise street canyon strengthens. The temperature increase is not large because relatively cold above-roof-level air comes into the spanwise street canyon. In the presence of both street-bottom and building-roof heating, the mean kinetic energy rather decreases in the spanwise street canyon. This is caused by the decrease in horizontal flow speed at the roof level, which results in the weakening of the mean flow circulation in the spanwise street canyon. It is found that the vorticity in the spanwise street canyon weakens. The temperature increase is relatively large compared with that in the street-bottom heating case, because relatively warm above-roof-level air comes into the spanwise street canyon.