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.     

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

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.     

UNCONDITIONAL STABLE DIFFERENCE METHODS WITH INTRINSIC PARALLELISM FOR SEMILINEAR PARABOLIC SYSTEMS OF DIVERGENCE TYPE

The general finite difference schemes with intrinsic parallelism for the boundary value problem of the semilinear parabolic system of divergence type with bounded measurable coefficients is studied. By the approach of the discrete functional analysis, the existence and uniqueness of the discrete vector solutions of the nonlinear difference system with intrinsic parallelism are proved. Moreover the unconditional stability of the general difference schemes with intrinsic parallelism justified in the sense of the continuous dependence of the discrete vector solution of the difference schemes on the discrete initial data of the original problems in the discrete W_2~(2,1) (Q△) norms. Finally the convergence of the discrete vector solutions of the certain difference schemes with intrinsic parallelism to the unique generalized solution of the original semilinear parabolic problem is proved.     

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.     

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.     

Convergence of the Cesàro mean for rational orthonormal bases

This paper deals with the convergence of the Cesaro mean for the rational orthonormal bases. Provided the set of zeroes of rational orthonormal bases is formed by a periodic repetition of the same finite sequence, the explicit expression of so-called block-Fejer kernel is available, and some properties of the block-Fejer kernel are discussed. Based on the convergence of the block-Cesaro mean, the convergence of Cesaro mean is also provided.     

基于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.     

WEIGHTED NONLINEAR REGRESSION

Minimization of the weighted nonlinear sum of squares of differences may be converted to the minimization of sum of squares. The Gauss-Newton method is recalled and the length of the step of the steepest descent method is determined by substituting the steepest descent direction in the Gauss-Newton formula. The existence of minimum is shown.     

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     

Noncommutative geometry and motives: The thermodynamics of endomotives

We combine aspects of the theory of motives in algebraic geometry with noncommutative geometry and the classification of factors to obtain a cohomological interpretation of the spectral realization of zeros of L-functions. The analogue in characteristic zero of the action of the Frobenius on ?-adic cohomology is the action of the scaling group on the cyclic homology of the cokernel (in a suitable category of motives) of a restriction map of noncommutative spaces. The latter is obtained through the thermodynamics of the quantum statistical system associated to an endomotive (a noncommutative generalization of Artin motives). Semigroups of endomorphisms of algebraic varieties give rise canonically to such endomotives, with an action of the absolute Galois group. The semigroup of endomorphisms of the multiplicative group yields the Bost-Connes system, from which one obtains, through the above procedure, the desired cohomological interpretation of the zeros of the Riemann zeta function. In the last section we also give a Lefschetz formula for the archimedean local L-factors of arithmetic varieties.     

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.     

Analysis of the Stress State of the Surface Layer of a Polymeric Mass upon Filling of Bulky Compression Molds

The stress state of the surface layer of a polymeric mass during filling of bulky compression molds is analyzed. It is shown that, at particular rheological characteristics of the mass, temperature, and filling rates, cracking of the surface layer occurs, which leads to defects in the finished products. A physical analysis of this process makes it possible to conclude that the cracks arise due to the normal stresses operating in the front region of the moving polymeric mass. It is found that, under certain flow conditions, areas with a pressure lower than the atmospheric one appear on the surface of the polymer. If the tensile stresses arising in these local regions are higher than the tensile strength of the mass, the continuity of the composition is broken in the direction determined by the greatest rate of the normal deformation. To confirm the reliability of the crack-formation mechanism proposed, the distribution of the pressure and normal stresses over the free surface is calculated based on a numerical method. These calculations show that, by comparing the stress level achieved in the front region with the tensile-strength characteristics of the polymeric composition, it is possible to predict, with a sufficient accuracy, the possibility of crack formation in the surface layer of such a mass under given flow conditions and thus to solve the question on flawless manufacturing of products.     

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.     

基于综合重要度的动车组多部件系统机会维修策略研究

针对我国动车组列车现行维修方式,提出基于综合重要度序列的动车组多部件系统机会维修策略,对提高系统可靠度贡献大的关键部件进行准时优先维修。建立部件综合重要度指数计算模型,并依据其对部件维修优先级进行排序。以维修总成本最低为目标计算单部件最优维修周期及时刻,以系统维修总成本最低为目标,以关键部件的维修时刻为系统停机时刻建立考虑重要度的多部件系统机会维修模型。算例选取某型动车组四级修时更换的四部件系统为研究对象,讨论机会维修里程窗的大小及其偏移量对维修效果的影响,对比结果表明,考虑综合重要度的机会维修策略能够在维修费用基本持平的条件下,保证对系统可靠性贡献大的关键部件的可靠性,进而保证系统的整体可靠性。     

Modelling and simulation of micro-well formation

Physico-chemical processes on the micro-scale require new modelling concepts because some effects become dominating that are negligible for macroscopic systems. This is illustrated by a new method for the production of micro-wells based on the placement of a small drop of toluene on a plate of polystyrene. After droplet evaporation, a micro-well is left. A mathematical model has been developed to understand the elementary processes of the micro-well formation. The model accounts for: (1) growth of the drop on the substrate, (2) evaporation process of the solvent, (3) dissolution of the substrate, (4) flow rate in the evaporating drop caused by the pinning effect, including the vertical velocity profile, and (5) increase in the concentration of dissolved material followed by precipitation. In the modelling and simulation process, it could be shown that the method of drop production also has a significant influence on the shape of the micro-wells.     

The Dirichlet divisor problem,traces and determinants for complex powers of the twisted bi-Laplacian

Estimating the counting function for the eigenvalues of the twisted bi-Laplacian leads to the Dirichlet divisor problem, which is then used to compute the trace of the heat semigroup and the Dixmier trace of the inverse of the twisted bi-Laplacian. The zeta function regularizations of the traces and determinants of complex powers of the twisted bi-Laplacian are computed. A formula for the zeta function regularizations of determinants of heat semigroups of complex powers of the twisted bi-Laplacian is given.     

Multi‐dimensional Lotka–Volterra systems for carcinogenesis mutations

In the paper we consider three classes of models describing carcinogenesis mutations. Every considered model is described by the system of (n+1) equations, and in each class three models are studied: the first is expressed as a system of ordinary differential equations (ODEs), the second—as a system of reaction–diffusion equations (RDEs) with the same kinetics as the first one and with the Neumann boundary conditions, while the third is also described by the system of RDEs but with the Dirichlet boundary conditions. The models are formulated on the basis of the Lotka–Volterra systems (food chains and competition systems) and in the case of RDEs the linear diffusion is considered. The differences between studied classes of models are expressed by the kinetic functions, namely by the form of kinetic function for the last variable, which reflects the dynamics of malignant cells (that is the last stage of mutations). In the first class the models are described by the typical food chain with favourable unbounded environment for the last stage, in the second one—the last equation expresses competition between the pre‐malignant and malignant cells and the environment is also unbounded, while for the third one—it is expressed by predation term but the environment is unfavourable. The properties of the systems in each class are studied and compared. It occurs that the behaviour of solutions to the systems of ODEs and RDEs with the Neumann boundary conditions is similar in each class; i.e. it does not depend on diffusion coefficients, but strongly depends on the class of models. On the other hand, in the case of the Dirichlet boundary conditions this behaviour is related to the magnitude of diffusion coefficients. For sufficiently large diffusion coefficients it is similar independently of the class of models, i.e. the trivial solution that is unstable for zero diffusion gains stability. Copyright © 2009 John Wiley & Sons, Ltd.     

ML Estimation of the MultivariatetDistribution and the EM Algorithm

Maximum likelihood estimation of the multivariatetdistribution, especially with unknown degrees of freedom, has been an interesting topic in the development of the EM algorithm. After a brief review of the EM algorithm and its application to finding the maximum likelihood estimates of the parameters of thetdistribution, this paper provides new versions of the ECME algorithm for maximum likelihood estimation of the multivariatetdistribution from data with possibly missing values. The results show that the new versions of the ECME algorithm converge faster than the previous procedures. Most important, the idea of this new implementation is quite general and useful for the development of the EM algorithm. Comparisons of different methods based on two datasets are presented.     

Nonlinear Electromagnetic Delay of Electromagnetic Signals Propagating in the Magnetic Meridian Plane of Pulsars and Magnetars

We analyze the propagation of electromagnetic waves in the magnetic meridian planes of neutron stars with a strong magnetic field in the framework of the parameterized post-Maxwellian electrodynamics of the vacuum. The origin of these electromagnetic waves is the curvature emission of X-rays and gamma rays from high-energy electrons in the vicinity of the magnetic poles of neutron stars. We show that in the case of a slowly varying intensity of X-ray and gamma-ray emission, the delay of the slow normal mode of electromagnetic waves relative to the fast mode results in a shift of the time dependence of the intensity of the detected radiation with one polarization relative to that of the radiation with the orthogonal polarization. In the case of single X-ray or gamma-ray pulses, the delay effect results in the polarization of the detected pulse varying during the pulse length, the leading edge of all pulses being polarized normally to the magnetic equator plane of the neutron star. We note that the modern level of the experimental technique, in principle, allows observing the manifestations of the delay effect for signals of different polarizations.