Intelligence in Scientific Computing

Combining numerical techniques with ideas from symbolic computation and with methods incorporating knowledge of science and mathematics leads to a new category of intelligent computational tools for scientists and engineers. These tools autonomously prepare simulation experiments from high-level specifications of physical models. For computationally intensive experiments, they automatically design special-purpose numerical engines optimized to perform the necessary computations. They actively monitor numerical and physical experiments. They interpret experimental data and formulate numerical results in qualitative terms. They enable their human users to control computational experiments in terms of high-level behavioral descriptions.

Variation of modal parameters of pile-supported structures in seismically liquefiable soils.

This paper presents experimental results that aimed to investigate the effects of soil liquefaction on the modal parameters (i.e. frequency and damping ratio) of pile-supported structures. The tests were carried out using the shaking table facility of the Bristol Laboratory for Advanced Dynamics Engineering (BLADE) at the University of Bristol (UK) whereby four pile-supported structures (two single piles and two pile groups) with and without superstructure mass were tested. The experimental investigation aimed to monitor the variation in natural frequency and damping of the four physical models at different degrees of excess pore water pressure generation and in full-liquefaction condition. The experimental results showed that the natural frequency of pile-supported structures may decrease considerably owing to the loss of lateral support offered by the soil to the pile. On the other hand, the damping ratio of structure may increase to values in excess of 20%. These findings have important design consequences: (a) for low-period structures, substantial reduction of spectral acceleration is expected; (b) during and after liquefaction, the response of the system may be dictated by the interactions of multiple loadings, that is, horizontal, axial and overturning moment, which were negligible prior to liquefaction; and (c) with the onset of liquefaction due to increased flexibility of pile-supported structure, larger spectral displacement may be expected, which in turn may enhance Pdelta effects and consequently amplification of overturning moment. Practical implications for pile design are discussed.

The dayside high-latitude trough under quiet geomagnetic conditions: Radio tomography and the CTIP model

Pryse, Sian; Dewis, K.L.; Middleton, H.R.; Balthazor, R.L., (2005) 'The dayside high-latitude trough under quiet geomagnetic conditions: Radio tomography and the CTIP model', Annales Geophysicae 23(4) pp.1199-1206 RAE2008

Modele fizyczne planetoid pojedynczych i podwójnych otrzymane na podstawie przeglądów fotometrycznych: przygotowanie do naukowego wykorzystania danych z kosmicznej misji Gaia

Wydział Fizyki

Modelling air quality in street canyons: a review

High pollution levels have been often observed in urban street canyons due to the increased traffic emissions and reduced natural ventilation. Microscale dispersion models with different levels of complexity may be used to assess urban air qualityand support decision-making for pollution control strategies and traffic planning. Mathematical models calculate pollutant concentrations by solving either analytically a simplified set of parametric equations or numerically a set of differential equations that describe in detail wind flow and pollutant dispersion. Street canyon models, which might also include simplified photochemistry and particle deposition–resuspension algorithms, are often nested within larger-scale urban dispersion codes. Reduced-scale physical models in wind tunnels may also be used for investigating atmospheric processes within urban canyons and validating mathematical models. A range of monitoring techniques is used to measure pollutant concentrations in urban streets. Point measurement methods (continuous monitoring, passive and active pre-concentration sampling, grab sampling) are available for gaseous pollutants. A number of sampling techniques (mainlybased on filtration and impaction) can be used to obtain mass concentration, size distribution and chemical composition of particles. A combination of different sampling/monitoring techniques is often adopted in experimental studies. Relativelysimple mathematical models have usually been used in association with field measurements to obtain and interpret time series of pollutant concentrations at a limited number of receptor locations in street canyons. On the other hand, advanced numerical codes have often been applied in combination with wind tunnel and/or field data to simulate small-scale dispersion within the urban canopy.

Combined vertex-based - cell-centred finite volume method for flows in complex geometries

CFD modelling of 'real-life' processes often requires solutions in complex three dimensional geometries, which can often result in meshes where aspects of it are badly distorted. Cell-centred finite volume methods, typical of most commercial CFD tools, are computationally efficient, but can lead to convergence problems on meshes which feature cells with high non-orthogonal shapes. The vertex-based finite volume method handles distorted meshes with relative ease, but is computationally expensive. A combined vertex-based - cell-centred (VB-CC) technique, detailed in this paper, allows solutions on distorted meshes that defeat purely cell-centred physical models to be employed in the solution of other transported quantities. The VB-CC method is validated with benchmark solutions for thermally driven flow and turbulent flow. An early application of this hybrid technique is to three-dimensional flow over an aircraft wing, although it is planned to use it in a wide variety of processing applications in the future.

Spin asymmetries in spatially resolved processes of ionization of heavy atoms by relativistic electrons

The triple-differential cross section for ionization of a heavy atom is shown to depend on the spin of the incident electron even if this is polarized entirely parallel or antiparallel to its direction of propagation, the atom is unpolarized, and the spins of the ejected electrons are not resolved. Quantitative predictions for the spin asymmetry are presented in a relativistic distorted-wave Born approximation. Simple physical models are introduced to understand both these results and further symmetry properties involving the reversal of a spatial momentum component also.

Design Methodology for Real-Time FPGA-Based Sound Synthesis

Explicit finite difference (FD) schemes can realise highly realistic physical models of musical instruments but are computationally complex. A design methodology is presented for the creation of FPGA-based micro-architectures for FD schemes which can be applied to a range of applications with varying computational requirements, excitation and output patterns and boundary conditions. It has been applied to membrane and plate-based sound producing models, resulting in faster than real-time performance on a Xilinx XC2VP50 device which is 10 to 35 times faster than general purpose and DSP processors. The models have developed in such a way to allow a wide range of interaction (by a musician) thereby leading to the possibility of creating a highly realistic digital musical instrument.

Thermophysical properties of ionic liquids

Ionic liquids have received significant interest from research groups and industry for a range of novel applications. Many of these require a thorough knowledge of the thermophysical properties of the pure fluids and their mixtures. Despite this need, the necessary experimental data for many properties are scarce and often inconsistent between the various sources. However, by using accurate data, predictive physical models can be developed which are highly useful, and some would consider essential, if ionic liquids are to realise their full potential. This is particularly true if one can use them to design new ionic liquids which maximise key desired attributes. This paper will review some of the recent advances in our understanding, prediction and correlation of selected ionic liquid physical properties.

Partial least squares model for dynamic rating of overhead line in wind intensive areas based on field measurements

This paper presents a statistical model for the thermal behaviour of the line model based on lab tests and field measurements. This model is based on Partial Least Squares (PLS) multi regression and is used for the Dynamic Line Rating (DLR) in a wind intensive area. DLR provides extra capacity to the line, over the traditional seasonal static rating, which makes it possible to defer the need for reinforcement the existing network or building new lines. The proposed PLS model has a number of appealing features; the model is linear, so it is straightforward to use for predicting the line rating for future periods using the available weather forecast. Unlike the available physical models, the proposed model does not require any physical parameters of the line, which avoids the inaccuracies resulting from the errors and/or variations in these parameters. The developed model is compared with physical model, the Cigre model, and has shown very good accuracy in predicting the conductor temperature as well as in determining the line rating for future time periods. 

Toward Characterization Of The Type IIP Supernova Progenitor Population: A Statistical Sample Of Light Curves From Pan-STARRS1

In recent years, wide-field sky surveys providing deep multi-band imaging have presented a new path for indirectly characterizing the progenitor populations of core-collapse supernovae (SN): systematic light curve studies. We assemble a set of 76 grizy-band Type IIP SN light curves from Pan-STARRS1, obtained over a constant survey program of 4 years and classified using both spectroscopy and machine learning-based photometric techniques. We develop and apply a new Bayesian model for the full multi-band evolution of each light curve in the sample. We find no evidence of a sub-population of fast-declining explosions (historically referred to as "Type IIL" SNe). However, we identify a highly significant relation between the plateau phase decay rate and peak luminosity among our SNe IIP. These results argue in favor of a single parameter, likely determined by initial stellar mass, predominantly controlling the explosions of red supergiants. This relation could also be applied for supernova cosmology, offering a standardizable candle good to an intrinsic scatter of 0.2 mag. We compare each light curve to physical models from hydrodynamic simulations to estimate progenitor initial masses and other properties of the Pan-STARRS1 Type IIP SN sample. We show that correction of systematic discrepancies between modeled and observed SN IIP light curve properties and an expanded grid of progenitor properties, are needed to enable robust progenitor inferences from multi-band light curve samples of this kind. This work will serve as a pathfinder for photometric studies of core-collapse SNe to be conducted through future wide field transient searches.

Multi-wavelength observations of Proxima Centauri

Aims: We report simultaneous observations of the nearby flare star Proxima Centauri with VLT/UVES and XMM-Newton over three nights in March 2009. Our optical and X-ray observations cover the star's quiescent state, as well as its flaring activity and allow us to probe the stellar atmospheric conditions from the photosphere into the chromosphere, and then the corona during its different activity stages. Methods: Using the X-ray data, we investigate variations in coronal densities and abundances and infer loop properties for an intermediate-sized flare. The optical data are used to investigate the magnetic field and its possible variability, to construct an emission line list for the chromosphere, and use certain emission lines to construct physical models of Proxima Centauri's chromosphere. Results: We report the discovery of a weak optical forbidden Fe xiii line at 3388 Å during the more active states of Proxima Centauri. For the intermediate flare, we find two secondary flare events that may originate in neighbouring loops, and discuss the line asymmetries observed during this flare in H i, He i, and Ca ii lines. The high time-resolution in the Hα line highlights strong temporal variations in the observed line asymmetries, which re-appear during a secondary flare event. We also present theoretical modelling with the stellar atmosphere code PHOENIX to construct flaring chromospheric models. Based on observations collected at the European Southern Observatory, Paranal, Chile, 082.D-0953A and on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member states and NASA.Full Table 6 is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/534/A133

Computing probable maximum loss in catastrophe reinsurance portfolios on multi-core and many-core architectures

In the reinsurance market, the risks natural catastrophes pose to portfolios of properties must be quantified, so that they can be priced, and insurance offered. The analysis of such risks at a portfolio level requires a simulation of up to 800 000 trials with an average of 1000 catastrophic events per trial. This is sufficient to capture risk for a global multi-peril reinsurance portfolio covering a range of perils including earthquake, hurricane, tornado, hail, severe thunderstorm, wind storm, storm surge and riverine flooding, and wildfire. Such simulations are both computation and data intensive, making the application of high-performance computing techniques desirable.

In this paper, we explore the design and implementation of portfolio risk analysis on both multi-core and many-core computing platforms. Given a portfolio of property catastrophe insurance treaties, key risk measures, such as probable maximum loss, are computed by taking both primary and secondary uncertainties into account. Primary uncertainty is associated with whether or not an event occurs in a simulated year, while secondary uncertainty captures the uncertainty in the level of loss due to the use of simplified physical models and limitations in the available data. A combination of fast lookup structures, multi-threading and careful hand tuning of numerical operations is required to achieve good performance. Experimental results are reported for multi-core processors and systems using NVIDIA graphics processing unit and Intel Phi many-core accelerators.

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.

Estudo hidráulico em modelo reduzido de descarregadores de cheias: utilização de conchas de rolo na dissipação de energia do escoamento

Dissertação para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização em Hidráulica