A non-uniform mesh scheme for compressible flow

A non-uniform mesh scheme is presented for the computation of compressible flows governed by the Euler equations of gas dynamics. The scheme is based on flux-difference splitting and represents an extension of a similar scheme designed for uniform meshes. The numerical results demonstrate that little, if any, spurious oscillation occurs as a result of the non-uniformity of the mesh; and importantly, shock speeds are computed correctly.

Flux difference splitting for open channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.

A numerical scheme for two-dimensional, open channel flows with non-rectangular geometries

An algorithm based on flux difference splitting is presented for the solution of two-dimensional, open channel flows. A transformation maps a non-rectangular, physical domain into a rectangular one. The governing equations are then the shallow water equations, including terms of slope and friction, in a generalized coordinate system. A regular mesh on a rectangular computational domain can then be employed. The resulting scheme has good jump capturing properties and the advantage of using boundary/body-fitted meshes. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

A molecular T-matrix approach to calculating low-energy electron diffusion intensities for ordered molecular adsorbates

We present a novel approach to calculating Low-Energy Electron Diffraction (LEED) intensities for ordered molecular adsorbates. First, the intra-molecular multiple scattering is computed to obtain a non-diagonal molecular T-matrix. This is then used to represent the entire molecule as a single scattering object in a conventional LEED calculation, where the Layer Doubling technique is applied to assemble the different layers, including the molecular ones. A detailed comparison with conventional layer-type LEED calculations is provided to ascertain the accuracy of this scheme of calculation. Advantages of this scheme for problems involving ordered arrays of molecules adsorbed on surfaces are discussed.

Integration of a 3D Variational data assimilation scheme with a coastal area morphodynamic model of Morecambe Bay

This paper describes the implementation of a 3D variational (3D-Var) data assimilation scheme for a morphodynamic model applied to Morecambe Bay, UK. A simple decoupled hydrodynamic and sediment transport model is combined with a data assimilation scheme to investigate the ability of such methods to improve the accuracy of the predicted bathymetry. The inverse forecast error covariance matrix is modelled using a Laplacian approximation which is calibrated for the length scale parameter required. Calibration is also performed for the Soulsby-van Rijn sediment transport equations. The data used for assimilation purposes comprises waterlines derived from SAR imagery covering the entire period of the model run, and swath bathymetry data collected by a ship-borne survey for one date towards the end of the model run. A LiDAR survey of the entire bay carried out in November 2005 is used for validation purposes. The comparison of the predictive ability of the model alone with the model-forecast-assimilation system demonstrates that using data assimilation significantly improves the forecast skill. An investigation of the assimilation of the swath bathymetry as well as the waterlines demonstrates that the overall improvement is initially large, but decreases over time as the bathymetry evolves away from that observed by the survey. The result of combining the calibration runs into a pseudo-ensemble provides a higher skill score than for a single optimized model run. A brief comparison of the Optimal Interpolation assimilation method with the 3D-Var method shows that the two schemes give similar results.

Source splitting via the point source method

We introduce a new algorithm for source identification and field splitting based on the point source method (Potthast 1998 A point-source method for inverse acoustic and electromagnetic obstacle scattering problems IMA J. Appl. Math. 61 119–40, Potthast R 1996 A fast new method to solve inverse scattering problems Inverse Problems 12 731–42). The task is to separate the sound fields uj, j = 1, ..., n of sound sources supported in different bounded domains G1, ..., Gn in from measurements of the field on some microphone array—mathematically speaking from the knowledge of the sum of the fields u = u1 + + un on some open subset Λ of a plane. The main idea of the scheme is to calculate filter functions , to construct uℓ for ℓ = 1, ..., n from u|Λ in the form We will provide the complete mathematical theory for the field splitting via the point source method. In particular, we describe uniqueness, solvability of the problem and convergence and stability of the algorithm. In the second part we describe the practical realization of the splitting for real data measurements carried out at the Institute for Sound and Vibration Research at Southampton, UK. A practical demonstration of the original recording and the splitting results for real data is available online.

A Time Efficient Optical Model for GATE Simulation of a LYSO Scintillation Matrix Used in PET Applications

A time efficient optical model is proposed for GATE simulation of a LYSO scintillation matrix coupled to a photomultiplier. The purpose is to avoid the excessively long computation time when activating the optical processes in GATE. The usefulness of the model is demonstrated by comparing the simulated and experimental energy spectra obtained with the dual planar head equipment for dosimetry with a positron emission tomograph ( DoPET). The procedure to apply the model is divided in two steps. Firstly, a simplified simulation of a single crystal element of DoPET is used to fit an analytic function that models the optical attenuation inside the crystal. In a second step, the model is employed to calculate the influence of this attenuation in the energy registered by the tomograph. The use of the proposed optical model is around three orders of magnitude faster than a GATE simulation with optical processes enabled. A good agreement was found between the experimental and simulated data using the optical model. The results indicate that optical interactions inside the crystal elements play an important role on the energy resolution and induce a considerable degradation of the spectra information acquired by DoPET. Finally, the same approach employed by the proposed optical model could be useful to simulate a scintillation matrix coupled to a photomultiplier using single or dual readout scheme.

A class of improved heteroskedasticity-consistent covariance matrix estimators

The heteroskedasticity-consistent covariance matrix estimator proposed by White (1980), also known as HC0, is commonly used in practical applications and is implemented into a number of statistical software. Cribari–Neto, Ferrari & Cordeiro (2000) have developed a bias-adjustment scheme that delivers bias-corrected White estimators. There are several variants of the original White estimator that also commonly used by practitioners. These include the HC1, HC2 and HC3 estimators, which have proven to have superior small-sample behavior relative to White’s estimator. This paper defines a general bias-correction mechamism that can be applied not only to White’s estimator, but to variants of this estimator as well, such as HC1, HC2 and HC3. Numerical evidence on the usefulness of the proposed corrections is also presented. Overall, the results favor the sequence of improved HC2 estimators.

On the relationship between a 2 x 2 matrix and second-order scalar spectral problems for integrable equations

A simple proof is given that a 2 x 2 matrix scheme for an inverse scattering transform method for integrable equations can be converted into the standard form of the second-order scalar spectral problem associated with the same equations. Simple formulae relating these two kinds of representation of integrable equations are established.

On the Baker-Akhiezer function in the AKNS scheme

Expressions for the Baker-Akhiezer function and their logarithmic space and time derivatives are derived in terms of the matrix elements of U - V matrices and 'squared basis functions'. These expressions generalize the well known formulas for the KdV equation case and establish links between different forms of the Whitham averaging procedure.

REDUCTION OF TODA LATTICE HIERARCHY TO GENERALIZED KDV HIERARCHIES AND THE 2-MATRIX MODEL

Toda lattice hierarchy and the associated matrix formulation of the 2M-boson KP hierarchies provide a framework for the Drinfeld-Sokolov reduction scheme realized through Hamiltonian action within the second KP Poisson bracket. By working with free currents, which Abelianize the second KP Hamiltonian structure, we are able to obtain a unified formalism for the reduced SL(M + 1, M - k) KdV hierarchies interpolating between the ordinary KP and KdV hierarchies. The corresponding Lax operators are given as superdeterminants of graded SL(M + 1, M - k) matrices in the diagonal gauge and we describe their bracket structure and field content. In particular, we provide explicit free field representations of the associated W(M, M - k) Poisson bracket algebras generalising the familiar nonlinear W-M+1 algebra. Discrete Backlund transformations for SL(M + 1, M - k) KdV are generated naturally from lattice translations in the underlying Toda-like hierarchy. As an application we demonstrate the equivalence of the two-matrix string model to the SL(M + 1, 1) KdV hierarchy.

Spectroscopic studies of the Eu3+ and Er3+ ions in the fluorozirconate LaZr2F11 matrix

An investigation by optical spectroscopy of the Eu3 + and Er3 + active ions in the crystallized fluorozirconate matrix LaZr2F11 is presented. The D-5(1) --> F-7(0-5) emission lines of Eu3 + are used to extract the F-7(0-5) energy scheme and the observed extinctions permit the deduction of irreducible representations (IRREPS) associated with corresponding sub-levels in the D-2 symmetry. The crystal field analysis was carried out on a 387 x 387 basis set, comprising the F-7, D-5(1,2,3) F-5(1,2), (5)G(1,2,3) and P-3(1,2,3,4,5,6) terms of the Eu-3 (+) 4f(6) configuration. The deviation and rms are 6.8 and 7.9 cm (-1), respectively for 38 levels and ten parameters. The experimental crystal field parameters are in good agreement with the ab-initio ones. Moreover, the relative intensities of the D-5(0) --> F-7(2,3,4) emissions are well reproduced by an 'ab-initio' calculation, except for three lines. The Er3 + ions introduced in LaZr2F11, microcrystals also lie in an unique crystallographic site. A total of 31 energy levels were recorded and the crystal field analysis led to 6.6 and 7.8 cm (-1) for the deviation and rms, respectively, for nine variable parameters taken into account. The experimental CF parameters for Er3 + and Eu3 + are very similar, which seems to show that the host lattice contracts around the smaller Er3 + ion. The informations given by both Eu3 + and Eu3 + optical probes in LaZr2F11 are very consistent with the structure previously determined for the isotypic PrZr2F11 fluoride. (C) 2000 Elsevier B.V. B.V. All rights reserved.

A W-matrix methodology for solving sparse network equations on multiprocessor computers

This paper describes a methodology for solving efficiently the sparse network equations on multiprocessor computers. The methodology is based on the matrix inverse factors (W-matrix) approach to the direct solution phase of A(x) = b systems. A partitioning scheme of W-matrix , based on the leaf-nodes of the factorization path tree, is proposed. The methodology allows the performance of all the updating operations on vector b in parallel, within each partition, using a row-oriented processing. The approach takes advantage of the processing power of the individual processors. Performance results are presented and discussed.

Magnetic and optical investigation of 40SiO(2)center dot 30Na(2)O center dot 1Al(2)O(3)center dot(29-x) B2O3 center dot xFe(2)O(3) glass matrix

Samples of 40SiO(2)center dot 30Na(2)O center dot 1Al(2)O(3)center dot(29 - x)B2O3 center dot xFe(2)O(3) (mol%), with 0.0 <= x <= 17.5, were prepared by the fusion method and investigated by electron paramagnetic resonance (EPR), optical absorption (OA) and Mossbauer spectroscopy (MS). The EPR spectra of the as-synthesized samples exhibit two well-defined EPR signals around g = 4.27 and g = 2.01 and a visible EPR shoulder around g = 6.4, assigned to isolated Fe3+ ion complexes (g = 4.27 and g = 6.4) and Fe3+-based clusters (g = 2.01). Analyses of both EPR line intensity and line width support the model picture of Fe3+-based clusters built in from two sources of isolated ions, namely Fe2+ and Fe3+; the ferrous ion being used to build in iron-based clusters at lower x-content (below about x = 2.5%) whereas the ferric ion is used to build in iron-based clusters at higher x-content (above about x = 2.5%). The presence of Fe2+ ions incorporated within the glass template is supported by OA data with a strong band around 1100 nm due to the spin-allowed E-5(g)-T-5(2g) transition in an octahedral coordination with oxygen. Additionally, Mossbauer data (isomer shift and quadrupole splitting) confirm incorporation of both Fe2+ and Fe3+ ions within the template, more likely in tetrahedral-like environments. We hypothesize that ferrous ions are incorporated within the glass template as FeO4 complex resulting from replacing silicon in non-bridging oxygen (SiO3O-) sites whereas ferric ions are incorporated as FeO4 complex resulting from replacing silicon in bridging-like oxygen silicate groups (SiO4). (C) 2012 Elsevier Masson SAS. All rights reserved.

Matrix Designs and Methods for Secure and Efficient Compressed Sensing

The idea of balancing the resources spent in the acquisition and encoding of natural signals strictly to their intrinsic information content has interested nearly a decade of research under the name of compressed sensing. In this doctoral dissertation we develop some extensions and improvements upon this technique's foundations, by modifying the random sensing matrices on which the signals of interest are projected to achieve different objectives. Firstly, we propose two methods for the adaptation of sensing matrix ensembles to the second-order moments of natural signals. These techniques leverage the maximisation of different proxies for the quantity of information acquired by compressed sensing, and are efficiently applied in the encoding of electrocardiographic tracks with minimum-complexity digital hardware. Secondly, we focus on the possibility of using compressed sensing as a method to provide a partial, yet cryptanalysis-resistant form of encryption; in this context, we show how a random matrix generation strategy with a controlled amount of perturbations can be used to distinguish between multiple user classes with different quality of access to the encrypted information content. Finally, we explore the application of compressed sensing in the design of a multispectral imager, by implementing an optical scheme that entails a coded aperture array and Fabry-Pérot spectral filters. The signal recoveries obtained by processing real-world measurements show promising results, that leave room for an improvement of the sensing matrix calibration problem in the devised imager.