Simulation Of Mass Transfer Research Articles

Thermal analysis of suspended single droplet evaporation measurements with a coupled lumped parameter model

The measurement data of single droplet evaporation experiments are often biased due to the extra heat input through the fiber suspension and the presence of thermal radiation in hot environments. This encumbers model validation for heat and mass transfer simulations of liquid droplets. In this paper, a thermal analysis of this measurement layout is presented with a coupled lumped parameter model, considering heat conduction through the suspension. The model was validated by experimental data from the literature and good agreements were found. The thermal analysis focused on fiber material and geometry, and thermal radiation properties. Calculations were performed on a broad range of ambient conditions for liquids with different volatility characteristics. Temporal squared droplet diameter- and temperature-profiles, furthermore, droplet stationary evaporation rate were used to characterize vaporization phenomena. The thermal balance of the droplet is dominated by the convective heat rate from the environment in the early stage of evaporation. The effect of heat conduction through the fiber becomes important at the end of the droplet lifetime when the droplet size is decreased. Temperature sensor suspension may seriously bias droplet temperature due to the larger thermal conductivity compared to quartz fiber. Large droplets in high-temperature environments show significant sensitivity to thermal radiation properties, which should be considered in measurements and model validation.

Numerical simulation of heat and mass transfer in heat exchange tubes based on action of interfacial slip

The liquid flow close to the interface has received much attention recently, that is determined by interfacial slip action. In heat exchange tube, the thermal diffusion in the flowing boundary layer is also related to the interfacial flow state. However, the study of convective heat transfer mostly ignores the interfacial slip and medium thermal diffusion. In this study, the variation law of heat transfer efficiency and the actual interfacial concentration under the action of interfacial slip and thermal diffusion in the heat exchange tubes are reported. We find that the heat transfer efficiency of no slip reaches the highest, compared to no slip, that of negative slip significantly reduces (at −0.99 mm, copper: −79.56%; aluminum alloy: −78.60%; carbon steel: −75.85%), and that of positive slip slightly reduces (at 0.58 mm, copper: −16.87%; aluminum alloy: −15.66%; carbon steel: −9.82%). As the slip length increases, the interfacial concentration of positive slip tends to the initial concentration, while the concentration changes rate of negative slip increases (at −0.99 mm, H2SO4: 110.29%, NaCl: 45.90%, C2H6O: −28.37%). These results show that the interfacial slip not only affects heat transfer performance but also interfacial concentration, and further affecting the metal corrosion behavior at the heat transfer interface.

Multirate mass transfer simulation of denitrification in a woodchip bioreactor

Denitrifying woodchip bioreactors (DWBs) have proven to be an efficient nature-based solution for nitrate removal. Modeling DWBs is required for improving their design and operation, but is hindered by the complexity of the modeled system where numerous chemical species and model parameters are needed. Reactions inside the woodchips are different from those at the edges, causing chemical localization (i.e., apparent simultaneous occurrence of incompatible reactions). We used the Multi Rate Mass Transfer (MRMT) approach to overcome these problems when simulating reactive transport processes in a DWB located at Kiruna, Sweden. Besides denitrification, other nitrogen-cycling processes (e.g., nitrification, dissimilatory nitrate reduction to ammonium, anammox) and alternative electron donors (e.g. oxygen, sulfate) were also considered. Biomass concentration is incorporated into the biochemical reaction rates, including growth and decay, to characterize microbial catalysis. We found that the MRMT model: 1) can account for the heterogeneity of the porous woodchips; 2) was capable of reproducing the nitrogen species evolution in the DWB with kinetic parameters from the literature; and 3) allows reproducing localized biochemical reactions (e.g., aerobic reactions on the woodchip edges, near the DWB entrance and anaerobic reactions inside); and 4) reproduces the full denitrification reactions sequence, but with the different reactions occurring in different portions of the woodchip (e.g., nitrate to nitrite near the edges and nitrite to nitrous oxide further inside). The latter observation suggests that increasing woodchip size may reduce the outflow of these undesired species.

Three-Dimensional Simulation of Heat and Mass Transfer by the Hybrid Lattice Boltzmann Method

Three-Dimensional Simulation of Heat and Mass Transfer by the Hybrid Lattice Boltzmann Method

Investigation of a falling film tube bank heat exchanger with baffle design for water recovery applications

Shell and tube heat exchangers (STHE) are essential thermal equipment and widely used in daily life. A novel thermosyphon system called falling-film thermosyphon (FFTS) is introduced and integrated into STHE system, resulting in a better thermal performance. In this study, a rectangular solid tube bank of FFTS bundles with a baffle design is studied. The numerical simulation for heat and mass transfer of the FFTS heat exchanger is developed to predict the condensation rate of the vapor in the flue gas, and a lab-scale prototype is also built up in COMSOL. The prediction is validated with the experimental data from references, and the model's accuracy is verified within 10%-12% error. Also, the Non-dominated Sorting Genetic Algorithm, version 2 (NSGA-II) is implemented to improve the thermal performance of rectangular tube banks in this paper. Several parameters, e.g., baffle number, tube number, and tube space, are optimized. As a result, compact configurations with more baffles are preferred to enhance the performance associated with a high-pressure drop correspondingly. The optimized layout for the lab-scale prototype can increase by 18 to 32% condensation with a pressure loss of less than 200 Pa.

CFD-PBE Simulation of Flow Dynamics and Mass Transfer in Two-Stage Countercurrent Mixer-Settler

Flow condition and mass transfer behavior of liquid–liquid extraction in a two-stage countercurrent extractor are investigated experimentally and numerically. CFD-PBE simulation is adopted to analyze the impact of operating parameters, including impeller speed, residence time, and organic-aqueous ratio of feed flow on volume fraction and droplet size of dispersed phase and volumetric mass transfer coefficient. Droplet size distribution is calculated by considering breakage and coalescence of dispersed phase. The allocation of dispersed phase between stages is found to be determined by impeller speed. Meanwhile, the allocation of dispersed phase between mixers and settlers is governed by residence time and organic-aqueous ratio of feed flow. Turbulence intensity, volume fraction, and droplet size distribution are revealed to be dominating parameters that affect the spatial nonuniform distribution of volumetric mass transfer coefficient. The mass transfer rate can be improved by increasing impeller speed in two stages. Higher residence time and organic-aqueous ratio of feed flow, however, have opposite effects on mass transfer rate in stage-1 and stage-2.

Toward multiscale simulations for solidification microstructure and microsegregation for selective laser melting of nickel-based superalloys

In this study, we carried out multiscale simulations that integrate a macro-scale finite element method simulation for mass and heat transfer of the transient molten pool and a micro-scale quantitative phase-field model for dendritic growth and solute distribution based on selective laser melting (SLM) experiments. Macroscopic simulations reveal an approximately inverse coupling of solidification velocity Vs and thermal gradient G near the solid–liquid interface at the tail end of the molten pool. It was found that the maximum flow velocity of the Inconel 718 alloy melt in the molten pool exceeds 2×10−2 m/s and the maximum negative pressure reaches −5.28×10−4 Pa. Through quantitative phase-field simulations, the map about G/Vs for each solidification morphology was obtained, and it was found that solidification velocity dependent solute partition coefficients k(Vs) can strongly influence the estimated stability regions. The power law relationship of the primary dendritic arm spacing (PDAS) λ1 versus solidification velocity λ1∝Vsα agrees with the previous experimental results for lower G, and the exponent α decreases with increasing G. The PDAS λ1 versus cooling rate R following λ1∝R−0.512 is consistent with previous investigations. We found that λ1 is a nonlinear function with G−0.5Vs−0.25, i.e., λ1∝(G−0.5Vs−0.25)α, with α ranging from 0.887 to 2.466 for various G. Comparison of the predicted λ1 with Trivedi's model reveals that λ1 is still nonlinear with G−0.25Vs−0.25.

Outdoor experimental validation for ultra-high concentrator photovoltaic with serpentine-based cooling system

With demand for renewable energy growing, concentrator photovoltaic thermal hybrids have great potential. Maximising concentration ratios through the deployment of multi-stage optics can yield high power outputs from multi-junction solar cells. To prevent damaging thermal stress and to enable extraction of thermal energy, a capable cooling system is necessary.The primary objective of this study is to maximise the effective concentration ratio over a solar cell and calibrate the system to optimise the energetic and exergetic efficiencies. The capability of the serpentine-based cooling system is investigated for each concentrator optic configuration. Originality is found in the presentation of the 3-stage optic, and the use of outdoor real-world experimental data to validate a computational model. This model uses both ray tracing, heat and mass transfer simulations to enhance the understanding of system operation and enable accurate prediction of performance under various conditions. Results show focal spot shape is more important than raw optical efficiency for electrical output, making the 3-stage optic superior to the other configurations in most regards. An effective concentration of over 1200 × is achieved. Higher exergetic efficiencies are consistently found in the double serpentine configuration, though variation does not exceed ±0.3% when only changing cooling system geometry.

Malachite green dye adsorption by jackfruit based activated carbon: Optimization, mass transfer simulation and surface area prediction

The exposure of basic dye towards human and environment can caused severe adverse effects. Hence, this study aimed (i) to optimize the adsorption of malachite green (MG) dye onto jackfruit peel based activated carbon (JPAC) by using response surface methodology (RSM) and (ii) to simulate the mass transfer process using Polymath mass transfer (PMT) model. The RSM revealed that the optimum activation temperature, activation time and potassium hydroxide (KOH) impregnation ratio, IR to be 782 °C, 1 h and 2.96 g/g which translated into 66.80 mg/g of MG adsorption and 34.52 % of JPAC's yield. The Freundlich isotherm model described the adsorption process the best while the Langmuir monolayer adsorption capacity was computed to be 376.33 mg/g. The kinetic data obeyed the PMT model the best where the rate constant, KPMT decreased from 1.20 to 0.41 h−1 as the MG concentration increased from 25 to 300 mg/L. Also, the PMT model was succeeded in predicting the surface area, aPMT of JPAC to be 807.77 m2/g and this value is comparable with the actual JPAC's mesopores surface area of 714.25 m2/g (error 11.58 %). The adsorption of MG onto JPAC was aided by several functional groups such as primary or secondary alcohol, non-bonded hydroxy group and methoxy through ion-dipole force and dipole-dipole force. Based on thermodynamic parameters, MG-JPAC adsorption system was endothermic in nature (∆H° = 37.12 kJ/mol), had an increment of randomness at solid-liquid interface (∆S° = 0.19 kJ/mol·K), spontaneous in nature (−20.11 kJ/mol) and physisorptionally governed (Ea = 9.46 kJ/mol).

Simulation of Heat and Mass Transfer on Blunt Bodies under Aerodynamic Heating of High-Speed Aircraft

Simulation of Heat and Mass Transfer on Blunt Bodies under Aerodynamic Heating of High-Speed Aircraft

A multilayer heat and mass transfer model for photoacoustics of aerosol particles (MHM-PA)

A diffusion multi-layer model is introduced which allows time-dependent simulations of heat and mass transfer for photoacoustic experiments on aerosol particles (MHM-PA model). It provides spatial and temporal information on the particle temperature, composition and size during photoacoustic modulation cycles, and thus enables accurate simulations of the photoacoustic amplitude and phase at different modulation frequencies. Simulations for micrometer-sized, mixed water-tetraethylene glycol droplets at different relative humidity and at different modulation frequencies between 4 and 100 kHz reveal a substantial loss of the information content of the photoacoustic signal above ∼30 kHz. The photoacoustic amplitude reduces by two orders of magnitude, causing sensitivity issues. The photoacoustic phase converges to the same constant value for all droplet sizes and humidities, meaning that useful information on heat and mass transfer is essentially lost in the phase.

Numerical Simulation of the 3D Simultaneous Heat and Mass Transfer in a Forced Convection Solar Drying System Integrated With Thermal Storage Material

Abstract The demand for quality dried products necessitates cost effective and innovative drying techniques that will improve its market value. The slow drying rate, weather dependency, and moisture reabsorption have been identified as the major challenges of solar drying operation. To address these shortcomings, hybrid solar drying systems have been recommended for the drying of various agricultural materials and other porous products. Designing a better drying system to accommodate thermal storage materials requires detailed analysis, which could be achieved through numerical simulation. Therefore, the numerical simulation of heat and mass transfer in a forced convection solar drying system integrated with black-coated firebrick sensible thermal storage materials (STSM) for the cocoa beans, locust beans, cereal grains, etc., was investigated under no-load conditions. The equations governing the fluid flow for a three-dimensional solar drying system were solved using the finite volume method with the aid of ansys, the computational fluid dynamics software to comprehend the dynamic and thermal behavior of the airflow within the dryer. The experimental maximum temperature values of 96.9 °C and 77.3 °C for the collector and drying chamber were in agreement with the simulated maximum collector and drying chamber temperatures of 116.9 °C and 80 °C respectively. The designed solar drying system with the incorporated STSM showed the capacity of raising the temperature of the air within the drying chamber to 3–37 °C above ambient temperature between 01:00 p.m. and 10:00 p.m. The agreement of the simulated dryer model with the experimental one is an indication that the developed dryer is suitable for drying cocoa, locust beans, fish, cereal grains, and some other agricultural products within an acceptable period based on the previous studies and therefore, the drying system is recommended to avoid the shortcomings associated with traditional/open sun drying.

Computational simulation of mass transfer in membranes using hybrid machine learning models and computational fluid dynamics

In this work, a modeling approach has been developed for describing component separation by selective membrane systems. The developed modeling strategy in this research is based on computational fluid dynamics (CFD) integrated with machine learning (ML) approach to reduce the computational costs and evaluate the possibility of ML model for integration to CFD models. The considered case study is a membrane-based molecular separation for removal of species from water. CFD was performed to solve mass transfer equations, and the data in form of concentration profiles have been used for ML computations. The data set in this research has about two thousand data vectors, each of which contains two inputs r and z and one output C which is the species concentration in the feed channel of membrane. Three models, Multi-Layer Perceptron (MLP), Support Vector Machines with RBF kernel (RBF-SVM), and Decision Tree (DT) were chosen for modeling since they were anticipated to perform well on the provided data by adjusting their hyper-parameters. After optimizing the models, their results were carefully checked with different criteria. All three models showed acceptable results with R2 score criteria higher than 0.9. MAPE decreased to 3.57 × 10−3 and improved to 3.6 with the RMSE criterion. We can introduce the best model in this research as the MLP model for description of the mass transfer in the membrane.

Numerical Simulation of Energy and Mass Transfer in a Magnetic Stirring Photocatalytic Reactor

Hydrogen production via photocatalytic water splitting is one of the promising solutions to energy and environmental issues. Understanding the relationship between hydrogen production in suspended photocatalytic reactions and various influencing factors is crucial for expanding the scale of the system. However, the complexity of physical and chemical factors involved in hydrogen production via photocatalytic water splitting makes systematic research of this technology challenging. In recent research, the simulated light source reactor has become a preferred study object due to its strong controllability. This paper presents a comprehensive energy and mass transfer model for the suspended photocatalytic reaction in a magnetically stirred reactor. The mutual impacts between the flow field, radiation field, and reaction field are analyzed. The simulation results show that the rotating speed of the stirring magneton in the reactor has a significant influence on the flow field. The rotation of the stirring magneton generates a vortex in the central axis area of the reactor, with the relationship between the depth of the vortex f(s) and the rotating speed of the magneton s described as f(s) = 0.27e0.0032s. The distribution of radiation within the reactor is influenced by both the incident radiation intensity and the concentration of the catalyst. The relationship between the penetration depth of radiation g(i) and the incident radiation intensity i is described as g(i) = 10.73ln(i) − 49.59. The relationship between the penetration depth of radiation h(c) and the particle concentration c is given as h(c) = −16.38ln(c) + 15.01. The radiation distribution in the reactor has a substantial impact on hydrogen production, which affects the concentration distribution law of hydrogen. The total amounts of hydrogen generated in the reactor are 1.04 × 10−3 mol and 1.35 × 10−3 mol when the reaction times are 1.0 s and 2.0 s, respectively. This study serves as a foundation for the future scaling of the system and offers theoretical guidance for the optimization of the photocatalytic reactor design and operating conditions.

Influence of aortic aneurysm on the local distribution of NO and O2 using image-based computational fluid dynamics

There is a pressing need to establish novel biomarkers to predict the progression of thoracic aortic aneurysm (TAA) dilatation. Aside from hemodynamics, the roles of oxygen (O2) and nitric oxide (NO) in TAA pathogenesis are potentially significant. As such, it is imperative to comprehend the relationship between aneurysm presence and species distribution in both the lumen and aortic wall. Given the limitations of existing imaging methods, we propose the use of patient-specific computational fluid dynamics (CFD) to explore this relationship. We have performed CFD simulations of O2 and NO mass transfer in the lumen and aortic wall for two cases: a healthy control (HC) and a patient with TAA, both acquired using 4D-flow magnetic resonance imaging (MRI). The mass transfer of O2 was based on active transport by hemoglobin, while the local variations of the wall shear stress (WSS) drove NO production. Comparing hemodynamic properties, the time-averaged WSS was considerably lower for TAA, while the oscillatory shear index and endothelial cell activation potential were notably elevated. O2 and NO showed a non-uniform distribution within the lumen and an inverse correlation between the two species. We identified several locations of hypoxic regions for both cases due to lumen-side mass transfer limitations. In the wall, NO varied spatially, with a clear distinction between TAA and HC. In conclusion, the hemodynamics and mass transfer of NO in the aorta exhibit the potential to serve as a diagnostic biomarker for TAA. Furthermore, hypoxia may provide additional insights into the onset of other aortic pathologies.

Numerical simulation and mathematical modeling for heat and mass transfer in MHD stagnation point flow of nanofluid consisting of entropy generation

The primary goal of this article is to explore the radiative stagnation point flow of nanofluid with cross-diffusion and entropy generation across a permeable curved surface. Moreover, the activation energy, Joule heating, slip condition, and viscous dissipation effects have been considered in order to achieve realistic results. The governing equations associated with the modeling of this research have been transformed into ordinary differential equations by utilizing appropriate transformation variable. The resulting system of equations was solved numerically by using Bvp4c built-in package in MATLAB. The impact of involved parameters have been graphically examined for the diverse features of velocity, temperature, and concentration profiles. Throughout the analysis, the volume fraction is assumed to be less than 5% while the Prandtl number is set to be 6. In addition, the entropy generation, friction drag, Nusselt, and Sherwood numbers have been plotted for describing the diverse physical aspects of the underlying phenomena. The major outcomes reveal that the curvature parameter reduces the velocity profile and skin friction coefficient whereas the magnetic parameter, temperature difference parameter, and radiation parameter intensify the entropy generation.

Transient simulation of heat and mass transfer in a parallel flow regenerative lime kiln with pulverized coal injection

Quicklime, which is an essential raw material in various industries, is primarily produced by thermally calcinating limestone in kilns. Among the various kiln designs, parallel flow regenerative (PFR) kilns with pulverized coal combustion are the predominant technology. The performance of the PFR kilns is determined by the interactions between the gas, limestone, quicklime, and coal particles inside the chamber. Additionally, the two chambers in the PFR kiln alternate after a short period during the operation. Consequently, the transfer behaviors between the various phases inside the chamber are in a transient state. A three-dimensional transient multiphase model was created to investigate the transient multiphase behavior. The entire operation cycle, which consists of calcination and regeneration processes, was simulated using a coupling method. The predicted CO2 concentration in the exhaust gas was compared with the industrial measurement data, and good agreement was obtained, with an error of less than 6%. Furthermore, variations in the gas and limestone temperatures and limestone concentrations with residence time during each period were obtained. During the calcination period, the gas temperature field stabilized after 100 s. After stabilization, the zone where intensive limestone decomposition occurred was located in the range 0–5 m from the lance tip. The maximum gas temperature in the calcination process was 1845 K, which was approximately 700 K higher than that of a PFR kiln fueled with natural gas. CO2 gradually accumulated in the middle part of the chamber; the average concentration of CO2 in the off-gas was 27.7%. In the regeneration period, the transfer process stabilization required 400 s. Furthermore, the maximum gas temperature was higher than that of a PFR kiln fueled with natural gas, with a difference of approximately 90 K. Contrary to the calcination period, CO2 accumulated in the peripheral region at the end of regeneration. The total CO2 emission is 1.2024 t CO2/t product for the investigated coal-fueled PFR kiln. These transient transfer characteristics, during the PFR operation, could aid the optimization of coal utilization and control of CO2 emissions in the future.

Mass transfer simulation on remazol brilliant blue R dye adsorption by optimized teak wood Based activated carbon

The reactive dye of remazol brilliant blue R (RBBR) can be toxic to aquatic ecosystems and humans. Hence, the objectives of this study were to optimize teak wood-based activated carbon (TWAC) via response surface methodology (RSM) to adsorb RBBR and to simulate the mass transfer process using Polymath software. The optimum conditions in preparing TWAC were discovered to be 470 Watt of radiation power, 6.40 min of radiation time, and 1.48 g/g of impregnation ratio (IR), resulting in RBBR uptakes of 83.97 mg/g and TWAC’s yield of 30.60%. Isotherm study revealed that the adsorption of RBBR onto TWAC was best described by the Langmuir model with maximum monolayer adsorption capacity, Qm of 337.55 mg/g. The kinetic data were best described by the Polymath Mass Transfer (PMT) model where the rate constant, kPTM was found to increase from 4.84 to 5.22 h−1 when RBBR initial concentration increased from 25 to 300 mg/L. The PMT model predicted the adsorption surface area, am to be 940.79 m2/g, which is highly accurate with the actual mesopores surface area of 983.25 m2/g. The RBBR molecules formed an attraction force with TWAC’s surface functional groups through hydrogen bond (cyclohexane), dipole–dipole force (polysulfides and terminal/medial alkyne), and ion–dipole force (phenol and hydroxy group). Thermodynamic parameters of ΔH°, ΔS°, ΔG°, and Ea were computed to be −3.55 kJ/mol (exothermic), 0.07 kJ/mol.K, −24.37 kJ/mol (spontaneous) and 21.23 kJ/mol (physisorption), respectively.

Numerical simulation on mass transfer in the bone lacunar-canalicular system under different gravity fields

The bone lacunar-canalicular system (LCS) is a unique complex 3D microscopic tubular network structure within the osteon that contains interstitial fluid flow to ensure the efficient transport of signaling molecules, nutrients, and wastes to guarantee the normal physiological activities of bone tissue. The mass transfer laws in the LCS under microgravity and hypergravity are still unclear. In this paper, a multi-scale 3D osteon model was established to mimic the cortical osteon, and a finite element method was used to numerically analyze the mass transfer in the LCS under hypergravity, normal gravity and microgravity and combined with high-intensity exercise conditions. It was shown that hypergravity promoted mass transfer in the LCS to the deep lacunae, and the number of particles in lacunae increased more significantly from normal gravity to hypergravity the further away from the Haversian canal. The microgravity environment inhibited particles transport in the LCS to deep lacunae. Under normal gravity and microgravity, the number of particles in lacunae increased greatly when doing high-intensity exercise compared to stationary standing. This paper presents the first simulation of mass transfer within the LCS with different gravity fields combined with high-intensity exercise using the finite element method. The research suggested that hypergravity can greatly promote mass transfer in the LCS to deep lacunae, and microgravity strongly inhibited this mass transfer; high-intensity exercise increased the mass transfer rate in the LCS. This study provided a new strategy to combat and treat microgravity-induced osteoporosis.

Non‐Isothermal Mass Transfer in Fluid Drops with Internal Circulation

AbstractThe influence of internal circulation on the heat and mass transfer within dispersed fluid drops in translation was studied numerically using the computational fluid dynamics simulation tool ANSYS FLUENT. The three phases identified for heat/mass transfer are: the initial and final diffusion‐dominated phases and a middle convection‐dominated phase. Non‐isothermal mass transfer simulations using temperature‐dependent diffusion coefficients were carried out to investigate the influence of heating and cooling of the drop on the solute transfer. Mass transfer is slower compared to heat transfer and the final diffusion‐dominated mass transfer occurs at the final temperature attained by the drop. Mass transfer is faster for drop heating than for drop cooling, as the final diffusion phase is rate controlling and significant solute transfer occurs at the final higher temperature.