Interface‐fitted moving mesh method for axisymmetric two‐phase flow in microchannels

A boundary‐fitted moving mesh scheme is presented for the simulation of two‐phase flow in two‐dimensional and axisymmetric geometries. The incompressible Navier‐Stokes equations are solved using the finite element method, and the mini element is used to satisfy the inf‐sup condition. The interface between the phases is represented explicitly by an interface adapted mesh, thus allowing a sharp transition of the fluid properties. Surface tension is modelled as a volume force and is discretized in a consistent manner, thus allowing to obtain exact equilibrium (up to rounding errors) with the pressure gradient. This is demonstrated for a spherical droplet moving in a constant flow field. The curvature of the interface, required for the surface tension term, is efficiently computed with simple but very accurate geometric formulas. An adaptive moving mesh technique, where smoothing mesh velocities and remeshing are used to preserve the mesh quality, is developed and presented. Mesh refinement strategies, allowing tailoring of the refinement of the computational mesh, are also discussed. Accuracy and robustness of the present method are demonstrated on several validation test cases. The method is developed with the prospect of being applied to microfluidic flows and the simulation of microchannel evaporators used for electronics cooling. Therefore, the simulation results for the flow of a bubble in a microchannel are presented and compared to experimental data.     

薄板结构加筋布局优化设计方法研究

提出一种两步优化策略，以加筋板结构的固有频率最大化为目标函数，以结构所受外载荷作用的最大静变形为约束条件，开展薄板结构加筋构件的布局优化设计研究。为了降低加筋布局优化的难度，提高优化设计的效率，将加筋等效为一系列弹性铰（点）支撑，以便快速获得加筋横向移动的灵敏度信息。在基本不改变结构重量的情形下，通过合理布局加筋位置，能显著改善结构的刚度分布，提高结构的整体承载能力。随后，小幅调整加筋的截面尺寸，以满足对结构最大变形的设计要求。最后，用两个算例验证了所提优化方法的可行性和有效性。     

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.     

Existence of strong minimizers for the Griffith static fracture model in dimension two

We consider the Griffith fracture model in two spatial dimensions, and prove existence of strong minimizers, with closed jump set and continuously differentiable deformation fields. One key ingredient, which is the object of the present paper, is a generalization to the vectorial situation of the decay estimate by De Giorgi, Carriero, and Leaci. This is based on replacing the coarea formula by a method to approximate $SB{D}^p$ functions with small jump set by Sobolev functions, and is restricted to two dimensions. The other two ingredients will appear in companion papers and consist respectively in regularity results for vectorial elliptic problems of the elasticity type and in a method to approximate in energy $GSB{D}^p$ functions by $SB{V}^p$ ones.     

MECHANICAL ANALYSIS OF TOROFLUX

流动环是一款魔术玩具,具有丰富的力学行为。制作流动环时应用了预应力储存技术,使得环内部储存了弹性弯曲变形能。一经外界扰动,流动环会迅速张开,进入一种新的稳定平衡状态;若给张开的流体环输入某种形式的能量,环会合拢,恢复原有的稳定平衡状态。本文从能量的角度对流动环张开与合拢等现象进行力学分析。     

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.     

Characteristic behaviour of plasma fluctuations inside a plasma bubble in presence of a magnetic field due to the formation of a potential well

Experiments were carried out to observe the effect of a magnetic field and grid biasing voltage in presence of a plasma bubble in a magnetized, filamentary discharge plasma system. A spherical mesh grid of 80% optical transparency was negatively biased and introduced into the plasma for creating a plasma bubble. Diagnostics via an electrical Langmuir probe and a hot emissive probe were extensively used for scanning the plasma bubble. Plasma floating potential fluctuations were measured at three different positions of the plasma bubble. The instability in the pattern showed the dynamic transition from periodic to chaotic for increasing magnetic fields. Time scale analysis using continuous wavelet transform was carried out to identify the presence of non‐linearity from the contour plots. The mechanisms of the low‐frequency instabilities along with the transition to chaos could be qualitatively explained. Non‐linear techniques such as fast Fourier transform, phase space plot, and recurrence plot were used to explore the dynamics of the system appearing during plasma fluctuations. In order to demonstrate the observed chaotic phenomena in this study, characteristics of chaos such as the Lyapunov exponent were obtained from experimental time series data. The experimentally observed potential structure is confirmed with numerical analysis based on fluid hydrodynamics.     

Dynamic Wetting Properties of Mesh Substrates with Tunable Water Adhesion

In nature, wetting phenomena are present nearly everywhere and are a source of inspiration for liquid transportation. A good understanding of the underlying dynamic phenomena that governs wettability is therefore extremely important for researchers involved in bio-inspired surfaces. Herein, we study the adhesive behavior with water of mesh substrates modified with structured copolymers in order to tune the surfaces from parahydrophobic states (high water adhesion) to superhydrophobic states (low water adhesion). Using the ejection test method (ETM), a new technique that consists of the ejection of water droplets deposited onto a substrate with the aid of a catapult system, we experimentally demonstrate that the elasticity of the mesh substrate can be exploited for efficient vertical actuation of droplets.     

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.     

Finite element modelling of cold drawing for high-precision tubes

Cold tube drawing is a metal forming process that allows manufacturers to produce high-precision tubes. The dimensions of the tube are reduced by pulling it through a conical converging die with or without inner tool. In this study, finite element modelling has been used to give a better understanding of the process.This paper presents a model that predicts the final dimensions of the tube with very high accuracy. It is validated thanks to experimental tests. Moreover, five studies are performed with this model, such as investigating the influence of the die angle on the drawing force or the influence of relative thickness on tube deformation.