Boundary value problems for third-order linear PDEs in time-dependent domains

We present the extension of a methodology to solve moving boundary value problems from the second-order case to the case of the third-order linear evolution PDE qt + qxxx = 0. This extension is the crucial step needed to generalize this methodology to PDEs of arbitrary order. The methodology is based on the derivation of inversion formulae for a class of integral transforms that generalize the Fourier transform and on the analysis of the global relation associated with the PDE. The study of this relation and its inversion using the appropriate generalized transform are the main elements of the proof of our results.

Necessary and sufficient conditions for realizability of point processes

We give necessary and sufficient conditions for a pair of (generali- zed) functions 1(r1) and 2(r1, r2), ri 2X, to be the density and pair correlations of some point process in a topological space X, for ex- ample, Rd, Zd or a subset of these. This is an infinite-dimensional version of the classical “truncated moment” problem. Standard tech- niques apply in the case in which there can be only a bounded num- ber of points in any compact subset of X. Without this restriction we obtain, for compact X, strengthened conditions which are necessary and sufficient for the existence of a process satisfying a further re- quirement—the existence of a finite third order moment. We general- ize the latter conditions in two distinct ways when X is not compact.

A high-order arbitrarily unstructured finite-volume model of the global atmosphere: tests solving the shallow-water equations

Simulations of the global atmosphere for weather and climate forecasting require fast and accurate solutions and so operational models use high-order finite differences on regular structured grids. This precludes the use of local refinement; techniques allowing local refinement are either expensive (eg. high-order finite element techniques) or have reduced accuracy at changes in resolution (eg. unstructured finite-volume with linear differencing). We present solutions of the shallow-water equations for westerly flow over a mid-latitude mountain from a finite-volume model written using OpenFOAM. A second/third-order accurate differencing scheme is applied on arbitrarily unstructured meshes made up of various shapes and refinement patterns. The results are as accurate as equivalent resolution spectral methods. Using lower order differencing reduces accuracy at a refinement pattern which allows errors from refinement of the mountain to accumulate and reduces the global accuracy over a 15 day simulation. We have therefore introduced a scheme which fits a 2D cubic polynomial approximately on a stencil around each cell. Using this scheme means that refinement of the mountain improves the accuracy after a 15 day simulation. This is a more severe test of local mesh refinement for global simulations than has been presented but a realistic test if these techniques are to be used operationally. These efficient, high-order schemes may make it possible for local mesh refinement to be used by weather and climate forecast models.

Blind equalization of nonminimum phase channels: higher order cumulant based algorithm

Higher order cumulant analysis is applied to the blind equalization of linear time-invariant (LTI) nonminimum-phase channels. The channel model is moving-average based. To identify the moving average parameters of channels, a higher-order cumulant fitting approach is adopted in which a novel relay algorithm is proposed to obtain the global solution. In addition, the technique incorporates model order determination. The transmitted data are considered as independently identically distributed random variables over some discrete finite set (e.g., set {±1, ±3}). A transformation scheme is suggested so that third-order cumulant analysis can be applied to this type of data. Simulation examples verify the feasibility and potential of the algorithm. Performance is compared with that of the noncumulant-based Sato scheme in terms of the steady state MSE and convergence rate.

A Proposed Modification to the Robert–Asselin Time Filter

The Robert–Asselin time filter is widely used in numerical models of weather and climate. It successfully suppresses the spurious computational mode associated with the leapfrog time-stepping scheme. Unfortunately, it also weakly suppresses the physical mode and severely degrades the numerical accuracy. These two concomitant problems are shown to occur because the filter does not conserve the mean state, averaged over the three time slices on which it operates. The author proposes a simple modification to the Robert–Asselin filter, which does conserve the three-time-level mean state. When used in conjunction with the leapfrog scheme, the modification vastly reduces the impacts on the physical mode and increases the numerical accuracy for amplitude errors by two orders, yielding third-order accuracy. The modified filter could easily be incorporated into existing general circulation models of the atmosphere and ocean. In principle, it should deliver more faithful simulations at almost no additional computational expense. Alternatively, it may permit the use of longer time steps with no loss of accuracy, reducing the computational expense of a given simulation.

Simultaneous analysis of v1, v4, 2v2, v2+v5, and 2v5 bands of 12CH3F

The Fourier-transform spectrum of CH3F from 2800 to 3100 cm−1, obtained by Guelachvili in Orsay at a resolution of about 0.003 cm−1, was analyzed. The effective Hamiltonian used contained all symmetry allowed interactions up to second order in the Amat-Nielsen classification, together with selected third-order terms, amongst the set of nine vibrational basis functions represented by the states ν1(A1), ν4(E), 2ν2(A1), ν2 + ν5(E), 2ν50(A1), and 2ν5±2(E). A number of strong Fermi and Coriolis resonances are involved. The vibrational Hamiltonian matrix was not factorized beyond the requirements of symmetry. A total of 59 molecular parameters were refined in a simultaneous least-squares analysis to over 1500 upper-state energy levels for J ≤ 20 with a standard deviation of 0.013 cm−1. Although the standard deviation remains an order of magnitude greater than the precision of the measurements, this work breaks new ground in the simultaneous analysis of interacting symmetric top vibrational levels, in terms of the number of interacting vibrational states and the number of parameters in the Hamiltonian.

Molecular structure of 2,5-dihydropyrrole (C4NH7), obtained by gas-phase electron diffraction and theoretical calculations

The structure of 2,5-dihydropyrrole (C4NH7) has been determined by gas-phase electron diffraction (GED), augmented by the results from ab initio calculations employing third-order Moller-Plesset (MP3) level of theory and the 6-311+G(d,p) basis set. Several theoretical calculations were performed. From theoretical calculations using MP3/6-311+G(d,p) evidence was obtained for the presence of an axial (63%) (N-H bond axial to the CNC plane) and an equatorial conformer (37%) (N-H bond equatorial to the CNC plane). The five-membered ring was found to be puckered with the CNC plane inclined at 21.8 (38)° to the plane of the four carbon atoms.

Investigation of the prototype silylene reaction, SiH2+H2O(and D2O): time-resolved gas-phase kinetic studies, isotope effects, RRKM calculations, and quantum chemical calculations of the reaction energy surface

Time-resolved kinetic studies of the reaction of silylene, SiH2, with H2O and with D2O have been carried out in the gas phase at 296 and at 339 K, using laser flash photolysis to generate and monitor SiH2. The reaction was studied over the pressure range 10-200 Torr with SF6 as bath gas. The second-order rate constants obtained were pressure dependent, indicating that the reaction is a third-body assisted association process. Rate constants at 339 K were about half those at 296 K. Isotope effects, k(H)/k(D), were small averaging 1.076 0.080, suggesting no involvement of H- (or D-) atom transfer in the rate determining step. RRKM modeling was undertaken based on a transition state appropriate to formation of the expected zwitterionic donoracceptor complex, H2Si...OH2. Because the reaction is close to the low pressure (third order) region, it is difficult to be definitive about the activated complex structure. Various structures were tried, both with and without the incorporation of rotational modes, leading to values for the high-pressure limiting (i.e., true secondorder) rate constant in the range 9.5 x 10(-11) to 5 x 10(-10) cm(3) molecule' s(-1). The RRKM modeling and mechanistic interpretation is supported by ab initio quantum calculations carried out at the G2 and G3 levels. The results are compared and contrasted with the previous studies.

Experimental and theoretical evidence for homogeneous catalysis in the gas-phase reaction of SiH2 with H2O (and D2O): a combined kinetic and quantum chemical study

Time-resolved kinetic studies of the reaction of silylene, SiH2, with H2O and with D2O have been carried out in the gas phase at 297 K and at 345 K, using laser flash photolysis to generate and monitor SiH2. The reaction was studied independently as a function of H2O (or D2O) and SF6 (bath gas) pressures. At a fixed pressure of SF6 (5 Torr), [SiH2] decay constants, k(obs), showed a quadratic dependence on [H2O] or [D2O]. At a fixed pressure of H2O or D2O, k(obs) Values were strongly dependent on [SF6]. The combined rate expression is consistent with a mechanism involving the reversible formation of a vibrationally excited zwitterionic donor-acceptor complex, H2Si...OH2 (or H2Si...OD2). This complex can then either be stabilized by SF6 or it reacts with a further molecule of H2O (or D2O) in the rate-determining step. Isotope effects are in the range 1.0-1.5 and are broadly consistent with this mechanism. The mechanism is further supported by RRKM theory, which shows the association reaction to be close to its third-order region of pressure (SF6) dependence. Ab initio quantum calculations, carried out at the G3 level, support the existence of a hydrated zwitterion H2Si...(OH2)(2), which can rearrange to hydrated silanol, with an energy barrier below the reaction energy threshold. This is the first example of a gas-phase-catalyzed silylene reaction.

On the momentum fluxes associated with mountain waves in directionally sheared flows

The direct impact of mountain waves on the atmospheric circulation is due to the deposition of wave momentum at critical levels, or levels where the waves break. The first process is treated analytically in this study within the framework of linear theory. The variation of the momentum flux with height is investigated for relatively large shears, extending the authors’ previous calculations of the surface gravity wave drag to the whole atmosphere. A Wentzel–Kramers–Brillouin (WKB) approximation is used to treat inviscid, steady, nonrotating, hydrostatic flow with directional shear over a circular mesoscale mountain, for generic wind profiles. This approximation must be extended to third order to obtain momentum flux expressions that are accurate to second order. Since the momentum flux only varies because of wave filtering by critical levels, the application of contour integration techniques enables it to be expressed in terms of simple 1D integrals. On the other hand, the momentum flux divergence (which corresponds to the force on the atmosphere that must be represented in gravity wave drag parameterizations) is given in closed analytical form. The momentum flux expressions are tested for idealized wind profiles, where they become a function of the Richardson number (Ri). These expressions tend, for high Ri, to results by previous authors, where wind profile effects on the surface drag were neglected and critical levels acted as perfect absorbers. The linear results are compared with linear and nonlinear numerical simulations, showing a considerable improvement upon corresponding results derived for higher Ri.

The sample autocorrelation function and the detection of long-memory processes

The detection of long-range dependence in time series analysis is an important task to which this paper contributes by showing that whilst the theoretical definition of a long-memory (or long-range dependent) process is based on the autocorrelation function, it is not possible for long memory to be identified using the sum of the sample autocorrelations, as usually defined. The reason for this is that the sample sum is a predetermined constant for any stationary time series; a result that is independent of the sample size. Diagnostic or estimation procedures, such as those in the frequency domain, that embed this sum are equally open to this criticism. We develop this result in the context of long memory, extending it to the implications for the spectral density function and the variance of partial sums of a stationary stochastic process. The results are further extended to higher order sample autocorrelations and the bispectral density. The corresponding result is that the sum of the third order sample (auto) bicorrelations at lags h,k≥1, is also a predetermined constant, different from that in the second order case, for any stationary time series of arbitrary length.

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.

First order k-th moment finite element analysis of nonlinear operator equations with stochastic data

We develop and analyze a class of efficient Galerkin approximation methods for uncertainty quantification of nonlinear operator equations. The algorithms are based on sparse Galerkin discretizations of tensorized linearizations at nominal parameters. Specifically, we consider abstract, nonlinear, parametric operator equations J(\alpha ,u)=0 for random input \alpha (\omega ) with almost sure realizations in a neighborhood of a nominal input parameter \alpha _0. Under some structural assumptions on the parameter dependence, we prove existence and uniqueness of a random solution, u(\omega ) = S(\alpha (\omega )). We derive a multilinear, tensorized operator equation for the deterministic computation of k-th order statistical moments of the random solution's fluctuations u(\omega ) - S(\alpha _0). We introduce and analyse sparse tensor Galerkin discretization schemes for the efficient, deterministic computation of the k-th statistical moment equation. We prove a shift theorem for the k-point correlation equation in anisotropic smoothness scales and deduce that sparse tensor Galerkin discretizations of this equation converge in accuracy vs. complexity which equals, up to logarithmic terms, that of the Galerkin discretization of a single instance of the mean field problem. We illustrate the abstract theory for nonstationary diffusion problems in random domains.