scholarly journals Evolution of a semidiscrete system modeling the scattering of acoustic waves by a piezoelectric solid

2018 ◽  
Vol 52 (2) ◽  
pp. 423-455 ◽  
Author(s):  
Thomas S. Brown ◽  
Tonatiuh Sánchez-Vizuet ◽  
Francisco-Javier Sayas

We consider a model problem of the scattering of linear acoustic waves in free homogeneous space by an elastic solid. The stress tensor in the solid combines the effect of a linear dependence of strains with the influence of an existing electric field. The system is closed using Gauss’s law for the associated electric displacement. Well-posedness of the system is studied by its reformulation as a first order in space and time differential system with help of an elliptic lifting operator. We then proceed to studying a semidiscrete formulation, corresponding to an abstract Finite Element discretization in the electric and elastic fields, combined with an abstract Boundary Element approximation of a retarded potential representation of the acoustic field. The results obtained with this approach improve estimates obtained with Laplace domain techniques. While numerical experiments illustrating convergence of a fully discrete version of this problem had already been published, we demonstrate some properties of the full model with some simulations for the two dimensional case.

Author(s):  
R. Becker ◽  
R. Koch ◽  
M. F. Modest ◽  
H.-J. Bauer

The present article introduces a new method to solve the radiative transfer equation (RTE). First, a finite element discretization of the solid angle dependence is derived, wherein the coefficients of the finite element approximation are functions of the spatial coordinates. The angular basis functions are defined according to finite element principles on subdivisions of the octahedron. In a second step, these spatially dependent coefficients are discretized by spatial finite elements. This approach is very attractive, since it provides a concise derivation for approximations of the angular dependence with an arbitrary number of angular nodes. In addition, the usage of high-order angular basis functions is straightforward. In the current paper the governing equations are first derived independently of the actual angular approximation. Then, the design principles for the angular mesh are discussed and the parameterization of the piecewise angular basis functions is derived. In the following, the method is applied to two-dimensional test cases which are commonly used for the validation of approximation methods of the RTE. The results reveal that the proposed method is a promising alternative to the well-established practices like the Discrete Ordinates Method (DOM) and provides highly accurate approximations. A test case known to exhibit the ray effect in the DOM verifies the ability of the new method to avoid ray effects.


2015 ◽  
Vol 8 (4) ◽  
pp. 582-604
Author(s):  
Zhengqin Yu ◽  
Xiaoping Xie

AbstractThis paper proposes and analyzes semi-discrete and fully discrete hybrid stress finite element methods for elastodynamic problems. A hybrid stress quadrilateral finite element approximation is used in the space directions. A second-order center difference is adopted in the time direction for the fully discrete scheme. Error estimates of the two schemes, as well as a stability result for the fully discrete scheme, are derived. Numerical experiments are done to verify the theoretical results.


2016 ◽  
Vol 309 ◽  
pp. 106-130 ◽  
Author(s):  
Matthew E. Hassell ◽  
Francisco-Javier Sayas

Author(s):  
Wenjun Cai ◽  
Huai Zhang ◽  
Yushun Wang

We propose a novel stable and efficient dissipation-preserving method for acoustic wave propagations in attenuating media with both correct phase and amplitude. Through introducing the conformal multi-symplectic structure, the intrinsic dissipation law and the conformal symplectic conservation law are revealed for the damped acoustic wave equation. The proposed algorithm is exactly designed to preserve a discrete version of the conformal symplectic conservation law. More specifically, two subsystems in conjunction with the original damped wave equation are derived. One is actually the conservative Hamiltonian wave equation and the other is a dissipative linear ordinary differential equation (ODE) system. Standard symplectic method is devoted to the conservative system, whereas the analytical solution is obtained for the ODE system. An explicit conformal symplectic scheme is constructed by concatenating these two parts of solutions by the Strang splitting technique. Stability analysis and convergence tests are given thereafter. A benchmark model in homogeneous media is presented to demonstrate the effectiveness and advantage of our method in suppressing numerical dispersion and preserving the energy dissipation. Further numerical tests show that our proposed method can efficiently capture the dissipation in heterogeneous media.


Author(s):  
Lonny L. Thompson

Abstract A Computational Structural Acoustics (CSA) capability for solving scattering, radiation, and other problems related to the acoustics of submerged structures has been developed by employing some of the recent algorithmic trends in Computational Fluid Dynamics (CFD), namely time-discontinuous Galerkin Least-Squares finite element methods. Traditional computational methods toward simulation of acoustic radiation and scattering from submerged elastic bodies have been primarily based on frequency domain formulations. These classical time-harmonic approaches (including boundary element, finite element, and finite difference methods) have been successful for problems involving a limited range of frequencies (narrow band response) and scales (wavelengths) that are large compared to the characteristic dimensions of the elastic structure. Attempts at solving large-scale structural acoustic systems with dimensions that are much larger than the operating wavelengths and which are complex, consisting of many different components with different scales and broadband frequencies, has revealed limitations of many of the classical methods. As a result, there has been renewed interest in new innovative approaches, including time-domain approaches. This paper describes recent advances in the development of a new class of high-order accurate and unconditionally stable space-time methods for structural acoustics which employ finite element discretization of the time domain as well as the usual discretization of the spatial domain. The formulation is based on a space-time variational equation for both the acoustic fluid and elastic structure together with their interaction. Topics to be discussed include the development and implementation of higher-order accurate non-reflecting boundary conditions based on the exact impedance relation through the. Dirichlet-to-Neumann (DtN) map, and a multi-field representation for the acoustic fluid based on independent pressure and velocity potential variables. Numerical examples involving radiation and scattering of acoustic waves are presented to illustrate the high-order accuracy achieved by the new methodology for CSA.


2016 ◽  
Vol 6 (4) ◽  
pp. 434-447 ◽  
Author(s):  
M. Mbehou ◽  
R. Maritz ◽  
P.M.D. Tchepmo

AbstractThis article is devoted to the study of the finite element approximation for a nonlocal nonlinear parabolic problem. Using a linearised Crank-Nicolson Galerkin finite element method for a nonlinear reaction-diffusion equation, we establish the convergence and error bound for the fully discrete scheme. Moreover, important results on exponential decay and vanishing of the solutions in finite time are presented. Finally, some numerical simulations are presented to illustrate our theoretical analysis.


2017 ◽  
Vol 27 (13) ◽  
pp. 2511-2556 ◽  
Author(s):  
Filip Rindler ◽  
Sebastian Schwarzacher ◽  
Endre Süli

Rate-independent systems arise in a number of applications. Usually, weak solutions to such problems with potentially very low regularity are considered, requiring mathematical techniques capable of handling nonsmooth functions. In this work, we prove the existence of Hölder-regular strong solutions for a class of rate-independent systems. We also establish additional higher regularity results that guarantee the uniqueness of strong solutions. The proof proceeds via a time-discrete Rothe approximation and careful elliptic regularity estimates depending in a quantitative way on the (local) convexity of the potential featuring in the model. In the second part of the paper, we show that our strong solutions may be approximated by a fully discrete numerical scheme based on a spatial finite element discretization, whose rate of convergence is consistent with the regularity of strong solutions whose existence and uniqueness are established.


2008 ◽  
Vol 8 (3) ◽  
pp. 294-308 ◽  
Author(s):  
A. PEDAS ◽  
E. TAMME

Abstract Approximations to a solution and its derivatives of a boundary value problem of an nth order linear Fredholm integro-differential equation with weakly sin-gular or other nonsmooth kernels have been determined. These approximations are piecewise polynomial functions on special graded grids. To find them, a fully discrete version of the Galerkin method has been constructed. This version is based on a dis-crete inner product concept and some suitable product integration techniques. Optimal global convergence estimates have been derived and a collection of numerical results of a test problem is given.


Sign in / Sign up

Export Citation Format

Share Document