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.     

On the Growth of Components of Meromorphic Solutions of Systems of Complex Differential Equations

This paper investigates the problem of the growth of the components of meromorphic solutions of a class of a system of complex algebraic differential equations, and generalized some of N. Toda＇s results concerning the growth of differential equations to the case of systems of differential equations. The paper considers the existence of admissible solutions of the system of differential equations.     

LOCAL AND GLOBAL HOPF BIFURCATIONS IN A DELAYED HUMAN RESPIRATORY SYSTEM

This paper considers a delayed human respiratory model. Firstly, the stability of the equilibrium of the model is investigated and the occurrence of a sequence of Hopf bifurcations of the model is proved. Secondly, the explicit algorithms which determine the direction of the Hopf bifurcations and the stability of the bifurcating periodic solutions are derived by applying the normal form method and the center manifold theory. Finally, the existence of the global periodic solutions is showed under some assumptions on the model.     

基于Legendre多项式的数值微分: 加权$L^2$空间中的收敛性分析

We consider the problem of estimating the derivative of a function f from its noisy version fδby using the derivatives of the partial sums of Fourier-Legendre series of f~δ. Instead of the observation L~2 space, we perform the reconstruction of the derivative in a weighted L~2 space. This takes full advantage of the properties of Legendre polynomials and results in a slight improvement on the convergence order. Finally, we provide several numerical examples to demonstrate the efficiency of the proposed method.     

拟线性可化约的双曲型方程组的奇性分析

In this paper we investigate the formation of singularities of hyperbolic systems. Employing the method of parametric coordinates and the existence of the solution of the blow-up system, we prove that the blow-up of classic solutions is due to the envelope of characteristics of the same family, analyze the geometric properties of the envelope of characteristics and estimate the blowup rates of the solution precisely.     

Time to most recent common ancestor for stationary continuous state branching processes with immigration

Motivated by sample path decomposition of the stationary continuous state branching process with immigration, a general population model is considered using the idea of immortal individual. We compute the joint distribution of the random variables： the time to the most recent common ancestor （MRCA）, the size of the current population, and the size of the population just before MRCA. We obtain the bottleneck effect as well. The distribution of the number of the oldest families is also established. These generalize the results obtained by Y. T. Chen and J. F. Delmas.     

Stability and Hopf bifurcation analysis of a predator-prey model with time delayed incomplete trophic transfer

A new nonlinear predator-prey model with incomplete trophic transfer is introduced. In this model,we assume that the rate of the trophic absorption of the predator is less than the rate of the conversion of consumed prey to predator in the Ivlev-type functional responses. The existence and uniqueness of the positive equilibrium of the model and the stability of the equilibrium of the model are studied under various conditions.Hopf bifurcation analysis of the delayed model is provided.     

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.     

广义中心-细鞍点的若干公式

In this paper, a problem of center-weak focus of a general system of degree n is transformed into a problem of generalized center-weak saddle. It offers: (1) the calculating formula of the saddle values of the former 4 orders in the general system of degree three; (2) some calculating formulae of the saddle values in the general homogeneous system of degree n. The problem of center-focus(saddle) plays an important role in the qualitative theory of ordinary differential equations. However, since the calculations are difficult, the results have not been satisfactory. On the basis of the author's work, Northeast. Math. J., 13(4)(1997), 487-494, the paper offers some calculating formulae of the saddle values.     

A local probability exponential inequality for the large deviation of an empirical process indexed by an unbounded class of functions and its application

A local probability exponential inequality for the tail of large deviation of an empirical process over an unbounded class of functions is proposed and studied. A new method of truncating the original probability space and a new symmetrization method are given. Using these methods, the local probability exponential inequalities for the tails of large deviations of empirical processes with non-i.i.d. independent samples over unbounded class of functions are established. Some applications of the inequalities are discussed. As an additional result of this paper, under the conditions of Kolmogorov theorem, the strong convergence results of Kolmogorov on sums of non-i.i.d. independent random variables are extended to the cases of empirical processes indexed by unbounded classes of functions, the local probability exponential inequalities and the laws of the logarithm for the empirical processes are obtained.     

基于模糊FNN-ELECTRE法的国有商业银行竞争力评价

随着金融业全方位开放,国有商业银行已成为国家经济命脉的核心,其竞争力强弱关系着国家经济的繁荣与衰退。本文基于模糊FNN-ELECTRE方法建立国有商业银行竞争力评价模型,以因素神经网络理论(FNN)与ELECTRE融合方法为基础,以现有商业银行竞争力评价指标为着眼,把国有商业银行竞争力的评价体系从现实竞争力和潜在竞争力两个因素抽取为一级指标,以规模、质量、业务结构、效率、成长性五个因素设为二级指标。从实证的角度对银行A、银行B、银行C、银行D、银行E五大国有银行进行分析验证,对其属性值进行和谐性与非和谐性检验得出,银行B竞争能力最强,银行E竞争能力最弱,五大国有银行均有改进空间。     

Wave parameter identification problem for ocean test structure data,part 2, discrete formulation

This paper deals with the solution of the wave parameter identification problem for ocean test structure data. A discrete formulation is assumed. An ocean test structure is considered, and wave elevation and velocities are assumed to be measured with a number of sensors. Within the frame of linear wave theory, a Fourier series model is chosen for the wave elevation and velocities. Then, the following problem is posed: Find the amplitudes of the various wave components of specified frequency and direction, so that the assumed model of wave elevation and velocities provides the best fit to the measured data. Here, the term best fit is employed in the least-square sense over a given time interval.At each time instant, the wave representation involves four indexes (frequency, direction, instrument, time); hence, four-dimensional arrays are required. This formal difficulty can be avoided by switching to an alternative representation involving only two indexes (frequency-direction, instrument-time); hence, standard vector-matrix notation can be used. Within this frame, optimality conditions are derived for the amplitudes of the assumed wave model.A characteristic of the wave parameter identification problem is that the condition number of the system matrix can be large. Therefore, the numerical solution is not an easy task and special procedures must be employed. Specifically, Gaussian elimination is avoided and advantageous use is made of the Householder transformation, in the light of the least-square nature of the problem and the discretized approach to the problem.Numerical results are presented. The effect of various system parameters (number of frequencies, number of directions, sampling time, number of sensors, and location of sensors) is investigated in connection with global or strong accuracy, local or weak accuracy, integral accuracy, and condition number of the system matrix.From the numerical experiments, it appears that the wave parameter identification problem has a unique solution if the number of directions is smaller than or equal to the number of sensors; it has an infinite number of solutions otherwise. In the case where a unique solution exists, the condition number of the system matrix increases as the size of the system increases, and this has a detrimental effect on the accuracy. However, the accuracy can be improved by proper selection of the sampling time and by proper choice of the number and location of the sensors.Generally speaking, the computations done for the discrete case exhibit better accuracy than the computations done for the continuous case (Ref. 5). This improved accuracy is a direct consequence of having used advantageously the Householder transformation and is obtained at the expense of increased memory requirements and increased CPU time.This work was supported by Exxon Production Research Company, Houston, Texas. This paper is based partly on Refs. 1–4.     

Formulation of linear problems and solution by a universal machine

Using the predicate language for ordered fields a class of problems referred to aslinear problems is defined. This class contains, for example, all systems of linear equations and inequalities, all linear programming problems, all integer programming problems with bounded variables, all linear complementarity problems, the testing of whether sets that are defined by linear inequalities are semilattices, all satisfiability problems in sentenial logic, the rank-computation of matrices, the computation of row-reduced echelon forms of matrices, and all quadratic programming problems with bounded variables. A single, one, algorithm, to which we refer as theUniversal Linear Machine, is described. It solves any instance of any linear problem. The Universal Linear Machine runs in two phases. Given a linear problem, in the first phase a Compiler running on a Turing Machine generates alinear algorithm for the problem. Then, given an instance of the linear problem, in the second phase the linear algorithm solves the particular instance of the linear problem. The linear algorithm is finite, deterministic, loopless and executes only the five ordered field operations — additions, multiplications, subtractions, divisions and comparisons. Conversely, we show that for each linear algorithm there is a linear problem which the linear algorithm solves uniquely. Finally, it is shown that with a linear algorithm for a linear problem, one can solve certain parametric instances of the linear problem.Research was supported in part by the National Science Foundation Grant DMS 92-07409, by the Department of Energy Grant DE-FG03-87-ER-25028, by the United States—Israel Binational Science Foundation Grant 90-00434 and by ONR Grant N00014-92-J1142.Corresponding author.     

A double-stage piecewise recursive filter for real-time applications

This paper examines the basic problem of estimating the state of a system described by a set of dynamical equations. The state estimation is performed by means of optimal filtering techniques in which the state observation is given by a set of nonlinear equations.The basic issue encountered in most engineering problems is the availability of more observable information than the processor can effectively process in real-time, while performing full optimal filtering on all available observable data. To alleviate this deficiency, an algorithm is developed in which the rate of information processing is kept at a high value, while the system statistics are evaluated at a much slower rate.The method utilizes a combination of fast and slow filtering loops, in which the observable data are processed at a high rate during the fast loop, while the system error covariance, gain computations, and all other system statistics are processed at a lower rate in the slow loop. Methods are provided to resolve any incompatibility in the system statistics resulting from this fast-loop/slow-loop processing combination.Applications describing problems pertaining to aircraft navigation are presented. Specifically, applications to aircraft navigation through a satellite network are studied. Appropriate simulation results from the above studies are shown.This work was performed under contract with the Department of the Air Force, Space and Missile Organization (SAMSO), Los Angeles, California, Contract No. F04-701-75-C-0180.The authors are indebted to Dr. C. Johnson, GPS System Engineer, Mr. G. Consolver, GPS Software Manager, and Mr. W. Riley, Communications/Navigation Department, Texas Instruments for their continuous encouragement in the course of this work. Also, they are indebted to Major M. Birnbaum, SAMSO, for his constructive criticism on the conceptual design of this work, and to Mr. A. Bierman, Aerospace Corporation, for kindly providing simulation data of aircraft trajectories. Finally, the authors would like to express their gratitude to Dr. N. Carlson, Intermetrics, for valuable consultation during the progress of this work.Portions of this paper were presented by the senior author at the GSP Conference, Plan 76, IEEE Position, Location, and Navigation Symposium, San Diego, California, 1976.     

The generalized hierarchical product of graphs

A generalization of both the hierarchical product and the Cartesian product of graphs is introduced and some of its properties are studied. We call it the generalized hierarchical product. In fact, the obtained graphs turn out to be subgraphs of the Cartesian product of the corresponding factors. Thus, some well-known properties of this product, such as a good connectivity, reduced mean distance, radius and diameter, simple routing algorithms and some optimal communication protocols, are inherited by the generalized hierarchical product. Besides some of these properties, in this paper we study the spectrum, the existence of Hamiltonian cycles, the chromatic number and index, and the connectivity of the generalized hierarchical product.     

考虑成本、等待时间和安全水平的分类安检模式研究

近年来,随着民航机场旅客吞吐量快速增加,机场安检部门的工作压力正不断加大,旅客等待时间也正在逐渐增加。针对该问题,引入了分类安检模式。然而,分类安检模式通常需要额外的投资成本,可能会产生较大的财务负担。因此,以成本为目标,综合考虑旅客等待时间和安全水平的要求,研究了分类安检模式和传统安检模式的比较以及最优决策性质等相关问题。首先,考虑安全水平和旅客等待时间的约束,分别为传统安检模式和分类安检模式建立了模型,并分析了其最优运营决策性质。其次,从成本的角度,对两种安检模式做了比较研究。研究结果表明,当旅客数量多、安检犯错容忍度较高或旅客等待敏感时,分类安检模式比传统安检模式表现更优;反之,传统安检模式更优。此外,在分类安检模式中,分配到高风险安检通道的旅客比例存在一个最优值,并且其通常在20%到30%之间取到。最后,通过数值分析验证了分类安检模式的最优服务配置。     

Representation of the solution to a model problem in semiconductor physics

We study a mixed type problem for the Poisson equation arising in the modeling of charge transport in semiconductor devices [V. Romano, 2D simulation of a silicon MESFET with a non-parabolic hydrodynamical model based on the maximum entropy principle, J. Comput. Phys. 176 (2002) 70-92; A.M. Blokhin, R.S. Bushmanov, A.S. Rudometova, V. Romano, Linear asymptotic stability of the equilibrium state for the 2D MEP hydrodynamical model of charge transport in semiconductors, Nonlinear Anal. 65 (2006) 1018-1038]. Unlike well-studied elliptic boundary-value problems in domains with smooth boundaries (see, for example, [O.A. Ladyzhenskaya, N.N. Uralceva, Linear and Quasilinear Elliptic Equations, Nauka, Moscow, 1973; D. Gilbarg, N.S. Trudinger, Elliptic Partial Differential Equations of Second Order, Springer-Verlag, Berlin, 1983]), our problem has two significant features: firstly, the boundary is not a smooth curve and, secondly, the type of boundary conditions is mixed (the Dirichlet condition is satisfied on the one part of the boundary whereas the Neumann condition on the other part). The well-posedness of the problem in Hölder and Sobolev spaces is proved. The representation of the solution to the problem is obtained in an explicit form.     

公众参与视角下的央地环境规制博弈分析

本文首先构建环境规制中中央政府和地方政府的两方演化博弈模型,并在此基础上将公众作为第三方参与主体,构建中央政府、地方政府和公众三方演化博弈模型,详细比较两方和三方博弈模型的区别,探究各个主体策略行为的影响因素。研究发现:(1)未有公众参与下,地方政府策略选择主要受地方政府积极执行成本、环境收益、经济损失,消极执行的环境政绩损失,中央政府监管力度、治理补贴和对地方政府的处罚等因素影响;中央政府的监管策略主要受到严格监管的成本以及对地方政府的治理补贴和处罚等因素影响。引入公众参与后,在央地两方博弈的基础上,地方政府环境规制执行策略的影响因素还增加了地方政府被举报后所受到的追加处罚,中央政府监管策略的影响因素还增加了中央政府监管力度、对地方政府的追加处罚以及中央政府的公信力损失。(2)未有公众参与下,中央政府严格监管的概率随地方政府积极执行概率的增大而减小。引入公众参与后,中央政府严格监管率随地方政府积极执行概率的增大而增大。说明在公众参与下,地方政府积极执行环境规制对中央政府严格监管产生的抑制作用转变成了促进作用。(3)地方政府积极执行的概率、中央政府严格监管的概率都随公众举报概率的增大而增大。说明公众参与不仅对地方政府承担环保责任具有促进作用,而且有利于促使中央政府落实环境治理政策。     

Rolling-ball method for estimating the boundary of the support of a point-process intensity

We suggest a generalisation of the convex-hull method, or ‘DEA’ approach, for estimating the boundary or frontier of the support of a point cloud. Figuratively, our method involves rolling a ball around the cloud, and using the equilibrium positions of the ball to define an estimator of the envelope of the point cloud. Constructively, we use these ideas to remove lines from a triangulation of the points, and thereby compute a generalised form of a convex hull. The radius of the ball acts as a smoothing parameter, with the convex-hull estimator being obtained by taking the radius to be infinite. Unlike the convex-hull approach, however, our method applies to quite general frontiers, which may be neither convex nor concave. It brings to these contexts the attractive features of the convex hull: simplicity of concept, rotation-invariance, and ready extension to higher dimensions. It admits bias corrections, which we describe and illustrate through implementation.