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.

An Improved Ce/Se Scheme For Multi-Material Elastic-Plastic Flows And Its Applications

In this paper, a new definition of SE and CE, which is based on the hexahedron mesh and simpler than Chang's original CE/SE method (the space-time Conservation Element and Solution Element method), is proposed and an improved CE/SE scheme is constructed. Furthermore, the improved CE/SE scheme is extended in order to solve the elastic-plastic flow problems. The hybrid particle level set method is used for tracing the interfaces of materials. Proper boundary conditions are presented in interface tracking. Two high-velocity impact problems are simulated numerically and the computational results are carefully compared with the experimental data, as well as the results from other literature and LS-DYNA software. The comparisons show that the computational scheme developed currently is clear in physical concept, easy to be implemented and high accurate and efficient for the problems considered. (C) 2008 Elsevier Ltd. All rights reserved.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

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.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

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.

Single-mode behaviour judgment of optical waveguides by imaginary-distance beam propagation method under perfectly matched layer boundary condition

Imaginary-distance beam propagation method under the perfectly matched layer boundary condition is applied to judge single-mode behaviour of optical waveguides, for the first time to our knowledge. A new kind of silicon-on-insulator-based rib structures with half-circle cross-section is presented. The single-mode behaviour of this kind of waveguide with radius 2mum is investigated by this method. It is single-mode when the slab height is not smaller than the radius.

New insight into the origin of twin and grain boundary in InP

X-ray reflectivity curves show bi-crystal (twin) characteristics. Defect segregations at the twin boundary can be seen, whereas stress is relaxed at the edge of the boundary. Relaxation of the stress resulted in the formation of twins and other defects. As a result of the formation of such defects, a defect-free and stress-free zone or low defect density and small stress zone is created around the defects. Stress, chemical stoichiometry deviation and non-homogeneous distribution of impurities are the key factors that cause twins in LEC InP crystal growth. (C) 1999 Elsevier Science Ltd. All rights reserved.

Interaction of linear waves with infinitely long horizontal cylinders studied by boundary element method

The two-dimensional problems concerning the interaction of linear water waves with cylinders of arbitrary shape in two-layer deep water are investigated by use of the Boundary Integral Equation method (BIEM). Simpler new expressions for the Green functions are derived, and verified by comparison of results obtained by BIEM with these by an analytical method. Examined are the radiation and scattering of linear waves by two typical configurations of cylinders in two-layer deep water. Hydrodynamic behaviors including hydrodynamic coefficients, wave forces, reflection and transmission coefficients and energies are analyzed in detail, and some interesting physical phenomena are observed.

Boundary conditions at the solid-liquid surface over the multiscale channel size from nanometer to micron

The boundary condition at the solid surface is one of the important problems for the microfluidics. In this paper we study the effects of the channel sizes on the boundary conditions (BC), using the hybrid computation scheme adjoining the molecular dynamics (MD) simulations and the continuum fluid mechanics. We could reproduce the three types of boundary conditions (slip, no-slip and locking) over the multiscale channel sizes. The slip lengths are found to be mainly dependent on the interfacial parameters with the fixed apparent shear rate. The channel size has little effects on the slip lengths if the size is above a critical value within a couple of tens of molecular diameters. We explore the liquid particle distributions nearest the solid walls and found that the slip boundary condition always corresponds to the uniform liquid particle distributions parallel to the solid walls, while the no-slip or locking boundary conditions correspond to the ordered liquid structures close to the solid walls. The slip, no-slip and locking interfacial parameters yield the positive, zero and negative slip lengths respectively. The three types of boundary conditions existing in "microscale" still occur in "macroscale". However, the slip lengths weakly dependent on the channel sizes yield the real shear rates and the slip velocity relative to the solid wall traveling speed approaching those with the no-slip boundary condition when the channel size is larger than thousands of liquid molecular diameters for all of the three types of interfacial parameters, leading to the quasi-no-slip boundary conditions.

A New Criterion Number for the Boundary Conditions at the Solid/Liquid Interface in Nanoscale

A simple, but important three-atom model was proposed at the solid/liquid interface, leading to a new criterion number, lambda, governing the boundary conditions (BCs) in nanoscale. The solid wall is considered as the face-centered-cubic (fcc) structure. The fluid is the liquid argon with the well-known LJ potential. Based on the concept, the two micro-systems have the same BCs if they have The same criterion number. The degree of the locking BCs is enhanced when lambda equals to 0.757. Such critical criterion number results in the substantial epitaxial ordering and one, two, or even three liquid layers are locked by the solid wall, depending on the coupling energy scale ratio of the solid and liquid atoms. With deviation from the critical criterion number, the flow approaches the slip BCs and there are little ordering structures within the liquid. Always at the same criterion number, the degree of the slip is decreased or the locking is enhanced with increasing the coupling energy scale ratio of the solid and liquid atoms. The above analysis is well confirmed by the molecular dynamics (MD) simulation. The slip length is well correlated in terms of the new criterion number. The future work is suggested to extend the present theory for other microstructures of the solid wall atoms and quasi-LJ potentials.