Parallelization strategies for resolved simulations of fluid-structure-particle interactions

Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.

Wall-modeling for large-eddy simulation of flows around an axisymmetric body using the diffuse-interface immersed boundary method

A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.

The spatial-temporal evolution of coherent structures in log law region of turbulent boundary layer

The spatial-temporal evolution of coherent structures （CS） is significant for turbulence control and drag re- duction. Among the CS, low and high speed streak structures show typical burst phenomena. The analysis was based on a time series of three-dimensional and three-component （3D-3C） velocity fields of the flat plate turbulent boundary layer （TBL） measured by a Tomographic and Time-resolved PIV （Tomo TRPIV） system. Using multi-resolution wavelet transform and conditional sampling method, we extracted the intrinsic topologies and found that the streak structures appear in bar-like patterns. Furthermore, we seized locations and velocity information of transient CS, and then calculated the propagation velocity of CS based on spatial-temporal cross-correlation scanning. This laid a foundation for further studies on relevant dynamics properties.

Derivation of a Nonlinear Reynolds Stress Model Using Renormalization Group Analysis and Two-Scale Expansion Technique

Adopting Yoshizawa＇s two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

A note on the numerical simulations of flow past a wavy square-section cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Kairmain vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

How far the wake of a wind farm can persist for?

With the increased penetration of wind energy in our nation’s energy portfolio, wind farms are placed in a way close to each other. Thus, their wakes have to be fully considered in the design and operation of a wind farm. In this study, we investigate the wake of a wind farm using large-eddy simulation with wind turbine rotor modelled by the actuator disk model. The simulated results show that the wake of a wind farm can persist for a long distance in its downstream. For the considered wind farm layout, the velocity in the wake recovers 95% of that of the undisturbed inflow at 55 rotor diameters downstream from its last row, suggesting that the wake of a wind farm should be fully considered in the optimal design and operation for its downstream wind farms.

An overview of fluids mixing in T-shaped mixers

In a T-shaped mixer,the two liquid streams in the inlet channels meet each other at the T-junction,and their liquid-liquid contacting face exhibits planar,swirling folds and the folds breaking to be chaos and turbulence,as the Reynolds number increases.The characteristic mixing scenario attracts long-time attention,given these mixings are of fundamental importance in fluid physics and also have been successfully used in engineering applications.The experimental and numerical studies of flow features and mixing characteristics in T-mixers are overviewed in this manuscript.This review introduces the experimental and numerical techniques in the studies,the flow and mixing characteristics in the corresponding regimes and application examples of the T-mixers at last,aiming at introducing fundamentals to researchers with initial interests on this topic.

Multi-Scale-Matching neural networks for thin plate bending problem

Physics-informed neural networks are a useful machine learning method for solving differential equations,but encounter challenges in effectively learning thin boundary layers within singular perturbation problems.To resolve this issue,multi-scale-matching neural networks are proposed to solve the singular perturbation problems.Inspired by matched asymptotic expansions,the solution is decomposed into inner solutions for small scales and outer solutions for large scales,corresponding to boundary layers and outer regions,respectively.Moreover,to conform neural networks,we introduce exponential stretched variables in the boundary layers to avoid semiinfinite region problems.Numerical results for the thin plate problem validate the proposed method.

A two-parameter model to predict fatigue life of high-strength steels in a very high cycle fatigue regime

In this paper, ultrasonic （20 kHz） fatigue tests were performed on specimens of a high-strength steel in very high cycle fatigue （VHCF） regime. Experimental results showed that for most tested specimens failed in a VHCF regime, a fatigue crack originated from the interior of specimen with a fish-eye pattern, which contained a fine granular area （FGA） centered by an inclusion as the crack origin. Then, a two-parameter model is proposed to predict the fatigue life of high-strength steels with fish-eye mode failure in a VHCF regime, which takes into account the inclusion size and the FGA size. The model was verified by the data of present experiments and those in the literature. Furthermore, an analytic formula was obtained for estimating the equivalent crack growth rate within the FGA. The results also indicated that the stress intensity factor range at the front of the FGA varies within a small range, which is irrespective of stress amplitude and fatigue life.

Self-propulsion of flapping bodies in viscous fluids: Recent advances and perspectives

Flapping-powered propulsion is used by many animals to locomote through air or water. Here we review recent experimental and numerical studies on self-propelled mechanical systems powered by a flapping motion. These studies improve our understanding of the mutual interaction between actively flapping bodies and surrounding fluids. The results obtained in these works provide not only new insights into biolocomotion but also useful information for the biomimetic design of artificial flyers and swimmers.

A large eddy simulation of flows around an underwater vehicle model using an immersed boundary method

A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.

Shedding Phenomenon of Ventilated Partial Cavitation around an Underwater Projectile

A new shedding phenomenon of ventilated partial cavitations is observed around an axisymmetric projectile in a horizontal launching experiment.The experiment system is established based on SHPB launching and high speed photography.A numerical simulation is carried out based on the homogeneous mixture approach,and its predicted evolutions of cavities are compared with the experimental results.The cavity breaks off by the interaction between the gas injection and the re-entry jet at the middle location of the projectile,which is obviously different from natural cavitation.The mechanism of cavity breaking and shedding is investigated,and the influences of important factors are also discussed.

Scaling law of resolved-scale isotropic turbulence and its application in large-eddy simulation

Eddy-damping quasinormal Markovian （EDQNM） theory is employed to calculate the resolved-scale spectrum and transfer spectrum, based on which we investigate the resolved-scale scaling law. Results show that the scaling law of the resolved-scale turbulence, which is affected by several factors, is far from that of the full-scale turbulence and should be corrected. These results are then applied to an existing subgrid model to improve its performance. A series of simulations are performed to verify the necessity of a fixed scaling law in the subgrid modeling.

Experimental study on spectrum and multi-scale nature of wall pressure and velocity in turbulent boundary layer

When using a miniature single sensor boundary layer probe, the time sequences of the stream-wise velocity in the turbulent boundary layer （TBL） are measured by using a hot wire anemometer. Beneath the fully developed TBL, the wall pressure fluctuations are attained by a microphone mechanism with high spatial resolution. Analysis on the statistic and spectrum properties of velocity and wall pressure reveals the relationship between the wall pressure fluctuation and the energy-containing structure in the buffer layer of the TBL. Wavelet transform shows the multi-scale natures of coherent structures contained in both signals of velocity and pressure. The most intermittent wall pressure scale is associated with the coherent structure in the buffer layer. Meanwhile the most energetic scale of velocity fluctuation at y＋ = 14 provides a specific frequency f9 ≈ 147 Hz for wall actuating control with Ret = 996.

Convection and correlation of coherent structure in turbulent boundary layer using tomographic particle image velocimetry

The present experimental work focuses on a new model for space-time correlation and the scale-dependencies of convection velocity and sweep velocity in turbulent boundary layer over a flat wail. A turbulent boundary layer flow at Reo = 2460 is measured by tomographic particle image velocimetry （tomographic PIV）. It is demonstrated that arch, cane, and hairpin vortices are dominant in the logarithmic layer. Hairpins and hairpin packets are responsible for the elongated low-momentum zones observed in the instantaneous flow field. The conditionally-averaged coherent structures systemically illustrate the key roles of hairpin vortice in the turbulence dynamic events, such as ejection and sweep events and energy transport. The space-time correlations of instantaneous streamwise fluctuation velocity are calculated and confirm the new elliptic model for the space-time correlation instead of Taylor hypothesis. The convection velocities derived from the space-time correlation and conditionally-averaged method both suggest the scaling with the local mean velocity in the logarithmic layer. Convection velocity result based on Fourier decomposition （FD） shows stronger scale- dependency in the spanwise direction than in streamwise direction. Compared with FD, the proper orthogonal decomposition （POD） has a distinct distribution of convection velocity for the large- and small-scales which are separated in light of their contributions of turbulent kinetic energy.

Wall-resolved large-eddy simulation of turbulent channel flows with rough walls

Turbulent flows over rough surfaces widely exist in nature and industry.Investigating its mechanism is of theoretical and practical significance.In this work we simulate the turbulent channel flow with rough walls using large-eddy simulation with rough elements resolved using the curvilinear immersed boundary method and compare the results obtained in this work with those in the paper by Yuan and Piomelli(J.Fluid Mech.,vol.760,pp.R1,2014),where the volume of fluid method was employed for modeling rough elements.The mean streamwise velocity profiles predicted by the two methods agree well with each other.Differences in Reynolds stresses and dispersive stresses are observed,which are attributed to the different approaches in dealing with the complex geometry of the rough surface.

A study on coherent structures and drag-reduction in the wall turbulence with polymer additives by TRPIV

An experimental measurement was performed us- ing time-resolved particle image velocimetry （TRPIV） to in- vestigate the spatial topological character of coherent struc- tures in wall-bounded turbulence of polymer additive solu- tion. The fully developed near-wall turbulent flow fields with and without polymer additives at the same Reynolds number were measured by TRPIV in a water channel. The compar- isons of turbulent statistics confirm that due to viscoelastic structure of long-chain polymers, the wall-normal velocity fluctuation and Reynolds shear stress in the near-wall region are suppressed significantly. Furthermore, it is noted that such a behavior of polymers is closely related to the decease of the motion of the second and forth quadrants, i.e., the ejection and sweep events, in the near-wall region. The spa- tial topological mode of coherent structures during bursts has been extracted by the new mu-level criteria based on locally averaged velocity structure function. Although the general shapes of coherent structures are unchanged by polymer additives, the fluctuating velocity, velocity gradient, velocity strain rate and vorticity of coherent structures during burst events are suppressed in the polymer additive solution com- pared with that in water. The results show that due to the polymer additives the occurrence and intensity of coherent structures are suppressed, leading to drag reduction.

Towards multi-fidelity simulation of flows around an underwater vehicle with appendages and propeller

Flow around a real-life underwater vehicle often happens at a high Reynolds number with flow structures at different scales from the boundary layer around a blade to that around the hull. This poses a great challenge for large-eddy simulation of an underwater vehicle aiming at resolving all relevant flow scales. In this work, we propose to model the hull with appendages using the immersed boundary method, and model the propeller using the actuator disk model without resolving the geometry of the blade. The proposed method is then applied to simulate the flow around Defense Advanced Research Projects Agency(DARPA) suboff. An overall acceptable agreement is obtained for the pressure and friction coefficients. Complex flow features are observed in the near wake of suboff. In the far wake, the core region is featured by a jet because of the actuator disk, surrounded by an annular region with velocity deficit due to the body of suboff.

Effects of laser power density on static and dynamic mechanical properties of dissimilar stainless steel welded joints

The mechanical properties of laser welded joints under impact loadings such as explosion and car crash etc. are critical for the engineering designs. The hardness, static and dynamic mechanical properties of AISI304 and AISI316L dissimilar stainless steel welded joints by CO2 laser were experimentally studied. The dynamic strainstress curves at the strain rate around 103 s-1 were obtained by the split Hopkinson tensile bar （SHTB）. The static mechanical properties of the welded joints have little changes with the laser power density and all fracture occurs at 316 L side. However, the strain rate sensitivity has a strong depen- dence on laser power density. The value of strain rate factor decreases with the increase of laser power density. The welded joint which may be applied for the impact loading can be obtained by reducing the laser power density in the case of welding quality assurance.

Numerical simulation of unsteady flows over a slow-flying bat

This letter describes numerical simulation of the unsteady flow over a slow-flying bat by using the immersed boundary method based on the measured bat wing geometry and kinematics. The main vortical structures around the bat flapping wings are identified, illuminating the lift-generating role of the leading- edge vortices generated mainly in the downstroke. Furthermore, the lift decomposition indicates that the vortex lift has the dominant contribution to the time-averaged lift and the lift associated with the fluid acceleration has the relatively moderate effect.