A numerical method for fractional integral with applications

IntroductionThefractionalcalculushasalonghistoryandthereareamassofworkstodiscussthefractionalderivativesandfractionalintegralswitharbitrary (realorcomplex)order[1- 3 ].Thefractionalcalculushasawideapplicationbackground ,especiallyinthefieldsofchemistry ,electromagnetics,materialscienceandmechanics.Forexample,Gement[4 ]proposedthefractionalderivativeconstitutivemodelsofaviscoelasticmaterialatfirst.Themodelshavereceivedincreasingattention[5 - 7].Onlyafewparametersarecontainedinthemodelsandthemo…     

Fluid-solid coupling mathematical model of contaminant transport in unsaturated zone and its asymptotical solution

IntroductionThetransportofcontaminantsinunsaturatedzonehascausedmuchattention .Inearly1960s,contaminationproblemsofsoilandgroundwaterhadbeenstudiedathomeandabroad[1].Andinrecentyears ,thetransformationandtransportationofcontaminantshavebeendeeplystudiedinthefieldsofhydrogeology ,petroleumengineering ,environmentalengineeringandsoon[2 ,3].Somecontaminanttransportmodelshavebeenpresentedsofar.Forexample ,Paker[4 ]etal.presentedaconstitutivemodelgoverningparametersofwater,gasandcontaminantswhenth…     

Analytical solution of a simply supported piezoelectric beam subjected to a uniformly distributed loading

IntroductionIntelligentstructureisakindofnewstructuremodelsandhasbeenreceivedmuchattentioninrecentyears.Usedaspiezoelectricsensorsandactuators,thiskindofintelligentstructureshasmanyadvantagessuchasthepromptnessofresponseandtheconvenienceforthesignalcontrol.Soitisfollowedwithinterestbothinthetheoreticalstudyandintheengineeringapplications[1].Forexample ,theresearchoffundamentalsolution[2 ,3]andellipsoidalinclusion[4 ]forthree_dimensionalpiezoelectricmaterialbyuseofFouriertransformation ;thestu…     

Shift of shock position for a class of nonlinear singularly perturbed problems

Thenonlinearsingularlyperturbedproblemisaveryattentiveobjectofstudyintheinternationalacademiccircles[1].Recently ,manyscholarsstudiedagreatdealofwork .Forexample,Boh啨[2 ],ButuzovandSmurov[3],O’Malley ,Jr.[4 ],Butuzov ,NefedovandSchneider[5 ]andKelley[6 ]andsoon .TheauthorconsideredaclassofsingularlyperturbednonlinearboundaryvalueproblemsfortheordinarydifferentialequationinRefs .[7～ 1 0 ] ,reactiondiffusionequationsinRefs.[1 1～ 1 3 ] ,singularlyperturbedproblemswithnonlocalboundarycond…     

On analysis of the steady flow in an irrectangular parallel-plate flow chamber

I.IntroductionForthestudyofrelationshipbetweencell'sgrowthandshearstress,thesubjectthatshouldbeworkedoutatfirstishowtosimulateshearstressenvironnlelltwhichsuitcell'sgrowthandofwhichtheshearstresscanbecalculatedconvelliently.NowthePPFCisoneofthemainapparatusesthatsimulatetheshearstressenvironmelltforcell'sgrowthinvitro.Forexample,tilerectangularPPFCisusedforthestudyoftherelationshipbetweentheendotheliumcells'injuryinartery,orplatelet'sadhesionordeformationalldatherosclerosisl3].Apartfromthe…     

Bifurcation in a nonlinear duffing system with multi-frequency external periodic forces

I.IntroductionMuchworkofDuffingsystemwithsinglefrequencyexternalf'orcchasbeendoned]'121'131.Also,muchworkofquasi-periodicsolutionwithparanletricandexterllalforces11asbeeddolle.Forexample,LuQishao['j'131studiedthelocalbifurcationandhoplybit'urcationofDLllTingsystemwithtwoparametricforces.Yagasakil']'t71'[8]sttldicdtilechaosphcllOlllcll('nof'DLll'Iillgsystelllwithparalnetricandexternalforces.ChenYtlshuillldWtlllgDcsllil"lstudiedtileroutetochaosfrondquasi-periodicsolutiollofDullingsystenlw…     

Hybrid finite analytic solutions of shallow water circulation

IntroductionThepredictionoftidalandstormsurgesisessentialfornavigation ,fishingandprotectionofthecoast.Furtherthedeterminationofhydrodynamiccirculationbecomesnecessaryforthemodelingofdispersionofpollutants,sedimenttransportandsalinitydistribution .Recently ,somekindsofnumericalmodelingtechniqueshavebeenextensivelydevelopedinthelasttwodecades.Inordertoimprovethewaterqualityinabays/lagoons ,manymethodshavebeenproposedbyresearchers.Forexample ,Komatsuetal.(2 0 0 0 ) [1] proposedanewmethodtocontr…     

IMPULSIVE CONTROL FOR THE STABILIZATION AND SYNCHRONIZATION OF LUR'E SYSTEMS

IntroductionInpractice ,thereexistmanyexamplesofimpulsivecontrolsystems (seeRefs.[1 ,2 ,3 ] ) .Recently ,impulsivecontrolhasbeenwidelyusedtostabilizeandsynchronizechaoticsystems(seeRefs.[4～ 1 5] ) .Forexample,thepredictivePoincarecontrol[4 ]andtheoccasionalproportionalfeedbackcontrol[5 ]aretwoimpulsivecontrolschemeswithvaryingimpulsiveintervals.Yangetal.inRefs.[7,1 1 ]studied ,respectively ,thestabilizationandsynchronizationofaclassofchaoticsystemscalledLorenzsystems .YangandChuainRef.[1 …     

On the calculation of energy release rate for viscoelastic cracked laminates

IntroductionGenerallyspeaking ,acrackinviscoelasticmaterialswillgrowinsomeunknownspeedeveniftheappliedloadisquasistatic,thusthemaindifficultyisinducedindervingtheenergyreleaserate .Knowledgeoftheconditiongoverningthedelaminationincompositelaminateswithviscoelasticlayersisofparamountimportanceinpracticalapplications.Forexample,inplasticencapsulatedICpackages,theinterfacialdelaminationbetweentheSilicondieandviscoelasticepoxymoldingcompoundunderthethermalloadingisthemainfailuremodeofthestructure…     

Plastic dynamic stability of a column under nonconservative forces

I.Introduction'Whenadeformableobjectmovinginamedium,thedirectionsofinteractiveforcessuchasfrictionalforceshouldvarywiththedeformationoftheobject.Thisisakindofnonconservativeforces.Forexample,thedirectionofgasfrictionalforceactingontherocketsandthefossiles…     

Quasiperiodic structures in the problem of fluid convection between horizontal planes

A steady flow induced by quasiperiodic temperature distributions at the boundaries is considered with reference to the example of fluid convection in a plane horizontal layer heated from below. A spectral-finite-difference method is employed to obtain convective vortex structures varying quasiperiodically along the horizontal coordinate. The properties of the quasiperiodic solutions, together with their spectra and integral characteristics (Nusselt number) are studied as functions of the Grashof number and the boundary conditions of the problem.     

Reconstruction of a quasi-instantaneous image of coherent structures from hotwire signals obtained by a multi-point simultaneous measurement system

The turbulent velocity components (u, v) at 11 points in a reciprocating oscillatory turbulent flow have been measured simultaneously by a set of eleven X-type hotwire probes located in a plane perpendicular to the mean flow. Using a conditional sampling technique and a new method of data analysis for the inverse estimation of flow fields called the “virtual plate/load and MASCON model”, a quasi-instantaneous three-dimensional image of coherent structures of turbulence was first reconstructed directly from the experimental velocity data. The quasi-instantaneous image was expressed in terms of the velocity components u, v, w and the vorticity components ω_x ω_y, ω_z and we found that the large-scale coherent structure was composed of a pair of counter-rotating fluid motions with asymmetry which was quite different from that of the ensemble-averaged one. Flow patterns induced by the large-scale structure have been clarified by perspective representations visualized by computer simulations that produce timelines and streaklines of fluid particle traces. Results showed that the new experimental method was applicable for investigating the three-dimensional feature of coherent structures including asymmetry.     

Numerical modeling of flow past a volumeless and thin rigid body using direct forcing immersed boundary method

The new capability has been added as the numerical method for modeling volumeless and thin rigid bodies to the direct forcing immersed boundary (DFIB) method. The DFIB approach is based on adding a virtual force to the Navier–Stokes equations of incompressible flow to account for the interaction between the fluid and structures. The volume of a solid function (VOS) identifies the stationary or moving solid structures in a given fluid domain. A new VOS-based algorithm was developed to identify thin, rigid structure boundary points in fluid flow and ensure that the fluid cannot cross through the boundary of a thin rigid structure while moving or stationary. The DFIB method was first validated in a three-dimensional (3D) turbulent flow over a circular cylinder. The large-eddy simulation simulated the turbulent flow scales. The proposed algorithm was tested using a 3D turbulent flow past a stationary and rotating Savonius wind turbine that functions as a thin, rigid body. The validation results showed that the selected DFIB approach, combined with the novel algorithm, could simulate a thin, volumeless, rigid structure that is stationary and rotating in incompressible turbulent flows. The current method is also applicable for two-way fluid-structure interaction problems.     

Fully resolved simulations of viscoelastic suspensions by an efficient immersed boundary-lattice Boltzmann method

An efficient immersed boundary-lattice Boltzmann method (IB-LBM) is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows. Stress LBM based on Giesekus and Oldroyd-B constitutive equation are used to model the viscoelastic stress tensor. A boundary thickening-based direct forcing IB method is adopted to solve the particle–fluid interactions with high accuracy for non-slip boundary conditions. A universal law is proposed to determine the diffusivity constant in a viscoelastic LBM model to balance the numerical accuracy and stability over a wide range of computational parameters. An asynchronous calculation strategy is adopted to further improve the computing efficiency. The method was firstly applicated to the simulation of sedimentation of a single particle and a pair of particles after good validations in cases of the flow past a fixed cylinder and particle migration in a Couette flow against FEM and FVM methods. The determination of the asynchronous calculation strategy and the effect of viscoelastic stress distribution on the settling behaviors of one and two particles are revealed. Subsequently, 504 particles settling in a closed cavity was simulated and the phenomenon that the viscoelastic stress stabilizing the Rayleigh–Taylor instabilities was observed. At last, simulations of a dense flow involving 11001 particles, the largest number of particles to date, were performed to investigate the instability behavior induced by elastic effect under hydrodynamic interactions in a viscoelastic fluid. The elasticity-induced ordering of the particle structures and fluid bubble structures in this dense flow is revealed for the first time. These simulations demonstrate the capability and prospects of the present method for aid in understanding the complex behaviors of viscoelastic particle suspensions.     

Simulation of a Viscous Flow Around an Oscillating Cylinder at Low Keulegan-Carpenter Numbers

A finite difference simulation method for a viscous flow around a circular cylinder sinusoidally oscillating at low Keulegan-Carpenter numbers is presented. Navier-Stokes equations in finite difference form are solved on a moving grid system, based on a time dependent coordinate transformation. Evolution with time of the flow structures induced by a circular cylinder performing sinusoidal oscillations in a fluid at rest, by means of stream lines, pressure contours and vortex shedding is studied in detail at Keulegan-Carpenter numbers, Kc = 9.4 and 14. The time dependent drag and lift are also explained.     

X-ray based measurements of the local solid phase content in a three-phase flow of a bubble column: statistical significance

Information regarding the flow properties in bubble columns with a three-phase flow is of great interest for research into its performance, as well as for the validation of computer models. In an earlier paper, we proposed an X-ray based particle tracking velocimetry (PTV), called XPTV. This method allows not only the measurement of the velocities of both solid and fluid phases in a three-phase flow, but also the assessment of the time-averaged local solid phase content. In this paper, we are concerned with the statistical significance of the obtained measurements.     

Numerical stabilities of loosely coupled methods for robust modeling of lightweight and flexible structures in incompressible and viscous flows

The growing interest to examine the hydroelastic dynamics and stabilities of lightweight and flexible materials requires robust and accurate fluid–structure interaction(FSI)models. Classically, partitioned fluid and structure solvers are easier to implement compared to monolithic methods;however, partitioned FSI models are vulnerable to numerical("virtual added mass") instabilities for cases when the solid to fluid density ratio is low and if the flow is incompressible.As a partitioned method, the loosely hybrid coupled(LHC)method, which was introduced and validated in Young et al.(Acta Mech. Sin. 28:1030–1041, 2012), has been successfully used to efficiently and stably model lightweight and flexible structures. The LHC method achieves its numerical stability by, in addition to the viscous fluid forces, embedding potential flow approximations of the fluid induced forces to transform the partitioned FSI model into a semi-implicit scheme. The objective of this work is to derive and validate the numerical stability boundary of the LHC. The results show that the stability boundary of the LHC is much wider than traditional loosely coupled methods for a variety of numerical integration schemes. The results also show that inclusion of an estimate of the fluid inertial forces is the most critical to ensure the numerical stability when solving for fluid–structure interaction problems involving cases with a solid to fluid-added mass ratio less than one.     

Wavy mode of the streaked flow around an oscillating cylinder in a stratified fluid at rest

The structures of the flow induced by a vertical circular cylinder performing transverse oscillations in a linearly stratified fluid at rest are investigated. The density gradient inhibits the onset of the three-dimensional instability. The instability appears as regularly spaced streaked flows along the cylinder axis. Each streaked flow follows a wavy path in the horizontal plane in a certain range of amplitude and frequency of the oscillations.     

明渠湍流边界层中颗粒的运动与分布

利用直接数值模拟、点球浸入边界法和颗粒离散元法相结合的方法, 模拟了颗粒在明渠湍流边界层中的运动, 并对颗粒的瞬时位置进行了Voronoi 分析, 定量研究了颗粒在湍流边界层中的运动和分布规律. 研究发现:颗粒的输运对湍流的统计特征有影响, 其运动与近壁区湍流拟序结构密切相关, 在"喷发"结构作用下被带离壁面, 在"扫掠" 结构和自身重力作用下回到壁面; 在湍流边界层中, 颗粒倾向于聚集在低流速带, 呈条带状分布;颗粒在大部分时间处于"簇"状态, 偶尔跳跃到"空" 状态, 但能够很快回到邻近低速区域.      

Natural convection induced in a fluid by the hot inner walls of a “submerged” channel

The jet flows induced around a submerged channel due to the hot inner channel walls and the flow inside the channel are calculated. The formation of high-and low-density regions at the inlet and outlet of the channel is detected. The dependence of the flow rate on the channel orientation relative to the gravity force is analyzed. The onset of coherent flow structures results in the development of unsteady oscillating flows. Natural convection in the fluid is studied using the JoinCAD/FEM program package. The regularized Oberbeck-Boussinesq equations are solved using a finite-element method with the same order of the approximating functions.