Uzawa iteration method for stokes type variational inequality of the second kind

In this paper,the Uzawa iteration algorithm is applied to the Stokes problem with nonlinear slip boundary conditions whose variational formulation is the variational inequality of the second kind.Firstly, the multiplier in a convex set is introduced such that the variational inequality is equivalent to the variational identity.Moreover,the solution of the variational identity satisfies the saddle-point problem of the Lagrangian functional ?.Subsequently,the Uzawa algorithm is proposed to solve the solution of the saddle-point problem. We show the convergence of the algorithm and obtain the convergence rate.Finally,we give the numerical results to verify the feasibility of the Uzawa algorithm.     

THE RELATIONSHIP BETWEEN PROJECTIVE GEOMETRIC AND RATIONAL QUADRATIC B-SPLINE CURVES

For two rational quadratic B-spline curves with same control vertexes, the cross ratio of four eollinear points are represented; which are any one of the vertexes, and the two points that the ray initialing from the vertex intersects with the corresponding segments of the twocurves, and the point the ray intersecting with the connecting line between the two neighboring vertexes. Different from rational quadratic Beeier curves, the value is generally related with the loeation of the ray, and the necessary and sufficient condition o5 the ratio being independent of the ray‘s loeation is showed. Alsn another cross ratio o5 the following four collinear points are suggested, i.e. one vertex, the points that the ray from the initlal vertex intersects respectivdy with the curve segmentt the line connecting the segments end points, and the line connecting the two neighboring vertexes. This cross ratio is concerned only whh the ray‘s location, butnot with the weights of the curve. Furthermore, the cross ratio is projective invariant under the projective transformation between the two segments.     

Boundary Shape Control of the Navier-Stokes Equations and Applications

In this paper, the geometrical design for the blade's surface in an impeller or for the profile of an aircraft, is modeled from the mathematical point of view by a boundary shape control problem for the Navier-Stokes equations. The objective function is the sum of a global dissipative function and the power of the fluid. The control variables are the geometry of the boundary and the state equations are the Navier-Stokes equations. The Euler-Lagrange equations of the optimal control problem are derived, which are an elliptic boundary value system of fourth order, coupled with the Navier-Stokes equations. The authors also prove the existence of the solution of the optimal control problem, the existence of the solution of the Navier-Stokes equations with mixed boundary conditions, the weak continuity of the solution of the Navier-Stokes equations with respect to the geometry shape of the blade's surface and the existence of solutions of the equations for the Gateaux derivative of the solution of the Navier-Stokes equations with respect to the geometry of the boundary.     

PROPERTIES OF SOLUTIONS OF n-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

Consider the n-dimensional incompressible Navier-Stokes equations ?/(?t)u-α△u +(u · ?)u + ?p = f(x, t), ? · u = 0, ? · f = 0,u(x, 0) = u0(x), ? · u0= 0.There exists a global weak solution under some assumptions on the initial function and the external force. It is well known that the global weak solutions become sufficiently small and smooth after a long time. Here are several very interesting questions about the global weak solutions of the Cauchy problems for the n-dimensional incompressible Navier-Stokes equations.· Can we establish better decay estimates with sharp rates not only for the global weak solutions but also for all order derivatives of the global weak solutions?· Can we accomplish the exact limits of all order derivatives of the global weak solutions in terms of the given information?· Can we use the global smooth solution of the linear heat equation, with the same initial function and the external force, to approximate the global weak solutions of the Navier-Stokes equations?· If we drop the nonlinear terms in the Navier-Stokes equations, will the exact limits reduce to the exact limits of the solutions of the linear heat equation?· Will the exact limits of the derivatives of the global weak solutions of the Navier-Stokes equations and the exact limits of the derivatives of the global smooth solution of the heat equation increase at the same rate as the order m of the derivative increases? In another word, will the ratio of the exact limits for the derivatives of the global weak solutions of the Navier-Stokes equations be the same as the ratio of the exact limits for the derivatives of the global smooth solutions for the linear heat equation?The positive solutions to these questions obtained in this paper will definitely help us to better understand the properties of the global weak solutions of the incompressible Navier-Stokes equations and hopefully to discover new special structures of the Navier-Stokes equations.     

An algorithm for decomposing a polynomial system into normal ascending sets

We present an algorithm to decompose a polynomial system into a finite set of normal ascending sets such that the set of the zeros of the polynomial system is the union of the sets of the regular zeros of the normal ascending sets.If the polynomial system is zero dimensional,the set of the zeros of the polynomials is the union of the sets of the zeros of the normal ascending sets.     

求解一种带有微分方程约束的条件变分问题的多目标最优化方法

Aiming at the conditional variational problem with the bi-objective functionals and the differential equational constraint in the optimal design of the electrostatic lenses, the conditional variational problem is transformed into multi-objective optimization problem by means of the spline function and the integral transformation. For solving the transformed problem, the analytic representation formula of the optimal solution about the original problem is obtained with regard to the voltage distribution and the electron trajectory. It will provide a new effective method for the design of the electrostatic lenses.     

The unconditional stable and convergent difference methods with intrinsic parallelism for quasilinear parabolic systems

A kind of the general finite difference schemes with intrinsic parallelism forthe boundary value problem of the quasilinear parabolic system is studied without assum-ing heuristically that the original boundary value problem has the unique smooth vectorsolution. By the method of a priori estimation of the discrete solutions of the nonlineardifference systems, and the interpolation formulas of the various norms of the discretefunctions and the fixed-point technique in finite dimensional Euclidean space, the exis-tence and uniqueness of the discrete vector solutions of the nonlinear difference systemwith intrinsic parallelism are proved. Moreover the unconditional stability of the generalfinite difference schemes with intrinsic parallelism is justified in the sense of the continu-ous dependence of the discrete vector solution of the difference schemes on the discretedata of the original problems in the discrete w_2~(2,1) norms. Finally the convergence of thediscrete vector solutions of the certain differe     

GLOBAL SOLUTION TO THE CAUCHY PROBLEM ON A UNIVERSE FIREWORKS MODEL

We prove existence and uniqueness of the global solution to the Cauchy problem on a universe fireworks model with finite total mass at the initial state when the ratio of the mass surviving the explosion, the probability of the explosion of fragments and the probability function of the velocity change of a surviving particle satisfy the corresponding physical conditions. Although the nonrelativistic Boltzmann-like equation modeling the universe fireworks is mathematically easy, this article leads rather theoretically to an understanding of how to construct contractive mappings in a Banach space for the proof of the existence and uniqueness of the solution by means of methods taken from the famous work by DiPerna ＆ Lions about the Boltzmann equation. We also show both the regularity and the time-asymptotic behavior of solution to the Cauchy problem.     

Mesoscopic Numerical Study on Flow Boiling Heat Transfer Performance in Channels With Multiple Rectangular Heaters北大核心CSCD

The flow boiling phenomenon in a channel with multiple rectangular heaters under a constant wall temperature was numerically studied with the lattice Boltzmann method. The effects of spacings between heaters, heater lengths and heater surface wettabilities on the bubble morphology, the bubble area and the heat flux on the heater surface, were studied. The results show that, the bubble growth rate increases with the spacing between heaters. The larger the bubble area is, the earlier the nucleated bubbles will leave the heater surface. The corresponding boiling heat transfer performance increases by 12% with the spacing between heaters growing from 250 lattices to 1 000 lattices. On the other hand, the longer the heater length is, the earlier the bubble will nucleate and leave the heater surface, and the better the boiling heat transfer performance will be. The boiling heat transfer performance increases by 13% with the heater length rising from 16 lattices to 22 lattices. In addition, the bubble nucleates later on the hydrophilic surface than on the hydrophobic surface. Compared with the hydrophilic surface, the hydrophobic surface retains residual bubbles after the leaving of bubbles from the heater. The average heat flux and the bubble area of the hydrophilic surface are less than those of the hydrophobic surface. With the contact angle changing from 77° to 120°, the heat transfer performance increases by 26%. Finally, the orthogonal test results indicate that, the wettability of the heat exchanger surface has the greatest influence on the flow boiling heat transfer performance, while the heater length has the least influence. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Global regularity for 3D magneto-hydrodynamics equations with only horizontal dissipation

We investigate the Cauchy problem for the 3D magneto-hydrodynamics equations with only horizontal dissipation for the small initial data. With the help of the dissipation in the horizontal direction and the structure of the system, we analyze the properties of the decay of the solution and apply these decay properties to get the global regularity of the solution. In the process, we mainly use the frequency decomposition in Green's function method and energy method.     

颤振分析中的特征值排序问题

颤振分析中判断颤振临界速度的重要依据是系统V-g和V-f图,即系统特征值随参数的变化曲线.在几乎所有商用软件及自编程序的输出结果中,有时会出现所谓的"窜支"现象,这给颤振临界速度和颤振穿越分支及耦合形式的判断带来很大不便.通过隐函数定理可以证明,除重特征值点以外,系统特征值连续依赖于系统参数变化.依据多元向量值函数连续性,建立对特征值的排列算法,给出系统特征根轨迹的正确曲线,再输出V-g和V-f图数据,从而避免"窜支"现象.编制应用程序,通过几个典型算例对算法进行了验证.该工作能够有效简化颤振分析的后处理工作,提高分析效率.     

Creep of polymer concrete in the nonlinear region

Two polyester-based polymer concretes with various volume content of diabase as an extender and aggregate are tested in creep under compression at different stress levels. The phenomenological and structural approaches are both used to analyze the experimental data. Common features of changes in the instantaneous and creep compliances are clarified, and a phenomenological creep model which accounts for the changes in the instantaneous compliance and in the retardation spectrum depending on the stress level is developed. It is shown that the model can be used to describe the experimental results of stress relaxation and creep under repeated loading. Modeling of the composite structure and subsequent solution of the optimization problem confirm the possibility of the existence of an interphase layer more compliant than the binder. A direct correlation between the interphase volume content and the instantaneous compliance of the composite is revealed. It is found that the distinction in nonlinearity of the viscoelastic behavior of the two polymer concretes under investigation can be due to the difference in their porosity. Submitted to the 11th International Conference on Mechanics of Composite Materials (Riga, June 11–15, 2000.) Translated from Mekhanika Kompozitnykh Materialov, Vol. 36, No. 2, pp. 147–164, 2000.     

On the hyperbolicity of the period-doubling fixed point

We give a new proof of the hyperbolicity of the fixed point for the period-doubling renormalization operator using the local dynamics near a semi-attractive fixed point (in a Banach space) and the theory of holomorphic motions. We also give a new proof of the exponential contraction of the Feigenbaum renormalization operator in the hybrid class of the period-doubling fixed point: our proof uses the non-existence of invariant line fields in the period-doubling tower (C. McMullen), the topological convergence (D. Sullivan), and a new infinitesimal argument.

     

考虑服务和退货的双渠道供应链定价问题研究

在考虑零售服务和消费者退货情况下,运用两阶段优化方法分别探究了集中式和分散式两种双渠道供应链的最优定价策略。研究发现零售服务有助于提高零售商的渠道权力,并始终使制造商获益,但只有当服务水平小于某一临界值时,才使制造商和零售商同时得到帕累托改进。零售价格和直销价格与退货率正相关,而与服务水平的关系受渠道结构的影响。服务水平的提高加剧了供应链的“双重边际化”效应,制造商通过生产与消费者需求相匹配的产品来降低退货风险可以减弱这种效应,从而提高供应链效率。     

树脂基复合材料粘接的优化设计

本文在试验数据和理论分析的基础上,提出树脂基复合材料单搭接的优化设计的原则,胶粘剂的选择、被粘物的设计,粘接长度的选择、胶层厚度的选择等.认为通过优化设计可以提高粘接强度,减轻产品的重量,提高其质量.     

Stress intensity factors for a finite plate with an inclined crack by boundary collocation

In this paper, we combine the Muskhelishvili＇s complex variable method and boundary collocation method, and choose a set of new stress function based on the stress boundary condition of crack surface, the higher precision and less computation are reached. This method is applied to calculating the stress intensity factor for a finite plate with an inclined crack. The influence of θ （the obliquity of crack） on the stress intensity factors, as well as the number of summation terms on the stress intensity factor are studied and graphically represented.     

拍卖机制设计在区块链中的应用与挑战

区块链是新一代信息技术的重要组成部分,是分布式网络、加密技术、智能合约等多种技术集成的新型数据库软件。过去的十多年,区块链技术在全球范围内产生广泛影响。如今的区块链技术,已从最初的关注于解决货币和支付的去中心化问题,转入到解决市场的去中心化问题。智能合约的出现使得基于区块链技术的去中心化金融进入高速发展状态,也涌现出区块链环境下的各类拍卖场景。本文首次从机制设计角度,以区块链交易费机制,非同质化代币(Non-Fungible Token,NFT)拍卖和矿工可提取价值(Miner-Extractable Value,MEV)交易位置拍卖为主要对象,总结和剖析近些年来区块链上特有的拍卖机制;并针对区块链特性,提出区块链上拍卖机制设计所面临的挑战和未来亟待解决的问题。     

Adiabatic divergence of the chaotic layer width and acceleration of chaotic and noise-induced transport

We show that, in spatially periodic Hamiltonian systems driven by a time-periodic coordinate-independent (AC) force, the upper energy of the chaotic layer grows unlimitedly as the frequency of the force goes to zero. This remarkable effect is absent in any other physically significant systems. It gives rise to the divergence of the rate of the spatial chaotic transport. We also generalize this phenomenon for the presence of a weak noise and weak dissipation. We demonstrate for the latter case that the adiabatic AC force may greatly accelerate the spatial diffusion and the reset rate at a given threshold.     

Numerical study on the dynamics of a jet-vortex interaction

The interaction of a vortex and a jet plays an important role for many industrial processes such as carbon black formation or combustion in diesel engines. The knowledge of physics of these phenomena is crucial for engineers, but also for scientists who wish to study this complex flow pattern. In this research, we numerically investigate cases where a cylindrical reactor is charged with gas injected through five ports. The first one is located along the main axis of the cylinder and this leads to the formation of the main jet. The other four ports are situated along the side walls such that the gas entering the cylinder tangentially through them causes a vortex to be formed. The objective of this paper is to show the fundamental physical phenomena and also how the initial and boundary conditions influence the results. It is shown how the process is influenced not only by modifying the velocity of injected gas, but also by the orientation and position of the vortex inlet nozzles.     

Peridynamics via finite element analysis

Peridynamics is a recently developed theory of solid mechanics that replaces the partial differential equations of the classical continuum theory with integral equations. Since the integral equations remain valid in the presence of discontinuities such as cracks, the method has the potential to model fracture and damage with great generality and without the complications of mathematical singularities that plague conventional continuum approaches. Although a discretized form of the peridynamic integral equations has been implemented in a meshless code called EMU, the objective of the present paper is to describe how the peridynamic model can also be implemented in a conventional finite element analysis (FEA) code using truss elements. Since FEA is arguably the most widely used tool for structural analysis, this implementation may hasten the verification of peridynamics and significantly broaden the range of problems that the practicing analyst might attempt. Also, the present work demonstrates that different subregions of a model can be solved with either the classical partial differential equations or the peridynamic equations in the same calculation thus combining the efficiency of FEA with the generality of peridynamics. Several example problems show the equivalency of the FEA and the meshless peridynamic approach as well as demonstrate the utility and robustness of the method for problems involving fracture, damage and penetration.