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.     

On the performance of adaptive mesh refinement computation

A finite difference scheme for the unsteady Euler equations using an adaptive mesh refinement (AMR) algorithm was applied to the time-dependent flowfield of shock diffraction problems. The effectiveness of the AMR algorithm was evaluated against a uniform mesh algorithm. Computational results showed that to obtain solutions with equivalent resolution, the AMR algorithm requires much less processing time, when compared with a uniform mesh algorithm.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.     

从流场干涉图定量计算冲击波传播速度和压力

研究了含冲击波流场干涉图的图像处理方法。首先利用快速条纹细化算法细化干涉条纹,并利用条纹追踪算法提取冲击波波阵面,然后根据配准算法把序列冲击波波阵面配准在同一幅画面上。根据上述算法,在PC-VISION100图像处理机上开发了应用软件。实验结果表明,该方法不但可用于从流场干涉图中提取冲击波波阵面,更重要的是可用于从含冲击波流场干涉图中定量计算冲击波的传播速度和压力分布。     

计算动力学中的伪弧长方法研究

动力学问题通常采用微分方程来描绘,但由于工程实际问题的复杂性,微分方程模型常伴随着解的不连续性、刚性或激波间断奇异性特点,传统方法很难求解,奇异性问题是计算动力学难点,同时也是国内外学者研究的热点.伪弧长数值算法是针对计算动力学中的奇异性问题所提出的,其基本思想为通过在解曲线上引入伪弧长参数,并增加一个约束方程,在伪弧长参数作用下,使得原始离散单元发生扭曲形变,从而达到消除或减弱奇异性的目的.本文首先介绍伪弧长方法求解定常对流-扩散方程的奇异性问题,并提出针对双曲守恒定律的局部伪弧长算法,其思想在于首先通过间断解的梯度变换来确定强间断所处位置,进而通过局部网格点重构以及数值修正来达到强间断处奇异性消除与降低的目的.针对高维问题,提出全局伪弧长方法,通过对整个计算区域内的网格点进行重构,使得所有网格点向奇异间断点处移动,从而降低间断点的影响域,达到降低奇异性的目的.重点讨论了三维全局伪弧长算法问题的计算难点,即三维空间网格扭曲大变形导致的数值算法不收敛,并提出在算法设计过程中采用分块重构与整体计算相结合的策略,实现了三维空间中的伪弧长数值算法,最后通过数值实验来验证伪弧长算法对于奇异性问题的有效性.     

脆性材料动态破坏过程的数值模拟

根据连续介质力学基本原理和离散元法基本思想，针对脆性材料动态破坏过程问题，推导出基于圆盘单元的正交各向异性离散元二维计算模型，并利用此模型计算了各向异性薄板在冲击载荷下的应力波传播问题。通过将计算结果与LS－DYNA程序相比较，表明了本文提出的离散元模型的精确性。另外，对钢弹侵彻下滑凝土圆板破坏过程这个典型的发生连续介质向非连续介质转化的动力学问题进行了数值模拟和动画显示，并同侵彻实验中出现的几种破坏形式进行了分析和对比，从而证明了本计算方法是计算和模拟脆性材料出现从连续介质向非连续介质转变的动态破坏问题的有力工具。     

Thermoelastic wave propagation in inhomogeneous media

Summary A novel approach to the modeling of thermoelastic wave propagation is presented, based on the thermodynamics of discrete systems. The first novelty includes the representation of integral balance laws for thermoelasticity in terms of contact quantities that describe the nonequilibrium state of elements. The next new aspect is a modification of the recently proposed wave-propagation algorithm, which is used as a tool for determining the contact quantities in a finite-volume scheme for the numerical simulation of two-dimensional thermoelastic wave propagation in inhomogeneous media. Such a modification is needed to provide the satisfaction of the thermodynamic consistency conditions between adjacent discrete elements. Results of computations for certain test problems show the efficiency and physical consistency of the algorithm. Received 11 November 1999; accepted for publication 2 May 2000     

The genetically optimized tunnel-entrance hood

A numerical procedure is investigated for optimizing the design of tunnel-entrance hoods used for controlling the compression wave generated when a high-speed train enters a tunnel. Long hoods are required for long tunnels and train speeds exceeding about 350 km/h. The hood must minimize the maximum pressure gradient across the compression wave-front by taking advantage of the pressure-release provided by open windows distributed along one or both of its walls. The compression wave produced by the train can be evaluated by means of a rapid computational scheme devised and validated against experiment. Optimization is achieved by representing a possible distribution of windows by a binary string. The individuals in an initial, random population of such strings are allowed to ‘mate’ and evolve by ‘natural selection’ through several generations towards an optimal configuration by application of a genetic algorithm. The genetically fittest hood is associated with the minimum possible maximum pressure gradient for prescribed values of the train speed and hood dimensions. The algorithm yields an optimal design from among a theoretically unlimited number of possibilities; it can also supply near-optimal, smoothly varying window distributions (or optimize the variation in width of a long slit-like window in the hood wall) satisfying additional constraints imposed by the designer.     

A NUMERICAL METHOD OF MODELING ELASTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA

The velocity-stress finite-difference method is adopted to simulate the elastic wave propa-gation in azimuthal anisotropic media.The difference grids are completely staggered in the numerical im-plementation.To reduce the computational work,the absorbin8 boundary conditions for anisotropic mediaare introduced first and the corner points are specially treated.Examples show that more accurate resultscan be obtained from the modeling algorithm,which cost much less computational time than the conven-tional methods.Therefore,the algorithm has broad application prospects in engineering.     

On a new algorithm of constructing solitary wave solutions for systems of nonlinear evolution equations in mathematical physics

ＩｎｔｒｏｄｕｃｔｉｏｎＳｏｌｖｉｎｇｎｏｎｌｉｎｅａｒｅｑｕａｔｉｏｎｓｉｓａｌｗａｙｓａｖｅｒｙｉｎｔｅｒｅｓｔｉｎｇｓｕｂｊｅｃｔｆｏｒｍａｔｈｅｍａｔｉｃｉａｎａｎｄｐｈｙｓｉｃｉａｎ ,ｉｎｐａｒｔｉｃｕｌａｒ,ｓｏｌｉｔａｒｙｗａｖｅｓｏｌｕｔｉｏｎｓｆｏｒｎｏｎｌｉｎｅａｒｅｑｕａｔｉｏｎｓａｒｅｏｆｂｏｔｈｔｈｅｏｒｅｔｉｃａｌａｎｄｐｒａｃｔｉｃａｌｉｍｐｏｒｔａｎｃｅ.Ｒｅｃｅｎｔｌｙ ,Ｙａｎ[1]ｏｂｔａｉｎｅｄａｔｒａｎｓｆｏｒｍａｔｉｏｎｄｉｒｅｃｔｌｙｆｒｏｍｔｈｅｆａｍｏ…     

各向异性介质中弹性波传问题的数值解法

通过运用速度－应力有限差分法研究方位各向异性介质中的弹性波传播问题，在计算实施过程中，使用了交错网格技术，为了减少计算量，首次引入了适用于各向异性体的吸收边界条件，并对角点处的吸收做特殊的处理，算例表明，该算法不仅具有较高的精度；与传统方法相比，计算时间也大为缩短，从而可望在实际中获得良好的应用前景。     

连续分层海洋中内波传播的一种数值模式

基于Euler方程,使用有限体积法建立了一种密度为连续分层情况下、适应水深变化的水域中内波传播的数值模式.为了使计算格式能够达到二阶精度,对流项的处理使用了TVD (total variation diminishing)格式.将SIMPLE算法引入连续分层海洋中内波的数值计算,为了简化计算并方便地适应多种TVD格式,在计算预估速度场时采用了显式格式,而没有采用传统的隐式格式;鉴于在原始的SIMPLE算法中没有涉及到由于密度扰动而引起的静水压力场的改变问题,给出了该问题的计算方法.因此改进了SIMPLE算法.出流边界的处理采用阻尼消波和Sommerfeld辐射条件相结合的方式,以使内波得到有效的衰减和释放.将等水深水域的数值解和理论解进行了比较,两者吻合较好;并对存在潜堤时数值计算的不同时刻密度变化的空间分布进行了详细的定性分析.计算结果表明,所建立的数值模式能有效地模拟内波的传播和变形.      

Free tides in two-dimensional uniform-depth basins

The problem of free tidal wave propagation in two-dimensional uniform-depth basins is considered. To solve the problem, a numerical algorithm, a version of the Babenko method, is developed. The algorithm makes it possible to solve the problem in a singly-connected domain admitting conformal mapping onto a circle. A solution can be obtained for both nonrotating and rotating basins. The first natural frequencies are determined for elliptic basins with different eccentricities and angular velocities and certain characteristic natural shapes demonstrating the distinctive features of tidal flows in elliptic basins are constructed. The time evolution of modes in a rotating elliptic basin is studied.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, 2004, pp. 119–130. Original Russian Text Copyright © 2004 by Ivanov.     

Contact problem for the flat crack under two normally incident shear H-waves with wave mode-shifting

The frictional contact interaction of the finite crack edges in an infinite plane is studied for the case of normal incidence of two harmonic shear H-waves with multiple mode-shifted frequencies. Boundary integral equation method and constrained optimization algorithm are used for the problem solution. The forces of contact interaction and displacement discontinuity are analyzed. Influence of the waves frequencies on the stress intensity factor for different normalized wave numbers and coefficients of friction is considered here.     

半空间弹性波双参数强非线性逆散射

半空间弹性波强散射源逆问题是一个实用而又难度很大的强非线性问题,而经典的Born迭代方法只适用于弱非线性问题．基于DBI方法（Distorted-Born Iterative方法）,提出了一种求解弹性波强非线性逆散射问题的迭代方法．在数值模拟运算时利用矩量法进行离散处理,并采用正则化原则避免求解病态矩阵方程。计算结果表明该方法具有较快的收敛速度及较高的精度。     

Simulation of thermoacoustic waves by a pressure-based algorithm for compressible flows

A modified SIMPLEC method which can solve compressible flows at low Mach number is introduced and used to study thermoacoustic waves induced by a rapid change of temperature at a solid wall and alternatingdirection flows generated by thermoacoustic effects in a tapered resonator.The results indicate that the algorithm adopted in this paper can be used for calculating compressible flows and thermoacoustic waves.It is found that the pressure and velocity in the resonator behave as standing waves,and the tapered resonator can suppress highfrequency harmonic waves as observed in a cylindrical resonator.     

A pressure-based Mach-uniform method for viscous fluid flows

A pressure-based, Mach-uniform method is developed by combining the SLAU2 numerical scheme and the higher temporal order pressure-based algorithm. This hybrid combination compensates the limitation of the SLAU2 numerical scheme in the low-Mach number regime and deficiencies of the pressure-based method in the high-Mach number regime. A momentum interpolation method is proposed to replace the Rhie-Chow interpolation for accurate shock-capturing and to alleviate the carbuncle phenomena. The momentum interpolation method is consistent in addition to preserving pressure–velocity coupling in the incompressible limit . The postulated pressure equation allows the algorithm to compute the subsonic flows without empirical scaling of numerical dissipation at low-Mach number computation. Several test cases involving a broad range of Mach number regimes are presented. The numerical results demonstrate that the present algorithm is remarkable for the calculation of viscous fluid flows at arbitrary Mach number including shock wave/laminar boundary layer interaction and aerodynamics heating problem.     

非理想情况下磁化球形重质气体物理爆炸的数值模拟

采用磁流体力学(magnetohydrodynamics, MHD)方程组对非理想情况下磁场对球形重质气体物理爆炸的作用过程进行数值模拟,同时,为了保证每一步中磁场散度为零,引入由12-solve CTU算法衍生出的CTU+CT算法。结果清晰地显示了重质气体团在磁场影响下的整个爆炸过程。在非理想情况下,爆炸过程中气体团界面产生液滴状结构,随着气体团被压缩,不稳定性现象最终得到抑制。从结果中看出,电阻和双极扩散效应会阻碍磁场对气体团的作用,同时,双极扩散效应还会增大磁压力的作用范围。     

A 2D implicit time‐marching algorithm for shallow water models based on the generalized wave continuity equation

This paper builds upon earlier work that developed and evaluated a 1D predictor–corrector time‐marching algorithm for wave equation models and extends it to 2D. Typically, the generalized wave continuity equation (GWCE) utilizes a three time‐level semi‐implicit scheme centred at k, and the momentum equation uses a two time‐level scheme centred at k+12. It has been shown that in highly non‐linear applications, the algorithm becomes unstable at even moderate Courant numbers. This work implements and analyses an implicit treatment of the non‐linear terms through the use of an iterative time‐marching algorithm in the two‐dimensional framework. Stability results show at least an eight‐fold increase in the maximum time step, depending on the domain. Studies also examined the sensitivity of the G parameter (a numerical weighting parameter in the GWCE) with results showing the greatest increase in stability occurs when 1?G/τ_max?10, a range that coincides with the recommended range to minimize errors. Convergence studies indicate an increase in temporal accuracy from first order to second order, while overall error is less than the original algorithm, even at higher time steps. Finally, a parallel implementation of the new algorithm shows that it scales well. Copyright © 2004 John Wiley & Sons, Ltd.     

基于P-M谱的波浪荷载数值模拟及其对跨海斜拉桥的冲击响应规律研究

针对深水桥梁在海洋中面对波浪的冲击问题,以某跨海斜拉桥为例,设定P-M谱为目标谱,依据Welch算法对海浪功率谱进行验证,基于以上理论采用fluent软件建立数值波浪水槽,根据重现期为10年、50年和100年分别计算对应的波高为6.4 m,7.2 m和9.6 m的波浪力,并将波浪力导入斜拉桥进行动力响应计算。综合分析大桥的动力响应, 结果表明,(1) 纵桥向波浪力是由横桥向波浪冲击桥墩出现绕射效应而产生,所以斜拉桥所受横桥向波浪力远大于纵桥向波浪力;(2) 随着波高的增大,斜拉桥的动力响应峰值逐渐增大,主跨跨中和塔顶动力响应峰值随波高的增长率最快,次边跨跨中动力响应峰值随波高的增长率次之,边跨跨中动力响应峰值随波高的增长率最慢;(3) 随着波浪荷载的增大,塔顶横向位移响应增幅显著,塔顶位移响应峰值最大达到0.0598 m。相关研究对大跨度深水桥梁的设计和动力分析有一定的参考意义。     

层状粘接复合结构界面力学参数声导波定征方法

在推导层状粘接复合结构良好粘接及存在弱界面、滑移界面和脱层等几种不同界面条件下声导波的广义频散方程的基础上,分析了界面径向与轴向力学参数对声导波传播特性的影响,进一步提出以频散特性为基础的超声导波定征方法和在最小二乘意义下的反向算法对粘接复合结构层间界面力学参数进行了估计,分析了影响估计准确性的各种因素,研究了超声导波定征方法对粘接复合结构层间力学参数的灵敏度及其在误差传递中的意义。