This paper presents an efficient method to calculate the displacement and stress fields in an isotropic elastic half-space having a hemispherical pit and being subject to gravity. The method is semi-analytical and takes advantage of the axisymmetry of the problem. The Boussinesq potentials are used to obtain an analytical solution in series form, which satisfies the equilibrium equations of elastostatics, traction-free boundary conditions on the infinite plane surface and decaying conditions at infinity. The boundary conditions on the free surface of the pit are then imposed numerically, by minimising a quadratic functional of surface elastic energy. The minimisation yields a symmetric and positive definite linear system of equations for the coefficients of the series, whose particular block structure allows its solution in an efficient and robust way. The convergence of the series is verified and the obtained semi-analytical solution is then evaluated, providing numerical results. The method is validated by comparing the semi-analytical solution with the numerical results obtained using a commercial finite element software.

This paper describes the applications of the method of fundamental solutions (MFS) as a mesh-free numerical method for the Stokes' first and second problems which prevail in the semi-infinite domain with constant and oscillatory velocity at the boundary in the fluid-mechanics benchmark problems. The time-dependent fundamental solutions for the semi-infinite problems are used directly to obtain the solution as a linear combination of the unsteady fundamental solution of the diffusion operator. The proposed numerical scheme is free from the conventional Laplace transform or the finite difference scheme to deal with the time derivative term of the governing equation. By properly placing the field points and the source points at a given time level, the solution is advanced in time until steady state solutions are reached. It is not necessary to locate and specify the condition at the infinite domain such as other numerical methods. Since the present method does not need mesh discretization and nodal connectivity, the computational effort and memory storage required are minimal as compared to the domain-oriented numerical schemes. Test results obtained for the Stokes' first and second problems show good comparisons with the analytical solutions. Thus the present numerical scheme has provided a promising mesh-free numerical tool to solve the unsteady semi-infinite problems with the space-time unification for the time-dependent fundamental solution.