This study investigates the phenomenon of targeted energy transfer (TET) from a linear oscillator to a nonlinear attachment behaving as a nonlinear energy sink for both transient and stochastic excitations. First, the dynamics of the underlying Hamiltonian system under deterministic transient loading is studied. Assuming that the transient dynamics can be partitioned into slow and fast components, the governing equations of motion corresponding to the slow flow dynamics are derived and the behaviour of the system is analysed. Subsequently, the effect of noise on the slow flow dynamics of the system is investigated. The Itô stochastic differential equations for the noisy system are derived and the corresponding Fokker–Planck equations are numerically solved to gain insights into the behaviour of the system on TET. The effects of the system parameters as well as noise intensity on the optimal regime of TET are studied. The analysis reveals that the interaction of nonlinearities and noise enhances the optimal TET regime as predicted in deterministic analysis.

Electromagnetic wave propagation in complex dispersive media is governed by the time dependent Maxwell's equations coupled to equations that describe the evolution of the induced macroscopic polarization. We consider “polydispersive” materials represented by distributions of dielectric parameters in a polarization model. The work focuses on a novel computational framework for such problems involving Polynomial Chaos Expansions as a method to improve the modeling accuracy of the Debye model and allow for easy simulation using the Finite Difference Time Domain (FDTD) method. Stability and dispersion analyzes are performed for the approach utilizing the second order Yee scheme in two spatial dimensions.

We study the discrete-time approximation of doubly reflected backward stochastic differential equations (BSDEs) in a multidimensional setting. As in Ma and Zhang (2005) or Bouchard and Chassagneux (2008), we introduce the discretely reflected counterpart of these equations. We then provide representation formulae which allow us to obtain new regularity results. We also propose an Euler scheme type approximation and give new convergence results for both discretely and continuously reflected BSDEs.

This paper introduces a method of generating real harmonizable multifractional Lévy motions (RHMLMs). The simulation of these fields is closely related to that of infinitely divisible laws or Lévy processes. In the case where the control measure of the RHMLM is finite, generalized shot-noise series are used. An estimation of the error is also given. Otherwise, the RHMLM Xh is split into two independent RHMLMs, Xε,1 and Xε,2. More precisely, Xε,2 is an RHMLM whose control measure is finite. It can then be rewritten as a generalized shot-noise series. The asymptotic behaviour of Xε,1 as ε → 0+ is further elaborated. Sufficient conditions to approximate Xε,1 by a multifractional Brownian motion are given. The error rate in terms of Berry-Esseen bounds is then discussed. Finally, some examples of simulation are given.

A new procedure that generates the transient solution of the first moment of the state of a Markovian queueing network with state-dependent arrivals, services, and routeing is developed. The procedure involves defining a partial differential equation that relates an approximate multivariate cumulant generating function to the intensity functions of the network. The partial differential equation then yields a set of ordinary differential equations which are numerically solved to obtain the first moment.

Let X = (X(t):t ≥ 0) be a Lévy process and X∊ the compensated sum of jumps not exceeding ∊ in absolute value, σ2(∊) = var(X∊(1)). In simulation, X - X∊ is easily generated as the sum of a Brownian term and a compound Poisson one, and we investigate here when X∊/σ(∊) can be approximated by another Brownian term. A necessary and sufficient condition in terms of σ(∊) is given, and it is shown that when the condition fails, the behaviour of X∊/σ(∊) can be quite intricate. This condition is also related to the decay of terms in series expansions. We further discuss error rates in terms of Berry-Esseen bounds and Edgeworth approximations.

In this paper, we provide a generalization of the classical Rao bound for orthogonal arrays, which can be applied to ordered orthogonal arrays and $\left( t,\,m,\,s \right)$-nets. Application of our new bound leads to improvements in many parameter situations to the strongest bounds (i.e., necessary conditions) for existence of these objects.

In an earlier paper [10], we studied a generalized Rao bound for ordered orthogonal arrays and $(T,\,M,\,S)$-nets. In this paper, we extend this to a coding-theoretic approach to ordered orthogonal arrays. Using a certain association scheme, we prove a MacWilliams-type theorem for linear ordered orthogonal arrays and linear ordered codes as well as a linear programming bound for the general case. We include some tables which compare this bound against two previously known bounds for ordered orthogonal arrays. Finally we show that, for even strength, the $\text{LP}$ bound is always at least as strong as the generalized Rao bound.

In this paper statistical properties of estimators of drift parameters for diffusion processes are studied by modern numerical methods for stochastic differential equations. This is a particularly useful method for discrete time samples, where estimators can be constructed by making discrete time approximations to the stochastic integrals appearing in the maximum likelihood estimators for continuously observed diffusions. A review is given of the necessary theory for parameter estimation for diffusion processes and for simulation of diffusion processes. Three examples are studied.

Given a parametric family of regenerative processes on a common probability space, we investigate when the derivatives (with respect to the parameter) are regenerative. We primarily consider sequences satisfying explicit, Lipschitz recursions, such as the waiting times in many queueing systems, and show that derivatives regenerate together with the original sequence under reasonable monotonicity or continuity assumptions. The inputs to our recursions are i.i.d. or, more generally, governed by a Harris-ergodic Markov chain. For i.i.d. input we identify explicit regeneration points; otherwise, we use coupling arguments. We give conditions for the expected steady-state derivative to be the derivative of the steady-state mean of the original sequence. Under these conditions, the derivative of the steady-state mean has a cycle-formula representation.