In blast furnace ironmaking, a mixture of iron ore pellets and sinter is charged in layers at the furnace top, with particle velocities reaching up to ∼10 m/s at the stock surface. The inherent differences in particle size, shape, and density between pellets and sinter pose challenges for maintaining a uniform mixture during this high-velocity charging, leading to segregation and uneven material distribution. This non-uniformity can negatively affect furnace efficiency and stability. Understanding segregation during charging is therefore crucial for optimizing the ironmaking process. The Discrete Element Method (DEM) can offer valuable insights, provided that the model parameters are calibrated and verified. This study presents a calibrated DEM model for a pellet–sinter mixture with a 50–50 mass ratio of both components. A novel high-velocity laboratory setup was used to simultaneously measure five different key performance indicators (KPIs) related to flow and packing behavior at various discharge heights, corresponding to different flow velocities. Calibration was performed at the highest flow velocity, representative of actual blast furnace conditions. The process involved creating response surface models for each KPI and using a multi-objective optimization approach with a desirability function to determine the model parameters. A step-wise calibration strategy was employed, first optimizing pellet and sinter interaction parameters individually, followed by calibration of the pellet–sinter interaction parameters. This approach proved effective, as the calibrated model accurately reproduced experimental data. Results also suggest that the calibration outcome is flow-invariant in this setup, with the model successfully predicting flow and packing behavior at lower discharge heights.

While the Discrete Element Method (DEM) provides high-fidelity insights into granular processes like high-energy ball milling, its computational cost can become prohibitive when exploring extensive parameter spaces required for scale-up. This limitation hinders the rapid design and optimization cycles crucial for emerging applications, like mechanochemistry. Surrogate modeling offers a promising path to overcome these computational barriers, yet existing approaches often struggle to accurately represent the complex, moving boundary conditions typical of milling equipment. In this work, we leverage a Signed Distance Function Graph Neural Network (SGN) surrogate tailored to the high-energy, moving-boundary regime of the Emax mill. Trained on DEM data, the SGN jointly predicts particle kinematics for recursive roll-out and a mechanochemistry-relevant global quantity, the global dissipated energy. The model exhibits strong generalization to unseen motion trajectories and moderate jar-shape edits without retraining, while operating with a timestep over 100x larger than required by DEM. In a CPU-only comparison, it achieves a minimum of 6.6× wall-clock speedup. This approach provides a powerful and promising technique for the simulation, analysis, and design optimization of high-energy ball milling equipment.

Segregation of the ferrous burden during blast furnace (BF) charging can cause uneven layer formation at the furnace throat, reducing bed permeability and disrupting gas–solid interaction. This study applies a discrete element method (DEM) model to the industrial-scale BF charging system (from the skip car to top hopper discharge) to examine segregation under real operating conditions. The model includes the full ferrous mixture (pellets, sinter, lump ore, and nut coke) and the real-scale geometries. A reference case representing current practice is analysed in detail and compared with systematically varied case studies. The results show that segregation generally decreases from the skip car to the top hopper due to partial remixing, but strong segregation is still observed. Lump ore and nut coke exhibit the strongest segregation, while pellets remain the least segregated. The order of pellets and sinter in the weighing bunkers strongly influences their segregation patterns, whereas variations in the sinter particle size distribution (PSD) and particle shape have only limited effects. The insights from this study provide a basis for developing practical strategies to mitigate segregation in industrial BF charging.

Powder flowability underlies reliable solids handling, influencing dosing accuracy and production stability. Wood powders are usually cohesive and susceptible to flow problems like bridging because of their irregular, fibrous particles that are hygroscopic and heterogeneous. Two lignocellulosic powders were tested: spruce (softwood) and poplar (hardwood). Their particle size distribution, particle shape, and density were measured experimentally. Crucially, the interparticle parameters that govern powder bulk behavior, which are the cohesion energy density (CED), rolling friction coefficient (μᵣ), and sliding friction coefficient (μₛ), are not directly measurable at the scale and morphological complexity of fibrous wood particles. Therefore, using the Discrete Element Method (DEM), (μₛ,μᵣ, CED) were identified as effective DEM parameters by inverse calibration against rotating drum tests. A novel calibration workflow was developed to compare DEM simulations with real rotating drum experiment indicators, which can be used for unconfined, dynamic flow. These indicators correspond to newly discovered macroscopic flow descriptors that are processed from the powder bed: average projected area Area¯, its fluctuation σArea, and the average surface profile irregularity r²¯. Wood particles were modeled as multi-sphere clumps with different sizes to balance realism and computational cost. The calibrated parameters were: spruce—μ_s=0.10, μ_r=0.367, CED=130 kJ/m³; poplar—μ_s=0.10, μ_r=0.772, CED=100 kJ/m³. Following a comprehensive results analysis, increasing CED and friction parameters deteriorates powder unconfined flowability by promoting agglomeration and particle interlocking. The resulting calibrated DEM inputs provide a baseline for predicting and improving the handling of wood powders in hoppers, feeders, and conveying screws.

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.

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.

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.

Sodium borohydride (NaBH4) is increasingly considered as an alternative fuel for maritime vessels due to its relatively high energy density. When stored in dry solid form, it is a granular material, similar to coal, starch, and iron ore. As NaBH4 is historically used in the chemical industry in aqueous solutions, virtually no details regarding its behaviour as a solid granular material are known. Therefore, after determining particle properties such as size, shape, and density, this study characterises granular NaBH4 in three flow regimes using three experimental setups. Ring shear tests are used for the quasi-static regime, ledge tests for the dense flow regime, and rotating drum tests characterise both dense and gaseous flow, depending on the rotational speed. Various operational conditions, including temperature, humidity, time consolidation, and handling stresses, are taken into account. Experimental results demonstrate that above a threshold temperature and humidity, NaBH4 readily absorbs moisture from ambient air but remains free-flowing for most scenarios. However, time consolidation can transform this free-flowing material into a very cohesive substance. While this cohesiveness is reversible, requiring minimal agitation, the transformation from free-flowing to cohesive is accelerated by elevated moisture contents and a reduced particle size. Additionally, handling stresses were found to have minimal effect on the flow behaviour and characteristics of NaBH4. These findings are ultimately used to derive implications for the design of handling and storage equipment for NaBH4, enabling its use as an alternative fuel for maritime vessels.

Clay is a notoriously challenging material to dredge. Due to its adhesion and plastic behaviour, it may clog the suction head and/or clay balls could form in the pipeline. This will raise difficulties in estimating the production, the required power and increase the risk of downtime. As this is an expensive risk for the dredging industry, the cutting process requires in depth research to achieve better understanding of the process, to prevent problems and to mitigate risks. Available literature on clay deformation and soil cutting has been reviewed. Important topics for the cutting process are the interaction between the clay sliding over the blade and the resulting macro deformation of the chip. Various cutting regimes can be distinguished, including the: Flow regime, Tear regime, Curling regime, etc. Additionally the best practices for soil bin experiments have been included. Review of the available literature and analysis of published models is used to design a soil bin experiment dedicated to test the process under conditions relevant for the dredging industry. The objective of the CHiPS project is to study cutting regime transitions for dimensionless parameter groups of soil properties and operating conditions. Transitions range from static traction problems on soft mud to grinding action on stiff clay. Preliminary results and analysis of these clay cutting experiments are presented. The test rig developed for the CHiPS project is functionally performing satisfactorily, but requires a stronger drive to test high-strength soils.

Locomotion over granular terrain poses significant challenges for autonomous robotic systems, particularly in coastal regions characterized by loose, shifting sands. To optimize the locomotion on these challenging terrains, a simulation-aided design approach was used to develop a soft, shape-adapting, wheeled locomotion system. A co-simulation framework combining the discrete element method (DEM) and multibody dynamics (MBD) is employed to simulate the locomotion of a wheeled robot on varying sandy soils, covering both dry and wet sandy soil conditions. A shape-adapting wheel design is proposed, incorporating soft, inflatable elements that enable the wheel to transform between lugged and circular configurations. A discretized flexbody approach is adopted to model the interactions between the sandy soil and the soft, flexible bodies of the shape-adapting wheel design. Simulation results demonstrate improved performance of the shape-adapting wheels across a variety of sandy terrains, including slopes and obstacles. Integrating softness into the wheel improves obstacle climbing performance, while a lugged wheel configuration performs particularly well on loose, dry sandy slopes. This DEM-MBD co-simulation further enables efficient evaluation of locomotion strategies without the need for extensive physical prototyping.

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.

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.

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.

Mechanochemical synthesis faces reproducibility and scale-up challenges due to complex parameter interactions. This study employs machine learning (ML) to predict NaBH4 regeneration yield, integrating chemical experimental data and DEM (Discrete Element Method) derived invariant mechanical descriptors (Ēn, Ēt, fcol/nball). Various algorithms were evaluated, including a two-step modeling strategy to isolate the dominant effect of milling time in our process. Results demonstrate that a two-step Gaussian Process Regression (GPR) model achieves good predictive performance (R2 = 0.83), significantly outperforming single-stage models and providing valuable uncertainty estimates. Tree-based ensembles (XGBoost, RF) also benefit from the two-step approach and can enhance interpretability. This work establishes a framework for using ML to optimize mechanochemical processes, reducing experimental cost and offering a method to link mechanical milling conditions to chemical outcomes, thereby enabling predictive mechanochemistry.

The growth of offshore wind farms is accelerating to meet the renewable energy target by 2030, driving the development of larger offshore wind turbines (OWTs) to boost energy capacity. To support these OWTs, large monopiles are being installed by using impact hammers, which in turn emit low-frequency underwater noise, posing challenges for traditional noise mitigation systems and increasing risks to marine life. To address this, a metamaterial-based cushion (meta-cushion) was proposed, embedding resonators to filter longitudinal waves associated with high underwater noise levels. While prior work has demonstrated the meta-cushion's noise attenuation potential, design guidelines are required for adaptation to various monopile installations. This paper introduces, for the first time, a design methodology for the meta-cushion, which based on the input parameters of the monopile system, it details the procedure for selecting the resonant elements contributing to the attenuation performance and their spatial arrangement on the cushion for enhancing mechanical performance. Such performance indicators are evaluated via finite element simulations and experimental modal analyses. The methodology concludes with a nondimensional study of the spiral resonator, which showed the best attenuation response in experiments, exploring its behavior under varying material and geometric parameters. This methodology enables the development of meta-cushions adaptable to monopile installations under any environmental conditions.

Waterborne transport is very important for moving freight and passengers globally. To make this transport more efficient, vessel design must adapt to changing missions, regulations and the occurrence of malfunctions. This paper presents the design of an intelligent decision-support framework to assist marine engineers and vessel operators in updating the system and control architecture of marine vessels before and during a mission. The connection between the system architecture and control design perspectives is enabled using a semantics-based technique. To this end, the multi-level vessel control system is described by a semantic database, a knowledge graph used to connect the components automatically, and quantitative service criteria. Considering the system architecture, the optimal modification is deduced using modularity and complexity criteria, originating from the field of network theory. On the control side, an intelligent automation supervisor is designed to make offline and online decisions regarding the energy deficit to execute a new mission and the active automation configuration during operation. For offline decisions, system architecture modifications are requested by the vessel designers to cover the energy deficit. During operation, switching between hardware and virtual sensors as well as switching between energy management controllers is implemented to handle the effects of sensor faults. The framework is successfully applied to a case study of a tugboat used to adapt to missions with different power requirements, while simulation results are used to indicate its application in supporting the decisions of vessel designers and human vessel operators.

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.