4,062 research outputs found
Near-optimal perfectly matched layers for indefinite Helmholtz problems
A new construction of an absorbing boundary condition for indefinite
Helmholtz problems on unbounded domains is presented. This construction is
based on a near-best uniform rational interpolant of the inverse square root
function on the union of a negative and positive real interval, designed with
the help of a classical result by Zolotarev. Using Krein's interpretation of a
Stieltjes continued fraction, this interpolant can be converted into a
three-term finite difference discretization of a perfectly matched layer (PML)
which converges exponentially fast in the number of grid points. The
convergence rate is asymptotically optimal for both propagative and evanescent
wave modes. Several numerical experiments and illustrations are included.Comment: Accepted for publication in SIAM Review. To appear 201
Adaptive Finite Element Method for Simulation of Optical Nano Structures
We discuss realization, properties and performance of the adaptive finite
element approach to the design of nano-photonic components. Central issues are
the construction of vectorial finite elements and the embedding of bounded
components into the unbounded and possibly heterogeneous exterior. We apply the
finite element method to the optimization of the design of a hollow core
photonic crystal fiber. Thereby we look at the convergence of the method and
discuss automatic and adaptive grid refinement and the performance of higher
order elements
Transparent boundary conditions based on the Pole Condition for time-dependent, two-dimensional problems
The pole condition approach for deriving transparent boundary conditions is
extended to the time-dependent, two-dimensional case. Non-physical modes of the
solution are identified by the position of poles of the solution's spatial
Laplace transform in the complex plane. By requiring the Laplace transform to
be analytic on some problem dependent complex half-plane, these modes can be
suppressed. The resulting algorithm computes a finite number of coefficients of
a series expansion of the Laplace transform, thereby providing an approximation
to the exact boundary condition. The resulting error decays super-algebraically
with the number of coefficients, so relatively few additional degrees of
freedom are sufficient to reduce the error to the level of the discretization
error in the interior of the computational domain. The approach shows good
results for the Schr\"odinger and the drift-diffusion equation but, in contrast
to the one-dimensional case, exhibits instabilities for the wave and
Klein-Gordon equation. Numerical examples are shown that demonstrate the good
performance in the former and the instabilities in the latter case
A 2.5D BEM-FEM using a spectral approach to study scattered waves in fluid–solid interaction problems
42nd International Conference on Boundary Elements and other Mesh Reduction Methods, BEM/MRM 2019; ITeCons-University of Coimbra, Coimbra; Portugal; 2 July 2019 through 4 July 2019.
- Publicado en WIT Transactions on Engineering Sciences, Volume 126, 2019, Pages 111-123This work presents a two-and-a-half dimensional (2.5D) spectral formulation based on the finite element method (FEM) and the boundary element method (BEM) to study wave propagation in acoustic and elastic waveguides. The analysis involved superposing two dimensional (2D) problems with different longitudinal wavenumbers. A spectral finite element (SFEM) is proposed to represent waveguides in solids with arbitrary cross-section. Moreover, the BEM is extended to its spectral formulation (SBEM) to study unbounded fluid media and acoustic enclosures. Both approaches use Lagrange polynomials as element shape functions at the Legendre–Gauss–Lobatto (LGL) points. The fluid and solid subdomains are coupled by applying the appropriate boundary conditions at the limiting interface. The proposed method is verified by means of a benchmark problem regarding the scattering of waves by an elastic inclusion. The convergence and the computational effort are evaluated for different h-p strategies. Numerical results show good agreement with the reference solution. Finally, the proposed method is used to study the pressure field generated by an array of elastic fluid-filled scatterers immersed in an acoustic mediumMinisterio de EconomĂa y Competitividad BIA2016-75042-C2-1-
Fourier-Splitting methods for the dynamics of rotating Bose-Einstein condensates
We present a new method to propagate rotating Bose-Einstein condensates
subject to explicitly time-dependent trapping potentials. Using algebraic
techniques, we combine Magnus expansions and splitting methods to yield any
order methods for the multivariate and nonautonomous quadratic part of the
Hamiltonian that can be computed using only Fourier transforms at the cost of
solving a small system of polynomial equations. The resulting scheme solves the
challenging component of the (nonlinear) Hamiltonian and can be combined with
optimized splitting methods to yield efficient algorithms for rotating
Bose-Einstein condensates. The method is particularly efficient for potentials
that can be regarded as perturbed rotating and trapped condensates, e.g., for
small nonlinearities, since it retains the near-integrable structure of the
problem. For large nonlinearities, the method remains highly efficient if
higher order p > 2 is sought. Furthermore, we show how it can adapted to the
presence of dissipation terms. Numerical examples illustrate the performance of
the scheme.Comment: 15 pages, 4 figures, as submitted to journa
Generalized thick strip modelling for vortex-induced vibration of long flexible cylinders
We propose a generalized strip modelling method that is computationally efficient for the VIV prediction of long flexible cylinders in three-dimensional incompressible flow. In order to overcome the shortcomings of conventional strip-theory-based 2D models, the fluid domain is divided into “thick” strips, which are sufficiently thick to locally resolve the small scale turbulence effects and three dimensionality of the flow around the cylinder. An attractive feature of the model is that we independently construct a three-dimensional scale resolving model for individual strips, which have local spanwise scale along the cylinder's axial direction and are only coupled through the structural model of the cylinder. Therefore, this approach is able to cover the full spectrum for fully resolved 3D modelling to 2D strip theory. The connection between these strips is achieved through the calculation of a tensioned beam equation, which is used to represent the dynamics of the flexible body. In the limit, however, a single “thick” strip would fill the full 3D domain. A parallel Fourier spectral/hp element method is employed to solve the 3D flow dynamics in the strip-domain, and then the VIV response prediction is achieved through the strip-structure interactions. Numerical tests on both laminar and turbulent flows as well as the comparison against the fully resolved DNS are presented to demonstrate the applicability of this approach
- …