Electromagnetic backscattering from one-dimensional drifting fractal sea surface I: Wave-current coupled model

To study the electromagnetic backscattering from a one-dimensional drifting fractal sea surface, a fractal sea surface wave-current model is derived, based on the mechanism of wave-current interactions. The numerical results show the effect of the ocean current on the wave. Wave amplitude decreases, wavelength and kurtosis of wave height increase, spectrum intensity decreases and shifts towards lower frequencies when the current occurs parallel to the direction of the ocean wave. By comparison, wave amplitude increases, wavelength and kurtosis of wave height decrease, spectrum intensity increases and shifts towards higher frequencies if the current is in the opposite direction to the direction of ocean wave. The wave-current interaction effect of the ocean current is much stronger than that of the nonlinear wave-wave interaction. The kurtosis of the nonlinear fractal ocean surface is larger than that of linear fractal ocean surface. The effect of the current on skewness of the probability distribution function is negligible. Therefore, the ocean wave spectrum is notably changed by the surface current and the change should be detectable in the electromagnetic backscattering signal.     

近海波-流相互作用的缓坡方程理论体系

在对近海一大类水波模型——缓坡方程的总体考察、分析和判断基础上,引进一个表征近海普遍波-流相互作用机制的算子,进而借助普适的Hamilton水波变分原理,建立了一个关联于时间的近海波-流相互作用缓坡方程理论体系,从形式到内容上达到了一种高度协调和统一.     

无旋性前进重力波传递在均匀流中的Lagrange解析解与试验验证Ⅰ.理论解析解

对于三维空间等深水中,无旋性自由表面周期性规则前进重力波传递在均匀流中的波流场,依质量守恒取一波长的流体质点的运动位移的波长平均高程,所得其标注参数恰为其在原静止水中的位置下,完全以Lagrange方式的参数控制式,解出此波流场至第三阶的全Lagrange形式解且得到检核验证;其中波流交互作用效应存在于Lagrange流速势中,使得波流场中的压力不受均匀流的影响.而Euler形式解所无法描述的流场特性,包括大于前进波周期的流体质点的运动周期,与其受前进波引起的质量传输速度、它们间的关系、及流体质点对其运动周期平均的高程与成因等,都说明是随流体质点所在的高程向下做指数函数样递减;而流体质点的三维空间螺旋曲线式的运动轨迹与烟线,其随均匀流的流向流速而变化的情况,例如其在均匀流于前进波波向有同向的流速分量时,是受流体质点恰在波谷断面处时的流速大小而变的形式,与其在均匀流于前进波波向有反向的流速分量时,则受流体质点恰在波峰断面处时的流速大小而变的形式,有很大不同的倒反形式甚至以封闭曲线形式呈现.最后,说明波流场变成稳定性运动流场时的特性,并证实其在无流时退化成纯前进波的情况.     

Particle image velocimetry measurements of wave-current interaction in a laboratory flume

The hydrodynamics of wave-current interaction is of interest to those concerned with marine and offshore structures. In particular the fluid loading characteristic may be radically altered in a sea state consisting of a mean current flow mixed with freely propagating gravity waves.

The present paper describes water flume experiments, using Particle Image Velocimetry (PIV), executed to examine hydrodynamics of wave-current interaction. A variety of wave and current conditions were investigated to determine the major influences on the combined flow.

This paper describes the experimental procedures used to obtain simultaneous measurements of the resulting wave velocity vectors over an extended region of the wave. It also describes how the directional ambiguity inherent in the basic PIV method was resolved by ‘pulse tagging’ technique.

Velocity vectors under waves at various phase points for different current and wave conditions are presented and compared in some cases with measurements derived using Laser Doppler Anemometry (LDA). The resulting velocity vectors are used to estimate how the mutual interaction, between wave and current, effect the calculation of structural loads using Morrison's equation.     

台风条件下一种新的浪流相互作用参数化方法在耦合模式中的应用

浪致辐射应力及浪致混合效应, 作为海洋上层两种重要的物理过程, 由于当前理论中关于其分层效应及计算公式中的误差还不是非常清楚, 目前还没有在浪流耦合模式中同时考虑. 本文将一种基于浪致辐射应力及浪致混合效应的参数化方案加入由浪(WAVEWATCH Ⅲ)和流模式(POM)组成的耦合模式中. 以台风"格美"(2006)为天气背景, 进行了四组数值试验. 结果表明: 浪致混合效应加强了混合层的湍流, 加大了动量及热量向下传递率, 但对表面流场的规律性改变不明显. 浪致辐射应力传递表面动量通量, 同时, 规律性地改变了水平流场. 当同时考虑两种物理过程时, 其效果叠加. 海流作为一种重要的动力因子, 对有效波高的高度和分布有显著影响, 使得动量转移, 有效波高增高. 模拟时次的第48 h, 流矢量差值的近惯性震荡以及顺时针旋转, 可能是台风靠近时浪流耦合效应加强的信号.     

伴随均匀流作用的有限水深三阶三色三向表面张力-重力波之完备对称解

着眼于"多色多向"这一典型的海洋表面波传播特征而以"三色三向"加以基本地刻画,并纳入普遍的波-流相互作用机理和丰富多样的表面张力波效应,由此给出了有限水深三阶三色三向波运动系统的一完备对称解.这是对目前业已存在的经典、现代"单色、多色多向波理论"的一次充分包容和集中反映. 关键词： 三色三向波 表面张力-重力波 波-流相互作用 完备对称解     

A primitive-variable Riemann method for solution of the shallow water equations with wetting and drying

A Riemann flux that uses primitive variables rather than conserved variables is developed for the shallow water equations with nonuniform bathymetry. This primitive-variable flux is both conservative and well behaved at zero depth. The unstructured finite-volume discretization used is suitable for highly nonuniform grids that provide resolution of complex geometries and localized flow structures. A source-term discretization is derived for nonuniform bottom that balances the discrete flux integral both for still water and in dry regions. This primitive-variable formulation is uniformly valid in wet and dry regions with embedded wetting and drying fronts. A fully nonlinear implicit scheme and both nonlinear and time-linearized explicit schemes are developed for the time integration. The implicit scheme is solved by a parallel Newton-iterative algorithm with numerically computed flux Jacobians. A concise treatment of characteristic-variable boundary conditions with source terms is also given. Computed results obtained for the one-dimensional dam break on wet and dry beds and for normal-mode oscillations in a circular parabolic basin are in very close agreement with the analytical solutions. Other results for a forced breaking wave with friction interacting with a sloped bottom demonstrate a complex wave motion with wetting, drying and multiple interacting wave fronts. Finally, a highly nonuniform, coastline-conforming unstructured grid is used to demonstrate an unsteady simulation that models an artificial coastal flooding due to a forced wave entering the Gulf of Mexico.     

Absorbing boundary conditions for the Euler and Navier–Stokes equations with the spectral difference method

Two absorbing boundary conditions, the absorbing sponge zone and the perfectly matched layer, are developed and implemented for the spectral difference method discretizing the Euler and Navier–Stokes equations on unstructured grids. The performance of both boundary conditions is evaluated and compared with the characteristic boundary condition for a variety of benchmark problems including vortex and acoustic wave propagations. The applications of the perfectly matched layer technique in the numerical simulations of unsteady problems with complex geometries are also presented to demonstrate its capability.     

A nodal discontinuous Galerkin finite element method for nonlinear elastic wave propagation

A nodal discontinuous Galerkin finite element method (DG-FEM) to solve the linear and nonlinear elastic wave equation in heterogeneous media with arbitrary high order accuracy in space on unstructured triangular or quadrilateral meshes is presented. This DG-FEM method combines the geometrical flexibility of the finite element method, and the high parallelization potentiality and strongly nonlinear wave phenomena simulation capability of the finite volume method, required for nonlinear elastodynamics simulations. In order to facilitate the implementation based on a numerical scheme developed for electromagnetic applications, the equations of nonlinear elastodynamics have been written in a conservative form. The adopted formalism allows the introduction of different kinds of elastic nonlinearities, such as the classical quadratic and cubic nonlinearities, or the quadratic hysteretic nonlinearities. Absorbing layers perfectly matched to the calculation domain of the nearly perfectly matched layers type have been introduced to simulate, when needed, semi-infinite or infinite media. The developed DG-FEM scheme has been verified by means of a comparison with analytical solutions and numerical results already published in the literature for simple geometrical configurations: Lamb's problem and plane wave nonlinear propagation.     

Unstructured Adaptive Grid Flow Simulations of Inert and Reactive Gas Mixtures

Unstructured adaptive grid flow simulation is applied to the calculation of high-speed compressible flows of inert and reactive gas mixtures. In the present case, the flowfield is simulated using the 2-D Euler equations, which are discretized in a cell-centered finite volume procedure on unstructured triangular meshes. Interface fluxes are calculated by a Liou flux vector splitting scheme which has been adapted to an unstructured grid context by the authors. Physicochemical properties are functions of the local mixture composition, temperature, and pressure, which are computed using the CHEMKIN-II subroutines. Computational results are presented for the case of premixed hydrogen–air supersonic flow over a 2-D wedge. In such a configuration, combustion may be triggered behind the oblique shock wave and transition to an oblique detonation wave is eventually obtained. It is shown that the solution adaptive procedure implemented is able to correctly define the important wave fronts. A parametric analysis of the influence of the adaptation parameters on the computed solution is performed.     

A large-scale wave-current coupled module with wave diffraction effect on unstructured meshes

Based on the extended mild-slope equation,a large-scale wave module is developed.By combining the eikonal equation and the modified wave action equation,the wave model can account for diffraction in most situations such as in the lee of islands and breakwaters,and using unstructured meshes provides great flexibility for modelling the wave in the complex geomorphology of barriers and islands,also allowing for refinement of the grid resolution within computationally important domains.The numerical implementation of the module is based on the explicit second-order upwind finite-volume schemes in geographic space,the Flux-Corrected Transport(FCT)algorithm in frequency space and the implicit Crank-Nicolson method in directional space.The three-dimensional hydrodynamic module is then modified to couple with the wave model,where the wave readily provides the depth-dependent radiation stress and the wave-induced turbulence coefficient for the current fields,and the wave propagation takes into account the current-induced advection,refraction and diffraction of wave energy and the effect of water level.The applicability of the proposed model to calculate Snell’s Law,wave transformation over the breakwaters and the elliptic shoal,wave propagation over the rip current field and the undertow on a sloping beach is evaluated.Numerical results show that the present model makes better predictions of the near-shore wave propagation and complex three-dimensional(3D)near-shore circulation driven by the waves,considering analytical solutions and experimental values.     

BGK方法在三维非结构网格上的初步应用

将BGK计算方法从二维拓展到三维并且应用于三维非结构网格,具有重要的理论价值和实用价值.采用旋转局部座标的方法,发展了一种针对三维非结构网格的BGK计算方法.在计算过程中,将最小二乘法应用于三维非结构网格的导数计算.对三维激波管和三维欠膨胀垂直冲击射流等两个算例进行了细致分析.这两个算例的计算结果表明,该方法在三维非结构网格上的初步应用是成功的 关键词： 气动BGK方法 三维 非结构网格     

一种求解Euler方程的BGK型非结构化自适应算法

本文构造了一种BGK型二阶非结构化网格自适应算法,用于求解Euler方程.若干一维、二维标准算例表明,本文所构造的算法在模拟复杂流动时具有很高的流场分辨率.     

含运动物体流场计算的自适应非结构二维网格生成方法

采用自动插点的Delaunay方法和局部网格重新生成方法,对含运动边界自适应非结构二维网格生成方法进行了系统研究,实现了带有非定常运动物体流场的数值计算,并进行了数值实验.     

Probing Protein Secondary Structure Using EPR: Investigating a Dynamic Region of Visual Arrestin

One key application of site-directed spin labeling electron paramagnetic resonance (EPR) spectroscopy is the determination of sequence-specific secondary structure in proteins. Regular secondary structure leads to a periodic variation in both side chain motion and solvent accessibility, two properties easily monitored by EPR techniques. Specifically, saturation recovery (SR) EPR spectroscopy has proven to be useful for making accessibility measurements for multiple protein structure populations by determining individual accessibilities and is, therefore, well suited to study the structure of proteins exhibiting multiple conformations in equilibrium. Here we employ both continuous wave and SR EPR spectroscopy in combination to examine the secondary structure of a short sequence showing conformational heterogeneity in visual rod arrestin. The EPR data presented here clearly distinguish between the unstructured loop and the helical structure formed in the crystallographic tetramer of visual arrestin and show that this region is unstructured in solution.     

Multi-dimensional limiting process for hyperbolic conservation laws on unstructured grids

The present paper deals with an efficient and accurate limiting strategy for the multi-dimensional hyperbolic conservation laws on unstructured grids. The multi-dimensional limiting process (MLP) which has been successfully proposed on structured grids is extended to unstructured grids. The basic idea of the proposed limiting strategy is to control the distribution of both cell-centered and cell-vertex physical properties to mimic multi-dimensional nature of flow physics, which can be formulated to satisfy so called the MLP condition. The MLP condition can guarantee high-order spatial accuracy and improved convergence without yielding spurious oscillations. Starting from the MUSCL-type reconstruction on unstructured grids followed by the efficient implementation of the MLP condition, MLP slope limiters on unstructured meshes are obtained.Thanks to its superior limiting strategy and maximum principle satisfying characteristics, the newly developed MLP on unstructured grids is quite effective in controlling numerical oscillations as well as accurate in capturing multi-dimensional flow features. Numerous test cases are presented to validate the basic features of the proposed approach.     

非结构网格上Level Set方程的间断有限元解法

针对间断有限元弱形式难于求解可压缩流场中Level Set方程的问题提出间断有限元强形式,从而在统-框架下解决Level Set方程在可压缩与不可压缩流场中的求解问题.通过非结构网格上采用Legendre-Gauss-Lobatto节点构造基函数,在复杂区域上可以达到任意高阶的精度.将若干-、二、三维算例与已有文献或解析解比较,验证方法追踪自由界面的有效性.结果表明,该方法适合各种情形下Level Set方程求解,易于在复杂区域的非结构网格上实施,精度高、分辨率高且具有高质量守恒性,既能避免重新初始化过程又方便向高维扩展.     

用非结构网格与欧拉方程计算复杂区域的二维流动

提出用Delaunay三角化方法生成非结构网格的一种过程。所生成的网格可用于复杂多连通域内的可压流计算。采用Euler方程和格心有限体积法,研制出程序,给出了算例。     

基于MULTI2D-Z程序的Z箍缩动态黑腔形成过程模拟

对辐射流体力学程序MULTI-2D进行改造,增加磁场演化方程程序模块,自洽地在运动方程模块中增加洛伦兹力,在能量方程模块中增加欧姆加热,将它改造成辐射磁流体力学程序MULTI2D-Z.验证了新增磁场程序模块的可靠性,并发现温度和密度的增大会抑制磁场的扩散,负径向速度梯度的流体对流也会抑制磁场的扩散.利用改造好的MULTI2D-Z程序模拟了峰值为8 MA的脉冲电流驱动的钨丝阵Z箍缩动态黑腔形成过程.得到了X光功率(约30 TW)和能量(约300 k J)、泡沫辐射温度(约120 eV)、箍缩轨迹等模拟结果.在动态黑腔形成过程中,磁场主要分布在钨主体等离子体中;辐射向内传播,烧蚀泡沫柱而使它膨胀;辐射热波在被撞击的泡沫柱中传播,其传播速度比物质温度传播得快,当辐射热波传播到中心轴时泡沫柱中的辐射场变得比较均匀,并且除了冲击波处外辐射温度与物质温度基本上没有分离.这些模拟结果可增强人们对磁场扩散和对流规律以及动态黑腔形成机制的理解,同时表明了MULTI2D-Z程序可成为Z箍缩及其应用的新的程序模拟工具.     

Large-eddy simulation of evaporating spray in a coaxial combustor

Large-eddy simulation of an evaporating isopropyl alcohol spray in a coaxial combustor is performed. The Favre-averaged, variable density, low-Mach number Navier-Stokes equations are solved on unstructured grids with dynamic subgrid scale model to compute the turbulent gas-phase. The original incompressible flow algorithm for LES on unstructured grids by [Mahesh et al., J. Comp. Phys. 197 (2004) 215–240] is extended to include density variations and droplet evaporation. An efficient particle-tracking scheme on unstructured meshes is developed to compute the dispersed phase. Experimentally measured droplet size distribution and size-velocity correlation near the nozzle exit are used as the inlet conditions for the spray. The predictive capability of the LES approach on unstructured grids together with Lagrangian droplet dynamics models to capture the droplet dispersion characteristics, size distributions, and the spray evolution is examined in detail. The mean and turbulent quantities for the gas and particle phases are compared to experimental data to show good agreement. It is shown that for low evaporation rates considered in the present study, a well resolved large-eddy simulation together with simple subgrid models for droplet evaporation and motion provides good agreement of the mean and turbulent quantities for the gas and droplet phases compared to the experimental data. This work represents an important first step to assess the predictive capability of the unstructured grid LES approach applied to spray vaporization. The novelty of the results presented is that they establish a baseline fidelity in the ability to simulate complex flows on unstructured grids at conditions representative of gas-turbine combustors.