An Exponentially Convergent Nonpolynomial Finite Element Method for Time-Harmonic Scattering from Polygons

In recent years nonpolynomial finite element methods have received increasing attention for the efficient solution of wave problems. As with their close cousin the method of particular solutions, high efficiency comes from using solutions to the Helmholtz equation as basis functions. We present and analyze such a method for the scattering of two-dimensional scalar waves from a polygonal domain that achieves exponential convergence purely by increasing the number of basis functions in each element. Key ingredients are the use of basis functions that capture the singularities at corners and the representation of the scattered field towards infinity by a combination of fundamental solutions. The solution is obtained by minimizing a least-squares functional, which we discretize in such a way that a matrix least-squares problem is obtained. We give computable exponential bounds on the rate of convergence of the least-squares functional that are in very good agreement with the observed numerical convergence. Challenging numerical examples, including a nonconvex polygon with several corner singularities, and a cavity domain, are solved to around 10 digits of accuracy with a few seconds of CPU time. The examples are implemented concisely with MPSpack, a MATLAB toolbox for wave computations with nonpolynomial basis functions, developed by the authors. A code example is included.

Time-resolved gas-phase kinetic study of the germylene addition reaction, GeH2 + C2D2

Time-resolved studies of germylene, GeH2, generated by laser. ash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate coefficients for its bimolecular reaction with C2D2. The reaction was studied in the gas phase, mainly at a total pressure of 1.3 kPa (in SF6 bath gas) at five temperatures in the range 298-558 K. Pressure variation measurements over the range 0.13-13 kPa ( SF6) at 298, 397 and 558 K revealed a small pressure dependence but only at 558 K. After correction for this, the second-order rate coefficients gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.96 +/- 0.05) + ( 6.16 +/- 0.37 kJ mol(-1))/RT ln 10 Comparison with the reaction of GeH2 + C2H2 (studied earlier) showed a similar behaviour with almost identical rate coefficients. The lack of a significant isotope effect is consistent with a rate-determining addition process and is explained by irreversible decomposition of the reaction intermediate to give Ge(P-3) + C2H4. This result contrasts with that for GeH2 + C2H4/C2D4 and those for the analogous silylene reactions. The underlying reasons for this are discussed.

Time-resolved gas-phase kinetic study of the germylene addition reaction, CH3C CCH3

Time-resolved kinetic studies of the reaction of germylene, GeH2, generated by laser. ash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reaction with 2-butyne, CH3C CCH3. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 300-556 K. The second order rate constants obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.46 +/- 10.06) + (5.16 +/- 10.47) kJ mol(-1)/ RT ln 10 Calculations of the energy surface of the GeC4H8 reaction system were carried out employing the additivity principle, by combining previous quantum chemical calculations of related reaction systems. These support formation of 1,2-dimethylvinylgermylene (rather than 2,3-dimethylgermirene) as the end product. RRKM calculations of the pressure dependence of the reaction are in reasonable agreement with this finding. The reactions of GeH2 with C2H2 and with CH3CRCCH3 are compared and contrasted.

Direct time-resolved study of the kinetics of the reaction of silylene with phenylsilane in the gas phase. Does SiH2 react with the aromatic ring?

Laser flash photolysis studies of silylene, SiH2, generated by the 193 nm laser flash photolysis phenylsilane, PhSiH3, have been carried out to obtain rate constants for its bimolecular reaction with PhSiH3 itself, in the gas phase. The reaction was studied in SF6 (mostly at 10 Torr total pressure) over the temperature range 298-595 K. The rate constants (also found to be pressure independent) gave the following Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-9.92 +/- 0.04) + (3.31 +/- 0.27) kJ mol(-1)/RT ln 10 Similar investigations of the reaction of silylene with benzene, C6H6, (295-410 K) gave data suggestive of the fact that SiH2 might be reacting with photochemical products of C6H6 as well as with C6H6 itself. However, in the latter system, apparent rate constants were sufficiently low to indicate that in the reaction of SiH2 with PhSiH3 addition to the aromatic ring was unlikely to be in excess of 3% of the total. Quantum chemical calculations of the energy surface for SiH2 + C6H6 indicate that 7-silanorcaradiene and 7-silacycloheptatriene are possible products but that PhSiH3 formation is unlikely. RRKM calculations suggest that 7-silanorcaradiene should be the initial product but that it cannot be collisionally stabilized under experimental conditions

An improved digital deconvolution equation

This paper first points out the important fact that the rectangle formulas of continuous convolution discretization, which was widely used in conventional digital deconvolution algorithms, can result in zero-time error. Then, an improved digital deconvolution equation is suggested which is equivalent to the trapezoid formulas of continuous convolution discretization and can overcome the disadvantage of conventional equation satisfactorily. Finally, a simulation in computer is given, thus confirming the theoretical result.

Generalised Dirichelt-to-Neumann map in time dependent domains

We study the heat, linear Schrodinger and linear KdV equations in the domain l(t) < x < ∞, 0 < t < T, with prescribed initial and boundary conditions and with l(t) a given differentiable function. For the first two equations, we show that the unknown Neumann or Dirichlet boundary value can be computed as the solution of a linear Volterra integral equation with an explicit weakly singular kernel. This integral equation can be derived from the formal Fourier integral representation of the solution. For the linear KdV equation we show that the two unknown boundary values can be computed as the solution of a system of linear Volterra integral equations with explicit weakly singular kernels. The derivation in this case makes crucial use of analyticity and certain invariance properties in the complex spectral plane. The above Volterra equations are shown to admit a unique solution.

The influence of the carbon surface chemical composition on Dubinin-Astakhov equation parameters calculated from SF(6) adsorption data-grand canonical Monte Carlo simulation

Using grand canonical Monte Carlo simulation we show, for the first time, the influence of the carbon porosity and surface oxidation on the parameters of the Dubinin-Astakhov (DA) adsorption isotherm equation. We conclude that upon carbon surface oxidation, the adsorption decreases for all carbons studied. Moreover, the parameters of the DA model depend on the number of surface oxygen groups. That is why in the case of carbons containing surface polar groups, SF(6) adsorption isotherm data cannot be used for characterization of the porosity.

A stress relaxation model for the viscoelastic solids based on the steady-state creep equation

Creep and stress relaxation are inherent mechanical behaviors of viscoelastic materials. It is considered that both are different performances of one identical physical phenomenon. The relationship between the decay stress and time during stress relaxation has been derived from the power law equation of the steady-state creep. The model was used to analyse the stress relaxation curves of various different viscoelastic materials (such as pure polycrystalline ice, polymers, foods, bones, metal, animal tissues, etc.). The calculated results using the theoretical model agree with the experimental data very well. Here we show that the new mathematical formula is not only simple but its parameters have the clear physical meanings. It is suitable to materials with a very broad scope and has a strong predictive ability.

Error analysis of Trefftz-discontinuous Galerkin methods for the time-harmonic Maxwell equations

In this paper, we extend to the time-harmonic Maxwell equations the p-version analysis technique developed in [R. Hiptmair, A. Moiola and I. Perugia, Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of the p-version, SIAM J. Numer. Anal., 49 (2011), 264-284] for Trefftz-discontinuous Galerkin approximations of the Helmholtz problem. While error estimates in a mesh-skeleton norm are derived parallel to the Helmholtz case, the derivation of estimates in a mesh-independent norm requires new twists in the duality argument. The particular case where the local Trefftz approximation spaces are built of vector-valued plane wave functions is considered, and convergence rates are derived.

On the long time behavior of free stochastic Schrödinger evolutions

We discuss the time evolution of the wave function which is the solution of a stochastic Schrödinger equation describing the dynamics of a free quantum particle subject to spontaneous localizations in space. We prove global existence and uniqueness of solutions. We observe that there exist three time regimes: the collapse regime, the classical regime and the diffusive regime. Concerning the latter, we assert that the general solution converges almost surely to a diffusing Gaussian wave function having a finite spread both in position as well as in momentum. This paper corrects and completes earlier works on this issue.

A time-domain probe method for three-dimensional rough surface reconstructions

The task of this paper is to develop a Time-Domain Probe Method for the reconstruction of impenetrable scatterers. The basic idea of the method is to use pulses in the time domain and the time-dependent response of the scatterer to reconstruct its location and shape. The method is based on the basic causality principle of timedependent scattering. The method is independent of the boundary condition and is applicable for limited aperture scattering data. In particular, we discuss the reconstruction of the shape of a rough surface in three dimensions from time-domain measurements of the scattered field. In practise, measurement data is collected where the incident field is given by a pulse. We formulate the time-domain fieeld reconstruction problem equivalently via frequency-domain integral equations or via a retarded boundary integral equation based on results of Bamberger, Ha-Duong, Lubich. In contrast to pure frequency domain methods here we use a time-domain characterization of the unknown shape for its reconstruction. Our paper will describe the Time-Domain Probe Method and relate it to previous frequency-domain approaches on sampling and probe methods by Colton, Kirsch, Ikehata, Potthast, Luke, Sylvester et al. The approach significantly extends recent work of Chandler-Wilde and Lines (2005) and Luke and Potthast (2006) on the timedomain point source method. We provide a complete convergence analysis for the method for the rough surface scattering case and provide numerical simulations and examples.

Power allocation strategies for distributed space-time codes in two-way relay networks

We study a two-way relay network (TWRN), where distributed space-time codes are constructed across multiple relay terminals in an amplify-and-forward mode. Each relay transmits a scaled linear combination of its received symbols and their conjugates,with the scaling factor chosen based on automatic gain control. We consider equal power allocation (EPA) across the relays, as well as the optimal power allocation (OPA) strategy given access to instantaneous channel state information (CSI). For EPA, we derive an upper bound on the pairwise-error-probability (PEP), from which we prove that full diversity is achieved in TWRNs. This result is in contrast to one-way relay networks, in which case a maximum diversity order of only unity can be obtained. When instantaneous CSI is available at the relays, we show that the OPA which minimizes the conditional PEP of the worse link can be cast as a generalized linear fractional program, which can be solved efficiently using the Dinkelback-type procedure.We also prove that, if the sum-power of the relay terminals is constrained, then the OPA will activate at most two relays.

On the stability and convergence of the finite section method for integral equation formulations of rough surface scattering

We consider the Dirichlet and Robin boundary value problems for the Helmholtz equation in a non-locally perturbed half-plane, modelling time harmonic acoustic scattering of an incident field by, respectively, sound-soft and impedance infinite rough surfaces.Recently proposed novel boundary integral equation formulations of these problems are discussed. It is usual in practical computations to truncate the infinite rough surface, solving a boundary integral equation on a finite section of the boundary, of length 2A, say. In the case of surfaces of small amplitude and slope we prove the stability and convergence as A→∞ of this approximation procedure. For surfaces of arbitrarily large amplitude and/or surface slope we prove stability and convergence of a modified finite section procedure in which the truncated boundary is ‘flattened’ in finite neighbourhoods of its two endpoints. Copyright © 2001 John Wiley & Sons, Ltd.

Time valued life cycle greenhouse gas emissions from buildings

The UK has adopted legally binding carbon reduction targets of 34% by 2020 and 80% by 2050 (measured against the 1990 baseline). Buildings are estimated to be responsible for more than 50% of greenhouse gas (GHG) emissions in the UK. These consist of both operational, produced during use, and embodied, produced during manufacture of materials and components, and during construction, refurbishments and demolition. A brief assessment suggests that it is unlikely that UK emission reduction targets can be met without substantial reductions in both Oc and Ec. Oc occurs over the lifetime of a building whereas the bulk of Ec occurs at the start of a building’s life. A time value for emissions could influence the decision making process when it comes to comparing mitigation measures which have benefits that occur at different times. An example might be the choice between building construction using low Ec construction materials versus building construction using high Ec construction materials but with lower Oc, although the use of high Ec materials does not necessarily imply a lower Oc. Particular time related issues examined here are: the urgency of the need to achieve large emissions reductions during the next 10 to 20 years; the earlier effective action is taken, the less costly it will be; future reduction in carbon intensity of energy supply; the carbon cycle and relationship between the release of GHG’s and their subsequent concentrations in the atmosphere. An equation is proposed, which weights emissions according to when they occur during the building life cycle, and which effectively increases Ec as a proportion of the total, suggesting that reducing Ec is likely to be more beneficial, in terms of climate change, for most new buildings. Thus, giving higher priority to Ec reductions is likely to result in a bigger positive impact on climate change and mitigation costs.

Instantaneous gelation in Smoluchowski’s coagulation equation revisited

We study the solutions of the Smoluchowski coagulation equation with a regularization term which removes clusters from the system when their mass exceeds a specified cutoff size, M. We focus primarily on collision kernels which would exhibit an instantaneous gelation transition in the absence of any regularization. Numerical simulations demonstrate that for such kernels with monodisperse initial data, the regularized gelation time decreasesas M increases, consistent with the expectation that the gelation time is zero in the unregularized system. This decrease appears to be a logarithmically slow function of M, indicating that instantaneously gelling kernels may still be justifiable as physical models despite the fact that they are highly singular in the absence of a cutoff. We also study the case when a source of monomers is introduced in the regularized system. In this case a stationary state is reached. We present a complete analytic description of this regularized stationary state for the model kernel, K(m1,m2)=max{m1,m2}ν, which gels instantaneously when M→∞ if ν>1. The stationary cluster size distribution decays as a stretched exponential for small cluster sizes and crosses over to a power law decay with exponent ν for large cluster sizes. The total particle density in the stationary state slowly vanishes as [(ν−1)logM]−1/2 when M→∞. The approach to the stationary state is nontrivial: Oscillations about the stationary state emerge from the interplay between the monomer injection and the cutoff, M, which decay very slowly when M is large. A quantitative analysis of these oscillations is provided for the addition model which describes the situation in which clusters can only grow by absorbing monomers.