Conditions aux limites absorbantes sur des frontières libres

To ensure the exit of vortices from an open computational domain, new boundary conditions for compressible Navier–Stokes equations, called ‘absorbing conditions’, are introduced. They use a layer around computational domain which eliminates nonphysical perturbations. They are compared to three kinds of boundary conditions based on characteristic relations. Computations of two test cases show the best behaviours of the new boundary conditions. These last ones involve a very small perturbation and can be used to solve aeroacoustic problems.     

含泡沫吸能层防护结构爆炸能量分布的数值模拟研究

采用ANSYS/LS-DYNA软件对复合结构各层中的爆破能量的分布进行了系统的数值模拟研究;由计算得到的复合防护结构中能量的分布规律显示,泡沫混凝土的反射、吸能作用非常明显,改变了防护结构内的能量分布状态,各层的能量分布发生了很大的改变,能量的分布向上转移,其最下层能量仅为无吸能层结构的14％。     

Radiative-convective heat transfer in a plane layer of a selectively absorbing mist

This paper deals with the thermal field in a plane layer of selectively absorbing gas which has been injected into a steady turbulent stream of high-temperature gas flowing around a porous plate. The boundary-value problem in terms of the energy equation reduces to a nonlinear integral equation in terms of a dimensionless temperature, and this equation is solved numerically by the Newton-Kantorovich method. The results are presented on graphs of temperature and thermal flux in the absorbing gas layer as functions of the space coordinate. Such a problem has been analyzed in [1] for the case of an injected gray gas.Translated from Zhurnal Prikladnoi Mekhaniki i Technicheskoi Fiziki, No. 3, pp. 179–182, May–June, 1972.     

Multi-directional and multi-time step absorbing layer for unbounded domain

Variant techniques are proposed for reproducing the elastic wave propagation in an unbounded medium such as the infinite elements, the absorbing boundary conditions or the perfect matched layers. Here, a simplified approach is adopted by considering absorbing layers characterized by the viscous Rayleigh matrix as studied by Semblat et al. [16] and Rajagopal et al. [14]. Here, further improvements to this procedure are provided. First, we start by establishing the strong form for the elastic wave propagation in a medium characterized by the Rayleigh matrix. This strong form will be used for deriving optimal conditions for damping out in the most efficient way the incident waves while minimizing the spurious reflected waves at the interface between the domain of interest and the Rayleigh damping layer. A procedure for designing the absorbing layer is proposed by targeting a performance criterion expressed in terms of logarithmic decrement of the wave amplitude in the layer thickness. Second, the GC subdomain coupling method, proposed by Combescure and Gravouil [9], is introduced for enabling the choice of any Newmark time integration schemes associated with different time steps depending on subdomains. When wave propagation is predicted by an explicit time integrator, the subdomain strategy is of great interest because it enables a different time integrator for the absorbing layer to be adopted. An external coupling software, based on the GC method, is used to carry out multi=time step explicit/implicit co-computations, making interact in time an explicit FE code (Europlexus) for the domain of interest, with an implicit FE code (Cast3m) handling the absorbing boundary layers. The efficiency of the approach is shown in 1D and 2D elastic wave propagation problems.     

Calculation of inviscid radiating gas flow past a blunt body

During hypersonic gas flow past a blunt body with a velocity on the order of the escape velocity or more, the gas radiation in the disturbed region behind the shock wave becomes the primary mechanism for aerodynamic heating and has a significant effect on the distribution of the gasdynamic parameters in the shock layer. This problem has been considered from different points of view by many authors. A rather complete review of these studies is presented in [1–4].In earlier studies [5, 6] the approximation of bulk emission was used. In this approximation, in order to account for the effect of radiative heat transfer a term is added in the energy equation which is equivalent to the body efflux, whose magnitude depends on the local thermodynamic state of the gas. However, the use of this assumption to solve the problem of inviscid flow past a blunt body leads to a singularity at the body [7, 8]. To eliminate the singularity, account is taken of the radiation absorption in a narrow wall layer [7], or the concept of a viscous and heat-conductive shock layer is used [8]. A further refinement was obtained by Rumynskii, who considered radiation selectivity and studied the flow of a radiating and absorbing gas in the vicinity of the forward stagnation point of a blunt body.In the present paper we study the distribution of the gasdynamic parameters in the shock layer over the entire frontal surface of a blunt body in a hypersonic flow of a radiating and absorbing gas with account for radiation selectivity.     

Numerical boundary conditions for unsteady transonic flow calculations

In calculations of transonic flows it is necessary to limit the domain of computation to a size that is manageable by computers. At the boundary of the computational domain, boundary conditions are required to ensure a unique solution. Since wave solutions exist in the unsteady transonic flow field, incorrect boundary conditions may result in spurious reflections from the computational boundary. This may introduce errors into the solution. To prevent the spurious reflections, absorbing boundary conditions are often used on the computational boundary. In this paper we describe a method to derive absorbing boudary conditions for transonic calculations. We demonstrate both theoretically and numerically that the use of the absorbing boundary conditions will reduce the spurious reflections in the calculation.     

Numerical simulation of stabilization of the boundary layer on a surface with a porous coating in a supersonic separated flow

Stability of a supersonic (M _∞ = 5.373) boundary layer with local separation in a compression corner with a passive porous coating partly absorbing flow perturbations is considered by solving two-dimensional Navier-Stokes equations numerically. The second mode of disturbances of a supersonic boundary layer is demonstrated to be the most important one behind the boundary-layer reattachment point. The possibility of effective stabilization of these disturbances behind the reattachment point with the use of porous coatings is confirmed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 2, pp. 39–47, March–April, 2007.     

二维稳态辐射声场的光滑有限元-完美匹配层解法

针对外场声学有限元计算精度偏低的问题, 将光滑有限元技术引入到二维稳态辐射声场预测中, 提出了光滑有限元-完美匹配层解法. 该解法采用完美匹配层截断声场计算域, 并将其离散为等参四边形单元, 采用指数吸收函数实现完美匹配层内参数坐标和笛卡尔坐标的映射关系, 采用光滑声压梯度技术计算辐射声场刚度矩阵, 将形函数梯度的域内积分转换为形函数域边界积分. 某汽车二维声腔辐射声场的数值分析结果表明, 与标准有限元-完美匹配层相比, 光滑有限元-完美匹配层解法在完美匹配层内的声波吸收效果更好, 在计算域内的数值计算精度更高, 具有良好的工程应用前景.      

Minimizing the thickness of an inhomogeneous layer at a given reflection coefficient for a monochromatic wave

The problem of wave absorption in an inhomogeneous impedance layer is a familiar one in acoustics and electrodynamics. In practice, impedance systems are designed with account for known optimum requirements, the most important of which is the requirement that the absorbing layer be of minimum thickness. The corresponding mathematical problems for sound absorbers in air were treated by various artificial approaches toward the end of the 1930's by Malyuzhinets, who attracted the attention of the authors, and by Svirskii, in a dissertation (Moscow State University, 1943). Similar (and more general) optimization problems can be studied systematically when they are treated as Mayer-Bolza variational problems. That point of view is adopted in this paper, in which the thickness of an inhomogeneous layer on which a plane monochromatic wave is incident normally is minimized.     

Berechnung der Vergrößerung der Austauschflächen und Ventilationsleistungen beim Übergang von Wasser- auf Luftkühlung

Through use of a simple equation for the economic optimization of heat exchangers used to reject heat to the environment, a functional relationship is developed for heat transfer area and required system ventilator power. Thermal resistance has been neglected between the heat exchanger and heat rejecting fluid. It is found that for a system in which the heat absorbing fluid is air, both the required heat transfer area and ventilator power are 76 times greater than in a system in which the heat absorbing fluid is water.     

Natural convection due to solar radiation over a non-absorbing plate with and without heat losses

Natural convection over a non-reflecting, non-absorbing, ideally transparent semi-infinite vertical flat plate due to absorption of incident radiation (solar radiation) is considered. The absorbed radiation acts as a distributed source which initiates buoyancy-driven flow and convection in the absorbing layer. The plate when heated by the absorbing fluid loses heat to the surroundings from its external side. Solution of the governing equations of the flow under these circumstances is non-similar because of both the heat source term in the energy equation and the temperature boundary condition at the plate. A local non-similar technique is used to obtain solutions for a wide range of the dimensionless distance along the plate and of the dimensionless loss coefficient to the surroundings. The results show that the temperature distribution has a maximum temperature in the depth of the fluid rather than on the plate. A new definition for a local heat transfer coefficient between the plate and the absorbing fluid is introduced which is based on the local maximum temperature rise in the fluid. A formula to calculate this heat transfer coefficient is given for the anticipated range of the loss coefficient.     

The Fokker–Planck Equation with Absorbing Boundary Conditions

We study the initial-boundary value problem for the Fokker–Planck equation in an interval with absorbing boundary conditions. We develop a theory of well-posedness of classical solutions for the problem. We also prove that the resulting solutions decay exponentially for long times. To prove these results we obtain several crucial estimates, which include hypoellipticity away from the singular set for the Fokker–Planck equation with absorbing boundary conditions, as well as the Hölder continuity of the solutions up to the singular set.     

Complex frequency shifted convolution PML for FDTD modelling of elastic waves

The perfectly matched layer (PML) is nowadays considered as the best optimum absorbing boundary condition available. However, the PML with the classical stretching tensor has certain limitations. Strangely, these limitations have rarely been addressed in elastic wave modelling. For example, substantial reflections occur when strong evanescent waves are propagating parallel to the interface. To circumvent problems like this, the complex frequency shifted stretching tensor has been introduced in electromagnetic modelling. In this paper we show that the convolution PML with this stretching tensor as used in electromagnetic modelling can be adapted for elastic wave modelling. Numerical results of a model where the presence of evanescent waves is predominant show that the PML based on the complex frequency shifted stretching tensor can improve the performance of the absorbing boundary layer considerably.     

Laser–induced thermocapillary deformation of a thin liquid layer

An equation that describes the profile of laser–induced thermocapillary deformation of the thin layer of an absorbing liquid is obtained. Numerical estimates of the deformation profile are in agreement with experimental data. It is shown that with increase in the radiation power, the deformation depth of the liquid layer increases and leads to its rupture.     

Numerical modeling of the stabilization of a supersonic flat-plate boundary layer by a porous coating

On the basis of a numerical solution of the two-dimensional Navier-Stokes equations, the stability and the receptivity of a supersonic (M_∞ = 6) boundary layer on a flat plate with a passive porous coating partially absorbing flow disturbances is studied. The results of direct numerical simulation are in good agreement with the data of the linear stability theory. The studies confirm the possibility of effectively stabilizing the second mode of the supersonic boundary layer using porous coatings.     

On radiative heat transfer in a plane layer of an absorbing medium

Recently there has arisen increased interest in the study of radiative heat transfer between geometrically simple systems, both as autonomous problems and as elements of more complex problems.Problems of this kind have been treated by many authors [1–111 who have considered gray, diffusely emitting and absorbing boundaries and gray nonscattering media. In most cases these investigations were restricted either to the derivation of approximate formulas for the net radiative flux, without an exact analysis of the temperature distribution in the layer [5–7], or to numerical computation [1–4], In the latter case, with the exception of [8], which contains a numerical analysis for the case of optical symmetry, no attempt was made to analyze the effect of the optical properties of the boundaries on the temperature field in the layer.These papers can be divided into two groups according to the method of analysis used. The first group includes papers based on the integral equations of radiative transfer, with the corresponding integral analytical methods [1, 2], Similar in nature are [3, 4] which use the slab method, applicable to electrical-analog computation, as well as a recent paper [8] based on probability methods.The second group of papers [5–7] is based on the so-called differential methods. Of particular interest is [7], which develops these methods to an advanced degree. In several papers the problem of radiative transfer is analyzed in conjunction with more complex problems (cf., e.g. [10, 11]).In the present work we shall attempt to carry out an approximate analytical study of problems connected with radiative heat transfer in a plane layer of an absorbing, emitting, nonscattering gray medium with temperature-independent optical properties. The layer is bounded by two parallel, diffusely emitting and diffusely reflecting, isothermal, gray planes.The paper presents the fundamental formulation of the problem, which consists in: (a) the determination of the net heat flux on the basis of given temperature distribution (direct formulation), and (b) the determination of the temperature distribution on the basis of given distribution of the net radiative heat source per unit volume and boundary temperatures (inverse formulation). The analysis is based on integral methods appropriate to the integral equations which represent the net total and hemispherical radiation flux densities [12].The author would like to thank S. S. Kutateladze for his interest in this work.     

黏弹性多层介质中SH波动的一种吸收边界条件

提出一种高精度时域吸收边界条件,与有限元法结合用于模拟瞬态标量SH波在达朗贝尔黏弹性多层介质中传播问题.建立时域吸收边界条件的过程是:首先将半无限域沿着竖向半离散得到半离散的位移方程以及人工边界处的力-位移关系,再通过引入模态分解,将物理空间下的量转化到模态空间,从而获得半无限域模态空间下的频域动力刚度;其次采用一种在全频范围内收敛的连分式精确逼近单层介质模态空间下标量形式的频域动力刚度,将标量连分式扩展为矩阵形式用来表示多层介质的频域动力刚度;最后通过引入辅助变量技术,将模态空间下基于连分式的频域动力刚度关系转化为时域吸收边界条件,进一步转换到物理空间后得到物理空间下的时域吸收边界条件.单层介质和五层介质的数值算例表明,建立的高精度时域吸收边界条件对于达朗贝尔黏弹性单层介质是精确且稳定的;对于达朗贝尔黏弹性多层介质,为了保证其高精度特性,需要将人工边界放置在距离感兴趣区域约为0.5倍无限域高度的位置处.     

Perfectly matched layer boundary condition for two-dimensional Euler equations in generalized coordinate system

In the past, perfectly matched layer (PML) equations have been constructed in Cartesian and spherical coordinates. In this article, the focus is on the development of a PML absorbing technique for treating numerical boundaries, especially those with unbounded domains, in a generalized coordinate system for a flow in an arbitrary direction. The PML equations for two-dimensional Euler equations are developed in split form through a space–time transformation involving a complex variable transformation with the application of a pseudo-mean-flow in the PML domain. A numerical solver is developed using conventional numerical schemes without employing any form of filtering or artificial dissipation to solve the governing PML equations for two-dimensional Euler equations in a generalized coordinate system. Physical domains of arbitrary shapes are considered and numerical simulations are carried out to validate and demonstrate the effectiveness of the PML as an absorbing boundary condition in generalized coordinates.     

Nonstationary (transient) radiative-conductive heat transfer in a plane layer of gray absorbing medium

The kinetics of the heating of a plane layer of gray absorbing medium by radiative-conductive heat transfer are considered. The nonstationary energy equation is reduced to a nonlinear integral equation by means of a Green's function, and this is solved numerically by the Newton method. The results of the solution are presented in the form of the temperature fields in the layer for various values of the defining parameters (optical thickness, radiative-conductive heat-transfer criterion, heat-transfer criterion at the boundaries).Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 156–159, January–February, 1972.     

Determining the boundary layers in the plane problem for three-layer strips. Part 1

We propose an exact solution of the problem on a boundary layer (a stress-strain state decreasing away from the boundary) for three-layer strips (rods) whose layers are made of different materials. We use the asymptotic integration method to obtain boundary eigenfunctions and a characteristic equation for the parameter describing the boundary layer decay rate. We study how the middle layer material affects the boundary layer extent.