Computational Modelling Group

A Fast Multipole Method for the Bessel potential

Marc Molinari, Simon Cox (Investigators), Neil O'Brien

The fast multipole method (FMM) proposed by Greengard and Rokhlin provides a method by which the O(N-squared) many-body problem can be reduced to O(N) complexity. In this project, a multipole method is developed to calculate the energy of a system of vortices in a high temperature superconductor, where the many-body interactions give rise to rich and complex physics. The method developed here is suitable for systems where the interactions are governed by a Bessel potential rather than the usual logarithmic potentials occurring in gravitational and electrostatic problems. We derive and apply vectorised forms of the Gegenbauer addition formulae in order to achieve the O(N) scaling associated with fast multipole methods.

A Fortran Based Mesh Viewer

Gabriel Amine-Eddine (Investigator)

During my final year as an undergraduate, I developed a fully functional software program for visualising geometries, grids and grid quality aspects, from a customised CFD software tool. It is designed to be used in conjunction with the HARTREE CFD code, and is currently in use by other fellow researchers within my supervisors team.

Ab initio simulations of chemical reactions on platinum nanoparticles

Chris-Kriton Skylaris (Investigator), Alvaro Ruiz-Serrano, Peter Cherry

•Use first principles calculations to study the relationship between shape and size of nanoparticle and the oxygen adsorption energy.

• Investigate the effect of high oxygen coverage on the catalytic activity of the nanoparticles.

An Investigation into the Cascade Effect of Mergers on the Global Financial Markets

Seth Bullock, Antonella Ianni (Investigators), Camillia Zedan

An investigation into the external effects that horizontal mergers have on the interconnected global markets.

Bioinformatic identification and physiological analysis of ethanol-related genes in C. elegans

Richard Edwards, Vincent O'Connor, Lindy Holden-Dye (Investigators), Ben Ient

Investigating the broad molecular, cellular and systems level impacts of acute and chronic ethanol in the nematode, Caenorhabditis elegans, as a model.

Cavity-Mediated Cooling

Peter Horak, Timothy Freegarde (Investigators), Andre Xuereb

Optical resonators enhance the interaction of light with matter while simultaneously acting as a temporal buffer. Both effects can be exploited to generate light-induced friction, or cooling, forces on atoms, molecules, or micromirrors. We investigate various aspects of these effects through numerical simulations, assisted by approximate analytical models, in this EPSRC and ESF sponsored project.

Continuously Tunable Optical Buffer

Peter Horak (Investigator)

The project aims to design, fabricate and test a novel integrated all-optical buffer device that is based on MEMS technology and provides a continuously tunable delay for optical pulses over a broad wavelength region. Such a device could play a crucial role in future packet-switched optical networks, photonic integrated circuits and coherent light based applications such as optically steered phase array antennas, LIDAR and optical coherence tomography.

This EPSRC funded project is a collaboration between the Optoelectronics Research Centre, Southampton, and University College London.

Control and Prediction of the Organic Solid State

Richard Boardman

This project aims to produce a computer technology for the prediction of the crystal structure(s) of an organic molecule, that could be used even prior to the synthesis of the compound.

Such a computational study could be done relatively quickly to predict the dangers and opportunities of the solid phases of a molecule under development. Our project will develop the methods of experimental screening for polymorphs and their characterisation, and hence the combination will provide a major new technology for aiding industrial formulation.

Coronary Artery Stent Design for Challenging Disease

Neil Bressloff (Investigator), Georgios Ragkousis

In this work, a method has been setup to (i) reconstruct diseased patient specific coronary artery segments; (ii) simulate the deployment of state of the start stents into these segments and (iii) assess the degree of stent malapposition. The aim now is to devise a stent delivery system that can mitigate this problem

DePuy Technology Partnership

Mark Taylor (Investigator), Adam Briscoe

This initiative concerns the transfer of knowledge between three key institutions (University of Southampton, University of Leeds and University of Hamburg) and DePuy International limited. The project is concerned with the ongoing advancement of technology used in orthopaedic devices.

Desiging Near-Capacity Quantum Error Correction Codes

Lajos Hanzo (Investigator), Zunaira Babar

Design efficient quantum error correction codes to correct the errors encountered in a quantum transmission; thus, increasing reliability and robustness of the future quantum systems.

Discrete ECogeomorphic Aeolian Landscape (DECAL) modelling

Joanna Nield (Investigator)

DECAL is a cellular automaton based model which incorporated mutual feedback processes between geomorphic forcing and ecological growth to investigate fundamental controls, self-organising and non-linear behaviour in semi-arid aeolian dune environments. This project explores landscape evolution and disturbance response, developing a phase-space in which dune fields can be quantified.

Electrostatic embedded energy calculations of proteins, using the ONETEP DFT code

Chris-Kriton Skylaris (Investigator), Stephen Fox, Chris Pittock

Calculating the energy of a biomolecule in solvent, using quantum mechanics (QM) is possible, but extremely challenging, even with linear-scaling QM methods like ONETEP. Using electrostatic embedding, a novel twist on the existing QM/MM method is used to calculate the binding energy of a small ligand to a solvated protein, increasing the accuracy and realism of our general project work.

Fluid Loads and Motions of Damaged Ships

Dominic Hudson, Ming-yi Tan (Investigators), Christian Wood, James Underwood, Adam Sobey

An area of research currently of interest in the marine industry is the effect of damage on ship structures. Research into the behaviour of damaged ships began in the mid nineties as a result of Ro-Ro disasters (e.g. Estonia in 1994). Due to the way the Estonia sank early research mainly focused on transient behaviour immediately after the damage takes place, the prediction of capsize, and of large lateral motions. Further research efforts, headed by the UK MoD, began following an incident where HMS Nottingham ran aground tearing a 50m hole from bow to bridge, flooding five compartments and almost causing the ship to sink just off Lord Howe Island in 2002. This project intends to answer the following questions:
“For a given amount of underwater damage (e.g. collision or torpedo/mine hit), what will be the progressive damage spread if the ship travels at ‘x’ knots? OR for a given amount of underwater damage, what is the maximum speed at which the ship can travel without causing additional damage?”

Fundamental Investigations of Cross-Coupled, Particle-Turbulence Interactions using a Pseudo Spectral DNS Code

Gabriel Amine-Eddine (Investigator), John Shrimpton

The behaviour of multiphase flows is of primary importance in many engineering applications. In the past, experimental observations have provided many researchers with the ability to understand and probe the phenomena and physical processes occurring in such flows. With advancements in modern day computational power, we now have the ability to gain an even greater wealth of knowledge, from what used to be a physical experiment, is now a virtual experiment, running across multiple computers in parallel architectures.

In this project, we simulate the full Navier-Stokes equations in a virtual experiment, and resolve to the best of degree, all possible scales of turbulence. We have the capability to track millions of computational particles in conjunction with the turbulence, and if the particles are charged, coupled to the turbulence, or if gravity is in the scenarios, we can examine the complex physical processes that occur in such a flow.

Provision has been made to simulate particles in conjunction with turbulence that has been subjected to deformations due to shear, strain, axi-symmetric contraction or expansion. Advancements in this code are soon to include the transport and coupling of scalar temperature between particles and the turbulence.

Currently, focus is on the coupling behaviour of poly-sized particle with the turbulence, and how such turbulence can be modelled accurately using stochastic Langevin methods.

Generic Operational Simulation of Civil Unmanned Air Vehicle Operations

Kenji Takeda, James Scanlan (Investigators)

This project creates a generic operational simulation of Unmanned Air Vehicle Operations. UAVs can be valued for their mission-suitability and compared against various configurations.

Graphical Simulation of Archaeological Environments

Graeme Earl (Investigator)

This project defines an emerging area of interest in physically accurate rendering within the Archaeological Computing Research Group. Sub-projects include analysis of Roman spaces at herculaneum, Neolithic buildings at Catalhoyuk and simulation of a range of artefacts.

Identification of novel Crustacean Pathogen Receptor Proteins

Richard Edwards, Chris Hauton, Timothy Elliott (Investigators), Oyindamola Lawal, Lloyd Mushambadzi

We are mining EST libraries (sequence fragments of expressed genes) for novel proteins that might play a role in the immune response of crustaceans.

Immunotherapy Research: Modelling MHC Class I Complex Assembly

Timothy Elliott, Jorn Werner (Investigators), Alistair Bailey

This project uses mathematical modelling and simulation to investigate mechanisms by which our cells process and present biological information that is used by our immune system to distinguish between healthy and diseased cells.

Integrated in silico prediction of protein-protein interaction motifs

Richard Edwards (Investigator), Nicolas Palopoli, Kieren Lythgow

Many vital protein-protein interactions are mediated by Short Linear Motifs (SLiMs) which are short proteins typically 5-15 amino acids long containing only a few positions crucial to function. This project integrates a number of leading computational techniques to predict novel SLiMs and add crucial detail to protein-protein interaction networks.

Interactome-wide prediction of short linear protein interaction motifs in humans

Richard Edwards (Investigator)

Short Linear Motifs (SLiMs) are important in many protein-protein interactions. In previous work, we have developed a computational tool, SLiMFinder, which places the interpretation of evidence for motifs within a statistical framework with high specificity, and subsequently enhanced sensitivity through application of conservation-based sequence masking. We are now applying these tools to a comprehensive set of human protein-protein interactions in order to predict novel human SLiMs in silico.

Measuring biomolecules - improvements to the spectroscopic ruler

Pavlos Lagoudakis, Tom Brown (Investigators), Jan Junis Rindermann, James Richardson

The spectroscopic ruler is a technique to measure the geometry of biomolecules on the nm scale by labeling them with pairs of fluorescent markers and measuring distance dependent non-radiative energy transfer between them. The remaining uncertainty in the application of the technique originates from the unknown orientation between the optical dipole moments of the fluorescent markers, especially when the molecule undergoes thermal fluctuations in physiological conditions. Recently we introduced a simulation based method for the interpretation of the fluorescence decay dynamics of the markers that allows us to retrieve both the average orientation and the extent of directional fluctuations of the involved dipole moments.

Meshless Methods for Photonic Crystal Modelling

Kamal Djidjeli, Marc Molinari, Simon Cox (Investigators), Neil O'Brien, Elizabeth Hart

We apply meshless methods to the problems of simulating photonic crystals. The meshless methods utilise compactly-supported radial basis functions (CSRBFs) and offer a promising alternative to the conventional plane-wave expansion method for calculating the band structure of photonic crystals.

Modelling Macro-Nutrient Release & Fate Resulting from Sediment Resuspension in Shelf Seas

Chris Wood

This study involves adapting a previously-published model to take into account the effect resuspension events (both natural and anthropogenic) may have on nutrient dynamics at the sediment-water interface, and hence produce better estimates for the total nutrient budgets for shelf seas.

Multi-objective design optimisation of coronary stents

Neil Bressloff, Georges Limbert (Investigators), Sanjay Pant

Stents are tubular type scaffolds that are deployed (using an inflatable balloon on a catheter), most commonly to recover the shape of narrowed (diseased) arterial segments. Despite the widespread clinical use of stents in cardiovascular intervention, the presence of such devices can cause adverse responses leading to fatality or to the need for further treatment. The most common unwanted responses of inflammation are in-stent restenosis and thrombosis. Such adverse biological responses in a stented artery are influenced by many factors, including the design of the stent. This project aims at using multi-objective optimisation techniques to find an optimum family of coronary stents which are more resistant to the processes of in-stent restenosis (IR) and stent thrombosis (ST).

MXL Project

Mark Taylor, Junfen Shi (Investigators)

‘MXL’ is short for “Enhanced patient safety by computational Modelling from clinically available X-rays to minimise the risk of overload and instability for optimised function and Longevity”. This is an international EU-funded project which the Bioengineering Sciences Research Group at Southampton is involved in. For more information, visit http://www.m-x-l.eu

nano-CMOS

Mark Zwolinski (Investigator), Michael Merrett

Modelling random device variations within systems using nano-CMOS technologies.

Nmag finite difference

Hans Fangohr (Investigator), Dmitri Chernyshenko, Matteo Franchin, Massoud Najafi

The goal of this project is to extends the finite element based micromagnetic simulation tool Nmag by the finite difference based extension Nmagfd and so to get an simulation tool where the user can easily switch between the used discretization method.

Nonlinear Optical Pulse Propagation

Peter Horak, Francesco Poletti (Investigators)

The work is concerned with the propagation of high-power short-pulse propagation in microstructured fibres or waveguides. Dispersion properties and optical nonlinearities are exploited for pulse shaping techniques in space, time, and frequency. Investigated microstructures include silica or soft-glass templates, gas-filled capillaries, and semiconductor-filled fibres, and optical wavelengths range from the X-ray to the mid-infrared regime.

OCCASION: Overcoming Capacity Constraints - A Simulation Integrated with Optimisation for Nodes

Tolga Bektas (Investigator)

OCCASION is a collaboration between TRG and the Schools of Mathematics and Management. The project's objective is to identify and investigate innovative methods of increasing the capacity of nodes (i.e. junctions and stations) on the railway network, without substantial investment in additional infrastructure. To this end, a state-of-the-art review of recent and ongoing work in this area will be conducted, followed by the development of tools to (i) assess existing levels of capacity utilisation at nodes, and (ii) investigate options for re-routeing and re-scheduling trains, with a view to reducing capacity utilisation levels. These tools will be used in combination to develop solutions delivering reduced levels of capacity utilisation, and thus increases in capacity and/or service reliability. Incremental changes to existing railway technologies (e.g. improved points) and operating practice (e.g. relaxations of the Rules of the Plan) will be investigated, as will concepts from other modes (e.g. road and air transport) and sectors (e.g. production scheduling).

Operational Simulation of the Solent Search-and-Rescue environment

James Scanlan, Kenji Takeda, Hans Fangohr (Investigators), Ben Schumann

This project aims to identify useful metrics for a proposed Search-and-Rescue UAV and test it virtually in a realistic environment.

Partial Discharge Signal Extraction using Spectral Methods

Condition based maintenance of 3 phase belted cables is increasing in demand as asset lifetimes approach their end...(some more info on the relevance of analysing belted cables).
The application of hard thresholding methods for Partial Discharge [PD] data extraction in high noise-ratio signals may produce series of PD trains with missing events and conversely correct events interlinked with noise. Subsequent analysis of the series will potentially be fraudulent and may lead to inaccurate conclusions. In this work, spectral methods for PD data extraction from very noisy environments are presented and a previously derived PD source classification technique is employed. Tests were conducted on a large set of wide bandwidth field data from three phase belted cables placed around London (UK) and Cyprus, and compared with accelerated aging experimental data. In the search of deeper insights, the deterministic origins of the produced sets of spike trains are explored and some linear and non-linear characteristics derived.

Prediction of orifice flow flooding rates through generic orifices

Dominic Hudson, Ming-yi Tan (Investigators), Christian Wood, Adam Sobey

This presearch concentrates on the modelling of compartment flooding rates following the occurrence of damage in a ship's side shell. Typical state of the art flooding models use Torricelli’s formula to calculate flooding rates using a constant co-efficient of discharge (Cd). Based on Bernoulli’s theorem, turbulence and viscosity effects are not included using a Cd independent of damage shape or size. Previous work indicates that this assumption over-simplifies the problem to an extent where the flooding rates used for calculation are in error. This project will use CFD validated by experiment to calculate flooding rates for a large number of cases from which a 'krigged' response surface will be generated. Validity of the subsequent response surface will be interrogated.

Real-time CFD for helicopter flight simulation

Kenji Takeda (Investigator), James Kenny

Project aims to show how real-time computational fluid dynamics (CFD) could be used to improve the realism of helicopter flight simulators.

Simulating Hydro-geomorphic Changes in European Climate Hotspots

John Dearing (Investigator), Ying Wang

This project will simulate the behaviour of hydro-geomorphological processes in a fluvial system over decadal timescales is an important basis for research on catchment environmental management, especially with regards climate changes and human impacts on fluvial system.

Simulation modelling of habitat permeability for mammalian wildlife

Patrick Doncaster, Jason Noble (Investigators), Angela Watkins

Using and integrating least-cost models and agent-based simulations to explore the way in which mammals interact with, and hence move, through fragmented landscapes.

Simulations investigating droplet diameter-charge models, for predicting electrostatically atomized dielectric liquid spray chracteristics

Gabriel Amine-Eddine (Investigator), John Shrimpton

Liquid sprays are atomized using electrostatic methods in many scientific, industrial and engineering applications. Due to jet and droplet breakup mechanisms, these spray plumes contain a range of drop diameters with a range of droplet charges. As a result of space charge repulsion forces between droplets, a wide range of inertial characteristics can be observed. Using an transient charged spray CFD code we have performed simulations to investigate charge-diameter relationship models for predicting the dynamics of poly-disperse and electrostatically charged hydrocarbon sprays.

A numerical procedure has been developed to allow for easy simulation and modelling of electrostatically atomised liquid sprays.

Currently work is underway, to extend the methodology developed in this study towards high-pressure sprays, where secondary atomization plays a dominant role within the spray dynamics and subsequent performance of the spray itself.

Spatial Mobility in the Formation of Agent-Based Economic Networks

Antonella Ianni, Seth Bullock (Investigators), Camillia Zedan

An investigation into the effect of spatial mobility on endogenous economic network formation.

Statistical model of the knee

Mark Taylor (Investigator), Francis Galloway, Prasanth Nair

Development of methods for large scale computational testing of a tibial tray incorporating inter-patient variability.

Tag based transcriptome analysis of gene expression in a promising green algae

Richard Edwards (Investigator), Andreas Johansson

We use SuperSAGE in combination with next-generation sequencing to compare differences in gene expression between selected mutants and the wild type of a green algae. The data in the form of millions of 26 bp tags representing short stretches of expressed genes, will be analysed to find patterns of variation in gene expression under different conditions.

The application of next-generation sequencing to unresolved familial disease

Andrew Collins, Sarah Ennis (Investigators), Jane Gibson, Reuben Pengelly

Next-generation sequencing (NGS) allows us to sequence individual patients cost-effectively, allowing us to enter a new era of genomic medicine. The level of genetic detail that we can access through these methods is unprecedented making it suitable for clinical molecular diagnostics.

The Endogenous Formation of Economic Networks

Antonella Ianni, Seth Bullock (Investigators), Camillia Zedan

An investigation into endogenous network formation using a simple agent-based approach.

The Role of Information in Price Discovery

Antonella Ianni, Seth Bullock (Investigators), Camillia Zedan

The recent economic crisis has highlighted a continued vulnerability and lack of understanding in the financial markets. In order to overcome this, many believe that current market models must be improved. Recently, a trend towards agent-based modelling has emerged. Viewing the economy as a complex system is beginning to be seen as key to explaining certain market characteristics that were originally considered anomalies.

One of the fundamental assumptions in economics is that of information efficiency: that the price of a stock reflects its worth, that all possible information about a security is publicly known, and that any changes to price take place instantaneously. In reality, however, this is not the case.

This project considers the use of agents in modelling economic systems and demonstrates the effect of information levels on price discovery using a simple market simulation.

Today's Computation Enabling Tomorrow's Seamless Communication

Lajos Hanzo (Investigator), Varghese Thomas

Radio Over Fibre (ROF) is a communication technique that aims to gainfully amalgamate the benefits of optical and wireless communication, while keeping the system cost low. This technique would support the next generation of wireless services.

Traveling and movement during European Late Prehistory

Patricia Murrieta Flores

This project has as main purpose to investigate through spatial analysis and computational modelling the variables and factors that influenced how humans traveled during prehistoric times.
One of the principal objectives will be to clarify the role that certain landscape elements (i.e megalithic monuments) played in terrestrial navigation and territorial definition.

This project is supported by CONACYT (Mexico) as a doctoral research by Patricia Murrieta-Flores under the supervision of Dr. David Wheatley (University of Southampton) and Dr. Leonardo Garcia Sanjuan (University of Seville, Spain). It also counts with the collaboration of Dr. Dimitrij Mlekuz (Gent University, Belgium).

Validation of a spatial-temporal soil water movement and plant water uptake model

Tiina Roose, Sevil Payvandi (Investigators), James Heppell

We develop a model that estimates the water saturation level within the soil at different depths, and the uptake of water by the root system. Data from Smethurst et al (2012) is used to validate our model and obtain a fully calibrated system for plant water uptake. When compared quantitatively to other models such as CROPWAT, our model achieves a better fit to the experimental data because of the simpler, first, second and third order terms present in the boundary condition, as opposed to complicated non-linear functions.

Whisky Code

Ian Hawke (Investigator)

A 3D finite volume code for simulating compact relativistic hydrodynamics.

Wind direction effects on urban flows

Zheng-Tong Xie, Ian Castro (Investigators), Jean Claus

Numerical simulations of turbulent air flow are conducted on Iridis to investigate the effects of different wind directions on the flow within and above an urban-like canopy.

µ-VIS Computed Tomography Centre

Ian Sinclair, Richard Boardman, Dmitry Grinev, Philipp Thurner, Simon Cox, Jeremy Frey, Mark Spearing, Kenji Takeda (Investigators)

A dedicated centre for computed tomography (CT) at Southampton, providing complete support for 3D imaging science, serving Engineering, Biomedical, Environmental and Archaeological Sciences. The centre encompasses five complementary scanning systems supporting resolutions down to 200nm and imaging volumes in excess of one metre: from a matchstick to a tree trunk, from an ant's wing to a gas turbine blade.