Hybrid algorithm for modeling of fluid-structure interaction in incompressible, viscous flows

The objective of this paper is to present and to validate a new hybrid coupling (HC) algorithm for modeling of fluid-structure interaction (FSI) in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled (LC) and tightly coupled (TC) algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible structures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incompressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen’s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical instability issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.     

3D contaminant migration model with consolidation dependent transport coefficients

Soil consolidation would induce variations of its transport coefficients such as hydraulic conductivity and diffusion coefficient. This paper presents a study of the influence of barrier consolidation on transport coefficients, and a 3D transport model based on mixture theory is proposed for describing the liners that involve circular defects in the geomembrane. The elastoplastic ALPHA model is revised by using the spatially mobilized plane (SMP) criterion for simulating the deformation of the soils. Then, the 3D model coupling the nonlinear consolidation and contaminant advection-diffusion is solved using the finite element software ABAQUS. The results show that the importance of reducing the defect size in the geomembrane and the liner porosity to control the contaminant concentration increase.     

The SPH approach to the process of container filling based on non-linear constitutive models

In this work,the transient free surface of container filling with non-linear constitutive equation’s fluids is numerically investigated by the smoothed particle hydrodynamics(SPH) method.Specifically,the filling process of a square container is considered for non-linear polymer fluids based on the Cross model.The validity of the presented SPH is first verified by solving the Newtonian fluid and OldroydB fluid jet.Various phenomena in the filling process are shown,including the jet buckling,jet thinning,splashing or spluttering,steady filling.Moreover,a new phenomenon of vortex whirling is more evidently observed for the Cross model fluid compared with the Newtonian fluid case.     

Assessment of composite failure and ultimate strength without experiment on composite

This paper attempts to estimate the ultimate strength of a laminated composite only based on its con- stituent properties measured independently. Three important issues involved have been systematically addressed, i.e., stress calculation for the constituent fiber and matrix materials, failure detection for the lamina and laminate upon the internal stresses in their constituents, and input data determination of the constituents from monolithic measurements. There are three important factors to influence the accuracy of the strength prediction. One is the stress concentration factor （SCF） in the matrix. Another is matrix plasticity. The third is thermal residual stresses in the constituents. It is these three factors, however, that have not been sufficiently well realized in the composite community. One can easily find out the elastic and strength parameters of a great many laminae and laminates in the current literature. Unfortunately, necessary information to determine the SCF, the matrix plasticity, and the thermal residual stresses of the composites is rare or incomplete. A useful design methodology is demonstrated in the paper.     

Crack buckling in soft gels under compression

Recent interest in designing soft gels with high fracture toughness has called for simple and robust methods to test fracture behavior. The conventional method of applying tension to a gel sample suffers from a difficulty of sample gripping. In this paper, we study a possible fracture mechanism of soft gels under uni-axial compression. We show that the surfaces of a pre-existing crack, oriented parallel to the loading axis, can buckle at a critical compressive stress. This buckling instability can open the crack surfaces and create highly concentrated stress fields near the crack tip, which can lead to crack growth. We show that the onset of crack buckling can be deduced by a dimensional argument com- bined with an analysis to determine the critical compression needed to induce surface instabilities of an elastic half space. The critical compression for buckling was verified for a neo-Hookean material model using finite element simulations.     

Modeling development on the meso-scale reacting transport phenomena in proton exchange membrane fuel cells

The catalyst layer (CL) of proton exchange mem-brane fuel cell (PEMFC) involves various particles and pores in meso-scale, which has an important effect on the mass, charge (proton and electron) and heat transport coupled with the electrochemical reactions. The coarse-grained molecular dynamics (CG-MD) method is employed as a meso-scale structure reconstruction technique to mimic the self-organization phenomena in the fabrication steps of a CL. The meso-scale structure obtained at the equilibrium state is further analyzed by molecular dynamic (MD) software to provide the necessary microscopic parameters for understanding of multi-scale and-physics processes in CLs. The primary pore size distribution (PSD) and active platinum (Pt) surface areas are also calculated and then compared with the experiments. In addition, we also highlight the implementation method to combine microscopic elementary kinetic reaction schemes with the CG-MD approaches to provide insight into the reactions in CLs. The concepts from CG modeling with particle hydrodynamics (SPH) and the problems on coupling of SPH with finite element modeling (FEM) methods are further outlined and discussed to understand the effects of the meso-scale transport phenomena in fuel cells.     

Experimental study of the boundary layer over an airfoil in plunging motion

This is an experimental study on the boundary layer over an airfoil under steady and unsteady conditions.It specifically deals with the effect of plunging oscillation on the laminar/turbulent characteristics of the boundary layer.The wind tunnel measurements involved surfacemounted hot-film sensors and boundary-layer rake.The experiments were conducted at Reynolds numbers of 0.42×10 6 to 0.84 × 10 6 and the reduced frequency was varied from 0.01 to 0.11.The results of the quasi-wall-shear stress as well as the boundary layer velocity profiles provided important information about the state of the boundary layer over the suction surface of the airfoil in both static and dynamic cases.For the static tests,boundary layer transition occurred through a laminar separation bubble.By increasing the angle of attack,disturbances and the transition location moved toward the leading edge.For the dynamic tests,earlier transition occurred with increasing rather than decreasing effective angle of attack.The mean angle of attack and the oscillating parameters significantly affected the state of the boundary layer.By increasing the reduced frequency,the boundary layer transition was promoted to the upstroke portion of the equivalent angle of attack,but the quasi skin friction coefficient was decreased.     

Effect of induced magnetic field on natural convection in vertical concentric annuli

In the present paper,we have considered the steady fully developed laminar natural convective flow in open ended vertical concentric annuli in the presence of a radial magnetic field.The induced magnetic field produced by the motion of an electrically conducting fluid is taken into account.The transport equations concerned with the considered model are first recast in the non-dimensional form and then unified analytical solutions for the velocity,induced magnetic field and temperature field are obtained for the cases of isothermal and constant heat flux on the inner cylinder of concentric annuli.The effects of the various physical parameters appearing into the model are demonstrated through graphs and tables.It is found that the magnitude of maximum value of the fluid velocity as well as induced magnetic field is greater in the case of isothermal condition compared with the constant heat flux case when the gap between the cylinders is less or equal to 1.70 times the radius of inner cylinder,while reverse trend occurs when the gap between the cylinders is greater than 1.71 times the radius of inner cylinder.These fields are almost the same when the gap between the cylinders is equal to 1.71 times the radius of inner cylinder for both the cases.It is also found that as the Hartmann number increases,there is a flattening tendency for both the velocity and the induced magnetic field.The influence of the induced magnetic field is to increase the velocity profiles.     

Design optimization of the S-frame to improve crashworthiness

In this paper, the S-frames, the front side rail structures of automobile, were investigated for crashworthihess. Various cross-sections including regular polygon, nonconvex polygon and multi-cell with inner stiffener sections were investigated in terms of energy absorption of S-frames. It was determined through extensive numerical simulation that a multi-celI S-frame with double vertical internal stiffeners can absorb more energy than the other configurations. Shape optimization was also carried out to improve energy absorption of the S-frame with a rectangular section. The center composite design of experiment and the sequential response surface method （SRSM） were adopted to construct the approximate design sub-problem, which was then solved by the feasible direction method. An innovative double S- frame was obtained from the optimal result. The optimum configuration of the S-frame was crushed numerically and more plastic hinges as well as shear zones were observed during the crush process. The energy absorption efficiency of the structure with the optimal configuration was improved compared to the initial configuration.     

Study of modeling unsteady blade row interaction in a transonic compressor stage part 1: code development and deterministic correlation analysis

The average-passage equation system (APES) provides a rigorous mathematical framework for accounting for the unsteady blade row interaction through multistage compressors in steady state environment by introducing deterministic correlations (DC) that need to be modeled to close the equation system.The primary purpose of this study is to provide insight into the DC characteristics and the influence of DC on the time-averaged flow field of the APES.In Part 1 of this two-part paper,firstly a 3D viscous unsteady and time-averaging flow CFD solver is developed to investigate the APES technique.Then steady and unsteady simulations are conducted in a transonic compressor stage.The results from both simulations are compared to highlight the significance of the unsteady interactions.Furthermore,the distribution characteristics of DC are studied and the DC at the rotor/stator interface are compared with their spatial correlations (SC).Lastly,steady and time-averaging (employing APES with DC) simulations for the downstream stator alone are conducted employing DC derived from the unsteady results.The results from steady and time-averaging simulations are compared with the time-averaged unsteady results.The comparisons demonstrate that the simulation employing APES with DC can reproduce the time-averaged field and the 3D viscous time-averaging flow solver is validated.     

Turbulent wake behind a T-shaped plate: Comparison with a cross-shaped plate

The wake behind T-shaped intersecting flat plates has been studied by direct numerical simulations and compared with the wake behind intersecting plates forming a cross. The Reynolds number based on the uniform inflow velocity and the plate width d was 1000. Similar to the cross-plate the vortex shedding was suppressed in a 4d wide intersection region with a substantial base suction pressure reduction. Shear-layer (K-H) instabilities have been observed and its characteristic frequency obtained. In contrast to the cross-plate, a main feature of the mean wake structure behind the T-plate is the formation of two symmetrically positioned swirling vortices close to the internal corners of the T. This was examined by considering pressure contours and the turbulent production terms of mean streamwise vorticity. In spite of some similarities, major features of the wake behind the T-plate turned out to be distinctly different from the wake behind a cross-plate configuration.     

Base drag reduction by control of the three-dimensional unsteady vortical structures

The present paper deals with the wake of a 2D body equipped with a drag reduction device. The device is a 3D trailing edge consisting of alternate segments of blunt base and spanwise cavity. The aerodynamic mechanisms acting on the near wake are studied in a water tunnel from schlieren observations by thermally marking large scale structures. The results show that the efficiency of the device is directly related to the presence of longitudinal vortices. An optimization of the shapes in subsonic compressible flow had led to a decrease of more than 40% of the total drag of the profile.     

Effects of a trapped vortex cell on a thick wing airfoil

The effects of a trapped vortex cell (TVC) on the aerodynamic performance of a NACA0024 wing model were investigated experimentally at Re = 10⁶ and 6.67×10⁵6.67\times 10^{5}. The static pressure distributions around the model and the wake velocity profiles were measured to obtain lift and drag coefficients, for both the clean airfoil and the controlled configurations. Suction was applied in the cavity region to stabilize the trapped vortex. For comparison, a classical boundary layer suction configuration was also tested. The drag coefficient curve of the TVC-controlled airfoil showed sharp discontinuities and bifurcative behavior, generating two drag modes. A strong influence of the angle of attack, the suction rate and the Reynolds number on the drag coefficient was observed. With respect to the clean airfoil, the control led to a drag reduction only if the suction was high enough. Compared to the classical boundary layer suction configuration, the drag reduction was higher for the same amount of suction only in a specific range of incidence, i.e., α = −2° to α = 6° and only for the higher Reynolds number. For all the other conditions, the classical boundary layer suction configuration gave better drag performances. Moderate increments of lift were observed for the TVC-controlled airfoil at low incidence, while a 20% lift enhancement was observed in the stall region with respect to the baseline. However, the same lift increments were also observed for the classical boundary layer suction configuration. Pressure fluctuation measurements in the cavity region suggested a very complex interaction of several flow features. The two drag modes were characterized by typical unsteady phenomena observed in rectangular cavity flows, namely the shear layer mode and the wake mode.     

Numerical simulation of vortex-induced drag of elastic swimmer models

We present numerical simulations of simplified models for swimming organisms or robots, using chordwise flexible elastic plates. We focus on the tip vortices originating from three-dimensional effects due to the finite span of the plate. These effects play an important role when predicting the swimmer's cruising velocity, since they contribute significantly to the drag force. First we simulate swimmers with rectangular plates of different aspect ratios and compare the results with a recent experimental study. Then we consider plates with expanding and contracting shapes. We find the cruising velocity of the contracting swimmer to be higher than the rectangular one, which in turn is higher than the expanding one. We provide some evidence that this result is due to the tip vortices interacting differently with the swimmer.     

Wall proximity effects on the effectiveness of upstream control rod

Two dimensional flow over a circular cylinder with an upstream control rod of same diameter is simulated in unbound condition and in wall bounded conditions. The cylinders are placed at various heights from the wall and the inter-distance between cylinders is also varied. The control rod is subjected to different rotation rates. It is found that, in unbound condition, rotating the control rod decreases the critical pitch length (S/D_cr) and increases the drag and Strouhal number of the main cylinder. In presence of plane wall, the shielding provided by the separated shear layers from the control rod in cavity regime is deteriorated due to deflection of shear layers which results in higher drag and large fluctuation of lift coefficient. However, in wake impingement regime, the binary vortices from the control rod are weakened due to diffusion of vorticity and hence, the main cylinder experiences a lower drag and small lift fluctuations than that of unbound condition. The critical height of vortex suppression (H/D_cr) is higher in cavity regime than that of wake impingement regime due to the single extended-bluff body like configuration. The rotation of control rod energizes the wall boundary layer and increases the critical height of vortex suppression. Increasing the rotational rate of control rod decreases the drag force and reduces the amplitude of lift fluctuation. Analysis of the wall shear stress distribution reveals that it suffers a sudden drop at moderate height where the normal Karman vortex shedding changes to irregular shedding consisting of single row of negative vortices. Modal structures obtained from dynamic mode decomposition (DMD) reveal that the flow structures behind the main cylinder are suppressed due to wall and the flow is dominated by the wake of control rod.     

Reducing the pressure drag of a D-shaped bluff body using linear feedback control

The pressure drag of blunt bluff bodies is highly relevant in many practical applications, including to the aerodynamic drag of road vehicles. This paper presents theory revealing that a mean drag reduction can be achieved by manipulating wake flow fluctuations. A linear feedback control strategy then exploits this idea, targeting attenuation of the spatially integrated base (back face) pressure fluctuations. Large-eddy simulations of the flow over a D-shaped blunt bluff body are used as a test-bed for this control strategy. The flow response to synthetic jet actuation is characterised using system identification, and controller design is via shaping of the frequency response to achieve fluctuation attenuation. The designed controller successfully attenuates integrated base pressure fluctuations, increasing the time-averaged pressure on the body base by 38%. The effect on the flow field is to push the roll-up of vortices further downstream and increase the extent of the recirculation bubble. This control approach uses only body-mounted sensing/actuation and input–output model identification, meaning that it could be applied experimentally.     

Drag reduction on the 25° slant angle Ahmed reference body using pulsed jets

This paper highlights steady and unsteady measurements and flow control results obtained on an Ahmed model with slant angle of 25° in wind tunnel. On this high-drag configuration characterized by a large separation bubble along with energetic streamwise vortices, time-averaged and time-dependent results without control are first presented. The influence of rear-end periodic forcing on the drag coefficient is then investigated using electrically operated magnetic valves in an open-loop control scheme. Four distinct configurations of flow control have been tested: rectangular pulsed jets aligned with the spanwise direction or in winglets configuration on the roof end and rectangular jets or a large open slot at the top of the rear slant. For each configuration, the influence of the forcing parameters (non-dimensional frequency, injected momentum) on the drag coefficient has been studied, along with their impact on the static pressure on both the rear slant and vertical base of the model. Depending on the type and location of pulsed jets actuation, the maximum drag reduction is obtained for increasing injected momentum or well-defined optimal pulsation frequencies.     

Spanwise effects on instabilities of compressible flow over a long rectangular cavity

The stability properties of two-dimensional (2D) and three-dimensional (3D) compressible flows over a rectangular cavity with length-to-depth ratio of \(L/D=6\) are analyzed at a free-stream Mach number of \(M_\infty =0.6\) and depth-based Reynolds number of \(Re_D=502\). In this study, we closely examine the influence of three-dimensionality on the wake mode that has been reported to exhibit high-amplitude fluctuations from the formation and ejection of large-scale spanwise vortices. Direct numerical simulation (DNS) and bi-global stability analysis are utilized to study the stability characteristics of the wake mode. Using the bi-global stability analysis with the time-averaged flow as the base state, we capture the global stability properties of the wake mode at a spanwise wavenumber of \(\beta =0\). To uncover spanwise effects on the 2D wake mode, 3D DNS are performed with cavity width-to-depth ratio of \(W/D=1\) and 2. We find that the 2D wake mode is not present in the 3D cavity flow with \(W/D=2\), in which spanwise structures are observed near the rear region of the cavity. These 3D instabilities are further investigated via bi-global stability analysis for spanwise wavelengths of \(\lambda /D=0.5{-}2.0\) to reveal the eigenspectra of the 3D eigenmodes. Based on the findings of 2D and 3D global stability analysis, we conclude that the absence of the wake mode in 3D rectangular cavity flows is due to the release of kinetic energy from the spanwise vortices to the streamwise vortical structures that develops from the spanwise instabilities.     

Mechanism of drag reduction by spanwise oscillating Lorentz force in turbulent channel flow

Numerical simulations and experimental research are both carried out to investigate the controlled effect of spanwise oscillating Lorentz force on a turbulent channel flow. The variations of the streaks and the skin friction drag are obtained through the PIV system and the drag measurement system, respectively. The flow field in the near-wall region is shown through direct numerical simulations utilizing spectral method. The experimental results are consistent with the numerical simulation results qualitatively, and both the results indicate that the streaks are tilted into the spanwise direction and the drag reduction utilizing spanwise oscillating Lorentz forces can be realized. The numerical simulation results reveal more detail of the drag reduction mechanism which can be explained, since the spanwise vorticity generated from the interaction between the induced Stokes layer and intrinsic turbulent flow in the near-wall region can make the longitudinal vortices tilt and oscillate, and leads to turbulence suppression and drag reduction.     

DNS of viscoelastic turbulent channel flow with rectangular orifice at low Reynolds number

Direct numerical simulations of turbulent viscoelastic-fluid flow in a channel with a rectangular orifice were performed to investigate the influence of viscoelasticity on turbulence statistics and turbulent structures downstream of the orifice. The geometry considered is periodic rectangular orifices with 1:2 expansion. The constitutive equation follows the Giesekus model, valid for polymer (or surfactant) solutions, which are generally capable of reducing the turbulent frictional drag in a smooth channel. The friction Reynolds number and the Weissenberg number were set to 100 and 20-30, respectively. A drag reduction of about 20% was achieved in the viscoelastic flows. The onset Reynolds number for the transition from a symmetric to an asymmetric state was found to be shifted to higher values than that for the Newtonian flow. In the viscoelastic flow, the turbulent kinetic energy was decreased and fewer turbulent eddies were observed, as the Kelvin-Helmholtz vortices were quickly damped. Away from the orifice, quasi-streamwise vortices in the viscoelastic flow were sustained for a longer period, accompanied by energy exchange from elastic energy of the viscoelastic fluid to kinetic energy.