Must thermodynamic functions be piecewise analytic?

There seems to be no general theoretical argument to support the idea that thermodynamic functions are piecewise analytic. We suggest that nonanalyticity may be associated with Gibbs states which are quasiperiodic under space translations, or have a more general nonperiodic (turbulent) behavior.     

Toward low‐noise synthetic turbulent inflow conditions for aeroacoustic calculations

The accuracy of boundary conditions for computational aeroacoustics is a well‐known challenge, due in part to the necessity of truncating the flow domain and replacing the analytical boundary conditions at infinity with numerical boundary conditions. In particular, the inflow boundary condition involving turbulent velocity or scalar fields is likely to introduce spurious waves into the domain, therefore degrading the flow behavior and deteriorating the physical acoustic waves. In this work, a method to generate low‐noise, divergence‐free, synthetic turbulence for inflow boundary conditions is proposed. It relies on the classical view of turbulence as a superposition of random eddies convected with the mean flow. Within the proposed model, the vector potential and the requirement that the individual eddies must satisfy the linearized momentum equations about the mean flow are used. The model is tested using isolated eddies convected through the inflow boundary and an experimental benchmark data for spatially decaying isotropic turbulence. Copyright © 2013 John Wiley & Sons, Ltd.     

Dynamic release process of pollutants during suspended sediment transport in aquatic system

Pollutants release is highly consistent with suspended sediment concentration （SSC） in water column, especially during re-suspension and transport events. The present research focuses on pollutant dynamic release from re-suspended sediment, especially the vertical distribution relationship between them. The sediment erosion experiments on a series of uniform flow are conducted in a circulate flume. Reactive tracer （phosphorus） is used as the contaminant in fine-grained sediments to identify the release characteristic length and time. Experimental results show that the flow condition near-bed depends on the sediment surface roughness. The region with high turbulent intensities corresponds to a high concentration sediment layer. In addition, the SSC decreases with the distance, water depth, and particle grain size. The sediment in a smaller grain size takes much more time to reach equilibrium concentration. Total phosphorus （TP） concentration changes along the water depth as SSC in the initial re-suspension stage, appearing in two obvious concentration regimes： the upper low-concentration layer and the high-concentration near-bottom layer. This layered phenomenon remains for about 3 hours until SSC distri- bution tends to be uniform. Longitudinal desorption plays an important role in long-way transport to reduce the amount of suspended sediment in water column.     

Birth–death process of local structures in defect turbulence described by the one-dimensional complex Ginzburg–Landau equation

Defect turbulence described by the one-dimensional complex Ginzburg–Landau equation is investigated and analyzed via a birth–death process of the local structures composed of defects, holes, and modulated amplitude waves (MAWs). All the number statistics of each local structure, in its stationary state, are subjected to Poisson statistics. In addition, the probability density functions of interarrival times of defects, lifetimes of holes, and MAWs show the existence of long-memory and some characteristic time scales caused by zigzag motions of oscillating traveling holes. The corresponding stochastic process for these observations is fully described by a non-Markovian master equation.     

Turbulence modelling and role of compressibility on oil spilling from a damaged double hull tank

The viscosity plays an important role, and a multiphase solver is necessary to numerically simulate the oil spilling from a damaged double hull tank (DHT). However, it is uncertain whether turbulence modelling is necessary, which turbulence model is suitable; and what the role of compressibility of the fluids is. This paper presents experimental and numerical investigations to address these issues for various cases representing different scenarios of the oil spilling, including grounding and collision. In the numerical investigations, various approaches to model the turbulence, including the large eddy simulation (LES), direct numerical simulation and the Reynolds average Navier–Stokes equation (RANS) with different turbulence models, are employed. Based on the investigations, it is suggested that the effective Reynolds numbers corresponding to both oil outflow and water inflow shall be considered when classifying the significance of the turbulence and selecting the appropriate turbulence models. This is confirmed by new lab tests considering the axial offset between the internal and the external holes on two hulls of the DHT. The investigations conclude for numerically simulating oil spilling from a damaged DHT that when the effective Re is smaller the RANS approaches should not be used and LES modelling should be employed; while when the effective Reynolds numbers is large, the RANS models may be used as they can give similar results to LES in terms of the height of the mixture in the ballast tank and discharge but costing much less CPU time. The investigation on the role of the compressibility of the fluid reveals that the compressibility of the fluid may be considerable in a small temporal‐spatial scale but plays an insignificant role on macroscopic process of the oil spilling. Copyright © 2016 John Wiley & Sons, Ltd.     

Discrete filter operators for large‐eddy simulation using high‐order spectral difference methods

The combination of a high‐order unstructured spectral difference (SD) spatial discretization scheme with sub‐grid scale (SGS) modeling for large‐eddy simulation is investigated with particular focus on the consistent implementation of a structural mixed model based on the scale similarity hypothesis. The difficult task of deriving a consistent formulation for the discrete filter within the SD element of arbitrary order led to the development of a new class of three‐dimensional constrained discrete filters. The discrete filters satisfy a set of selected criteria and are completely local within the SD element. Their weights can be automatically computed at run time from the number of solution points within each element and the expected filter cutoff length scale. The novel discrete filters can be applied to any SGS model involving explicit filtering and to a broad class of high‐order discontinuous finite element numerical schemes. The code is applied to the computation of turbulent channel flows at three Reynolds numbers, namely Re_τ = 180, 395, and 590 (based on the friction velocity u_τ and channel half‐width δ). Results from computations with and without the SGS model are compared against results from direct numerical simulation. The numerical experiments suggest that the results are sensitive to the use of the SGS model, even when a high‐order numerical scheme is used, especially when the grid resolution is kept relatively low and mostly in terms of resolved Reynolds stresses. Results obtained using existing filters based on the projection of the solution over lower‐order polynomial bases are also shown and demonstrate that these filters are inadequate for SGS modeling purposes, mostly because of their inability to enforce the selected cutoff length scale with sufficient accuracy. The use of the similarity mixed formulation proved to be particularly accurate in reproducing SGS interactions, confirming that its well‐known potential can be realized in conjunction with state‐of‐the‐art high‐order numerical schemes.Copyright © 2012 John Wiley & Sons, Ltd.     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

Closure modeling in bridging regions of variable-resolution (VR) turbulence computations

Optimal utilization of computational resources mandates spatio-temporal variation in resolution for computing complex engineering flows. Closure modeling in regions bridging between different resolutions is rendered difficult due to changing interactions between resolved and unresolved fields. We develop a closure model for the bridging region based on energy conservation principles. Then we proceed to provide a proof of concept in decaying isotropic turbulence with temporally varying resolution. The simplicity of the flow permits a thorough examination of various aspects of the proposed closure not feasible in more complex flows. The results demonstrate the potential promise of the approach, but more validation studies need to be performed.While the present development is in the context of partially averaged Navier–Stokes (PANS) method, the closure principle should apply for other variable-resolution (VR) approaches.     

Energy flow along the medium-induced parton cascade

We discuss the dynamics of parton cascades that develop in dense QCD matter, and contrast their properties with those of similar cascades of gluon radiation in vacuum. We argue that such cascades belong to two distinct classes that are characterized respectively by an increasing or a constant (or decreasing) branching rate along the cascade. In the former class, of which the BDMPS, medium-induced, cascade constitutes a typical example, it takes a finite time to transport a finite amount of energy to very soft quanta, while this time is essentially infinite in the latter case, to which the DGLAP cascade belongs. The medium induced cascade is accompanied by a constant flow of energy towards arbitrary soft modes, leading eventually to the accumulation of the initial energy of the leading particle at zero energy. It also exhibits scaling properties akin to wave turbulence. These properties do not show up in the cascade that develops in vacuum. There, the energy accumulates in the spectrum at smaller and smaller energy as the cascade develops, but the energy never flows all the way down to zero energy. Our analysis suggests that the way the energy is shared among the offsprings of a splitting gluon has little impact on the qualitative properties of the cascades, provided the kernel that governs the splittings is not too singular.     

湍流大气传输中跟踪抖动对远场影响的仿真研究

分析了跟踪抖动对湍流大气传输远场光斑的影响。基于麦克斯韦电磁场理论,采用大气相干长度对大气湍流进行描述,推导了发射光束因跟踪抖动导致光轴偏离的远场表达式。在此基础上,利用相位屏法模拟抖动引起的倾斜相位和大气折射率起伏引起的相位调制,并采用低频补偿的功率谱反演法对传输过程进行了数值仿真。分析了不同跟踪抖动、湍流强度条件下远场光斑质心脱靶量的变化,以及不同尺寸模拟目标的回波概率。分析结果表明,在传输距离为10 km时,强湍流造成的远场光斑脱靶量可达几十μrad;当跟踪抖动较大时,湍流强弱对脱靶量影响差别很小。最后,对一定尺寸的模拟目标,从探测回波概率的角度给出了发射系统跟踪抖动量的控制范围。