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.     

Exact and approximate numerical solutions of steady-state distributions arising in the queueGI/G/1

In this paper we first obtain, in a unified way, closed-form analytic expressions in terms of roots of the so-called characteristic equation (c.e.), and then discuss the exact numerical solutions of steady-state distributions of (i) actual queueing time, (ii) virtual queueing time, (iii) actual idle time, and (iv) interdeparture time for the queueGI/R/1, whereR denotes the class of distributions whose Laplace-Stieltjes transforms (LSTs) are rational functions (ratios of a polynomial of degree at mostn to a polynomial of degreen). For the purpose of numerical discussions of idle- and interdeparture-time distributions, the interarrival-time distribution is also taken to belong to the classR. It is also shown that numerical computations of the idle-time distribution ofR/G/1 queues can be done even ifG is not taken asR. Throughout the discussions it is assumed that the queue discipline is first-come-first-served (FCFS). For the tail of the actual queueing-time distribution ofGI/R/1, approximations in terms of one or more roots of the c.e. are also discussed. If more than one root is used, they are taken in ascending order of magnitude. Numerical aspects have been tested for a variety of complex interarrival- and service-time distributions. The analysis is not restricted to generalized distributions with phases such as Coxian-n (C _n ), but also covers nonphase type distributions such as uniform (U) and deterministic (D). Some numerical results are also presented in the form of tables and figures. It is expected that the results obtained from the present study should prove to be useful not only to practitioners, but also to queueing theorists who would like to test the accuracies of inequalities, bounds or approximations.     

Expected hitting and cover times of random walks on some special graphs

We give general bounds (and in some cases exact values) for the expected hitting and cover times of the simple random walk on some special undirected connected graphs using symmetry and properties of electrical networks. In particular we give easy proofs for an N–1H_N-1 lower bound and an N² upper bound for the cover time of symmetric graphs and for the fact that the cover time of the unit cube is Φ(NlogN). We giver a counterexample to a conjecture of Freidland about a general bound for hitting times. Using the electric approach, we provide some genral upper and lower bounds for the expected cover times in terms of the diameter of the graph. These bounds are tight in many instances, particularly when the graph is a tree. © 1994 John Wiley & Sons, Inc.     

圈上的多重懒惰随机游走

本文考虑了n个定点的圈上的多重懒惰随机游走.利用偶和方法证明了其最大相遇时的期望的阶数为h_(max)×log n,其中h_(max)为圈上的一简单随机游走的最大击中时.     

Total Completion Time in a Two-machine Flowshop with Deteriorating Tasks

This paper deals with two-machine flowshop problems with deteriorating tasks, i.e. tasks whose processing times are a nondecreasing function that depend on the length of the waiting periods. We consider the so-called Restricted Problem. This problem can be defined as follows: for a given permutation of tasks, find an optimal placement on two machines so that the total completion time is minimized. We will show that the Restricted Problem is nontrivial. We give some properties for the optimal placement and we propose an optimal placement algorithm.      

A randomized first-passage problem for drifted Brownian motion subject to hold and jump from a boundary

We study an inverse first-passage-time problem for Wiener process X(t) subject to hold and jump from a boundary c. Let be given a threshold S > X(0) ≥ c, and a distribution function F on [0, +∞). The problem consists in finding the distribution of the holding time at c and the distribution of jumps from c, so that the first-passage time of X(t) through S has distribution F.     

An occupation time identity for reflecting Brownian motion with drift

Summary This note is about an occupation time identity derived in [14] for reflecting Brownian motion with drift ]]>]]>]]>]]>]]>]]>]]>]]>]]>]]>]]>]]>]]>-\mu<0,$ RBM($-\mu$), for short. The identity says that for RBM($-\mu$) in stationary state ]]>(I^{+}_t, I^{-}_t) \rr (t-G_t,D_t-t),\qquad t\in \mathbb{R},$$ where $G_t$ and $D_t$ denote the starting time and the ending time, respectively, of an excursion from 0 to 0 (straddling $t$) and $I^{+}_t$ and $I^{-}_t$ are the occupation times above and below, respectively, of the observed level at time $t$ during the excursion. Due to stationarity, the common distribution does not depend on $t.$ In fact, it is proved in [9] that the identity is true, under some assumptions, for all recurrent diffusions and stationary processes. In the null recurrent diffusion case the common distribution is not, of course, a probability distribution. The aim of this note is to increase understanding of the identity by studying the RBM($-\mu$) case via Ray--Knight theorems.     

Brown运动关于球的末遇时首出时及停留时的矩

设x=(X