Frequency-explicit stability estimates for time-harmonic elastodynamic problems in nearly incompressible materials

We consider time-harmonic elastodynamic problems in heterogeneous media.cWe focus on scattering problems in the high-frequency regime and incnearly incompressible media, where the the angular frequency $\omega$ and ratio of the Lam\'e parameters $\lambda/\mu$ may both be large. We derive stability estimates controlling the norm of the solution by the norm of the right-hand side up to a fully-explicit constant. Crucially, under natural assumptions on the domain and coefficients, this constant increases linearly with $\omega$ and is uniform in the ratio $\lambda/\mu$

Finite element approximation of electromagnetic fields using nonfitting meshes for Geophysics

We analyze the use of nonfitting meshes for simulating the propagation of electromagnetic waves inside the earth with applications to borehole logging. We avoid the use of parameter homogenization and employ standard edge finite element basis functions. For our geophysical applications, we consider a 3D Maxwell’s system with piecewise constant conductivity and globally constant permittivity and permeability. The model is analyzed and discretized using both the Eand H-formulations. Our main contribution is to develop a sharp error estimate for both the electric and magnetic fields. In the presence of singularities, our estimate shows that the magnetic field approximation is converging faster than the electric field approximation. As a result, we conclude that error estimates available in the literature are sharp with respect to the electric field error but provide pessimistic convergence rates for the magnetic field in our geophysical applications. Another surprising consequence of our analysis is that nonfitting meshes deliver the same convergence rate as fitting meshes to approximate the magnetic field. Our theoretical results are numerically illustrated via 2D experiments. For the analyzed cases, the accuracy loss due to the use of nonfitting meshes islimited, even for high conductivity contrasts

A postprocessing technique for a discontinuous Galerkin discretization of time-dependent Maxwell's equations

We present a novel postprocessing technique for a discontinuous Galerkin (DG) discretization of time-dependent Maxwell's equations that we couple with an explicit Runge-Kutta time-marching scheme. The postprocessed electromagnetic field converges one order faster than the unprocessed solution in the H(curl)-norm. The proposed approach is local, in the sense that the enhanced solution is computed independently in each cell of the computational mesh, and at each time step of interest. As a result, it is inexpensive to compute, especially if the region of interest is localized, either in time or space. The key ideas behind this postprocessing technique stem from hybridizable discontinuous Galerkin (HDG) methods, which are equivalent to the analyzed DG scheme for specific choices of penalization parameters. We present several numerical experiments that highlight the superconvergence properties of the postprocessed electromagnetic field approximation

An ensemble approach to assess hydrological models’ contribution to uncertainties in the analysis of climate change impact on water resources

Over the recent years, several research efforts investigated the impact of climate change on water resources for different regions of the world. The projection of future river flows is affected by different sources of uncertainty in the hydro-climatic modelling chain. One of the aims of the QBic3 5 project (Que´bec-Bavarian International Collaboration on Climate Change) is to assess the contribution to uncertainty of hydrological models by using an ensemble of hydrological models presenting a diversity of structural complexity (i.e. lumped, semi distributed and distributed models). The study investigates two humid, mid-latitude catchments with natural flow conditions; one located in 10 Southern Que´bec (Canada) and one in Southern Bavaria (Germany). Daily flow is simulated with four different hydrological models, forced by outputs from regional climate models driven by a given number of GCMs’ members over a reference (1971–2000) and a future (2041–2070) periods. The results show that the choice of the hydrological model does strongly affect the climate change response of selected hydrological indicators, especially those related to low flows. Indicators related to high flows seem less sensitive on the choice of the hydrological model. Therefore, the computationally less demanding models (usually simple, lumped and conceptual) give a significant level of trust for high and overall mean flows

On the need for bias correction in regional climate scenarios to assess climate change impacts on river runoff

In climate change impact research, the assessment of future river runoff as well as the catchment scale water balance is impeded by different sources of modeling uncertainty. Some research has already been done in order to quantify the uncertainty of climate 5 projections originating from the climate models and the downscaling techniques as well as from the internal variability evaluated from climate model member ensembles. Yet, the use of hydrological models adds another layer of incertitude. Within the QBic3 project (Qu´ebec-Bavaria International Collaboration on Climate Change) the relative contributions to the overall uncertainty from the whole model chain (from global climate 10 models to water management models) are investigated using an ensemble of multiple climate and hydrological models. Although there are many options to downscale global climate projections to the regional scale, recent impact studies tend to use Regional Climate Models (RCMs). One reason for that is that the physical coherence between atmospheric and land-surface 15 variables is preserved. The coherence between temperature and precipitation is of particular interest in hydrology. However, the regional climate model outputs often are biased compared to the observed climatology of a given region. Therefore, biases in those outputs are often corrected to reproduce historic runoff conditions from hydrological models using them, even if those corrections alter the relationship between temperature and precipitation. So, as bias correction may affect the consistency between RCM output variables, the use of correction techniques and even the use of (biased) climate model data itself is sometimes disputed among scientists. For those reasons, the effect of bias correction on simulated runoff regimes and the relative change in selected runoff indicators is explored. If it affects the conclusion of climate change analysis in 25 hydrology, we should consider it as a source of uncertainty. If not, the application of bias correction methods is either unnecessary in hydro-climatic projections, or safe to use as it does not alter the change signal of river runoff. The results of the present paper highlight the analysis of daily runoff simulated with four different hydrological models in two natural-flow catchments, driven by different regional climate models for a reference and a future period. As expected, bias correction of climate model outputs is important for the reproduction of the runoff regime of the 5 past regardless of the hydrological model used. Then again, its impact on the relative change of flow indicators between reference and future period is weak for most indicators with the exception of the timing of the spring flood peak. Still, our results indicate that the impact of bias correction on runoff indicators increases with bias in the climate simulations

A controllability method for Maxwell's equations

We propose a controllability method for the numerical solution of time-harmonic Maxwell's equations in their first-order formulation. By minimizing a quadratic cost functional, which measures the deviation from periodicity, the controllability method determines iteratively a periodic solution in the time domain. At each conjugate gradient iteration, the gradient of the cost functional is simply computed by running any time-dependent simulation code forward and backward for one period, thus leading to a non-intrusive implementation easily integrated into existing software. Moreover, the proposed algorithm automatically inherits the parallelism, scalability, and low memory footprint of the underlying time-domain solver. Since the time-periodic solution obtained by minimization is not necessarily unique, we apply a cheap post-processing filtering procedure which recovers the time-harmonic solution from any minimizer. Finally, we present a series of numerical examples which show that our algorithm greatly speeds up the convergence towards the desired time-harmonic solution when compared to simply running the time-marching code until the time-harmonic regime is eventually reached

Area distribution and the average shape of a L\'evy bridge

We consider a one dimensional L\'evy bridge x_B of length n and index 0 < \alpha < 2, i.e. a L\'evy random walk constrained to start and end at the origin after n time steps, x_B(0) = x_B(n)=0. We compute the distribution P_B(A,n) of the area A = \sum_{m=1}^n x_B(m) under such a L\'evy bridge and show that, for large n, it has the scaling form P_B(A,n) \sim n^{-1-1/\alpha} F_\alpha(A/n^{1+1/\alpha}), with the asymptotic behavior F_\alpha(Y) \sim Y^{-2(1+\alpha)} for large Y. For \alpha=1, we obtain an explicit expression of F_1(Y) in terms of elementary functions. We also compute the average profile < \tilde x_B (m) > at time m of a L\'evy bridge with fixed area A. For large n and large m and A, one finds the scaling form = n^{1/\alpha} H_\alpha({m}/{n},{A}/{n^{1+1/\alpha}}), where at variance with Brownian bridge, H_\alpha(X,Y) is a non trivial function of the rescaled time m/n and rescaled area Y = A/n^{1+1/\alpha}. Our analytical results are verified by numerical simulations.Comment: 21 pages, 4 Figure

On h h -transforms of one-dimensional diffusions stopped upon hitting zero

For a one-dimensional diffusion on an interval for which 0 is the regular-reflecting left boundary, three kinds of conditionings to avoid zero are studied. The limit processes are $h$-transforms of the process stopped upon hitting zero, where $h$'s are the ground state, the scale function, and the renormalized zero-resolvent. Several properties of the $h$-transforms are investigated

On the exchange of intersection and supremum of sigma-fields in filtering theory

We construct a stationary Markov process with trivial tail sigma-field and a nondegenerate observation process such that the corresponding nonlinear filtering process is not uniquely ergodic. This settles in the negative a conjecture of the author in the ergodic theory of nonlinear filters arising from an erroneous proof in the classic paper of H. Kunita (1971), wherein an exchange of intersection and supremum of sigma-fields is taken for granted.Comment: 20 page