Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.

Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

Analytical calculations of Eulerian and Lagrangian time correlations in turbulent shear flows

The Taylor series expansion method is used to analytically calculate the Eulerian and Lagrangian time correlations in turbulent shear flows. The short-time behaviors of those correlation functions can be obtained from the series expansions. Especially, the propagation velocity and sweeping velocity in the elliptic model of space-time correlation are analytically calculated and further simplified using the sweeping hypothesis and straining hypothesis. These two characteristic velocities mainly determine the space-time correlations.

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]

Receptivity to free-stream disturbance waves for blunt cone axial symmetry hypersonic boundary layer

Based on high-order compact upwind scheme, a high-order shock-fitting finite difference scheme is studied to simulate the generation of boundary layer disturbance waves due to free-stream waves. Both steady and unsteady flow solutions of the receptivity problem are obtained by resolving the full Navier-Stokes equations. The interactions of bow-shock and free-stream disturbance are researched. Direct numerical simulation (DNS) of receptivity to free-stream disturbances for blunt cone hypersonic boundary layers is performed.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

Wave loading on floating platforms by internal solitary waves

Morison's equation is used for estimating internal solitary wave-induced forces exerted on SPAR and semi-submersible platforms. And the results we got have also been compared to ocean surface wave loading. It is shown that Morison's equation is an appropriate approach to estimate internal wave loading even for SPAR and semi-submersible platforms, and the internal solitary wave load on floating platforms is comparable to surface wave counterpart. Moreover, the effects of the layers with different thickness on internal solitary wave force are investigated.

An algorithm for coarse particle sedimentation simulation by Stokesian dynamics

Coarse Particle sedimentation is studied by using an algorithm with no adjustable parameters based on stokesian dynamics. Only inter-particle interactions of hydrodynamic force and gravity are considered. The sedimentation of a simple cubic array of spheres is used to verify the computational results. The scaling and parallelism with OpenMP of the method are presented. Random suspension sedimentation is investigated with Mont Carlo simulation. The computational results are shown in good agreement with experimental fitting at the lower computational cost of O(N In N).

Effects of rainfall infiltration on the stability of soil slopes

Slope failure due to rainfall is a common geotechnical problem. The mechanics of rainfall induced landslides involves the interaction of a number of complex hydrologic and geotechnical factors. This study attempts to identify the influence of some of these factors on the stability of soil slope including rainfall intensity, hydraulic conductivity and the strength parameters of soil.

Comparisons of static, quasi-static and dynamic 3D porous media scale network models for two-phase immiscible flow in porous media

A dynamic 3D pore-scale network model is formulated for investigating the effect of interfacial tension and oil-water viscosity during chemical flooding. The model takes into account both viscous and capillary forces in analyzing the impact of chemical properties on flow behavior or displacement configuration, while the static model with conventional invasion percolation algorithm incorporates the capillary pressure only. From comparisons of simulation results from these models. it indicates that the static pore scale network model can be used successfully when the capillary number is low. With the capillary increases due to the enhancement of water viscosity or decrease of interfacial tension, only the quasi-static and dynamic model can give insight into the displacement mechanisms.

Lateral motion and departure of vapor bubbles in nucleate pool boiling on thin wires in microgravity

A space experiment on bubble behavior and heat transfer in subcooled pool boiling phenomenon has been performed utilizing the temperature-controlled pool boiling (TCPB) device both in normal gravity in the laboratory and in microgravity aboard the 22(nd) Chinese recoverable satellite. The fluid is R113 at 0.1 MPa and subcooled by 26 degrees C nominally. A thin platinum wire of 60 mu m in diameter and 30mm in length is simultaneously used as heater and thermometer. Only the lateral motion and the departure of discrete vapor bubbles in nucleate pool boiling are reported and analyzed in the present paper. A scale analysis on the Marangoni convection surrounding a bubble in the process of subcooled nucleate pool boiling leads to formulas of the characteristic velocity of the lateral motion and its observability. The predictions consist with the experimental observations. Considering the Marangoni effect, a new qualitative model is proposed to reveal the mechanism underlying the bubble departure processes and a quantitative agreement can also be acquired.

Effects of temperature on the threshold of phosphorus for algal blooms

Algal blooms, worsening marine ecosystems and causing great economic loss, have been paid much attention to for a long time. Such environmental factors as light penetration, water temperature, and nutrient concentration are crucial in blooms processes. Among them, only nutrients can be controlled. Therefore, the threshold of nutrients for algal blooms is of great concern. To begin with, a dynamic eutrophication model has been constructed to simulate the algal growth and phosphorus cycling. The model encapsulates the essential biological processes of algal growth and decay, and phosphorus regeneration due to algal decay. The nutrient limitation is based upon commonly used Monod's kinetics. The effects of temperature and phosphorus limitation are particularly addressed. Then, we have endeavored to elucidate the threshold of phosphorus at different temperature for algal blooms. Based on the numerical simulation, the isoquant contours of change rate of alga as shown in the figure are obtained, which obviously demonstrate the threshold of nutrient at an arbitrary reasonable temperature. The larger the change rate is, the more rapidly the alga grows. If the phosphorus concentration at a given temperature remains larger than the threshold the algal biomass may increase monotonically, leading to the algal blooming. With the rising of temperature, the threshold is apparently reduced, which may explain why likely red tide disasters occur in a fine summer day. So, high temperature and sufficient phosphorus supply are the major factors which result in algal growth and blowout of red tide.

Generalized Viscoplastic Modeling of Debris Flow

Various concepts have been proposed or used in the development of rheological models for debris flow. The earliest model developed by Bagnold was based on the concept of the “dispersive” pressure generated by grain collisions. Bagnold’s concept appears to be theoretically sound, but his empirical model has been found to be inconsistent with most theoretical models developed from non-Newtonian fluid mechanics. Although the generality of Bagnold’s model is still at issue, debris-flow modelers in Japan have generally accepted Takahashi’s formulas derived from Bagnold’s model. Some efforts have recently been made by theoreticians in non-Newtonian fluid mechanics to modify or improve Bagnold’s concept or model. A viable rheological model should consist both of a rate-independent part and a rate-dependent part. A generalized viscoplastic fluid (GVF) model that has both parts as well as two major rheological properties (i.e., the normal stress effect and soil yield criterion) is shown to be sufficiently accurate, yet practical, for general use in debris-flow modeling. In fact, Bagnold’s model is found to be only a particular case of the GVF model. Analytical solutions for (steady) uniform debris flows in wide channels are obtained from the GVF model based on Bagnold’s simplified assumption of constant grain concentration.

Modelling of submarine landslides and generated water waves

The slide of unstable sedimentary bodies and their hydraulic effects are studied by numerical means. A two-dimensional fluid mechanics model based on Navier-Stokes equations has been developed considering the sediments and water as a mixture. Viscoplastic and diffusion laws for the sediments have been introduced into the model. The numerical model is validated with an analytical solution for a Bingham flow. Laboratory experiments consisting in the slide of gravel mass have been carried out. The results of these experiments have shown the importance of the sediment rheology and the diffusion. The model parameters are adjusted by trial and error to match the observed “sandflow”.

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.