Soft robotics requires structural systems capable of performing complex and programmable deformations to adapt to unstructured or dynamic environments. Shape memory materials (SMMs) offer a promising solution owing to their shape memory effect and stimulus-responsive adaptability. However, actuators relying on a single type of SMM are often constrained by nonlinear actuation behavior and limited stiffness variation, which restrict their ability to achieve coordinated, multifunctional responses. Addressing these challenges, this study introduces a hybrid programmable morphing structure that integrates a shape memory polymer (SMP) and a shape memory alloy (SMA) to realize cooperative actuation and adaptive stiffness variation within a single unit. In the proposed configuration, the SMA springs act as thermally activated actuators that generate deformation. The SMP cylindrical core employs its shape memory effect to realize reversible shape locking and serves as a thermal switch that enables controlled stiffness variation through temperature regulation. A coupled numerical model was established to describe the cooperative behavior between the SMA and SMP components, and the numerical results were validated through experimental testing. The agreement between simulations and experiments confirms the feasibility and repeatability of the proposed design. The structure achieves a maximum bending angle of 55° under dual-SMA actuation and 42° under single-SMA actuation, while maintaining any intermediate shape during thermal cycling. Furthermore, the hybrid system demonstrates a reversible six-fold increase in stiffness and a motion range extending up to three times its original length, representing a significant improvement over conventional single-material soft actuator. Moreover, the proposed hybrid structure offers a flexible strategy for programmable morphing and demonstrates scalable applicability in practical applications, such as adaptive grasping, reconfigurable locomotion, and environmental exploration. In conclusion, this work provides a feasible and generalizable framework for integrating multiple SMM into programmable morphing structures which can be applied into multifunctional soft robotic systems.

Sodium borohydride (NaBH₄) is a promising hydrogen carrier for maritime applications due to its high gravimetric and volumetric energy densities compared to compressed or liquefied hydrogen. As a solid granular material, NaBH₄ can be stored under atmospheric pressure and room temperature, eliminating the need for extreme pressures, or cryogenic conditions. Although NaBH₄ is hygroscopic and can become cohesive when exposed to humid environments, when stored and handled sufficiently dry it remains free-flowing, which is essential for reliable conveying, storage, and discharge operations. Designing equipment for these processes requires an accurate understanding of NaBH₄’s free-flowing behaviour, motivating the need for modelling tools. To address this, we model granular, free-flowing NaBH₄ using the Discrete Element Method. Granular flow behaviour depends strongly on the flow regime—defined by the applied shear and confining pressure through the inertial number—and practical handling systems typically operate in the dense flow regime; consequently, this study focuses on dense-flow conditions. A standard calibration–verification–validation methodology is applied, using a ledge test and rotating drum for calibration and verification, and an inclined surface test for validation. In parallel, we evaluate the inertial number across all setups using a velocity-based approach to characterise both the global flow regime and locally occurring flow regimes within the flowing layer. Rather than assigning a single characteristic value per setup, we demonstrate that distinct segments exhibit different inertial numbers, indicating that the inertial number is spatially dependent.

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.

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.

This study introduces a computational framework for modelling raw chicken breast fillets using the Discrete Element Method (DEM), aimed at providing a baseline efficient simulation model for large-scale poultry handling processes. A bonded multi-sphere meta-particle representation was developed and calibrated through mechanical testing of raw fillets. Compression experiments yielded a Young’s modulus of approximately 48.6 kPa, which informed the stiffness properties of the DEM sub-particle assembly. Numerical Design of Experiments (DoEs) highlighted the need for an unbalanced ratio between normal and shear bond stiffness to ensure correct damping behaviour and preserve realistic flexibility. The framework was validated using a full-scale hopper–conveyor discharge experiment, demonstrating the model’s ability to reproduce key physical behaviours such as large deformations, curling during discharge, and the transition between jammed and free-flow regimes. The simulation closely matched the measured discharge rate, with all chicken fillets discharged within 4 s at a 6 cm gate opening height. The proposed model required approximately 9 mins to simulate a 10-second industrial-scale process, underscoring the model’s practical suitability for simulation-aided design and optimisation of poultry processing 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.

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.

We compare the influence of tangential (shear) and normal (compressive) stress events on the mechanochemical regeneration of sodium borohydride NaBH₄ from hydrated sodium metaborate [Figure presented] and magnesium hydride MgH₂. Discrete element method (DEM) mechanical descriptors are used to design experiments that either maintain the mill at a constant rotational speed or maintain a constant total dissipation power, thereby separating stress distribution from net power input. Under constant power operation, a tangential rich regime achieves a record 94% conversion yield with 37.5% shorter milling time, 40% lower ball-to-powder ratio, and 34% reduced speed. However, this high yield requires such a substantial power consumption that the converted mass per Watt drops to only 0.090 gW⁻¹, below both balanced (0.113 gW⁻¹) and normal-bias (0.108 gW⁻¹) cases. By contrast, a tangential bias at half the power (3 W) still delivers 84% yield and peaks at 0.185 gW⁻¹, illustrating the often disregarded trade-off between absolute conversion and energetic productivity in mechanochemistry. Specific yield (conversion per Watt) likewise peaks at 0.28 W⁻¹ and declines linearly with fill ratio (R²>0.99). Mechanochemical energy leverage analysis reveals that, at most, 1.7–3.7% of input mechanical work is theoretically recoverable on an enthalpy basis, 2.1–4.4% on a Gibbs free energy basis, and 4–8.7% when considering the fuel value of all available hydrogen. Our mill-agnostic framework provides a transferable blueprint for cross-platform optimization of mechanochemical processes.

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.

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.

Micro-scale mechanisms, such as inter-particle and particle-fluid interactions, govern the behaviour of granular systems. While particle-scale simulations provide detailed insights into these interactions, their computational cost is often prohibitive. At a recent Lorentz Center Workshop on “Machine Learning for Discrete Granular Media”, researchers explored how machine learning approaches can aid the development of constitutive laws and efficient data-driven surrogates for granular materials while also addressing uncertainty quantification. Attended by researchers from both the granular materials (GM) and machine learning (ML) communities, the workshop brought the ML community up to date with GM challenges. This position paper emerged from the workshop discussions. In this position paper, we define granular materials and identify seven key challenges that characterise their distinctive behaviour across various scales and regimes–ranging from gas-like to fluid-like and solid-like. Addressing these challenges is essential for developing robust and efficient models for the digital twinning of granular systems in various industrial applications. To showcase the potential of ML to the GM community, we present classical and emerging machine/deep learning techniques that have been, or could be, applied to granular materials. We reviewed sequence-based learning models for path-dependent constitutive behaviour, followed by encoder-decoder type models for representing high-dimensional data in reduced spaces. We then explore graph neural networks and recent advances in neural operator learning. The latter captures the emerging field evolution of interacting particles via efficient latent space representation. Lastly, we discuss model-order reduction and probabilistic learning techniques for high-dimensional parameterised systems, both of which are crucial for quantifying and incorporating uncertainties arising from physics-based and data-driven models. We present a typical workflow aimed at unifying data structures and modelling pipelines and guiding readers through the selection, training, and deployment of ML surrogates for granular material simulations. Finally, we illustrate the workflow’s practical use with two representative examples, focusing on granular materials in solid-like and fluid-like regimes.