非旋转性前进波Euler与Lagrange解间的转换

对等深水中非旋转性的前进重力波动场,以求得的Euler与Lagrange两种形式至第三阶的解,按照同一流体质点在相同时间与位置处其流速唯一与质量守恒性及在自由表面水位处Euler形式解与Lagrange形式解为同一值的特性,来推导二者可相互转换.由连续的Taylor级数展开,考虑波动场中各流体质点的运动轨迹与运动周期,将已知的Euler形式解转换成完全未知的Lagrange形式解,解决了以往成果中出现含时间的不合理的共振项,以及无法得到与Euler系统不同的Lagrange形式的流体质点运动频率与平均运动 关键词： 非旋转性前进波 Euler-Lagrange转换 质点运动轨迹 质点运动频率     

The shallow water equations in Lagrangian coordinates

Recent advances in the collection of Lagrangian data from the ocean and results about the well-posedness of the primitive equations have led to a renewed interest in solving flow equations in Lagrangian coordinates. We do not take the view that solving in Lagrangian coordinates equates to solving on a moving grid that can become twisted or distorted. Rather, the grid in Lagrangian coordinates represents the initial position of particles, and it does not change with time. We apply numerical methods traditionally used to solve differential equations in Eulerian coordinates, to solve the shallow water equations in Lagrangian coordinates. The difficulty with solving in Lagrangian coordinates is that the transformation from Eulerian coordinates results in solving a highly nonlinear partial differential equation. The non-linearity is mainly due to the Jacobian of the coordinate transformation, which is a precise record of how the particles are rotated and stretched. The inverse Jacobian must be calculated, thus Lagrangian coordinates cannot be used in instances where the Jacobian vanishes. For linear (spatial) flows we give an explicit formula for the Jacobian and describe the two situations where the Lagrangian shallow water equations cannot be used because either the Jacobian vanishes or the shallow water assumption is violated. We also prove that linear (in space) steady state solutions of the Lagrangian shallow water equations have Jacobian equal to one. In the situations where the shallow water equations can be solved in Lagrangian coordinates, accurate numerical solutions are found with finite differences, the Chebyshev pseudospectral method, and the fourth order Runge–Kutta method. The numerical results shown here emphasize the need for high order temporal approximations for long time integrations.     

考虑谐波互作用的行波管欧拉非线性理论模型

基于拉格朗日体系的考虑谐波互作用的理论模型,将离散的粒子近似处理为流体,得到电子相位的连续分布函数.对电子相位连续分布函数进行傅里叶一阶展开,并结合贝塞尔母函数关系式,建立了考虑谐波互作用的欧拉非线性理论模型.应用考虑谐波互作用的欧拉非线性理论模型对一支L波段空间行波管和一支C波段空间行波管进行大信号分析,并与拉格朗日理论模型进行对比.结果表明:在增益1dB压缩点之前,考虑谐波互作用的欧拉非线性理论模型与拉格朗日理论模型十分符合,增益最大误差不超过4%.考虑谐波互作用的欧拉非线性理论模型能够有效的对增益1 dB压缩点之前的谐波进行分析.仿真结果验证了考虑谐波互作用的欧拉非线性理论模型的正确性和有效性.考虑谐波互作用的欧拉非线性理论不但提供了一个谐波快速计算模型,而且为后续研究行波管谐波的产生机理与抑制方法奠定了基础.     

非旋转性自由表面重力驻波的Euler与Lagrange方法解间的转换

对于等深水中的非旋转性重力驻波流场,本文用Euler与Lagrange两种方法求得其至三阶的解,根据同一粒流体质点在相同时间与位置处其流速值为唯一与质量守恒及在自由表面水位的Euler形式解与Lagrange形式解相同等特性,来推导其间互可转换.由一系列连续的Taylor级数展开,在考虑波动场中各流体质点的运动轨迹与运动周期条件下,将已知的Euler解转换成完全未知的Lagrange形式解.接着再将所得的Lagrange解转换成对应的Euler形式,均可得到完全相同的结果.由此可得知,在考虑波动场各流体质 关键词： 重力驻波 Euler与Lagrange解间的转换 质点运动轨迹     

Poincaré引力规范场和1/2自旋粒子的运动

木文以Poinaré群作为引力规范群,在有挠率和曲率的空间中,讨论了当引力拉氏量包含场强的线性项与二次项时体系的运动方程,指出球对称真空静引力场方程在“宏观”极限下可以得到Schwarzchild解,因此它与目前关于广义相对论的实验验证是一致的,但在“微观”极限下,方程预示着一种新的短程作用,讨论了自旋1/2的粒子作为检测粒子在这种球对称真空静场中的运动,指出运动方程只与仿射联络的黎曼部分有关,并和广义相对论的相应方程具有同样的形式。 关键词：     

Nonlinear,nondispersive wave equations: Lagrangian and Hamiltonian functions in the hodograph transformation

The hodograph transformation is generally used in order to associate a system of linear partial differential equations to a system of nonlinear (quasilinear) differential equations by interchanging dependent and independent variables. Here we consider the case when the nonlinear differential system can be derived from a Lagrangian density and revisit the hodograph transformation within the formalism of the Lagrangian-Hamiltonian continuous dynamical systems.Restricting to the case of nondissipative, nondispersive one-dimensional waves, we show that the hodograph transformation leads to a linear partial differential equation for an unknown function that plays the role of the Lagrangian in the hodograph variables. We then define the corresponding hodograph Hamiltonian and show that it turns out to coincide with the wave amplitude. i.e., with the unknown function of the independent variables to be solved for in the initial nonlinear wave equation.     

Material point method applied to multiphase flows

The particle-in-cell method (PIC), especially the latest version of it, the material point method (MPM), has shown significant advantage over the pure Lagrangian method or the pure Eulerian method in numerical simulations of problems involving large deformations. It avoids the mesh distortion and tangling issues associated with Lagrangian methods and the advection errors associated with Eulerian methods. Its application to multiphase flows or multi-material deformations, however, encounters a numerical difficulty of satisfying continuity requirement due to the inconsistence of the interpolation schemes used for different phases. It is shown in Section 3 that current methods of enforcing this requirement either leads to erroneous results or can cause significant accumulation of errors. In the present paper, a different numerical method is introduced to ensure that the continuity requirement is satisfied with an error consistent with the discretization error and will not grow beyond that during the time advancement in the calculation. This method is independent of physical models. Its numerical implementation is quite similar to the common method used in Eulerian calculations of multiphase flows. Examples calculated using this method are presented.     

General relativity and quantum mechanics in five dimensions

In 5D, I take the metric in canonical form and define causality by null-paths. Then spacetime is modulated by a factor equivalent to the wave function, and the 5D geodesic equation gives the 4D Klein-Gordon equation. These results effectively show how general relativity and quantum mechanics may be unified in 5D.     

Derivation of Dirac's Equation from the Evans Wave Equation

The Evans wave equation [1] of general relativity is expressed in spinor form, thus producing the Dirac equation in general relativity. The Dirac equation in special relativity is recovered in the limit of Euclidean or flat spacetime. By deriving the Dirac equation from the Evans equation it is demonstrated that the former originates in a novel metric compatibility condition, a geometrical constraint on the metric vector qused to define the Einstein metric tensor. Contrary to some claims by Ryder, it is shown that the Dirac equation cannot be deduced unequivocally from a Lorentz boost in special relativity. It is shown that the usually accepted method in Clifford algebra and special relativity of equating the outer product of two Pauli spinors to a three-vector in the Pauli basis leads to the paradoxical result X = Y = Z = 0. The method devised in this paper for deriving the Dirac equation from the Evans equation does not use this paradoxical result.     

Six-dimensional Dirac equation

The Dirac equation in six-dimensional relativity (three space and three time) is considered and shown to correspond to particles which have spatial spin-1/2 and temporal spin-1/2. Explicit forms of the spinor transformation are found. Plane wave solutions are obtained and their properties are given in terms of spatial and temporal spins and helicities. An expression is found for the charge conjugation operator.     

Two-fluid equations for phonons without momentum

In the customary microscopic derivation of the two-fluid equations of liquid He II explicit use is made of the assumption that an elementary excitation of wave vector k carries a momentum ?k. In this paper it is shown that phonons in a liquid can be defined as carrying momentum ?k (Eulerian phonons) or zero momentum (Lagrangian phonons). A careful analysis—in particular of the concept “velocity of the bearer fluid”—shows that the two-fluid equations turn out to be the same in either case.     

Solution of the scalar wave equation over very long distances using nonlinear solitary waves: Relation to finite difference methods

The linear wave equation represents the basis of many linear electromagnetic and acoustic propagation problems. Features that a computational model must have, to capture large scale realistic effects (for over the horizon or “OTH” radar communication, for example), include propagation of short waves with scattering and partial absorption by complex topography. For these reasons, it is not feasible to use Green’s Function or any simple integral method, which neglects these intermediate effects and requires a known propagation function between source and observer. In this paper, we describe a new method for propagating such short waves over long distances, including intersecting scattered waves. The new method appears to be much simpler than conventional high frequency schemes: Lagrangian “particle” based approaches, such as “ray tracing” become very complex in 3-D, especially for waves that may be expanding, or even intersecting. The other high frequency scheme in common use, the Eikonal, also has difficulty with intersecting waves.Our approach, based on nonlinear solitary waves concentrated about centroid surfaces of physical wave features, is related to that of Whitham [1], which involves solving wave fronts propagating on characteristics. Then, the evolving electromagnetic (or acoustic) field can be approximated as a collection of propagating co-dimension one surfaces (for example, 2-D surfaces in three dimensions). This approach involves solving propagation equations discretely on an Eulerian grid to approximate the linear wave equation. However, to propagate short waves over long distances, conventional Eulerian numerical methods, which attempt to resolve the structure of each wave, require far too many grid cells and are not feasible on current or foreseeable computers. Instead, we employ an “extended” wave equation that captures the important features of the propagating waves. This method is first formulated at the partial differential equation (PDE) level, as a wave equation with an added “confining” term that involves both a positive and a negative dissipation. Once we have the stable PDE, the discrete formulation is simply a multidimensional PDE with (stable) perturbations caused by the discretization. The resulting discrete solution can then be low order and very simple and yet remain stable over arbitrarily long times. When discretized and solved on an Eulerian grid, this new method allows far coarser grids than required by conventional resolution considerations, while still accounting for the effects of varying atmospheric and topographic features. An important point is that the new method is in the same form as conventional discrete wave equation methods. However, the conventional solution eventually decays, and only the “intermediate asymptotic” solution can be used. Simply by adding an extra term, we show that a nontrivial true asymptotic solution can be obtained. A similar solitary wave based approach has been used successfully in a different problem (involving “Vorticity Confinement”), for a number of years.     

Comparisons of Lagrangian and Eulerian PDF methods in simulations of non-premixed turbulent jet flames with moderate-to-strong turbulence-chemistry interactions

Transported probability density function (PDF) methods have been applied widely and effectively for modelling turbulent reacting flows. In most applications of PDF methods to date, Lagrangian particle Monte Carlo algorithms have been used to solve a modelled PDF transport equation. However, Lagrangian particle PDF methods are computationally intensive and are not readily integrated into conventional Eulerian computational fluid dynamics (CFD) codes. Eulerian field PDF methods have been proposed as an alternative. Here a systematic comparison is performed among three methods for solving the same underlying modelled composition PDF transport equation: a consistent hybrid Lagrangian particle/Eulerian mesh (LPEM) method, a stochastic Eulerian field (SEF) method and a deterministic Eulerian field method with a direct-quadrature-method-of-moments closure (a multi-environment PDF-MEPDF method). The comparisons have been made in simulations of a series of three non-premixed, piloted methane–air turbulent jet flames that exhibit progressively increasing levels of local extinction and turbulence-chemistry interactions: Sandia/TUD flames D, E and F. The three PDF methods have been implemented using the same underlying CFD solver, and results obtained using the three methods have been compared using (to the extent possible) equivalent physical models and numerical parameters. Reasonably converged mean and rms scalar profiles are obtained using 40 particles per cell for the LPEM method or 40 Eulerian fields for the SEF method. Results from these stochastic methods are compared with results obtained using two- and three-environment MEPDF methods. The relative advantages and disadvantages of each method in terms of accuracy and computational requirements are explored and identified. In general, the results obtained from the two stochastic methods (LPEM and SEF) are very similar, and are in closer agreement with experimental measurements than those obtained using the MEPDF method, while MEPDF is the most computationally efficient of the three methods. These and other findings are discussed in detail.