Non-linear free surface flows and an application of the Orlanski boundary condition

The focus of this paper is the analysis of spatially two-dimensional non-linear free surface problems. The critical aspects of the problem concern the treatment of the non-linear free surface, the body boundary condition for large motions and the imposition of suitable radiation conditions. To address such complexities, time domain simulation was chosen as the method of analysis. With the use of a finite domain for simulation, a major concern is with the radiation condition to be applied at the open or truncation boundary. For the two-dimensional problem at hand, no theoretical radiation conditions are known to exist. An extension of the Orlanski open boundary condition, based on phase velocity determination at the free surface, is proposed. Three categories of problems were analysed using numerical simulation-namely, freely moving steep waves, waves over a submerged body and forced body motion. Simulation results have been compared with linear theory and experiments.     

A formal finite element approach for open boundaries in transport and diffusion ground-water problems

In the application of the finite element method to diffusion and convection-dispersion equations over a ground-water domain, the Galerkin technique was used to incorporate Neumann (or second-type) and Cauchy (or third-type) boundary conditions. While mass movement through open boundaries is a priori unknown, these boundaries are usually treated as a zero Neumann condition at some far distance from the domain of interest. Nevertheless, cheaper and better solutions can be obtained if these unknown conditions are adequately incorporated in the weak formulation and in the transient solution schemes (open boundary condition). Theoretical and numerical proofs are given of the equivalences between this approach and a ‘well-posed’ problem in a semi-infinite domain with a zero Neumann condition at a boundary placed at infinity. Transport and diffusion equations were applied in one dimension to show the numerical performances and limitations of this procedure for some linear and non-linear problems. No a priori limitations are foreseen in order to find similar solutions in two or three dimensions. Thus the spatial discretization in the proximity of open boundaries could be drastically reduced to the domain of interest.     

Méthode couplée CVFE/BE pour les calculs acoustiques en domaine non borné

A new method based on the flux conservation for acoustic computations in unbounded domain is presented. The domain is divided in outer and inner regions by a fictitious boundary. The inner problem is solved by a control volume finite element method while the outer flow is treated with a boundary element method. The coupling of the two sub-domains is performed using an appropriate expression for the outgoing fluxes satisfying the Sommerfeld condition. A two-dimensional configuration with a circular interface is considered to show that the method is effective. Results of this method are compared to an analytical solution and to results from the literature.     

Development of a numerical wave tank for analysis of nonlinear and irregular wave field

A numerical wave-absorption filter has been developed for an open boundary condition in the analysis of nonlinear and irregular wave evolution. The filter is composed of a simulated sponge layer and Sommerfeld's radiation condition at the outer edge of the layer. The wave-absorption characteristics of the filter have been investigated by applying the linear potential theory and a two-dimensional nonlinear boundary element model. In both cases, the filter is found to he applicable for a wide range of wave parameters. In order to realize an idealized “numerical wave tank”, the present model also incorporates a nonreflective wave generator in the computational domain composed of a series of vertically aligned point sources. Numerous numerical experiments demonstrate that the present approach is effective in generating an arbitrary wave profile without reflection not only at the open boundaries but also at the wave generator.     

Effects of radiation and heat source on MHD flow of a viscoelastic liquid and heat transfer over a stretching sheet

We study the MHD flow and also heat transfer in a viscoelastic liquid over a stretching sheet in the presence of radiation. The stretching of the sheet is assumed to be proportional to the distance from the slit. Two different temperature conditions are studied, namely (i) the sheet with prescribed surface temperature (PST) and (ii) the sheet with prescribed wall heat flux (PHF). The basic boundary layer equations for momentum and heat transfer, which are non-linear partial differential equations, are converted into non-linear ordinary differential equations by means of similarity transformation. The resulting non-linear momentum differential equation is solved exactly. The energy equation in the presence of viscous dissipation (or frictional heating), internal heat generation or absorption, and radiation is a differential equation with variable coefficients, which is transformed to a confluent hypergeometric differential equation using a new variable and using the Rosseland approximation for the radiation. The governing differential equations are solved analytically and the effects of various parameters on velocity profiles, skin friction coefficient, temperature profile and wall heat transfer are presented graphically. The results have possible technological applications in liquid-based systems involving stretchable materials.     

The non-linear dynamic response of the Euler-Bernoulli beam with an arbitrary number of switching cracks

In this study the non-linear dynamic response of the Euler-Bernoulli beam in presence of multiple concentrated switching cracks (i.e. cracks that are either fully open or fully closed) is addressed. The overall behaviour of such a beam is non-linear due to the opening and closing of the cracks during the dynamic response; however, it can be regarded as a sequence of linear phases each of them characterised by different number and positions of the cracks in open state. In the paper the non-linear response of the beam with switching cracks is evaluated by determining the exact modal properties of the beam in each linear phase and evaluating the corresponding time history linear response through modal superposition analysis. Appropriate initial conditions at the instant of transition between two successive linear phases have been considered and an energy control has been enforced with the aim of establishing the minimum number of linear modes that must be taken into account in order to obtain accurate results. Some numerical applications are presented in order to illustrate the efficiency of the proposed approach for the evaluation of the non-linear dynamic response of beams with multiple switching cracks. In particular, the behaviour under different boundary conditions both for harmonic loading and free vibrations has been investigated.     

A complete boundary integral formulation for steady compressible inviscid flows governed by non-linear equations

A complete boundary integral formulation for steady compressible inviscid flows governed by non-linear equations is established by using the specific mass flux as a dependent variable. Thus, the dimensionality of the problem to be solved is reduced by one and the computational mesh to be generated is needed only on the boundary of the domain. It is shown that the boundary integral formulation developed in this paper is equivalent of the results of distributions of the fundamental solutions of the Laplacian operator equation with a different order along the boundaries of the domain. Hence, we have succeeded in establishing the fundamental-solution method for compressible inviscid flows governed by non-linear equations.     

Discrete non‐local absorbing boundary condition for exterior problems governed by Helmholtz equation

The finite element method is employed to approximate the solutions of the Helmholtz equation for water wave radiation and scattering in an unbounded domain. A discrete, non‐local and non‐reflecting boundary condition is specified at an artificial external boundary by the DNL method, yielding an equivalent problem that is solved in a bounded domain. This procedure formulates a boundary value problem in a bounded region by imposing a relation in the discrete medium between the nodal values at the two last layers. For plane geometry, this relation can be found by straightforward eigenvalue decomposition. For circular geometry, the plane condition is applied at the external layer and this condition is condensed through a structured annular region, resulting in a condition at an inner radius. Exterior problems with a bounded internal physical obstacle are considered. It is well‐known that these kind of problems are well‐posed, and have a unique solution. Numerical studies based on standard Galerkin methodology examine the dependence of the DNL condition with respect to the circular annular region width. The DNL condition is compared with local boundary conditions of several orders. Numerical examples confirm the important improvement in accuracy obtained by the DNL method over standard conditions. Copyright © 1999 John Wiley & Sons, Ltd.     

A coupled fluid-elasticity model for the wave forcing of an ice-shelf

A mathematical model for predicting the vibrations of ice-shelves based on linear elasticity for the ice-shelf motion and potential flow for the fluid motion is developed. No simplifying assumptions such as the thinness of the ice-shelf or the shallowness of the fluid are made. The ice-shelf is modelled as a two-dimensional elastic body of an arbitrary geometry under plane-strain conditions. The model is solved using a coupled finite element method incorporating an integral equation boundary condition to represent the radiation of energy in the infinite fluid. The solution is validated by comparison with thin-beam theory and by checking energy conservation. Using the analyticity of the resulting linear system, we show that the finite element solution can be extended to the complex plane using interpolation of the linear system. This analytic extension shows that the system response is governed by a series of singularities in the complex plane. The method is illustrated through time-domain simulations as well as results in the frequency domain.     

应用边界积分法求圆形夹杂问题的解析解

边界元方法作为一种数值方法, 在各种科学工程问题中得到了广泛的应用.本文参考了边界元法的求解思路, 从Somigliana等式出发, 利用格林函数性质,得到了一种边界积分法, 使之可以用来寻求弹性问题的解析解.此边界积分法也可以从Betti互易定理得到. 应用此新方法, 求解了圆形夹杂问题.首先设定夹杂与基体之间完美连接, 将界面处的位移与应力按照傅里叶级数展开,根据问题的对称性与三角函数的正交性来简化假设, 减少待定系数的个数.其次选择合适的试函数(试函数满足位移单值条件以及无体力的线弹性力学问题的控制方程),应用边界积分法, 求得界面处的位移与应力的值. 然后再求解域内位移与应力.得到了问题的精确解析解, 当夹杂弹性模量为零或趋向于无穷大时,退化为圆孔或刚性夹杂问题的解析解. 求解过程表明,若问题的求解区域包含无穷远处时, 所取的试函数应满足无穷远处的边界条件.若求解区域包含坐标原点, 试函数在原点处位移与应力应是有限的.结果表明了此方法的有效性.      

Impact of thermal radiation on MHD slip flow over a flat plate with variable fluid properties

The influence of partial slip, thermal radiation and temperature dependent fluid properties on the hydro-magnetic fluid flow and heat transfer over a flat plate with convective surface heat flux at the boundary and non-uniform heat source/sink is studied. The transverse magnetic field is assumed as a function of the distance from the origin. Also it is assumed that the fluid viscosity and the thermal conductivity vary as an inverse function and linear function of temperature respectively. Using the similarity transformation, the governing system of non-linear partial differential equations are transformed into similarity non-linear ordinary differential equations and are solved numerically using symbolic software MATHEMATICA 7.0. The numerical values obtained within the boundary layer for the dimensionless velocity, temperature, skin friction coefficient and the Nusselt number are presented through graphs and tables for several sets of values of the parameters. The effects of various physical parameters on the flow and heat transfer characteristics are discussed from the physical point of view.     

物面引射对滑移边界条件的影响

讨论壁面有引射的滑移边界条件的提法．利用Ｃｈａｐｍａｎ Ｅｎｓｋｏｇ速度分布函数,通过分析Ｋｎｕｄｓｅｎ层外缘和壁面处的质量、动量和能量等通量守恒,得到了有壁面引射的多组元气体有催化反应的壁面滑移边界条件方程组．利用粘性激波层方程和所得的边界条件对钝头体绕流驻点区的流场进行了计算,讨论了物面引射对边界上诸量及流场的影响     

浅埋圆形孔洞附近的半圆形凸起对SH波的散射

采用"契合"的思想，给出了地下孔洞与地面上的半圆形凸起地形对SH波散射问题的 解答. 将整个求解区域分割成两部分来处理. 其一为包括半圆形凸起地形在内的一个圆形区 域I，其余为区域II. 在区域I和II中分别构造位移解，并在两个区域的"公共边界"上 实施"契合". 在区域I中构造一个上半部边界应力为零，而其余部分位移、应力任意的驻 波解，在区域II中构造出半圆形凹陷和浅埋圆孔的散射波，且要求它满足水平界面上应力 为零的约束条件. 然后再通过移动坐标，满足"公共边界"的"契合"条件和地下孔洞的边 界条件，建立起求解该问题的无穷代数方程组. 最后，给出了分析例题和数值结果，并 对其进行了讨论.     

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%.     

流固耦合介质轴对称动力问题解法的改进

用直接求解常微分方程组解文［１］所得的控制方程，减少了传递矩阵计算工作量，避免了子阵求逆，使问题的求解得到了简化     

On the crack face boundary conditions in electromechanical fracture and an experimental protocol for determining energy release rates

The fracture mechanics of electromechanical materials has been investigated for well over a decade, yet there still exists controversy over the appropriate crack face boundary conditions for non-conducting cracks. In this paper an experimental protocol for measuring the energy release rate in a non-linear reversible electromechanical body is proposed and summarized. The potential results from the proposed experimental approach are capable of shedding light on the true physical nature of the conditions prevailing at the crack surface and in the space within the crack. The experimental procedure is simulated numerically for a linear piezoelectric specimen in a four point bending configuration subjected to electrical loading perpendicular to the crack. The focus of these investigations is on a comparison between the commonly used exact crack face boundary condition and the recently proposed energetically consistent boundary conditions. To perform the numerical calculation with a wide range of electrical and mechanical loadings, two efficient finite element formulations are presented for the general analysis of crack problems with non-linear crack face boundary conditions. Methods for the numerical determination of the crack tip energy release rate and the simulation of the experimental method for obtaining the total energy release rate are developed. Numerical results for the crack tip and total energy release rate are given for both the exact and energetically consistent boundary conditions. It is shown that the crack tip energy release rate calculated under energetically consistent boundary conditions is equal to the total energy release rate generated from the simulated experimental method. When the exact boundary conditions are used, there is no such agreement.     

On the immersed boundary‐lattice Boltzmann simulations of incompressible flows with freely moving objects

For simulating freely moving problems, conventional immersed boundary‐lattice Boltzmann methods encounter two major difficulties of an extremely large flow domain and the incompressible limit. To remove these two difficulties, this work proposes an immersed boundary‐lattice Boltzmann flux solver (IB‐LBFS) in the arbitrary Lagragian–Eulerian (ALE) coordinates and establishes a dynamic similarity theory. In the ALE‐based IB‐LBFS, the flow filed is obtained by using the LBFS on a moving Cartesian mesh, and the no‐slip boundary condition is implemented by using the boundary condition‐enforced immersed boundary method. The velocity of the Cartesian mesh is set the same as the translational velocity of the freely moving object so that there is no relative motion between the plate center and the mesh. This enables the ALE‐based IB‐LBFS to study flows with a freely moving object in a large open flow domain. By normalizing the governing equations for the flow domain and the motion of rigid body, six non‐dimensional parameters are derived and maintained to be the same in both physical systems and the lattice Boltzmann framework. This similarity algorithm enables the lattice Boltzmann equation‐based solver to study a general freely moving problem within the incompressible limit. The proposed solver and dynamic similarity theory have been successfully validated by simulating the flow around an in‐line oscillating cylinder, single particle sedimentation, and flows with a freely falling plate. The obtained results agree well with both numerical and experimental data. Copyright © 2016 John Wiley & Sons, Ltd.     

Effect of electric boundary conditions on crack propagation in ferroelectric ceramics

In this paper, the effect of electric boundary conditions on Mode I crack propagation in ferroelectric ceramics is studied by using both linear and nonlinear piezoelectric fracture mechanics. In linear analysis, impermeable cracks under open circuit and short circuit are analyzed using the Stroh formalism and a rescaling method. It is shown that the energy release rate in short circuit is larger than that in open circuit. In nonlinear analysis, permeable crack conditions are used and the nonlinear effect of domain switching near a crack tip is considered using an energy-based switching criterion proposed by Hwang et al.（Acta Metal. Mater.,1995）. In open circuit, a large depolarization field induced by domain switching makes switching much more diffcult than that in short circuit. Analysis shows that the energy release rate in short circuit is still larger than that in open circuit, and is also larger than the linear result. Consequently,whether using linear or nonlinear fracture analysis, a crack is found easier to propagate in short circuit than in open circuit, which is consistent with the experimental observations of Kounga Njiwa et al.（Eng. Fract. Mech., 2006）.