Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

A new type of boundary integral-equation for plane problems of elasticity including rotational forces

Based on the idea proposed by Hu [Scientia Sinica Series A XXX, 385-390 (1987)], a new type of boundary integral equation for plane problems of elasticity including rotational forces is derived and its boundary element formulation is presented. Numerical results for a rotating hollow disk are given to demonstrate the accuracy of the new type of boundary integral equation.

Interaction of linear waves with infinitely long horizontal cylinders studied by boundary element method

The two-dimensional problems concerning the interaction of linear water waves with cylinders of arbitrary shape in two-layer deep water are investigated by use of the Boundary Integral Equation method (BIEM). Simpler new expressions for the Green functions are derived, and verified by comparison of results obtained by BIEM with these by an analytical method. Examined are the radiation and scattering of linear waves by two typical configurations of cylinders in two-layer deep water. Hydrodynamic behaviors including hydrodynamic coefficients, wave forces, reflection and transmission coefficients and energies are analyzed in detail, and some interesting physical phenomena are observed.

A combined BIE-FE method for the Stokes equations

A numerical algorithm for the biharmonic equation in domains with piecewise smooth boundaries is presented. It is intended for problems describing the Stokes flow in the situations where one has corners or cusps formed by parts of the domain boundary and, due to the nature of the boundary conditions on these parts of the boundary, these regions have a global effect on the shape of the whole domain and hence have to be resolved with sufficient accuracy. The algorithm combines the boundary integral equation method for the main part of the flow domain and the finite-element method which is used to resolve the corner/cusp regions. Two parts of the solution are matched along a numerical ‘internal interface’ or, as a variant, two interfaces, and they are determined simultaneously by inverting a combined matrix in the course of iterations. The algorithm is illustrated by considering the flow configuration of ‘curtain coating’, a flow where a sheet of liquid impinges onto a moving solid substrate, which is particularly sensitive to what happens in the corner region formed, physically, by the free surface and the solid boundary. The ‘moving contact line problem’ is addressed in the framework of an earlier developed interface formation model which treats the dynamic contact angle as part of the solution, as opposed to it being a prescribed function of the contact line speed, as in the so-called ‘slip models’. Keywords: Dynamic contact angle; finite elements; free surface flows; hybrid numerical technique; Stokes equations.

A well-posed integral equation formulation for three-dimensional rough surface scattering

We consider the problem of scattering of time-harmonic acoustic waves by an unbounded sound-soft rough surface. Recently, a Brakhage Werner type integral equation formulation of this problem has been proposed, based on an ansatz as a combined single- and double-layer potential, but replacing the usual fundamental solution of the Helmholtz equation with an appropriate half-space Green's function. Moreover, it has been shown in the three-dimensional case that this integral equation is uniquely solvable in the space L-2 (Gamma) when the scattering surface G does not differ too much from a plane. In this paper, we show that this integral equation is uniquely solvable with no restriction on the surface elevation or slope. Moreover, we construct explicit bounds on the inverse of the associated boundary integral operator, as a function of the wave number, the parameter coupling the single- and double-layer potentials, and the maximum surface slope. These bounds show that the norm of the inverse operator is bounded uniformly in the wave number, kappa, for kappa > 0, if the coupling parameter h is chosen proportional to the wave number. In the case when G is a plane, we show that the choice eta = kappa/2 is nearly optimal in terms of minimizing the condition number.

Approximations for the generalized temperature integral: a method based on quadrature rules

The generalized temperature integral I(m, x) appears in non-isothermal kinetic analysis when the frequency factor depends on the temperature. A procedure based on Gaussian quadrature to obtain analytical approximations for the integral I(m, x) was proposed. The results showed good agreement between the obtained approximation values and those obtained by numerical integration. Unless other approximations found in literature, the methodology presented in this paper can be easily generalized in order to obtain approximations with the maximum of accurate.

Integral transform method for laminar heat transfer convection of herschel-bulkley fluids within concentric annular ducts

ABSTRACT: Related momentum and energy equations describing the heat and fluid flow of Herschel-Bulkley fluids within concentric annular ducts are analytically solved using the classical integral transform technique, which permits accurate determination of parameters of practical interest in engineering such as friction factors and Nusselt numbers for the duct length. In analyzing the problem, thermally developing flow is assumed and the duct walls are subjected to boundary conditions of first kind. Results are computed for the velocity and temperature fields as well as for the parameters cited above with different power-law indices, yield numbers and aspect ratios. Comparisons are also made with previous work available in the literature, providing direct validation of the results and showing that they are consistent.

Estudo comparativo das técnicas de elementos finitos e equação integral na modelagem eletromagnética bidimensional

Os métodos numéricos de Elementos Finitos e Equação Integral são comumente utilizados para investigações eletromagnéticas na Geofísica, e, para essas modelagens é importante saber qual algoritmo é mais rápido num certo modelo geofísico. Neste trabalho são feitas comparações nos resultados de tempo computacional desses dois métodos em modelos bidimensionais com heterogeneidades condutivas num semiespaço resistivo energizados por uma linha infinita de corrente (com 1000Hz de freqüência) e situada na superfície paralelamente ao "strike" das heterogeneidades. Após a validação e otimização dos programas analisamos o comportamento dos tempos de processamento nos modelos de corpos retangulares variandose o tamanho, o número e a inclinação dos corpos. Além disso, investigamos nesses métodos as etapas que demandam maior custo computacional. Em nossos modelos, o método de Elementos Finitos foi mais vantajoso que o de Equação Integral, com exceção na situação de corpos com baixa condutividade ou com geometria inclinada.