Search

This paper presents the first classical Convolutional Neural Network (CNN) that can be applied directly to data from unstructured finite element meshes or control volume grids. CNNs have been hugely influential in the areas of image classification and image compression, both of which typically deal with data on structured grids. Unstructured meshes are frequently used to solve partial differential equations and are particularly suitable for problems that require the mesh to conform to complex geometries or for problems that require variable mesh resolution. Central to our approach are space-filling curves, which traverse the nodes or cells of a mesh tracing out a path that is as short as possible (in terms of numbers of edges) and that visits each node or cell exactly once. The space-filling curves (SFCs) are used to find an ordering of the nodes or cells that can transform multi-dimensional solutions on unstructured meshes into a one-dimensional (1D) representation, to which 1D convolutional layers can then be applied. Although developed in two dimensions, the approach is applicable to higher dimensional problems. To demonstrate the approach, the network we choose is a convolutional autoencoder (CAE), although other types of CNN could be used. The approach is tested by applying CAEs to data sets that have been reordered with a space-filling curve. Sparse layers are used at the input and output of the autoencoder, and the use of multiple SFCs is explored. We compare the accuracy of the SFC-based CAE with that of a classical CAE applied to two idealised problems on structured meshes, and then apply the approach to solutions of flow past a cylinder obtained using the finite-element method and an unstructured mesh., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Roadside concentrations of harmful pollutants such as NO x are highly variable in both space and time. This is rarely considered when assessing pedestrian and cyclist exposures. We aim to fully describe the spatio-temporal variability of exposures of pedestrians and cyclists travelling along a road at high resolution. We evaluate the value added of high spatio-temporal resolution compared to high spatial resolution only. We also compare high resolution vehicle emissions modelling to using a constant volume source. We highlight conditions of peak exposures, and discuss implications for health impact assessments. Using the large eddy simulation code Fluidity we simulate NO x concentrations at a resolution of 2 m and 1 s along a 350 m road segment in a complex real-world street geometry including an intersection and bus stops. We then simulate pedestrian and cyclist journeys for different routes and departure times. For the high spatio-temporal method, the standard deviation in 1 s concentration experienced by pedestrians (50.9 μg.m -3 ) is nearly three times greater than that predicted by the high-spatial only (17.5 μg.m -3 ) or constant volume source (17.6 μg.m -3 ) methods. This exposure is characterised by low concentrations punctuated by short duration, peak exposures which elevate the mean exposure and are not captured by the other two methods. We also find that the mean exposure of cyclists on the road (31.8 μg.m -3 ) is significantly greater than that of cyclists on a roadside path (25.6 μg.m -3 ) and that of pedestrians on a sidewalk (17.6 μg.m -3 ). We conclude that ignoring high resolution temporal air pollution variability experienced at the breathing time scale can lead to a mischaracterization of pedestrian and cyclist exposures, and therefore also potentially the harm caused. High resolution methods reveal that peaks, and hence mean exposures, can be meaningfully reduced by avoiding hyper-local hotspots such as bus stops and junctions., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

Vaccination strategy is crucial in fighting the COVID-19 pandemic. Since the supply is still limited in many countries, contact network-based interventions can be most powerful to set an efficient strategy by identifying high-risk individuals or communities. However, due to the high dimension, only partial and noisy network information can be available in practice, especially for dynamic systems where contact networks are highly time-variant. Furthermore, the numerous mutations of SARS-CoV-2 have a significant impact on the infectious probability, requiring real-time network updating algorithms. In this study, we propose a sequential network updating approach based on data assimilation techniques to combine different sources of temporal information. We then prioritise the individuals with high-degree or high-centrality, obtained from assimilated networks, for vaccination. The assimilation-based approach is compared with the standard method (based on partially observed networks) and a random selection strategy in terms of vaccination effectiveness in a SIR model. The numerical comparison is first carried out using real-world face-to-face dynamic networks collected in a high school, followed by sequential multi-layer networks generated relying on the Barabasi-Albert model emulating large-scale social networks with several communities., Competing Interests: Conflict of interestThe authors declare that there is no conflict of interest., (© The Author(s) 2023.)

Underground railway systems are recognised spaces of increased personal pollution exposure. We studied the number-size distribution and physico-chemical characteristics of ultrafine (PM 0.1 ), fine (PM 0.1 - 2.5 ) and coarse (PM 2.5 - 10 ) particles collected on a London underground platform. Particle number concentrations gradually increased throughout the day, with a maximum concentration between 18:00 h and 21:00 h (local time). There was a maximum decrease in mass for the PM 2.5 , PM 2.5 - 10 and black carbon of 3.9, 4.5 and ~ 21-times, respectively, between operable (OpHrs) and non-operable (N-OpHrs) hours. Average PM 10 (52 μg m -3 ) and PM 2.5 (34 μg m -3 ) concentrations over the full data showed levels above the World Health Organization Air Quality Guidelines. Respiratory deposition doses of particle number and mass concentrations were calculated and found to be two- and four-times higher during OpHrs compared with N-OpHrs, reflecting events such as train arrival/departure during OpHrs. Organic compounds were composed of aromatic hydrocarbons and polycyclic aromatic hydrocarbons (PAHs) which are known to be harmful to health. Specific ratios of PAHs were identified for underground transport that may reflect an interaction between PAHs and fine particles. Scanning transmission electron microscopy (STEM) chemical maps of fine and ultrafine fractions show they are composed of Fe and O in the form of magnetite and nanosized mixtures of metals including Cr, Al, Ni and Mn. These findings, and the low air change rate (0.17 to 0.46 h -1 ), highlight the need to improve the ventilation conditions., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

We propose a novel use of generative adversarial networks (GANs) (i) to make predictions in time (PredGAN) and (ii) to assimilate measurements (DA-PredGAN). In the latter case, we take advantage of the natural adjoint-like properties of generative models and the ability to simulate forwards and backwards in time. GANs have received much attention recently, after achieving excellent results for their generation of realistic-looking images. We wish to explore how this property translates to new applications in computational modelling and to exploit the adjoint-like properties for efficient data assimilation. We apply these methods to a compartmental model in epidemiology that is able to model space and time variations, and that mimics the spread of COVID-19 in an idealised town. To do this, the GAN is set within a reduced-order model, which uses a low-dimensional space for the spatial distribution of the simulation states. Then the GAN learns the evolution of the low-dimensional states over time. The results show that the proposed methods can accurately predict the evolution of the high-fidelity numerical simulation, and can efficiently assimilate observed data and determine the corresponding model parameters., Competing Interests: Conflict of interestThe authors declare that they have no conflict of interest., (© The Author(s) 2022.)

We propose the use of reduced order modeling (ROM) to reduce the computational cost and improve the convergence rate of nonlinear solvers of full order models (FOM) for solving partial differential equations. In this study, a novel ROM-assisted approach is developed to improve the computational efficiency of FOM nonlinear solvers by using ROM's prediction as an initial guess. We hypothesize that the nonlinear solver will take fewer steps to the converged solutions with an initial guess that is closer to the real solutions. To evaluate our approach, four physical problems with varying degrees of nonlinearity in flow and mechanics have been tested: Richards' equation of water flow in heterogeneous porous media, a contact problem in a hyperelastic material, two-phase flow in layered porous media, and fracture propagation in a homogeneous material. Overall, our approach maintains the FOM's accuracy while speeding up nonlinear solver by 18-73% (through suitable ROM-assisted FOMs). More importantly, the proximity of ROM's prediction to the solution space leads to the improved convergence of FOMs that would have otherwise diverged with default initial guesses. We demonstrate that the ROM's accuracy can impact the computational efficiency with more accurate ROM solutions, resulting in a better cost reduction. We also illustrate that this approach could be used in many FOM discretizations (e.g., finite volume, finite element, or a combination of those). Since our ROMs are data-driven and non-intrusive, the proposed procedure can easily lend itself to any nonlinear physics-based problem., (© 2022. The Author(s).)

The outbreak of the coronavirus disease 2019 (COVID-19) has now spread throughout the globe infecting over 150 million people and causing the death of over 3.2 million people. Thus, there is an urgent need to study the dynamics of epidemiological models to gain a better understanding of how such diseases spread. While epidemiological models can be computationally expensive, recent advances in machine learning techniques have given rise to neural networks with the ability to learn and predict complex dynamics at reduced computational costs. Here we introduce two digital twins of a SEIRS model applied to an idealised town. The SEIRS model has been modified to take account of spatial variation and, where possible, the model parameters are based on official virus spreading data from the UK. We compare predictions from one digital twin based on a data-corrected Bidirectional Long Short-Term Memory network with predictions from another digital twin based on a predictive Generative Adversarial Network. The predictions given by these two frameworks are accurate when compared to the original SEIRS model data. Additionally, these frameworks are data-agnostic and could be applied to towns, idealised or real, in the UK or in other countries. Also, more compartments could be included in the SEIRS model, in order to study more realistic epidemiological behaviour., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2021 Elsevier B.V. All rights reserved.)

The paper presents a novel method for numerically modelling fluid-structure interactions. The method consists of solving the fluid-dynamics equations on an extended domain, where the computational mesh covers both fluid and solid structures. The fluid and solid velocities are relaxed to one another through a penalty force. The latter acts on a thin shell surrounding the solid structures. Additionally, the shell is represented on the extended domain by a non-zero shell-concentration field, which is obtained by conservatively mapping the shell mesh onto the extended mesh. The paper outlines the theory underpinning this novel method, referred to as the immersed-shell approach. It also shows how the coupling between a fluid- and a structural-dynamics solver is achieved. At this stage, results are shown for cases of fundamental interest., (© 2015 The Author(s) Published by the Royal Society. All rights reserved.)

Research into the use of unstructured mesh methods in oceanography has been growing steadily over the past decade. The advantages of this approach for domain representation and non-uniform resolution are clear. However, a number of issues remain, in particular those related to the computational cost of models produced using unstructured mesh methods compared with their structured mesh counterparts. Mesh adaptivity represents an important means to improve the competitiveness of unstructured mesh models, where high resolution is only used when and where necessary. In this paper, an optimization-based approach to mesh adaptivity is described where emphasis is placed on capturing anisotropic solution characteristics. Comparisons are made between the results obtained with uniform isotropic resolution, isotropic adaptive resolution and fully anisotropic adaptive resolution.