两微极磁热弹性固体分界面上波的反射与折射

给出了磁场、热场和弹性场多场耦合作用下微极广义热弹性固体的一般控制方程.该方 程既包含了磁场、热场和弹性场的耦合作用,又在其广义热传导方程中涵盖了耦合热弹理论 (C-D)及其5类推广(L-S理论,G-L理论,G-N(II,III)理论和C-T理论).运用该微极广义磁热 弹性控制方程,研究了在定常磁场作用下, 具有均匀初始温度的两理想接触微极弹性介质平面分界面上磁热弹性波的反射和折射现象.给出了分别在缺少磁场、热场作用或不同广义热传 导理论下反射或折射热波、纵向位移波、耦合横向和微旋转波与入射纵向位移波的振幅比随 入射角变化的关系曲线.对缺少磁、热和微极性以及热松弛时间时对应的反射、折射系数进 行了对比.结果表明磁、热和微极性以及热松弛时间对振幅比均有不同程度的影 响,与磁、热和微极性一样,热松弛时间对不同类型波的影响能力差别明显,但对同 一类型的反射波和折射波的影响相似.     

Reflection and refraction of attenuated waves at boundary of elastic solid and porous solid saturated with two immiscible viscous fluids

The propagation of elastic waves is studied in a porous solid saturated with two immiscible viscous fluids.The propagation of three longitudinal waves is represented through three scalar potential functions.The lone transverse wave is presented by a vector potential function.The displacements of particles in different phases of the aggregate are defined in terms of these potential functions.It is shown that there exist three longitudinal waves and one transverse wave.The phenomena of reflection and refraction due to longitudinal and transverse waves at a plane interface between an elastic solid half-space and a porous solid half-space saturated with two immiscible viscous fluids are investigated.For the presence of viscosity in pore-fluids,the waves refracted to the porous medium attenuate in the direction normal to the interface.The ratios of the amplitudes of the reflected and refracted waves to that of the incident wave are calculated as a nonsingular system of linear algebraic equations.These amplitude ratios are used to further calculate the shares of different scattered waves in the energy of the incident wave.The modulus of the amplitude and the energy ratios with the angle of incidence are computed for a particular numerical model.The conservation of the energy across the interface is verified.The effects of variations in non-wet saturation of pores and frequencies on the energy partition are depicted graphically and discussed.     

Wave reflection and refraction in triclinic crystalline media

Summary In this paper, the reflection and refraction of a plane wave at an interface between two half-spaces composed of triclinic crystalline material is considered. It is shown that due to incidence of plane wave three types of waves, namely quasi-P (qP), quasi-SV (qSV) and quasi-SH (qSH), will be generated governed by the propagation condition involving the acoustic tensor. A simple procedure has been presented for the calculation of all the three phase velocities of the quasi waves. It has been established that the direction of particle motion is neither parallel nor perpendicular to the direction of propagation. Relations are established between directions of motion and propagation, respectively. The expressions for reflection and refraction coefficients of qP, qSV and qSH waves are obtained. Numerical results of reflection and refraction coefficients are presented for different types of anisotropic media and for different types of incident waves. Graphical representations have been made for incident qP waves, and for incident qSV and qSH waves numerical data are presented in tables.The work was completed while the author was visiting the University of Kaiserslautern, Department of Geomathematics as Visiting Professor. The Author is grateful to Professor Dr. W. Freeden for providing DAAD fellowship and all the facilities for conducting research, as well as to Dr. V.Michel for various discussions about the research work and also for all kinds of help during his stay at Kaiserslautern, Germany. This award is very gratefully acknowledged.     

A compact numerical algorithm for solving the time‐dependent mild slope equation

The mild slope equation has been widely used to describe combined wave refraction and diffraction. In this study, a new numerical algorithm is developed to solve the time‐dependent mild slope equation in a second‐order hyperbolic form. The numerical algorithm is based on a compact and explicit finite difference method that is second‐order accurate in both time and space. The algorithm has the similar structure to the leap‐frog method but is constructed on three time levels for the second‐order time derivative term. The numerical model has the capability of simulating transient wave motion by correctly predicting the speed of wave energy propagation, which is important for the real‐time forecast of the arrival time of storm waves generated in the far field. The model is validated against analytical solution for wave shoaling and experimental data for combined wave refraction and diffraction over a submerged elliptic shoal on a slope (Coastal Eng. 1982; 6 :255). Lastly, the realistic scale Homma's island (Geophys. Mag. 1950; 21 :199) is studied with the use of various wave periods of T = 720s, T = 120 s, and T = 24 s. These wave periods correspond to long, intermediate, and short waves for the given topography, respectively. Comparisons are made between numerical results and existing analytical solutions in terms of the wave amplification around the island, which serves as the indicator for the potential wave runup. Excellent agreements are obtained. The model runs on a PC (Pentium IV 1.8GHz) and the computer capacity allows the computation of a mesh system up to 3000 × 3000, which is equivalent to about 150 × 150 waves or a large area of 540km × 540km for a wave train with the period of T = 60 s. Copyright 2004 John Wiley & Sons, Ltd.     

固-固有限周期声子晶体中SH波全反射隧穿的机理

利用边界条件推导出SH波在一维固-固有限周期声子晶体中的转移矩阵,并用倏逝波的概念和多波束干涉理论,从理论上解释了一维固-固有限周期声子晶体全反射隧穿效应产生的物理机理,得出了全反射隧穿峰带波长满足的解析公式。利用解析公式和转移矩阵研究了一维固-固有限周期声子晶体中SH波的全反射隧穿效应,结果发现两种方法得出的结论是吻合的。     

两层流体中水波在垂直薄板上的反射与透射

研究在两层流体中表面波模态和内波模态的波浪与半潜式刚性垂直薄板 相互作用的问题. 基于特征函数展开理论,建立了两种模态入射波作用下,半潜式刚性垂直 薄板的反射与透射能量的计算方法,证明了对每一种模态的入射波,另一种模态波浪的反射 与透射能量是相等的. 对水面漂浮和座底半潜式薄板的反射与透射能量,以及作用在其上的 水平波浪力进行了数值计算分析,表明在某个频率范围内,流体的分层效应对这些水动力 量的影响是不可忽视的. 特别地,当薄板的一端位于两层流体的内界面上时,两种模态波浪 的能量转化是最大的.     

A numerical study of oblique shock wave reflections over wedges in an ideal quantum gas

The various oblique shock wave reflection patterns generated by a moving incident shock on a planar wedge using an ideal quantum gas model are numerically studied using a novel high resolution quantum kinetic flux splitting scheme. With different incident shock Mach numbers and wedge angles as flow parameters, four different types of reflection patterns, namely, the regular reflection, simple Mach reflection, complex Mach reflection and the double Mach reflection as in the classical gas can be classified and observed. Both Bose–Einstein and Fermi–Dirac gases are considered.      

Spurious currents in finite element based level set methods for two‐phase flow

A study of spurious currents in continuous finite element based simulations of the incompressible Navier–Stokes equations for two‐phase flows is presented on the basis of computations on a circular drop in equilibrium. The conservative and the standard level set methods are used. It is shown that a sharp surface tension force, expressed as a line integral along the interface, can give rise to large spurious currents and oscillations in the pressure that do not decrease with mesh refinement. If instead a regularized surface tension representation is used, exact force balance at the interface is possible, both for a fully coupled discretization approach and for a fractional step projection method. However, the numerical curvature calculation introduces errors that cause spurious currents. Different ways to extend the curvature from the interface to the whole domain are discussed and investigated. The impact of using different finite element spaces and stabilization methods is also considered. Copyright © 2011 John Wiley & Sons, Ltd.     

分层流体中内孤立波在台阶上的反射和透射

基于匹配渐近展开和格林函数的方法,研究了两层流体系统中内孤立波在台阶地形上透射、 反射及其分裂的演化特征. 通过保角变换和求解奇异Fredholm积分方程,获得了反映地形 效应对Boussinesq方程影响的约化边界条件,藉此建立了KdV演化方程的``初值'问题, 根据散射反演理论获得了反射波和透射波的解析表达式. 分析结果表明：上下流体层的厚度 比、密度比以及台阶高度对于反射和透射波振幅及其分裂具有显著的影响. 尤其当上层流体 厚度小于下层厚度时,由于存在临界点,在其附近反射波的幅值随台阶高度的演化由单调增 变为单调减,透射波的幅值由单调减变为单调增;上台阶的反射波与入射波反相,其最大幅 值可达到入射波的数倍;此外,下台阶反射波也可发展为单支孤立波,它区别于单层流体中 反射波仅为衰减的振荡波列.     

Influence of a rigid skeleton pore structure on wave propagation in a fluid-filling porous medium

This work concerns an analysis of the influence of a rigid skeleton pore structure on wave propagation in a fluid-filling porous medium. The analysis is based on the continuum theory of a deformable porous medium in which the pore structure is described by two macroparameters. Considerations comprise two questions: the influence of the pore structure on wave-propagation velocity analysed for the quasilinear case and the role of structure in the reflection-refraction wave phenomenon in fluid at the contact surface of two porous media. It has been shown that the pore structure reduces the velocity of wave and together with the angle of incidence it defines the reflection-refraction wave phenomenon.     

A comparison of four recent numerical schemes giving high resolution of shock wave and concentrated vortex

ＡＣＯＭＰＡＲＩＳＯＮＯＦＦＯＵＲＲＥＣＥＮＴＮＵＭＥＲＩＣＡＬＳＣＨＥＭＥＳＧＩＶｉＮＧＨＩＧＨＲＥＳＯＬＵＴＩＯＮＯＦＳＨＯＣＫＷＡＶＥＡＮＤＣＯＮＣＥＮＴＲＡＴＥＤＶＯＲＴＥＸＨｕａｎｇＤｕｎ（黄敦）（Ｄｅｐａｒｔ．ｏｆＭａｔｈ．ＰｅｋｉｎｇＵ...     

带隔板装药爆轰波马赫反射理论研究和数值模拟

基于三波理论和Whitham方法对带隔板装药爆轰波相互作用后发生的正规反射和非正规反射进行了理论分析,给出了爆轰波发生马赫反射时临界入射角和马赫杆增长角等参数的变化规律,提出了马赫杆高度的计算模型。基于凝聚炸药爆轰Jones-Wilkins-Lee(JWL)模型和冲击起爆的Lee-Tarver模型,利用有限元计算软件对带隔板装药爆轰波的传播过程进行了数值模拟。结果表明,发生马赫反射后,随着爆轰波的传播,马赫杆的高度不断增加。数值模拟结果与理论计算结果吻合较好,说明本文中采用的理论模型和数值模拟方法能够较准确地描述带隔板装药爆轰波马赫反射的传播过程。     

Interaction of two Mach reflections over concave double wedges-analytical model

A simplified analytical model for solving the wave configuration resulting when two triple points join together over a concave double wedge is developed. The model equations are solved analytically. Predictions of the model are compared to experimental results and fairly good agreement is obtained. It is believed that the present model can be used to better understand complex wave interaction phenomena. Received 2 May 1995 / Accepted 18 December 1997     

Specific features of incident and reflected blast waves

On the basis of numerical modeling specific features of shock wave reflections were analyzed. It was found, that after diaphragm rupture self-modeling pressure and velocity distributions nearby the shock front establish. But in some special cases the temperature behind the shock front can rise. This peculiarity should be taken into account when performing experiments with high reactive gaseous mixtures. The temperature on the shock front and the velocity gradient behind it are uniform in the case of strong blast wave reflections. This effect is observed in the zone with an elevated temperature profile behind the incident blast wave. The reflected triangular waves conserve a quasi-self-modeling character by pressure. Typical experiments were carried out to verify the theoretical predictions. The effects of reflected wave acceleration in the case of triangular waves and the self-similar character of the pressure profiles were observed.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

An adaptive artificial viscosity for the displacement shallow water wave equation

The numerical oscillation problem is a difficulty for the simulation of rapidly varying shallow water surfaces which are often caused by the unsmooth uneven bottom,the moving wet-dry interface, and so on. In this paper, an adaptive artificial viscosity(AAV) is proposed and combined with the displacement shallow water wave equation(DSWWE) to establish an effective model which can accurately predict the evolution of multiple shocks effected by the uneven bottom and the wet-dry interface. The effec...     

地下结构多尺度动力分析方法

近年地下结构震害频发,针对目前地下结构抗震研究中实际存在的多尺度动力问题,提出了地下结构多尺度动力分析方法,该方法不需要任何附加的过滤和阻尼就能有效地消除高频波的虚假反射. 以桥域耦合理论为基础,引入拉格朗日乘子将不同尺度区域之间的约束关系,通过能量势函数隐含到动力方程中,推导出不同尺度域的动力控制方程;基于中心差分法,提出了用于地下结构多尺度分析的动力显式算法,以求解所建立的多尺度动力耦合体系. 以实际工程为应用实例,通过与传统的位移耦合方法的对比分析,说明了该多尺度方法用于地下结构动力分析的可行性以及对消除高频波虚假反射的有效性.      

近水面水下爆炸的数值研究

近水面水下爆炸后形成的冲击波很快与空气一水的自由界面相互作用，产生折射激波和反射稀疏波，这些波在形成初期是规则的，但很快转变为不规则波。伴随冲击波与自由界面的相互作用，水面要上升直至被炸开。这些复杂的物理现象对数值计算提出了很高的要求。本文采用近期发展起来的位标函数方法以及高精度的ＮＮＤ格式来数值模拟冲击波、自由界面的运动及其相互作用，获得了成功，近水面水下爆炸时的一系列复杂的物理现象均得到了合理的再现。与现有的跟踪自由界面运动的方法相比较，本文方法更简洁，更实用。     

Pressure wave propagation in liquid-filled tubes of viscoelastic material

The propagation of small-amplitude waves in a thick-walled long viscoelastic tube of variable cross-section, filled with a viscous incompressible fluid, is considered with account for wave reflection at the tube end in application to arterial pulse wave propagation. A solution is obtained in the form of expansions in a small parameter. The effect of the coefficient of wave reflection at the tube end and the wall material parameters on the fluid volume flow-rate and the tube wall displacement is investigated. It is shown that the volume flow-rate phase spectrum characteristics depend only slightly on the wall properties and can be used in clinical diagnostics for finding the reflection coefficient from pressure and flow-rate records.     

Detonation Initiation by Shock Reflection from an Orifice Plate

Results from an experimental investigation of the interaction of a “non-ideal” shock wave and a single obstacle are reported. The shock wave is produced ahead of an accelerated flame in a 14 cm inner-diameter tube partially filled with orifice plates. The shock wave interacts with a single larger blockage orifice plate placed 15–45 cm after the last orifice plate in the flame acceleration section of the tube. Experiments were performed with stoichiometric ethylene–oxygen mixtures with varying amounts of nitrogen dilution at atmospheric pressure and temperature. The critical nitrogen dilution was found for detonation initiation. It is shown that detonation initiation occurs if the chemical induction time based on the reflected shock state is shorter than the time required for an acoustic wave to traverse the orifice plate upstream surface, from the inner to the outer diameter. The similarity between the present results and those obtained from previous investigators looking at detonation initiation by ideal shock reflection produced in a shock tube indicates that the phenomenon is not sensitive to the detailed structure of the shock front but only on the average shock strength.This paper is based on work that was presented at the 20th International Colloquium on the Dynamics of Explosions and Reactive Systems, Montreal, Canada, July 31–August 5, 2005.     

Nonreflecting boundary conditions based on nonlinear multidimensional characteristics

Nonlinear characteristic boundary conditions based on nonlinear multidimensional characteristics are proposed for 2‐ and 3‐D compressible Navier–Stokes equations with/without scalar transport equations. This approach is consistent with the flow physics and transport properties. Based on the theory of characteristics, which is a rigorous mathematical technique, multidimensional flows can be decomposed into acoustic, entropy, and vorticity waves. Nonreflecting boundary conditions are derived by setting corresponding characteristic variables of incoming waves to zero and by partially damping the source terms of the incoming acoustic waves. In order to obtain the resulting optimal damping coefficient, analysis is performed for problems of pure acoustic plane wave propagation and arbitrary flows. The proposed boundary conditions are tested on two benchmark problems: cylindrical acoustic wave propagation and the wake flow behind a cylinder with strong periodic vortex convected out of the computational domain. This new approach substantially minimizes the spurious wave reflections of pressure, density, temperature, and velocity as well as vorticity from the artificial boundaries, where strong multidimensional flow effects exist. The numerical simulations yield accurate results, confirm the optimal damping coefficient obtained from analysis, and verify that the method substantially improves the 1‐D characteristics‐based nonreflecting boundary conditions for complex multidimensional flows. Copyright © 2009 John Wiley & Sons, Ltd.