Simulation of Seawater Intrusion in Coastal Confined Aquifer Using a Point Collocation Method Based Meshfree Model

Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of fresh water from the coastal aquifers. The focus of the present study is on the development of a numerical model based on Meshfree (MFree) method to study the seawater intrusion problem. For the simulation of seawater intrusion problem, widely used models are based on Finite Difference (FDM) and Finite Element (FEM) Methods, which demand well defined grids/meshes and considerable pre-processing efforts. Here, MFree Point Collocation Method (PCM) based on the Radial Basis Function (RBF) is proposed for the simulation. Diffusive interface approach with density-dependent dispersion and solution of flow and solute transport is adopted. These equations are solved using PCM with appropriate boundary conditions. The developed model has been verified with Henry’s problem, and found to be satisfactory. Further the model has been applied to another established problem and an attempt is made to examine the influence of important system parameters including pumping and recharge on the seawater intrusion. The PCM based MFree model is found computationally efficient as preprocessing is avoided when compared to other numerical methods.

Radial point collocation method (RPCM) for solving convection-diffusion problems

In this paper, Radial point collocation method （RPCM）, a kind of meshfree method, is applied to solve convectiondiffusion problem. The main feature of this approach is to use the interpolation schemes in local supported domains based on radial basis functions. As a result, this method is local and hence the system matrix is banded which is very attractive for practical engineering problems. In the numerical examination, RPCM is applied to solve non-linear convection-diffusion 2D Burgers equations. The results obtained by RPCM demonstrate the accuracy and efficiency of the proposed method for solving transient fluid dynamic problems. A fictitious point scheme is adopted to improve the solution accuracy while Neumann boundary conditions exist. The meshfree feature of the nresent method is verv attractive in solving comnutational fluid nroblems.

AHermitian C^(2) Differential Reproducing Kernel Interpolation Meshless Method for the 3D Microstructure-Dependent Static Flexural Analysis of Simply Supported and Functionally Graded Microplates

This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.

NONLINEAR DYNAMIC RESPONSE AND DYNAMIC BUCKLING OF SHALLOW SPHERICAL SHELLS WITH CIRCULAR HOLE

In this paper, the nonlinear equations of motion for shallow spherical shells with axisymmetric deformation including transverse shear are derived. The nonlinear static and dynamic response and dynamic buckling of shallow spherical shells with circular hole on elastically restrained edge are investigated. By using the orthogonal point collocation method for space and Newmarh-β scheme for time, the displacement functions are separated and the nonlinear differential equations are replaced by linear algebraic equations to seek solutions. The numerical results are presented for different cases and compared with available data.

Axisymmetrical Nonlinear Bending and Buckling of a Circular Plate Under Large Load

With the terms of the exact series solution taken as trial functions, the method of point collocation was used to calculate the large deflection of a circular plate. The axisymmetrical bending formulae were developed for the calculation of a circular plate subjected to polynomial distributed loads, a concentrated load at the center, uniform radial forces and moments along the edge or their combinations. The support may be elastic. The buckling load was calculated. Under action of uniformly distributed load, central load or their compound load, solutions were compared with those obtained by other methods. Buckling beyond critical thrust was compared with that calculated by the power series method. The method presented in this paper has advantages of wide convergent range, high precision and short computing time. Moreover, the computing time is nearly independent of the complexity of the loads.