One of the most promising energy carriers for transport applications are hydrogen-based energy carriers. NaBH₄ is a hydrogen energy carrier and produces hydrogen bubbles when it is dissolved in water. The formation of hydrogen bubbles hinders experimental measurements of the thermodynamic and transport properties of aqueous NaBH₄ solutions at elevated temperatures. Accurate knowledge of these properties is essential for the NaBH₄ hydrolysis reactor modeling and design. Molecular dynamics (MD) simulations provide the option to study the thermodynamic and transport properties of NaBH₄ aqueous solutions without hindering hydrogen bubble formation. In this work, a new force field is developed for BH₄^-, namely, the Delft force field of BH₄^- (DFF/BH₄^-), which, combined with additional force fields, can accurately describe experimental densities and viscosities of 0 to 5 m (mol salt/kg water) NaBH₄, 0 to 3 m NaB(OH)₄, and 1 m NaOH aqueous solutions at 295 K within 1.8% and 10.8% maximum deviation, respectively. Empirical fitting correlations are created for densities, viscosities, and self-diffusivities obtained from the MD simulations of 0 to 5 m NaBH₄, 0 to 5 m NaB(OH)₄, and 0 to 1 m NaOH aqueous solutions at 295-363 K for NaBH₄ hydrolysis reactor modeling and design purposes.

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.

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.

Hydrogen carriers are attractive alternative fuels for the shipping sector. They are zero-emission, have high energy densities, and are safe, available, and easy to handle. Sodium borohydride, potassium borohydride, dibenzyltoluene, n-ethylcarbazole, and ammoniaborane are hydrogen carriers with high theoretical energy densities. The energy density is paramount to implementing hydrogen carriers as a high energy density enables compact and lightweight storage. The effective energy density depends on integrating heat and masses with energy converters. This combination defines the energy efficiency and, thus, the energy density of the system. This paper addresses the effective energy density of the hydrogen carriers, including the dehydrogenation process. Using a 0D model, we combined the five carriers with two types of fuel cells, namely proton exchange membrane (PEM) and solid oxide fuel cells (SOFC), an internal combustion engine and a gas turbine. N-ethylcarbazole and dibenzyltoluene offer medium energy densities, reaching almost 4 MJ/kg. However, the effective energy density of sodium borohydride and ammoniaborane is very high, up to 15 MJ/kg, including the energy converter. This is similar to the energy density of marine diesel oil combined with an internal combustion engine. Thus, we conclude hydrogen carriers are alternative fuels that deserve more attention because of their strong potential to make shipping zero-emission.

Electrochemical reactors, such as water electrolyzers, CO₂ electrolyzers, fuel cells, and flow batteries, will be essential in electrifying industry as part of the global transition towards a defossilized and sustainable economy. These technologies require further optimization to enhance efficiency and reduce costs for widespread adoption. Hydrodynamics and mass transfer at electrode–electrolyte interfaces significantly affect electrochemical conversion reactions by influencing the reactant availability and pH in the local reaction environment. 3D electrodes, such as flow-through foams and suspension electrodes, hold a great advantage over 2D electrodes as they moderate pH changes and reactant depletion by spreading the current over a larger electrode area and electrolyte volume. We study the diffusion boundary layer in operando around a single mm-sized particle, representing an element of a 3D electrode. We visualize the local and transient pH with Fluorescence Lifetime Imaging Microscopy (FLIM) during H₂O reduction at various current densities and electrolyte flow velocities at a resolution down to 9 μm and 2 Hz. In addition, we apply an intermittent current to investigate how long the capacitive electric double layer of a suspension electrode particle can maintain an electrochemical reaction during their time of non-contact with a current collector, mimicking applications with Faradaic charge transfer (i.e. flow batteries, microbial fuel cells, capacitance-based electrolyzers). We demonstrate that the diffusion boundary layer is not symmetrical, but depend on the direction of the electric field, the current density and the flow conditions. The substantial pH gradients and boundary layer formation at the scale of hundreds of micrometers underline the importance of controlling flow in or around electrodes, making 3D electrodes an important asset for creating suitable reaction conditions in mass transport-limited electrochemical conversions.

– Hydrogen carriers, such as liquid organic hydrogen carriers (LOHCs) and borohydrides, are promising zero-emission alternative fuels for ships. Bringing these hydrogen carriers on board, however, creates new challenges. A major challenge is their spill behaviour. Knowing the spill behaviour is paramount to avoid large-scale environmental disasters. This paper investigates the spill behaviour of four hydrogen carriers (and their conjugates): sodium borohydride, ammonia borane, dibenzyltoluene, and n-ethylcarbazole. The hydrogen carriers were all dissolved in artificial seawater to test their behaviour. Sodium borohydride reacts with seawater, as it also reacts with pure water. However, contrary to expectations, it reacts faster with seawater than regular water. The reaction mechanism behind this is unknown. Ammonia borane does not visibly react with normal water or with seawater. Dibenzyltoluene sinks and forms tiny bubbles which are easily perturbed. Unfortunately, perhydro dibenzyltoluene could not be tested due to technical problems. N-ethylcarbazole breaks up into smaller pieces and predominantly stays afloat, likely due to the surface tension of water. Perhydro n-ethylcarbazole floats but is barely visible in seawater due to its transparency. Preventive measures must be established to avoid large-scale spills if these substances are utilised on ships, as they are likely challenging to clean up.

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.

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).

Hydrogen carriers are attractive alternative fuels for the shipping sectors. They are zero-emission, have high energy densities, and are safe, available, and easy to handle. Sodium borohydride, potassium borohydride, dibenzyltoluene, n-ethylcarbazole, and ammoniaborane are interesting hydrogen carriers, with high theoretical energy densities. The exact energy density of these hydrogen carriers depends on the integration of heat and mass with the energy converters. This combination defines the energy efficiency and, thus, the energy density of the system. Using a 0D model, we combined the five carriers with two types of fuel cells (PEM and SOFC), an internal combustion engine and a gas turbine. This resulted in 20 combinations. Despite the limitations of the 0D model and the occasional difficulty of validating input values, this model still produces exciting findings, which are valuable for further research. For the dehydrogenation of both dibenzyltoluene and n-ethylcarbazole, an external hydrogen burner is required if no waste heat resources from the integrated system are available. For the borohydrides, on the other hand, energy integration is essential for reducing cooling power. Dehydrogenation produces substantial energy, but only a fraction of this energy can be used for internal preheating. Dehydrogenation of ammoniaborane produces less energy. Among all hydrogen carriers, both ammoniaborane and sodium borohydride provide energy densities comparable to that of marine diesel oil. In particular, ammoniaborane possesses a remarkably high energy density. Thus, we conclude, that hydrogen carriers are attractive alternative fuels that deserve more attention, including their potential performance for hydrogen imports.

CO2 conversion is an important part of the transition towards clean fuels and chemicals. However, low solubility of CO2 in water and its slow diffusion cause mass transfer limitations in aqueous electrochemical CO2 reduction. This significantly limits the partial current densities towards any desired CO2-reduction product. We propose using flowable suspension electrodes to spread the current over a larger volume and alleviate mass transfer limitations, which could allow high partial current densities for CO2 conversion even in aqueous environments. To identify the requirements for a well-performing suspension electrode, we use a transmission line model to simulate the local electric and ionic current distributions throughout a channel and show that the electrocatalysis is best distributed over the catholyte volume when the electric, ionic and charge transfer resistances are balanced. In addition, we used electrochemical impedance spectroscopy to measure the different resistance contributions and correlated the results with rheology measurements to show that particle size and shape impact the ever-present trade-off between conductivity and flowability. We combine the modelling and experimental results to evaluate which carbon type is most suitable for use in a suspension electrode for CO2 reduction, and predict a good reaction distribution throughout activated carbon and carbon black suspensions. Finally, we tested several suspension electrodes in a CO2 electrolyzer. Even though mass transport limitations should be reduced, the CO partial current densities are capped at 2.8 mA cm−2, which may be due to engineering limitations. We conclude that using suspension electrodes is challenging for sensitive reactions like CO2 reduction, and may be more suitable for use in other electrochemical conversion reactions suffering from mass transfer limitations that are less affected by competing reactions and contaminations.

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).

Oscillatory flow in confined spaces is central to understanding physiological flows and rational design of synthetic periodic-actuation based micromachines. Using theory and experiments on oscillating flows generated through a laser-induced cavitation bubble, we associate the dynamic bubble size (fluid velocity) and bubble lifetime to the laser energy supplied—a control parameter in experiments. Employing different channel cross-section shapes, sizes and lengths, we demonstrate the characteristic scales for velocity, time and energy to depend solely on the channel geometry. Contrary to the generally assumed absence of instability in low Reynolds number flows (<1000), we report a momentary flow distortion that originates due to the boundary layer separation near channel walls during flow deceleration. The emergence of distorted laminar states is characterized using two stages. First the conditions for the onset of instabilities is analyzed using the Reynolds number and Womersley number for oscillating flows. Second the growth and the ability of an instability to prevail is analyzed using the convective time scale of the flow. Our findings inform rational design of microsystems leveraging pulsatile flows via cavitation-powered microactuation.

Background: Electrohydrodynamic (EHD) drying relies on the generation of a corona wind that is created with a high electric potential between an emitter and a ground electrode. The impinging corona wind enhances convective heat and mass transport close to the sample, which is positioned on the ground electrode. Although EHD drying is reported to have large potential for milder and energy efficient drying, most studies have only evaluated its potential at the lab scale and for specific applications. Scope and approach: This review focusses on discussing the underlying physical phenomena that may be expected to influence the performance of EHD drying processes. Specifically, we discuss how corona wind changes due to discharge behaviour as it is influenced by moisture release during drying, how the microstructural properties of the product could influence electromigration in the drying material and how EHD drying could be efficiently and safely scaled up for biobased & food products. We conclude with an outlook to the research that still is required to concretize the potential of EHD drying in practice. Key findings and conclusions: In EHD drying, electrical discharge induces airflow which enhances external convective drying. Studies coupling physical phenomena are still lacking, although some possible couplings are mentioned in chapter 2. Internal mass transfer is enhanced due to presence of an electrical potential over the product, which can cause several electrically-driven transport phenomena as discussed in chapter 3. Several drying setups have been explored in industry, of which the characteristics are summarized in chapter 4.

The local conditions inside a gas diffusion electrode (GDE) pore, especially in the electrical double layer (EDL) region, influence the charge transfer reactions and the selectivity of desired CO₂ER products. Most GDE computational models ignore the EDL or are limited in their applicability at high potentials. In this work, we present a continuum model to describe the local environment inside a catalytic pore at varying potentials, electrolyte concentrations and pore diameters. The systems studied in this work are based on an Ag catalyst in contact with KHCO₃ solution. Our study shows that steric effects dominate the local environment at high cathodic potentials (≪−25 mV vs pzc at the OHP), leading to a radial drop of CO₂ concentration. We also observe a drop in pH value within 1 nm of the reaction plane due to electrostatic repulsion and attraction of OH⁻ and H⁺ ions, respectively. We studied the influence of pore radii (1-10 nm) on electric field and concentrations. Pores with a radius smaller than 5 nm show a higher mean potential, which lowers the mean CO₂ concentration. Pores with a favourable local environment can be designed by regulating the ratio between the pore radius and Debye length.

Electrochemical oxygen reduction is a promising and sustainable alternative to the current industrial production method for hydrogen peroxide (H₂O₂), which is a green oxidant in many (emerging) applications in the chemical industry, water treatment, and fuel cells. Low solubility of O₂ in water causes severe mass transfer limitations and loss of H₂O₂ selectivity at industrially relevant current densities, complicating the development of practical-scale electrochemical H₂O₂ synthesis systems. We tested a flow-by and flow-through configuration and suspension electrodes in an electrochemical flow cell to investigate the influence of electrode configuration and flow conditions on mass transfer and H₂O₂ production. We monitored the H₂O₂ production using Cu-tmpa (tmpa = tris(2-pyridylmethyl)amine) as a homogeneous copper-based catalyst in a pH-neutral phosphate buffer during 1 h of catalysis and estimated the limiting current density from CV scans. We achieve the highest H₂O₂ production and a 15-20 times higher geometrical limiting current density in the flow-through configuration compared to the flow-by configuration due to the increased surface area and foam structure that improved mass transfer. The activated carbon (AC) material in suspension electrodes, which have an even larger surface area, decomposes all produced H₂O₂ and proves unsuitable for H₂O₂ synthesis. Although the mass transfer limitations seem to be alleviated on the microscale in the flow-through system, the high O₂ consumption and H₂O₂ production cause challenges in maintaining the initially reached current density and Faradaic efficiency (FE). The decreasing ratio between the concentrations of the O₂ and H₂O₂ in the bulk electrolyte will likely pose a challenge when proceeding to larger systems with longer electrodes. Tuning the reactor design and operating conditions will be essential in maximizing the FE and current density.

In flow-by capacitive deionization (CDI) brackish water flows between two electrodes that capacitively remove salt. We assume low inlet concentrations so “salt shocks” appear in the electrodes and the process becomes diffusion-limited. For unit charge efficiency, a simplified model is derived consisting of two coupled partial differential equations. We obtain approximations, and exact solutions in terms of the Lambert W function, for the salt concentration as a function of time and space and for the equilibrium charge-voltage relation. These surprisingly simple solutions compare well with the results from comprehensive two-dimensional simulations. Useful analytical expressions are obtained for optimal geometrical and operational parameters that maximize the productivity and minimize the specific energy losses. By making cells much thinner the productivity can be increased an order of magnitude compared to typical values in the literature. The optimal electrode is found to be roughly six times thinner than the spacer. The associated pressure drop is around 0.4 bar per 1 mM of inlet salt concentration, making our recommendations practically feasible only for relatively low concentrations. The obtained model and analytical expressions provide useful guidance to strongly improve the design process.