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.

Simulation and Experimental Study of a Novel Negative-Pressure Flapper-Nozzle Mechanism

As the manufacturing industry evolves, the significance of control valve positioners in chemical production escalates. The flapper-nozzle system, the heart of control valve positioners, directly influences the linearity of system control. Presently, studies on the flapper-nozzle system primarily focus on dynamic system modeling and computational fluid dynamics simulations. However, traditional flapper-nozzle mechanisms often fail to achieve linear control objectives. This paper proposes a novel negative-pressure nozzle structure to tackle this issue, combining computational fluid dynamics and experimental methods, and considering gas compressibility during high-speed flow. Both simulation and experimental results suggest that the new structure improves the supply air pressure and broadens the linear pressure output range of the system, showing significant potential for practical applications.

Prediction of Hemodynamic-Related Hemolysis in Carotid Stenosis and Aiding in Treatment Planning and Risk Stratification Using Computational Fluid Dynamics

Atherosclerosis affects human health in many ways, leading to disability or premature death due to ischemic heart disease, stroke, or limb ischemia. Poststenotic blood flow disruption may also play an essential role in artery wall impairment linked with hemolysis related to shear stress. The maximum shear stress in the atherosclerotic plaque area is the main parameter determining hemolysis risk. In our work, a 3D internal carotid artery model was built from CT scans performed on patients qualified for percutaneous angioplasty due to its symptomatic stenosis. The obtained stenosis geometries were used to conduct a series of computer simulations to identify critical parameters corresponding to the increase in shear stress in the arteries. Stenosis shape parameters responsible for the increase in shear stress were determined. The effect of changes in the carotid artery size, length, and degree of narrowing on the change in maximum shear stress was demonstrated. Then, a correlation for the quick initial diagnosis of atherosclerotic stenoses regarding the risk of hemolysis was developed. The developed relationship for rapid hemolysis risk assessment uses information from typical non-invasive tests for treated patients. Practical guidelines have been developed regarding which stenosis shape parameters pose a risk of hemolysis, which may be adapted in medical practice.

A hemodynamic analysis of fenestrated physician-modified endograft repair for complicated aortic dissections involving the visceral arteries

OBJECTIVE

The aim of this study is to perform patient-specific hemodynamic simulations of the patients with complicated aortic dissection underwent Physician-modified endograft (PMEG) and evaluate the treatment outcome.

METHOD

12 patient-specific models were reconstructed from computed tomography angiography (CTA) data of 6 patients with complicated aortic dissection before and after the PMEG. 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. Hemodynamic indicators such as time-averaged wall shear stress (TAWSS), relative residence time (RRT) and endothelial cell activation potential (ECAP) were obtained to evaluate the postoperative effect of PMEG.

RESULTS

Comparing with the preoperative models, the flow rates of most visceral arteries were increased in the postoperative models (PSMA = 0.012, PRRA = 0.013, and PLRA = 0.005). Pressure and TAWSS in visceral regions were significantly reduced (PP = 0.003 and PTAWSS = 0.017). With the false lumens (FL) covered by the stent grafts, the average TAWSS level increased in the regions of postoperative abdominal aorta (P = 0.002), and the average RRT and ECAP values decreased significantly (PRRT = 0.02 and PECAP = 0.003).

CONCLUSION

This study shows that PMEG, as a new technique for the treatment of complicated aortic dissection involving the distal tears in the visceral region, can effectively restore the abnormal blood supply of the visceral arteries, reduce the risk of aortic rupture, the formation of aortic dissection aneurysm (ADA), and thrombosis. This corresponds well with clinical retrospective studies and 1-year follow-up outcomes. The findings of this study are of great significance for the development of PMEG.