Visualizing diffusion tensor fields on streamsurfaces with merging ellipsoids

Streamsurfaces in diffusion tensor fields are used to represent structures with primarily planar diffusion. So far, however, no effort has been made on the visualization of the anisotropy of diffusion on them, although this information is very important to identify the problematic regions of these structures. We propose two methods to display this anisotropy information. The first one employs a set of merging ellipsoids, which simultaneously characterize the local tensor details- anisotropy- on them and portray the shape of the streamsurfaces.The weight between the streamsurfaces continuity and the discrete local tensors can be interactively adjusted by changing some given parameters. The second one generates a dense LIC(line integral convolution) texture of the two tangent eigenvector fields along the streamsurfaces firstly, and then blends in some color mapping indicating the anisotropy information. For high speed and high quality of texture images, we confine both the generation and the advection of the LIC texture in the image space. Merging ellipsoids method reveals the entire anisotropy information at discrete points by exploiting the geometric attribute of ellipsoids, and thus suits for local and detailed examination of the anisotropy; the texture-based method gives a global representation of the anisotropy on the whole streamsurfaces with texture and color attributes.To reveal the anisotropy information more efficiently, we integrate the two methods and use them at two different levels of details.     

Some classes of Kenmotsu manifolds with respect to semi-symmetric metric connection

In this paper, we study conharmonic curvature tensor in Kenmotsu manifolds with respect to semi-symmetric metric connection and also characterize conharmonically flat, conharmonically semi-symmetric and φ-conharmonically flat Kenmotsu manifolds with respect to semi-symmetric metric connection.     

ENERGY ESTIMATES FOR DELAY DIFFUSION-REACTION EQUATIONS

In this paper we consider nonlinear delay diffusion-reaction equations with initial and Dirichlet boundary conditions. The behaviour and the stability of the solution of such initial boundary value problems （IBVPs） are studied using the energy method. Simple numerical methods are considered for the computation of numerical approximations to the solution of the nonlinear IBVPs. Using the discrete energy method we study the stability and convergence of the numerical approximations. Numerical experiments are carried out to illustrate our theoretical results.     

On pricing of corporate securities in the case of jump-diffusion

Structural models of credit risk are known to present vanishing spreads at very short maturitiesThis shortcoming, which is due to the diffusive behavior assumed for asset values, can be circumvented by considering discontinuities of the jump type in their evolution over timeIn this paper, we extend the pricing model for corporate bond and determine the default probability in jump-diffusion model to address this issueTo make the problem clearly,we first investigate the case that the firm value follows a geometric Brownian motion under similar assumptions to those in Black and Scholes(1973), Briys and de Varenne(1997), i.e, the default barrier is KD(t, T) and the recovery rate is(1- ω), where D(t, T) is the price of zero coupon default free bond and ω is a constant(0 ω≤ 1)By changing the numeraire, we obtain the closed-form solution for both the price of bond and default probabilityFurther, we consider the case of jump-diffusion and suppose that a firm will go bankruptcy if its value Vt ≤ KD(t, T)and at the same time, the bondholder will receive(1- ω)Vt KBy introducing the Green function of PDE with absorbing boundary and converting the problem to an II-type Volterra integral equation, we get the closed-form expressions in series form for bond price and corresponding default probabilityNumerical results are presented to show the impact of different parameters to credit spread of bond.     

VARIATIONAL DISCRETIZATION FOR OPTIMAL CONTROL GOVERNED BY CONVECTION DOMINATED DIFFUSION EQUATIONS

In this paper, we study variational discretization for the constrained optimal control problem governed by convection dominated diffusion equations, where the state equation is approximated by the edge stabilization Galerkin method. A priori error estimates are derived for the state, the adjoint state and the control. Moreover, residual type a posteriori error estimates in the L^2-norm are obtained. Finally, two numerical experiments are presented to illustrate the theoretical results.     

A Class of Nonlinear Degenerate Parabolic Equations Not in Divergence Form

This paper is devoted to the study of a class of singular nonlinear diffusion problem. The existence and uniqueness of solutions are obtained. Moreover, some properties of solutions such as blow-up property etc. are also discussed.     

Necessary Maximum Principle of Stochastic Optimal Control with Delay and Jump Diffusion

In this paper, we have studied the necessary maximum principle of stochastic optimal control problem with delay and jump diffusion.     

EXISTENCE OF PERIODIC SOLUTIONS TO A SYSTEM WITH FUNCTIONAL RESPONSE ON TIME SCALES

This paper investigates the existence of periodic solutions of a three-species food-chain diffusive system with Beddington-DeAngelis functional responses and time delays in a two-patch environment on time scales. By using a continuation theorem based on coincidence degree theory, we obtain sufficient criteria for the existence of periodic solutions for the system. Moreover, when the time scale T is chosen as R or Z, the existence of the periodic solutions of the corresponding continuous and discrete models follows. Therefore, the method is unified to provide the existence of the desired solutions for continuous differential equations and discrete difference equations.     

A POSTERIORI ENERGY-NORM ERROR ESTIMATES FOR ADVECTION-DIFFUSION EQUATIONS APPROXIMATED BY WEIGHTED INTERIOR PENALTY METHODS

We propose and analyze a posteriori energy-norm error estimates for weighted interior penalty discontinuous Galerkin approximations of advection-diffusion-reaction equations with heterogeneous and anisotropic diffusion. The weights, which play a key role in the analysis, depend on the diffusion tensor and are used to formulate the consistency terms in the discontinuous Galerkin method. The error upper bounds, in which all the constants are specified, consist of three terms： a residual estimator which depends only on the elementwise fluctuation of the discrete solution residual, a diffusive flux estimator where the weights used in the method enter explicitly, and a non-conforming estimator which is nonzero because of the use of discontinuous finite element spaces. The three estimators can be bounded locally by the approximation error. A particular attention is given to the dependency on problem parameters of the constants in the local lower error bounds. For moderate advection, it is shown that full robustness with respect to diffusion heterogeneities is achieved owing to the specific design of the weights in the discontinuous Galerkin method, while diffusion anisotropies remain purely local and impact the constants through the square root of the condition number of the diffusion tensor. For dominant advection, it is shown, in the spirit of previous work by Verfiirth on continuous finite elements, that the local lower error bounds can be written with constants involving a cut-off for the ratio of local mesh size to the reciprocal of the square root of the lowest local eignevalue of the diffusion tensor.     

A NOTE ON n-PERINORMAL OPERATORS

The study of operators satisfying σja(T) = σa(T) is of significant interest. Does σja(T) = σa(T) for n-perinormal operator T ∈ B(H)ff This question was raised by Mecheri and Braha [Oper. Matrices 6(2012), 725–734]. In the note we construct a counterexample to this question and obtain the following result: if T is a n-perinormal operator in B(H), then σja(T)\{0} = σa(T)\{0}. We also consider tensor product of n-perinormal operators.     

PREDUAL SPACES FOR Q SPACES

To find the predual spaces Pα（R^n） of Qα（R^n） is an important motivation in the study of Q spaces. In this article, wavelet methods are used to solve this problem in a constructive way. First, an wavelet tent atomic characterization of Pα（Rn） is given, then its usual atomic characterization and Poisson extension characterization are given. Finally, the continuity on Pα of Calderon-Zygmund operators is studied, and the result can be also applied to give the Morrey characterization of Pα（Rn）.     

STABILIZED FEM FOR CONVECTION-DIFFUSION PROBLEMS ON LAYER-ADAPTED MESHES

The application of a standard Galerkin finite element method for convection-diffusion problems leads to oscillations in the discrete solution, therefore stabilization seems to be necessary. We discuss several recent stabilization methods, especially its combination with a Galerkin method on layer-adapted meshes. Supercloseness results obtained allow an improvement of the discrete solution using recovery techniques.     

A TWO-SCALE HIGHER-ORDER FINITE ELEMENT DISCRETIZATION FOR SCHROEDINGER EQUATION

In this paper, a two-scale higher-order finite element discretization scheme is proposed and analyzed for a Schroedinger equation on tensor product domains. With the scheme, the solution of the eigenvalue problem on a fine grid can be reduced to an eigenvalue problem on a much coarser grid together with some eigenvalue problems on partially fine grids. It is shown theoretically and numerically that the proposed two-scale higher-order scheme not only significantly reduces the number of degrees of freedom but also produces very accurate approximations.     

FINITE ELEMENTS WITH LOCAL PROJECTION STABILIZATION FOR INCOMPRESSIBLE FLOW PROBLEMS

In this paper we review recent developments in the analysis of finite element methods for incompressible flow problems with local projection stabilization （LPS）. These methods preserve the favourable stability and approximation properties of classical residual-based stabilization （RBS） techniques but avoid the strong coupling of velocity and pressure in the stabilization terms. LPS-methods belong to the class of symmetric stabilization techniques and may be characterized as variational multiscale methods. In this work we summarize the most important a priori estimates of this class of stabilization schemes developed in the past 6 years. We consider the Stokes equations, the Oseen linearization and the NavierStokes equations. Furthermore, we apply it to optimal control problems with linear（ized） flow problems, since the symmetry of the stabilization leads to the nice feature that the operations ＂discretize＂ and ＂optimize＂ commute.     

AN EFFICIENT METHOD FOR MULTIOBJECTIVE OPTIMAL CONTROL AND OPTIMAL CONTROL SUBJECT TO INTEGRAL CONSTRAINTS

We introduce a new and efficient numerical method for multicriterion optimal control and single criterion optimal control under integral constraints. The approach is based on extending the state space to include information on a ＂budget＂ remaining to satisfy each constraint; the augmented Hamilton-Jacobi-Bellman PDE is then solved numerically. The efficiency of our approach hinges on the causality in that PDE, i.e., the monotonicity of characteristic curves in one of the newly added dimensions. A semi-Lagrangian ＂marching＂ method is used to approximate the discontinuous viscosity solution efficiently. We compare this to a recently introduced ＂weighted sum＂ based algorithm for the same problem [25]. We illustrate our method using examples from flight path planning and robotic navigation in the presence of friendly and adversarial observers.     

DIRECT MINIMIZATION FOR CALCULATING INVARIANT SUBSPACES IN DENSITY FUNCTIONAL COMPUTATIONS OF THE ELECTRONIC STRUCTURE

In this article, we analyse three related preconditioned steepest descent algorithms, which are partially popular in Hartree-Fock and Kohn-Sham theory as well as invariant subspace computations, from the viewpoint of minimization of the corresponding functionals, constrained by orthogonality conditions. We exploit the geometry of the admissible manifold, i.e., the invariance with respect to unitary transformations, to reformulate the problem on the Grassmann manifold as the admissible set. We then prove asymptotical linear convergence of the algorithms under the condition that the Hessian of the corresponding Lagrangian is elliptic on the tangent space of the Grassmann manifold at the minimizer.     

Cartesian closed categories of FƵ-domains

A subset system Z assigns to each partially ordered set P a certain collection Z（P） of subsets. In this paper, a new kind of subset systems called directable subset systems is introduced. For a directable subset system Z, the concepts of FZ-way-below relation and FZ-domain are introduced. The well-known Scott topology is naturally generalized to the Z-level and the resulting topology is called FZ-Scott topology, and the continuous functions with respect to this topology are characterized by preserving the suprema of directed Z-sets. Then, we mainly consider a generalization of the cartesian closedness of the categories DCPO of directed complete posets, BF of bifinite domains and FS of FS-domains to the Z-level. Corresponding to them, it is proved that, for a suitable subset system Z, the categories FZCPO of Z-complete posets, FSFZ of finitely separated FZ-domains and BFFZ of bifinite FZ-domains are all cartesian closed. Some examples of these categories are given.     

RECENT PROGRESS ON SPHERICAL HARMONIC APPROXIMATION MADE BY BNU RESEARCH GROUP --In memory of Professor Sun Yongsheng

As early as in 1990, Professor Sun Yongsheng, suggested his students at Beijing Normal University to consider research problems on the unit sphere. Under his guidance and encouragement his students started the research on spherical harmonic analysis and approximation. In this paper, we incompletely introduce the main achievements in this area obtained by our group and relative researchers during recent 5 years （2001-2005）. The main topics are： convergence of Cesaro summability, a.e. and strong summability of Fourier-Laplace series; smoothness and K-functionals; Kolmogorov and linear widths.     

FINITE ELEMENT APPROXIMATIONS FOR SCHROEDINGER EQUATIONS WITH APPLICATIONS TO ELECTRONIC STRUCTURE COMPUTATIONS

In this paper, both the standard finite element discretization and a two-scale finite element discretization for SchrSdinger equations are studied. The numerical analysis is based on the regularity that is also obtained in this paper for the Schroedinger equations. Very satisfying applications to electronic structure computations are provided, too.     

Projections,Birkhoff Orthogonality and Angles in Normed Spaces

Let X be a Minkowski plane, i.e., a real two dimensional normed linear space. We use projections to give a definition of the angle Aq（x, y） between two vectors x and y in X, such that x is Birkhoff orthogonal to y if and only if Aq（x,y）=π/2. Some other properties of this angle are also discussed.