On the spectral problem for anyons

We consider the spectral problem resulting from the Schrödinger equation for a quantum system ofn2 indistinguishable, spinless, hard-core particles on a domain in two dimensional Euclidian space. For particles obeying fractional statistics, and interacting via a repulsive hard core potential, we provide a rigorous framework for analysing the spectral problem with its multi-valued wave functions.Partially supported by the Mathematical Sciences Research Institute, Berkeley California, under NSF Grant # DMS 8505550Partially supported under NSF Grant no. DMR-9101542     

Exponential convergence of a spectral projection of the KdV equation

It is shown that a spectral approximation of the Korteweg–de Vries equation converges exponentially fast to the true solution if the Fourier basis is used and if the solution is analytic in a fixed strip about the real axis. Computations are carried out which show that the exponential convergence rate can be achieved in practice.     

A numerical method for simulating the dynamics of 3D axisymmetric vesicles suspended in viscous flows

We extend [Shravan K. Veerapaneni, Denis Gueyffier, Denis Zorin, George Biros, A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D, Journal of Computational Physics 228(7) (2009) 2334–2353] to the case of three-dimensional axisymmetric vesicles of spherical or toroidal topology immersed in viscous flows. Although the main components of the algorithm are similar in spirit to the 2D case—spectral approximation in space, semi-implicit time-stepping scheme—the main differences are that the bending and viscous force require new analysis, the linearization for the semi-implicit schemes must be rederived, a fully implicit scheme must be used for the toroidal topology to eliminate a CFL-type restriction and a novel numerical scheme for the evaluation of the 3D Stokes single layer potential on an axisymmetric surface is necessary to speed up the calculations. By introducing these novel components, we obtain a time-scheme that experimentally is unconditionally stable, has low cost per time step, and is third-order accurate in time. We present numerical results to analyze the cost and convergence rates of the scheme. To verify the solver, we compare it to a constrained variational approach to compute equilibrium shapes that does not involve interactions with a viscous fluid. To illustrate the applicability of method, we consider a few vesicle-flow interaction problems: the sedimentation of a vesicle, interactions of one and three vesicles with a background Poiseuille flow.     

An eigen-based high-order expansion basis for structured spectral elements

We present an eigen-based high-order expansion basis for the spectral element approach with structured elements. The new basis exhibits a numerical efficiency significantly superior, in terms of the conditioning of coefficient matrices and the number of iterations to convergence for the conjugate gradient solver, to the commonly-used Jacobi polynomial-based expansion basis. This basis results in extremely sparse mass matrices, and it is very amenable to the diagonal preconditioning. Ample numerical experiments demonstrate that with the new basis and a simple diagonal preconditioner the number of conjugate gradient iterations to convergence has essentially no dependence or only a very weak dependence on the element order. The expansion bases are constructed by a tensor product of a set of special one-dimensional (1D) basis functions. The 1D interior modes are constructed such that the interior mass and stiffness matrices are simultaneously diagonal and have identical condition numbers. The 1D vertex modes are constructed to be orthogonal to all the interior modes. The performance of the new basis has been investigated and compared with other expansion bases.     

On equivalent diffuse conditions for an internal layer in multislab atmospheric radiation

This paper reports on new discrete ordinates conditions for efficiently solving a set of multislab atmospheric radiation problems characterized by an optically stationary internal layer, i.e. an internal layer whose optical (absorption/scattering) properties and optical thickness do not change from one problem to another in the set. The discrete ordinates conditions reported here are founded in a recently developed spectral nodal method for solving multislab atmospheric radiation problems with anisotropic scattering. We suitably use the optically discretized equations of our recently developed spectral nodal method to derive discrete ordinates diffuse conditions, which model the response—diffuse radiation leaving the layer—of an internal layer to an anisotropic inner source and to diffuse radiation that is incident upon the layer at top and bottom. These conditions can be used to replace an optically stationary internal layer in multislab atmospheric radiation computations, while saving computer resources and without degrading the numerical results.     

On the spectral and conservation properties of nonlinear discretization operators

Following the study of Pirozzoli [1], the objective of the present work is to provide a detailed theoretical analysis of the spectral properties and the conservation properties of nonlinear finite difference discretizations. First, a Nonlinear Spectral Analysis (NSA) is proposed in order to study the statistical behavior of the modified wavenumber of a nonlinear finite difference operator, for a large set of synthetic scalar fields with prescribed energy spectrum and random phase. Second, the necessary conditions for local and global conservation of momentum and kinetic energy are derived and verified for nonlinear discretizations. Because the nonlinear mechanisms result in a violation of the energy conservation conditions, the NSA is used to quantify the energy imbalance. Third, the effect of aliasing errors due to the nonlinearity is analyzed. Finally, the theoretical observations are verified for two simple, thought relevant, numerical simulations.     

An analysis of the incompressible viscous flows problem by meshless method

Generally the incompressible viscous flow problem is described by the Navier--Stokes equation. Based on the weighted residual method the discrete formulation of element-free Galerkin is inferred in this paper. By the step-by-step computation in the field of time, and adopting the least-square estimation of the-same-order shift, this paper has calculated both velocity and pressure from the decoupling independent equations. Each time fraction Newton--Raphson iterative method is applied for the velocity and pressure. Finally, this paper puts the method into practice of the shear-drive cavity flow, verifying the validity, high accuracy and stability.     

An analytic and numerical solution with spectral Green's function method for transport equation in spherical geometry

Since in many cases curvilinear geometry is more appropriate than cartesian geometry for precise modeling of the complex systems for reactor calculation, we have developed the spectral Green's function (SGF) method which is employed to obtain angular and scalar flux distributions in heterogeneous sphere geometry with isotropic scattering. In this study, we showed that the neutron transport problems of homogeneous spheres could be reduced to the solution of plane geometry equation.Finally, some results are discussed and compared with those already obtained by diamond difference scheme to test the accuracy of the results. The agreement is satisfactory. SGF method is very suitable for the numerical solution of the neutron transport equation with isotropic scattering.     

Conservation of mass and momentum of the least-squares spectral collocation scheme for the Stokes problem

From the literature it is known that spectral least-squares schemes perform poorly with respect to mass conservation and compensate this lack by a superior conservation of momentum. This should be revised, since the here presented new least-squares spectral collocation scheme leads to an outstanding performance with respect to conservation of momentum and mass. The reasons can be found in using only a few elements, each with high polynomial degree, avoiding normal equations for solving the overdetermined linear systems of equations and by introducing the Clenshaw–Curtis quadrature rule for imposing the average pressure to be zero. Furthermore, we combined the transformation of Gordon and Hall (transfinite mapping) with our least-squares spectral collocation scheme to discretize the internal flow problems.     

On the problem of reality conditions for the Ashtekar formulation

In this paper we discuss the reality conditions for Lorentzian and Euclidean gravity in the Ashtekar formulation by introducing a double conformal transformation.We generalize Marugan‘s results and demonstrate that the values of the double conformal factor have to be either real or double complex numbers.Either Lorentzian or Enclidean gravitational theory is up to the different values of the double conformal factor.Furthermore,the reality conditions of Lorentzian and Euclidean gravitational theory can be expressed in a unified way be use of the double complex function method.     

Characteristic boundary conditions for simulations of compressible reacting flows with multi-dimensional,viscous and reaction effects

A generalized formulation of the characteristic boundary conditions for compressible reacting flows is proposed. The new and improved approach resolves a number of lingering issues of spurious solution behaviour encountered in turbulent reacting flow simulations in the past. This is accomplished (a) by accounting for all the relevant terms in the determination of the characteristic wave amplitudes and (b) by accommodating a relaxation treatment for the transverse gradient terms with the relaxation coefficient properly determined by the low Mach number asymptotic expansion. The new boundary conditions are applied to a comprehensive set of test problems including: vortex-convection; turbulent inflow; ignition front propagation; non-reacting and reacting Poiseuille flows; and counterflow cases. It is demonstrated that the improved boundary conditions perform consistently superior to existing approaches, and result in robust and accurate solutions with minimal acoustic wave interactions at the boundary in hostile turbulent combustion simulation conditions.     

A unifying lifting collocation penalty formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids

Recently a new high-order formulation for 1D conservation laws was developed by Huynh using the idea of “flux reconstruction”. The formulation was capable of unifying several popular methods including the discontinuous Galerkin, staggered-grid multi-domain method, or the spectral difference/spectral volume methods into a single family. The extension of the method to quadrilateral and hexahedral elements is straightforward. In an attempt to extend the method to other element types such as triangular, tetrahedral or prismatic elements, the idea of “flux reconstruction” is generalized into a “lifting collocation penalty” approach. With a judicious selection of solution points and flux points, the approach can be made simple and efficient to implement for mixed grids. In addition, the formulation includes the discontinuous Galerkin, spectral volume and spectral difference methods as special cases. Several test problems are presented to demonstrate the capability of the method.     

High-order compact schemes for incompressible flows: A simple and efficient method with quasi-spectral accuracy

In this paper, a finite difference code for Direct and Large Eddy Simulation (DNS/LES) of incompressible flows is presented. This code is an intermediate tool between fully spectral Navier–Stokes solvers (limited to academic geometry through Fourier or Chebyshev representation) and more versatile codes based on standard numerical schemes (typically only second-order accurate). The interest of high-order schemes is discussed in terms of implementation easiness, computational efficiency and accuracy improvement considered through simplified benchmark problems and practical calculations. The equivalence rules between operations in physical and spectral spaces are efficiently used to solve the Poisson equation introduced by the projection method. It is shown that for the pressure treatment, an accurate Fourier representation can be used for more flexible boundary conditions than periodicity or free-slip. Using the concept of the modified wave number, the incompressibility can be enforced up to the machine accuracy. The benefit offered by this alternative method is found to be very satisfactory, even when a formal second-order error is introduced locally by boundary conditions that are neither periodic nor symmetric. The usefulness of high-order schemes combined with an immersed boundary method (IBM) is also demonstrated despite the second-order accuracy introduced by this wall modelling strategy. In particular, the interest of a partially staggered mesh is exhibited in this specific context. Three-dimensional calculations of transitional and turbulent channel flows emphasize the ability of present high-order schemes to reduce the computational cost for a given accuracy. The main conclusion of this paper is that finite difference schemes with quasi-spectral accuracy can be very efficient for DNS/LES of incompressible flows, while allowing flexibility for the boundary conditions and easiness in the code development. Therefore, this compromise fits particularly well for very high-resolution simulations of turbulent flows with relatively complex geometries without requiring heavy numerical developments.     

Adaptive rational spectral methods for the linear stability analysis of nonlinear fourth-order problems

This paper presents the application of adaptive rational spectral methods to the linear stability analysis of nonlinear fourth-order problems. Our model equation is a phase-field model of infiltration, but the proposed discretization can be directly extended to similar equations arising in thin film flows. The sharpness and structure of the wetting front preclude the use of the standard Chebyshev pseudo-spectral method, due to its slow convergence in problems where the solution has steep internal layers. We discuss the effectiveness and conditioning of the proposed discretization, and show that it allows the computation of accurate traveling waves and eigenvalues for small values of the initial water saturation/film precursor, several orders of magnitude smaller than the values considered previously in analogous stability analyses of thin film flows, using just a few hundred grid points.     

基于谱估计方法的等离子体湍流三波耦合数据处理算法

基于数字谱分析技术求解波耦合方程,进而计算与三波相互作用相关联的线性耦合系数和能量转移,以此开发了数据处理程序用于研究HL-2A 装置等离子体边缘湍流中的非线性能量传递过程。介绍了算法设计和开发的主要思想。应用该程序对与带状流相关的一次放电的实验数据进行了数据处理研究。结果表明,带状流是由等离子体湍流的能量逆级联所驱动的。     

On the analytical solution of viscous fluid flow past a flat plate

In this Letter, we propose a simple approach using HAM to obtain accurate totally analytical solution of viscous fluid flow over a flat plate. First, we show that the solution obtained using HPM is not a reliable one; moreover, we show that HPM is only a special case of HAM and its basic assumptions are restrictive rather than useful. We set ?=−1 for the case of comparison of our results to those obtained using HPM. Afterwards, we introduce an extra auxiliary parameter and a straightforward approach to find best values of this auxiliary parameter which plays a prominent role in the frame of our solution and makes it more convergent in comparison to previous works.     

时间-频率联合分析消除太赫兹光谱干涉

基于Fourier变换的传统太赫兹时域光谱分析技术从不同长度的太赫兹信号中产生了不一致的光谱结果,增加太赫兹时域采样长度又会造成光谱的干涉,不利于太赫兹光谱测量的研究和应用。应用小波变换方法对太赫兹波谱作时间—频率联合分析,将不同时刻的太赫兹波谱展开到二维的时间-频率平面,消除了光谱分析的不一致性和光谱干涉的影响。