爱因斯坦引力场Arnowitt-Deser-Misner约束方程的两种推导方法

给出两种不同方法,分别导出爱因斯坦引力理论中著名的Arnowitt-Deser-Misner(ADM)约束方程.其一是在具有洛伦兹号差的时空中,构造一个单参数引力场作用量,由此导出单参数ADM约束方程.该参数取某特定值时对应的就是熟知的ADM约束方程.其二是将二重复函数理论运用于爱因斯坦引力场的哈密顿形式表述中,得到引力场ADM约束的二重化形式,从而也能将通常的ADM约束作为其特殊情况包含其中.此外,这两种方法还能统一地表述具有不同时空号差(洛伦兹号差和欧几里得号差)的洛伦兹引力理论和欧几里得引力理论 关键词： Arnowitt-Deser-Misner约束方程 哈密顿表述 时空号差 引力场作用量     

Maslov’s asymptotic methods in problems of the theory of optical lattices

Application of Maslov’s asymptotic methods to equations that arise in the theory of optical lattices was considered. The occurrence of a small parameter in the Schrödinger equations with potentials of three-dimensional and controlled optical lattice was investigated and the conditions of application of Maslov’s asymptotic methods for solving these equations was determined. Consideration of different conditions imposed on the parameters in these potentials led to the use of two different methods for solving the arising equations: Maslov’s method of a complex germ and Maslov’s operator-valued method of a complex germ. Expressions that can be used to calculate the desired characteristics of the atomic systems under consideration in optical lattices were derived.     

Free vibration of three-dimensional anisotropic layered composite nanoplates based on modified couple-stress theory

Based on the modified couple-stress theory, three-dimensional analytical solutions of free vibration of a simply supported, multilayered and anisotropic composite nanoplate are derived by solving an eigenvalue system and using the propagator matrix method. By expanding the solutions of the extended displacements in terms of two-dimensional Fourier series, the final governing equations of motion with modified couple-stress effect are reduced to an eigenvalue system of ordinary differential equations. Analytical expressions for the natural frequencies and mode shapes of multilayered anisotropic composite plates with modified couple-stress effect are then derived via the propagator matrix method. Numerical examples are carried out for homogeneous thick-plates and sandwich composite plates to show the effect of the non-local parameter in different layers and stacking sequence on the mode shapes. The present solutions can serve as benchmarks to various thick-plate theories and numerical methods, and could be further useful for designing layered composite structures involving small scale.     

Adjoint-based sensitivity analysis of flames

Simulation of chemically reacting flows using detailed chemistry introduces a large number of chemistry model parameters. While not all significantly affect the target outcomes of a simulation, the parameters that do are not always known a priori. In order to improve simulations for specified target outcomes, termed quantities of interest (QoIs), the sensitivity of these QoIs to the model parameters are needed. However, evaluating the sensitivities is computationally expensive, especially for complex fuels that may involve many parameters. For these simulations, the forward sensitivity method requires the solution of an additional number of governing equations proportional to the number of parameters. Here, an adjoint sensitivity approach is formulated where the computational cost scales as the number of QoIs and not the number of parameters. Specifically, adjoint equations are derived for laminar, incompressible, variable density reacting flow and applied to hydrogen flame simulations. From the solution of the corresponding adjoint equations, sensitivity of the QoIs to chemistry model parameters is calculated. The one-dimensional simulation results show that the adjoint sensitivity results closely match those of forward sensitivity methods, thus providing validation of the adjoint method. The two-dimensional simulation results indicate the most sensitive parameters for two QoIs, flame tip temperature and NOx emission. For these tests, the adjoint method reduces computational expense compared to forward sensitivity methods by a factor proportional to the number of QoIs over the number of parameters, here 2/172. Such savings can be more drastic for cases that involve complex fuels, such as combustion of jet fuel, requiring thousands of chemistry model parameters. Further, this sensitivity information can be used in development of experiments by pointing out which are the critical chemistry model parameters.     

Invariant manifold methods for metabolic model reduction

After the decay of transients, the behavior of a set of differential equations modeling a chemical or biochemical system generally rests on a low-dimensional surface which is an invariant manifold of the flow. If an equation for such a manifold can be obtained, the model has effectively been reduced to a smaller system of differential equations. Using perturbation methods, we show that the distinction between rapidly decaying and long-lived (slow) modes has a rigorous basis. We show how equations for attracting invariant (slow) manifolds can be constructed by a geometric approach based on functional equations derived directly from the differential equations. We apply these methods to two simple metabolic models. (c) 2001 American Institute of Physics.     

A computation of the reflection and transmission coefficients for dielectric coated waveguides

A previously developed variational technique for finding the approximate solution to the electromagnetic field inside waveguides of varying shapes and containing non-uniform dielectric and magnetic materials is applied to the specific case of an axisymmetric waveguide containing three layers of different dielectric materials. The magnetic permeability is taken equal to unity. The method enables the problem to be reduced to a two point boundary value problem for a pair of second order, linear, ordinary differential equations. The values of the reflection and transmission coefficients obtained by this method are in good agreement with those derived from solving the partial differential equations for the field.     

Connection between the immersing and transfer-matrix methods in one-dimensional scattering problems

We show that application of the immersing and transfer-matrix methods to one-dimensional problems of particles scattering leads to the system of two linear equations for the functions F and Φ expressed by means of the transmission and reflection amplitudes. The expressions of these functions are derived. The offered method is illustrated by the finding of transmission and reflection coefficients for the potential barrier with a constant height. The developed method can be applied in solving the quasi-one-dimensional and two-dimensional problems of scattering.     

Asymptotic solutions for Mathieu instability under random parametric excitation and nonlinear damping

A theoretical analysis is presented of the response of a lightly and nonlinearly damped mass-spring system in which the spring constant contains a small randomly fluctuating component. Damping is represented by a combination of linear and nonlinear power-law damping. System response to some initial disturbance at time zero is described by a sinusoidal wave whose amplitude and phase vary slowly and randomly with time. Leading order formulations for the equations of amplitude and phase are obtained through the application of methods of stochastic averaging of Stratonovich. The equations of amplitude and phase are given in two versions: Fokker-Planck equations for transient probability and Langevin equations for response in the time-domain. Solutions in closed-form of these equations are derived by methods of mathematical and theoretical physics involving higher transcendental functions. They are used to study the behavior of system response for ever increasing time applying asymptotic methods of analysis such as the method of steepest descent or saddle-point method. It is found that system behavior depends on the power density of the parametric excitation at twice the natural frequency and on the magnitude and form of the damping. Depending on these parameters different types of system behavior are found to be possible: response which decays exponentially to zero, response which leads to a stationary state of random behavior, and response which can either grow unboundedly or which approaches zero in a finite time.     

Granular matter and the time-dependent viscous eikonal equation

Deposition of granular matter under gravity can be described by the well-known two-layer model for a standing and a rolling layer. Matter from sources enters the rolling layer which flows along the gradient of the standing layer and finally enters the standing layer via interaction of the two layers. From this system of two coupled hyperbolic partial differential equations a time-dependent viscous eikonal equation is derived as a limiting case for weak sources, a thin rolling layer and fast convection of the rolling layer. This equation, supplied with boundary conditions, describes the deposition of dry sand from evenly distributed sources onto a flat table with a vertical rim of variable height. The stationary problem can also be seen as an application of the method of vanishing viscosity to the eikonal equation. For certain types of interaction between the two layers the resulting eikonal equation can be transformed into a linear equation. This transformation yields additional insight into the problem.     

Calculation of self-sputtering spectra of thin films by the discrete flow method

Self-sputtering of plane-parallel layers is described theoretically on the basis of the invariant immersion method combined with the discrete flow technique. Differential equations are derived for forward and backward self-sputtering functions, and exact solutions to these equations are obtained. Energy distributions of self-sputtered ions are calculated; total self-sputtering coefficients are compared to the results of computer experiment. The dependences of the total self-sputtering coefficients on the film thickness and the initial energy of ions are analyzed.     

Determination of the Electromagnetic Lagrangian from a System of Poisson Brackets

The Lagrangian and Hamiltonian formulations of electromagnetism are reviewed and the Maxwell equations are obtained from the Hamiltonian for a system of many electric charges. It is shown that three of the equations which were obtained from the Hamiltonian, namely the Lorentz force law and two Maxwell equations, can be obtained as well from a set of postulated Poisson brackets. It is shown how the results derived from these brackets can be used to reconstruct the original Lagrangian for the theory aided by some reasoning based on physical concepts.     

On one-port characterization of noise sources in ducts by using external loads

Equivalent acoustic source characterization of duct-borne fluid machinery noise is often undertaken by interpolating the results of two-microphone pressure measurements with different external acoustic loads over a linear one-port source model. If the source is time-invariant, the one-port source characteristics can be determined by using only two external loads. This is well known as the two-load method. An extension of the two-load method for time-variant sources is also available and known as the multiple-load method. In these methods the source is treated as a ‘black-box’. This paper addresses the problem of one-port source characterization when the linear operations inherent in the ‘black-box’ are known explicitly. The equations governing the explicit one-port source models are derived and the source characteristics are shown to be measurable using only few acoustic loads. It is not the purpose of this paper to discuss the application of these models to any specific fluid machinery; however, of particular interest are the explicit source models that require only two loads. Numerical results are presented to show some features of such time-invariant and time-variant explicit one-port source models.     

Random vibration analysis of stochastic time-varying systems

The analysis of the free and forced vibration of a randomly time-varying system is the subject matter of this paper. This is a complicated problem which has received relatively little discussion in the literature. Herein two methods are presented, apart from the digital simulation technique, of finding the response moments. The first one is a series technique which can be considered as a generalization of the well known Galerkin method. The second method belongs to the class of closure techniques. Upon presuming some of the joint distributions to be Gaussian, equations are derived for the first two response moments. It is shown further that the non-Gaussian output density can be approximately predicted by a simple transformation. Detailed numerical results are obtained and compared with computer simulated response statistics. It is demonstrated that the methods developed here are highly efficient. In particular it is found that the Gaussian closure approximation has a wide range of application.     

A critical evaluation of closure methods via two simple dynamic systems

Two simple dynamic systems with cubic nonlinearity and additive Gaussian white noise are used to assess the performance and the usefulness of closure methods in nonlinear random vibration. One of the systems has a single potential well while the other has two potential wells. It is shown that the performance of closure methods is determined by the structure of the moment equations rather than the way in which these equations are closed. For the system with one potential well, any closure method is satisfactory. For the system with two potential wells, closure methods can be inaccurate irrespective of the closure level. It is also shown that moment equations can be augmented with moment inequalities to solve approximately the infinite hierarchy of moment equations.     

A higher order dynamic theory for viscoelastic plates and layered composites

By using a new technique proposed by the first author [1] approximate theories are developed for the dynamic response of viscoelastic plates and layered composites. The originality of the new technique lies in the fact that it permits the approximate theory to satisfy correctly the lateral boundary conditions of a plate, or the interface (continuity) conditions of a layered composite. This, in turn, enables the approximate theory to describe accurately the geometric dispersion of waves propagating in a plate or layered composite. The approximate equations of a single viscoelastic plate are first derived by making use of the new technique. To develop the approximate theory for viscoelastic layered composites made of two alternating layers it is noted that the approximate equations of a single plate already established also hold in each layer of the composite. The theory is completed by adding the continuity conditions to these equations and using a smoothing operation. The equations thus obtained constitute a continuum (homogeneous) model (CM) which simplifies the determination of the dynamic response of viscoelastic layered composites when the number of layers in the composite is large. The proposed approximate theories are open to improvement in the sense that their region of validity in the wave number-frequency plane can be enlarged as much as one wishes by increasing the orders of the theories and continuity conditions.     

The refined theory of deep rectangular beams based on general solutions of elasticity

The problem of deducing one-dimensional theory from two-dimensional theory for a homogeneous isotropic beam is investigated. Based on elasticity theory, the refined theory of rectangular beams is derived by using Papkovich-Neuber solution and Lur’e method without ad hoc assumptions. It is shown that the displacements and stresses of the beam can be represented by the angle of rotation and the deflection of the neutral surface. Based on the refined beam theory, the exact equations for the beam without transverse surface loadings are derived and consist of two governing differential equations: the fourth-order equation and the transcendental equation. The approximate equations for the beam under transverse loadings are derived directly from the refined beam theory and are almost the same as the governing equations of Timoshenko beam theory. In two examples, it is shown that the new theory provides better results than Levinson’s beam theory when compared with those obtained from the linear theory of elasticity.     

Refraction of active waves in reaction—diffusion media

The refraction of active waves is analyzed for a stable—metastable reaction—diffusion system consisting of two regions with different diffusion coefficients. The equations governing the evolution of wavefronts are derived by means of an asymptotic perturbation method for boundary layers. These equations describe non-stationary refraction near the steady state regime. It is shown that the dynamics of wavefronts separates into that in the region near the boundary and that far from the boundary.     

On continuation of inviscid vortex patches

This investigation concerns solutions of the steady-state Euler equations in two dimensions featuring finite area regions with constant vorticity embedded in a potential flow. Using elementary methods of the functional analysis we derive precise conditions under which such solutions can be uniquely continued with respect to their parameters, valid also in the presence of the Kutta condition concerning a fixed separation point. Our approach is based on the Implicit Function Theorem and perturbation equations derived using shape-differentiation methods. These theoretical results are illustrated with careful numerical computations carried out using the Steklov–Poincaré method which show the existence of a global manifold of solutions connecting the point vortex and the Prandtl–Batchelor solution, each of which satisfies the Kutta condition.