Boundary interpolation in the ball

We solve a boundary interpolation problem in the reproducing kernel Hilbert space of functions analytic in the unit ball of with reproducing kernel 1/(1−∑₁^Nz_kw_k^*). We introduce the notion of Brune factor (or Blaschke–Potapov factor of the third kind) in this setting.     

A stable and accurate projection method on a locally refined staggered mesh

In this paper, a robust projection method on a locally refined mesh is proposed for two‐ and three‐dimensional viscous incompressible flows. The proposed method is robust not only when the interface between two meshes is located in a smooth flow region but also when the interface is located in a flow region with large gradients and/or strong unsteadiness. In numerical simulations, a locally refined mesh saves many grid points in regions of relatively small gradients compared with a uniform mesh. For efficiency and ease of implementation, we consider a two‐level blocked structure, for which both of the coarse and fine meshes are uniform Cartesian ones individually. Unfortunately, the introduction of the two‐level blocked mesh results in an important but difficult issue: coupling of the coarse and fine meshes. In this paper, by properly addressing the issue of the coupling, we propose a stable and accurate projection method on a locally refined staggered mesh for both two‐ and three‐dimensional viscous incompressible flows. The proposed projection method is based on two principles: the linear interpolation technique and the consistent discretization of both sides of the pressure Poisson equation. The proposed algorithm is straightforward owing to the linear interpolation technique, is stable and accurate, is easy to extend from two‐ to three‐dimensional flows, and is valid even when flows with large gradients cross the interface between the two meshes. The resulting pressure Poisson equation is non‐symmetric on a locally refined mesh. The numerical results for a series of exact solutions for 2D and 3D viscous incompressible flows verify the stability and accuracy of the proposed projection method. The method is also applied to some challenging problems, including turbulent flows around particles, flows induced by impulsively started/stopped particles, and flows induced by particles near solid walls, to test the stability and accuracy. Copyright © 2010 John Wiley & Sons, Ltd.     

Free surface Stokes flows obstructed by multiple obstacles

Gravity‐driven Stokes flow down an inclined plane over and around multiple obstacles is considered. The flow problem is formulated in terms of a boundary integral equation and solved using the boundary element method. A Hermitian radial basis function (RBF) is used for the interpolation of the free surface, generation of the unit normal and curvature, and to prescribe the far‐field conditions. For flow over an obstacle, hemispheres are taken. For flow around an obstacle, circular cylinders are modelled and the contact angle condition on the obstacle/free surface intersection specified using the RBF formulation. Explicit profiles are produced for flow over and around two obstacles placed in various locations relative to one another. Interaction due to two obstacles is given by comparisons made with the profiles for flow over and around individual obstacles. In general, when the obstacles are separated by a sufficiently large distance the flow profiles are identical to a single obstacle analysis. For flow over and around two obstacles in‐line with the incident flow, effects of the governing parameters are examined, with variations in plane inclination angle, Bond number, obstacle size, and in the case of obstacles intersecting the free surface, static contact angle is considered. Finally flows over and around three obstacles are modelled. Copyright © 2009 John Wiley & Sons, Ltd.     

An algorithm for constrained convex optimization

This paper presents a feasible direction algorithm for the minimization of a pseudoconvex function over a smooth, compact, convex set. We establish that each cluster point of the generated sequence is an optimal solution of the problem without introducing anti-jamming procedures. Each iteration of the algorithm involves as subproblems only one line search for a zero of a continuously differentiable convex function and one univariate function minimization on a compact interval.     

Hilbert function spaces and the Nevanlinna-Pick problem on the polydisc

In [4], Agler and McCarthy (2005) used Hilbert function spaces to study the uniqueness properties of the Nevanlinna-Pick problem on the bidisc. In this work we give a geometric procedure for constructing a Nevanlinna-Pick problem on Dn with a specified set of uniqueness. On the way to establishing this procedure, we prove a result about Hilbert function spaces and partially answer a question posed by Agler and McCarthy.     

Support-type properties of generalized convex functions

Chebyshev systems induce in a natural way a concept of convexity. The functions convex in this sense behave in many aspects similarly to ordinary convex functions. In this paper support-type properties are investigated. Using osculatory interpolation, the existence of support-like functions is established for functions convex with respect to Chebyshev systems. Unique supports are determined. A characterization of the generalized convexity via support properties is presented.     

Coherent trace-class operators in Hilbert couples

In this paper it is proved that the ideal of trace-class operators acting in a couple of Hilbert spaces coincides with the ideal of coherent trace-class operators. A new formula is derived for theK-functional in the couple of algebras of all bounded linear operators acting in a Hilbert couple, and new interpolation theorems are proved for trace class operators.Translated fromMatematicheskie Zametki, Vol. 63, No. 6, pp. 866–872, June, 1998.This research was partially supported by the International Science Foundation and the Russian Government under grant JD 7100.     

基于插值的全矢量有限差分法求解阶跃光波导本征模式

给出了一种适合求解阶跃折射率分布光波导本征模式的全矢量有限差分法.该方法充分考虑了电场分量在折射率分布跃变处的不连续性, 在拉格朗日插值的基础上得到了一种精度较高的差分格式.为了验证其有效性,对一种典型结构、折射率阶跃分布的波导作了计算,并与其它方法的结果作了比较.最后对不理想刻蚀情况做了模拟计算.     

Time shifted aliasing error upper bounds for truncated sampling cardinal series

Time shifted aliasing error upper bound extremals for the sampling reconstruction procedure are fully characterized. Sharp upper bounds are found on the aliasing error of truncated cardinal series and the corresponding extremals are described for entire functions from certain specific Lp, p>1, classes. Analogous results are obtained in multidimensional regular sampling. Truncation error analysis is provided in all cases considered. Moreover, sharpness of bounding inequalities, convergence rates and various sufficient conditions are discussed.     

Estimate of the truncation error of finite volume discretization of the Navier–Stokes equations on colocated grids

A methodology is proposed for the calculation of the truncation error of finite volume discretizations of the incompressible Navier–Stokes equations on colocated grids. The truncation error is estimated by restricting the solution obtained on a given grid to a coarser grid and calculating the image of the discrete Navier–Stokes operator of the coarse grid on the restricted velocity and pressure field. The proposed methodology is not a new concept but its application to colocated finite volume discretizations of the incompressible Navier–Stokes equations is made possible by the introduction of a variant of the momentum interpolation technique for mass fluxes where the pressure part of the mass fluxes is not dependent on the coefficients of the linearized momentum equations. The theory presented is supported by a number of numerical experiments. The methodology is developed for two‐dimensional flows, but extension to three‐dimensional cases should not pose problems. Copyright © 2005 John Wiley & Sons, Ltd.