On atomistic-to-continuum couplings without ghost forces in three dimensions

In this paper we construct energy based numerical methods free of ghost forces in three dimen- sional lattices arising in crystalline materials. The analysis hinges on establishing a connection of the coupled system to conforming finite elements. Key ingredients are: (i) a new representation of discrete derivatives related to long range interactions of atoms as volume integrals of gradients of piecewise linear functions over bond volumes, and (ii) the construction of an underlying globally continuous function representing the coupled modeling method

Benchmarking two simulation models for underwater and atmospheric sound propagation

Computational wave propagation models are widely used in underwater and atmospheric sound propagation simulation. In most realistic cases the physical domains involved are irregular. We have developed finite element techniques, applied to general irregular meshes and coupled with discrete, artificial absorbing boundary conditions of nonlocal type, for the Helmholtz equation and its 'standard' parabolic approximation. The physical domain is axially symmetric, with several fluid layers of variable acoustic properties. Boundaries and interfaces of general topography are allowed. The resulting models are referred to as the FENL and CNP1-NL models, respectively. We present results of the FENL model for underwater acoustic applications related to object identification and of the CNP1-NL for atmospheric sound propagation over an irregular terrain. © 2005 Elsevier Ltd. All rights reserved

Coupled mode and finite element approximations of underwater sound propagation problems in general stratified environments

We compare the results of a coupled mode method with those of a finite element method and also of COUPLE on two test problems of sound propagation and scattering in cylindrically symmetric, underwater, multilayered acoustic waveguides with range-dependent interface topographies. We observe, in general, very good agreement between the results of the three codes. In some cases in which the frequency of the harmonic point source is such that an eigenvalue of the local vertical problem remains small in magnitude and changes sign several times in the vicinity of the interface nonhomogeneity, the discrepancies between the results of the three codes increase, but remain small in absolute terms. © 2008 IMACS