In this study we investigate the mechanochemical regeneration of sodium borohydride (NaBH₄) from a system comprising hydrated sodium metaborate ( [Formula presented] ) and magnesium hydride (MgH₂). We explore the individual and joint impact of key operational parameters (rotational speed, milling time, ball-to-powder ratio (BPR), and molar ratio) on the regeneration yield. Furthermore, a method for quantifying chemical conversion is introduced relying only on water and thus, offering environmental benefits. This approach additionally facilitates the production and storage of a “ready-to-use” NaBH₄ solution with minimal losses at room temperature. Notably, a yield of 90% is achieved, with a 20% reduction in rotational speed compared to prior literature. This research contributes to sustainable hydrogen storage and presents practical advancements in mechanochemical processes.

Mechanical metamaterials are architected structures with unique functionalities, such as negative Poisson's ratio and negative stiffness, which are widely employed for absorbing energy of quasi-static and impact loads, giving improved mechanical response. Acoustic/elastic metamaterials, their dynamic counterparts, rely on frequency-dependent properties of their microstructure elements, including mass density and bulk modulus, to control the propagation of waves. Although such metamaterials introduced significant contribution for solving independently static and dynamic problems, they were facing certain resistance to their use in real-world engineering problems, mainly because of a lack of integrated systems possessing both mechanical and vibration attenuation performance. Advances in manufacturing processes and material and computational science now enable the creation of hybrid mechanical metamaterials, offering multifunctionality in terms of simultaneous static and dynamic properties, giving them the ability of controlling waves while withstanding the applied loading conditions. Exploring towards this direction, this review paper introduces the hybrid mechanical metamaterials in terms of their design process and multifunctional properties. We emphasize the still remaining challenges and how they can be potentially implemented as engineering solutions.

Monopiles are the dominant foundation type for offshore wind turbines, accounting for approximately 80% of the installed capacity. Installing offshore monopile foundations on seabeds susceptible to scour erosion requires monopiles to penetrate several pre-installed scour protection rock layers before securing them into the seabed. The accurate prediction of the pile penetration resistance is crucial to ensure successful monopile installations. To complement, and potentially reduce the dependence on the costly and labour-intensive experimental small-scale penetration tests, a numerical model has been developed using the Discrete Element Method (DEM) that captures the discrete nature of interactions between rocks and piles and predicts the resistance during the penetration process. The developed DEM model includes armour and filter rocks represented by multispheres and sand particles represented by spheres. A multistage calibration, verification and validation DEM modelling framework is proposed and examined with small-scale penetration tests conducted using plates and piles in a double-layer scour protection configuration. The sand material model is calibrated and verified using penetrometer tests and the rock material models are calibrated and verified using a plate penetration test. The DEM model with three verified materials predicts the penetration resistance well in small-scale pile penetration tests and proves the validity of the proposed framework. The DEM model presented in this paper facilitates the modelling in areas where traditional continuum-based numerical methods give less accurate predictions and provide insights that are difficult or nearly impossible to obtain through experimental methods.

Conventional scour protection installation around monopiles for offshore wind farms involves placement of smaller filter layer rocks and larger armor layer rocks in two separate operations, requiring multiple visits of the rock dumping vessel to the site which increases costs, time spent offshore by the vessels, and consequently carbon emissions. The joint industry project Optimizing Pile Installation through Scour Protection (OPIS), established in 2023, helps reaching the target of energy transmission by streamlining the pile and scour protection installations, reducing costs and carbon emissions associated with offshore wind developments. The project investigates the technical feasibility of pile installation through coarse rock to enhance the feasibility of the operation. A series of small and medium scale laboratory experiments has been conducted, considering different scour protection designs such as single- and double-layer systems, and different rock densities (high and normal density rocks). This paper delves into the physical modelling aspects of the OPIS research project. Furthermore, this contribution elaborates on the design and intricacies of the scaled laboratory tests, providing in-depth insights into the design and implementation of laboratory setups, along with a detailed account of experimental procedures. Moreover, preliminary results and challenges encountered during the experiments are discussed, and innovative solutions devised to overcome specific challenges are highlighted.

The use of optimization procedures for designing acoustic/elastic metamaterials (A/E MMs) has gained significant interest since they enable the efficient attainment of unique functionalities often contradicting. When it comes to vibration attenuation caused by mechanical stress waves, such as impact loads, the dynamic properties of A/E MMs are optimized so that their wave-control ability is maximized. However, the mechanical performance of A/E MMs during the propagation of such waves is normally not evaluated into the design optimization stages. This may compromise not only the load-bearing capacity of MMs, but also their ability in attenuating vibrations. To prevent such effects, we propose a design strategy that incorporates the stress analysis in the early design phase of A/E MMs subjected to an impact load. The effective mass density approach is applied, from which the vibration attenuation is identified at frequency ranges where the resonator moves out-of-phase in relation to the applied excitation. Regarding to the A/E MM mechanical behavior, maximum von Mises stress is calculated through the transient analysis of a unit cell array subjected to a dynamic load. A Pareto front shows a trade-off behavior between the A/E MM functionalities. With that, we emphasize the importance of incorporating the mechanical performance into the design stage of A/E MMs for vibration attenuation of structures undergoing high impact loads, such as installation of foundations by impact hammering. This brings A/E MMs closer to real applications involving energy filtering at specific frequencies from transient loads, designed in an optimized and efficient way.

High-energy ball milling is a versatile method utilized in mechanochemical reactions and material transformations. Understanding and characterizing the relevant mechanical variables is crucial for the optimization and up-scaling of these processes. To achieve this, the present study delves into differentiating the contributions of normal and tangential interactions during high-energy collisions. Using Discrete Element Method (DEM) simulations, we characterize how operational parameters influence these energy dissipation modes, emphasizing the significance of tangential interactions. Our analysis also reveals how different operational parameters such as ball size, fill ratio, and rotational speed affect the mechanical action inside the milling jar giving rise to multiple operating zones where different modes of energy dissipation can thrive. Finally, we present master curves that generalize findings across a wide range of configurations, offering a tool for characterizing and predicting mechanochemical processes beyond the presented cases. These results provide a robust framework for improving mechanochemical reaction efficiency, and equipment design.

To reduce global emissions, hydrogen is increasingly considered as an energy carrier for renewable energy storage. However, traditional storage methods for hydrogen such as compression or liquefaction require high pressures, extremely low temperatures, and still result in a low volumetric energy density. As a solution, sodium borohydride (NaBH₄) is proposed as an alternative method to store hydrogen. NaBH₄ is a granular material that can be stored using ambient temperature and pressure, and has a relatively high volumetric and gravimetric energy density compared to traditional hydrogen storage. This paper explores the application of NaBH₄ as a fuel in the maritime industry, and elaborates on how the use of NaBH₄ leads to a circular bunkering (refuelling) process. By using hydrolysis to extract hydrogen from NaBH₄ during vessel operation, a so called spent fuel remains and needs to be stored on the vessel until next port call. Additionally, examples of various bunkering equipment that can be used to design the circular bunkering process of NaBH₄ are presented. Moreover, it explains how design of bunkering equipment depends on the mechanical characteristics of the fuel and spent fuel. The main finding of this work is that NaBH₄ is a promising solution for a sustainable future. Before NaBH₄ can be used as a fuel, vessels and ports need to be adapted to facilitate circular bunkering with such a novel solid-state energy carrier.

Calibration of discrete element method (DEM) models is crucial for the realistic simulation of granular materials. However, it remains a challenging task, especially for multi-component mixtures due to their higher complexity and larger number of parameters involved. This study presents a systematic and computationally efficient calibration framework designed to address these challenges, focusing on pellet-sinter mixtures, as a representative case of two-component mixtures commonly used in blast furnace steelmaking. The framework integrates sensitivity analysis, machine learning-based surrogate modelling with adaptive sampling, and genetic algorithm-driven optimisation techniques to minimise the number of required DEM simulations. Using this approach, we achieved a high-accuracy surrogate model (R² = 0.95) for seven DEM parameters with only 110 data points, highlighting the efficiency and robustness of the framework. These parameters were successfully calibrated with a relative error of less than 2 %. Moreover, the calibrated parameters for the base case (i.e., 50–50 pellet-sinter mass ratio) remained valid across different mass ratios and layering orders, eliminating the need for recalibration. Overall, the proposed framework offers a reliable, cost-effective, and adaptable solution for DEM calibration of two-component mixtures. Its flexibility and efficiency make it a promising tool for extending to more complex systems, facilitating the development of DEM models for a wide range of industrial applications involving granular mixtures.

The International Maritime Organization aims to achieve full decarbonization by 2050 in response to climate change. This ambitious goal demands well-defined strategies guided by techno-economic assessments. The complexity of global shipping systems makes predicting long-term maritime trade patterns challenging, necessitating scenario-building rather than precise forecasts. Investigating shipping demand scenarios is crucial due to the uncertainty brought by the energy transition and its role as the primary driver of shipping emissions. This paper improves the representation of maritime shipping in Integrated Assessment Models (IAMs) by examining the impacts of climate targets on future shipping demand. A novel econometric model, grounded in advanced gravity theory and integrated with machine-learning algorithms, is proposed to estimate the elasticities of variables in bilateral seaborne trade. By coupling this model with the WITCH IAM, we explore various scenarios, providing deeper insights into trade patterns and their implications. The results show that stricter climate policies and higher carbon taxes reduce GDP due to higher abatement costs, higher fuel prices, and therefore reduced seaborne trade, especially for oil products and containerized cargo. Early adoption of carbon taxes in Europe may shift oil production and consumption patterns, temporarily boosting seaborne trade. Sub-Saharan Africa could experience significant demand growth due to economic and population increases.

The offshore wind industry has grown significantly over the past decade with the increasing demand for clean energy, and over 80% of offshore wind turbines have monopile foundations. The installation of these foundations requires monopiles to penetrate several pre-installed scour protection rock layers before reaching the required penetration depth. With the increasing sizes of both monopiles and the extent of scour protection layers, various challenges arise during the monopile’s driving and self-weight penetration process, such as misalignment and early pile refusal. Lab scale testing and empirical modelling are time-consuming and normally not general enough to be applied to different scenarios. Within the Optimising Pile Installation through Scour Protection (OPIS) project, the high fidelity discrete element method (DEM) has been deployed to consider both the discrete nature and shapes of the scour protection rocks. A double-layer (armour & filter) scour protection with a sand layer underneath was simulated using the DEM. Detailed rock dynamics and the armour rocks dragging down during penetration are revealed and particle-geometry rolling friction is identified as the dominant parameter. The effects of the sand layer thickness and pile wall thickness on penetration resistance are also addressed.

Although DEM (Discrete Element Method) was introduced >40 years ago, there are still various challenges in applying it with a certain level of confidence. To develop a realistic material model, calibration, validation, particle shape representation and scaling are well-known challenges that have been studied by various researchers over the past two decades. In addition, once a realistic DEM material model is developed and published in line with open science principles, it is unknown to what extent the calibrated values can be used across software packages. Although a benchmark between 8 open source codes was recently published, it did not cover cohesive materials, rolling friction models, realistic calibrated material behaviour, and commercial software that is widely used in industries.

The aim of this study is to investigate the portability of input parameters between different codes and to identify if software independent calibration is possible. We propose a framework to assist DEM users in industry and academia to reach a software independent DEM simulation when required.

To compare the results between software packages, the framework considers the following aspects: 1) identical implementation of contact models, 2) single contact modelling, 3) bulk level simulation, and 4) identical post-processing. This framework is demonstrated for two contact models (Hertz-Mindlin and Edinburgh Elasto Plastic Adhesive) with rolling friction model and two commercial software packages (EDEM and PFC) but is applicable for any other (open-source) software. Moreover, two realistic calibrated material models are included to cover a range of material characteristics, from free-flowing incompressible materials to cohesive compressible materials. This study shows that individual particle level simulations give identical results, bulk level simulations show differences for both contact models. In general, users should use parameter values from other software packages with caution, especially where critical or sensitive applications are modelled. This paper also highlights the use of novel computation techniques, such as the GPU engine, to achieve practical computation times when modelling industrial applications.

Greenhouse gas emissions drive global warming, posing significant risks to ecosystems and human society. The maritime sector contributes approximately 3% of global emissions, leading the International Maritime Organization (IMO) to target a 40% reduction in emissions by 2030 and 70% by 2050, compared to 2008 levels. This reduction can be achieved by improving vessel efficiency or adopting alternative fuels such as hydrogen. Traditionally, hydrogen is stored as a gas, liquid, or cryocompressed, but these methods have drawbacks, including high pressure and energy demands. A promising alternative is sodium borohydride (NaBH4), a solid hydrogen carrier that offers high volumetric energy density and safe storage at ambient conditions. When used in maritime vessels, NaBH4 reacts with water to produce hydrogen and a byproduct, which is stored and later regenerated onshore. Both NaBH4 and its spent fuel are granular materials, requiring specific handling equipment. This paper aims to identify the mechanical characteristics of NaBH4 and its spent fuel to design appropriate storage and handling systems using the Discrete Element Method (DEM).

Sodium borohydride (NaBH4) is considered as an alternative fuel for the maritime industry [1]. In contrast to conventional fuels, NaBH4 is a granular material. To use a simulation-supported design for assessing the feasibility of equipment designs for storing and handling this material, its mechanical characteristics are required. These are then used to calibrate and verify simulations using the Discrete Element Method (DEM). However, as this is a novel application for this material, virtually no bulk characteristics are known yet. Therefore extensive testing has been done to extract required mechanical characteristics, such as cohesion, adhesion, internal friction, wall friction, and the Angle of Repose (AoR). These experiments showed that NaBH4 is initially free-flowing, but an increase in the moisture content because of an increase in relative humidity leads to an increase in cohesion, effectively reducing the flowability of the bulk material. Furthermore, our experimental results showed plastic deformation of individual NaBH4 particles.
This work focuses on capturing both the free-flowing and cohesive behaviour of NaBH4 in DEM. To this end, the two-step calibration approach introduced by Grima [2] is adopted and adjusted. First, the free-flowing behaviour is calibrated using a non-cohesive contact model, Hertz-Mindlin (HM). Second, the cohesive material is calibrated using the appropriate cohesive parameters of the Edinburgh Elasto-Plastic Adhesion contact model (EEPA), while the (calibrated) non-cohesive parameters are kept constant. The novelty of this work is the use of EEPA for the second calibration step, which allows the modelling of both the cohesive behaviour and the plastic deformation of the individual particles in the bulk material.


[1] M.C. van Benten, J.T. Padding, and D.L. Schott. “Towards Hydrogen-Fuelled Marine Vessels using Solid Hydrogen Carriers”. In: The 14th International Conference on Bulk Materials Storage, Handling and Transportation. Wollongong, Australia, July 2023.
[2] A. Grima. “Quantifying and modelling mechanisms of flow in cohesionless and cohesive granular materials”. In: University of Wollongong Thesis Collection 1954- 2016 (Jan. 2011).

Bed permeability is a crucial factor in blast furnace performance which depends on the material distribution achieved through charging. Since a homogeneous bed of pellet and sinter is recommended, it is crucial to understand whether segregation of the pellet-sinter mixture occurs during charging. The Discrete Element Method is useful in this regard; however, simulations of pellet-sinter mixture charging currently lack credibility since pellet-sinter interaction parameters have not yet been calibrated and validated. Determining pellet-sinter interaction parameters will require significant efforts, so it is useful to know whether mixture segregation and the resulting bed permeability are sensitive to these parameters. In this work, we investigate to what extent the restitution coefficient, sliding friction coefficient and rolling friction coefficient between pellet and sinter affect segregation during bed formation and the resulting permeability in terms of porosity using a simplified charging setup. The investigation is done for different mixture compositions and flow velocities, and analysis settings including sample size and sampling directions. We conclude that all parameters affect segregation and porosity, regardless of the composition and velocity. Hence, all mixture parameters including the interaction parameters between the components must be carefully calibrated when developing a model for predicting permeability.