Nonlinear stability of overcompresive shock waves in a rotationally invariant system of viscous conservation laws

This paper proves that certain non-classical shock waves in a rotationally invariant system of viscous conservation laws posses nonlinear large-time stability against sufficiently small perturbations. The result applies to small intermediate magnetohydrodynamic shocks in the presence of dissipation.Research supported by Deutsche ForschungsgemeinschaftResearch supported in part by NSF Grant DMS 90-0226 and Army Grant DAAL 03-91-G-0017     

Time-Asymptotic Behavior of Wave Propagation Around a Viscous Shock Profile

We study the nonlinear stability of shock waves for viscous conservation laws. Our approach is based on a new construction of a fundamental solution for a linearized system around a shock profile. We obtain, for the first time, the pointwise estimates of nonlinear wave interactions across a shock wave. Our results apply to all ranges of weak shock waves and small perturbations. In particular, our results reduce to the time-asymptotic behavior of constant state perturbation, uniformly as the strength of the shock wave tends to zero. The research of the first author was partially supported by NSC Grant 96-2628-M-001-011 and NSF Grant DMS-0709248. The research of the second author was partially supported by NSF Grant DMS-0207154 and UAB Advance Program, sponsored by NSF.     

Nonlinear stability of rarefaction waves for compressible Navier Stokes equations

It is shown that expansion waves for the compressible Navier-Stokes equations are nonlinearly stable. The expansion waves are constructed for the compressible Euler equations based on the inviscid Burgers equation. Our result shows that Navier-Stokes equations and Euler equations are time-asymptotically equivalent on the level of expansion waves. The result is proved using the energy method, making essential use of the expansion of the underlining nonlinear waves and the specific form of the constitutive eqution for a polytropic gas.Supported in part by NSF Grant DMS-87-03971 and Army Grant DAAL03-87-K-0063Supported in part by Army Grant DAAL03-87-K-0063     

Stability of shock waves for the Broadwell equations

For the Broadwell model of the nonlinear Boltzmann equation, there are shock profile solutions, i.e. smooth traveling waves that connect two equilibrium states. For weak shock waves, we prove asymptotic (in time) stability with respect to small perturbations of the initial data. Following the work of Liu [7] on shock wave stability for viscous conservation laws, the method consists of analyzing the solution as the sum of a shock wave, a diffusive wave, a linear hyperbolic wave and an error term. The diffusive and linear hyperbolic waves are approximate solutions of the fluid dynamic equations corresponding to the Broadwell model. The error term is estimated using a variation of the energy estimates of Kawashima and Matsumura [6] and the characteristic energy method of Liu [7].Research supported by the Office of Naval Research through grant N00014-81-0002 and by the National Science Foundation through grant NSF-MCS-83-01260Research supported by the National Science Foundation through grant DMS-84-01355     

An Integral Spectral Representation of the Propagator for the Wave Equation in the Kerr Geometry

We consider the scalar wave equation in the Kerr geometry for Cauchy data which is smooth and compactly supported outside the event horizon. We derive an integral representation which expresses the solution as a superposition of solutions of the radial and angular ODEs which arise in the separation of variables. In particular, we prove completeness of the solutions of the separated ODEs. This integral representation is a suitable starting point for a detailed analysis of the long-time dynamics of scalar waves in the Kerr geometry. Research supported by NSERC grant # RGPIN 105490-2004. Research supported in part by the NSF, Grant No. DMS-010-3998. Research supported in part by the NSF, Grant No. 33-585-7510-2-30.     

Shock profile solutions of the Boltzmann equation

Shock waves in gas dynamics can be described by the Euler Navier-Stokes, or Boltzmann equations. We prove the existence of shock profile solutions of the Boltzmann equation for shocks which are weak. The shock is written as a truncated expansion in powers of the shock strength, the first two terms of which come exactly from the Taylor tanh (x) profile for the Navier-Stokes solution. The full solution is found by a projection method like the Lyapunov-Schmidt method as a bifurcation from the constant state in which the bifurcation parameter is the difference between the speed of soundc ₀ and the shock speeds.Research supported in part by the National Science Foundation, the Army Research Office, the Air Force of Scientific Research, the Office of Naval Research, and the Department of Energy     

A Rigorous Treatment of Energy Extraction from a Rotating Black Hole

The Cauchy problem is considered for the scalar wave equation in the Kerr geometry. We prove that by choosing a suitable wave packet as initial data, one can extract energy from the black hole, thereby putting supperradiance, the wave analogue of the Penrose process, into a rigorous mathematical framework. We quantify the maximal energy gain. We also compute the infinitesimal change of mass and angular momentum of the black hole, in agreement with Christodoulou’s result for the Penrose process. The main mathematical tool is our previously derived integral representation of the wave propagator. Research supported in part by the Deutsche Forschungsgemeinschaft. Research supported by NSERC grant # RGPIN 105490-2004. Research supported in part by the Humboldt Foundation and the National Science Foundation, Grant No. DMS-0603754. Research supported in part by the NSF, Grant No. 33-585-7510-2-30.     

A Hybrid Lattice Boltzmann Flux Solver for Simulation of Viscous Compressible Flows

In this paper, a hybrid lattice Boltzmann flux solver (LBFS) is proposed for simulation of viscous compressible flows. In the solver, the finite volume method is applied to solve the Navier-Stokes equations. Different from conventional Navier-Stokes solvers, in this work, the inviscid flux across the cell interface is evaluated by local reconstruction of solution using one-dimensional lattice Boltzmann model, while the viscous flux is still approximated by conventional smooth function approximation. The present work overcomes the two major drawbacks of existing LBFS [28–31], which is used for simulation of inviscid flows. The first one is its ability to simulate viscous flows by including evaluation of viscous flux. The second one is its ability to effectively capture both strong shock waves and thin boundary layers through introduction of a switch function for evaluation of inviscid flux, which takes a value close to zero in the boundary layer and one around the strong shock wave. Numerical experiments demonstrate that the present solver can accurately and effectively simulate hypersonic viscous flows.     

Ground state representation of the infinite one-dimensional Heisenberg ferromagnet

In its ground state representation, the infinite, spin 1/2 Heisenberg chain provides a model for spin wave scattering, which entails many features of the quantum mechanicalN-body problem. Here, we give a complete eigenfunction expansion for the Hamiltonian of the chain in this representation, forall numbers of spin waves. Our results resolve the questions of completeness and orthogonality of the eigenfunctions given by Bethe for finite chains, in the infinite volume limit.Research supported in part by NSF Grant No. MCS-76-05857Research supported in part by NSF Grant No. MCS-74-07313-A02     

激波冲击火焰的流动拓扑研究

为深入研究激波冲击火焰现象的内在机制,采用二维带化学反应的Navier-Stokes方程对现象进行数值研究,通过对速度梯度张量特征方程的分析证明Okubo-Weiss函数适用于可压缩流动,并重点分析火焰区的流动拓扑特性.结果表明,波后火焰区内Okubo-Weiss函数积分量基本守恒,但在火焰区内部和表面具有截然不同的流动状态,且火焰发展基本不受流场可压缩性的影响;波后火焰区的流动拓扑分类主要以焦点和鞍点为主,意味着流场中变形占主导.     

Decay of Solutions of the Wave Equation in the Kerr Geometry

We consider the Cauchy problem for the massless scalar wave equation in the Kerr geometry for smooth initial data compactly supported outside the event horizon. We prove that the solutions decay in time in L ^∞ _loc. The proof is based on a representation of the solution as an infinite sum over the angular momentum modes, each of which is an integral of the energy variable ω on the real line. This integral representation involves solutions of the radial and angular ODEs which arise in the separation of variables. Research supported in part by the Deutsche Forschungsgemeinschaft. Research supported by NSERC grant #RGPIN 105490-2004. Research supported in part by the NSF, Grant No. DMS-010-3998. Research supported in part by the NSF, Grant No. 33-585-7510-2-30. An erratum to this article is available at .     

Bleustein-Gulyaev waves in a functionally graded piezoelectric material layered structure

This work presents a theoretical study of the propagation behavior of Bleustein-Gulyaev waves in a layered structure consisting of a functionally graded piezoelectric material (FGPM) layer and a transversely isotropic piezoelectric substrate. The influence of the graded variation of FGPM coefficients on the dispersion relations of Bleustein-Gulyaev waves in the layered structure is investigated. It is demonstrated that, for a certain frequency range of Bleustein-Gulyaev waves, the mechanical perturbations of the particles are restricted in the FPGM layer and the phase velocity is independent of the electrical boundary conditions at the free surface. Results presented in this study can not only provide further insight on the electromechanical coupling behavior of surface waves in FGPM layered structures, but also lend a theoretical basis for the design of high-performance surface acoustic wave (SAW) devices. Supported by the National Natural Science Foundation of China (Grant No. 10632060), the National Basic Research Program of China (Grant No. 2006CB601202), the National 111 Project of China (Grant No. B06024), and the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20070698064)     

可压剪切层线性稳定性特征直接数值模拟

基于伪随机数生成技术促白噪声扰动,以高精度迎风／对称紧致混合差分算法求解二维／三维非定常可压Navier－Stokes方程,揭示了可压自由剪切层初始剪切过程中扰动的线性演化特征,以及该过程对扰动波数和方向的内在选择性,验证了所用算法的有效性,表明线性理论同数值模拟相结合是可压剪切层研究的合理途径之一。     

Guided wave propagation in multilayered piezoelectric structures

A general formulation of the method of the reverberation-ray matrix (MRRM) based on the state space formalism and plane wave expansion technique is presented for the analysis of guided waves in multilayered piezoelectric structures. Each layer of the structure is made of an arbitrarily anisotropic piezoelectric material. Since the state equation of each layer is derived from the three-dimensional theory of linear piezoelectricity, all wave modes are included in the formulation. Within the framework of the MRRM, the phase relation is properly established by excluding exponentially growing functions, while the scattering relation is also appropriately set up by avoiding matrix inversion operation. Consequently, the present MRRM is unconditionally numerically stable and free from computational limitations to the total number of layers, the thickness of individual layers, and the frequency range. Numerical examples are given to illustrate the good performance of the proposed formulation for the analysis of the dispersion characteristic of waves in layered piezoelectric structures. Supported by the National Natural Science Foundation of China (Grant Nos. 10725210 and 10832009), the Specialized Research Fund for the Doctoral Program of Higher Education (Grant No. 20060335107), the National Basic Research Program of China (Grant No. 2009CB623204), and the Scientific Research Foundation for Tsuiying Talents of Lanzhou University     

The validity of nonlinear geometric optics for weak solutions of conservation laws

The method of weakly nonlinear geometric optics is one of the main formal perturbation techniques used in analyzing nonlinear wave motion for hyperbolic systems. The tacit assumption in using such perturbation methods is that the corresponding solutions of the hyperbolic system remain smooth; since shock waves typically form in such solutions, these assumptions are rarely satisfied in practice. Nevertheless, in a variety of applied contexts, these methods give qualitatively reliable answers for discontinuous weak solutions. Here we give a rigorous proof for the validity of nonlinear geometric optics for general weak solutions of systems of hyperbolic conservation laws in a single space variable. The methods of proof do not mimic the formal construction of weakly nonlinear asymptotics but instead rely on structural symmetries of the approximating equations, stability estimates for intermediate asymptotic times, and the rapid decay in variation of weak solutions for large asymptotic times.Partially supported by NSF Grant No. DMS-8301135Partially supported by NSF Grant No. MCS-81-02360 and ARO Grant No. 483964-25530     

A Switch Function-Based Gas-Kinetic Scheme for Simulation of Inviscid and Viscous Compressible Flows

In this paper, a switch function-based gas-kinetic scheme (SF-GKS) is presented for the simulation of inviscid and viscous compressible flows. With the finite volume discretization, Euler and Navier-Stokes equations are solved and the SF-GKS is applied to evaluate the inviscid flux at cell interface. The viscous flux is obtained by the conventional smooth function approximation. Unlike the traditional gas-kinetic scheme in the calculation of inviscid flux such as Kinetic Flux Vector Splitting (KFVS), the numerical dissipation is controlled with a switch function in the present scheme. That is, the numerical dissipation is only introduced in the region around strong shock waves. As a consequence, the present SF-GKS can well capture strong shock waves and thin boundary layers simultaneously. The present SF-GKS is firstly validated by its application to the inviscid flow problems, including 1-D Euler shock tube, regular shock reflection and double Mach reflection. Then, SF-GKS is extended to solve viscous transonic and hypersonic flow problems. Good agreement between the present results and those in the literature verifies the accuracy and robustness of SF-GKS.     

On representation spaces in geometric quantization

The structure of the space of wave functions in the representation given by a complete strongly admissible polarization of the phase space is investigated. The conditions that the wave functions should be covariant constant along the real part of the polarization define the Bohr-Sommerfeld set of the representation containing the supports of all wave functions. There is a natural scalar product for the wave functions defined on the Bohr-Sommerfeld set. It is shown, for a real polarization, that the resulting Hilbert space of wave functions is not trivial if and only if the Bohr-Sommerfeld set is not empty.Partially supported by the National Research Council, Grant No. A8091.     

Preferred shapes and modes of internal motion in a four-boson system

A four-body quantum-mechanical system is studied and through suitable manipulations of its wave functions, dominant shapes and modes of internal motion are ascribed to a number of its low-lying 0⁺ states.Research supported in part by NSF Grant No. PHY-8712229 and Grant No. PHY-8945627     

Preferred shapes and modes of internal motion in a five-boson system

A five-body system is studied using harmonic-oscillator basis functions and through suitably-defined shape-density functions which are based on its wave functions, we are able to ascribe dominant shapes and modes of internal motion to a number of its low-lying 0⁺ states.Research supported in part by NSF Grant No. PHY-8712229 and Grant No. PHY-8945627     

Surface instability of the liquid turbulent flows in the open inclined channels

The flow of a viscous liquid layer in an open inclined channel under the turbulent mode is considered in this paper. To describe turbulent viscosity, the Van Driest model is used. The spectrum of characteristic values of the problem on linear stability of a plane-parallel flow is studied numerically. Parameters of the maximal growth waves are found out, the surface tension effect is studied, and theoretical results are compared with experimental data. The work was financially supported by the Russian Foundation for Basic Research (Grant No. 05-08-33585a).