ADAPTIVE FINITE ELEMENT METHOD FOR HIGH-SPEED FLOW-STRUCTURE INTERACTION

An adaptive finite element method for high-speed flow-structure interaction is pre- sented.The cell-centered finite element method is combined with an adaptive meshing technique to solve the Navier-Stokes equations for high-speed compressible flow behavior.The energy equation and the quasi-static structural equations for aerodynamically heated structures are solved by applying the Galerkin finite element method.The finite element formulation and computational procedure are de- scribed.Interactions between the high-speed flow,structural heat transfer,and deformation are studied by two applications of Mach 10 flow over an inclined plate,and Mach 4 flow in a channel.

Length and stiffness effects of the attached flexible plate on the flow over a traveling wavy foil

Flow over a traveling wavy foil attached with a flexible plate has been numerically investigated using the lattice Boltzmann method combined with the immersed boundary method. The influence of the flexibility and length of the caudal fin on the locomotion of swimming fish through this simplified model, whereas the fish body is modeled by the undulating foil and the caudal fin by the plate passively flapping as a consequence of fluid-structure interaction. It is found that the plate flexibility denoted by the bending stiffness, as well as the length ratio of tail and body, plays an important role in improving thrust generation and propulsive efficiency. It is also revealed that there exists a parameter region of the plate length and stiffness, in which positive propeller efficiency can be achieved. The effect of the passively flapping flexible plate on the pressure field and the vortex production on the wake is further discussed. It is found that when the length ratio of caudal fin and body is greater than 0.2, a reverse von Kármán vortex street occurs when the bending stiffness is about greater than 1.0, and a great thrust is generated as a result of a large pressure difference occurring across the flexible plate. This work provides physical insight into the role of the caudal fin in fish swimming and may inspire the design of robotic fish.

Dynamical separation of rigid bodies in supersonic flow

A computational investigation of the unsteady separation behavior of rigid bodies in Mach-4 flow is carried out. Two rigid bodies, a sphere and a cube, initially stationary, centroid axially aligned, are released and thereafter fly freely according to the aerodynamic forces experienced. During the separation process, the smaller cube can experience different types of movement and our principal interest here is the non-dimensional transverse velocity of it. The separation behavior is investigated for interactions between a sphere and a cube with different mass ratio and a constant initial distance between them. The qualitative separation behavior and the final transverse velocity of the small body are found to vary strongly with the mass ratio but less sensitive to the initial distance between the two bodies. At a critical mass ratio for a given distance, the smaller body transit from entrainment within the flow region bounded by the larger body's shock to expulsion and the accumulated transverse velocity of the small body is close to maximum. This phenomenon is the so-called ‘shock-wave surfing' phenomenon noted by Laurence & Deiterding for two spheres at hypersonic Mach numbers. Then we investigate the separation behavior of a sphere interaction with a rotary cube and with a non-rotary cube for a given mass ratio and different distance between them. The rotary is found to increase the likelihood of ‘surfing'. Only at a certain initial distance for a given mass ratio the rotary effect of cube can be neglectable.

Enhanced External Counterpulsation Treatment Inhibitting Advanced Atherosclerotic Plaque Progression by Augmenting the Plaque Wall Stress： An in vivo FSI Study Based on Animal Experiment

Enhance extemal counterpulsation （ECP） procedure has exhbited itself to be an effective therapy for the m anagem entof ischem ic card iovascu lar diseases, However, considering that EECP significantly increases the acute diastolic pressure, whether it will intervene in the chronic progression of advanced plaque causing great concern in clilical applkation, but yet rein ains elusive presently. In the current paper, a flu id-structure interface （FSI） num erical model of artery with p iaque corn ponent w as developed based on in vivo hem odynam ic m easurem entperfotm ed h a porcine model, to caku late the m echanical stresses of the plaque before and during EECP, and h lum to assess the potential effects of long-term EECP treatm ent on plaque progression. The resu Its show that E E C P augm ented the wall shear stress （WSS） and plaque w all stress （PWS） over the card lac cycles, aswell as the spacial oscillatory of W SS （WSSG ）. Durhg EECP treatm ent, the PW S level respectively raised 6.82% and 6.07% in two simulation cases. The currentpilot study suggests that E E C P treatm entre ay p lay a positive effect on inh biting the conthued plaque progression by hcreashg the PW S level over the card iac cycles. Meanwhile, the plaque morphology should be taken into consideration while m aking patient- specific plan for Ion g- term E E C P treatment in clinic.

Effects of material damping on flag flutter

Complex interactions of plates with ambient fluid are common in daily lives,e.g.flags flapping in wind,aerofoils oscillating in flow.Recently,the feasibility to harvest energy using the flutter motion has been demonstrated.The objectives of this study are to systematically explore the effects of the material damping on flag flutter,and then to study the energy interchange between the fluid and the flag.In this study,a two-dimensional model was developed.Three dimensionless parameters govern the system,i.e.the mass ratio between the structure and the fluid,the dimensionless fluid velocity and the dimensionless material damping.Results show that the critical velocity increases with the increase of the material damping.The oscillation frequency of the flag decreases with the increase of the material damping,and the time-averaged energy dissipation rate initially increases and then decreases.The increase of the material damping causes the transition of the system from a higher frequency oscillating state to a lower frequency oscillating state,and from a chaotic state to a periodic state.

Experimental study on the vortex structure and path instability of freely falling annular disks

Bodies freely falling in steady water or air are common scenes encountered in various scientific and engineering fields, including the flapping flight of birds and the reentry of a space shuttle. In this work, the freely falling annular thin disks with small dimensionless moments of inertia f and Reynolds number Re are investigated experimentally in a water tank. We use stereo- scopic vision to record the position and orientation of the disks. The flow structure behind the disks is studied by applying fluorescent dye visualization and PIV method. Varying the geometry dimensionless parameter （the inner to outer diameter ratio η and I＊） of the disks reveals two new falling patterns. When ηcritl=0.6〈η〈ηcrit2=0.8, the disks show a random lateral vibration while falling. For high ηcrit2〉0.8, the circular vortex loops shed frequently from the disk, which causes a lengthways vibration superimposed onto straight vertical motion. We also observe another two falling patterns reported previously： hula-hoop and helical motion. By comparing the wake structure of the two motions, we find that the vortex layer twists more violently in the hula-hoop motion, which is the reason for the different trajectory between them. Further research on flow field reveals that the torque on the disk that causes the vibration is due to the formation, elongation and shedding of the vortex.

A new three-dimensional finite-volume non-hydrostatic shock-capturing model for free surface flow

In this paper a new finite-volume non-hydrostatic and shock-capturing three-dimensional model for the simulation of wave-structure interaction and hydrodynamic phenomena（wave refraction, diffraction, shoaling and breaking） is proposed. The model is based on an integral formulation of the Navier-Stokes equations which are solved on a time dependent coordinate system： a coordinate transformation maps the varying coordinates in the physical domain to a uniform transformed space. The equations of motion are discretized by means of a finite-volume shock-capturing numerical procedure based on high order WENO reconstructions. The solution procedure for the equations of motion uses a third order accurate Runge-Kutta（SSPRK） fractional-step method and applies a pressure corrector formulation in order to obtain a divergence-free velocity field at each stage. The proposed model is validated against several benchmark test cases.