Simulating free surface problems using Discrete Least Squares Meshless method

2010 ◽  
Vol 39 (3) ◽  
pp. 461-470 ◽  
Author(s):  
G. Shobeyri ◽  
M.H. Afshar
2014 ◽  
Vol 2014 ◽  
pp. 1-16 ◽  
Author(s):  
J. Amani ◽  
A. Saboor Bagherzadeh ◽  
T. Rabczuk

The node moving and multistage node enrichment adaptive refinement procedures are extended in mixed discrete least squares meshless (MDLSM) method for efficient analysis of elasticity problems. In the formulation of MDLSM method, mixed formulation is accepted to avoid second-order differentiation of shape functions and to obtain displacements and stresses simultaneously. In the refinement procedures, a robust error estimator based on the value of the least square residuals functional of the governing differential equations and its boundaries at nodal points is used which is inherently available from the MDLSM formulation and can efficiently identify the zones with higher numerical errors. The results are compared with the refinement procedures in the irreducible formulation of discrete least squares meshless (DLSM) method and show the accuracy and efficiency of the proposed procedures. Also, the comparison of the error norms and convergence rate show the fidelity of the proposed adaptive refinement procedures in the MDLSM method.


2016 ◽  
Vol 24 (11) ◽  
pp. 2149-2164 ◽  
Author(s):  
Majid Darehmiraki ◽  
Mohammad Hadi Farahi ◽  
Sohrab Effati

We use a hybrid local meshless method to solve the distributed optimal control problem of a system governed by parabolic partial differential equations with Caputo fractional time derivatives of order α ∈ (0, 1]. The presented meshless method is based on the linear combination of moving least squares and radial basis functions in the same compact support, this method will change between interpolation and approximation. The aim of this paper is to solve the system of coupled fractional partial differential equations, with necessary and sufficient conditions, for fractional distributed optimal control problems using a combination of moving least squares and radial basis functions. To keep matters simple, the problem has been considered in the one-dimensional case, however the techniques can be employed for both the two- and three-dimensional cases. Several test problems are employed and results of numerical experiments are presented. The obtained results confirm the acceptable accuracy of the proposed method.


Geofluids ◽  
2018 ◽  
Vol 2018 ◽  
pp. 1-13 ◽  
Author(s):  
Jing-En Xiao ◽  
Cheng-Yu Ku ◽  
Chih-Yu Liu ◽  
Wei-Chung Yeih

A novel boundary-type meshless method for modeling geofluid flow in heterogeneous geological media was developed. The numerical solutions of geofluid flow are approximated by a set of particular solutions of the subsurface flow equation which are expressed in terms of sources located outside the domain of the problem. This pioneering study is based on the collocation Trefftz method and provides a promising solution which integrates the T-Trefftz method and F-Trefftz method. To deal with the subsurface flow problems of heterogeneous geological media, the domain decomposition method was adopted so that flux conservation and the continuity of pressure potential at the interface between two consecutive layers can be considered in the numerical model. The validity of the model is established for a number of test problems. Application examples of subsurface flow problems with free surface in homogenous and layered heterogeneous geological media were also carried out. Numerical results demonstrate that the proposed method is highly accurate and computationally efficient. The results also reveal that it has great numerical stability for solving subsurface flow with nonlinear free surface in layered heterogeneous geological media even with large contrasts in the hydraulic conductivity.


Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. S83-S94 ◽  
Author(s):  
Yikang Zheng ◽  
Yibo Wang ◽  
Xu Chang

Free-surface-related multiples can provide extra illumination of the subsurface and thus can be usefully included in migration procedures. However, most multiple migration approaches require separation of primaries and free-surface-related multiples or at least prediction of multiples in advance, which is time consuming and prone to errors. The data-to-data migration (DDM) method migrates free-surface-related multiples by forward and backward propagating the recorded full data (containing primaries and free-surface-related multiples). For DDM, there is no need to predict or separate multiples, but the migration results suffer from the crosstalk generated by crosscorrelations of undesired seismic events, e.g., primaries and second-order free-surface-related multiples. We have developed least-squares DDM (LSDDM) for marine data to eliminate the crosstalk generated by DDM. In each iteration, the forward-propagated primaries and free-surface-related multiples are crosscorrelated with the backward-propagated primary and free-surface-related multiple residuals to form the reflectivity gradient. We use a three-layer model and the Marmousi model for numerical tests. The results validate that LSDDM can provide a migrated image with higher signal-to-noise ratio and more balanced amplitudes than DDM. The LSDDM approach might be valuable for general subsurface imaging for marine seismic data when the migration velocity is accurate, and the acquired data have sufficient recording time.


Sign in / Sign up

Export Citation Format

Share Document