Modeling elastic wave propagation in fluid-filled boreholes drilled in nonhomogeneous media: BEM – MLPG versus BEM-FEM coupling

The efficiency of two coupling formulations, the boundary element method (BEM)-meshless local Petrov–Galerkin (MLPG) versus the BEM-finite element method (FEM), used to simulate the elastic wave propagation in fluid-filled boreholes generated by a blast load, is compared. The longitudinal geometry is assumed to be invariant in the axial direction (2.5D formulation). The material properties in the vicinity of the borehole are assumed to be nonhomogeneous as a result of the construction process and the ageing of the material. In both models, the BEM is used to tackle the propagation within the fluid domain inside the borehole and the unbounded homogeneous domain. The MLPG and the FEM are used to simulate the confined, damaged, nonhomogeneous, surrounding borehole, thus utilizing the advantages of these methods in modeling nonhomogeneous bounded media. In both numerical techniques the coupling is accomplished directly at the nodal points located at the common interfaces. Continuity of stresses and displacements is imposed at the solid–solid interface, while continuity of normal stresses and displacements and null shear stress are prescribed at the fluid–solid interface. The performance of each coupled BEM-MLPG and BEM-FEM approach is determined using referenced results provided by an analytical solution developed for a circular multi-layered subdomain. The comparison of the coupled techniques is evaluated for different excitation frequencies, axial wavenumbers and degrees of freedom (nodal points).Ministerio de Economía y Competitividad BIA2013-43085-PCentro Informático Científico de Andalucía (CICA

Fluctuating surface-current formulation of radiative heat transfer: theory and applications

We describe a novel fluctuating-surface current formulation of radiative heat transfer between bodies of arbitrary shape that exploits efficient and sophisticated techniques from the surface-integral-equation formulation of classical electromagnetic scattering. Unlike previous approaches to non-equilibrium fluctuations that involve scattering matrices---relating "incoming" and "outgoing" waves from each body---our approach is formulated in terms of "unknown" surface currents, laying at the surfaces of the bodies, that need not satisfy any wave equation. We show that our formulation can be applied as a spectral method to obtain fast-converging semi-analytical formulas in high-symmetry geometries using specialized spectral bases that conform to the surfaces of the bodies (e.g. Fourier series for planar bodies or spherical harmonics for spherical bodies), and can also be employed as a numerical method by exploiting the generality of surface meshes/grids to obtain results in more complicated geometries (e.g. interleaved bodies as well as bodies with sharp corners). In particular, our formalism allows direct application of the boundary-element method, a robust and powerful numerical implementation of the surface-integral formulation of classical electromagnetism, which we use to obtain results in new geometries, including the heat transfer between finite slabs, cylinders, and cones

Numerical methods for computing Casimir interactions

We review several different approaches for computing Casimir forces and related fluctuation-induced interactions between bodies of arbitrary shapes and materials. The relationships between this problem and well known computational techniques from classical electromagnetism are emphasized. We also review the basic principles of standard computational methods, categorizing them according to three criteria---choice of problem, basis, and solution technique---that can be used to classify proposals for the Casimir problem as well. In this way, mature classical methods can be exploited to model Casimir physics, with a few important modifications.Comment: 46 pages, 142 references, 5 figures. To appear in upcoming Lecture Notes in Physics book on Casimir Physic

A Space-time Smooth Artificial Viscosity Method For Nonlinear Conservation Laws

We introduce a new methodology for adding localized, space-time smooth, artificial viscosity to nonlinear systems of conservation laws which propagate shock waves, rarefactions, and contact discontinuities, which we call the $C$-method. We shall focus our attention on the compressible Euler equations in one space dimension. The novel feature of our approach involves the coupling of a linear scalar reaction-diffusion equation to our system of conservation laws, whose solution $C(x,t)$ is the coefficient to an additional (and artificial) term added to the flux, which determines the location, localization, and strength of the artificial viscosity. Near shock discontinuities, $C(x,t)$ is large and localized, and transitions smoothly in space-time to zero away from discontinuities. Our approach is a provably convergent, spacetime-regularized variant of the original idea of Richtmeyer and Von Neumann, and is provided at the level of the PDE, thus allowing a host of numerical discretization schemes to be employed. We demonstrate the effectiveness of the $C$-method with three different numerical implementations and apply these to a collection of classical problems: the Sod shock-tube, the Osher-Shu shock-tube, the Woodward-Colella blast wave and the Leblanc shock-tube. First, we use a classical continuous finite-element implementation using second-order discretization in both space and time, FEM-C. Second, we use a simplified WENO scheme within our $C$-method framework, WENO-C. Third, we use WENO with the Lax-Friedrichs flux together with the $C$-equation, and call this WENO-LF-C. All three schemes yield higher-order discretization strategies, which provide sharp shock resolution with minimal overshoot and noise, and compare well with higher-order WENO schemes that employ approximate Riemann solvers, outperforming them for the difficult Leblanc shock tube experiment.Comment: 34 pages, 27 figure

Thermoacoustic instability - a dynamical system and time domain analysis

This study focuses on the Rijke tube problem, which includes features relevant to the modeling of thermoacoustic coupling in reactive flows: a compact acoustic source, an empirical model for the heat source, and nonlinearities. This thermo-acoustic system features a complex dynamical behavior. In order to synthesize accurate time-series, we tackle this problem from a numerical point-of-view, and start by proposing a dedicated solver designed for dealing with the underlying stiffness, in particular, the retarded time and the discontinuity at the location of the heat source. Stability analysis is performed on the limit of low-amplitude disturbances by means of the projection method proposed by Jarlebring (2008), which alleviates the linearization with respect to the retarded time. The results are then compared to the analytical solution of the undamped system, and to Galerkin projection methods commonly used in this setting. This analysis provides insight into the consequences of the various assumptions and simplifications that justify the use of Galerkin expansions based on the eigenmodes of the unheated resonator. We illustrate that due to the presence of a discontinuity in the spatial domain, the eigenmodes in the heated case, predicted by using Galerkin expansion, show spurious oscillations resulting from the Gibbs phenomenon. By comparing the modes of the linear to that of the nonlinear regime, we are able to illustrate the mean-flow modulation and frequency switching. Finally, time-series in the fully nonlinear regime, where a limit cycle is established, are analyzed and dominant modes are extracted. The analysis of the saturated limit cycles shows the presence of higher frequency modes, which are linearly stable but become significant through nonlinear growth of the signal. This bimodal effect is not captured when the coupling between different frequencies is not accounted for.Comment: Submitted to Journal of Fluid Mechanic

Stationary bumps in a piecewise smooth neural field model with synaptic depression

We analyze the existence and stability of stationary pulses or bumps in a one–dimensional piecewise smooth neural field model with synaptic depression. The continuum dynamics is described in terms of a nonlocal integrodifferential equation, in which the integral kernel represents the spatial distribution of synaptic weights between populations of neurons whose mean firing rate is taken to be a Heaviside function of local activity. Synaptic depression dynamically reduces the strength of synaptic weights in response to increases in activity. We show that in the case of a Mexican hat weight distribution, there exists a stable bump for sufficiently weak synaptic depression. However, as synaptic depression becomes stronger, the bump became unstable with respect to perturbations that shift the boundary of the bump, leading to the formation of a traveling pulse. The local stability of a bump is determined by the spectrum of a piecewise linear operator that keeps track of the sign of perturbations of the bump boundary. This results in a number of differences from previous studies of neural field models with Heaviside firing rate functions, where any discontinuities appear inside convolutions so that the resulting dynamical system is smooth. We also extend our results to the case of radially symmetric bumps in two–dimensional neural field models