Computational Fluid Dynamics Model Research Articles

Effect of Hydrodynamic Shear Stress on Algal Cell Fate in 3D Extrusion Bioprinting

The 3D bioprinting of aquatic photosynthetic organisms holds potential for applications in biosensing, wastewater treatment, and biofuel production. While algae cells can be immobilized in bioprinted cell‐friendly matrices, there is a knowledge gap regarding the thresholds of hydrodynamic shear stress that affect the cells’ functionality and viability during bioprinting. This study examines the effect of hydrodynamic shear stress on the fate of Chlamydomonas reinhardtii cells. Computational fluid dynamics models based on the Navier–Stokes equations are developed to numerically predict the shear stresses experienced by the cells during extrusion. Parallelly, cell culture experiments are conducted to evaluate the functionality, growth rates, and viability of algae cells within bioprinted constructs. By correlating cell culture and simulation results, the causal link between shear stress in the nozzle and cell viability and function has been characterized. The findings highlight that cell viability and function are significantly impacted by process factors. Notably, algae cell function is more sensitive to shear stress than cell viability. Functional impairments occur at maximum shear stresses around 5 kPa, while viability remains unaffected. Beyond 14 kPa, both functionality and viability decline significantly and irreversibly. The results emphasize the importance of assessing viability and function after bioprinting, rather than just viability.

How Irregular Geometry and Flow Waveform Affect Pulsating Arterial Mass Transfer.

Alzheimer's disease is a progressive degenerative condition that has various levels of effect on one's memory. It is thought to be caused by a buildup of protein in small fluid-filled spaces in the brain called perivascular spaces (PVS). The PVS often takes on the form of an annular region around arteries and is used as a protein-clearing system for the brain. To analyze the modes of mass transfer in the PVS, a digitized scan of a mouse brain PVS segment was meshed and used for computational fluid dynamics (CFD) studies. Tandem analyses were then carried out and compared between the mouse PVS section and a cylinder with commensurate dimensionless parameters and hydraulic resistance. The geometry pair was used to first validate the CFD model and then assess mass transfer in various advection states: no-flow, constant flow, sinusoidal flow, sinusoidal flow with zero net solvent flux, and an anatomically correct asymmetrical periodic flow. Two mass transfer situations were considered, one being a protein build-up and the other being a protein blend-down using a multitude of metrics. Bulk arterial solute transport was found to be advection-controlled. The consideration of temporal evolution and trajectories of contiguous protein bolus volumes revealed that flow pulsation was beneficial at bolus break-up and that additional local wall curvature-based geometry irregularities also were. Using certain measures, local solute peak concentration blend-down appeared to be diffusion-dominated even for high Peclet numbers; however, bolus size evolution analyses showed definite advection support.

Computational study on the impact of gasoline-ethanol blending on autoignition and soot/NOx emissions under low-load gasoline compression ignition conditions

In the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine are performed to investigate the impact of gasoline-ethanol blending on autoignition, nitrogen oxide (NOx), and soot emissions under low-load conditions. In order to represent the test gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF)+ethanol (ETPRF) surrogate (with 10% ethanol by volume; E10) is employed. A three-dimensional (3D) engine CFD model employing finite-rate chemistry with a skeletal kinetic mechanism (including NOx sub-mechanism), adaptive mesh refinement (AMR), and hybrid method of moments (HMOM) is adopted to capture the in-cylinder combustion phenomena and soot/NOx emissions. The engine CFD model is validated against experimental data for three gasoline-ethanol blends: E10, E30 and E100, with varying ethanol content by volume. Model validation is carried out for a broad range of start-of-injection (SOI) timings (−21, −27, −36, and −45 crank angle degrees (°CA) after top-dead-center (aTDC)) with respect to in-cylinder pressure, heat release rate, combustion phasing, NOx and soot emissions. For relatively later injection timings (−21 and −27 °CA aTDC), E30 yields higher amount of soot than E10; while the trend reverses for early injection cases (−36 and −45 °CA aTDC ). On the other hand, E100 yields the lowest amount of soot among all fuels irrespective of SOI timing. Further, E10 shows a non-monotonic trend in soot emissions with SOI timing: SOI-36>SOI-45>SOI-21>SOI-27, while soot emissions from E30 exhibit monotonic decrease with advancing SOI timing. NOx emissions from various fuels follow a trend of E10>E30>E100. On the other hand, NOx emissions increase as SOI timing is advanced for all fuels, with an anomaly for E10 and E100 where NOx decreases when SOI is advanced beyond −36 °CA aTDC. Detailed analysis of the numerical results is performed to investigate the soot/NOx emission trends and elucidate the impact of chemical composition and physical properties on autoignition and emissions characteristics.

Characterization of the dynamic viscosity of cell cultures and its effect on mixing performance in a spinner flask bioreactor

Computational fluid dynamics (CFD) models have been developed to simulate cell culturing bioreactors but assume water-like viscosity properties due to significant data gaps. This study characterized the dynamic viscosity of HEK-293 cell cultures and evaluated its effect on mixing performance in a spinner flask bioreactor. Viscosity measurements indicated that the cell culture media, media with microcarriers, and cell cultures presented shear thinning behaviors within the measured shear rate range of 1–100 s−1. The viscosity also increased with the microcarrier concentrations and growth of cell culture. The CFD model, incorporating dynamic viscosity data, showed that shear stress and Kolmogorov length profiles are significantly influenced by microcarrier concentrations and cell culture growth. Higher microcarrier concentrations led to higher average shear stress and Kolmogorov values. The cultured HEK-293 cells after seven days of growth also had higher average shear stress and Kolmogorov values than at the day of seeding, indicating an impact caused by the cells’ metabolism and biomass. Overall, the results indicated that assuming water-like properties underestimates shear stress and Kolmogorov length scales, especially at zones of lower shear rates due to the observed shear thinning behavior. Thus, careful monitoring of dynamic viscosity of cell cultures and proper control of mixing parameters are critical to deliver the desired mixing conditions for optimized cell growth especially during scale-up production operations.

Evaluation navigation controlled gate of aging spillway on cavitation damage.

High-speed flow of clean water or water with sediment released from aging spillways may cause abrasion and cavitation on the concrete surface gradually. The occurrence of irregularities on the concrete surface can exacerbate the erosion problem. Which might jeopardize the safety of dams constantly, hence the rehabilitation efforts become urgent tasks in dam safety projects. This study employs a 3D Computational Fluid Dynamics (CFD) model to quantitatively analyze the cavitation risk on the aging concrete surface of the Chay 5 spillway in Ha Giang, Vietnam, under various operation scenarios. There are two standards used to measure cavitation: the cavitation index (σ) which indicates the danger due to the drop of pressure in rapid flow, and the new gasification index (β) which takes into consideration the formation and collapse of bubbles behind asperities. Three extreme flood cases may not result in potential cavitation because both σ and β exceed critical thresholds. Regarding the six controlled gate scenarios with normal water level, the σ profiles are approximated 1,0 showing a low likelihood of cavitation damage while the β values are smaller than 0.8, indicating a considerable risk of cavitation. Besides, the opening height of 100 cm poses the greatest risk of creating severe cavitation erosion in the concave area and slope portion. The flip bucket experienced the most vulnerable cavitation when the opening height is 400 cm. In addition, an approach to spillway surface rehabilitation involving specialized mortars has been presented. For aging conveyance structure, gasification index (β) takes into account irregularities surface, providing a more comprehensive assessment of the likelihood of cavitation damage than cavitation index (σ). After rehabilitation with anti-shrinkage high abrasion resistance mortar, the entire spillway surface is smooth. This allows for reducing the cavitation risk and improvement of life service thereof.

Comparison of experimental and simulated separation performance in capillary tube-in-manifold devices

A metal tube-in-manifold packed bed capillary column device, designed to overcome common limitations associated with capillary LC separations, is described. Experimental results of initial packing tests with sub-3 μm core-shell particles demonstrated efficiencies greater than 47,000 plates/m for a separation performed using the column device. Computational fluid dynamics (CFD) modeling of the multicomponent separation used for this work was validated against experimental LC results and the optimized model was able to effectively predict component peak retention times. However, the accuracy of predicted efficiencies requires further refinement. The tube-in-manifold design demonstrates that packed capillary columns with cylindrical cross-sectional channel geometry and ultrahigh pressure, low dead volume fluidic connections are achievable.

Can the infection risk in elevators be negligible? A comparative study of airborne infection probability in elevators and conference rooms

People in crowded and poorly ventilated elevators are at risk of respiratory infection. However, due to the short duration of the elevator ride, the transmission of respiratory diseases in elevators does not get enough attention. To evaluate the infection risk, this study investigated the airborne transmission of respiratory diseases in the hospital elevator by comparison to the conference room. A validated computational fluid dynamics (CFD) model was adopted to simulate airflow, temperature and tracer gas dispersion in the elevator and conference room. We used Wells-Riley model to evaluate the infection probability of the susceptible persons. In addition, the influences of source patient posture and ventilation strategy on the transmission of tracer gas were analyzed. The results showed that the infection probability in the elevator with 5 min (average 2.70%) was higher than that in the conference room with 50 min (average 1.77%). The effects of source patient posture and ventilation strategy on the infection probability in the elevator were more significant than those in the conference room. Local air circulation could gather the tracer gas inside a confined space in the elevator and led to a high infection probability. The infection risk of respiratory diseases in the elevator was non-negligible.

Multi-objective optimization of chemical reaction characteristics of selective catalytic reduction in denitrification of diesel engine using ELM-MOPSO methodology

This study systematically investigated the effects of structural parameters in a two-stage selective catalytic reduction (SCR) catalyst on ammonia storage, NO conversion efficiency, and ammonia slip using computational fluid dynamics (CFD) simulations. The rank ordering of the correlations between structural parameters and performance metrics was determined by gray relational analysis (GRA). Both the catalyst diameter and the upstream catalyst porosity had correlations exceeding 0.9, making them key parameters affecting SCR performance. An extreme learning machine (ELM) prediction model was trained using data obtained from CFD simulations. The prediction performance of this ELM model was then evaluated using statistical metrics. Finally, the optimized Pareto frontier diagram was obtained through the multi-objective particle swarm optimization (MOPSO) algorithm, and the optimal results were selected for further CFD modeling analysis. The results showed that the upstream NO conversion efficiency of the SCR system was improved by 7.7 %, the downstream NO conversion efficiency was improved by 18.67 %, the ammonia slip was reduced by 71.97 %, and the upstream ammonia storage capacity was improved by 6.83 % compared with the original model.

Estimation of Turbulent Mixing Factor and Study of Turbulent Flow Structures in Pressurized Water Reactor Sub Channel by Direct Numerical Simulation

Abstract Subchannel analysis codes are presently a requirement for design and safety analysis of nuclear reactors. Among the crucial inputs for these codes, the turbulent mixing factor holds particular significance. However, acquiring this factor through experimental means proves to be a challenging endeavor, primarily due to the necessity for precise pressure equilibrium between subchannels. Consequently, this requirement leads to the undertaking of expensive and intricate experiments for each new reactor or in cases where there are modifications in fuel bundle design. The need for direct numerical simulation (DNS) stems from the challenges and costs involved in experimental techniques, and the uncertainties due to empiricism in computational fluid dynamics (CFD) models. In this study, DNS has been conducted across six Reynolds numbers, ranging from 17,640 to 1.176 × 105, in the geometry of a pressurized water reactor (PWR) subchannel. The resulting turbulent flow structures have been computed and their dynamics are examined. Furthermore, this paper presents a methodology for directly calculating the turbulent mixing factor from the fluctuating velocity field obtained from DNS data. The turbulent mixing process has been scrutinized in-depth, and a correlation for the turbulent mixing factor is developed. It is noted that most of the mixing occurs in the near-wall region. The study suggests different mixing factors for mass and momentum mixing. This paper aims to provide a comprehensive insight into the turbulent mixing phenomenon.

Experimental and CFD Analysis of Hydrodynamics in Dual-Impeller Crystallizer at Different Off-Bottom Clearances

Producing tailor-made crystals demands knowledge of the influence of hydrodynamics and nucleation kinetics. In this paper, the hydrodynamic conditions in a dual-impeller crystallizer and their influence on the key nucleation parameters of the batch cooling crystallization of borax at different impeller off-bottom clearances were investigated. Two different impeller configurations were used—a dual pitched-blade turbine (2 PBT) and a dual straight-blade turbine (2 SBT). Hydrodynamics was analyzed in depth based on the developed computational fluid dynamics model. The experimental results on mixing time and power input were used to validate the numerical model. The results show that the properties of the final product are affected by the impeller position in both dual-impeller configurations. An increase in the impeller off-bottom clearance in both systems results in a decrease in the mean crystal size. The hydrodynamic conditions generated at C/D = 1 in the 2 PBT impeller system and at C/D = 0.6 in the 2 SBT impeller system favored an earlier onset of nucleation compared to other impeller positions. It was found that the eddy dissipation rate and the Kolmogorov length scale correlate highly with the mean crystal size, suggesting that the size is affected by the shear stress in the vessel, rather than the overall convective flow.

The Impact of Flow-Thermal Characteristics in Ship-Board Solid Oxide Fuel Cells

Hydrogen is increasingly recognized as a clean and reliable energy vector for decarbonization in the future. In the marine sector, marine solid oxide fuel cells (SOFCs) that employ hydrogen as an energy source have already been developed. In this study, a multi-channel plate-anode-loaded SOFC was taken as the research object. A three-dimensional steady-state computational fluid dynamics (CFD) model for anode-supported SOFC was established, which is based on the mass conservation, energy conservation, momentum conservation, electrochemical reactions, and charge transport equations, including detailed geometric shapes, model boundary condition settings, and the numerical methods employed. The polarization curves calculated from the numerical simulation were compared with experimental results from the literature to verify the model’s accuracy. The curved model was applied by enlarging the flow channels or adding blocks. Numerical calculations were employed to obtain the current density, temperature distribution, and component concentration distribution under the operating conditions of the SOFC. Subsequently, the distribution patterns of various physical parameters during the SOFC operation were analyzed. Compared to the classical model, the temperature of the curved model was reduced by 1.3%, and the velocities of the cathode and anode were increased by 4.9% and 5.0%, respectively, with a 2.42% enhancement in performance. The findings of this study provide robust support for research into and the application of marine SOFCs, and offer they insights into how we may achieving “dual carbon” goals.

A Validated Thermal Computational Fluid Dynamics Model of Wine Warming in a Glass

Oenophiles are aware that the temperature at the time of drinking can profoundly shape wine’s sensory attributes. Wine is usually served and drunk below room temperature but warms up after pouring due to heat exchange with warmer surroundings. This study investigates how quickly wine warms up in a wine glass and identifies the relevant heating effects. A numerical simulation using conjugate heat transfer is established, representing the complex multi-physical process. Experiments are conducted to validate the simulation. It is shown that the simulation must take into account thermal conduction, convection, and even radiation to provide accurate results. Without simulating radiation and convection of the room air, the predicted temperature is off by 66.3% or 3.3 °C. As warming is independent of the alcohol content, the simulation results are valid for non-sparkling wine types with moderate sugar levels within the considered configuration. A parameter study investigated the temperature increase over time depending on the ambient temperature and the initial wine temperature for 150 mL wine in a medium-sized red wine glass. The results can provide information on preparing a wine to obtain the desired drinking temperature.

Estimating daily road traffic pollution in Hong Kong using CFD modelling: Validation and application

Road traffic is a major source of urban air pollution, exposing residents to harmful levels of pollutants, particularly in densely populated urban centers. In this study, Computational Fluid Dynamics (CFD) modelling is applied to investigate the impact of road traffic emissions on air pollution within a dense urban environment. The model, applied in Mong Kok, Hong Kong, incorporates high-resolution traffic activity, meteorological, and background concentration data. It demonstrates strong agreement with roadside measurements, deviating by 18.8 % for NOx, 15.5 % for CO, and 0.7 % for PM10, with potential underestimation linked to unmodeled NOx-O3 photochemistry during midday hours (12:00–17:00). Results reveal that local road traffic significantly increases background NOx levels (41 % average increase), highlighting the importance of dispersion modelling for accurate pollution assessments. Traffic-induced increases in CO and PM10 background levels are less pronounced. Two applications demonstrate the model's utility. Firstly, analysis of daily pollution distribution reveals that approximately 10 % of the studied area experiences moderate to very high pollution levels during peak hours. Secondly, an investigation into pollution sources within a leisure park identifies heavy-duty and light-duty vehicles as the primary NOx polluters. This modelling approach provides a valuable decision-making tool for policymakers and local authorities. It enables the assessment of urban design changes, strategic placement of sensitive structures, and future traffic fleet composition choices to mitigate pollution exposure. As air quality regulations tighten, this study offers a robust method for guiding impactful mitigation strategies.

Computational investigation of combustion and emission characteristics of HCCI engine fueled with Carbon-Free NH3 –H2O2 blend

A computational investigation is provided of the utilization of blends of NH3 and H2O2 as carbonless fuels for HCCI combustion. A Computational Fluid Dynamics (CFD) model is developed and validated to assess the engine’s performance under varying proportions of NH3 and H2O2, equivalence ratios, and inlet temperatures. The heat release rate profiles demonstrated dual peaks, which were rationalized as an initial fuel decomposition followed by a slow thermal explosion and then a fast depletion of the fuel after the generation of a pool of OH radicals through decomposition of H2O2. Findings suggest that an increase in mole fraction of H2O2 advances combustion. Increasing the H2O2 volume fraction also decreases combustion duration, while increasing the IMEP. In-cylinder peak pressure and Maximum Pressure Rise Rate (MPRR) increase with H2O2 mole fraction, leading to reduced thermal efficiency, particularly at higher equivalence ratios, and an increased danger of uncontrolled combustion (knocking). Additionally, higher inlet mixture temperatures were shown to result in an advancement in the start of combustion and a reduction in combustion duration, with a corresponding increase in peak pressure. This points to the possibility of potentially controlling the notoriously difficult-to-control HCCI combustion through intake temperature for NH3/H2O2 fuel blends. However, increased NOx emissions are observed for increased intake temperature and in any case the NOx emissions are orders of magnitude above current regulatory limits, which means that practical engine operation will only be possible with an aftertreatment scheme based on the de-NOx activity of NH3.

Examination of Thermal Dispersion and Airflow within a Refrigerator

A blast freezer, characterized by its capability to diminish the core temperature of cooked food from 100 °C to -18 °C within 270 minutes, constitutes a critical component in this preservation process. This study endeavors to model a blast freezer system employing Computational Fluid Dynamics (CFD) methodologies, subsequently validating the CFD analysis through empirical investigations. The pressure-based k-ε turbulence model is employed to solve the Navier-Stokes and energy equations. The ensuing analyses encompass airflow assessments and temperature evaluations for unloaded and fully loaded blast freezers. Results gleaned from experiments and analyses indicate a temperature escalation within the cabin as it approaches the enclosure walls. Maximum velocities of 31.1 m/s and 26.9 m/s are recorded for unloaded and fully loaded freezers. The average disparity between the CFD and experimental models is computed as -0.7 °C, signifying a close alignment between the simulated and actual outcomes.

Computational fluid dynamics investigation of pesticide spraying by agricultural drones

Precise control over pesticide spraying by agricultural drones requires fundamental research of multiphase flow, heat and mass transfer of liquid droplets in downwash flow. In this study, water was used as a carrier agent in pesticide formulations, in which the liquid droplets exhibit a Newtonian fluid behavior. A general computational fluid dynamics (CFD) model that characterizes propeller motion coupled with droplet evaporation and deposition was developed. The downwash flow generated by the propeller rotation was simulated using a multiple reference frame method along with the SST-k-ω turbulence model. Spraying of liquid droplets was simulated using a discrete phase model. The flow model validation was conducted by comparing the simulated velocities with experimental data, and the spray model validation was completed by a comparison of droplet depositions from CFD predictions and experiments. The results indicated that the droplet deposition fluctuates with an increase of propeller speed, and that both temperature and relative humidity affect the droplet deposition. In addition, the relation between the number of particles and the impact velocity was identified quantitatively, impact velocity with respect to particle diameter was predicted, and the effects of surface inclination, surface roughness, and crosswind speed on droplet deposition were analyzed. Moreover, an in-depth discussion about the intrinsic relation between downwash flow and droplet spraying was conducted.

Investigation on the influence of particle size and porosity on solid phase fluidization and heat transfer in carbon fibre porous media

AbstractIn this paper, the 3D pore structure was reconstructed and solid phase fluidization in porous media was investigated. Based on the two fluid model, a computational fluid dynamics (CFD) model of gas–solid fluidization was established and the influence of particle size and porosity were investigated. When porosity was constant, during the fluidization process, the particle velocity and solid concentration standard deviation gradually decreased, and the bed height increased. At the same time, the smaller the particle size was, the smaller the solid concentration standard deviation was, the faster the bed height increased, and the more the particle temperature decreased. Based on the fixed value of particle size, when studying the effect of porosity on fluidization, it was found that with the enhancement of solid‐phase fluidization, particle velocity, solid concentration standard deviation, and particle temperature decreased, and bed height increased. Moreover, when porosity was large, particle velocity decreased rapidly, and solid concentration standard deviation reached a smaller minimum value, which took a longer time. Within the same time, the particle temperature also decreased less.

Hydrogen enrichment in methanol SI engine at varying injection timing during compression stroke

This study investigates the effect of hydrogen enrichment in a direct-injected methanol-fuelled SI engine, with combustion and emission performance analysed by varying the methanol injection timing from 150° to 60° CA bTDC. A three-dimensional computational fluid dynamic model of the hydrogen-enriched methanol engine was developed to predict the engine performance, in-cylinder, and emission characteristics. The numerical model was solved using three-dimensional Reynolds-Averaged Navier Stokes, the RNG k-epsilon model for turbulence, the O'Rourke and Amsden sub-model for heat transfer, the extended Zeldovich mechanism for nitric oxide emissions, and the Hiroyasu-NSC model for soot emissions. The results indicated that retarding the injection timing resulted in a decrease in the indicated specific CO and soot emissions along with a rise in the indicated specific NOx emissions. Furthermore, hydrogen enrichment with methanol enhanced the hydroxyl radical concentration, reduced the combustion duration, reduced also the indicated specific CO and soot emissions while increasing the indicated specific NOx emissions. The study indicates that hydrogen enrichment could extend the late injection timing limit of methanol by enhancing fuel-air mixing, which improves and controls the combustion process more effectively.

CFD modeling of CO2 fixation by microalgae cultivated in a lab scale photobioreactor

The green microalga Chlorella vulgaris has been studied for efficient fixation of carbon dioxide, CO2, from elevated concentrations of in-feed gas flows. In this work, a Computational Fluid Dynamics (CFD) model for algae cultivation including the dependence on CO2 concentration in liquid has been derived based on cultivation experiments of Chlorella vulgaris. The experiments were performed in a laboratory scale cylindrical stirred tank reactor with a range of enriched CO2 concentrations in-feed (0.04–50 %). The model provided by Béchet et al. (2015) for algae cultivation considering dependence on local light, algae concentration and temperature has been implemented for CFD and modified for comprising CO2 dependence. The model has been verified with experimental results for the algae productivity showing the proper trends for dissolved CO2 concentrations and algae concentrations. Thus, the developed CFD model can serve as a tool for evaluating photobioreactor performance for CO2 capture

Calculating the local and overall view factors of a multi-segment human model

In non-uniform environments, such as spaces with radiant ceilings or floors, or with hot or cold windows, view factors (VFs) of a person towards these surfaces are needed for calculating the heat exchanges between the person and the surroundings. Multi-body segmented human models can simulate the heat exchanges of individual body parts with the environment. Simulating heat transfer at the local body parts is important for assessing comfort in non-uniform environments. These models require VFs for each individual body part. Accurately determining these VFs is time-consuming, often requiring methods such as ray-tracing techniques. A fast calculation method to accurately determine the view factor (VF) for each segment has yet to be investigated. This study employed a 16-segment human model into various architectural radiation scenarios. The surface-to-surface model in computational fluid dynamics was utilized to simulate VFs for the 16 body segments in relation to radiant surfaces at different locations. Through a series of simulations, comprehensive formulas for VFs were derived for each individual at various positions within a room, encompassing 16 body segments and computing parameters for 24 building surfaces. The accuracy of these general formulas was validated by comparing the whole-body VF with results from several available studies. Additionally, this study developed a visualization tool for calculating VF and mean radiant temperatures with different room space dimension and embedded radiant surfaces. The users can use this online tool to easily calculate VFs and mean radiant temperatures for each body parts, and for the whole-body. The tool offers valuable insights for predicting thermal comfort in asymmetric radiant environments and facilitating control over building environments.