Least-squares Galerkin procedures for parabolic integro-differential equations

Two least-squares Galerkin finite element schemes are formulated to solve parabolic integro-differential equations. The advantage of this method is that it is not subject to the LBB condition. The convergence analysis shows that the least-squares mixed element schemes yield the approximate solution with optimal accuracy in H(div;Ω)×H¹(Ω) and (L²(Ω))²×L²(Ω), respectively.     

Asymptotic and Partial Asymptotic Hankel Operators on $$H^2(\mathbb{D}^n)$$ H 2 ( D n ) )

In this paper, we generalize the concept of asymptotic Hankel operators on H²(D) to the Hardy space H²(Dⁿ) (over polydisk) in terms of asymptotic Hankel and partial asymptotic Hankel operators and investigate some properties in case of its weak and strong convergence. Meanwhile, we introduce i^th-partial Hankel operators on H²(Dⁿ) and obtain a characterization of its compactness for n > 1. Our main results include the containment of Toeplitz algebra in the collection of all strong partial asymptotic Hankel operators on H²(Dⁿ). It is also shown that a Toeplitz operator with symbol φ is asymptotic Hankel if and only if φ is holomorphic function in L^∞(Tⁿ).     

Cohomology of the Universal Enveloping Algebras of Certain Bigraded Lie Algebras

Let p be an odd prime and q = 2(p-1).Up to total degree t-s max{(5p~3+ 6p~2+ 6 p +4)q-10,p~4q},the generators of H~(s,t)(U(L)),the cohomology of the universal enveloping algebra of a bigraded Lie algebra L,are determined and their convergence is also verified.Furthermore our results reveal that this cohomology satisfies an analogous Poinare duality property.This largely generalizes an earlier classical results due to J.P.May.     

On the relation between interior critical points and parameters for a class of nonlinear problems with Neumann–Robin boundary conditions

We consider boundary value problem

where   0, λ > 0 are parameters and f  C²[0, ∞) such that f(0) < 0. In this paper we study for the cases p  (0, β) and p  (β, θ) (p is the value of the solution at x = 0 and β, θ are such that f(β) = 0, , the relation between λ and the number of interior critical points of the positive solutions of the above system.     

On the extreme eigenvalues of hermitian (block) toeplitz matrices

We are concerned with the behavior of the minimum (maximum) eigenvalue λ₀⁽ⁿ⁾ (λ_n⁽ⁿ⁾) of an (n + 1) × (n + 1) Hermitian Toeplitz matrix T_n(ƒ) where ƒ is an integrable real-valued function. Kac, Murdoch, and Szegö, Widom, Parter, and R. H. Chan obtained that λ₀⁽ⁿ⁾ — min ƒ = O(1/n^2k) in the case where ƒ C^2k, at least locally, and ƒ — inf ƒ has a zero of order 2k. We obtain the same result under the second hypothesis alone. Moreover we develop a new tool in order to estimate the extreme eigenvalues of the mentioned matrices, proving that the rate of convergence of λ₀⁽ⁿ⁾ to inf ƒ depends only on the order ρ (not necessarily even or integer or finite) of the zero of ƒ — inf ƒ. With the help of this tool, we derive an absolute lower bound for the minimal eigenvalues of Toeplitz matrices generated by nonnegative L¹ functions and also an upper bound for the associated Euclidean condition numbers. Finally, these results are extended to the case of Hermitian block Toeplitz matrices with Toeplitz blocks generated by a bivariate integrable function ƒ.     

Stable limits of empirical processes of moving averages with infinite variance

The paper obtains a functional limit theorem for the empirical process of a stationary moving average process X_t with i.i.d. innovations belonging to the domain of attraction of a symmetric -stable law, 1<<2, with weights b_j decaying as j^−β, 1<β<2/. We show that the empirical process (normalized by N^1/β) weakly converges, as the sample size N increases, to the process c_x⁺L⁺+c_x⁻L⁻, where L⁺,L⁻ are independent totally skewed β-stable random variables, and c_x⁺,c_x⁻ are some deterministic functions. We also show that, for any bounded function H, the weak limit of suitably normalized partial sums of H(X_s) is an β-stable Lévy process with independent increments. This limiting behavior is quite different from the behavior of the corresponding empirical processes in the parameter regions 1/<β<1 and 2/<β studied in Koul and Surgailis (Stochastic Process. Appl. 91 (2001) 309) and Hsing (Ann. Probab. 27 (1999) 1579), respectively.     

Choosability conjectures and multicircuits

This paper starts with a discussion of several old and new conjectures about choosability in graphs. In particular, the list-colouring conjecture, that ch′=χ′ for every multigraph, is shown to imply that if a line graph is (a : b)-choosable, then it is (ta : tb)-choosable for every positive integer t. It is proved that ch(H²)=χ(H²) for many “small” graphs H, including inflations of all circuits (connected 2-regular graphs) with length at most 11 except possibly length 9; and that ch″(C)=χ″(C) (the total chromatic number) for various multicircuits C, mainly of even order, where a multicircuit is a multigraph whose underlying simple graph is a circuit. In consequence, it is shown that if any of the corresponding graphs H² or T(C) is (a : b)-choosable, then it is (ta : tb)-choosable for every positive integer t.     

Irreducible components of L2 functions on hopf bundles

This paper provides an explicit decomposition of the L² function space on the unit sphereS^dn-1 for d = 1,2 and 4, into irreducible representations under the action of the Lie Groups K= SO(n) × SO(1)S(U(n) × U(1)), and Sp(n) × Sp(l), respectively. The decomposition is realized as the eigenspaces of the Laplacian acting on homogeneous polynomials over the reals, complex numbers and quaternions. For the quaternionic case, an additional differential operator that commutes with the Laplacian is used to find the decomposition.     

Phase Retrieval of Real-valued Functions in Sobolev Space

The Sobolev space H^s(R^d) with s > d/2 contains many important functions such as the bandlimited or rational ones. In this paper we propose a sequence of measurement functions {ϕ_j,k^γ} ⊆ H^-s(R^d) to the phase retrieval problem for the real-valued functions in H^s(R^d). We prove that any real-valued function f ∈ H^s(R^d) can be determined, up to a global sign, by the phaseless measurements {|<f, ϕ_j,k^γ>|}. It is known that phase retrieval is unstable in infinite dimensional spaces with respect to perturbations of the measurement functions. We examine a special type of perturbations that ensures the stability for the phase-retrieval problem for all the real-valued functions in H^s(R^d) ∩ C¹(R^d), and prove that our iterated reconstruction procedure guarantees uniform convergence for any function f ∈ H^s(R^d) ∩ C¹(R^d) whose Fourier transform f is L¹-integrable. Moreover, numerical simulations are conducted to test the efficiency of the reconstruction algorithm.     

Computational analysis of supersonic turbulent boundary layers over rough surfaces using the k–ω and the stress–ω models

The influence of surface roughness in the prediction of the mean flow and turbulent properties of a high-speed supersonic (M = 2.7, Re/m = 2 × 10⁷) turbulent boundary layer flow over a flat plate is numerically investigated. In particular, the performance of the k–ω and stress–ω turbulence models is evaluated against the available experimental data. Even though the performance of these models have been proven satisfactory in the computation of incompressible boundary layer flow over rough surfaces, their validity for high-speed compressible has not been investigated yet. It is observed from this study that, for smooth surface, both k–ω and stress–ω models perform very well in predicting the mean flow and turbulence quantities in supersonic flow. For rough surfaces, both models matched the experimental data fairly well for lower roughness heights but performed unsatisfactorily for higher roughness conditions. Overall the performance of the k–ω model is better than the stress–ω model. The stress–ω model does not show any strong advantages to make up for the extra computational cost associated with it. The predictions indicate that the ω boundary conditions at the wall in both models, especially the stress–ω model, need to be refined and reconsidered to include the geometric factor for supersonic flow over surfaces with large roughness values.     

New constructions of divisible designs

A construction is given for a (p^2a(p+1),p²,p^2a+1(p+1),p^2a+1,p^2a(p+1)) (p a prime) divisible difference set in the group H×Z²_pa+1 where H is any abelian group of order p+1. This can be used to generate a symmetric semi-regular divisible design; this is a new set of parameters for λ₁≠0, and those are fairly rare. We also give a construction for a (p^a−1+p^a−2+…+p+2,p^a+2, p^a(p^a+p^a−1+…+p+1), p^a(p^a−1+…+p+1), p^a−1(p^a+…+p²+2)) divisible difference set in the group H×Z_p2×Z^a_p. This is another new set of parameters, and it corresponds to a symmetric regular divisible design. For p=2, these parameters have λ₁=λ₂, and this corresponds to the parameters for the ordinary Menon difference sets.     

二阶双曲方程的全离散格式下的混合元超收敛分析

通过在空间方向上使用双线性元和最低阶的Nedelec元(即Q₁₁+Q₀₁×Q₁₀)以及在时间方向上使用二阶精度的数值逼近格式,得到了在矩形网格上二阶双曲方程全离散混合元格式下的对原始变量的L^∞(H¹)和流量的L^∞((L²)²)的超逼近和超收敛的误差结果.在分析过程中,巧妙地使用了上述混合单元对在矩形网格上的特有的高精度积分恒等式和精确解的投影和插值之间的在H¹范数意义下的超逼近的估计.最后,给出一些数值结果来验证理论分析的正确性.     

The thermal equilibrium state of semiconductor devices

The thermal equilibrium state of two oppositely charged gases confined to a bounded domain , m = 1,2 or m = 3, is entirely described by the gases' particle densities p, n minimizing the total energy (p, n). it is shown that for given P, N > 0 the energy functional admits a unique minimizer in {(p, n) ε L²(Ω) x L ²(Ω) : p, n ≥ 0, _Ωp = P, _Ωn = N} and that p, n ε C(Ω) ∩ L^∞(Ω).

The analysis is applied to the hydrodynamic semiconductor device equations. These equations in general possess more than one thermal equilibrium solution, but only the unique solution of the corresponding variational problem minimizes the total energy. It is equivalent to prescribe boundary data for electrostatic potential and particle densities satisfying the usual compatibility relations and to prescribe V^e and P, N for the variational problem.     

Large 2P3-free graphs with bounded degree

Let ex^* (D;H) be the maximum number of edges in a connected graph with maximum degree D and no induced subgraph H; this is finite if and only if H is a disjoint union of paths. If the largest component of such an H has order m, then ex^*(D; H) = O(D²ex^*(D; P_m)). Constructively, ex^*(D;qP_m) = Θ(gD²ex^*(D;P_m)) if q>1 and m> 2(Θ(gD²) if m = 2). For H = 2P₃ (and D 8), the maximum number of edges is if D is even and

Full-size image

if D is odd, achieved by a unique extremal graph.     

Numerical approximation of the product of the square root of a matrix with a vector

Given an n×n symmetric positive definite matrix A and a vector , two numerical methods for approximating are developed, analyzed, and computationally tested. The first method applies a Newton iteration to a specific nonlinear system to approximate while the second method applies a step-control method to numerically solve a specific initial-value problem to approximate . Assuming that A is first reduced to tridiagonal form, the first method requires O(n²) operations per iteration while the second method requires O(n) operations per iteration. In contrast, numerical methods that first approximate A^1/2 and then compute generally require O(n³) operations per iteration.     

Fractional Fourier transform on R2 and an application

We focus on the L^p(R²) theory of the fractional Fourier transform (FRFT) for 1 ≤ p ≤ 2. In L¹(R²), we mainly study the properties of the FRFT via introducing the two-parameter chirp operator. In order to get the point-wise convergence for the inverse FRFT, we introduce the fractional convolution and establish the corresponding approximate identities. Then the well-defined inverse FRFT is given via approximation by suitable means, such as fractional Gauss means and Able means. Furthermore, if the signal F_α,βf is received, we give the process of recovering the original signal f with MATLAB. In L²(R²), the general Plancherel theorem, direct sum decomposition, and the general Heisenberg inequality for the FRFT are obtained.     

Non blow-up criterion for the 3-D Magneto-hydrodynamics equations in the limiting case

In this paper, we prove that suitable weak solution(u, b) of the 3-D MHD equations can be extended beyond T if u∈L~∞(0,T; L~3(R~3)) and the horizontal components b_h of the magnetic field satisfies the well-known Ladyzhenskaya-Prodi-Serrin condition, which improves the corresponding regularity criterion by Mahalov-Nicolaenko-Shilkin.     

Stability properties of non-negative solutions to a non-autonomous p-Laplacian equation

We study the stability of non-negative stationary solutions of

where Δ_p denotes the p-Laplacian operator defined by Δ_pz = div(z^p−2z); p > 2, Ω is a bounded domain in R^N(N  1) with smooth boundary where [0,1],h:∂Ω→R⁺ with h = 1 when  = 1, λ > 0, and g:Ω×[0,∞)→R is a continuous function. If g(x, u)/u^p−1 be strictly increasing (decreasing), we provide a simple proof to establish that every non-trivial non-negative solution is unstable (stable).     

Numerical method for solving 3D inverse scattering problems

Let q(x) L²(D), D R³ is a bounded domain, q = 0 outside D, q is real-valued. Assume that A(\Gj;\t';,\Gj;,k) A(\Gj;\t';,\Gj), the scattering amplitude, is known for all \Gj;|t',\Gj; S², S² is the unit sphere, an d a fixed k \r>0. These data determine q(x) uniquely and a numerical method is given for computing q(x).     

The smallest degree sum that yields potentially kCℓ-graphic sequences

Gould et al. (Combinatorics, Graph Theory and Algorithms, Vol. 1, 1999, pp. 387–400) considered a variation of the classical Turán-type extremal problems as follows: For a given graph H, determine the smallest even integer σ(H,n) such that every n-term graphic sequence π=(d₁,d₂,…,d_n) with term sum σ(π)=d₁+d₂++d_nσ(H,n) has a realization G containing H as a subgraph. In this paper, for given integers k and ℓ, ℓ7 and 3kℓ, we completely determine the smallest even integer σ(_kC_ℓ,n) such that each n-term graphic sequence π=(d₁,d₂,…,d_n) with term sum σ(π)=d₁+d₂++d_nσ(_kC_ℓ,n) has a realization G containing a cycle of length r for each r, krℓ.