Effect of Nasal Inhalation on Drug Particle Deposition and Size Distribution in the Upper Airway: With Soft Mist Inhalers

Delivery of drugs to the lungs is commonly achieved using nasal and/or oral breathing-assisted techniques. The route of inhalation can substantially change the fate of inhaled droplets. The Respimat® Soft Mist™ Inhaler (SMI) is a commercially available efficient inhaler with 40-60% effectiveness. In the present study, we used computational fluid dynamics (CFD) with a custom setup to investigate the effect of a combined oral/nasal inhalation route on the SMI's regional droplet deposition, size distribution, and flow field. Our setup used a modified induction port (MIP) to mimic nasal inhalation inside the human respiratory tract. Six different oral/nasal flow rate ratios inside the MIP were applied (total flow rate of 30 l/min). An overall good agreement was achieved between simulation outcomes and in vitro results. Our results confirmed that the combined inhalation route affects the flow field, altering the MIP's droplet deposition and size distribution. The lowest depositional loss, mainly in the mouth area, was observed at oral/nasal flow rate ratios of O/N = 1 and O/N = 2 with 3% and 7.7% values, respectively. Droplets with a 2-5 µm diameter range showed the highest droplet mass inside the MIP at all combined flow rates. We observed less intense vortexes followed by a lower level of turbulent kinetic energy at the oral/nasal ratio of 1. Increasing the relative humidity (RH) at oral/nasal flow rate ratios of 0.07, 1, and 14 led to an increase in droplet deposition at the outlet of the MIP.

Investigation of direct contact membrane distillation (DCMD) performance using CFD and machine learning approaches

Direct Contact Membrane Distillation (DCMD) is emerging as an effective method for water desalination, known for its efficiency and adaptability. This study delves into the performance of DCMD by integrating two powerful analytical tools: Computational Fluid Dynamics (CFD) and Artificial Neural Networks (ANN). The research thoroughly examines the impact of various factors, such as inlet temperatures, velocities, channel heights, salt concentration, and membrane characteristics, on the process's efficiency, specifically calculating the water vapor flux. A rigorous validation of the CFD model aligns well with established studies, ensuring reliability. Subsequently, over 1000 data points reflecting variations in input factors are utilized to train and validate the ANN. The training phase demonstrated high accuracy, with near-zero mean squared errors and R2 values close to one, indicating a strong predictive capability. Further analysis post-ANN training shed light on key relationships: higher membrane porosity boosts water vapor flux, whereas thicker membranes reduce it. Additionally, it was detailed how salt concentration, channel dimensions, inlet temperatures, and velocities significantly influence the distillation process. Finally, a mathematical model was proposed for water vapor flux as a function of key input factors. The results highlighted that salt mole fraction and hot water inlet temperature have the most effect on the water vapor flux. This comprehensive investigation contributes to the understanding of DCMD and emphasizes the potential of combining CFD and ANN for optimizing and innovating water desalination technology.

Influence of building spatial patterns on wind environment and air pollution dispersion inside an industrial park based on CFD simulation

The "spatial pattern-wind environment-air pollution" within building clusters is closely interconnected, where different spatial pattern parameters may have varying degrees of impact on the wind environment and pollutant dispersion. Due to the complex spatial structure within industrial parks, this complexity may lead to the accumulation and retention of air pollutants within the parks. Therefore, to alleviate the air pollution situation in industrial parks in China and achieve the circular transformation and construction of parks, this study takes Hefei Circular Economy Demonstration Park as the research object. The microscale Fluent model in computational fluid dynamics (CFD) is used to finely simulate the wind flow field and the diffusion process of pollutants within the park. The study analyzes the triad relationship and influence mechanism of "spatial pattern-wind environment-air pollution" within the park and studies the influence of different spatial pattern parameters on the migration and diffusion of pollutants. The results show a significant negative correlation between the content of pollutants and wind speed inside the industrial park. The better the wind conditions, the higher the air quality. The spatial morphology parameters of the building complex are the main influences on the condition of its internal wind environment. Building coverage ratio and degree of enclosure have a significant negative correlation with wind conditions. Maintaining them near 0.23 and 0.37, respectively, is favorable to the quality of the surrounding environment. Moreover, the average height of the building is positively correlated with the wind environment condition. The rate of transport and dissipation of pollutants gradually increases as the average building height reaches 16 m. Therefore, a reasonable building planning strategy and arrangement layout can effectively improve the wind environment condition inside the park, thus alleviating the pollutant retention situation. The obtained results serve as a theoretical foundation for optimizing morphological structure design within urban industrial parks.

Enhancing Inhalation Drug Delivery: A Comparative Study and Design Optimization of a Novel Valved Holding Chamber

This paper presents an innovative approach to the design optimization of valved holding chambers (VHCs), crucial devices for aerosol drug delivery. We present the design of an optimal cylindrical VHC body and introduce a novel valve based on particle impaction theory. The research combines computational simulations and physical experiments to assess the performance of various VHCs, with a special focus on the deposition patterns of medication particles within these devices. The methodology incorporates both experimental and simulation approaches to validate the reliability of the simulation. Emphasis is placed on the deposition patterns observed on the VHC walls and the classification of fine and large particles for salbutamol sulfate particles. The study reveals the superior efficacy of our valve design in separating particles compared to commercially available VHCs. In standard conditions, our valve design allows over 95% of particles under 7 μm to pass through while effectively filtering those larger than 8 μm. The optimized body design accomplishes a 60% particle mass flow fraction at the outlet and an average particle size reduction of 58.5%. When compared numerically in terms of size reduction, the optimal design outperforms the two commercially available VHCs selected. This study provides valuable insights into the optimization of VHC design, offering significant potential for improved aerosol drug delivery. Our findings demonstrate a new path forward for future studies, aiming to further optimize the design and performance of VHCs for enhanced pulmonary drug delivery.

Determination of Significant Three-Dimensional Hemodynamic Features for Postembolization Recanalization in Cerebral Aneurysms Through Explainable Artificial Intelligence

BACKGROUND

Recanalization poses challenges after coil embolization in cerebral aneurysms. Establishing predictive models for postembolization recanalization is important for clinical decision making. However, conventional statistical and machine learning (ML) models may overlook critical parameters during the initial selection process.

METHODS

In this study, we automated the identification of significant hemodynamic parameters using a PointNet-based deep neural network (DNN), leveraging their three-dimensional spatial features. Further feature analysis was conducted using saliency mapping, an explainable artificial intelligence (XAI) technique. The study encompassed the analysis of velocity, pressure, and wall shear stress in both precoiling and postcoiling models derived from computational fluid dynamics simulations for 58 aneurysms.

RESULTS

Velocity was identified as the most pivotal parameter, supported by the lowest P value from statistical analysis and the highest area under the receiver operating characteristic curves/precision-recall curves values from the DNN model. Moreover, visual XAI analysis showed that robust injection flow zones, with notable impingement points in precoiling models, as well as pronounced interplay between flow dynamics and the coiling plane, were important three-dimensional features in identifying the recanalized aneurysms.

CONCLUSIONS

The combination of DNN and XAI was found to be an accurate and explainable approach not only at predicting postembolization recanalization but also at discovering unknown features in the future.

Scale resolving simulations of the effect of glottis motion and the laryngeal jet on flow dynamics during respiration

BACKGROUND AND OBJECTIVE

The movement of the respiratory walls has a significant impact on airflow through the respiratory tract. The majority of computational fluid dynamics (CFD) studies assume a static geometry which may not provide a realistic flow field. Furthermore, many studies use Reynolds Averaged Navier-Stokes (RANS) turbulence models that do not resolve turbulence structure. Combining the application of advanced scale-resolving turbulence models with moving respiratory walls using CFD will provide detailed insights into respiratory flow structures.

METHODS

This study simulated a complete breathing cycle involving inhalation and exhalation in a nasal cavity to trachea geometry that incorporated moving glottis walls. A second breathing cycle was simulated with static glottis walls for comparison. A recently developed hybrid RANS-LES turbulence model, the Stress-Blended Eddy Simulation (SBES), was incorporated to resolve turbulent flow structures in fine detail for both transient simulations. Transient results were compared with steady-state RANS simulations for the same respiratory geometry.

RESULTS

Glottis motion caused substantial effects on flow structure through the complete breathing cycle. Significant flow structure and velocity variations were observed due to glottal motion, primarily in the larynx and trachea. Resolved turbulence structures using SBES showed an intense mixing section in the glottis region during inhalation and in the nasopharynx during expiration, which was not present in the RANS simulations.

CONCLUSION

Transient simulations of a realistic breathing cycle uncovered flow structures absent in simulations with a constant flow rate. Furthermore, the incorporation of glottis motion impacted airflow characteristics that suggest rigid respiratory walls do not accurately describe respiratory flow. Future research in respiratory airflow should be conducted using transient scale-resolving models in conjunction with moving respiratory walls to capture flow structures in detail.

Computational Analysis of the Effects of Fiber Deformation on the Microstructure and Permeability of Blood Oxygenator Bundles

Mechanical loads on the polymeric fibers of oxygenating bundles are commonly present due to bundle press-fitting during device assembly and blood pressure load. However, computational fluid dynamics (CFD) simulations for fiber bundle optimization neglect possible changes in microstructure due to such deformations. The aim of this study is to investigate the impact of fiber deformability on bundle microstructure and fluid dynamics mainly in terms of permeability. Fibers from commercial mats typically used for blood oxygenators were mechanically tested and based on these experimental data, a material model was developed to simulate the structural deformations the fibers undergo under press-fitting and blood pressure loads. Then, CFD simulations were performed on deformed bundle repetitive units to investigate permeability under varying loading conditions. The effects of different bundle geometric parameters on the variation of bundle permeability due to press-fitting were evaluated. Bundle press-fitting results in significant changes in microstructure that are reflected in a bundle permeability more than halved for a 15% press-fitting. This impact on permeability is present in all the simulated fiber bundles and becomes more pronounced as the pitch between fibers and thus bundle porosity decreases. Instead, the analyses on pressurized bundle show only small deformations caused by pressure load, with permeability changes below 1%. While blood pressure effects could be neglected, bundle press-fitting turns out to have a significant impact on bundle microstructure and permeability. Neglecting such microstructure variations during CFD simulations could also lead to incorrect assessment of the local fluid dynamics within the bundle.

Optimization of up-flow anaerobic sludge blanket second-layer influent distribution structure

ABSTRACTA parametric model of the second-layer influent distributor is proposed to increase the reaction efficiency of the up-flow anaerobic sludge blanket. The impacts on the flow efficiency of the main parameters, including the length of distribution pipe, that is, the cylinder radius r1, the eccentricity of the nozzle e, and the number of nozzles n1, are investigated. The optimal parameter combination of the second-layer influent distributoris obtained by single factor analysis and orthogonal analysis. Then the two-layer distributor combination model, including the bottom influent distributor and the second-layer influent distributor, is established and the simulations are conducted to study the effect of the two-layer distributor on flow inside the reactor. The simulation results show that the proposed two-layer distributor can achieve a higher mixture efficiency than the single distributor.

Flow patterns in ascending aortic aneurysms: Determining the role of hypertension using phase contrast magnetic resonance and computational fluid dynamics

Thoracic aortic aneurysm (TAA) is a local dilation of the thoracic aorta. Although universally used, aneurysm diameter alone is a poor predictor of major complications such as rupture. There is a need for better biomarkers for risk assessment that also reflect the aberrant flow patterns found in TAAs. Furthermore, hypertension is often present in TAA patients and may play a role in progression of aneurysm. The exact relation between TAAs and hypertension is poorly understood. This study aims to create a numerical model of hypertension in the aorta by using computational fluid dynamics. First, a normotensive state was simulated in which flow and resistance were kept unaltered. Second, a hypertensive state was modeled in which blood inflow was increased by 30%. Third, a hypertensive state was modeled in which the proximal and peripheral resistances and capacitance parameters from the three-element Windkessel boundary condition were adjusted to mimic an increase in resistance of the rest of the cardiovascular system. One patient with degenerative TAA and one healthy control were successfully simulated at hypertensive states and were extensively analyzed. Furthermore, three additional TAA patients and controls were simulated to validate our method. Hemodynamic variables such as wall shear stress, oscillatory shear index, endothelial cell activation potential (ECAP), vorticity and helicity were studied to gain more insight on the effects of hypertension on flow patterns in TAAs. By comparing a TAA patient and a control at normotensive state at peak-systole, helicity and vorticity were found to be lower in the TAA patient throughout the entire domain. No major changes in flow and flow derived quantities were observed for the TAA patient and control when resistance was increased. When flow rate was increased, regions with high ECAP values were found to reduce in TAA patients in the aneurysm region which could reduce the risk of thrombogenesis. Thus, it may be important to assess cardiac output in patients with TAA.

The Impact of Nozzle Opening Thickness on Flow Characteristics and Primary Electron Beam Scattering in an Environmental Scanning Electron Microscope

This paper describes the methodology of combining experimental measurements with mathematical-physics analyses in the investigation of flow in the aperture and nozzle. The aperture and nozzle separate the differentially pumped chamber from the specimen chamber in an environmental scanning electron microscope (ESEM). Experimental measurements are provided by temperature and pressure sensors that meet the demanding conditions of cryogenic temperature zones and low pressures. This aperture maintains the required pressure difference between the chambers. Since it separates the large pressure gradient, critical flow occurs on it and supersonic gas flow with the characteristic properties of critical flow in the state variables occurs behind it. As a primary electron beam passes through the differential pumped chamber and the given aperture, the aperture is equipped with a nozzle. The shape of the nozzle strongly influences the character of the supersonic flow. The course of state variables is also strongly influenced by this shape; thus, it affects the number of collisions the primary beam's electrons have with gas molecules, and so the resulting image. This paper describes experimental measurements made using sensors under laboratory conditions in a specially created experimental chamber. Then, validation using mathematical-physical analysis in the Ansys Fluent system is described.

Numerical study of drop dynamics for inkjet based 3D printing of pharmaceutical tablets

Interest in 3D printing has been growing rapidly especially in pharmaceutical industry due to its multiple advantages such as manufacturing versatility, personalization of medicine, scalability, and cost effectiveness. Inkjet based 3D printing gained special attention after FDA's approval of Spritam® manufactured by Aprecia pharmaceuticals in 2015. The precision and printing efficiency of 3D printing is strongly influenced by the dynamics of ink/binder jetting, which further depends on the ink's fluid properties. In this study, Computational Fluid Dynamics (CFD) has been utilized to study the drop formation process during inkjet-based 3D printing for piezoelectric and thermal printhead geometries using Volume of Fluid (VOF) method. To develop the CFD model commercial software ANSYS-Fluent was used. The developed CFD model was experimentally validated using drop watcher setup to record drop progression and drop velocity. During the study, water, Fujifilm model fluid, and Amitriptyline drug solutions were evaluated as the ink solutions. The drop properties such as drop volume, drop diameter, and drop velocity were examined in detail in response to change ink solution properties such as surface tension, viscosity, and density. A good agreement was observed between the experimental and simulation data for drop properties such as drop volume and drop velocity.

Numerical performance of CO2 accumulation and droplet dispersion from a cough inside a hospital lift under different ventilation strategies

The impact of mechanical ventilation on airborne diseases is not completely known. The recent pandemic of COVID-19 clearly showed that additional investigations are necessary. The use of computational tools is an advantage that needs to be included in the study of designing safe places. The current study focused on a hospital lift where two subjects were included: a healthy passenger and an infected one. The elevator was modelled with a fan placed on the middle of the ceiling and racks for supplying air at the bottom of the lateral wall. Three ventilation strategies were evaluated: a without ventilation case, an upwards-blowing exhausting fan case and a downwards-blowing fan case. Five seconds after the elevator journey began, the infected person coughed. For the risk assessment, the CO2 concentration, droplet removal performance and dispersion were examined and compared among the three cases. The results revealed some discrepancies in the selection of an optimal ventilation strategy. Depending on the evaluated parameter, downward-ventilation fan or no ventilation strategy could be the most appropriate approach.

Intrinsic dissolution rate modeling for the pharmacopoeia apparatus rotating disk compared to flow channel method

For a solid understanding of drug characteristics, in vitro measurement of the intrinsic dissolution rate is important. Hydrodynamics are often emphasized as the decisive parameter influencing the dissolution. In this study, experiments and computational fluid dynamic (CFD) simulations showed that the mixing behavior in the rotating disc apparatus causes an inhomogeneous flow field and a systematic error in the calculation of the intrinsic dissolution rate. This error is affected by both the experimental time and the velocity. Due to the rotational movement around the tablet center, commonly utilized in pharmacopeia methods, a broad variance is present with regard to the impact of fluid velocity on individual particles of the specimen surface. As this is significantly reduced in the case of uniform overflow, the flow channel is recommended for investigating the dissolution behavior. It is shown that rotating disc measurements can be compared with flow channel measurements after adjusting the measured data for the rotating disc based on a proposed, representative Reynolds number and a suggested apparatus-dependent correction factor. Additionally, modeling the apparatus-independent intrinsic dissolution rate for different temperatures in the rotating disc apparatus is possible using the adapted Levich's equation.

Using shock generator for the fuel mixing of the extruded single 4-lobe nozzle at supersonic combustion chamber

The importance of the fuel injection configuration on the propulsion efficiency of high-speed vehicles is apparent. In this article, the use of an annular extruded 4-lobe nozzle for the injection of fuel jet in a supersonic combustor of a scramjet engine in the existence of a shock generator is examined. The main aim of this study is to obtain the efficient jet arrangement for efficient fuel mixing inside the engine of hypersonic vehicles. A numerical approach is used to model the supersonic air stream and cross-jet flow with the SST turbulence model. The role of nozzle altitude and internal air jet on the fuel mixing of the hydrogen within the high-speed domain are disclosed. The importance of the horseshoe vortex and counter-rotating vortex on the fuel distribution is also presented. Our results show that the usage of a coaxial jet instead of an annular jet would increase fuel mixing by more than 40% in the combustion chamber.

Design, control, aerodynamic performances, and structural integrity investigations of compact ducted drone with co-axial propeller for high altitude surveillance

Compact multi-rotor unmanned aerial vehicles (UAVs) can be operated in many challenging environmental conditions. In case the UAV requires certain considerations in designing like lightweight, efficient propulsion system and others depending upon the application, the hybrid UAV comes into play when the usual UAV types cannot be sufficient to meet the requirements. The propulsion system for the UAV was selected to be coaxial rotors because it has a high thrust-to-weight ratio and to increase the efficiency of the propulsion system, a unique propeller was proposed to achieve higher thrust. The proposed propeller was uniquely designed by analyzing various airfoil sections under different Reynolds's number using X-Foil tool to obtain the optimum airfoil section for the propellers. Since the design with duct increases efficiency, the Hybrid UAV presented in this paper has the modified novel convergent-divergent (C-D)-based duct which is a simplified model of a conventional C-D duct. The yawing and rolling maneuverings of the UAV could be achieved by the thrust vectoring method so that the design is simpler from a structural and mechanical perspective. The use of UAVs has risen in recent years, especially compact UAVs, which can be applied for applications like surveillance, detection and inspection, and monitoring in a narrow region of space. The design of the UAV is modeled in CATIA, and its further performance enactment factors are picked from advanced computational simulations relayed bottom-up approach. The predominant computational fluid dynamics (CFD) and fluid structure interaction (FSI) investigations are imposed and optimized through Computational Analyses using Ansys Workbench 17.2, which includes analysis of structural behaviour of various alloys, CFRP and GFRP based composite materials. From the structural analysis Titanium alloy came out to be the best performing materials among the others by having lower total deformation and other parameters such as normal and equivalent stress. The dynamics control response is obtained using MATLAB Simulink. The validations are carried out on the propeller using a thrust stand for CFD and on the duct through a high-jet facility for structural outcomes to meet the expected outcome.

Modelling and simulation of fuel tank with increased capacity and improving its location for better stability of three wheeled vehicle

In the automotive industry, the fuel tank is an essential component of the vehicle. It is; designed as an integral part of the fuel system to transport the fuel and deliver it to the engine via the fuel filter and fuel pump. In Ethiopia, it is observed that due to the smaller fuel tank of three-wheeled vehicles, more frequent visits to fuel stations result in a lot of time being wasted waiting in long queues at fuel stations. In addition, the current location of the fuel tank causes some stability problems. Subsequently, the method, location, installation techniques and stability analysis of the three-wheeled vehicle are carried out. This study includes a new design of fuel tanks for three-wheeled vehicles in which the capacity is doubled and the rollover stability; of the vehicle is maintained. The SOLIDWORKS software and the CFD software ANSYS FLUENT for the sloshing analysis were used for the modeling of this work. The CFD results of the comparative sloshing analysis of partially filled (50% level) fuel tanks with and without baffles are performed; and conclusions are drawn. It is concluded that a baffle plate for the fuel tank offers the greatest advantage in reducing the effects of sloshing. The newly developed part is mounted on the underside of the passenger seat. A bottom-mounted fuel tank optimizes the vehicle's center of gravity, as the height of the vehicle's center of gravity is lowered by 30 mm. This increases the rollover safety of the vehicle. According to the analysis, the mileage of the filled existing and new tanks is 136 km and 272 km respectively. To refuel once, the driver waits an average of 5 h and works about two days for the existing tank and four days for the new tank.

Numerical modeling and CFD simulation of diffuser augmented dual vertical axis hydrokinetic Banki-Michell turbine

Hydrokinetic Banki turbines present an affordable, technically feasible, environmentally friendly technology. Their construction without requiring more expensive structures like diversion weirs, canals, forebay, and penstock, makes their initial investment much lower than commonly used horizontal Banki turbine of the same capacity. The possibility to install in the existing canals for Ultra Low Head applications is the additional motivating factor for this research. The system studied includes two Banki runners without internal shafts mounted vertically side by side surrounded by nozzle and diffuser structures. In the first scenario, Nozzle and then the Nozzle-diffuser augmented structures were separately studied to enhance the output of the runner for ultra-low head application, and the effects of each on the speed, pressure, and power output were analyzed. For the case of commonly used Banki, without nozzle and diffuser augmentation the speed for Ultra Low Head was minimum and determined to be 344 rpm, which is far below the recommended value of 800 rpm for safe operation at a flow rate of 1 m^3/s. In view of this, in the present study the enhanced speed on account of improvement was found to be 850 rpm and 1025 rpm for the design without and with diffuser assemblies respectively. Besides, the performance is seen to be improved by 7.6% with the diffuser as compared with the one without diffuser assembly. Detailed simulation results are presented and discussed: 3D ANSYS-FLUENT optimization result provided optimum number of blades for each runner to be 19 and with the optimum throat width in both cases as 202 mm. On account of the lack of any results reported so far for this innovative geometry, validation of the simulated results was carried out with reported results for the dual horizontal axis Banki turbines with good agreement.

Parametric study of an empty diffuser geometric parameters and shape for a wind turbine using CFD analysis

This research examines the optimal 3D geometric parameters and shape of empty diffusers to enhance the mass flow rate of the HAWT rotor plane using a detailed parametric study. Previous works have investigated the use of diffusers to augment turbine power output; however, different curvature surfaces and the effects of all associated angles have not been considered for a thorough evaluation. This work mainly focuses on analyzing the effect of opening angles (2° to 22°), inlet shroud angles (8° to 24°), flange height ratios, flange angles (0° and 15°), and shape of the diffuser as well as flanges on velocity, pressure at the diffuser entry, and through the diffuser section at a wind speed of 4.5 m/s. At an inlet-shroud angle of 24° and an opening angle of 8°, with a diffuser flange height-throat diameter ratio of 0.3, the system achieved an 82.9% increase in flowrate. The diffuser with an inlet shroud-side lower stepped flange showed an optimum velocity of 9.12 m/s (maximum) and 8.2 m/s (average), resulting in a 102.66% and 82.2% increase in velocity, respectively. The percentage increase in velocity of the present study is 92.61%, compared with the previous maximum increase in rate of 53.8%, and then an increase in velocity of 38.81% was obtained. The optimum speed occurred at 0.175 m from the inlet section of the diffuser, indicating where the turbine should best be located. The CFD results from this work were validated with experimental data from the literature, showing a good agreement between the two. Integrated diffuser-turbine system simulation and experimental work with field tests are recommended as a way forward.

Design of experiment approach to identify the dominant geometrical feature of left coronary artery influencing atherosclerosis

Background and Objective. Coronary artery geometry heavily influences local hemodynamics, potentially leading to atherosclerosis. Consequently, the unique geometrical configuration of an individual by birth can be associated with future risk of atherosclerosis. Although current researches focus on exploring the relationship between local hemodynamics and coronary artery geometry, this study aims to identify the order of influence of the geometrical features through systematic experiments, which can reveal the dominant geometrical feature for future risk assessment.Methods. According to Taguchi's method of design of experiment (DoE), the left main stem (LMS) length (lLMS), curvature (kLMS), diameter (dLMS) and the bifurcation angle between left anterior descending (LAD) and left circumflex (LCx) artery (αLAD-LCx) of two reconstructed patient-specific left coronary arteries (LCA) were varied in three levels to create L9 orthogonal array. Computational fluid dynamic (CFD) simulations with physiological boundary conditions were performed on the resulting eighteen LCA models. Average helicity intensity (h2) and relative atheroprone area (RAA) of near-wall hemodynamic descriptors were analyzed.Results. The proximal LAD (LADproximal) was identified to be the most atheroprone region of the left coronary artery due to higherh2,large RAA of time averaged wall shear stress (TAWSS < 0.4 Pa), oscillatory shear index (OSI ∼ 0.5) and relative residence time (RRT > 4.17 Pa-1). In both patient-specific cases, based onh2and TAWSS,dlmsis the dominant geometric parameter while based on OSI and RRT,αLAD-LCxis the dominant one influencing hemodynamic condition in proximal LAD (p< 0.05). Based on RRT, the rank of the geometrical factors is:αLAD-LCx>dLMS>lLMS>kLMS, indicating thatαLAD-LCxis the most dominant geometrical factor affecting hemodynamics at proximal LAD which may influence atherosclerosis.Conclusion. The proposed identification of the rank of geometrical features of LCA and the dominant feature may assist clinicians in predicting the possibility of atherosclerosis, of an individual, long before it will occur. This study can further be translated to be used to rank the influence of several arterial geometrical features at different arterial locations to explore detailed relationships between the arterial geometrical features and local hemodynamics.

Simulation on two-phase refrigerant compression in the cylinder of rotary compressors using CFD method

The two-phase compression process in the rotary compressor often occurs, such as defrosting and startup processes, which has a significant impact on the performance and reliability of air conditioning systems. In this paper, the CFD simulations predicting the two-phase refrigerant compression process in the compressor cylinder are conducted using the commercial software ANSYS Fluent. The dynamic mesh for the fluid domain and phase change model for the refrigerant are considered in the simulation. Effects of initial liquid volume fraction, refrigerant type and compressor type on the two-phase compression characteristics using R290 as refrigerant are carried out. Variations of the pressure, temperature, gas fraction distribution and evaporation rate in the cylinder are discussed. The results show that most liquid accumulates near the leakage gap and the bottom of the compression chamber during the two-phase compression process. The peak pressure during the two-phase compression decreases with the increase of the liquid volume fraction. The evaporation rate of R32 in the cylinder is much higher than that of R290. The maximum pressure of the reciprocating compressor is 2.26 times higher than that of the rotary compressor.

Tailoring left ventricular assist device cannula implantation using coupled multi-scale multi-objective optimization

BACKGROUND

The frequent occurrence of thromboembolic cerebral events continues to limit the widespread implementation of Ventricular Assist Devices (VAD) despite continued advancements in VAD design and anti-coagulation treatments. Recent studies point to the optimal positioning of the outflow graft (OG) as a potential mitigator of post implantation thromboembolism.

OBJECTIVE

This study aims to examine the tailoring of the OG implantation orientation with the goal of minimizing the number of thrombi reaching the cerebral vessels by means of a formal shape optimization scheme incorporated into a multi-scale hemodynamics analysis.

METHODS

A 3-D patient-specific computational fluid dynamics model is loosely coupled in a two-way manner to a 0-D lumped parameter model of the peripheral circulation. A Lagrangian particle-tracking scheme models and tracks thrombi as non-interacting solid spheres. The loose coupling between CFD and LPM is integrated into a geometric shape optimization scheme which aims to optimize an objective function that targets a drop in cerebral embolization, and an overall reduction in particle residence times.

RESULTS

The results elucidate the importance of OG anastomosis orientation and placement particularly in the case that studied particle release from the OG, as a fivefold decrease in cerebral embolization was observed between the optimal and non-optimal implantations. Another case considered particle release from the ventricle and aortic root walls, in which optimal implantation was achieved with a shallow insertion angle. Particle release from all three origins was investigated in the third case, demonstrating that the optimal configurations were generally characterized by VAD flow directed along the central lumen of the aortic arch. Because optimal configurations depended on the anatomic origin of the thrombus, it is important to determine, in clinical studies, the most likely sites of thrombus formation in VAD patients.

Flow characterization of Maquet and Bio-Medicus multi-stage drainage cannulae during venoarterial extracorporeal membrane oxygenation

BACKGROUND

Drainage cannulae extract blood from a patient during venoarterial extracorporeal membrane oxygenation (VA ECMO), a treatment that temporarily supports patients undergoing severe heart and/or lung dysfunction. Currently, the two most commonly used multi-stage drainage cannulae are manufactured by Maquet and Bio-Medicus, but their designs vary in many aspects which impacts the generated flow dynamics. Therefore, this study aimed to use computational fluid dynamics (CFD) to explore the flow characteristics of the aforementioned cannulae and their impact on complications such as thrombosis.

METHODS

The Maquet and Bio-Medicus cannulae were 3D modelled within a patient-specific geometry of the venous vasculature taken from a computed tomography scan of a patient undergoing VA ECMO. A drainage flow rate of 4 L/min was assigned to each cannula. Lastly, a stress blended eddy simulation turbulence model was employed to resolve bulk flow turbulence.

RESULTS

The proximal row of side holes in both cannulae generated high intensity counter-rotating vortices, thus generating supraphysiological shear. These proximal rows were also responsible for the majority of flow extraction in both cannulae (>1.6 L/min). Despite identical simulation settings, each cannulae had differing impacts on global flow dynamics. For instance, the Bio-Medicus model produced a total stagnant blood volume of 25.6 ml, compared to 17.8 ml the Maquet cannula, thereby increasing the risk of thrombosis.

CONCLUSIONS

Overall, our results demonstrate that differences in design clearly impact flow dynamics and risk of complications. Therefore, further work in optimizing cannula design may be beneficial to prevent harmful flow characteristics.

Operator learning for urban water clarification hydrodynamics and particulate matter transport with physics-informed neural networks

Computational fluid dynamics (CFD) can be a powerful tool for higher-fidelity water infrastructure planning and design. Despite decades of development and demonstration over a wide range of water systems such as clarification basins, activated sludge processes, ozone contactors, etc., CFD remains primarily used in academic research, with limited application in civil and environmental engineering practice. This limitation is contributed by its higher computational cost and demand for specialized user skills. This, however, need not be the case, if a robust and efficient surrogate model can be developed from CFD simulations and independently deployed for engineering purposes. Leveraging the emerging scientific machine learning (ML) techniques of physics-informed ML and operator learning, this study develops a composite neural network (CPNN) for learning the flow hydrodynamics and particulate matter (PM) transport and fate in clarification systems. The CPNN consists of a deep operator network (DeepONet) as an encoder and a physics-informed neural network (PINN) as a decoder. In contrast to common "black box" and lumped ML approaches, the developed CPNN directly incorporates physics principles into its architecture. Furthermore, the CPNN is designed for process-resolved and operator learning, enabling it to predict spatial hydrodynamics and PM concentration distribution (i.e., contours) across different basin geometrics and loading conditions. Compared to CFD simulation, the developed CPNN model has significantly higher computational efficiency (∼ milliseconds) while demonstrating robust predictive capability. For predicting basin hydrodynamics across 10,000 test cases, the trained CPNN model achieves an R2 above 0.8 for 66.4% of cases and an R2 above 0.4 for 89.2% of cases. A similar performance is also demonstrated by the CPNN in predicting basin PM concentration. Further investigation reveals that basin geometrics that trigger bi-modal flow solutions can be particularly challenging for ML. Additionally, this study visualizes the dependency of basin hydrodynamics and PM concentration on basin geometrics and loading conditions, providing valuable insights for optimizing basin configuration. Lastly, the potentials and benefits of web-based applications, e.g., DeepXtorm, as a user-friendly interface for the developed CPNN model is discussed. This study represents the initial step toward achieving real-time higher-fidelity water infrastructure planning, design, optimization, and regulation.

Engineering innovations in medicine and biology: Revolutionizing patient care through mechanical solutions

The overlap between mechanical engineering and medicine is expanding more and more over the years. Engineers are now using their expertise to design and create functional biomaterials and are continually collaborating with physicians to improve patient health. In this review, we explore the state of scientific knowledge in the areas of biomaterials, biomechanics, nanomechanics, and computational fluid dynamics (CFD) in relation to the pharmaceutical and medical industry. Focusing on current research and breakthroughs, we provide an overview of how these fields are being used to create new technologies for medical treatments of human patients. Barriers and constraints in these fields, as well as ways to overcome them, are also described in this review. Finally, the potential for future advances in biomaterials to fundamentally change the current approach to medicine and biology is also discussed.

A general-purpose tool for modeling multifunctional thin porous media (POREnet): From pore network to effective property tensors

POREnet, a novel approach to model effective properties of thin porous media, TPM, is presented. The methodology allows the extraction of local effective property tensors by volume averaging from discrete pore networks, PNs, built on the tessellated continuum space of a TPM. The gradient theorem is used to describe 3D transport in bulk tessellated space, providing an appropriate metric to normalize network fluxes. Implemented effective transport properties include diffusivity, permeability, solid-phase conductivity, and entry capillary pressure and contact angle under two-phase conditions, considering multi-component materials with several solid phases and local contact resistances. Calculated property tensors can be saved on 3D image stacks, where interfacial and sub-CV scale features can be added before exporting data to CFD meshes for simulation. Overall, POREnet provides a general-purpose, versatile methodology for modeling TPM in an ample range of conditions within a single CFD framework. Among other advantages, coupling of PN and continuum models at TPM-channel interfaces is simplified, interfacial contact resistances can be included using robin boundary conditions, and transient multiphysics simulations can be implemented more easily using CFD. The code is tested against a miscellaneousness of examples extracted from electrochemical applications.

The Effect of Mechanical Circulatory Support on Blood Flow in the Ascending Aorta: A Combined Experimental and Computational Study

Percutaneous mechanical circulatory support (MCS) devices are designed for short-term treatment in cases of acute decompensated heart failure as a bridge to transplant or recovery. Some of the known complications of MCS treatments are related to their hemodynamics in the aorta. The current study investigates the effect of MCS on the aortic flow. The study uses combined experimental and numerical methods to delineate complex flow structures. Particle image velocimetry (PIV) is used to capture the vortical and turbulent flow characteristics in a glass model of the human aorta. Computational fluid dynamics (CFD) analyses are used to complete the 3D flow in the aorta. Three specific MCS configurations are examined: a suction pump with a counterclockwise (CCW) rotating impeller, a suction pump with a clockwise (CW) rotating impeller, and a discharge pump with a straight jet. These models were examined under varying flow rates (1-2.5 L/min). The results show that the pump configuration strongly influences the flow in the thoracic aorta. The rotating impeller of the suction pump induces a dominant swirling flow in the aorta. The swirling flow distributes the incoming jet and reduces the turbulent intensity near the aortic valve and in the aorta. In addition, at high flow rates, the local vortices formed near the pump are washed downstream toward the aortic arch. Specifically, an MCS device with a CCW rotating impeller induces a non-physiological CCW helical flow in the descending aorta (which is opposite to the natural helical flow), while CW swirl combines better with the natural helical flow.

Study of cooling experiment and simulation for edible oil storage

This paper proposes a refrigerant cooling method using an inner tube in a storage tank to improve the cooling performance and thermal uniformity during the storing of edible oil. With a prototype of an oil tank in Central Grain Reserve of Zhenjiang, the experimental oil tank was built in a scale of 50:1. Both natural and manual cooling experiments were carried out for the experimental tank. The manual cooling process involved two supplying modes for the refrigerant tube (top and bottom) and four different refrigerant temperatures (10 ℃, 12 ℃, 14 ℃, 16 ℃). The experimental results show that, compared with natural cooling, manual cooling can effectively reduce the temperature difference and thermal stratification between upper and lower layers. The temperature difference is 6.79 ℃, 1.93 ℃, and 3.67 ℃ for the natural cooling, manual top supplying, and manual bottom supplying mode, respectively. Furthermore, for the two manual modes, the cooling efficiency of bottom supplying is 21.4% higher than that of the top supplying, and the average oil temperature drops by 0.8-1 ℃. Based on experimental results, different working conditions (20, 40, and 60 ml/s) were simulated to determine the optimal flow rate for bottom supplying mode. The simulation results indicate that the low flow rate (20 ml/s) corresponds to the best thermal uniformity, and the maximum temperature has no obvious change under different flow rate conditions. Therefore, it is not necessary to increase the flow rate to improve cooling efficiency considering the rising energy consumption.

Agglomeration of particulate matter in chimneys using acoustic flow

The emission of micrometer-sized particulate matter on an industrial scale is causing increasing environmental concern about air pollution. Numerous industries and research communities need help to reduce micrometer-sized pollutants in the atmosphere. The current research investigates the acoustic agglomeration of particulate matter through a combination of experimental and numerical methods. Acoustic agglomeration is a process that involves using acoustic waves to influence the movement of particles in the air. Acoustic agglomeration operates by facilitating particle collision and simplifying the formation of agglomerates that are later removed through filtration. This article is focused on research on acoustic pre-processing with the aim of reducing atmospheric pollution caused by toxic combustion products. The capture of fine silica particles with diameters ranging from 0.3 to 10 μm, emitted through the chimneys of industrial enterprises, can be considered a significant technological innovation. The experimental part of the current research is conducted using a newly developed experimental bench. The assembly comprises the following key components: a wind tunnel, a particle dosing device, the agglomeration camera, and a particle concentration measurement device on the edge of the bench. A loudspeaker was used to evaluate the effect of sound pressure in the frequency range of 500-3000 Hz. A comprehensive CFD study of the particles was conducted, which included analysis of the boundary layer, facilitating a better understanding of the behavior of the particles and its potential to agglomerate. An experimental study of particle agglomeration, using an acoustic field with a frequency range of 500-3000 Hz, demonstrated the effectiveness of particle agglomeration of different diameters. The efficacy of particle agglomeration is up to 80 % when the sound pressure values were 129-135 dB; the highest efficiency was found at excitation frequencies of 1500 and 3000 Hz, respectively.

Mathematical modelling of a single tethered aerostat using longitudinal stability derivatives

Lighter-than-air (LTA) aerial vehicles such as airships and aerostats can be found in various strategic and commercial applications, primarily due to their capability to hover and stealth. The mathematical model of these vehicles helps in understanding their complex dynamics and designing and developing proper stabilisation systems for them. Stability derivatives have been used for developing mathematical models for heavier-than-air aerial vehicles since their introduction. This paper presents a methodology to develop a mathematical model of an aerostat based on stability derivatives. One of the major contributions of this study is the estimation of aerostat's added mass terms expressed as longitudinal stability derivatives due to acceleration of the longitudinal motion variables. A longitudinally decoupled linear mathematical model of a single-tethered aerostat using stability derivatives is investigated in this study. A computational fluid dynamics (CFD)-based analysis of the 3D model of the vehicle is used to obtain the stability derivatives. The methodology presented considers the aerostat and tether models separately before coupling them to create the full model. The stability derivative analysis is carried out using ANSYS Fluent, and the coupled tethered aerostat model is investigated using MATLAB 2020. The negative pitch angle of the aerostat is caused by the selection of the pitching centre as the aerostat centre of volume instead of the tether confluence point. The tension force on the tether, which is proportional to the wind velocity, and aerostat velocity components are found to be stabilised within 200-400 s.

A Novel Miniature and Selective CMOS Gas Sensor for Gas Mixture Analysis-Part 4: The Effect of Humidity

This is the fourth part of a study presenting a miniature, combustion-type gas sensor (dubbed GMOS) based on a novel thermal sensor (dubbed TMOS). The TMOS is a micromachined CMOS-SOI transistor, which acts as the sensing element and is integrated with a catalytic reaction plate, where ignition of the gas takes place. The GMOS measures the temperature change due to a combustion exothermic reaction. The controlling parameters of the sensor are the ignition temperature applied to the catalytic layer and the increased temperature of the hotplate due to the released power of the combustion reaction. The solid-state device applies electrical parameters, which are related to the thermal parameters. The heating is applied by Joule heating with a resistor underneath the catalytic layer while the signal is monitored by the change in voltage of the TMOS sensor. Voltage, like temperature, is an intensive parameter, and one always measures changes in such parameters relative to a reference point. The reference point for both parameters (temperature and voltage) is the blind sensor, without any catalytic layer and hence where no reaction takes place. The present paper focuses on the study of the effect of humidity upon performance. In real life, the sensors are exposed to environmental parameters, where humidity plays a significant role. Humidity is high in storage rooms of fruits and vegetables, in refrigerators, in silos, in fields as well as in homes and cars. This study is significant and innovative since it extends our understanding of the performance of the GMOS, as well as pellistor sensors in general, in the presence of humidity. The three main challenges in simulating the performance are (i) how to define the operating temperature based on the input parameters of the heater voltage in the presence of humidity; (ii) how to measure the dynamics of the temperature increase during cyclic operation at a given duty cycle; and (iii) how to model the correlation between the operating temperature and the sensing response in the presence of humidity. Due to the complexity of the 3D analysis of packaged GMOS, and the many aspects of humidity simultanoesuly affecting performane, advanced simulation software is applied, incorporating computational fluid dynamics (CFD). The simulation and experimental data of this study show that the GMOS sensor can operate in the presence of high humidity.

Infusion Simulation of Graphene-Enhanced Resin in LCM for Thermal and Chemo-Rheological Analysis

The present numerical study proposes a framework to determine the heat flow parameters-specific heat and thermal conductivity-of resin-graphene nanoplatelets (GNPs) (modified) as well as non-modified resin (with no GNPs). This is performed by evaluating the exothermic reaction which occurs during both the filling and post-filling stages of Liquid Composite Moulding (LCM). The proposed model uses ANSYS Fluent to solve the Stokes-Brinkman (momentum and mass), energy, and chemical species conservation equations to a describe nano-filled resin infusion, chemo-rheological changes, and heat release/transfer simultaneously on a Representative Volume Element (RVE). The transient Volume-of-Fluid (VOF) method is employed to track free-surface propagation (resin-air interface) throughout the computational domain. A User-Defined Function (UDF) is developed together with a User-Defined Scaler (UDS) to incorporate the heat generation (polymerisation), which is added as an extra source term into the energy equation. A separate UDF is used to capture intra-tow (microscopic) flow by adding a source term into the momentum equation. The numerical findings indicate that the incorporation of GNPs can accelerate the curing of the resin system due to the high thermal conductivity of the nanofiller. Furthermore, the model proves its capability in predicting the specific heat and thermal conductivity of the modified and non-modified resin systems utilising the computed heat of reaction data. The analysis shows an increase of ∼15% in the specific heat and thermal conductivity due to different mould temperatures applied (110-170 °C). This, furthermore, stresses the fact that the addition of GNPs (0.2 wt.%) improves the resin-specific heat by 3.68% and thermal conductivity by 58% in comparison to the non-modified thermoset resin. The numerical findings show a satisfactory agreement with and in the range of experimental data available in the literature.

On the comparison between pre- and post-surgery nasal anatomies via computational fluid dynamics

Nasal breathing difficulties (NBD) are widespread and difficult to diagnose; the failure rate of their surgical corrections is high. Computational fluid dynamics (CFD) enables diagnosis of NBD and surgery planning, by comparing a pre-operative (pre-op) situation with the outcome of virtual surgery (post-op). An equivalent comparison is involved when considering distinct anatomies in the search for the functionally normal nose. Currently, this comparison is carried out in more than one way, under the implicit assumption that results are unchanged, which reflects our limited understanding of the driver of the respiratory function. The study describes how to set up a meaningful comparison. A pre-op anatomy, derived via segmentation from a CT scan, is compared with a post-op anatomy obtained via virtual surgery. State-of-the-art numerical simulations for a steady inspiration carry out the comparison under three types of global constraints, derived from the field of turbulent flow control: a constant pressure drop (CPG) between external ambient and throat, a constant flow rate (CFR) through the airways and a constant power input (CPI) from the lungs can be enforced. A significant difference in the quantities of interest is observed depending on the type of comparison. Global quantities (flow rate, pressure drop and nasal resistance) as well as local ones are affected. The type of flow forcing affects the outcome of the comparison between pre-op and post-op anatomies. Among the three available options, we argue that CPG is the least adequate. Arguments favouring either CFR or CPI are presented.

Computational fluid dynamics-based design and in vitro characterization of a novel pediatric pump-lung

BACKGROUND

Although extracorporeal membrane oxygenation (ECMO) has been used to provide temporary support for pediatric patients suffering severe respiratory or cardiac failure since 1970, ECMO systems specifically designed for pediatric patients, particularly for long-term use, remain an unmet clinical need. We sought to develop a new pediatric ECMO system, that is, pediatric pump-lung (PPL), consisting of a unique cylinder oxygenator with an outside-in radial flow path and a centrifugal pump.

METHODS

Computational fluid dynamics was used to analyze the blood fluid field for optimized biocompatible and gas exchange performances in terms of flow characteristics, hemolysis, and gas transfer efficiency. Ovine blood was used for in vitro hemolysis and gas transfer testing.

RESULTS

Both the computational and experimental data showed that the pressure drop through the PPL's oxygenator is significantly low, even at a flow rate of more than 3.5 L/min. The PPL showed better hemolysis performance than a commercial ECMO circuit consisting of the Quadrox-iD pediatric oxygenator and the Rotaflow pump at a 3.5 L/min flow rate and 250 mm Hg afterload pressure. The oxygen transfer rate of the PPL can reach over 200 mL/min at a flow rate of 3.5 L/min.

CONCLUSIONS

The PPL has the potential to provide adequate blood pumping and excellent respiratory support with minimal risk of hemolysis for a wide range of pediatric patients.

CFD-based bioreactor model with proportional-integral-derivative controller functionality for dissolved oxygen and pH

A physics-based model for predicting cell culture fluid properties inside a stirred tank bioreactor with embedded PID controller logic is presented. The model evokes a time-accurate solution to the fluid velocity field and overall volumetric mass transfer coefficient, as well as the ongoing effects of interfacial mass transfer, species mixing, and aqueous chemical reactions. The modeled system also includes a direct coupling between process variables and system control variables via embedded controller logic. Satisfactory agreement is realized between the model prediction and measured bioreactor data in terms of the steady-state operating conditions and the response to setpoint changes. Simulation runtimes are suitable for industrial research and design timescales.

CFD evaluation of hydrophobic feedstock bench-scale fermenters for efficient high agitation volumetric mass transfer

A new biomanufacturing platform combining intracellular metabolic engineering of the oleaginous yeast Yarrowia lipolytica and extracellular bioreaction engineering provides efficient bioconversion of plant oils/animal fats into high-value products. However, predicting the hydrodynamics and mass transfer parameters is difficult due to the high agitation and sparging required to create dispersed oil droplets in an aqueous medium for efficient yeast fermentation. In the current study, commercial computational fluid dynamic (CFD) solver Ansys CFX coupled with the MUSIG model first predicts two-phase system (oil/water and air/water) mixing dynamics and their particle size distributions. Then, a three-phase model (oil, air, and water) utilizing dispersed air bubbles and a polydispersed oil phase was implemented to explore fermenter mixing, gas dispersion efficiency, and volumetric mass transfer coefficient estimations (kL a). The study analyzed the effect of the impeller type, agitation speed, and power input on the tank's flow field and revealed that upward-pumping pitched blade impellers (PBI) in the top two positions (compared to Rushton-type) provided advantageous oil phase homogeneity and similar estimated kL a values with reduced power. These results show good agreement with the experimental mixing and kL a data.

The influence of blood velocity and vessel geometric parameters on wall shear stress

Vascular geometry was proposed to be one risk factor of atherosclerosis (AS). When developing this hypothesis, the discussion of geometry-wall shear stress (WSS) has often been included. However, further exploration on how various geometric parameters were affecting WSS was needed. The purpose of this study was to investigate the influence degree of vessel geometric parameters and blood velocity on WSS. A computational fluid dynamics (CFD) analyses of the vertebral and basilar arteries (VA and BA, respectively) was used. Twenty patients with no plaques or vessel wall thickening at the VA and BA were included. CFD analyses using both specific vessel models and flow conditions measured by ultrasound Doppler were performed. Subsequently, CFD results were post-processed with multiple linear regression to investigate numerical correlations between geometrical and flow parameters and WSS. The results of the multiple linear regression analysis further demonstrated that the BA proximal velocity was the most influential factor positively influencing BA WSS. The lower the WSS was, the stronger the influence brought by BA average diameter would be. The regression demonstrated that the contributions brought by average diameter and proximal velocity in lower WSS regions were lower than that in higher WSS regions. Tortuosity was only positively correlated with 97.5th WSS percentile, and vessel length and curvature showed no correlation with WSS. This study quantified the influence degree of BA morphology and flow velocity on WSS, which may have practical implications for predicting hemodynamic risks.

Experimental and numerical analysis of the thermal performance of pebble solar thermal collector

In this work, pebbles of higher specific heat than the conventional absorber materials like aluminium or copper are proposed as a absorber in the solar flat plate collector. The proposed collector are integrated into the building design and constructed with masonry. Tests were conducted by varying the operating parameters which influence its performance, like the flow rate of the heat-absorbing medium, and the tilt of the collector using both coated and uncoated pebbles. The maximum temperature difference that could be measured for a conventional absorber was approximately 8 °C for a flow rate of 0.6 L/min. While for a coated and uncoated absorber, it was 7 °C and 5.5 °C respectively. This difference decreased with an increase in flow rates from 0.6 L/min to 1.2 L/min. For all the flow rates, it was observed that the average difference in efficiency between the coated and the conventional absorber collector is 5.82 %, while the difference between the coated and uncoated absorber collector is 15.68 %. Thus, it is very much evident that by replacing the conventional absorber with the proposed coated pebble absorber, the overall loss in efficiency is just 5.82 %, but the advantages are enormous. Along with the experimental study, numerical analysis was also carried out with CFD modeling. The numerical results agreed well with experimental results with the least error. Therefore, CFD simulation can be further used to optimize the design of the collector.

Study of aerosol dispersion and control in dental practice

OBJECTIVES

In this study, we investigated the dispersion patterns of aerosols and droplets in dental clinics and developed a suction device to evaluate its effectiveness in reducing aerosols during dental procedures.

MATERIALS AND METHODS

Firstly, the continuous images of oral aerosols and droplets were photographed with a high-speed camera, and the trajectories of these particles were recognized and processed by Image J to determine key parameters affecting particle dispersion: diffusion velocity, distance, and angle. Secondly, based on the parameter data, the flow field of aerosol particles around the oral cavity was simulated using computational fluid dynamics (CFD), and the flow field under adsorption conditions was simulated to demonstrate the aerodynamic characteristics and capture efficiencies of the single-channel and three-channel adsorption ports at different pressures. Finally, according to the simulated data, a three-channel suction device was developed, and the capture efficiency of the device was tested by the fluorescein tracer method.

RESULTS

The dispersion experimental data showed that aerosol particles' maximum diffusion velocity, distance, and angle were 6.2 m/s, 0.55 m, and 130°, respectively. The simulated aerosol flow-field distribution was consistent with the aerosol dispersion patterns. The adsorption simulation results showed that the outlet flow rate of single-channel adsorption was 184.5 L/s at - 350 Pa, and the aerosol capture efficiency could reach 79.4%. At - 350 Pa and - 150 Pa, the outlet flow rate of three-channel adsorption was 228.9 L/s, and the capture efficiency was 99.23%. The adsorption experimental data showed that the capture efficiency of three-channel suction device was 97.71%.

CONCLUSIONS

A three-channel suction device was designed by simulations and experiments, which can capture most aerosols in the dental clinic and prevent them from spreading.

CLINICAL RELEVANCE

Using three-channel suction devices during oral treatment effectively reduces the spread of oral aerosols, which is essential to prevent the spread of epidemics and ensure the health and safety of patients and dental staff.

The Slump Flow of Cementitious Pastes: Simulation vs. Experiments

Understanding the transient properties of cementitious pastes is crucial for construction materials engineering. Computational modeling, particularly through Computational Fluid Dynamics (CFD), offers a promising avenue to enhance our understanding of these properties. However, there are several numerical uncertainties that affect the accuracy of the simulations using CFD. This study focuses on evaluating the accuracy of CFD simulations in replicating slump flow tests for cementitious pastes by determining the impact of the numerical setup on the simulation accuracy and evaluates the transient, viscosity-dependent flows for different viscous pastes. Rheological input parameters were sourced from rheometric tests and Herschel-Bulkley regression of flow curves. We assessed spatial and temporal convergence and compared two regularization methods for the rheological model. Our findings reveal that temporal and spatial refinements significantly affected the final test results. Adjustments in simulation setups effectively reduced computational errors to less than four percent compared to experimental outcomes. The Papanastasiou regularization was found to be more accurate than the bi-viscosity model. Employing a slice geometry, rather than a full three-dimensional cone mesh, led to accurate results with decreased computational costs. The analysis of transient flow properties revealed the effect of the paste viscosity on the time- and shear-dependent flow progress. The study provides an enhanced understanding of transient flow patterns in cementitious pastes and presents a refined CFD model for simulating slump flow tests. These advancements contribute to improving the accuracy and efficiency of computational analyses in the field of cement and concrete flow, offering a benchmark for prospective analysis of transient flow cases.

A hemodynamic analysis of energy loss in abdominal aortic aneurysm using three-dimension idealized model

Objective: The aim of this study is to perform specific hemodynamic simulations of idealized abdominal aortic aneurysm (AAA) models with different diameters, curvatures and eccentricities and evaluate the risk of thrombosis and aneurysm rupture. Methods: Nine idealized AAA models with different diameters (3 cm or 5 cm), curvatures (0° or 30°) and eccentricities (centered on or tangent to the aorta), as well as a normal model, were constructed using commercial software (Solidworks; Dassault Systemes S.A, Suresnes, France). Hemodynamic simulations were conducted with the same time-varying volumetric flow rate extracted from the literature and 3-element Windkessel model (3 EWM) boundary conditions were applied at the aortic outlet. Several hemodynamic parameters such as time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), endothelial cell activation potential (ECAP) and energy loss (EL) were obtained to evaluate the risk of thrombosis and aneurysm rupture under different conditions. Results: Simulation results showed that the proportion of low TAWSS region and high OSI region increases with the rising of aneurysm diameter, whereas decreases in the curvature and eccentric models of the corresponding diameters, with the 5 cm normal model having the largest low TAWSS region (68.5%) and high OSI region (40%). Similar to the results of TAWSS and OSI, the high ECAP and high RRT areas were largest in the 5 cm normal model, with the highest wall-averaged value (RRT: 5.18 s, ECAP: 4.36 Pa-1). Differently, the increase of aneurysm diameter, curvature, and eccentricity all lead to the increase of mean flow EL and turbulent EL, such that the highest mean flow EL (0.82 W · 10-3) and turbulent EL (1.72 W · 10-3) were observed in the eccentric 5 cm model with the bending angle of 30°. Conclusion: Collectively, increases in aneurysm diameter, curvature, and eccentricity all raise mean flow EL and turbulent flow EL, which may aggravate the damage and disturbance of flow in aneurysm. In addition, it can be inferred by conventional parameters (TAWSS, OSI, RRT and ECAP) that the increase of aneurysm diameter may raise the risk of thrombosis, whereas the curvature and eccentricity appeared to have a protective effect against thrombosis.

Fluid dynamic analysis in predicting the recanalization of intracranial aneurysms after coil embolization - A study of spatiotemporal characteristics

Purpose

Hemodynamics play a key role in the management of cerebral aneurysm recanalization after coil embolization; however, the most reliable hemodynamic parameter remains unknown. Previous studies have explored the use of both spatiotemporally averaged and maximal definitions for hemodynamic parameters, based on computational fluid dynamics (CFD) analysis, to build predictive models for aneurysmal recanalization. In this study, we aimed to assess the influence of different spatiotemporal characteristics of hemodynamic parameters on predictive performance.

Methods

Hemodynamics were simulated using CFD for 66 cerebral aneurysms from 65 patients. We evaluated 14 types of spatiotemporal definitions for two hemodynamic parameters in the pre-coiling model and five in virtual post-coiling model (VM) created by cutting the aneurysm from the pre-coiling model. A total of 91 spatiotemporal hemodynamic features were derived and utilized to develop univariate predictor (UP) and multivariate logistic regression (LR) models. The model's performance was assessed using two metrics: the area under the receiver operating characteristic curve (AUROC) and the area under the precision-recall curve (AUPRC).

Results

Different spatiotemporal hemodynamic features exhibited a wide range of AUROC values ranging from 0.224 to 0.747, with 22 feature pairs showing a significant difference in AUROC value (P-value <0.05), despite being derived from the same hemodynamic parameter. PDave,q1 was identified as the strongest UP with AUROC/AUPRC values of 0.747/0.385, yielding sensitivity and specificity value of 0.889 and 0.614 at the optimal cut-off value, respectively. The LR model further improved the prediction performance, having AUROC/AUPRC values of 0.890/0.903. At the optimal cut-off value, the LR model achieved a specificity of 0.877, sensitivity of 0.719, outperforming the UP model.

Conclusion

Our research indicated that the characteristics of hemodynamic parameters in terms of space and time had a significant impact on the development of predictive model. Our findings suggest that LR model based on spatiotemporal hemodynamic features could be clinically useful in predicting recanalization after coil embolization in patients, without the need for invasive procedures.

Calibration method of the k-ω SST turbulence model for wind turbine performance prediction near stall condition

The present study intends to enhance the accuracy of the k-ω SST turbulence model for numerical wind turbine simulation in stall condition. In order to achieve this purpose a calibration approach is proposed, and is applied to NREL S826 NTNU wind rotor. This method consists in adjusting the two turbulence model coefficients: a1 and β*, which are found to be 0.8 and 0.45 respectively. Then the calculated power and thrust coefficients are compared to the experimental results. The power coefficient results revealed that the simulation relative error at the range of tip speed ratio between 3 to 6 where light stall occurs, is reduced from 17.89 % to 4.58 % by application of calibration. The effects of calibration on flow behaviour are implemented thereafter, by analysing pressure, and skin friction coefficients distribution along the blade. The limiting streamlines on the blade suction side are examined for more flow behaviour understanding. The effects on turbulent kinetic energy around the blade are also highlighted. The main important conclusions that can be made is that calibration reduce the separation zone on the blade suction side, and limits the vortex shedding strength, leading to improve the rotor efficiency and hence to improve the model accuracy.

Accelerated Computational Fluid Dynamics Simulations of Microfluidic Devices by Exploiting Higher Levels of Abstraction

The design of microfluidic devices is a cumbersome and tedious process that can be significantly improved by simulation. Methods based on Computational Fluid Dynamics (CFD) are considered state-of-the-art, but require extensive compute time-oftentimes limiting the size of microfluidic devices that can be simulated. Simulation methods that abstract the underlying physics on a higher level generally provide results instantly, but the fidelity of these methods is usually worse. In this work, a simulation method that accelerates CFD simulations by exploiting simulation methods on higher levels of abstraction is proposed. Case studies confirm that the proposed method accelerates CFD simulations by multiple factors (often several orders of magnitude) while maintaining the fidelity of CFD simulations.

Numerical Investigation of Dimensionless Parameters in Carangiform Fish Swimming Hydrodynamics

Research into how fish and other aquatic organisms propel themselves offers valuable natural references for enhancing technology related to underwater devices like vehicles, propellers, and biomimetic robotics. Additionally, such research provides insights into fish evolution and ecological dynamics. This work carried out a numerical investigation of the most relevant dimensionless parameters in a fish swimming environment (Reynolds Re, Strouhal St, and Slip numbers) to provide valuable knowledge in terms of biomechanics behavior. Thus, a three-dimensional numerical study of the fish-like lambari, a BCF swimmer with carangiform kinematics, was conducted using the URANS approach with the k-ω-SST transition turbulence closure model in the OpenFOAM software. In this study, we initially reported the equilibrium Strouhal number, which is represented by St∗, and its dependence on the Reynolds number, denoted as Re. This was performed following a power-law relationship of St∝Re(-α). We also conducted a comprehensive analysis of the hydrodynamic forces and the effect of body undulation in fish on the production of swimming drag and thrust. Additionally, we computed propulsive and quasi-propulsive efficiencies, as well as examined the influence of the Reynolds number and Slip number on fish performance. Finally, we performed a vortex dynamics analysis, in which different wake configurations were revealed under variations of the dimensionless parameters St, Re, and Slip. Furthermore, we explored the relationship between the generation of a leading-edge vortex via the caudal fin and the peak thrust production within the motion cycle.

Optimal design and development of a fast steering robot inspired by scallops

The improvement of the steering performance of jet robots is challenging due to single inflexible jet aperture. Scallops provide a potential solution with hard shells and a double-hole jet propulsion, which are expected to achieve fast steering movement under water. Inspired by scallops, a bionic propulsion dynamic mesh is proposed in this article, and a three-dimensional computational model of scallops is established. We further calculated the scallop propulsion mechanism under the swing of shells with different shapes. The coupling of simultaneous swing of two shells and their coupling with velum are presented, revealing the unique movement mechanism of Bivalvia. Based on this, the advantages of the double-hole jet propulsion are applied to develop a scallop robot with excellent steering capabilities. Experiments are conducted to verify the steering performance of the scallop robot.

Investigation of ultrasound-assisted starch acetylation by single- and dual- frequency ultrasound based on rheology modelling, non-isothermal reaction kinetics, and flow/acoustic simulation

To achieve wheat starch acetylation (AC) with a high degree of substitution (DS), the acetylation process was carried out using various ultrasonication frequencies, including 25 kHz, 40 kHz, and 25 + 40 kHz. In the second step, wheat starch's ultrasound-assisted acetylation (UAA) is simulated using various approaches including the rheology models, non-isothermal reaction kinetics, and flow/acoustic modelling. The computational fluid dynamics (CFD) simulation solves the non-linear acoustic governing equation to determine the flow field and the amount of delivered ultrasound energy. The acetylated starch increased peak and final viscosity, with the highest values observed for the 25 + 40 kHz frequency than other single frequencies (25 kHz and 40 kHz). The viscosity of the starch is specified based on the experimental data using Herschel-Bulkley, power law, and Casson rheology models. According to differential scanning calorimetry (DSC) analysis, the gelatinization parameters and enthalpy of gelatinization (ΔHgel), were found to be lower in acetylated starches at the frequency of 25 + 40 kHz compared to those at frequencies of 25 kHz and 40 kHz, as well as native starches (NS). Moreover, the gelatinization process is examined by implementing the non-isothermal reaction kinetics to obtain the activation energy and reaction order. Based on the results obtained, implementing sonication at 25 kHz reduces the activation energy by 70.3 % compared to native starch. However, the same parameter is obtained to be 69.9 % and 67.1 % for the application of 40 and 25 + 40 kHz transducers, respectively. Additionally, during the sonication treatment, the yield shear stress increases between 24.1 and 31.8 %, based on the applied frequency. Morphology analysis determined by scanning electron microscopy (SEM) revealed that the surfaces and small granules underwent more damage in acetylated starches at frequencies of 25 kHz and 40 kHz. However, in acetylated starches at 25 + 40 kHz, the larger granules were more affected than the smaller ones.

Development of an effective two-equation turbulence modeling approach for simulating aerosol deposition across a range of turbulence levels

Pharmaceutical aerosol systems present a significant challenge to computational fluid dynamics (CFD) modeling based on the need to capture multiple levels of turbulence, frequent transition between laminar and turbulent flows, anisotropic turbulent particle dispersion, and near-wall particle transport phenomena often within geometrically complex systems over multiple time scales. Two-equation turbulence models, such as the k - ω family of approximations, offer a computationally efficient solution approach, but are known to require the use of near-wall (NW) corrections and eddy interaction model (EIM) modifications for accurate predictions of aerosol deposition. The objective of this study was to develop an efficient and effective two-equation turbulence modeling approach that enables accurate predictions of pharmaceutical aerosol deposition across a range of turbulence levels. Key systems considered were the traditional aerosol deposition benchmark cases of a 90-degree bend (R e = 6,000 ) and a vertical straight section of pipe (R e = 10,000 ), as well as a highly complex case of direct-to-infant (D2I) nose-to-lung pharmaceutical aerosol delivery from an air-jet dry powder inhaler (DPI) including a patient interface and infant nasal geometry through mid-trachea (500 < R e < 7,000 ). Of the k - ω family of models, the low Reynolds number (LRN) shear stress transport (SST) approach was determined to provide the best agreement with experimental aerosol deposition data in the D2I system, based on an improved simulation of turbulent jet flow that frequently occurs in DPIs. Considering NW corrections, a new correlation was developed to quantitatively predict best regional values of the y + l i m i t , within which anisotropic NW turbulence is approximated. Considering EIM modifications, a previously described drift correction approach was implemented in pharmaceutical aerosol simulations for the first time. Considering all model corrections and modifications applied to the D2I system, regional relative errors in deposition fractions between CFD predictions and new experimental data were improved from 19-207% (no modifications) to 2-15% (all modifications) with a notable decrease in computational time (up to ∼15%). In conclusion, the highly efficient two-equation k - ω models with physically realistic corrections and modifications provided a viable, efficient and accurate approach to simulate the transport and deposition of pharmaceutical aerosols in complex airway systems that include laminar, turbulent and transitional flows.

Numerical evaluation and experimental validation of fluid flow behavior within an organ-on-a-chip model

BACKGROUND AND OBJECTIVE

By combining biomaterials, cell culture, and microfluidic technology, organ-on-a-chip (OoC) platforms have the ability to reproduce the physiological microenvironment of human organs. For this reason, these advanced microfluidic devices have been used to resemble various diseases and investigate novel treatments. In addition to the experimental assessment, numerical studies of biodevices have been performed aiming at their improvement and optimization. Despite considerable progress in numerical modeling of biodevices, the validation of these computational models through comparison with experimental assays remains a significant gap in the current literature. This step is critical to ensure the accuracy and reliability of numerical models, and consequently enhance confidence in their predictive results. The aim of the present work is to develop a numerical model capable of reproducing the fluid flow behavior within an OoC, for future investigations, encompassing the geometry optimization.

METHODS

In this study, the validation of a numerical model for an OoC microfluidic device was undertaken. This comprised both quantitative and qualitative assessments of trace microparticles flowing through a physical OoC model. High-speed microscopy images of the flow, using a blood analog fluid, were analyzed and compared with the numerical simulations run using the Ansys Fluent software. For a qualitative analysis, the particles' paths through the inlet and bifurcations were observed whereas, for a quantitative analysis, the particle velocities were measured. Furthermore, oxygen transport was simulated and evaluated for different Reynolds numbers.

RESULTS

In both qualitative and quantitative analyses, the results predicted by the numerical model and the ones outputted by the experimental model were in good agreement. These findings underscore the capability and potential of the developed numerical model. The examination of oxygen transport at various vertical positions within the organoid has revealed that for lower positions, oxygen transport predominantly occurs through diffusion, leading to a symmetric distribution of oxygen. Contrastingly, the convection phenomenon becomes more evident in the upper region of the organoid.

CONCLUSIONS

The successful validation of the numerical model against experimental data shows its accuracy and reliability in simulating the fluid flow within the OoC, which consequently can expedite the OoC design process by reducing the need for prototypes' fabrication and costly laboratory experiments.

Experimental and numerical investigation of the stenosed coronary artery taken from the clinical setting and modeled in terms of hemodynamics

The study was carried out to investigate the effect of the artery with different pulse values and stenosis rates on the pressure drop, the peristaltic pump outlet pressure, fractional flow reserve (FFR) and most importantly the amount of power consumed by the peristaltic pump. For this purpose, images taken from the clinical environment were produced as models (10 mm inlet diameter) with 0% and 70% percent areal stenosis rates (PSR) on a three-dimensional (3D) printer. In the experimental system, pure water was used as the fluid at 54, 84, 114, 132, and 168 bpm pulse values. In addition, computational fluid dynamics (CFD) analyzes of the test region were performed using experimental boundary conditions with the help of ANSYS-Fluent software. The findings showed that as PSR increases in the arteries, the pressure drop in the stenosis region increases and this amount increases dramatically with increasing effort. An increase of approximately 40% was observed in the pump outlet pressure value from 54 bpm to 168 bpm in the PSR 0% model and 51% increase in the PSR 70% model. It has been observed that the pump does more work to overcome the increased pressure difference due to increased pulse rate and PSR. With the effect of contraction, the power consumption of the pump increased from 9.2% for 54 bpm to 13.8% for 168 bpm. In both models, the Wall Shear Stress (WSS) increased significantly. WSS increased abruptly in the stenosis and arcuate regions, while sudden decreases were observed in the flow separation region.

Computational Study of Aerodynamic Effects of the Dihedral and Angle of Attack of Biomimetic Grids Installed on a Mini UAV

In this paper, a numerical analysis of a biomimetic unmanned aerial vehicle (UAV) is presented. Its wings feature three grids at the tip similar to the primary feathers of birds in order to modify the lift distribution over the wing and help in reducing the induced drag. Numerical analysis using computational fluid dynamics (CFD) is presented to analyze the aerodynamic effects of the changes in dihedral and angle of attack (with respect of the rest of the wing) of these small grids at the tip. The aerodynamic performances (lift, drag, and efficiency) and rolling capabilities are obtained under different flight conditions. The effects of changing the dihedral are small. However, the change in the grid angle of attack increases aerodynamic efficiency by up to 2.5 times when the UAV is under cruise flight conditions. Changes to the angle of attack of the grids also provide increased capabilities for rolling. Finally, boundary values of the pressure coefficient and non-dimensional velocity contours are presented on the surfaces of the UAV, in order to relate the aerodynamic results to the aerodynamic patterns observed over the wing.