An analysis of subgrid-resolved scale interactions with use of results from direct numerical simulations

Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range; By decomposing the subgrid energy transfer and nonlinear interaction into 'forward' and 'backward' groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.

Teracluster LSSC-II - Its Designing Principles and Applications in Large Scale Numerical Simulations

The teracluster LSSC-II installed at the State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences is one of the most powerful PC clusters in China. It has a peek performance of 2Tflops. With a Linpack performance of 1.04Tflops, it is ranked at the 43rd place in the 20th TOP500 List (November 2002), 51st place in the 21st TOP500 List (June 2003), and the 82nd place in the 22nd TOP500 List (November 2003) with a new Linpack performance of 1.3Tflops. In this paper, we present some design principles of this cluster, as well as its applications in some largescale numerical simulations.

Large Scale Numerical Simulations For Multi-Phase Fluid Dynamics With Moving Interfaces

This short communication presents our recent studies to implement numerical simulations for multi-phase flows on top-ranked supercomputer systems with distributed memory architecture. The numerical model is designed so as to make full use of the capacity of the hardware. Satisfactory scalability in terms of both the parallel speed-up rate and the size of the problem has been obtained on two high rank systems with massively parallel processors, the Earth Simulator (Earth simulator research center, Yokohama Kanagawa, Japan) and the TSUBAME (Tokyo Institute of Technology, Tokyo, Japan) supercomputers.

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 Karman 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.

Numerical simulations of interrupted and conventional microchannel heat sinks

We provide three-dimensional numerical simulations of conjugate heat transfer in conventional and the newly proposed interrupted microchannel heat sinks. The new microchannel heat sink consists of a set of separated zones adjoining shortened parallel microchannels and transverse microchambers. Multi-channel effect, physical property variations, and axial thermal conduction are considered. It is found that flow rate variations in different channels can be neglected, while heat received by different channels accounts for 2% deviations from the averaged value when the heat flux at the back surface of the silicon chip reaches 100 W/cm(2). The computed hydraulic and thermal boundary layers are redeveloping in each separated zone due to shortened flow length for the interrupted microchannel heat sink. The periodic thermal developing flow is responsible for the significant heat transfer enhancement. Two effects influence pressure drops across the newly proposed microchannel heat sink. The first one is the pressure recovery effect in the microchamber, while the second one is the head loss when liquid leaves the microchamber and enters the next zone. The first effect compensates or suppresses the second one, leading to similar or decreased pressure drop than that for the conventional microchannel heat sink, with the fluid Prandtl number larger than unity.

Two-dimensional numerical simulations of shock waves in micro convergent-divergent nozzles

The steady two-dimensional Navier-Stokes equations with the slip wall boundary conditions were used to simulate the supersonic flow in micro convergent-divergent nozzles. It is observed that shock waves can take place inside or outside of the micronozzles under the earth environment. For the over-expanded flows, there is a boundary layer separation point, downstream of which a wave interface separates the viscous boundary layer with back air flow and the inviscid core flow. The oblique shock wave is followed by the bow shock and shock diamond. The viscous boundary layer thickness relative to the whole nozzle width on the exit plane is increased but attains the maximum value around of 0.5 and oscillates against this value with the continuous increasing of the nozzle upstream pressures. The viscous effect either changes the normal shock waves outside of the nozzle for the inviscid flow to the oblique shock waves inside the nozzle, or transfers the expansion jet flow without shock waves for the inviscid flow to the oblique shock waves outside of the nozzle. 

Physical mechanism and numerical simulations of surface layer temperature inversion in tropical ocean

The one-dimensional Kraus-Turner mixed layer model improved by Liu is developed to consider the effect of salinity and the equations of temperature and salinity under the mixed layer. On this basis, the processes of growth and death of surface layer temperature inversion is numerically simulated under different environmental parameters. At the same time, the physical mechanism is preliminarily discussed combining the observations at the station of TOGA-COARE 0 degrees N, 156 degrees E. The results indicate that temperature inversion sensitively depends on the mixed layer depth, sea surface wind speed and solar shortwave radiation, etc., and appropriately meteorological and hydrological conditions often lead to the similarly periodical occurrence of this inversion phenomenon.

Factors influencing the climatological mixed layer depth in the South China Sea: numerical simulations

The mixed layer depth (MLD) in the upper ocean is an important physical parameter for describing the upper ocean mixed layer. We analyzed several major factors influencing the climatological mixed layer depth (CMLD), and established a numerical simulation in the South China Sea (SCS) using the Regional Ocean Model System (ROMS) with a high-resolution (1/12A degrees x1/12A degrees) grid nesting method and 50 vertical layers. Several ideal numerical experiments were tested by modifying the existing sea surface boundary conditions. Especially, we analyzed the sensitivity of the results simulated for the CMLD with factors of sea surface wind stress (SSWS), sea surface net heat flux (SSNHF), and the difference between evaporation and precipitation (DEP). The result shows that of the three factors that change the depth of the CMLD, SSWS is in the first place, when ignoring the impact of SSWS, CMLD will change by 26% on average, and its effect is always to deepen the CMLD; the next comes SSNHF (13%) for deepening the CMLD in October to January and shallowing the CMLD in February to September; and the DEP comes in the third (only 2%). Moreover, we analyzed the temporal and spatial characteristics of CMLD and compared the simulation result with the ARGO observational data. The results indicate that ROMS is applicable for studying CMLD in the SCS area.

FEM Simulations and Optimization about Residual Stresses in Coating Structures with Functionally Graded Materials Layer

An elasto-plastic finite element method is developed to predict the residual stresses of thermal spraying coatings with functionally graded material layer. In numerical simulations, temperature sensitivity of various material constants is included and mix

Simulations and Experiments of Stochastic Characteristics for Collective Short Fatigue Cracks in Steels

Stochastic characteristics prevail in the process of short fatigue crack progression. This paper presents a method taking into account the balance of crack number density to describe the stochastic behaviour of short crack collective evolution. The results from the simulation illustrate the stochastic development of short cracks. The experiments on two types of steels show the random distribution for collective short cracks with the number of cracks and the maximum crack length as a function of different locations on specimen surface. The experiments also give the variation of total number of short cracks with fatigue cycles. The test results are consistent with numerical simulations.

Direct Numerical Simulation of Supersonic Turbulent Boundary Layer Flow

Direct numerical simulations of a spatially evolving supersonic flat-plate turbulent boundary layer flow with free Mach number M = 2.25 and Reynolds number Re = 365000/in are performed. The transition process from laminar to turbulent flow is obtained by solving the three-dimensional compressible Navier-Stokes, equations, using high-order accurate difference schemes. The obtained statistical results agree well with the experimental and theoretical data. From the numerical results it can be seen that the transition process under the considered conditions is the process which skips the Tolimien-Schlichting instability and the second instability through the instability of high gradient shear layer and becomes of laminar flow breakdown. This means that the transition process is a bypass-type transition process. The spanwise asymmetry of the disturbance locally upstream imposed is important to induce the bypass-type transition. Furthermore, with increasing the time disturbance frequency the transition will delay. When the time disturbance frequency is large enough, the transition will disappear.