低渗透多孔介质渗流动边界模型的解析与数值解

考虑启动压力梯度的低渗透多孔介质非达西渗流模型属于强非线性动边界问题, 分别利用相似变量变换方法和基于空间坐标变换的有限差分方法, 对内边界变压力情况下、考虑启动压力梯度的一维低渗透多孔介质非达西渗流动边界模型进行了精确解析与数值求解研究. 研究结果表明:该动边界模型存在唯一的精确解析解, 且所求得的精确解析解可严格验证数值解的正确性;且当启动压力梯度值趋于零时, 非达西渗流动边界模型的精确解析解将退化为达西渗流情况下的精确解析解. 由求解结果作出的非零无因次启动压力梯度下的地层压力分布曲线表现出紧支性特点, 其与达西渗流模型的有显著不同. 因此, 研究低渗透多孔介质中非稳态渗流问题时, 应该考虑动边界的影响. 研究内容完善了低渗透多孔介质的非达西渗流力学理论, 为低渗透油气藏开发的试井解释与油藏数值模拟技术提供了理论基础.      

一种考虑薄壁散射效应的声学计算模型

采用薄壁边界元/FW-H理论混合方法建立了考虑薄壁声学散射效应的数值计算模型.这种声学计算模型可以预测存在薄壁如风扇机匣、蜗壳等条件下的声波的传播及散射问题.计算模型的建立主要包含噪声源的计算和声源的传播两方面:首先建立FW-H的频域方程,并采用计算流体力学方法计算流场,通过流场数据计算气动噪声源;然后采用薄壁面边界元法计算固壁对声波的散射,并计算声波在固壁散射后的声场分布.数值计算结果和实验结果及经典的叶轮机管道风扇噪声理论进行了对比,结果表明,这种计算模型与理论计算结果及实验结果吻合较好,可以准确的预测机匣壁的散射效应对声源传播的影响.     

二维含裂纹结构与声耦合问题研究

应用分形有限元方法结合边界元方法研究了二维含裂纹结构和声耦合问题.采用二级分形有限元方法对含裂纹的弹性结构体进行离散处理,这样可以使得自由度数大大地减少;无限大外域声场的计算使用边界元方法,可以自动满足无穷远辐射条件.数值仿真算例结果表明:结构声耦合系统的共振频率随着裂纹深度的增加而下降;裂纹附近的声场所受的影响较为明显.     

敷设多孔吸声材料声腔特征值分析的径向积分边界元法

由于Helmholtz方程的基本解是频率的函数,因此传统边界元法在处理声场特征值问题时具有天生的缺陷。本文采用Laplace方程基本解生成积分方程,通过径向积分法将在此过程中产生的域积分项转化为边界积分。此方法克服了传统边界元法系数矩阵对频率的依赖,同时克服了特解积分法对特解的依赖,并通过对表面声导纳的多项式逼近,将敷设多孔吸声材料声腔特征值问题转化为矩阵多项式,从而避免了复杂的非线性求解。通过数值算例验证了算法的有效性。     

A Conceptual Study of Cavity Aeroacoustics Control Using Porous Media Inserts

In this study, an integrated flow simulation and aeroacoustics prediction methodology is applied to testing a sound control technique using porous inserts in an open cavity. Large eddy simulation (LES) combined with a three-dimensional Ffowcs Williams–Hawkings (FW–H) acoustic analogy is employed to predict the flow field, the acoustic sources and the sound radiation. The Darcy pressure – velocity law is applied to conceptually mimic the effect of porous media placed on the cavity floor and/or rear wall. Consequently, flow in the cavity could locally move in or out through these porous walls, depending on the local pressure differences. LES with “standard” subgrid-scale models for compressible flow is carried out to simulate the flow field covering the sound source and near fields, and the fully three-dimensional FW–H acoustic analogy is used to predict the sound field. The numerical results show that applying the conceptual porous media on cavity floor and/or rear wall could decrease the pressure fluctuations in the cavity and the sound pressure level in the far field. The amplitudes of the dominant oscillations (Rossiter modes) are suppressed and their frequencies are slightly modified. The dominant sound source is the transverse dipole term, which is significantly reduced due to the porous walls. As a result, the sound pressure in the far field is also suppressed. The preliminary study reveals that using porous-inserts is a promising technology for flow and sound radiation control.     

Structural–acoustic sensitivity analysis of radiated sound power using a finite element/ discontinuous fast multipole boundary element scheme

The complete interaction between the structural domain and the acoustic domain needs to be considered in many engineering problems, especially for the acoustic analysis concerning thin structures immersed in water. This study employs the finite element method to model the structural parts and the fast multipole boundary element method to model the exterior acoustic domain. Discontinuous higher‐order boundary elements are developed for the acoustic domain to achieve higher accuracy in the coupling analysis. Structural–acoustic design sensitivity analysis can provide insights into the effects of design variables on radiated acoustic performance and thus is important to the structural–acoustic design and optimization processes. This study is the first to formulate equations for sound power sensitivity on structural surfaces based on an adjoint operator approach and equations for sound power sensitivity on arbitrary closed surfaces around the radiator based on the direct differentiation approach. The design variables include fluid density, structural density, Poisson's ratio, Young's modulus, and structural shape/size. A numerical example is presented to demonstrate the accuracy and validity of the proposed algorithm. Different types of coupled continuous and discontinuous boundary elements with finite elements are used for the numerical solution, and the performances of the different types of finite element/continuous and discontinuous boundary element coupling are presented and compared in detail. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamic interaction of a spherical radiator in a fluid-filled cylindrical borehole within a poroelastic formation

Harmonic acoustic radiation from a modally oscillating spherical source positioned at the center of a fluid-filled cylindrical cavity embedded within a fluid-saturated porous elastic formation is studied in an exact manner. The formulation utilizes the Biot theory of dynamic poroelasticity along with the cylindrical to spherical wave-field transformations, and the pertinent boundary conditions to obtain a closed-form series solution. The analytical results are illustrated with a numerical example in which the spherical source, with its polar axis oriented along the main axis of a water-filled borehole and embedded within a water-saturated Ridgefield sandstone formation, is excited in vibrational modes of various orders. The magnitude of the reflected component of acoustic pressure along the axis of the borehole for a pulsating (n = 0), an oscillating (n = 1), and also a multipole (n = 0–3) spherical source as a function of the excitation frequency is calculated and discussed for representative values of the parameters characterizing the system. Special attention is paid to the effects of source excitation frequency, size, surface velocity profile, and internal impedance as well as borehole interface permeability condition on the reflected pressure magnitudes. Limiting cases are considered and fair agreements with well-known solutions are obtained.     

流固耦合地震波动问题的显式谱元模拟方法

流固耦合地震波动问题主要研究由流体和固体构成的复杂系统中地震波传播特性及其规律. 传统模拟方法中一般以声波方程、弹性波方程的数值解分别描述理想流体和弹性固体中的波动, 并实时地处理两种不同性质介质之间的相互耦合作用, 数值格式复杂且限制数值模拟精度与计算效率. 本文采用谱元法结合多次透射公式人工边界条件实现了一种流固耦合地震波动问题的高阶显式数值计算方法. 该方法利用了流固耦合问题统一计算框架,可将饱和多孔介质的Biot波动方程分别退化为理想流体的声波方程和弹性固体的弹性波方程. 通过P波垂直入射的水平成层理想流体-饱和多孔介质-弹性固体场地模型、P波斜入射的不规则层状界面以及任意形状界面的理想流体-饱和多孔介质-弹性固体场地模型等三个算例, 与传递函数法解析解以及集中质量有限元法计算结果进行对比分析, 证明了本文方法的正确性与有效性. 数值模拟结果表明, 本文方法相较传统有限元法可以少得多的节点数量获得更高的数值精度, 并且在较宽的频率范围内都能可靠地模拟出流固耦合系统的动力响应, 充分体现出本文方法兼顾高精度、计算效率和复杂场地建模灵活的特点.      

Harmonic radiation from a liquid-filled spherical acoustic lens with an internal eccentric spherical source

Radiation of sound from a modally vibrating shell-encapsulated (eccentric) spherical source is analyzed in an exact manner using the classical method of separation of variables. The proposed model is a realistic idealization of a spherical acoustic lens with focal point inside the lens when used as a sound projector. The analytical results are illustrated with a numerical example in which the modal acoustic radiation impedance load on the source and the radiated far-field pressure are evaluated for representative values of the parameters characterizing the system. Numerical results clearly illustrate that in addition to frequency, surface velocity distribution and eccentricity of the source, the dynamic interaction of the encapsulating shell can be of great consequence in sound radiation.     

A radiation and convection fluxmeter for high temperature applications

Heat flux is an essential parameter for the diagnostic of thermal systems. In high temperature industrial environment, there are difficulties in measuring incident radiation heat flux as well as in differentiating between the convective and radiative components of heat flux on the heat transfer surface. A new method for heat flux measurement is being developed using a porous sensing element. The gas stream flowing through the porous element is used to measure the heat received by the sensor surface exposed to the hot gas environment. A numerical model of sensor with appropriate boundary condition has been developed in order to perform analysis of possible options regarding its design. The analysis includes: geometry of element, physical parameters of gas and solid and gas flow rate through the porous element. For the optimal selection of parameters, an experimental set-up was designed, including the sensor element with respective cooling and monitoring systems and a high temperature radiation source. The experimental set-up was used to obtain calibration curves for a number of sensors. The linear dependency of the heat flux and respective temperature difference of the gas were verified. The accuracy analysis of the sensor reading has proved high linearity of the calibration curve and accuracy of ±5%.     

冲击声场的时域边界元分析

本文用时域边界元技术从理论上探求机械结构在冲击力作用下的瞬态声辐射规律,预报声场特性,为冲击噪声的研究提供理论分析方法和计算手段。     

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.     

细分曲面边界元法的黏附吸声材料结构拓扑优化分析

等几何分析采用样条基函数构造几何模型和实施变量近似,实现了计算机辅助设计和辅助工程的无缝连接,并已广泛应用于弹性力学、电磁场和位势问题等领域.然而直接采用等几何方法难以构造复杂模型,限制了该方法在大规模实际工程问题上的应用.细分曲面法可用于克服这一问题,该方法对传统模型的离散网格进行细分和拟合操作,构造出极限光滑曲面,连续性更高,对复杂结构的适用性更强.该方法主要有以下优点:(1)适用于任意拓扑结构;(2)数值计算稳定;(3)实施简单;(4)局部细化与连续性控制.由于该方法在复杂结构模型构造方面具有较强的灵活性和便利性,已被广泛应用于航空航天、汽车、动画、游戏制作等建模领域.将细分曲面法与边界元法相结合进行结构声学分析,几何场与物理场均采用箱样条基函数进行插值近似.以黏附吸声材料结构的声散射问题为例,建立吸声材料分布拓扑优化数学模型,并采用移动渐进线算法进行设计变量更新,最终获得最优材料分布.      

广义有限差分法在含阻抗边界空腔声学分析中的应用

声学分析在噪声控制、室内隔音等工程计算中有着重要的作用. 由于现实生活中的声学模型往往伴随着吸声材料, 因此分析含阻抗边界条件的声学问题显得十分必要. 广义有限差分法是一种新型区域型无网格数值离散方法, 该方法基于多元函数泰勒级数展开式和加权最小二乘拟合, 将控制方程中未知参量的各阶偏导数表示为相邻节点函数值的线性组合. 本文首次将广义有限差分法应用于含阻抗边界条件空腔声学问题的分析中, 建立了空腔声场问题的广义有限差分法数值离散格式. 与传统算法相比, 所建立的数值模型具有无需网格剖分和数值积分、计算精度高、适用于大规模声学分析等优点. 通过具有解析解的经典算例, 研究了总节点数目和局部支撑点数目对数值结果的影响, 得到了最大计算频率与节点间距之间关系的经验公式. 此外, 将广义有限差分法应用于无解析解的二维和三维复杂声学模型, 并与COMSOL Multiphysics软件所得的有限元结果进行了比较分析. 数值实验表明, 该算法是一种高效、精确、稳定、收敛的数值模拟方法, 在含阻抗边界空腔声学分析中具有广阔的应用前景.      

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

Experimental study and numerical modeling of brittle fracture of carbonate rock under uniaxial compression

The brittle carbonate rock taken from the Tarim Oilfield is tested in laboratory under uniaxial compression. The acoustic emission (AE) is used to monitor the microcracking activity in rock during the experiment. Moreover, the 3D tomograms of carbonate rock after uniaxial compression are obtained by using CT imaging technology, which indicates that microcracks mutually interconnect and eventually form macroscopic fractures after failure. The PFC2D is used to model the behavior of brittle rock including microcracks propagation. The stress–strain curve and cracks distribution in rock model are obtained from the PFC simulation. The numerical results agree with the experimental test well.     

The application of three-dimensional hydroelastic analysis of ship structures in Pekeris hydro-acoustic waveguide environment简

The hydroelastic analysis and sonoelastic analysis methods are incorporated with the Green＇s function of the Pekeris ocean hydro-acoustic waveguide model to produce a three-dimensional sonoelastic analysis method for ships in the ocean hydro-acoustic environment. The seabed condition is represented by a penetrable boundary of prescribed density and sound speed. This method is employed in this paper to predict the vibration and acoustic radiation of a 1 500 t Small Water Area Twin Hull （SWATH） ship in shallow sea acoustic environment. The wet resonant frequencies and radiation sound source levels are predicted and compared with the measured results of the ship in trial.     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

Development of numerical techniques for near‐field aeroacoustic computations

This work is concerned with the development of a numerical scheme capable of producing accurate simulations of sound propagation in the presence of a mean flow field. The method is based on the concept of variable decomposition, which leads to two separate sets of equations. These equations are the linearised Euler equations and the Reynolds‐averaged Navier–Stokes equations. This paper concentrates on the development of numerical schemes for the linearised Euler equations that leads to a computational aeroacoustics (CAA) code. The resulting CAA code is a non‐diffusive, time‐ and space‐staggered finite volume code for the acoustic perturbation, and it is validated against analytic results for pure 1D sound propagation and 2D benchmark problems involving sound scattering from a cylindrical obstacle. Predictions are also given for the case of prescribed source sound propagation in a laminar boundary layer as an illustration of the effects of mean convection. Copyright © 1999 John Wiley & Sons, Ltd.     

A fractal study of sound propagation characteristics in roughened porous materials

The influence of rough elements of porous materials on sound propagation characteristics in roughened porous materials is investigated. Rough surface topography on the roughened porous materials is simulated by the fractal geometry theory, in which relative roughness is defined clearly as a function of the fractal dimension, porosity, average diameter and diameter ratio of the hexagonal elements on pore wall of the porous materials. Based on the sound propagation model of the smooth porous materials, a sound propagation model of the roughened porous materials is built. The effective density and bulk modulus, acoustic impedance and propagation constant, flow resistivity and sound absorption coefficient of the roughened porous materials are derived and discussed as a function of the relative roughness. It is demonstrated that the sound absorption coefficient of the roughened porous materials is improved as the relative roughness increases. The model predictions for the sound absorption coefficient of the roughened porous materials are well agreed with that of the existing test results.