Search

A method for nonlinear material modeling and design using statistical learning is proposed to assist in the mechanical analysis of structural materials. Conventional computational homogenization schemes are proven to underperform in analyzing the complex nonlinear behavior of such microstructures with finite deformations. Also, the higher computational cost of the existing homogenization schemes inspires the inception of a data-driven multiscale computational homogenization scheme. In this paper, a statistical nonlinear homogenization scheme is discussed to mitigate these issues using the Gaussian Process Regression technique. A data-driven model is trained for different strain states of microscale unit cells. In the macroscale, nonlinear response of the macroscopic structure is analyzed, for which the stresses and material responses are predicted by the trained surrogate model.

Wind-fed models offer a unique way to form predictive models of the accretion flow surrounding Sagittarius A*. We present 3D wind-fed magnetohydrodynamic (MHD) and general relativistic magnetohydrodynamic (GRMHD) simulations spanning the entire dynamic range of accretion from parsec scales to the event horizon. We expand on previous work by including non-zero black hole spin and dynamically evolved electron thermodynamics. Initial conditions for these simulations are generated from simulations of the observed Wolf–Rayet stellar winds in the Galactic Centre. The resulting flow tends to be highly magnetized (β ≈ 2) with an ∼r−1 density profile independent of the strength of magnetic fields in the winds. Our simulations reach the magnetically arrested disc (MAD) state for some, but not all cases. In tilted flows, standard and normal evolution (SANE) jets tend to align with the angular momentum of the gas at large scales, even if that direction is perpendicular to the black hole spin axis. Conversely, MAD jets tend to align with the black hole spin axis. The gas angular momentum shows similar behaviour: SANE flows tend to only partially align while MAD flows tend to fully align. With a limited number of dynamical free parameters, our models can produce accretion rates, 230 GHz flux, and unresolved linear polarization fractions roughly consistent with observations for several choices of electron heating fraction. Absent another source of large-scale magnetic field, winds with a higher degree of magnetization (e.g. where the magnetic pressure is 1/100 of the ram pressure in the winds) may be required to get a sufficiently large rotation measure with consistent sign.

Virus infection dynamics is governed by the processes on multiple scales: on the whole organism level, tissue level, and intracellular level. In this paper, we develop a multi-scale multi-compartment model of HIV infection in a simplified setting and the computational methods for numerical realization of the model. The multiscale model describes the processes from various scales and of different nature (cell motility, virus diffusion, intracellular virus replication). Intracellular replication model is based on a Markov chain with time-inhomogeneous propensities that depend on the extracellular level of virions. Reaction diffusion equations used to model free virion diffusion in the lymphoid tissue have moving sources, which are determined by the positions of the infected cells (immune cell motility model) and the rate of virion secretion from them (intracellular model). Immune cell motility model parameterizes the intercellular interaction forces, friction and the stochastic force of active cell motility. Together, this allows for a proper description of the intracellular stochasticity that propagates across multiple scales. A hybrid discrete-continuous stochastic-deterministic algorithm for simulation of the multiscale model based on the uniformization Monte Carlo method is implemented.

Based on a multiscale approach, the in-home built Fortran code SIMFC allows for high-temperature fuel cell simulation from material intrinsic properties to system overall operation.

The healthcare industry is different from other industries–patient data are sensitive, their storage needs to be handled with care and in compliance with regulative, while prediction accuracy needs to be high. This fast expansion in medical image modalities and data collection leads to generation of so called “Big Data” which is time-consuming to be analyzed by medical experts. This paper provides an insight into the Big Data from the aspect of its role in multiscale modelling. Special attention is paid to the workflow, starting from medical image processing all the way to creation of personalized models and their analysis. A review of literature regarding Big Data in healthcare is provided and two proposed solutions are described–carotid artery ultrasound image processing and 3D reconstruction, and drug testing on personalized heart models. Related to the carotid artery ultrasound image processing, the starting point is ultrasound images, which are segmented using convolutional neural network U-net, while segmented masks were further used in 3D reconstruction of geometry. Related to the drug testing on personalized heart model, similar approach was proposed, images were used in creation of personalized 3D geometrical model that is used in computational modelling to determine pressure in the left ventricle before and after drug testing. All the aforementioned methodologies are complex, include Big Data analysis and should be performed using servers or high-performance computing. Future development of Big Data applications in healthcare domains offers a lot of potential due to new data standards, rapid development of research and technology, as well as strong government incentives.

Our study provides a global assessment of water balance components accounting for the uncertainty in globally available precipitation products. This assessment is carried out consistently using a multiscale modelling framework established over more than 6000 GRDC river basins at various spatial resolutions (daily time step, period 1990-2019). The framework is based on the mesoscale Hydrologic Model (mHM; [1,2,3]) equipped with the multiscale parameter regionalization (MPR) scheme. All basins share the same parameterization and are driven with four different state-of-the-art meteorological products: ERA5 reanalysis [4], MSWEP [5], and deterministic EM-Earth v2 [6]. Additionally, the hydrological simulations are benchmarked against E-OBS [7] over Europe and against locally interpolated 1km gridded rain gauge dataset over Germany.Our results show that EM-Earth clearly exhibits the best streamflow performance across North and South America and Asia with respect to the other products. The MSWEP is the best product in Africa, where the overall model’s performance is rather poor. In Australia, MSWEP and EM-Earth have comparable skills. In Europe, the differences get narrower, although slightly better performance is seen for EM-Earth, with a median performance of 0.5 of daily KGE (N basins=2000), mainly due to correcting the under-catch error of rain gauges which is not considered, e.g. in E-OBS. Furthermore, when we zoom into a subset consisting of medium-sized 200 German basins, the in-house high-resolution meteorologic product clearly overperformed all global products, mainly due to better captured temporal correlation and smaller biases. On the other hand, ERA5 leads to very strong positive biases over German basins using standard parameterization. Finally, our study contributes to discussions on objective quantification of the optimal spatial resolution of hydrological studies. [1] https://doi.org/10.5281/zenodo.5119952[2] https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2008WR007327 [3] https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2012WR012195[4] https://rmets.onlinelibrary.wiley.com/doi/10.1002/qj.3803[5] https://hess.copernicus.org/articles/21/589/2017/hess-21-589-2017.html[6] https://journals.ametsoc.org/view/journals/bams/103/4/BAMS-D-21-0106.1.xml[7] https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2017JD028200

Accurate designing of electrical machines made of customized materials, such as soft magnetic composites (SMC), is essential when aiming for higher energy efficiency and minimized amount of material per device. To simulate the electromagnetic behavior of SMC materials, multiscale methods are needed. In this paper, a novel multiscale approach using an EH EH EH-conforming formulation and multiharmonic method is presented and applied to SMC materials. Based on the simulation results, the prospects of the approach are discussed in detail.

The M4F project brings together the fusion and fission materials communities working on the prediction of radiation damage production and evolution and their effects on the mechanical behaviour of irradiated ferritic/martensitic (F/M) steels. It is a multidisciplinary project in which several different experimental and computational materials science tools are integrated to understand and model the complex phenomena associated with the formation and evolution of irradiation induced defects and their effects on the macroscopic behaviour of the target materials. In particular the project focuses on two specific aspects: (1) To develop physical understanding and predictive models of the origin and consequences of localised deformation under irradiation in F/M steels; (2) To develop good practices and possibly advance towards the definition of protocols for the use of ion irradiation as a tool to evaluate radiation effects on materials. Nineteen modelling codes across different scales are being used and developed and an experimental validation programme based on the examination of materials irradiated with neutrons and ions is being carried out. The project enters now its 4th year and is close to delivering high-quality results. This paper overviews the work performed so far within the project, highlighting its impact for fission and fusion materials science. Peer Reviewed L. Malerba a,*, M.J. Caturla b, E. Gaganidze c, C. Kaden d, M.J. Konstantinovi´c e, P. Olsson f, C. Robertson g, D. Rodney h, A.M. Ruiz-Moreno i, M. Serrano a, J. Aktaa c, N. Anento j, S. Austin i, A. Bakaev e, J.P. Balbuena b, F. Bergner d, F. Boioli k, M. Boleininger l, G. Bonny e, N. Castin e, J.B. J. Chapman l, P. Chekhonin d, M. Clozel m, B. Devincre k, L. Dupuy g, G. Diego a, S.L. Dudarev l, C.-C. Fu g, R. Gatti k, L. G´el´ebart g, B. G´omez-Ferrer n, D. Gonçalves g, C. Guerrero a, P.M. Gueye n, P. H¨ahner i, S.P. Hannula o, Q. Hayat p, M. Hern´andez-Mayoral a, J. Jagielski m, N. Jennett p, F. Jim´enez a, G. Kapoor d, A. Kraych h, T. Khvan e,q, L. Kurpaska m, A. Kuronen r, N. Kvashin j, O. Libera s, P.-W. Ma l, T. Manninen o, M.-C. Marinica g, S. Merino a, E. Meslin g, F. Mompiou t, F. Mota a, H. Namburi s, C.J. Ortiz a, C. Pareige n, M. Prester u, R.R. Rajakrishnan t, M. Sauzay g, A. Serra j, I. Simonovski i, F. Soisson g, P. Sp¨atig v, D. Tanguy h, D. Terentyev e, M. Trebala o, M. Trochet g, A. Ulbricht d, M.Vallet g, K. Vogel d, T. Yalcinkaya w, J. Zhao r a Centro de Investigaciones Energ´eticas, Medioambientales y Tecnol´ogicas (CIEMAT), Madrid, Spain b Universidad de Alicante, San Vicente del Raspeig, Spain c Karlsruher Institut für Technologie (KIT), Karlsruhe, Germany d Helmholtz-Zentrum Dresden-Rossendorf Ev (HZDR), Rossendorf, Germany e Studiecentrum voor Kernenergie / Centre d’´Etude de l’´Energie Nucl´eaire (SCK CEN), Mol, Belgium f KTH Royal Institute of Technology, Stockholm, Sweden g Universit´e Paris-Saclay, Commissariat `a l’´Energie Atomique et aux ´Energies Alternatives (CEA), Gif-sur-Yvette, France h Institut Lumi`ere Mati`ere (ILM), Centre National de la Recherche Scientifique, Lyon, France i Joint Research Centre (JRC)- European Commission, Petten, the Netherlands j Universitat Polit`ecnica de Catalunya, Barcelona, Spain k Laboratoire d’Etude des Microstructures (LEM), Centre National de la Recherche Scientifique, Chˆatillon, France l United Kingdom Atomic Energy Authority (UKAEA), Culham, UK m Narodowe Centrum Badan Jadrowych (NCBJ), Swierk, Poland n Normandie Univ, UNIROUEN, INSA Rouen, CNRS, Groupe de Physique des Mat´eriaux, 76000 Rouen, France o Aalto University, Espoo, Finland p Coventry University, UK q Universit´e de Li`ege, Belgium r Helsingin Yliopisto, Helsinki, Finland s Centrum Vyzkumu ˇReˇz S.R.O. (CVR), ˇReˇz, Czech Republic t Centre pour l’´Elaboration Elaboration de Mat´eriaux et pour l’´Etude des Structures (CEMES), Centre National de la Recherche Scientifique, Toulouse, France u Institut za Fiziku, Zagreb, Croatia v Paul Scherrer Institut (PSI), Villingen, Switzerland w Middle East Technical University (METU), Ankara, Turkey

The atomistic structure of calcium-silicate-hydrate (C-S-H) has been studied by applying force field simulations and quantum mechanical methods. Large number of crystallographically equivalent structures were used to account for the disordered nature of C-S-H. An extensive set of periodic structures of C-S-H (with a variation of chemical composition) has been generated and optimised. Special attention has been paid to two crystallographic models for tobermorite, which are used for the construction of C-S-H models. Re-optimisation of atomic structures, initially optimised with ReaxFF force field, by quantum chemical methods (semiempirical and Density functional theory) have been performed. Although the main structural characteristics of C-S-H remain the same in all calculations, the balance between the amount of OH groups and water in the structure is different. The obtained structures (optimised with a different level of theory) are combined into an on-line database and can be used for a future simulation of C-S-H materials.

Mammalian cells produce up to 80% of the commercially available therapeutic proteins, with Chinese Hamster Ovary (CHO) cells being the primary production host. Manufacturing involves a train of reactors, the last of which is typically run in fed-batch mode, where cells grow and produce the required protein. The feeding strategy is decided a priori, from either past operations or the design of experiments and rarely considers the current state of the process. This work proposes a Model Predictive Control (MPC) formulation based on a hybrid kinetic-stoichiometric reactor model to provide optimal feeding policies in real-time, which is agnostic to the culture, hence transferable across CHO cell culture systems. The benefits of the proposed controller formulation are demonstrated through a comparison between an open-loop simulation and closed-loop optimization, using a digital twin as an emulator of the process.

Nasl. z nasl. zaslona. Opis vira z dne 16. 5. 2023. Bibliografija: str. 16-18. Abstract.

Food processing and consumption involves multiple contacts between biological fluids and solid materials of processing devices of which steel is one of the most common. Due to complexity of these interactions, it is difficult to identify the main control factors in the formation of undesirable deposits on the device surfaces that may affect safety and efficiency of the processes. Mechanistic understanding of biomolecule-metal interactions involving food proteins could improve management of these pertinent industrial processes and consumer safety in food industry and beyond. In this work, we perform a multiscale study of formation of protein corona on iron surfaces and nanoparticles in contact with cow milk proteins. By calculating the binding energies of proteins with the substrate we quantify the adsorption strength and rank proteins by the adsorption affinity. We use a multiscale method involving all-atom and coarse-grained simulations as well as generation of 3D structures of the proteins for this purpose. Finally, using the adsorption energy results, we predict the composition of protein corona on iron nanoparticles and flat surfaces via a competitive adsorption model.

As global power generation is currently relying on fossil fuel-based power plants, more anthropogenic CO2 is being released into the atmosphere. During the transition period to alternative energy sources, carbon capture and storage seems to be a promising solution. Chemical-looping combustion (CLC) is an energy conversion technology designed for combustion of fossil fuel with advantageous carbon capture capabilities. In this work, a 1D computational fluid dynamics (CFD) multiscale model was developed to study the reduction step in a syngas-based CLC system and was validated using literature data (R=0.99). In order to investigate mass transfer effects, flow rate and particle dimension studies were carried out. Sharper mass transfer rates were seen at lower flow rates and smaller granule sizes due to suppression of diffusion limitations. In addition, a 3D CFD particle model was developed to investigate in depth the reduction within an ilmenite particle, with focus on heat transfer effects. Minor differences of 1 K were seen when comparing temperature changes predicted by the two models during the slightly exothermic reduction reaction with syngas.

Prediction of material performance in fusion reactor environments relies on computational modelling, and will continue to do so until the first generation of fusion power plants come on line and allow long-term behaviour to be observed. In the meantime, the modelling is supported by experiments that attempt to replicate some aspects of the eventual operational conditions. In 2019, a group of leading experts met under the umbrella of the IEA to discuss the current position and ongoing challenges in modelling of fusion materials and how advanced experimental characterisation is aiding model improvement. This review draws from the discussions held during that workshop. Topics covering modelling of irradiation-induced defect production and fundamental properties, gas behaviour, clustering and segregation, defect evolution and interactions are discussed, as well as new and novel multiscale simulation approaches, and the latest efforts to link modelling to experiments through advanced observation and characterisation techniques. MRG, SLD, and DRM acknowledge funding by the RCUK Energy Programme [grant number EP/T012250/1]. Part of this work has been carried out within the framework of the EUROFusion Consortium and has received funding from the Euratom research and training programme 2014–2018 and 2019–2020 under grant Agreement No. 633053. The views and opinions expressed herein do not necessarily reflect those of the European Commission. JRT acknowledges funding from the US Department of Energy (DOE) through grant DE-SC0017899. ZB, LY,BDW, and SJZ acknowledge funding through the US DOE Fusion Energy Sciences grant DE-SC0006661ZB, LY and BDW also were partially supported from the US DOE Office of Science, Office of Fusion Energy Sciences and Office of Advanced Scientific Computing Research through the Scientific Discovery through Advanced Computing (SciDAC) project on Plasma-Surface Interactions. JMa acknowledges support from the US-DOEs Office of Fusion Energy Sciences (US-DOE), project DE-SC0019157. Pacific Northwest National Laboratory is operated by Battelle Memorial Institute for the US Department of Energy (DOE) under contract DE-AC05-76RL01830. YO and YZ were supported as part of the Energy Dissipation to Defect Evolution (EDDE), an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Sciences under contract number DE-AC05-00OR22725. TS and TT are supported by JSPS KAKENHI Grant Number 19K05338.

Solid oxide fuel cells are high-efficiency renewable energy devices and considered one of the most promising net-zero carbon energy technologies. Numerical modelling is a powerful tool for the virtual design and optimisation of the next-generation solid oxide fuel cells but needs to tackle issues for incorporating the multi-scale character of the cell and further improving the accuracy and computational efficiency. While most of solid oxide fuel cell models were developed based on closed source platforms which limit the freedom of customisation in numerical discretization schemes and community participation. Here, an open source multi-physics and multiscale platform for advanced SOFC simulations consisting of both cell- and pore-scale performance models was developed using OpenFOAM. The modelling aspects are elucidated in detail, involving the coupling of various physical equations and the implementation of the pore-scale electrode in the performance model. The entire platform was carefully validated against experimental data and the other numerical models which were implemented in commercial software ANSYS Fluent and based on the lattice Boltzmann method. The cell-scale model is subsequently employed to study the effects of different fuels, i.e., pure hydrogen and different ratios of pre-reformed methane gas under various operating temperatures. It is found that the cell-scale model reasonably predicts the effects of these parameters on the cell performance, aligning well with the Fluent model. This study further identified the size of representative element volume with respect to the current density for the anode via the pore-scale model where the realistic microstructures reconstructed by a Xe plasma focused ion beam–scanning electron microscopy are employed as computational domains. It is found that a volume element size of 1243 voxels is sufficient to yield the representative current density of the whole. All these numerical investigations show that OpenFOAM is a potential multi-physics and multi-scale computational platform that is capable of accurately predicting both cell-scale and pore-scale performance and spatial information of solid oxide fuel cells. The developed models are also made public in GitHub to inspire community to further develop around it.

Cancer cell invasion is recognised as one of the hallmarks of cancer and involves several inner-related multiscale processes that ultimately contribute to its spread into the surrounding tissue. In order to gain a deeper understanding of the tumour invasion process, we pay special attention to the interacting dynamics between the cancer cell population and various constituents of the surrounding tumour microenvironment. To that end, we consider the key role that ECM plays within the human body tissue, providing not only structure and support to surrounding cells, but also acting as a platform for cells communication and spatial movement. There are several other vital structures within the ECM, however we are going to focus primarily on fibrous proteins, such as fibronectin. These fibres play a crucial role in tumour progression, enabling the anchorage of tumour cells to the ECM. In this work we consider the two-scale dynamic cross-talk between cancer cells and a two component ECM (consisting of both a fibre and a non-fibre phase). To that end, we incorporate the interlinked two-scale dynamics of cells-ECM interactions within the tumour support that contributes simultaneously both to cell-adhesion and to the dynamic rearrangement and restructuring of the ECM fibres. Furthermore, this is embedded within a multiscale moving boundary approach for the invading cancer cell population, in the presence of cell-adhesion at the tissue scale and cell-scale fibre redistribution activity and leading edge matrix degrading enzyme molecular proteolytic processes. The overall modelling framework will be accompanied by computational results that will explore the impact on cancer invasion patterns of different levels of cell adhesion in conjunction with the continuous ECM fibres rearrangement., Comment: 44 pages, 17 figures

This paper aims at broadening the perspective in life cycle assessment methodology, exploiting multiscale modelling towards the generation of life cycle inventory data during an early-stage product design. Our approach involves the usage of molecular modelling techniques, such as electronic, atomistic or mesoscale models, in combination with continuum models, such as process simulation or finite element methods, to provide data for the generation of life cycle inventories of complex materials and production processes. In particular, each simulation is performed at a specific length and time scale through dedicated software, passing information from the lower to the upper scale. A practical application of the proposed approach will be illustrated for property predictions of nanostructured polymer systems used for manufacturing a marine engine cover. The material formulations capable of fulfilling the marine thermomechanical constraints were defined using in-silico techniques, allowing for the generation of materials life cycle inventories, which have been compared to identify the most environmentally friendly option.

A better understanding of the time-dependent biomechanical behaviour of the biphasic hip articular cartilage (AC) under physiological loadings is important to understand the onset of joint pathology and guide the clinical treatment. Current computational studies for the biphasic hip AC were usually limited to short-term duration or using elaborate loading. The present study aimed to develop a multiscale computational modelling to investigate the long-term biphasic behaviour of the hip AC under physiological loadings over multiple gait cycles. Two-scale computational modelling including a musculoskeletal model and a finite element model of the natural hip was created. These two models were then combined and used to investigate the biphasic behaviour of hip AC over 80 gait cycles. The results showed that the interstitial fluid pressure in the AC supported over 89% of the loading during gait. When the contact area was located at the AC centre, the contact pressure and fluid pressure increased over time from the first cycle to the 80th cycle, while when the contact area approached the edge, these pressures decreased first dramatically and then slowly over time. The peak stresses and strains in the solid matrix of the AC remained at a low level and increased over time from the first cycle to the 80th cycle. This study demonstrated that the long-term temporal variations of the biphasic behaviour of hip AC under physiological loadings are significant. The methodology has potentially important implications in the biomechanical studies of human cartilage and supporting the development of cartilage substitution.

The influence of the inherent fabric anisotropy of sand and its evolution on the footing bearing capacity and failure mechanism has been modelled and analysed with a multiscale approach that couples the finite-element method (FEM) and the discrete-element method (DEM). In this approach, the material constitutive responses are captured by DEM simulations on representative volume elements (RVEs) assigned to the FEM Gauss points. The inherent fabric anisotropy of sand is reflected by the bedding plane angle α of the RVEs composed of elliptic particles, where α measures the average angle between the long axes of the particles and the horizontal direction. The bearing capacity is found to decrease with an increasing α and the reduction can reach as high as 25%. Asymmetric deformation has been observed in the α other than 0° or 90° cases, which is mainly attributable to the different relative principal directions between the stress and the fabric at the two sides of the domain. Further multiscale analyses reveal that the load-bearing structure of the material expressed by the contact-normal-based fabric tensor is more fragile under the footing loading at one side of the domain and in larger α cases, which offers a micromechanical explanation to the observation of the asymmetric deformation and the decreased bearing capacity with increased α.

In recent years, multiscale modelling approach has begun to receive an overwhelming appreciation as an appropriate technique to characterize the complexity of infectious disease systems. In this study, we develop an embedded multiscale model of paratuberculosis in ruminants at host level that integrates the within-host scale and the between-host. A key feature of embedded multiscale models developed at host level of organization of an infectious disease system is that the within-host scale and the between-host scale influence each other in a reciprocal (i.e., both) way through superinfection, that is, through repeated infection before the host recovers from the initial infectious episode. This key feature is demonstrated in this study through a multiscale model of paratuberculosis in ruminants. The results of this study, through numerical analysis of the multiscale model, show that superinfection influences the dynamics of paratuberculosis only at the start of the infection, while the MAP bacteria replication continuously influences paratuberculosis dynamics throughout the infection until the host recovers from the initial infectious episode. This is largely because the replication of MAP bacteria at the within-host scale sustains the dynamics of paratuberculosis at this scale domain. We further use the embedded multiscale model developed in this study to evaluate the comparative effectiveness of paratuberculosis health interventions that influence the disease dynamics at different scales from efficacy data.

Haematococcus pluvialis can produce significant amounts of industrially important compounds belonging to lipids and starch classes, including various specific pigments such as β-carotene, lutein and astaxanthin, as well as lipids, carbohydrates and proteins. Their production can vary depending on environmental stress conditions like nutrient starvation. However, stress conditions lead also to undesired phenomena such as cell lysis, which is likely to be related to products loss. The microorganism develops towards smaller single cell volumes during the growth process, and eventually, more likely towards lysis when fission (i.e. cell division) slows down. The lysis process takes place simultaneously with nutrient depletion, so both growth and lysis are linked to the change of environmental conditions. In this work, we develop a novel multiscale segregated-structured model based on Population Balance Equations (PBEs) to describe the photoautotrophic growth of H.pluvialis, in particular cell growth, and lysis, making possible the description of the relationship between cell volume/transition, cell loss, and metabolic product availability. Cell volume is the internal coordinate of the population balance model, and its link with intrinsic concentrations is also presented. The model parameters are fitted against experimental data, extensive sensitivity analysis is performed and the model predictive capabilities are tested in terms of cell density distributions, as well as 0th and 1st order moments.

Light-weight fabric membranes have gained increasing popularity over the past years due to their tailorable structural and material performances. These tailorable properties include stretch forming and deep drawing formability that exhibits excellent stretchability and drapeability properties of textiles and textile composites. Since the inception of computerised numerical control for three-dimensional textile-manufacturing machines, technical textiles paved their way to numerous applications, certainly not limited to; aerospace, biomedical, civil engineering, defence, marine and medical industries. Digital interlooping and digital interlacing technology in additive manufacturing greatly advanced the manufacturing processes of textiles. In this work, we consider two branches of technical fabrics, namely plain-woven and weft-knitted. Multiscale modelling is the tool of choice for homogenising periodic structures and has been used extensively to model and analyse the mechanical behaviour of woven and knitted fabrics. But there is a plethora of literature discussing the demerits of such conventional multiscale modelling. These demerits include higher computational costs, rigid numerical models, ineffcient algorithmic computations and inability to incorporate geometric nonlinearities. We propose a data-driven nonlinear multiscale modelling technique to analyse the complex mechanical behaviour of plain-woven and weft-knitted fabrics with a neat extension to fabric material designing. We show how the integration of statistical learning techniques mitigates the weaknesses of conventional multiscale modelling. Moreover, we discuss the avenues that will open in many potential fields with regard to material modelling, structural engineering and textile industries. In the proposed data-driven nonlinear computational homogenisation technique, we effi ciently integrate the microscale and macroscale using Gaussian Process Regression (GPR) statistical learning technique. In the microscale, representative volume elements (RVEs) are modelled using nite deformable isogeometric spatial rods and deformation is homogenised using periodic boundary conditions. This nite deformable rod is profi cient in handling large deformations, rod-to-rod contacts, arbitrary cross-section de finitions and follower loads. Respecting the principle of separation of scales, we construct response databases by applying different homogenised strain states to the RVEs and recording the respective incremental volume-averaged energy values. We use GPR to learn a model using a 5-fold cross-validation technique by optimising the log marginal likelihood. In the macroscale, textiles are modelled as nonlinear orthotropic membranes for which the stresses and material constitutive relations are predicted by the trained GPR model. This coupling between GPR and membrane models is achieved through a systematic and seamless nite element integration using C++ and Python environments. A neat extension to material designing is also discussed with potentials to extend the work into other related fi elds., Cambridge trust and Trinity Hall scholarship

The present study aims to predict the performance of hydrostatic journal bearings during steady-state operation through multiscale modelling and analysis. An innovative multiscale modelling approach is presented for numerical modelling and analysis of the hydrostatic journal bearing as configured and designed, particularly for the precision press machine. With the challenging environment for journal bearings under hydrostatic pressures and surface features, elastic modulus of oil and viscosity-temperature characteristics are often neglected while against heavy-loading and precision engineering requirement. This research takes into account of oil elastic modulus and viscosity-temperature characteristics to investigate the intrinsic relationship between the pressures within the hydrostatic bearing and its surface texture, feature size and locations being considered. This paper concludes with a further discussion on the potential and application of the approach in precision engineering.

Propane (C3H8) and butane (C4H10) are short straight-chain alkane molecules that are difficult to convert catalytically. Analogous to propane, butane can be dehydrogenated to butenes (also known as...

Acetohydroxamic acid (AHA) is a small organic acid with a wide variety of industrial, biological, and pharmacological applications. A deep fundamental molecular level understanding of the mechanisms responsible for the radical-induced reactions of AHA in these environments is necessary to predict and control their behaviour and elucidate their interplay with other attendant chemical species, for example, the oxidative degradation products of AHA. To this end, we present a comprehensive, multiscale computer model for interrogating the radical-induced degradation of AHA in acidic aqueous solutions. Model predictions were critically evaluated by a systematic experimental radiation chemistry investigation, leveraging time-resolved electron pulse irradiation techniques for the measurement of new radical reaction rate coefficients, and steady-state gamma irradiations for the identification and quantification of AHA degradation products: acetic acid, hydroxylamine, nitrous oxide, and molecular hydrogen, with formic acid and methane as minor products. Excellent agreement was achieved between calculation and experiment, indicating that this fundamental model can accurately predict the degradation pathways of AHA under irradiation in acidic aqueous solutions.

Molecular modelling applications in metalloenzyme design are still scarce due to a series of challenges. In top of that, the simulations of metal-mediated binding and the identification of catalytic competent geometries requiring both large conformational exploration and simulation of fine electronic properties. Here, we demonstrate how the incorporation of new tools in multiscale strategies, namely substrate diffusion exploration, allow taking a step further. As a showcase the enantioselective profiles of the most outstanding variants of an artificial Rh2-based cyclopropanase (GSH, HFF and RFY) developed by Lewis and co-workers have been rationalized. DFT calculations on the free-cofactor-mediated process identify the carbene insertion and the cyclopropanoid formation as crucial events, being the latter the enantiodetermining step, which displays up to 8 competitive orientations easily altered by the protein environment. The key intermediates of the reaction were docked into the protein scaffold showing that some mutated residues have direct interaction with the cofactor and/or the co-substrate. These interactions take form of a direct coordination of Rh in GSH and HFF and a strong hydrophobic patch with the carbene moiety in RFY. Posterior molecular dynamics sustain that the cofactor induces global re-arrangements of the protein. Finally, massive exploration of substrate diffusion, based on the GPathFinder approach, defines this event as the origin of the enantioselectivity in GSH and RFY. For HFF, fine molecular dockings suggest that it is likely related to local interactions upon diffusion. This work shows how modelling of long-range mutations on the catalytic profiles of metalloenzymes may be unavoidable and software simulating substrate diffusion should be applied.

A multiscale mathematical model describing the genesis and ecology of algal-bacterial photogranules and the metals biosorption on their solid matrix within a sequencing batch reactor (SBR) is presented. The granular biofilm is modelled as a spherical free boundary domain with radial symmetry and a vanishing initial value. The free boundary evolution is governed by an ODE accounting for microbial growth, attachment and detachment phenomena. The model is based on systems of PDEs derived from mass conservation principles. Specifically, two systems of nonlinear hyperbolic PDEs model the growth of attached species and the dynamics of free adsorption sites; and two systems of quasi-linear parabolic PDEs govern the diffusive transport and conversion of nutrients and metals. The model is completed with systems of impulsive ordinary differential equations (IDEs) describing the evolution of dissolved substrates, metals, and planktonic and detached biomasses within the granular-based SBR. All main phenomena involved in the process are considered in the mathematical model. Moreover, the dual effect of metal presence on the formation process of photogranules is accounted: metal stimulates the production of EPS by sessile species and negatively affects the metabolic activities of microbial species. To describe the effects related to metal presence, a stimulation term for EPS production and an inhibition term for metal are included in all microbial kinetics. The model is used to examine the role of the microbial species and EPS in the adsorption process, and the effect of metal concentration and adsorption proprieties of biofilm components on the metal removal. Numerical results show that the model accurately describes the photogranules evolution and ecology and confirm the applicability of algal-bacterial photogranules systems for metal-rich wastewater treatment., Comment: 43 pages, 18 figures, preprint version

The complexity of the thermo-hydro-mechanical processes in expansive clays is such that coupled models are needed to explain their behaviour under variable hydraulic, mechanic, and thermal conditions. The effects of temperature changes on the behaviour of bentonite are being studied, in particular, with respect to their role in the chemical and hydraulic insulation of radioactive waste. Results demonstrate that thermal effects can be significant. These effects can also occur owing to different processes, such as seasonal thermal variations in the shallow (possibly clay-rich) soil layers along mountain slopes. Very limited research in this field has shown that some landslides might be activated by thermal excursions. In the proposed project, experiments and modelling at various scales have been performed to better account for thermo-hydro-mechanical coupling and for the control of temperature on the shear strength, and hence on the stability of slopes. An extensive variety of geotechnical parameters are subjected to changes under thermal action. For instance, we need to mention the shear strength. To develop a more complete framework to account for temperature in slope stability analysis, a multi-scale approach has been defined to improve our understanding on how soil behaves under thermal changes....

How morphogenesis depends on cell properties is an active direction of research. Here, we focus on mechanical models of growing plant tissues, where microscopic (sub)cellular structure is taken into account. In order to establish links between microscopic and macroscopic tissue properties, we perform a multiscale analysis of a model of growing plant tissue with subcellular resolution. We use homogenization to rigorously derive the corresponding macroscopic tissue scale model. Tissue scale mechanical properties are computed from all microscopic structural and material properties, taking into account deformation by the growth field. We then consider case studies and numerically compare the detailed microscopic model and the tissue-scale model, both implemented using finite element method. We find that the macroscopic model can be used to efficiently make predictions about several configurations of interest. Our work will help making links between microscopic measurements and macroscopic observations in growing tissues.

Propane (C3H8) and butane (C4H10) are short straight-chain alkane molecules that are difficult to convert catalytically. Analogous to propane, butane can be dehydrogenated to butenes (also known as butylenes) or butadiene, which are used industrially as raw materials when synthesizing various chemicals (plastics, rubbers, etc.). In this study, we present results of detailed first-principles-based multiscale modelling of butane dehydrogenation, consisting of three size- and time-scales. The reaction is modelled over Cr2O3(0001) chromium oxide, which is commonly used in the industrial setting. A complete 108-step reaction pathway of butane (C4H10) dehydrogenation was studied, yielding 1-butene (CH2CHCH2CH3) and 2-butene (CH3CHCHCH3), 1-butyne (CHCCH2CH3) and 2-butyne (CH3CCCH3), butadiene (CH2CHCHCH2), butenyne (CH2CHCCH), and ultimately butadiyne (CHCCCH). We include cracking and coking reactions (yielding C1, C2, and C3 hydrocarbons) in the model to provide a thorough description of catalyst deactivation as a function of the temperature and time. Density functional theory calculations with the Hubbard U model were used to study the reaction on the atomistic scale, resulting in the complete energetics and first-principles kinetic parameters for the dehydrogenation reaction. They were cast in a kinetic model using mean-field microkinetics and kinetic Monte Carlo simulations. The former was used to obtain gas equilibrium conditions in the steady-state regime, which were fed in the latter to provide accurate surface kinetics. A full reactor simulation was used to account for the macroscopic properties of the catalytic particles: their loading, specific surface area, and density and reactor parameters: size, design, and feed gas flow. With this approach, we obtained first-principles estimates of the catalytic conversion, selectivity to products, and time dependence of the catalyst activity, which can be paralleled to experimental data. We show that 2-butene is the most abundant product of dehydrogenation, with selectivity above 90% and turn-over frequency above 10–3 s–1 at T = 900 K. Butane conversion is below 5% at such low temperature, but rises above 40% at T > 1100 K. Activity starts to drop after ∼6 h because of surface poisoning with carbon. We conclude that the dehydrogenation of butane is a viable alternative to conventional olefin production processes.

The cane toad (Bufo marinus) is an invasive species in Australia that has a negative impact on native species. Control methods such as trapping, fencing, and water exclusion have been devised to contain the spread of cane toads and reduce their ecological impact. However, implementing these interventions is expensive, and estimating the likely impact of a proposed intervention on spread at a large spatial scale, comprising one or more control methods, is challenging due to the lack of large-scale data and the computational cost of modelling a large number of toads.To address this challenge, we developed a multiscale model which uses individual-level data on cane toad behaviour to estimate the impact of trapping, fencing, and water exclusion when applied at scale in the Pilbara region in north-western Australia. Compared to previous work, our model allows us to explore more complex combinations and tradeoffs of control methods by utilising data sources at different scales.Our results suggest that exclusion of toads from water points is the most effective method for containing the spread of cane toads, and that trapping and fencing alone are unlikely to be sufficient. However, trapping and fencing are still useful supplementary measures in scenarios where exclusion cannot be broadly applied to a large number of water points.Synthesis and applications. Our analyses highlight the importance of limiting access to sheltering and breeding sites in invasive species control. Furthermore, this study illustrates the value of multiscale computational models for exploring scenarios where parameters and calibration data are available at the scale of individuals and small groups, but management questions are framed at a much larger scale.

Acetohydroxamic acid (AHA) has been proposed for inclusion in advanced, single-cycle, used nuclear fuel reprocessing solvent systems for the reduction and complexation of plutonium and neptunium ions. For this application, a detailed description of the fundamental degradation of AHA in dilute aqueous nitric acid is required. To this end, we present a comprehensive, multiscale computer model for the coupled radiolytic and hydrolytic degradation of AHA in aqueous sodium nitrate and nitric acid solutions. Rate coefficients for the reactions of AHA and hydroxylamine (HA) with the oxidizing nitrate radical were measured for the first time using electron pulse radiolysis and used as inputs for the kinetic model. The computer model results are validated by comparison to experimental data from steady-state gamma ray irradiations, for which the agreement is excellent. The presented model accurately predicts the yields of the major degradation products of AHA: acetic acid, HA, nitrous oxide, and molecular hydrogen.