A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

Effects of wall temperature on boundary layer stability over a blunt cone at Mach 7.99

Effects of wall temperature on stabilities of hypersonic boundary layer over a 7-degree half-cone-angle blunt cone are studied by using both direct numerical simulation (DNS) and linear stability theory (LST) analysis. Four isothermal wall cases with Tw/T0= 0.5, 0.7, 0.8 and 0.9, as well as an adiabatic wall case are considered. Results of both DNS and LST indicate that wall temperature has significant effects on the growth of disturbance waves. Cooling the surface accelerates unstable Mack II mode waves and decelerates the first mode (Tollmien–Schlichting mode) waves. LST results show that growth rate of the most unstable Mack II mode waves for the cases of cold wall Tw/T0=0.5 and 0.7 are about 45% and 25% larger than that for the adiabatic wall, respectively. Numerical results show that surface cooling modifies the profiles of rdut/dyn and temperature in the boundary layers, and thus changes the stability haracteristic of the boundary layers, and then effects on the growth of unstable waves. The results of DNS indicate that the disturbances with the frequency range from about 119.4 to 179.1 kHz, including the most unstable Mack modes, produce strong mode competition in the downstream region from about 11 to 100 nose radii. And adiabatic wall enhances the amplitudes of disturbance according to the results of DNS, although the LST indicates that the growth rate of the disturbance of cold wall is larger. That because the growth of the disturbance does not only depend on the development of the second unstable mode.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Euler Solution Using Adaptive Cartesian Grid With A Gridless Boundary Treatment

A quadtree-based adaptive Cartesian grid generator and flow solver were developed. The grid adaptation based on pressure or density gradient was performed and a gridless method based on the least-square fashion was used to treat the wall surface boundary condition, which is generally difficult to be handled for the common Cartesian grid. First, to validate the technique of grid adaptation, the benchmarks over a forward-facing step and double Mach reflection were computed. Second, the flows over the NACA 0012 airfoil and a two-element airfoil were calculated to validate the developed gridless method. The computational results indicate the developed method is reasonable for complex flows.

Acoustic Calculation for Supersonic Turbulent Boundary Layer Flow

An approach which combines direct numerical simulation (DNS) with the Lighthill acoustic analogy theory is used to study the potential noise sources during the transition process of a Mach 2.25 flat plate boundary layer. The quadrupole sound sources due to the flow fluctuations and the dipole sound sources due to the fluctuating surface stress are obtained. Numerical results suggest that formation of the high shear layers leads to a dramatic amplification of amplitude of the fluctuating quadrupole sound sources. Compared with the quadrupole sound source, the energy of dipole sound source is concentrated in the relatively low frequency range.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.

Experimental Study Of Vortex-Induced Vibrations Of A Cylinder Near A Rigid Plane Boundary In Steady Flow

In this study, the vortex-induced vibrations of a cylinder near a rigid plane boundary in a steady flow are studied experimentally. The phenomenon of vortex-induced vibrations of the cylinder near the rigid plane boundary is reproduced in the flume. The vortex shedding frequency and mode are also measured by the methods of hot film velocimeter and hydrogen bubbles. A parametric study is carried out to investigate the influences of reduced velocity, gap-to-diameter ratio, stability parameter and mass ratio on the amplitude and frequency responses of the cylinder. Experimental results indicate: (1) the Strouhal number (St) is around 0.2 for the stationary cylinder near a plane boundary in the sub-critical flow regime; (2) with increasing gap-to-diameter ratio (e (0)/D), the amplitude ratio (A/D) gets larger but frequency ratio (f/f (n) ) has a slight variation for the case of larger values of e (0)/D (e (0)/D > 0.66 in this study); (3) there is a clear difference of amplitude and frequency responses of the cylinder between the larger gap-to-diameter ratios (e (0)/D > 0.66) and the smaller ones (e (0)/D < 0.3); (4) the vibration of the cylinder is easier to occur and the range of vibration in terms of V (r) number becomes more extensive with decrease of the stability parameter, but the frequency response is affected slightly by the stability parameter; (5) with decreasing mass ratio, the width of the lock-in ranges in terms of V (r) and the frequency ratio (f/f (n) ) become larger.