Pseudotensors in asymptotically curvilinear coordinates

We show how to calculate pseudotensor-based conserved quantities for isolated systems in general relativity, in a way which allows an arbitrary asymptotic behavior of the coordinate system used. Our method is a generalization of that given by Persides [1], and allows the asymptotic evaluation of energy, momentum, and angular momentum in any coordinate system. We carry out the calculation for the Schutz-Sorkin gravitational Noether operator, which is a pseudotensorial operator on vector fields that reduces to the familiar pseudotensors for particular choices of the fields.     

Quantization in generalized coordinates III-Lagrangian formulation

It is shown that if one incorporates the generalized coordinate quantum velocitiesQ ¹ as given byQ ¹=l[H,Q ¹](h=1) into the generalized classical Lagrangian for a free particle (the total energy),L=1/2Q ¹ g _tk Q ^k one does not obtain (no matter what ordering of the operatorsq ^l ,q ^k andg _lkwe choose the correct quantum Lagrangian operator which is a transformation from -1/2V² to generalized coordinates (Gruber, 1971, 1972).q ^l as given byq ^l=i[H,q ^l] turns out to be the Hermitian part of a more generaiized operator which we call the total generalized velocity operator similar to the notation in ear previous articles (Gruber, 1971, 1972). This total velocity operator really determines the fundamental structure governing our system in the Lagrangian formulation. We show that ft is through the total velocity operator that we make the transition from classical to quantum mechanics and through our procedure we arrive at the correct quantum Lagrangian operator.     

正交曲线坐标系中加速度的矢量求法

利用矢量法,给出一种求解正交曲线坐标系中加速度的简捷方法。     

Stabilized non-dissipative approximations of Euler equations in generalized curvilinear coordinates

We discuss stabilization strategies for finite-difference approximations of the compressible Euler equations in generalized curvilinear coordinates that do not rely on explicit upwinding or filtering of the physical variables. Our approach rather relies on a skew-symmetric-like splitting of the convective derivatives, that guarantees preservation of kinetic energy in the semi-discrete, low-Mach-number limit. A locally conservative formulation allows efficient implementation and easy incorporation into existing compressible flow solvers. The validity of the approach is tested for benchmark flow cases, including the propagation of a cylindrical vortex, and the head-on collision of two vortex dipoles. The tests support high accuracy and superior stability over conventional central discretization of the convective derivatives. The potential use for DNS/LES of turbulent compressible flows in complex geometries is discussed.     

Formulations of the 3+1 evolution equations in curvilinear coordinates

Following Brown (Phys Rev D79:104029, 2009), in this paper we give an overview of how to modify standard hyperbolic formulations of the 3+1 evolution equations of General Relativity in such a way that all auxiliary quantities are true tensors, thus allowing for these formulations to be used with curvilinear sets of coordinates such as spherical or cylindrical coordinates. After considering the general case for both the Nagy–Ortiz–Reula and the Baumgarte–Shapiro–Shibata–Nakamura (BSSN) formulations, we specialize to the case of spherical symmetry and also discuss the issue of regularity at the origin. Finally, we show some numerical examples of the modified BSSN formulation at work in spherical symmetry.     

BiGlobal stability analysis in curvilinear coordinates of massively separated lifting bodies

A methodology based on spectral collocation numerical methods for global flow stability analysis of incompressible external flows is presented. A potential shortcoming of spectral methods, namely the handling of the complex geometries encountered in global stability analysis, has been dealt with successfully in past works by the development of spectral-element methods on unstructured meshes. The present contribution shows that a certain degree of regularity of the geometry may be exploited in order to build a global stability analysis approach based on a regular spectral rectangular grid in curvilinear coordinates and conformal mappings. The derivation of the stability linear operator in curvilinear coordinates is presented along with the discretisation method. Unlike common practice to the solution of the same problem, the matrix discretising the eigenvalue problem is formed and stored. Subspace iteration and massive parallelisation are used in order to recover a wide window of its leading Ritz system. The method is applied to two external flows, both of which are lifting bodies with separation occurring just downstream of the leading edge. Specifically the flow configurations are a NACA 0015 airfoil, and an ellipse of aspect ratio 8 chosen to closely approximate the geometry of the airfoil. Both flow configurations are at an angle of attack of 18° with a Reynolds number based on the chord length of 200. The results of the stability analysis for both geometries are presented and illustrate analogous features.     

New interpolation technique for the CIP method on curvilinear coordinates

We propose a new interpolation technique for the CIP method applied to curvilinear coordinates. The CIP method can hardly maintain third-order accuracy on curvilinear coordinates. The reason for the degeneracy in accuracies has not been discussed in detail. This paper reveals the problems of the CIP method on curvilinear coordinates and presents an improved CIP method to solve the advection equation accurately. The features of the presented method are: (1) the metric computation on the upwind stencil is defined in the same manner as in the advection phase of the CIP method; and (2) gradient values in the physical domain in the computation on the curvilinear coordinates are used. Various test problems show that the improved CIP method has approximate third-order accuracy.     

Employment of curvilinear coordinates in ab initio calculations of insulators using pseudopotentials

The standard ab initio scheme for calculating the structure of crystals using nonlocal pseudopotentials is modified for use in curvilinear coordinates. A method for solving the Poisson equation for the Coulomb potential in a curved space in the k representation is found. It is shown in the example of calculations for crystals of insulators having an NaCl structure that the employment of a curved space permits a very significant decrease in the required size of the basis set. Fiz. Tverd. Tela (St. Petersburg) 41, 241–246 (February 1999)     

Canonical quantization of classical mechanics in curvilinear coordinates. Invariant quantization procedure

In the paper is presented an invariant quantization procedure of classical mechanics on the phase space over flat configuration space. Then, the passage to an operator representation of quantum mechanics in a Hilbert space over configuration space is derived. An explicit form of position and momentum operators as well as their appropriate ordering in arbitrary curvilinear coordinates is demonstrated. Finally, the extension of presented formalism onto non-flat case and related ambiguities of the process of quantization are discussed.     

正交曲线坐标系在变化地形声场计算中的应用*

变化地形声场计算是越来越经常碰到的情况。变化地形往往给声场计算带来计算量增大、求解精度下降的问题。坐标系选择是声场计算中非常重要的问题,本文研究基于地形的合适正交曲线坐标系选取规则。在此坐标系中求解Helmholtz方程,能够极大简化求解过程、提高计算精度。结合简正波理论,从新的角度展示了正交曲线坐标系在某些典型变化地形线源声场计算问题中的应用。结果表明,本文方法在变化地形处理方面,能克服传统阶梯步进法的不足。     

Coupling of Dirichlet-to-Neumann boundary condition and finite difference methods in curvilinear coordinates for multiple scattering

The applicability of the Dirichlet-to-Neumann technique coupled with finite difference methods is enhanced by extending it to multiple scattering from obstacles of arbitrary shape. The original boundary value problem (BVP) for the multiple scattering problem is reformulated as an interface BVP. A heterogenous medium with variable physical properties in the vicinity of the obstacles is considered. A rigorous proof of the equivalence between these two problems for smooth interfaces in two and three dimensions for any finite number of obstacles is given. The problem is written in terms of generalized curvilinear coordinates inside the computational region. Then, novel elliptic grids conforming to complex geometrical configurations of several two-dimensional obstacles are constructed and approximations of the scattered field supported by them are obtained. The numerical method developed is validated by comparing the approximate and exact far-field patterns for the scattering from two circular obstacles. In this case, for a second order finite difference scheme, a second order convergence of the numerical solution to the exact solution is easily verified.     

A three-dimensional solid-liquid two-phase turbulence model with the effect of vegetation in non-orthogonal curvilinear coordinates

A three-dimensional k-ε-Ap solid-liquid two-phase two-fluid model with the effect of vegetation is solved numerically with a finite-volume method on an adaptive grid to study water-sediment movements and bed evolution in vegetated channels. The additional drag force and additional turbulence generation due to vegetation are added to the relevant control equations for simulating the interaction between vegetation and flow. The flow structure and the bed-topography changes in a 60° partly vegetated channel be...     

Natural curvilinear coordinates for ideal MHD equations. Non-stationary flows with constant total pressure

Equations of magnetohydrodynamics (MHD) in the natural curvilinear system of coordinates where trajectories and magnetic lines play a role of coordinate curves are reduced to the non-linear vector wave equation coupled with the incompressibility condition in the form of the generalized Cauchy integral. The symmetry group of obtained equation, equivalence transformation, and group classification with respect to the constitutive equation are calculated. New exact solutions with functional arbitrariness describing non-stationary incompressible flows with constant total pressure are given by explicit formulae. The corresponding magnetic surfaces have the shape of deformed nested cylinders, tori, or knotted tubes.