Fluid-solid Interface Research Articles

Custom lens design for ultrasonic inspection of immersed anisotropic parts

Inspection of anisotropic parts using immersion ultrasonic waves poses resolution challenges caused by directionally dependent wave properties. Currently, focused lenses with spherical or cylindrical geometries are the standard method for creating convergent, small cross-section, and high amplitude beams within solids. These lenses form a conical beam that uniformly refracts at the fluid-solid interface of the immersed part when the solid is isotropic, producing the desired focusing behavior. However, for anisotropic solids, refraction becomes directionally dependent, which distorts the beam and causes non-uniform focal profiles and, in some cases, multiple foci. In this work, the forward wave propagation for an immersed anisotropic solid is modeled and exploited to develop custom lens geometries. These custom lenses utilize inverted slowness profiles to produce conical beams within the anisotropic solid. By incorporating these lenses, the focused beam behavior is improved and becomes more circular, higher in amplitude, and smaller in the cross-sectional area. These lenses have the potential to improve inspection resolution in anisotropic solids, an effort beneficial to the additive, aerospace, and electronics industries.

Coupling of multicomponent transport models in particle-resolved fluid-solid simulations

An accurate and self-consistent methodology for mass transport of multi-component mixtures in multi phase media is a necessity for a proper description of complex physical and chemical processes in reactors such as catalytic packed beds. In this regard, a novel methodology has been developed to describe and couple underlying transport phenomena in fluid and porous media as well as at the solid-fluid interface. The methodology is symmetric as it treats all components in a mixture equally. The Maxwell-Stefan equations are symmetrically formulated, discretized conservatively and coupled with a compressible flow solver for the fluid part. The Dusty Gas Model is applied inside porous media by developing a self-consistent and robust numerical formulation. A ghost-cell Immersed Boundary Method is used to capture the physics at the solid-fluid interface with the implementation of a novel symmetric non-singular mass flux formulation. Several test cases are established to demonstrate the accuracy and robustness of the newly developed symmetric methodology in this paper. These test cases can be used as benchmark for the future development of symmetric methodologies for multicomponent systems in multi phase media.

Modeling of thermal radiation in a participating media with conjugate heat transfer

AbstractModeling a combination of thermal radiation and conjugate heat transfer in a three‐dimensional rectangular domain which has a participating media CO2 flowing through is done numerically in OpenFOAM. The rectangular duct has a vertical step (facing forward to the inlet) which is located at a distance from the inlet (the distance is same as the height of the inlet section). The domain is divided into two regions (namely solid and fluid). Carbon dioxide, a highly absorbing fluid with extinction, is used here as the participating medium. The ability of the code is verified to analyze the thermal radiation in a participating media with conjugate heat transfer. The study was carried out for a constant Reynolds number 250 and a contraction ratio of 0.5. The study focused primarily on the importance of adding thermal radiation on to thermal analysis and the reason behind the Nusselt number variation on different regions of solid–fluid interface. It also discussed the effect of radiative properties, such as optical thickness and linear scattering albedo, on the average convective Nusselt Number.

Thermal boundary conditions at the fluid–solid interface in the case of a conducting body: a novel thermal lattice Boltzmann analysis

Thermal boundary conditions at the fluid–solid interface in the case of a conducting body: a novel thermal lattice Boltzmann analysis

Important of Slip Effects in Non-Newtonian Nanofluid Flow with Heat Generation for Enhanced Heat Transfer Devices

In various fields such as engineering, nanotechnology, and biomedical sciences, the study of non-Newtonian nanofluid flow with heat generation is becoming increasingly important. However, it is challenging to accurately model such flows due to their complex behavior and slip effects at the fluid-solid interface. This research investigates the impact of first and second-order slip conditions on the flow and heat transfer properties of a non-Newtonian nanofluid using a power law model to describe the fluid's non-Newtonian behavior and numerical methods to solve the resulting equations. To determine the influence of various parameters such as slip parameters, Brinkman number, power law index, and Eckert number on the velocity, temperature, and concentration profiles, which this study examines. The study shows that slip parameters significantly determine the flow and heat transfer properties of non-Newtonian nanofluids, the study also reveals that slip parameters are a crucial factor in understanding the flow and heat transfer characteristics of nanofluids, with the second-order slip condition having a greater impact on velocity and temperature profiles than the first-order slip condition. These findings are valuable for developing and optimizing heat transfer devices that involve non-Newtonian nanofluids with heat generation, which is essential for technological advancements in today's industry.

Microfluidic organ chip of fluid–solid dynamic curved interface

Dynamic curved interfaces are fundamental and ubiquitous structures in biological systems. However, replicating the structure and function associated with these interfaces for mechanobiology and drug screening is challenging. Here, we develop a dynamic curvature-enabled microfluidic organ chip of two fluid–solid dynamic curved interfaces. One interface effectively integrates adjustable biomechanics, and the other controls drug release with open microfluidics. The fluid–solid interface sensed by the cells can modulate the residual stress, stiffness, strain of the solid phase, and the flow shear stress of the fluid phase. Using the chip, we investigate the mechanotransductive responses of endothelial and epithelial cells, including Piezo1, Ca2+, and YAP, and reveal that the response of the endothelium to combined dynamic cyclic strain and flow shear stress is different from separate stimulation and also disparate from the epithelium. Furthermore, direct and high-efficiency drug release to cells is realized by constructing the other fluid–solid interface on the back side of cells, where drugs are encapsulated within cross-linked alginate hydrogel in the open microfluidic channel. Then, we replicate object-specific and location-specific biomechanical environments within carotid bifurcation and prove the effectiveness of drug delivery. Our design exemplifies dynamic curved biological interfaces with controlled mechanical environments and holds potential for patient-specific medicine.

An accurate and scalable direction-splitting solver for flows laden with non-spherical rigid bodies—Part 2: Moving rigid bodies

The Direction Splitting algorithm is a computational method that facilitates the solution of the Navier–Stokes equation by decomposing a multi-dimensional Poisson equation into a series of one-dimensional equations. This approach allows for efficient and scalable numerical computations. The incorporation of rigid bodies in the computational domain is made possible by modifying the diffusion and divergence stencils near the fluid–solid interface. In our recent study (Goyal and Wachs, 2023), we implemented stencil corrections for fixed non-spherical rigid bodies and demonstrated the scalability and accuracy of the DS algorithm for a wide range of test cases. This study presents an extension of the Direction Splitting solver to effectively compute the dynamics of moving non-spherical rigid bodies in the presence of fluid flow. The interaction between rigid bodies and fluid can be categorized as either one-way or two-way coupled. The stencils resolving the fluid–solid interface are updated at every time-iteration to account for changes in particle position and orientation. We demonstrate that the adaptation of the diffusion and divergence stencil in the literature (Goyal and Wachs (2023); Keating and Minev (2013)) is inadequate for accurately predicting the dynamics of moving rigid bodies in a two-way coupled simulation, especially at O(102) Reynolds number. Hence, we implement a conservative cut-cell method for the discretization of the divergence term, and we also take into account the proximity of the fluid–solid interface in the advection discretization. The results of multiple test cases demonstrate that the suggested stencil modifications effectively capture the dynamics of fluid flow in the presence of freely moving rigid bodies. However, the cut-cell discretization of the divergence term constitutes a significant computing overhead that impacts the fast computing property of the Direction Splitting algorithm.

A conjugated heat and mass transfer model to implement reaction in particle-resolved CFD simulations of catalytic fixed bed reactors

Modelling catalytic fixed bed reactors with a small tube-to-particle diameter ratio requires a detailed description of the interactions between fluid flow, intra-particle transport, and the chemical reaction(s) within the catalyst. Particle-resolved computational fluid dynamics (PRCFD) simulations are the most promising approach to predict the behaviour of these reactors accurately, since they take into account the local packed bed structure explicitly. In this work, a conjugated heat and mass transfer model for use in PRCFD simulations is presented in order to couple the fluid flow through the fixed bed with transport and reaction in the porous catalyst, while guaranteeing the no-slip boundary condition at the fluid–solid interface. For this purpose, the solutions of the solid and fluid domain are computed separately and are coupled by calculation and updating the boundary condition at the particle surface. Owing to the consideration of secondary gradients, the developed transfer model is also valid for unstructured calculation meshes containing non-orthogonal cells at the fluid–solid interface. Such meshes are often used to resolve complex geometries, such as a packed bed, in a computationally efficient manner. The coupling approach is validated using cases for which an analytical solution or literature correlations derived from experimental data are available. The simulation results of a short catalytic packed bed with rings catalysing the partial oxidation of n-butane to maleic anhydride exemplify the potential of PRCFD involving reactions to analyse the catalyst performance in great detail.

NG-DPSM: A neural green-distributed point source method for modelling ultrasonic field emission near fluid-solid interface using physics informed neural network

Distributed point source method (DPSM) is a collocation point-based semi-analytical method to model the scattered ultrasonic fields in isotropic and anisotropic materials. It requires the evaluation of Green's function solutions and its rigorous differentiation for complex geometry problems having a large set of distributed point sources. DPSM is much faster than the finite element method (FEM) however, it is analytically demanding as it requires differentiation of the Green's function. The intricacy associated with differentiation operation hinders the potential application of DPSM in large-scale automation and real-time structural health monitoring. The current paper introduces machine-learning-based strategies within DPSM and yields the proposed Neural Green-DPSM (NG-DPSM) framework. A physics-informed neural network (PINN) is incorporated in the DPSM framework for evaluating the Green's function and its gradients. To train the PINN, displacement Green's function solutions of wave propagation in isotropic solid is used as the governing physics. Once the PINN is trained for displacement Green's function solutions, the NG-DPSM framework leverages the trained network and its automatic differentiation capability to predict displacement and stress fields annihilating the rigorous differentiation of Green's functions. The accuracy and efficacy of NG-DPSM are demonstrated using two numerical experiments. First, the ultrasonic fields are evaluated for the problem geometry of a plate immersed in fluid excited at different angles of incidence. Further, the efficacy of the proposed approach is demonstrated for wave scattering through a circular hole in the plate. The results show a strong agreement between the ultrasonic fields computed using both NG-DPSM and conventional DPSM.

Predicting fluid–structure interaction with graph neural networks

We present a rotation equivariant, quasi-monolithic graph neural network framework for the reduced-order modeling (ROM) of fluid–structure interaction systems. With the aid of an arbitrary Lagrangian–Eulerian (ALE) formulation, the system states are evolved temporally with two sub-networks. The movement of the mesh is reduced to the evolution of several coefficients via complex-valued proper orthogonal decomposition (POD), and the prediction of these coefficients over time is handled by a single multi-layer perceptron (MLP). A finite element-inspired hypergraph neural network is employed to predict the evolution of the fluid state based on the state of the whole system. The structural state is implicitly modeled by the movement of the mesh on the solid–fluid interface; hence, it makes the proposed framework quasi-monolithic. The effectiveness of the proposed framework is assessed on two prototypical fluid–structure systems, namely, the flow around an elastically mounted cylinder and the flow around a hyperelastic plate attached to a fixed cylinder. The proposed framework tracks the interface description and provides stable and accurate system state predictions during roll-out for at least 2000 time steps and even demonstrates some capability in self-correcting erroneous predictions. The proposed framework also enables direct calculation of the lift and drag forces using the predicted fluid and mesh states, in contrast to existing convolution-based architectures. The proposed reduced-order model via the graph neural network has implications for the development of physics-based digital twins concerning moving boundaries and fluid–structure interactions.

Numerical study of gas–surface interface effects due to transpiration in a hypersonic flow over a blunt body

This paper investigates the impact of transpiration on a hypersonic flow over a cylinder, considering different degrees of rarefaction. The study analyzes the interaction between freestream argon gas flow at Mach 5 and transpiring argon gas at the fluid–solid interface at a velocity of 10 m/s. Freestream Knudsen numbers considered are 0.002, 0.01, 0.05, and 0.25, spanning from a continuum to rarefied regime. Flow simulations utilize the open-source direct simulation Monte Carlo solver, Stochastic PArallel Rarefied-gas Time-accurate Analyzer. The influence of transpiration on flow and surface properties is examined by comparing non-transpiration and transpiration cases. At all regimes, transpiration increases the normal shock stand-off distance, while a comparison of flow properties along the stagnation line reveals a reduction in the velocity and an increase in the post-shock temperature with transpiration. Surface heat flux comparison indicates that transpiring gas reduces heat flux on the cylinder's upstream-facing front surface at all Knudsen numbers. However, at Kn∞ = 0.25, a shift occurs, and surface heat flux starts increasing locally from the top/bottom point on the cylinder surface through the rear face of the cylinder. Furthermore, a test for the validity of the continuum-based blowing correction correlation function reveals the failure of the empirical model, even in the continuum regime at Kn∞ = 0.002, casting doubt on its applicability to vehicles with curvilinear blunt-body shapes. Additionally, a sensitivity analysis demonstrates that transpiring gas with a number density an order of magnitude higher than the freestream reduces stagnation peak heat flux by nearly 30%, while transpiring gas with a temperature two times higher than the freestream shows a ∼13% reduction.

Implementation of a level-set-based volume penalization method for solving fluid flows around bluff bodies

A volume penalization-based immersed boundary technique is developed and thoroughly validated for fluid flow problems, specifically flow over bluff bodies. The proposed algorithm has been implemented in an open source field operation and manipulation (OpenFOAM), a computational fluid dynamics solver. The immersed boundary method offers the advantage of inserting a complex solid object inside a Cartesian grid system, and therefore, the governing equations can be applied to such a simpler grid arrangement. For capturing the fluid–solid interface more accurately, the grid is refined near the solid surface using topoSetDict and refineMeshDict utilities in OpenFOAM. In order to avoid any numerical oscillation and to compute the gradients accurately near the interface, the present volume penalization method (VPM) is integrated with a signed distance function, which is also referred to as a level-set function. Benchmark problems, such as flows around a cylinder and a sphere, are considered and thoroughly validated with the results available in the literature. For the flow over a stationary cylinder, the Reynolds number is varied so that it covers from a steady two-dimensional flow to an unsteady three-dimensional flow. The capability of the present solver has been further verified by considering the flow past a vibrating cylinder in the cross-stream direction. In addition, a flow over a sphere, which is inherently three-dimensional due to its geometrical shape, is validated in both steady and unsteady regimes. The results obtained by the present VPM show good agreement with those obtained by a body-fitted grid using the same numerical scheme as that of the VPM, and also with those reported in the literature. The present results indicate that the VPM-based immersed boundary technique can be widely applicable to scientific and engineering problems involving flow past stationary and moving bluff bodies of arbitrary geometry.

Mechanism of boundary bubble drag reduction of Couette flow in nano-confined domain

Bubble drag reduction technology is of great significance in improving the propulsion efficiency of underwater vehicle and reducing the comprehensive energy consumption during navigation. Bubble drag reduction is a highly effective method of reducing the frictional resistance encountered by large ships and underwater vehicles during navigation. It exhibits excellent stability in drag reduction, and has advantages such as environmental friendliness, adaptability to various flow environments, and suitability for all underwater components of ships. Therefore, it is greatly significant to conduct in-depth research on bubble drag reduction and its underlying mechanism. In this work, the flow characteristics and the boundary bubble drag reduction mechanism of gas-liquid Couette flow in parallel wall nanochannels are studied by molecular dynamics method, and the influences of surface wettability, wall roughness, and gas concentration on boundary slip velocity and bubble drag reduction effect are analyzed. The results indicate that the bubble drag reduction effect is enhanced with the increase of boundary slip velocity. In the gas-liquid two-phase flow region, with the increase of shear velocity, the lateral deformation of boundary adsorbed bubble and boundary slip velocity increase, thus enhancing the bubble drag reduction effect. The increase of solid-gas interaction strength and gas concentration can lead to the enrichment of gas atoms near the wall, improve the bubble spreading characteristics on the wall, and thus increase the slip velocity of the solid-liquid interface. The wall roughness can change the spreading characteristics of bubble, affect the boundary slip velocity, and then change the drag reduction effect of the fluid-solid interface. As the rib height increases, gas atoms accumulate in the grooves between ribs and the adsorption quantity of gas atoms on the upper surface of the rib decreases, which leads to the decrease of the boundary slip velocity of the solid-liquid interface and ultimately reduces the drag reduction effect. The research results will provide important theoretical guidance for implementing the boundary drag reduction technology in large ships and underwater vehicles.

A NUMERICAL APPROACH FOR THE STUDY OF HEAT GENERATION IN THE PRESENCE OF THERMAL BOUNDARY LAYER FOR A FLAT PLATE

In this study, we investigate the laminar boundary layer flow in two dimensions, steadiness, and incompressibility around a moving vertical flat plate in a uniform free stream of fluid with a convective surface boundary condition. The similarity transformation technique has been applied to convert the governing nonlinear partial differential equation into two nonlinear ordinary differential equations. By combining the finite difference method with the shooting technique, the problem is solved numerically. We present a tabular and graphical representation of the variation in dimensionless temperature and fluid-solid interface characteristics for various values of the Prandtl number. As a special case of the problem, a comparison between the current result and the previously published result demonstrates a good agreement.

Experimental and Numerical Study of Novel Vortex Bladeless Wind Turbine with an Economic Feasibility Analysis and Investigation of Environmental Benefits

This study combines experimental and numerical evaluations of Vortex Bladeless Wind Turbines (VBWTs) to understand their potential in renewable energy generation. The methodology employs Two-Way Fluid–Solid Interface (FSI) simulations, alongside real-world data, providing important insights into the turbine’s vibration dynamics and flow interactions during operation. Key findings include identifying optimal vibration frequencies and amplitudes that enhance energy harvesting and a clear advantage in power-generation estimations shown by one of the models used. The study reveals possible applications of VBWT in various settings like airport runways, highways, and buildings, indicating a promising avenue for incorporating such renewable-energy solutions. Discussions on the economic feasibility and environmental benefits of VBWT deployment are also presented, suggesting a need for further research and optimization in this area. A conceptual generator design and business model are introduced as part of a broader discussion on technology integration and energy storage. The research in this study encompasses experimental and numerical analysis, to achieve a broader understanding of the workings of a VBWT, realizing the feasibility of using such systems in lower-wind-speed conditions and upscaling to higher-wind-speed cases.

Anisotropic continuum framework of coupled gas flow – adsorption – deformation in sedimentary rocks

AbstractSolid deformation is always a crucial factor of gas transport in sedimentary rocks. While previous studies always adopt the assumption of isotropic poroelastic deformation, anisotropic poroelastoplastic deformation is rarely considered, despite anisotropy being a ubiquitous property of natural sedimentary rocks. In this work, an anisotropic poromechanical model is established to analyze the matrix porosity and apparent permeability evolutions during the process of gas migration. Using a thermodynamic formulation that treats the fluid–solid interface as an independent phase, we derive a rate form for matrix porosity and obtain the new dissipation function that contains three parts: dissipations from solid deformation, gas adsorption, and fluid flow. For gas adsorption, we justify the rationality of the adopted model; for fluid flow, the updated porosity can be substituted into sophisticated apparent permeability models for full‐scale analysis; and for solid deformation, a recently developed constitutive model appropriate for rocks exhibiting transverse isotropy in both the elastic and plastic responses is adopted in this work. Through the novel stress‐point simulation incorporating two effective stress measures and adsorption strain, new patterns of apparent permeability are obtained, which fit the experimental data quite well and cannot be reproduced from the assumption of isotropic poroelasticity. The advantages of our poromechanical model include thermodynamic consistency and the ability to employ finite‐element‐based formulation. Finally, an initial‐boundary value problem of gas production considering anisotropic plasticity is conducted, and the effects of the bedding plane and different adsorption models are highlighted.

Multi-objective topology optimization of heat transfer surface using level-set method and adaptive mesh refinement in OpenFOAM

The present study proposes a new efficient and robust algorithm for multi-objectives topology optimization of heat transfer surfaces to achieve heat transfer enhancement with a less pressure drop penalty based on a continuous adjoint approach. It is achieved with a customized OpenFOAM solver, which is based on a volume penalization method for solving a steady and laminar flow around iso-thermal solid objects with arbitrary geometries. The fluid-solid interface is captured by a level-set function combined with a newly proposed robust reinitialization scheme ensuring that the interface diffusion is always kept within a single local grid spacing. Adaptive mesh refinement is applied in near-wall regions automatically detected by the level-set function to keep high resolution locally, thereby reduces the overall computational cost for the forward and adjoint analyses. The developed solver is first validated in a drag reduction problem of a flow around a two-dimensional cylinder at the Reynolds numbers of 10 and 40 by comparing reference data. Then, the proposed scheme is extended to heat transfer problems in a two-dimensional flow at the Prandtl number of 0.7 and 6.9. Finally, three-dimensional topology optimization for multi-objective problems is considered for cost functionals with different weights for the total drag and heat transfer. Among various solutions obtained on the Pareto front, 4.0% of heat transfer enhancement with 12.6% drag reduction is achieved at the Reynolds number of 10 and the Prandtl number of 6.9. Moreover, the optimization of a staggered pin-fin array demonstrates that the optimal shapes and arrangement of the fins strongly depend on the number of rows from the inlet. Specifically, the pin-fins in the first and third rows extended in the upstream direction further enhance heat transfer, while the fins in the second row vanish to reduce pressure loss.

A model-based simulation framework for coupled acoustics, elastodynamics, and damage with application to nano-pulse lithotripsy

We develop a model for solid objects surrounded by a fluid that accounts for the possibility of acoustic pressures giving rise to damage on the surface of the solid. The propagation of an acoustic pressure in the fluid domain is modeled by the acoustic wave equation. On the other hand, the response of the solid is described by linear elastodynamics coupled with a gradient damage model, one that is based on a cohesive-type phase-field description of fracture. The interaction between the acoustic pressure and the deformation and damage of the solid are represented by transmission conditions at the fluid–solid interface. The resulting governing equations are discretized using a finite-element/finite-difference method that pays particular attention to the spatial and temporal scales that need to be resolved. Results from model-based simulations are provided for a benchmark problem as well as for recent experiments in nano-pulse lithotripsy. A parametric study is performed to illustrate how damage develops in response to the driving force (magnitude and location of the acoustic source) as a function of the fracture resistance of the solid. The results are shown to be qualitatively consistent with experimental observations for the location and size of the damage fields on the solid surface. A study of limiting cases also suggests that both the threshold for damage and the critical fracture energy are important to consider in order to capture the transition from damage initiation to complete localization. A low-cycle fatigue model is proposed that degrades the fracture resistance of the solid as a function of accumulated tensile strain energy, and it is shown to be capable of capturing damage localization in simulations of multi-pulse nano-pulse lithotripsy.

Flow field analysis and particle erosion of tunnel‐slope systems under coupling between runoff and fast (slow) seepage

AbstractThe presence of particles on the surface of a tunnel slope renders it susceptible to erosion by water flow, which is a major cause of soil and water loss. In this study, a nonlinear mathematical model and a mechanical equilibrium model are developed to investigate the distribution of flow fields and particle motion characteristics of tunnel slopes, respectively. The mathematical model of flow fields comprises three parts: a runoff region, a highly permeable soil layer, and a weakly permeable soil layer. The Navier‒Stokes equation controls fluid motion in the runoff region, while the Brinkman‐extended Darcy equation governs fast and slow seepage in the highly and weakly permeable soil layers, respectively. Analytical solutions are derived for the velocity profile and shear stress expression of the model flow field under the boundary condition of continuous transition of velocity and stress at the fluid‒solid interface. The shear stress distribution shows that the shear stress at the tunnel‐slope surface is the largest, followed by the shear stress of the soil interface, indicating that particles in these two locations are most vulnerable to erosion. A mechanical equilibrium model of sliding and rolling of single particles is established at the fluid‒solid interface, and the safety factor of particle motion (sliding and rolling) is derived. Sensitivity analysis shows that by increasing the runoff depth, slope angle, and soil permeability, the erosion of soil particles will be aggravated on the tunnel‐slope surface, but by increasing the particle diameter, particle‐specific gravity, and particle stacking angle, the erosion resistance ability of the tunnel‐slope surface particles will be enhanced. This study can serve as a reference for the analysis of surface soil and water loss in tunnel‐slope systems.

New trends in modelling of breakthrough curves to remove pollutants using adsorption on advanced monoliths geometries

This work proposes a rigorous mathematical model capable of reproducing the adsorption process in dynamic regime on advanced monoliths geometries. For this, four bed geometries with axisymmetric distribution of channels and similar solid mass were proposed. In each geometry a different distribution of channels was suggested, maintaining constant the bed dimensions of 15 cm high and 5 cm radius. The mathematical modeling includes mass and momentum transfer phenomena, and it was solved with the COMSOL Multiphysics software using mass transfer parameters published in the literature. The overall performance of the column was evaluated in terms of breakthrough (CA/CA0 = 0.1) and saturation times (CA/CA0 = 0.9). The mass and velocity distributions obtained from the proposed model show good physical consistency with what is expected in real systems. In addition, the model proved to be easy to solve given the short convergence times required (2–4 h). Modifications were made to the bed geometry to achieve a better use of the adsorbent material which reached up to 80%. The proposed bed geometries allow obtaining different mixing distributions, in such a way that inside the bed a thinning of the boundary layer is caused, thus reducing diffusive effects at the adsorbent solid-fluid interface, given dissipation rates of about 323 × 10−11 m2/s3. The bed geometry composed of intersecting rings deployed the best performance in terms of usage of the material adsorbent, and acceptable hydrodynamical behavior inside the channels (maximum fluid velocity = 35.4 × 10−5 m/s and drop pressure = 0.19 Pa). Based on these results, it was found that it is possible to reduce diffusional effects and delimit the mass transfer zone inside the monoliths, thus increasing the efficiency of adsorbent fixed beds.