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.

Adding social ties to the Schelling model

Seth Bullock, Sally Brailsford (Investigators), Elisabeth zu-Erbach-Schoenberg

The Schelling model is an abstract model for segregation in
a spatially arranged population. We extended the traditional model by the addition of a dynamic social network. The social network influences the spatial dynamics of agents moving on the grid by changing the agents’ evaluation of their neighbourhood. In turn, the spatial arrangement influences the change of the social network.

Advanced modelling for two-phase reacting flow

Edward Richardson (Investigator)

Engine designers want computer programs to help them invent ways to use less fuel and produce less pollution. This research aims to provide an accurate and practical model for the injection and combustion of liquid fuel blends.

Aerofoil noise

Richard Sandberg (Investigator)

High-performance computing is used to identify noise sources on aerofoils.

Antimicrobial Peptide and E. coli Membrane Interactions

Syma Khalid (Investigator), Thomas Piggot, Nils Berglund

Antimicrobial peptides (AMPs) are known to disrupt the membranes of bacterial cells such as E. coli. I work on investigating the nature of these interactions using molecular dynamics (MD) simulations.

Automated selection of suitable atmospheric calibration sites for satellite imagery

Robin Wilson (Investigator)

Ground calibration targets (GCTs) play a vital role in atmospheric correction of satellite sensor data in the optical region, but selecting suitable targets is a subjective and time- consuming task. This project is developing methods to automatically select suitable GCTs, using a combination of remotely sensed multispectral and topographic data.

B-meson coupling with relativistic heavy quarks

Jonathan Flynn (Investigator), Patrick Fritzsch, Dirk Broemmel

We non-perturbatively compute the coupling between B* and B pi meson states relying on relativistic heavy quarks and domain wall light fermions. The coupling is of importance for an effective description of hadronic heavy meson decays.

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.

BioSimGrid

Jonathan Essex, Hans Fangohr (Investigators), Richard Boardman, Syma Khalid, Steven Johnston

The aim of the BioSimGrid project is to make the results of large-scale computer simulations of biomolecules more accessible to the biological community. Such simulations of the motions of proteins are a key component in understanding how the structure of a protein is related to its dynamic function.

Cloud Computing for Planetary Defense

Hugh Lewis, Kenji Takeda (Investigators), Steven Johnston

We demonstrate how a cloud-based computing architecture can be used for planetary defense and space situational awareness (SSA). We show how utility compute can facilitate both a financially economical and highly scalable solution for space debris and near-earth object impact analysis. As we improve our ability to track smaller space objects, and satellite collisions occur, the volume of objects being tracked vastly increases, increasing computational demands. Propagating trajectories and calculating conjunctions becomes increasingly time critical, thus requiring an architecture which can scale with demand. The extension of this to tackle the problem of a future near-earth object impact is discussed, and how cloud computing can play a key role in this civilisation-threatening scenario.

Complexity in Modelling Electric Marine Propulsive Devices

Suleiman Sharkh, Neil Bressloff, Hans Fangohr (Investigators), Aleksander Dubas

This project involves the simulation of turbulent flow around a marine rim-driven thruster and the complex interaction of flow features involved through computational fluid dynamics.

Computational electromagnetic modelling of 3D photonic structures

Marc Molinari, Darren Bagnall, Simon Cox (Investigators), Asa Asadollahbaik, Elizabeth Hart

Nano-structured materials can provide very specific and often very special optical effects which can be exploited for a large range of optical applications including wavelength filters, LEDs, micro-lasers, HDTV, solar-cell coatings, optical high-Q fibres, diffraction gratings, polarisation devices, optical switches, etc. This research in “Computational Electromagnetic Modelling of 3D Photonic Structures” aims to address the need for accurate and fast three-dimensional modelling, simulation and analysis processes in the photonics industry. A FEM/FDTD software suite will be developed to simulate Maxwell’s field equations and thin-film quantum effects (plasmons) in the visible and near-infrared EM frequency spectrum. The results obtained from running the software on suitable compute clusters will then be compared to the analysis results of experimentally manufactured materials. We will investigate structures occurring in nature such as iridescent butterfly wings, white/black reflecting beetle shells, etc., and aim to optimise artificially designed structures with periodic, quasi-periodic and random configurations.

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.

Designer 3D Magnetic Mesostructures

Hans Fangohr (Investigator), Matteo Franchin, Andreas Knittel

A new electrodeposition self-assembly method allows for the growth of well defined mesostructures. This project's aim is to use this method in order to fabricate supraconducting and ferromagnetic mesostructures. Numerical methods based on well-established models are used in order to characterise the grown structures.

Development of a novel Navier-Stokes solver (HiPSTAR)

Richard Sandberg (Investigator)

Development of a highly efficient Navier-Stokes solver for HPC.

Dynamag: computational magnonics

Hans Fangohr, Atul Bhaskar (Investigators), Matteo Franchin, Andreas Knittel

Analytical treatment of long range magneto-dipole interactions is a bottle-neck of magnonics and more generally of the theory of spin waves in non-uniform media. This project develops a theoretical framework for analysis of magnonic phenomena in magnetic nano-structures, including isolated nano-elements, arrays of those, and extended magnonic crystals. The DYNAMAG project is funded by the EU FP7 and the DST of India.

Dynamics of interacting magnetic nanoparticles

Thomas Fischbacher (Investigator), Maximilian Albert

The project aims at extending the micromagnetic simulation framework 'nmag' developed at the University of Southampton to enable it to handle dynamic geometries. The extended framework will then be used to study systems such as interacting magnetic nanoparticles.

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?”

Fluid Structure Interactions of Yacht Sails

Stephen Turnock (Investigator), Daniele Trimarchi

The research is the main subject of the PhD topic. It regards the application of fluid structure interaction techniques to the domain of yacht sails simulation

Fracturing of small social networks

Seth Bullock, Sally Brailsford (Investigators), Elisabeth zu-Erbach-Schoenberg

A connected social network is a very important factor for the success of groups and organisations. We investigate which factors make a group more resistant to the effects of disagreements which commonly happen in small social networks.

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.

Gravitational waves from neutron stars

Ian Hawke (Investigator)

Gravitational waves, once detected, will give information about the extremes of space and time. Compact objects such as neutron stars are perfect locations for generating such waves.

Hadronic structure on the computer

Jonathan Flynn (Investigator), Dirk Broemmel, Thomas Rae, Ben Samways

In experiments at the Large Hadron Collider (LHC) at CERN, Geneva, the interactions that occur between the colliding particles (protons in this case) can be factorised into a simple scattering between two constituent particles, called quarks, followed by a hadronisation process, which describes the dynamics of forming the bound proton states. Quarks are particles within the proton that bind to form composite particles (hadrons) such as a proton. The scattering process can be computed relatively easily, but hadronisation is intrinsically non-perturbative and hard to calculate. Lattice QCD (computer simulation of QCD on a discrete space-time lattice) provides our only known first-principles and systematically-improvable method to address problems like hadronisation. This project uses Iridis to extract parton distribution amplitudes which are experimentally inaccessible, but needed to describe the quark structure of hadrons.

How far can we stretch the MARTINI?

Syma Khalid (Investigator), Ric Gillams

To date, coarse-grained lipid models have generally been parameterised to ensure the correct prediction of structural properties of membranes, such as the area per lipid and the bilayer thickness. The work described here explores the extent to which coarse-grained models are able to predict correctly bulk properties of lipids (phase behaviour) as well as the mechanical properties, such as lateral pressure profiles and stored elastic stress in bilayers. Such an evaluation is crucial for understanding the predictive capabilities of coarse-grained models.

Hybrid quantum and classical free energy methods in computational drug optimisation

Jonathan Essex, Chris-Kriton Skylaris (Investigators), Christopher Cave-Ayland

This work is based around the application of thermodynamics and quantum mechanics to the field of computational drug design and optimisation. Through the application of these theories the calculation of the physical properties of drug-like molecules is possible and hence some predictive power for their pharmaceutical activity in vivo can be obtained.

Hybrid RANS/LES methods

Richard Sandberg (Investigator), Markus Weinmann

Novel hybrid RANS/LES methods are developed for more accurate and efficient simulation of flow over complex geometries.

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.

Is fine-scale turbulence universal?

Richard Sandberg (Investigator)

Complementary numerical simulations and experiments of various canonical flows will try to answer the question whether fine-scale turbulence is universal.

Jet noise

Richard Sandberg (Investigator), Neil Sandham

Direct numerical simulations are used to investigate jet noise.

Kaon to two pion decays in lattice QCD

Jonathan Flynn (Investigator), Elaine Goode, Dirk Broemmel

We calculate kaon decay amplitudes on the lattice so we may compare the Standard Model to experiment.

Laminar to Turbulent Transition in Hypersonic Flows

Neil Sandham, Heinrich Luedeke

Understanding of laminar to turbulent transition in hypersonic boundary-layer flows is crucial for re-entry vehicle design and optimization. The boundary-layer state directly affects the temperatures on the vehicle surface and its viscous drag. Therefore transition has to be considered to correctly compensate for drag and to properly design the thermal protection system.
For the proposed study, in order to obtain a clear understanding of the transition process, the configuration is kept as simple as possible by varying only a minimum number of parameters affecting transition on a simple test geometry such as a swept ramp at different sweep angles. To investigate the influence of such sweep angles on the transition process in the hypersonic regime, Direct Numerical Simulations (DNS) of the turbulent flow field are carried out on the Iridis cluster.

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.

Molecular Fragments in Inhibitor Design

Jonathan Essex (Investigator), Michael Bodnarchuk

Fragment-Based Drug Discovery (FBDD) has emerged as an important tool in the drug discovery process. Instead of screening entire drug molecules, FBDD screens molecular fragments; constituents which make up drug molecules. A computational approach to identifying fragment binding is currently being sought which also yield binding free energy estimation.

Multi-scale simulations of bacterial outer-membrane proteins

Syma Khalid (Investigator), Jamie Parkin

Using Iridis to run multiple simulations, I aim to simulate the outer membrane proteins of Pseudomonas aeruginosa, using X-ray crystal structures of proteins only recently resolved by Bert van den Berg, University of Massachusetts. By modelling the proteins in a realistic P. aeruginosa outer membrane, I am to gain insight into the binding of these proteins to specific substrates and their function.

Multiscale modelling of biological membranes

Jonathan Essex (Investigator), Mario Orsi

Biological membranes are complex and fascinating systems, characterised by proteins floating in a sea of lipids. Biomembranes, besides being the fundamental structures employed by nature to encapsulate cells, play crucial roles in many phenomena indispensable for life, such as growth, energy storage, and in general information transduction via neural activity. In this project, we develop and apply multiscale computational models to simulate biological membranes and obtain molecular-level insights into fundamental structures and phenomena.

Multiscale Simulation of Cellular Calcium Signalling

Hans Fangohr, Jonathan Essex (Investigators), Dan Mason

Calcium ions play a vitally important role in signal transduction and are key to many cellular processes including muscle contraction and cell apoptosis (cell death). This importance has made calcium an active area in biomedical science and mathematical modelling.

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

Nmag - computational micromagnetics

Hans Fangohr, Thomas Fischbacher (Investigators), Matteo Franchin, Andreas Knittel, Maximilian Albert, Dmitri Chernyshenko, Massoud Najafi, Richard Boardman

Nmag is a micromagnetic simulation package based on the general purpose multi-physics library nsim. It is developed by the group of Hans Fangohr and Thomas Fischbacher in the School of Engineering Sciences at the University of Southampton and released under the GNU GPL.

Non-Perturbative Renormalisation on the Lattice

Jonathan Flynn (Investigator), Dirk Broemmel, Thomas Rae

In this project we compute renormalisation factors for various physical observables in a non-perturbative lattice framework. Renormalisation hereby arises due to a fundamental scale dependence of the physical processes.

OMSys Towards a system model of a bacterial outer membrane

Syma Khalid (Investigator)

Many bacteria have an outer membrane which is the interface between the cell and its environment. The components of this membrane are well studied at an individual level, but there is a need to model and understand the outer membrane as a whole. In this project we aim to develop such a model of a bacterial outer membrane, linking computer simulations of the component molecules through to a more "systems biology" approach to modelling the outer membrane as a whole. Such an approach to modelling an OM must be multi-scale i.e. it must embrace a number of levels ranging from atomistic level modelling of e.g. the component proteins through to higher level "agent-based" modelling of the interplay of multiple components within the outer membrane as a whole. The different levels of description will be integrated to enable predictive modelling in order to explore the roles of outer membrane changes in e.g. antibiotic resistance.

Population24/7: space-time specific population surface modelling

Samantha Cockings, David Martin, Samuel Leung (Investigators)

Project funded by Economic and Social Research Council to compute time-specific geographical representations of population distribution.

Probing the oligomeric state and interaction surface of Fukutin Transmembrane Domain in lipid bilayer via Molecular Dynamics simulations

Nils Berglund

Fukutin Transmembrane Domain (FK1TMD) is localised to the endoplasmic reticulum or Golgi Apparatus within the cell where it is believed to function as a glycosyltransferase. Its localisation within the cell is thought to be mediated by the interaction of its N-terminal transmembrane domain with the lipid bilayers surrounding these compartments, each of which possess a distinctive lipid composition. Studies have revealed that the N-terminal transmembrane domain of FK1TMD exists as dimer within dilauroylphosphatidylcholine bilayers and this interaction is driven by interactions between a characteristic TXXSS motif. Furthermore residues close to the N-terminus that have previously been shown to play a key role in the clustering of lipids are shown to play a key role in anchoring the protein in the membrane.

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.

Scalability of Energy Efficient Routing Algorithms in WSN

Davide Zilli

This project compares two broad classes of routing algorithms for Wireless Sensor Networks, message flooding and single path, by means of a simulation model. In particular, we want to understand how the two scale in terms of energy efficiency on large networks of sensors.

Self-Force and Black Hole Inspirals

Sam Dolan (Investigator)

We use IRIDIS to compute the self-force acting on a solar-mass black hole orbiting a supermassive black hole.

Separation of timescales in models of complex networks

Seth Bullock (Investigator), Elisabeth zu-Erbach-Schoenberg, Connor McCabe

In many real-world systems several processes act on the system state. The way these processes interact can have implications for the resulting system state. We investigate how separation of the timescales of two processes influences the system's equilibrium state.

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.

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.

Study of global instability in separated flows at high Mach number

Neil Sandham, Zhiwei Hu (Investigators), Kangping Zhang

Flow instability is observed when extending two-dimensional (2D) stable flow into three-dimensional (3D). Development of instability varies along different spanwise length. Thresholds are also discovered for the flow studied to become instable.

Supersonic axisymmetric wakes

Richard Sandberg (Investigator)

Direct numerical simulations are used to shed more light on structure formation and evolution in supersonic wakes.

The application of automated pattern metrics to surface moisture influences on modelled dune field development

Robin Wilson, Joanna Nield (Investigators)

Areas of sand dunes (known as dunefields) develop complex patterns over time. These are influenced by both the past and present environmental conditions, including surface moisture, vegetation distribution and human impact. This project develops a method of automated pattern analysis which allow the patterns produced by a large number of sand dune evolution simulations (performed using the DECAL model) to be quantified over time.

The effect of roughness upon turbulent supersonic flows

Neil Sandham (Investigator), Christopher Tyson

Understanding the interaction between surface roughness and supersonic air flows are crucial in the design of re-entry vehicles such as the space shuttle. Numerical simulations of these flows has been conducted in order to examine the relationship in order to attempt to achieve a much clearer understanding of the behaviour.

Transition to turbulence in high-speed boundary layers

Neil Sandham (Investigator), Nicola De Tullio

This work is focused on the numerical simulation of hypersonic transition to turbulence in boundary layers. We use direct numerical simulations of the Navier-Stokes equations to analyse the effects of different flow conditions and external disturbances on the transition process. The main objective is to gain insight into the different aspects of transition to turbulence at high speeds, which can lead to the design of new transition models and transition control techniques for high-speed flows.

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).

Water Molecules in Protein Binding Sites

Jonathan Essex (Investigator), Michael Bodnarchuk

Water molecules are commonplace in protein binding sites, although the true location of them can often be hard to predict from crystallographic methods. We are developing tools which enable the location and affinity of water molecules to be found.

Wave-based discontinuous Galerkin methods

Gwenael Gabard (Investigator), Greg Kennedy

Wave-based computational methods are developed to model sound propagation in moving inhomogeneous media.

Whisky Code

Ian Hawke (Investigator)

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

µ-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.