An hp‐adaptive unstructured finite volume solver for compressible flows

The idea of hp‐adaptation, which has originally been developed for compact schemes (such as finite element methods), suggests an adaptation scheme using a mixture of mesh refinement and order enrichment based on the smoothness of the solution to obtain an accurate solution efficiently. In this paper, we develop an hp‐adaptation framework for unstructured finite volume methods using residual‐based and adjoint‐based error indicators. For the residual‐based error indicator, we use a higher‐order discrete operator to estimate the truncation error, whereas this estimate is weighted by the solution of the discrete adjoint problem for an output of interest to form the adaptation indicator for adjoint‐based adaptations. We perform our adaptation by local subdivision of cells with nonconforming interfaces allowed and local reconstruction of higher‐order polynomials for solution approximations. We present our results for two‐dimensional compressible flow problems including subsonic inviscid, transonic inviscid, and subsonic laminar flow around the NACA 0012 airfoil and also turbulent flow over a flat plate. Our numerical results suggest the efficiency and accuracy advantages of adjoint‐based hp‐adaptations over uniform refinement and also over residual‐based adaptation for flows with and without singularities.     

Mesh adaptation using C1 interpolation of the solution in an unstructured finite volume solver

The goal of this paper is to show the effectiveness of a newly developed estimate of the truncation error calculated based on C¹ interpolation of the solution weighted by the adjoint solution as the adaptation indicator for an unstructured finite volume solver. We will show that adjoint‐based mesh adaptation based on the corrected functional using the new developed truncation error estimate is capable of adapting the mesh to improve the accuracy of the functional and the convergence rate. Both discrete and continuous adjoint solutions are used for adaptation. Results are significantly better with new truncation error estimate than with previously used estimates.     

The Method of R-Functions in the Solution of Elastic Problems on the Basis of Reissner's Mixed Variational Principle

A method is presented for solving boundary-value elastic problems on the basis of the variational–structural method of R-functions and Reissner's mixed variational principle. A mathematical formulation is given to problems on the deformation of elastic bodies under mixed boundary conditions and bodies interacting with smooth rigid dies. Solutions satisfying all the boundary conditions are proposed. For undetermined components of these solutions, the resolving equations are derived and their properties are studied. A posteriori estimation of numerical solutions is made. As examples, solutions are found to a problem on the stress–strain state of a short cylinder and to a contact problem on a cylinder interacting with a smooth die. A numerical method of solving such problems is analyzed for convergence, and the accuracy of the solutions is estimated.     

Simulation of two‐phase flow with moving immersed boundaries

A two‐dimensional multi‐phase model for immiscible binary fluid flow including moving immersed objects is presented. The fluid motion is described by the incompressible Navier–Stokes equation coupled with a phase‐field model based on van der Waals' free energy density and the Cahn–Hilliard equation. A new phase‐field boundary condition was implemented with minimization of the free energy in a direct way, to specifically improve the physical behavior of the contact line dynamics for moving immersed objects. Numerical stability and execution time were significantly improved by the use of the new boundary condition. Convergence toward the analytical solution was demonstrated for equilibrium contact angle, the Lucas–Washburn theory and Stefan's problem. The proposed model may be used for multi‐phase flow problems with moving boundaries of complex geometry, such as the penetration of fluid into a deformable, porous medium. Copyright © 2010 John Wiley & Sons, Ltd.     

Adjoint‐based airfoil optimization with discretization error control

In this article, we develop a new airfoil shape optimization algorithm based on higher‐order adaptive DG methods with control of the discretization error. Each flow solution in the optimization loop is computed on a sequence of goal‐oriented h‐refined or hp‐refined meshes until the error estimation of the discretization error in a flow‐related target quantity (including the drag and lift coefficients) is below a prescribed tolerance. Discrete adjoint solutions are computed and employed for the multi‐target error estimation and adaptive mesh refinement. Furthermore, discrete adjoint solutions are employed for evaluating the gradients of the objective function used in the CGs optimization algorithm. Furthermore, an extension of the adjoint‐based gradient evaluation to the case of target lift flow computations is employed. The proposed algorithm is demonstrated on an inviscid transonic flow around the RAE2822, where the shape is optimized to minimize the drag at a given constant lift and airfoil thickness. The effect of the accuracy of the underlying flow solutions on the quality of the optimized airfoil shapes is investigated. Copyright © 2014 John Wiley & Sons, Ltd.     

Application of the angular superposition method to the contact problem on the compression of an elastic cylinder

The angular superposition method is used to construct an approximate solution of the contact problem on the compression of an elastic cylinder by two rigid plates. The solution thus obtained has a closed-form analytic expression and can be used in the entire domain of the cylinder cross-section. We analyze the absolute error, which takes the largest value near the points of contact between the plates and the cylinder, where the boundary conditions are discontinuous. According to the von Mises criterion, when moving into the depth of the cylinder from the contact site along the symmetry axis, the second invariant J ₂ of the stress deviator tensor first decreases and then, after attaining a minimum, increases and attains the largest value at a small depth, which agrees with Johnson’s photoelastic experiments and Dinnik’s computations. We present the graphs of the displacement and normal stress distributions over the contact site, the dependence of the compressing force on the displacements of rigid plates, and the dependence of the invariant J ₂ on the coordinate along the symmetry axis. If 640 computation points are chosen on the cylinder boundary and the Hertz law for the normal pressure on the contact site is used, then the error in the approximate solution near the endpoint of the contact site is approximately 55%, and if the proposed two-parameter normal law is used, then the error is of the order of 4%. On the free lateral surface of the cylinder boundary, we find the critical pointM*, which separates the cylinder contraction and extension parts.The contact problems are the most difficult problems, and their solution is complicated by the discontinuous boundary conditions [1–5]. In [6], the contact problem is solved by the Fourier method, which can be used only for bodies of classical shapes. In such cases, the problem can be reduced to solving coupled integral equations [7]. The interaction between the bandage and a cylindrical body is considered in [2, 6, 7]. In [8], the possibility of using the finite element method is investigated in the case of contact problems for a differential wheel with roughness of the contacting surfaces taken into account. In [9, 10], the method of homogeneous solutions is used to consider contact problems for a finite-dimensional elastic cylinder loaded on its end surfaces. Note that only error estimates are given in the literature cited above; the absolute error over the entire domain of the elastic body is not studied, although this is one of the important characteristics of the obtained approximate solution. A sufficiently complete survey of the literature in the field of contact interactions of elastic bodies is given in [3–5].In what follows, we propose to solve contact problems by the angular superposition method [11]. This method can be used for bodies of nonclassical shapes, which can be multiply connected, and the friction on the contact site can be taken into account. In the present paper, as a first example of applied character, we show how this method can be used in the simplest case. The multiple connectedness and the curvilinearity of the shape of the body, as well as taking into account the friction on the boundary, do not create new essential difficulties in this method.     

Domain versus boundary computation of flow sensitivities with the continuous adjoint method for aerodynamic shape optimization problems

This paper considers the computation of flow sensitivities that arise in the context of design optimization. The scheme is based on the solution of a continuous adjoint problem, for which two complementary, although analytically equivalent, approaches have been routinely used for some time now, yielding expressions for the sensitivities that contain, respectively, boundary and domain integrals. These concepts are clarified in a unified framework and their equivalence at the continuous level is demonstrated through appropriate algebraic manipulations. Equivalence at the discrete level is assessed through numerical testing for various aerodynamic shape‐optimization problems. Copyright © 2011 John Wiley & Sons, Ltd.     

The thermoelastic Hertzian contact problem

A numerical solution is obtained for the steady-state thermoelastic contact problem in which heat is conducted between two elastic bodies of dissimilar materials at different temperatures with arbitrary quadratic profiles. Thermoelastic deformation causes the initially elliptical contact area to be reduced in size and to become more nearly circular as the temperature difference is increased. There is also a small but identifiable deviation from exact ellipticity at intermediate temperature differences. An approximate analytical solution is obtained, based on approximating the contact area by an ellipse.     

A hybrid immersed‐boundary and multi‐block lattice Boltzmann method for simulating fluid and moving‐boundaries interactions

A numerical method is developed for modelling the interactions between incompressible viscous fluid and moving boundaries. The principle of this method is introducing the immersed‐boundary concept in the framework of the lattice Boltzmann method, and improving the accuracy and efficiency of the simulation by refining the mesh near moving boundaries. Besides elastic boundary with a constitutive law, the method can also efficiently simulate solid moving‐boundary interacting with fluid by employing the direct forcing technique. The method is validated by the simulations of flow past a circular cylinder, two cylinders moving with respect to each other and flow around a hovering wing. The versatility of the method is demonstrated by the numerical studies including elastic filament flapping in the wake of a cylinder and fish‐like bodies swimming in quiescent fluid. Copyright © 2006 John Wiley & Sons, Ltd.     

Simultaneous solution of the Navier-Stokes and elastic membrane equations by a finite element method

Fluid flow through a significantly compressed elastic tube occurs in a variety of physiological situations. Laboratory experiments investigating such flows through finite lengths of tube mounted between rigid supports have demonstrated that the system is one of great dynamical complexity, displaying a rich variety of self-excited oscillations. The physical mechanisms responsible for the onset of such oscillations are not yet fully understood, but simplified models indicate that energy loss by flow separation, variation in longitudinal wall tension and propagation of fluid elastic pressure waves may all be important. Direct numerical solution of the highly non-linear equations governing even the most simplified two-dimensional models aimed at capturing these basic features requires that both the flow field and the domain shape be determined as part of the solution, since neither is known a priori. To accomplish this, previous algorithms have decoupled the solid and fluid mechanics, solving for each separately and converging iteratively on a solution which satisfies both. This paper describes a finite element technique which solves the incompressible Navier-Stokes equatikons simultaneously with the elastic membrane equations on the flexible boundary. The elastic boundary position is parametized in terms of distances along spines in a manner similar to that which has been used successfully in studies of viscous free surface flows, but here the membrane curvature equation rather than the kinematic boundary condition of vanishing normal velocity is used to determine these diatances and the membrane tension varies with the shear stresses exerted on it by the fluid motions. Bothy the grid and the spine positions adjust in response to membrane deformation, and the coupled fluid and elastic equations are solved by a Newton-Raphson scheme which displays quadratic convergence down to low membrane tensions and extreme states of collapse. Solutions to the steady problem are discussed, along with an indication of how the time-dependent problem might be approached.     

Remarks on the numerical solution of the adjoint quasi‐one‐dimensional Euler equations

We examine the numerical solution of the adjoint quasi‐one‐dimensional Euler equations with a central‐difference finite volume scheme with Jameson‐Schmidt‐Turkel (JST) dissipation, for both the continuous and discrete approaches. First, the complete formulations and discretization of the quasi‐one‐dimensional Euler equations and the continuous adjoint equation and its counterpart, the discrete adjoint equation, are reviewed. The differences between the continuous and discrete boundary conditions are also explored. Second, numerical testing is carried out on a symmetric converging–diverging duct under subsonic flow conditions. This analysis reveals that the discrete adjoint scheme, while being manifestly less accurate than the continuous approach, gives nevertheless more accurate flow sensitivities. Copyright © 2011 John Wiley & Sons, Ltd.     

大变形中摩擦接触问题的数值模拟及应用

大变形的摩擦接触是复杂的非线性问题 ,本文介绍了一种处理摩擦接触问题的数值方法。采用接触单元技术模拟接触界面 ,基于弹塑性理论形式的非经典Coulomb摩擦定律及罚函数方法建立了摩擦接触的增量本构关系。结合大变形的增量分析格式给出了积分摩擦接触本构方程的回映方法。这种处理摩擦接触问题的方法计算简单、使用方便。给出的计算实例及应用实例说明了方法的精度与稳定性     

Boundary element method for moving and rolling contact of 2D elastic bodies with defects

A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes. The project supported by the National Natural Science Foundation of China (19772025)     

Control of hydrothermal waves in a thermocapillary flow using a gradient‐based control strategy

The problem of controlling the hydrothermal waves in a thermocapillary flow is addressed using a gradient‐based control strategy. The state equations are the two‐dimensional unsteady incompressible Navier–Stokes and energy equations under the Boussinesq approximation. The modeled problem is the ‘open boat’ process of crystal growth, the flow which is driven by Marangoni and buoyancy effects. The control is a spatially and temporally varying heat flux boundary condition at the free surface. The control that minimizes the hydrothermal waves is found using a conjugate gradient method, where the gradient of the objective function with respect to the control variables is obtained from solving a set of adjoint equations. The effectiveness of choices of the parameters governing the control algorithm is examined. Almost complete suppression of the hydrothermal waves is obtained for certain choices of the parameters governing the control algorithm. The numerical issues involved with finding the control using the optimizer are discussed, and the features of the resulting control are analyzed with the goal of understanding how it affects the flow.Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical control of two‐dimensional shock waves in dual solution domain by instant temperature disturbances

The search for the temperature disturbance causing transition between regular and Mach reflections in the dual solution domain is addressed in an optimization statement. The gradient of the discrepancy between the current and target flow fields was calculated using adjoint equations. The control was determined by gradient‐based optimization. The flow field simulation is verified via a posteriori error estimates using the solution of an additional adjoint problem. Copyright © 2012 John Wiley & Sons, Ltd.     

Viscous flows in corner regions: Singularities and hidden eigensolutions

Numerical issues arising in computations of viscous flows in corners formed by a liquid–fluid free surface and a solid boundary are considered. It is shown that on the solid a Dirichlet boundary condition, which removes multivaluedness of velocity in the ‘moving contact‐line problem’ and gives rise to a logarithmic singularity of pressure, requires a certain modification of the standard finite‐element method. This modification appears to be insufficient above a certain critical value of the corner angle where the numerical solution becomes mesh‐dependent. As shown, this is due to an eigensolution, which exists for all angles and becomes dominant for the supercritical ones. A method of incorporating the eigensolution into the numerical method is described that makes numerical results mesh‐independent again. Some implications of the unavoidable finiteness of the mesh size in practical applications of the finite‐element method in the context of the present problem are discussed. Copyright © 2009 John Wiley & Sons, Ltd.     

移动接触下求解二维闭合裂纹的双重迭代边界元法

使用子域边界元法对受移动接触弹性体作用下的二维闭合裂纹问题进行了数值计算。由于两弹性体的接触界面和裂纹表面的接触范围的大小和接触状态事先是未知的 ,对此 ,在两个接触表面同时采用迭代的方法进行了求解。在裂纹的每个裂尖上都采用了四分之一的奇异单元以保证裂尖位移场和应力场奇异性的满足。用我们编制的二维裂纹问题程序对一些中心裂纹问题进行了计算 ,计算结果与经典断裂力学的理论值比较吻合。在无摩擦的条件下 ,对一些具有不同角度且受移动接触弹性体作用下的闭合裂纹问题进行了数值计算 ,得到了一些耦合作用下的应力强度因子的计算结果     

A simple and efficient error analysis for multi‐step solution of the Navier–Stokes equations

A simple error analysis is used within the context of segregated finite element solution scheme to solve incompressible fluid flow. An error indicator is defined based on the difference between a numerical solution on an original mesh and an approximated solution on a related mesh. This error indicator is based on satisfying the steady‐state momentum equations. The advantages of this error indicator are, simplicity of implementation (post‐processing step), ability to show regions of high and/or low error, and as the indicator approaches zero the solution approaches convergence. Two examples are chosen for solution; first, the lid‐driven cavity problem, followed by the solution of flow over a backward facing step. The solutions are compared to previously published data for validation purposes. It is shown that this rather simple error estimate, when used as a re‐meshing guide, can be very effective in obtaining accurate numerical solutions. Copyright © 2002 John Wiley & Sons, Ltd.     

The solution of two‐dimensional free‐surface problems using automatic mesh generation

A new method is described for the iterative solution of two‐dimensional free‐surface problems, with arbitrary initial geometries, in which the interior of the domain is represented by an unstructured, triangular Eulerian mesh and the free surface is represented directly by the piecewise‐quadratic edges of the isoparametric quadratic‐velocity, linear‐pressure Taylor–Hood elements. At each time step, the motion of the free surface is computed explicitly using the current velocity field and, once the new free‐surface location has been found, the interior nodes of the mesh are repositioned using a continuous deformation model that preserves the original connectivity. In the event that the interior of the domain must be completely remeshed, a standard Delaunay triangulation algorithm is used, which leaves the initial boundary discretisation unchanged. The algorithm is validated via the benchmark viscous flow problem of the coalescence of two infinite cylinders of equal radius, in which the motion is due entirely to the action of capillary forces on the free surface. This problem has been selected for a variety of reasons: the initial and final (steady state) geometries differ considerably; in the passage from the former to the latter, large free‐surface curvatures—requiring accurate modelling—are encountered; an analytical solution is known for the location of the free surface; there exists a large body of literature on alternative numerical simulations. A novel feature of the present work is its geometric generality and robustness; it does not require a priori knowledge of either the evolving domain geometry or the solution contained therein. Copyright © 1999 John Wiley & Sons, Ltd.     

预紧方式对弹流润滑下角接触球轴承内部力学特性的影响

油膜厚度预测在评估弹流润滑(EHL)下角接触球轴承的性能和耐久性方面发挥着重要的作用. 耦合拟静力学理论和自旋下椭圆接触弹流模型,以干接触角接触球轴承拟静力学分析方法为基础,建立了定压和定位预紧方式下考虑弹流润滑和钢球自旋运动的角接触球轴承的拟静力学分析模型. 采用快速傅里叶变换(FFT)计算椭圆接触的弹性变形,运用Gauss-Seidel迭代方法求解Reynolds方程,得到自旋弹流模型的完全数值解,将其代入轴承拟静力学模型中迭代,得到轴承内部接触载荷、三维接触压力及三维膜厚分布. 对采用不同预紧方式的SKF7210型角接触球轴承进行分析,结果表明:富油润滑下,当轴承转速从0增大到15 000 r/min时,定压预紧时内圈轴向位移减小17.83%,而定位预紧时内圈承受的轴向载荷增大23.17%;定压预紧方式下球与内外滚道间膜厚均略大于定位预紧. 此外,不同预紧方式下,外圈上的中心膜厚大于内圈10%. 与干接触相比,定压下考虑弹流润滑内圈上接触载荷略大0.64%.