Ultra-Clean-Air (UCA) operating theatres aim to minimise surgical instrument contamination and wound infection through high flow rates of ultra-clean air, reducing the presence of Microbe Carrying Particles (MCPs). This study investigates the airflow patterns and ventilation characteristics of a UCA operating theatre (OT) under standard ventilation system operating conditions, considering both empty and partially occupied scenarios. Utilising a precise computational model, quasi-Direct Numerical Simulations (qDNS) were conducted to delineate flow velocity profiles, energy spectra, distributions of turbulent kinetic energy, energy dissipation rate, local Kolmogorov scales, and pressure-based coherent structures. These results were also complemented by a tracer gas decay analysis following ASHRAE standard guidelines. Simulations showed that contrary to the intended laminar regime, the OT’s geometry inherently fosters a predominantly turbulent airflow, sustained until evacuation through the exhaust vents, and facilitating recirculation zones irrespective of occupancy level. Notably, the occupied scenario demonstrated superior ventilation efficiency, a phenomenon attributed to enhanced kinetic energy induced by the additional obstructions. The findings underscore the critical role of UCA-OT design in mitigating MCP dissemination, highlighting the potential to augment the design to optimise airflow across a broader theatre spectrum, thereby diminishing recirculation zones and consequently reducing the propensity for Surgical Site Infections (SSIs). The study advocates for design refinements to harness the turbulent dynamics beneficially, steering towards a safer surgical environment.

LAMMPS is a powerful simulator originally developed for molecular dynamics that, today, also accounts for other particle-based algorithms such as DEM, SPH, or Peridynamics. The versatility of this software is further enhanced by the fact that it is open-source and modifiable by users. This property suits particularly well Discrete Multiphysics and hybrid models that combine multiple particle methods in the same simulation. Modifying LAMMPS can be challenging for researchers with little coding experience. The available material explaining how to modify LAMMPS is either too basic or too advanced for the average researcher. In this work, we provide several examples, with increasing level of complexity, suitable for researchers and practitioners in physics and engineering, who are familiar with coding without been experts. For each feature, step by step instructions for implementing them in LAMMPS are shown to allow researchers to easily follow the procedure and compile a new version of the code. The aim is to fill a gap in the literature with particular reference to the scientific community that uses particle methods for (discrete) multiphysics.

The presence of spherical solid particles immersed in an incompressible turbulent flow was numerically investigated from the perspective of the particle mass fraction (PMF or ${\varphi }_m$), a measure of the particle-to-fluid mass ratio. Although a number of different changes have been reported to be obtained by the presence of solid particles in incompressible turbulent flows, the present study reports the findings of varying ${\varphi }_m$ in the the turbulent behaviour of the flow, including aspects such as: turbulent statistics, skin-friction coefficient, and the general dynamics of a particle-laden flow. For this purpose, a particle-laden turbulent channel flow transporting solid particles at three different friction Reynolds numbers, namely $R{e}_{\tau }=180$, 365, and 950, with a fixed particle volume fraction of ${\varphi }_v={10}^{-3}$, was employed as conceptual flow model and simulated using large eddy simulations. The value adopted for ${\varphi }_v$ allowed the use of a two-way coupling approach between the particles and the flow or carrier phase. Three different values of ${\varphi }_m$ were explored in this work ${\varphi }_m\approx 1,2.96$, and $12.4$. Assessment of the effect of ${\varphi }_m$ was performed by examining changes of mean velocity profiles, velocity fluctuation profiles, and a number of other relevant turbulence statistics. Our results show that attenuation of turbulence activity of the carrier phase is attained, and that such attenuation increases with ${\varphi }_m$ at fixed Reynolds numbers and ${\varphi }_v$. For the smallest Reynolds number case considered, flows carrying particles with higher ${\varphi }_m$ exhibited lower energy requirements to sustain constant fluid mass flow rate conditions. By examining the flow velocity field, as well as instantaneous velocity components contours, it is shown that the attenuation acts even on the largest scales of the flow dynamics, and not only at the smaller levels. These findings reinforce the concept of a selective stabilising effect induced by the solid particles, particularly enhanced by high values of ${\varphi }_m$, which could eventually be exploited for improvement of energetic efficiency of piping or equivalent particles transport systems.

In this work, we present a mechanobiochemical model for two-dimensional cell migration which couples mechanical properties of the cell cytosol with biochemical processes taking place near or on the cell plasma membrane. The modelling approach is based on a recently developed mathematical formalism of evolving bulk- surface partial differential equations of reaction-diffusion type. We solve these equations using finite element methods within a moving-mesh framework derived from the weak formulation of the evolving bulk-surface PDEs. In the present work, the cell cytosol interior (bulk) dynamics are coupled to the cell membrane (surface) dynamics through non-homogeneous Dirichlet boundary conditions. The modelling approach exhibits both directed cell migration in response to chemical cues as well as spontaneous migration in the absence of such cues. Asa by-product, the approach shows fundamental characteristics associated with single cell migration such as: (i) cytosolic and membrane polarisation, (ii) actin dependent protrusions, and (iii) continuous shape deformation of the cell during migration. Cell migration is an ubiquitous process in life that is mainly triggered by the dynamics of the actin cytoskeleton and therefore is driven by both mechanical and biochemical processes. It is a multistep process essential for mammalian organisms and is closely linked to a vast diversity of processes; from embryonic development to cancer invasion. Experimental, theoretical and computational studies have been key to elucidate the mechanisms underlying cell migration. On one hand, rapid advances in experimental techniques allow for detailed experimental measurements of cell migration pathways, while, on the other, computational approaches allow for the modelling, analysis and understanding of such observations. The bulk-surface mechanobiochemical modelling approach presented in this work, set premises to study single cell migration through complex non-isotropic environments in two- and three-space dimensions.

Multiphysics models have become a key tool in understanding the way different phenomenon are related in bone remodeling and various approaches have been proposed, yet, to the best of the author’s knowledge there is no model able to link a cell population model with a mechanical stimulus model using a discrete approach, which allows for an easy implementation. This article couples two classical models, the cell population model from Komarova and the Nackenhorst model in a 2D domain, where correlations between the mechanical loading and the cell population dynamics can be established, furthermore the effect of different paracrine and autocrine regulators is seen on the overall density of a portion of trabecular bone. A discretization is performed using frame 1D finite elements, representing the trabecular structure. The Nackenhorst model is implemented by using the finite element method to calculate the strain energy as the main mechanical stimulus that determines the bone mass density evolution in time. This density is normalized to be added to the bone mass percentage proposed by the Komarova model, where coupling terms have been added as well that guarantee a stable response. In the simulations, the equations were solved employing the finite element method with a user subroutine implemented in ABAQUS (2017) and by applying a direct formulation. The methodology presented can model the cell dynamics occurring in bone remodelling in accordance with the asynchronous nature of this process, yet allowing to differentiate zones with higher density, the main trabecular groups are obtained for the proximal femur. Finally, the model is tested in pathological cases, such as osteoporosis and osteopetrosis, yielding results similar to the pathology behavior. Furthermore, the discrete modelling technique is shown to be of use in this particular application.