DNS and LES of turbulent channel flow with hydrophobic surface

Hydrophobic surface benefits for drag reduction. Min and Kim[1] do the first Direct Numerical Simulation on drag reduction in turbulent channel flow. And Fukagata and Kasagi[2] make some theoretical analysis based on Dean[3]'s formula and some observations in the DNS results. Using their theory, they conclude that drag reduction is possible in large Reynolds number. Both Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed in our research. How the LES behaving in the turbulent channel flow with hydrophobic surface is examined. Original Smagorinsky model and its Dynamical model are used in LES. The slip velocities predicted by LES using Dynamical model are in good agreement with DNS as shown in the Figure. Although the percentage of drag reduction predicted by LES shows some discrepancies, it is in the error limit for industrial flow. First order and second order moments of LES are also examined and compared with DNS's results. The first-order moments is calculated well by LES. But there are some discrepancies of second-order moments between LES and DNS. [GRAPHICS]

Multi-pipe string electromagnetic detection tool and its applications

The MID-K, a new kind of multi-pipe string detection tool is introduced. This tool provides a means of evaluating the condition of in-place pipe string, such as tubing and casino. It is capable of discriminating the defects of the inside and outside, and estimating the thickness of tubing and casing. It is accomplished by means of a low frequency eddy current to detect flaws on the inner surface and a magnetic flux leakage to inspect the full thickness. The measurement principle, the technology and applications are presented in this paper.

Net phytoplankton and zooplankton in the New York Bight, January 1976 to February 1978, with comments on the effects of wind, Gulf Stream eddies, and slope water intrusions

Results are given of monthly net phytoplankton and zooplankton sampling from a 10 m depth in shelf, slope, and Gulf Stream eddy water along a transect running southeastward from Ambrose Light, New York, in 1976, 1977, and early 1978. Plankton abundance and temperature at 10 m and sea surface salinity at each station are listed. The effects of atmospheric forcing and Gulf Stream eddies on plankton distribution and abundance arc discussed. The frequency of Gulf Stream eddy passage through the New York Bight corresponded with the frequency of tropical-subtropical net phytoplankton in the samples. Gulf Stream eddies injected tropical-subtropical zooplankton onto the shelf and removed shelfwater and its entrained zooplankton. Wind-induced offshore Ekman transport corresponded generally with the unusual timing of two net phytoplankton maxima. Midsummer net phytoplankton maxima were recorded following the passage of Hurricane Belle (August 1976) and a cold front (July 1977). Tropical-subtropical zooplankton which had been injected onto the outer shelf by Gulf Stream eddies were moved to the inner shelf by a wind-induced current moving up the Hudson Shelf Valley. (PDF file contains 47 pages.)

Migrations of yellowfin tuna tagged off the southern coast of Mexico in 1960 and 1969

ENGLISH: The map method, the Jones method, the variance-covariance method, and the Skellam method were used to study the migrations of tagged yellowfin tuna released off the southern coast of Mexico in 1960 and 1969. The first three methods are all useful, and each presents information which is complementary to that presented by the others. The Skellam method, as used in this report, is less useful. The movements of the tagged fish released in 1960 appeared to have been strongly directed, but this was probably caused principally by the distribution of the fishing effort. The effort was much more widely distributed in 1970, and the movements of the fish released in 1969 appeared to have been much less directed. The correlation coefficients derived from the variance-covariance method showed that it was not random, however. The small fish released in the Acapulco and 10°N-100°W areas in 1969 migrated to the Manzanillo area near the beginning of February 1970. The medium and large fish released in the same areas in the same year tended to migrate to the southeast throughout the first half of 1970, however. SPANISH: El método de mapas, el de Jones, el de la variancia-covariancia y el de Skellam fueron empleados para estudiar las migraciones del atún aleta amarilla marcado y liberado frente a la costa meridional de México en 1960 y 1969. Los tres primeros métodos son todos útiles, y cada uno presenta información que complementa la presentada por los otros. El método de Skellam, conforme se usa en este informe, es menos útil. Parece que los desplazamientos de los peces marcados y liberados en 1960 hubieran sido fuertemente orientados, pero ésto probablemente fue causado principalmente por la distribución del esfuerzo de pesca. El esfuerzo se distribuyó más extensamente en 1970, y parece que los desplazamientos de los peces liberados en 1969 fueran menos orientados. Los coeficientes de correlación derivados del método variancia-covariancia indicaron, sin embargo, que no eran aleatorios. Los peces pequeños liberados en las áreas de Acapulco y los 10°N-100°W en 1969 migraron al área de Manzanillo a principios de febrero 1970. Los peces medianos y grandes liberados en las mismas áreas en el mismo año tuvieron, sin embargo, la tendencia a desplazarse al sudeste durante el primer semestre de 1970. (PDF contains 64 pages.)

A kinematic subgrid scale model for large-eddy simulation of turbulence generated sound

In the hybrid approach of large-eddy simulation (LES) and Lighthill’s acoustic analogy for turbulence-generated sound, the turbulence source fields are obtained using an LES and the turbulence-generated sound at far fields is calculated from Lighthill’s acoustic analogy. As only the velocity fields at resolved scales are available from the LES, the Lighthill stress tensor, serving as a source term in Lighthill’s acoustic equation, has to be evaluated from the resolved velocity fields. As a result, the contribution from the unresolved velocity fields is missing in the conventional LES. The sound of missing scales is shown to be important and hence needs to be modeled. The present study proposes a kinematic subgrid-scale (SGS) model which recasts the unresolved velocity fields into Lighthill’s stress tensors. A kinematic simulation is used to construct the unresolved velocity fields with the imposed temporal statistics, which is consistent with the random sweeping hypothesis. The kinematic SGS model is used to calculate sound power spectra from isotropic turbulence and yields an improved result: the missing portion of the sound power spectra is approximately recovered in the LES.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

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.

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included groundbased tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.