Power-of-two sampling in redundancy systems: The impact of assignment constraints

A classical sampling strategy for load balancing policies is power-of-two, where any server pair is sampled with equal probability. This does not cover practical settings with assignment constraints which force non-uniform sampling. While intuition suggests that non-uniform sampling adversely impacts performance, this was only supported through simulations, and rigorous statements have remained elusive. Building on product-form distributions for redundancy systems, we prove the stochastic dominance of uniform sampling for a four-server system as well as arbitrary-size systems in light traffic.     

Negative thermal expansion in NbF_3 and NbOF_2:A comparative theoretical study

Thermal expansion control is always an obstructive factor and challenging in high precision engineering field. Here,the negative thermal expansion of Nb F_3 and Nb OF_2 was predicted by first-principles calculation with density functional theory and the quasi-harmonic approximation(QHA). We studied the total charge density, thermal vibration, and lattice dynamic to investigate the thermal expansion mechanism. We found that the presence of O induced the relatively strong covalent bond in Nb OF_2, thus weakening the transverse vibration of F and O in Nb OF_2, compared with the case of Nb F_3.In this study, we proposed a way to tailor negative thermal expansion of metal fluorides by introducing the oxygen atoms.The present work not only predicts two NTE compounds, but also provides an insight on thermal expansion control by designing chemical bond type.     

A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

体积应变对立方钛酸铅电子结构和光学性质的影响

本文利用第一性原理方法计算并分析了体积应变(-11%~11%)对立方顺电相PbTiO₃的结构、稳定性、电子结构和光学性质的影响。研究发现体积应变后PbTiO₃形成焓增大,稳定性下降,其中压应变对其稳定性的影响比拉应变大。当受到拉伸应变时,立方PbTiO₃由直接带隙半导体变为间接带隙半导体,且带隙随应变增大呈先增大后降低的趋势。在发生压应变时,从复介电函数、复折射率及吸收系数的分析结果可知,在自然光照下PbTiO₃的光吸收能力仅在个别波段有所增大,但总体呈减弱趋势,当产生拉伸应变时,介电峰、吸收峰红移,表明PbTiO₃在可见光范围内光吸收能力增强,并且当应变增大到11%时,PbTiO₃的吸收能力远高于本征立方相。     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

A cost‐effective curvature calculation approach for interfacial flows on unstructured meshes

We present a simple and cost‐effective curvature calculation approach for simulations of interfacial flows on structured and unstructured grids. The interface is defined using volume fractions, and the interface curvature is obtained as a function of the gradients of volume fractions. The gradient computation is based on a recently proposed gradient recovery method that mimicks the least squares approach without the need to solve a system of equations and is quite easy to implement on arbitrary polygonal meshes. The resulting interface curvature is used in a continuum surface force formulation within the framework of a well‐balanced finite‐volume algorithm to simulate multiphase flows dominated by surface tension. We show that the proposed curvature calculation is at least as accurate as some of the existing approaches on unstructured meshes while being straightforward to implement on any mesh topology. Numerical investigations also show that spurious currents in stationary problems that are dependent on the curvature calculation methodology are also acceptably low using the proposed approach. Studies on capillary waves and rising bubbles in viscous flows lend credence to the ability of the proposed method as an inexpensive, robust, and reasonably accurate approach for curvature calculation and numerical simulation of multiphase flows.     

An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.     

以时温等效方法研究尼龙1010应力松弛行为及使用寿命预测

以尼龙材料的应力松弛行为作为研究对象, 考察初始应变为1.0%, 2.8%和5.1%的尼龙1010样品在温度区间293353 K的松弛曲线, 采用时间-温度等效叠加方法得到了松弛模量主曲线, 计算出叠加过程中的表观活化能、 松弛过程中的活化体积和应力辅助功. 结果表明, 整个松弛过程中的表观活化能和应力辅助功表现出相同的变化趋势, 体现出松弛过程中克服运动单元位垒的过程. 当293323 K区间的松弛曲线叠加时, 随着初始应变的增加, 表观活化能和应力辅助功均逐渐降低, 有助于聚合物内部的运动单元越过能垒发生松弛, 与松弛过程中的应力辅助热活化理论相一致; 当333353 K区间的松弛曲线叠加时, 不同初始应变样品的表观活化能均为260 kJ/mol, 应力辅助功均为60 MPa·nm³, 说明松弛过程中克服运动单元的能垒与应力作用无关. 根据松弛主曲线, 计算出了尼龙1010在1.0%, 2.8%和5.1% 3种形变下, 长时间范围内应力衰减与时间的关系, 为预测实际使用过程中的应力松弛行为提供了依据.     

基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.     

An adaptive multilevel wavelet framework for scale‐selective WENO reconstruction schemes

We put forth a dynamic computing framework for scale‐selective adaptation of weighted essential nonoscillatory (WENO) schemes for the simulation of hyperbolic conservation laws exhibiting strong discontinuities. A multilevel wavelet‐based multiresolution procedure, embedded in a conservative finite volume formulation, is used for a twofold purpose. (i) a dynamic grid adaptation of the solution field for redistributing grid points optimally (in some sense) according to the underlying flow structures, and (ii) a dynamic minimization of the in built artificial dissipation of WENO schemes. Taking advantage of the structure detection properties of this multiresolution algorithm, the nonlinear weights of the conventional WENO implementation are selectively modified to ensure lower dissipation in smoother areas. This modification is implemented through a linear transition from the fifth‐order upwind stencil at the coarsest regions of the adaptive grid to a fully nonlinear fifth‐order WENO scheme at areas of high irregularity. Therefore, our computing algorithm consists of a dynamic grid adaptation strategy, a scale‐selective state reconstruction, a conservative flux calculation, and a total variation diminishing Runge‐Kutta scheme for time advancement. Results are presented for canonical examples drawn from the inviscid Burgers, shallow water, Euler, and magnetohydrodynamic equations. Our findings represent a novel direction for providing a scale‐selective dissipation process without a compromise on shock capturing behavior for conservation laws, which would be a strong contender for dynamic implicit large eddy simulation approaches.