Study on connections of the MFS, Trefftz method, indirect BIEM and invariant MFS in the three-dimensional Laplace problems containing spherical boundaries

Following the success of a study on the method of fundamental solutions using an image concept [13], we extend to solve the three-dimensional Laplace problems containing spherical boundaries by using the three approaches. The case of eccentric sphere for the Laplace problem is considered. The optimal locations for the source distribution to include the foci in the MFS are also examined by using the image concept in the 3D problems. Whether a free constant is required or not in the MFS is also studied. The error distribution is discussed after comparing with the analytical solution derived by using the bispherical coordinates. Besides, the relationship between the Trefftz bases and the singularity in the MFS for the three-dimensional Laplace problems is also addressed. It is found that one source of the MFS contains several interior and exterior Trefftz sets through a degenerate kernel. On the contrary, one single Trefftz base can be superimposed by some lumped sources in the MFS through an indirect BIEM. Based on this finding, the relationship between the fictitious boundary densities of the indirect BIEM and the singularity strength in the MFS can be constructed due to the fact that the MFS is a lumped version of an indirect BIEM.     

Image solutions for boundary value problems without sources

In this paper, we employ the image method to solve boundary value problems in domains containing circular or spherical shaped boundaries free of sources. two and threeD problems as well as symmetric and anti-symmetric cases are considered. By treating the image method as a special case of method of fundamental solutions, only at most four unknown strengths, distributed at the center, two locations of frozen images and one free constant, need to be determined. Besides, the optimal locations of sources are determined. For the symmetric and anti-symmetric cases, only two coefficients are required to match the two boundary conditions. The convergence rate versus number of image group is numerically performed. The differences of the image solutions between 2D and 3D problems are addressed. It is found that the 2D solution in terms of the bipolar coordinates is mathematically equivalent to that of the simplest MFS with only two sources and one free constant. Finally, several examples are demonstrated to see the validity of the image method for boundary value problems.     

A dynamical model of storm surges along island coasts

A time-independent dynamical model of storm surge along island coasts using orthogonal curvilinear coordinates is presented. The curved annulus between an island coast and an arbitrary deep-water boundary is mapped conformally onto a rectangular image. Two configurations of island coasts are investigated; circular and elliptic coasts. The corresponding coordinates are circular polar and elliptic respectively. The linearized vertically-integrated equations of motion are used to model storm surges with two assumptions: (i) bottom stress is proportional to horizontal transport, and (ii) storm forces are shear stresses on water surface. Analytical solutions are presented for three dynamical cases: (i) a constant-depth basin acted upon by a uniform storm stress, (ii) variable-depth basin acted upon by a uniform-direction variable-magnitude stress, and. (iii) a basin with closed depth contours acted upon by vortex-shaped storm stress. The obtained solutions clarify the relative importance of the various parameters and variables that affect surge height distribution along island coasts. These solutions may be used to test a time-dependent, numerical dynamical storm model.     

Analytical estimation of liquid film thickness in two-phase annular flow using electrical resistance measurement

This paper presents an analytical solution to estimate the liquid film thickness in two-phase annular flow through a circular pipe using electrical resistance tomography. Gas–liquid flow with circular gas core surrounded by a liquid film is considered. Conformal mapping is employed to obtain the analytic solution for annular flow with an eccentric circular gas core. The liquid film thickness for an arbitrary annular flow is estimated by comparing the resistance values for concentric and eccentric annular flows. The film thickness estimation has a good performance when the normalized distance between the gas core center and the flow center is less than 0.2 and the void fraction is greater than 0.4, the estimated error of the normalized thickness is less than 0.04.     

Computing numerically the Green’s function of the half-plane Helmholtz operator with impedance boundary conditions

In this article we compute numerically the Green’s function of the half-plane Helmholtz operator with impedance boundary conditions. A compactly perturbed half-plane Helmholtz problem is used to motivate this calculation, by treating it through integral equation techniques. These require the knowledge of the calculated Green’s function, and lead to a boundary element discretization. The Green’s function is computed using the inverse Fourier operator of its spectral transform, applying an inverse FFT for the regular part, and removing the singularities analytically. Finally, some numerical results for the Green’s function and for a benchmark resonance problem are shown.     

The space fractional diffusion equation with Feller’s operator

This paper investigates the space fractional diffusion equation with fractional Feller’s operator. The Green’s function is obtained by using Fourier transform, and the analytical solutions of some space fractional diffusion equations with initial (or initial and boundary) condition are obtained in terms of Green’s function. In addition, numerical simulations are discussed. The results indicate that the effect range of skewness parameter θ has more effect on probability density than that of parameter α. The results also explain the property of the skewness and long tail in the asymmetry diffusion process.     

Analytical modelling for a three-dimensional hydrofoil with winglets operating beneath a free surface

The current study focuses on establishing a theoretical lifting surface model for predicting the hydrodynamic loads acting on the three-dimensional hydrofoil with winglets, which is considerably influenced by the proximity to the free surface through finding the three-dimensional Green’s function for the planar and vertical horseshoe vortices operating below a free surface. The hydrofoil surface is decomposed into a finite number of elements along the span direction and the chord directions, each of which can then be represented by a horseshoe vortex. The linearized free surface boundary condition is applied to analyze the influence of the free surface on the hydrofoil as well as the winglets. The thickness problem is considered using the source distribution among the hydrofoil and winglets surfaces and the analytical Green’s function that satisfies the linearized free surface boundary condition is used. As a sample application, numerical examples were conducted to show the performance of the hydrodynamic characteristics for the hydrofoil with winglets as a function of the Froude number. It was concluded that there are significant efficiency benefits from using winglets inside the free surface proximity effect. These results are substantiated by the comparison with the available published data.     

A Green’s function model for the analysis of laser heating of materials

An analytical model based on Green’s function method is developed to analyze the temperature distribution and heated regions in a material irradiated by a high-energy laser beam. The model is multi-dimensional, transient and incorporates different types of beam characteristics and boundary conditions. The multi-dimensional integration formulas in the Green’s function solution equation are evaluated using an adaptive numerical integration algorithm. A parametric study is conducted to show the effect of various laser beam parameters and material properties on the laser heating process.     

On the Green’s function for the Helmholtz operator in an impedance circular cylindrical waveguide

This paper addresses the problem of finding a series representation for the Green’s function of the Helmholtz operator in an infinite circular cylindrical waveguide with impedance boundary condition. Resorting to the Fourier transform, complex analysis techniques and the limiting absorption principle (when the undamped case is analyzed), a detailed deduction of the Green’s function is performed, generalizing the results available in the literature for the case of a complex impedance parameter. Procedures to obtain numerical values of the Green’s function are also developed in this article.     

Robertson-Stiff流体在偏心环空中的流动

石油钻井和完井过程中,经常出现泥浆、水泥浆和完井液在偏心环空中的流动。Robertson-Stif流体在描述这些流体的流变性时有很高的精度。文中根据Robertson-Stif流体的流变方程和偏心环空环隙间距的几何关系,利用相似原理,求解了该流体在偏心环形空间中作层流轴向流动时适于工程应用的二元速度分布,以及平均流速、流量和压降的表达式。研究表明:液体在偏心环形空间流动时与在同心环空中流动的最大差异,在于流经偏心环空不同间隙处流速的较大差异。文中还把Robertson-Stif流体的计算结果与采用宾汉模式和幂律模式的计算结果进行了对比。     

The Trefftz method using fundamental solutions for biharmonic equations

In this paper, the Trefftz method of fundamental solution (FS), called the method of fundamental solution (MFS), is used for biharmonic equations. The bounds of errors are derived for the MFS with Almansi’s fundamental solutions (denoted as the MAFS) in bounded simply connected domains. The exponential and polynomial convergence rates are obtained from highly and finitely smooth solutions, respectively. The stability analysis of the MAFS is also made for circular domains. Numerical experiments are carried out for both smooth and singularity problems. The numerical results coincide with the theoretical analysis made. When the particular solutions satisfying the biharmonic equation can be found, the method of particular solutions (MPS) is always superior to the MFS and the MAFS, based on numerical examples. However, if such singular particular solutions near the singular points do not exist, the local refinement of collocation nodes and the greedy adaptive techniques can be used for seeking better source points. Based on the computed results, the MFS using the greedy adaptive techniques may provide more accurate solutions for singularity problems. Moreover, the numerical solutions by the MAFS with Almansi’s FS are slightly better in accuracy and stability than those by the traditional MFS. Hence, the MAFS with the AFS is recommended for biharmonic equations due to its simplicity.     

A discrete fourth-order Lidstone problem with parameters

Various existence, multiplicity, and nonexistence results for nontrivial solutions to a nonlinear discrete fourth-order Lidstone boundary value problem with dependence on two parameters are given, using a symmetric Green’s function approach. An existence result is also given for a related semipositone problem, thus relaxing the usual assumption of nonnegativity on the nonlinear term.     

Three-dimensional dynamic Green’s functions in transversely isotropic tri-materials

An analytical derivation of the elastodynamic fundamental solutions for a transversely isotropic tri-material full-space is presented by means of a complete representation using two displacement potentials. The complete set of three-dimensional point-load, patch-load, and ring-load Green’s functions for stresses and displacements are given, for the first time, in the complex-plane line-integral representations. The formulation includes a complete set of transformed stress-potential and displacement-potential relations in the framework of Fourier expansions and Hankel integral transforms, that is useful in a variety of elastodynamic as well as elastostatic problems. For the numerical computation of the integrals, a robust and effective methodology is laid out. Selected numerical results for point-load and patch-load Green’s functions are presented to portray the dependence of the response on layering, the frequency of excitation, and type of loading.     

Inverse shape and surface heat transfer coefficient identification

In this paper, an inverse geometric problem for the modified Helmholtz equation arising in heat conduction in a fin is considered. This problem which consists of determining an unknown inner boundary of an annular domain and possibly its surface heat transfer coefficient from one or two pairs of boundary Cauchy data (boundary temperature and heat flux) is solved numerically using the meshless method of fundamental solutions (MFS). A nonlinear unconstrained minimisation of the objective function is regularised when noise is added to the input boundary data. The stability of the numerical results is investigated for several test examples with respect to noise in the input data and various values of the regularisation parameters.     

任意分布横向载荷下圆板混合边值问题的一种解析方法

通过引入边界混合函数的概念,选用极坐标下双调和方程的特解序列为试函数,本文提出了一种求解圆板混合边值问题的解析方法.最后通过三个实例计算说明了该方法具有较好的收敛性和稳定性.     

Numerical analysis of the method of fundamental solution for harmonic problems in annular domains

In this study, we investigate the application of the method of fundamental solutions (MFS) to the Dirichlet problem for Laplace's equation in an annular domain. We examine the properties of the resulting coefficient matrix and its eigenvalues. The convergence of the method is proved for analytic boundary data. An efficient matrix decomposition algorithm using fast Fourier transforms (FFTs) is developed for the computation of the MFS approximation. We also tested the algorithm numerically on several problems confirming the theoretical predictions. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Thermal analysis of moving induction heating of a hollow cylinder with subsequent spray cooling: Effect of velocity,initial position of coil,and geometry

In this paper, temperature analysis of the complete process of moving induction heat treatment is performed using numerical methods. A non-linear and transient magneto-thermal coupled problem with a moving coil which is considered as moving heat source, is investigated by an efficient finite-element procedure. A vertical hollow circular cylinder is heated by the moving coil at a given velocity along it, and the heated parts then quenched by a moving water–air spray. The effects of natural convection with air on the both inner and outer surfaces of cylinder, and also radiation of outer surface of cylinder with ambient are taken into account. For quenching of work-piece, a specific kind of atomized spray cooling which utilizes a mixture of water and air with different mass fractions is used. This procedure includes moving boundary conditions, temperature-dependent properties, and change in magnetic permeability of specified alloy at the Curie temperature. Obtained numerical results have been verified by comparison with analytical solutions using Green’s function methods. Also, the effect of velocity, initial position of inductor and inner to outer radius ratio on temperature distribution are investigated.     

The Helmholtz equation with impedance in a half-plane

This Note gives answers to the uniqueness and existence questions for solutions of the Helmholtz equation in an half-plane with an impedance or mixed boundary condition. We deal with unbounded domains which boundaries are unbounded too. The radiation conditions are different from the ones that we found in an usual exterior problem due to the appearance of surface waves. We first compute and study the half-plane Green's function to see how the solutions behave at infinity, and second obtain integral representation for these solutions. To cite this article: M. Duran et al., C. R. Acad. Sci. Paris, Ser. I 340 (2005).     

Elastic fields in two imperfectly bonded half-planes with a thermal inclusion of arbitrary shape

A general method is presented for the rigorous solution of Eshelby’s problem concerned with an arbitrary shaped inclusion embedded within one of two dissimilar elastic half-planes in plane elasticity. The bonding between the half-planes is considered to be imperfect with the assumption that the interface imperfections are uniform. Using analytic continuation, the basic boundary value problem is reduced to a set of two coupled nonhomogeneous first-order differential equations for two analytic functions defined in the lower half-plane which is free of the thermal inclusion. Using diagonalization, the two coupled differential equations are decoupled into two independent nonhomogeneous first-order differential equations for two newly defined analytic functions. The resulting closed-form solutions are given in terms of the constant imperfect interface parameters and the auxiliary function constructed from the conformal mapping which maps the exterior of the inclusion onto the exterior of the unit circle. The method is illustrated using several examples of an imperfect interface. In particular, when the same degree of imperfection is realized in both the normal and tangential directions between the two half-planes, a thermal inclusion of arbitrary shape in the upper half-plane does not cause any mean stress to develop in the lower half-plane. Alternatively, when the imperfect interface parameters are not equal, then a nonzero mean stress will be induced in the lower half-plane by the thermal inclusion of arbitrary shape in the upper half-plane. Detailed results are presented for the mean stress and the interfacial normal and shear stresses caused by a circular and elliptical thermal inclusion, respectively. Results from these calculations reveal that the imperfect bonding condition has a significant effect on the internal stress field induced within the inclusion as well as on the interfacial normal and shear stresses existing between the two half-planes especially when the inclusion is near the imperfect interface.     

Three-dimensional layerwise-finite element free vibration analysis of thick laminated annular plates on elastic foundation

Using a three-dimensional layerwise-finite element method, the free vibration of thick laminated circular and annular plates supported on the elastic foundation is studied. The Pasternak-type formulation is employed to model the interaction between the plate and the elastic foundation. The discretized governing equations are derived using the Hamilton’s principle in conjunction with the layerwise theory in the thickness direction, the finite element (FE) in the radial direction and trigonometric function in the circumferential direction, respectively. The fast rate of convergence of the method is demonstrated and to verify its accuracy, comparison studies with the available solutions in the literature are performed. The effects of the geometrical parameters, the material properties and the elastic foundation parameters on the natural frequency parameters of the laminated thick circular and annular plates subjected to various boundary conditions are presented.