Baran Aydın^1,2, Utku Kânoğlu³

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.2, pp. 149-156, 2007, DOI:10.3970/cmes.2007.021.149

Abstract We developed analytical solutions to the wind set-down and the wind set-down relaxation problems. The response of the ocean to the wind blowing over a long-narrow and linearly sloping shallow basin is referred to as wind set-down. The shoreline exhibits oscillatory behavior when the wind calms down and the resulting problem is referred to as wind set-down relaxation. We use an existing hodograph-type transformation that was introduced to solve the nonlinear shallow-water wave equations analytically for long wave propagation and obtain an explicit-transform analytical solution for wind set-down. For the wind set-down relaxation, the nonlinear shallow-water wave equations are solved… More >

H.B. Chen¹, D.J. Fu¹, P.Q. Zhang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.20, No.2, pp. 85-98, 2007, DOI:10.3970/cmes.2007.020.085

Abstract Researches today show that, both approximation and dispersion errors are encountered by classical Galerkin FEM solutions for Helmholtz equation governing the harmonic wave propagation, which leads to numerical inaccuracies especially for high wave number cases. In this paper, Local Boundary Integral Equation Method (LBIEM) is firstly implemented to solve the boundary value problem of Helmholtz equation. Then the regularized LBIE is proposed to overcome the singularities of the boundary integrals in the LBIEM. Owing to the advantages of the Moving Least Square Approximation (MLSA), the frequency-dependent basis functions modified by the harmonic wave propagation solutions are easily adopted instead of… More >

Q.W. Ma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.3, pp. 223-234, 2007, DOI:10.3970/cmes.2007.018.223

Abstract Two methods have been recently developed by the author and his group: one called MLPG_R (Meshless Local Petrov-Galerkin method based on Rankine source solution) and the other called QALE-FEM (Quasi Arbitrary Lagrangian-Eulerian Finite Element Method). The former is a meshless method developed from a general MLPG (Meshless Local Petrov-Galerkin) method and is more computationally efficient than the general one when applied to modelling nonlinear water waves. The later is a mesh-based method similar to a conventional finite element method (FEM) when discretizing the governing equations but different from the conventional one in managing the mesh. In this paper, they are… More >

A.N. Guz¹, O.V. Menshykov^1,2, V.V. Zozulya³, I.A. Guz²

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.3, pp. 205-214, 2007, DOI:10.3970/cmes.2007.017.205

Abstract The contact interaction of opposite faces of an elliptical crack is studied for the case of a normal time-harmonic shear wave loading. The distribution of stress intensity factors (shear modes II and III) as functions of the wave number and the friction coefficient is investigated. The results are compared with those obtained for an elliptical crack without allowance for the contact interaction. More >

Ivan Christov¹, C.I. Christov², P.M. Jordan³

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.1, pp. 47-54, 2007, DOI:10.3970/cmes.2007.017.047

Abstract Two widely-used weakly-nonlinear models of acoustic wave propagation --- the inviscid Kuznetsov equation (IKE) and the Lighthill--Westervelt equation (LWE) --- are investigated numerically using a Godunov-type finite-difference scheme. A reformulation of the models as conservation laws is proposed, making it possible to use the numerical tools developed for the Euler equations to study the IKE and LWE, even after the time of shock-formation. It is shown that while the IKE is, without qualification, in very good agreement with the Euler equations, even near the time of shock formation, the same cannot generally be said for the LWE. More >

D. Soares Jr.^1,2, W.J. Mansur¹, D.L. Lima³

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.1, pp. 19-34, 2007, DOI:10.3970/cmes.2007.017.019

Abstract The present paper discussion is concerned with the development of robust and efficient algorithms to model propagation of interacting acoustic and elastic waves. The paper considers acoustic-elastic, acoustic-acoustic and elastic-elastic partitioned analyses of coupled systems; however, the focus here is the acoustic-elastic coupling considering finite elements and the acoustic-acoustic coupling considering finite elements and finite differences (other coupling procedures can be implemented analogously). One important feature of the algorithms presented is that they allow considering different time-steps for different sub-domains; so it is possible to substantially improve efficiency, accuracy and stability of the central difference time integration algorithm employed here.… More >

Q.W. Ma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.2, pp. 193-210, 2005, DOI:10.3970/cmes.2005.009.193

Abstract Recently, the MLPG (Meshless Local Petrov-Galerkin Method) method has been successfully extended to simulating nonlinear water waves [Ma, (2005)]. In that paper, the author employed the Heaviside step function as the test function to formulate the weak form over local sub-domains, acquiring an expression in terms of pressure gradient. In this paper, the solution for Rankine sources is taken as the test function and the local weak form is expressed in term of pressure rather than pressure gradient. Apart from not including pressure gradient, velocity gradient is also eliminated from the weak form. In addition, a semi-analytical technique is developed… More >

J. Bielak¹, O. Ghattas², E.-J. Kim³

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.2, pp. 99-112, 2005, DOI:10.3970/cmes.2005.010.099

Abstract We present a parallel octree-based finite element method for large-scale earthquake ground motion simulation in realistic basins. The octree representation combines the low memory per node and good cache performance of finite difference methods with the spatial adaptivity to local seismic wavelengths characteristic of unstructured finite element methods. Several tests are provided to verify the numerical performance of the method against Green's function solutions for homogeneous and piecewise homogeneous media, both with and without anelastic attenuation. A comparison is also provided against a finite difference code and an unstructured tetrahedral finite element code for a simulation of the 1994 Northridge… More >

S.G. Bardenhagen¹, J.E. Guilkey², K.M. Roessig³, J.U. Brackbill⁴, W.M. Witzel⁵, J.C.Foster⁶

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.4, pp. 509-522, 2001, DOI:10.3970/cmes.2001.002.509

Abstract Contact between deformable bodies is a difficult problem in the analysis of engineering systems. A new approach to contact has been implemented using the Material Point Method for solid mechanics, Bardenhagen, Brackbill, and Sulsky (2000a). Here two improvements to the algorithm are described. The first is to include the normal traction in the contact logic to more appropriately determine the free separation criterion. The second is to provide numerical stability by scaling the contact impulse when computational grid information is suspect, a condition which can be expected to occur occasionally as material bodies move through the computational grid. The modifications… More >

S.-P. Zhu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 227-236, 2001, DOI:10.3970/cmes.2001.002.227

Abstract Problems defined on infinite domains must be treated on a finite computational domain. The treatment of the artificially placed boundaries (usually referred to as open boundaries) of such domain truncations can be quite subtle; an over truncation would normally result in large, undesirable reflection of signals back to the computational domain whereas an under truncation would imply an injudicious use of computational resources. In particular, problems occur when strongly nonlinear free surface waves generated in a numerical wave tank are passing through such an open boundary.
In this paper, some recent numerical test results of an innovative treatment of… More >