Inequalities relating to Lp-version of Petty＇s conjectured projection inequality

Petty＇s conjectured projection inequality is a famous open problem in the theory of convex bodies. In this paper, it is shown that an inequality relating to Lp-version of the Petty＇s conjectured projection inequality is developed by using the notions of the Lp-mixed volume and the Lp-dual mixed volume, the relation of the Lp-projection body and the geometric body Г-pK, the Bourgain-Milman inequality and the Lp-Bnsemann-Petty inequality. In addition, for each origin-symmetric convex body, by applying the Jensen inequality and the monotonicity of the geometric body Г-pK, the reverses of Lp-version of the Petty＇s conjectured projection inequality and the Lp-Petty projection inequality are given, respectively.     

Inequalities relating to L_p-version of Petty's conjectured projection inequality

Petty's conjectured projection inequality is a famous open problem in the theory of convex bodies.In this paper,it is shown that an inequality relating to L_p-version of the Petty's conjectured projection inequality is developed by using the notions of the L_p-mixed volume and the L_p-dual mixed volume,the relation of the L_p-projection body and the geometric bodyΓ_(-p)K,the Bourgain-Milman inequality and the L_p-Busemann- Petty inequality.In addition,for each origin-symmetric convex body,by applying the Jeusen inequality and the monotonicity of the geometric bodyΓ_(-p)K,the reverses of L_p-version of the Perry's conjectured projection inequality and the L_p-Petty projection inequality are given,respectively.     

线性区间有限元静力控制方程的组合解法

区间有限元的静力控制方程常被归结为区间方程组来求解。但实际上两者并不等价。本文根据不确定结构有限元分析的力学背景，直接从问题的基本参量的不确定性出发，将基本区间参量的边界组合与求解区间方程组的有关解法相结合，提出了线性区间有限元静力控制方程的两种组合解法－参量边界全组合法和组合迭代法。可以以较小的计算量获得或逼近位移和应力区间的准确界限。且不受基本参量变化范围的限制。算例分析表明文中方法是实用和可行的。     

High-efficiency improved symmetric successive over-relaxation preconditioned conjugate gradient method for solving large-scale finite element linear equations

Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering （CAE） and computer aided manufacturing （CAM）. This paper presents a high-efficiency improved symmetric successive over-relaxation （ISSOR） preconditioned conjugate gradient （PCG） method, which maintains lelism consistent with the original form. Ideally, the by 50% as compared with the original algorithm. the convergence and inherent paralcomputation can It is suitable for be reduced nearly high-performance computing with its inherent basic high-efficiency operations. By comparing with the numerical results, it is shown that the proposed method has the best performance.     

High-effciency improved symmetric successive over-relaxation preconditioned conjugate gradient method for solving large-scale finite element linear equations

Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering (CAE) and computer aided ...     

A linear solver based on algebraic multigrid and defect correction for the solution of adjoint RANS equations

A new solution approach for discrete adjoint Reynolds‐averaged Navier–Stokes equations is suggested. The approach is based on a combination of algebraic multigrid and defect correction. The efficient interplay of these two methods results in a fast and robust solution technique. Algebraic multigrid deals very well with unstructured grids, and defect correction completes it in such a way that a required second‐order accuracy of the solution is achieved. The performance of the suggested solution approach is demonstrated on a number of representative benchmarks.Copyright © 2014 John Wiley & Sons, Ltd.     

New nonconforming finite element method for solving transient Naiver-Stokes equations

For transient Naiver-Stokes problems, a stabilized nonconforming finite element method is presented, focusing on two pairs inf-sup unstable finite element spaces, i.e., pNC/pNC triangular and pNQ/pNQ quadrilateral finite element spaces. The semi- and full-discrete schemes of the stabilized method are studied based on the pressure projection and a variational multi-scale method. It has some attractive features： avoiding higher-order derivatives and edge-based data structures, adding a discrete velocity term only on the fine scale, being effective for high Reynolds number fluid flows, and avoiding increased computation cost. For the full-discrete scheme, it has second-order estimations of time and is unconditionally stable. The presented numerical results agree well with the theoretical results.     

线性蠕变问题中两种摄动半解析计算格式

将线性蠕变的控制方程做摄动展开,用加权余量法和二类变量的变分原理,引入有限元线和Ｈａｍｉｌｔｏｎ混合状态元技术,建立了两种求解缄性蠕变问题的摄动增解析计算格式。     

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.     

Existence and comparison results of solutions for nonlinear random integral and differential equations relative to weak topology

In this paper, we prove several existence theorems of random solutions to nonlinear random Volterra integral equations under the weak topology of Banach spaces. Then, as applications, we obtain the existence theorems of weak random solutions to random differential equations. Existence of extremal random solutions and a random comparison theorem for these random equations are also obtained. Our theorems improve and extend the corresponding results in [4,5,10,11,12]. Projects supported by the Science Fund of the Chinese Academy of Sciences.     

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.     

非定常不可压N-S方程的最小二乘算子分裂有限元法数值求解

采用最小二乘算子分裂有限元法求解非定常不可压N-S(Navier-Stokes)方程,即在每个时间层上采用算子分裂法将N-S方程分裂成扩散项和对流项,这样既能考虑对流占优特点又能顾及方程的扩散性质。扩散项是一个抛物型方程,时间离散采用向后差分格式,空间离散采用标准Galerkin有限元法。对流项的时间项采用后向差分格式,非线性部分用牛顿法进行线性化处理,再用最小二乘有限元法进行空间离散,得到对称正定的代数方程组系数矩阵。采用Re=1000的方腔流对该算法的有效性进行检验,表明其具有较高的精度,能够很好地捕捉流场中的涡结构。同时,对圆柱层流绕流进行了数值研究,通过流线图、压力场、阻力系数、升力系数及斯特劳哈数等结果的分析与对比,表明本文算法对于模拟圆柱层流绕流是准确和可靠的。     

Quasi-optimal complexity of adaptive finite element method for linear elasticity problems in two dimensions

This paper introduces an adaptive finite element method (AFEM) using the newest vertex bisection and marking exclusively according to the error estimator without special treatment of oscillation. By the combination of the global lower bound and the localized upper bound of the posteriori error estimator, perturbation of oscillation, and cardinality of the marked element set, it is proved that the AFEM is quasi-optimal for linear elasticity problems in two dimensions, and this conclusion is verified by the numerical examples.     

稀薄流非线性模型方程离散速度坐标法有限差分解

从一般非线性Bo ltzm ann方程出发,发展并实现了一套适于大范围K nudsen数稀薄流问题数值模拟的统一算法。采用BGK模型和Shakov模型近似碰撞项,进而引入两个二速度无量纲简化分布函数,通过关于分子速度第三分量取矩积分,将三速度单一模型方程变换为二速度微分方程组。基于G auss-H erm ite积分公式和正交多项式G auss积分公式,借助离散速度坐标法消除简化模型方程对分子速度空间的连续依赖性,从相空间到物理空间得到一组带源项双曲守恒离散方程,并给出其显式和隐式二阶迎风TVD有限差分解。以二维圆柱A r气体超声速绕流算例,验证了数值算法的有效性,比较分析了漫反射和镜面反射两种气体分子壁面反射模型的计算结果。     

A discontinuous Galerkin method/HLLC solver for the Euler equations

This paper proposes a fully three‐dimensional non‐linear Euler methodology for solving aerodynamic and acoustic problems in the presence of strong shocks and rarefactions. It uses a discontinuous Galerkin method (DGM) within the element, and a Riemann solver (HLLC) at the boundaries to propagate rarefactions while preserving the entropy condition and capturing shocks with no spurious oscillations. This approach is thought to marry the best aspects of finite element and finite volume methods, achieving conservation while not requiring the solution of a large matrix. Examples in which shock and rarefaction waves are well captured are presented and the propagation of acoustic pulses is well demonstrated. Copyright © 2003 John Wiley & Sons, Ltd.     

Characteristic Galerkin method for convection-diffusion equations and implicit algorithm using precise integration

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability. The project sponsored by the State Scientific and Technological Commission of China through “China State Key Project: the Theory and Methodology for Scientific and Engineering Computations with Large Scale”, the National Natural Science Foundation of China and the European Commission Research Project CI1^*CT94-0014.     

附加自由阻尼结构系统的有限元分析和拓扑优化设计方法研究

给出了一种适用于车身等复杂附加阻尼结构系统的有限元建模和动态特性分析方法,包括附加自由阻尼薄壁结构和附加阻尼材料粘弹性特性的有限元建模方法以及附加阻尼结构动态特性的有限元分析方法。在此基础上,又以强迫振动响应最小为优化设计目标,给出了一个附加阻尼结构的拓扑优化设计方法,包括优化设计问题的列式和适用的求解算法。通过对一个简单薄壁构件和一个车身地板上的附加阻尼结构的拓扑优化设计,验证了提出的有限元分析和拓扑优化设计方法的有效性。     

H 1 space-time discontinuous finite element method for convection-diffusion equations

An H~1 space-time discontinuous Galerkin (STDG) scheme for convectiondiffusion equations in one spatial dimension is constructed and analyzed. This method is formulated by combining the H~1 Galerkin method and the space-time discontinuous finite element method that is discontinuous in time and continuous in space. The existence and the uniqueness of the approximate solution are proved. The convergence of the scheme is analyzed by using the techniques in the finite difference and finite element methods. An optimal a-priori error estimate in the L~∞ (H~1 ) norm is derived. The numerical exper- iments are presented to verify the theoretical results.     

Low order nonconforming mixed finite element method for nonstationary incompressible Navier-Stokes equations

This paper studies a low order mixed finite element method (FEM) for nonstationary incompressible Navier-Stokes equations. The velocity and pressure are approximated by the nonconforming constrained Q ₁ ^rot element and the piecewise constant, respectively. The superconvergent error estimates of the velocity in the broken H ¹-norm and the pressure in the L ²-norm are obtained respectively when the exact solutions are reasonably smooth. A numerical experiment is carried out to confirm the theoretical results.     

形成不同拉压允许应力Michell桁架的有限元方法

提出了形成具有不同拉压允许应力材料的Michell桁架的有限元方法。采用正交异性材料模型,以结点位置杆件的密度和方向作为设计变量,由结点位置的值插值得到有限单元内部的材料性质,材料在设计域内连续分布,这样避免了单元之间材料性质的突变,因此从根本上避免了单元铰接、棋盘格现象以及单元依赖性等数值计算不稳定性问题。提出了一种基于优化准则法的迭代算法,用以形成不同拉压允许应力材料的Michell桁架。由有限元分析得到结点位置处的主应力方向作为杆件方向,根据主应力方向的应变和材料拉压允许应力调整杆件密度。最后由结点位置的杆件方向形成Michell桁架中连续分布杆件。