Multiphysics Software Research Articles

Thermal and irreversibility analysis of non-Newtonian fluid within trapezoidal cavity containing heated cylinder

In this article, the mixed convection transport of the Carreau fluid in the presence of a magnetic field and a heated circular cylinder within a trapezoidal cavity is examined. Finite element analysis is used via COMSOL Multiphysics software to simulate the two-dimensional flow. The heat equation is formulated by incorporating the effects of viscous dissipation and Joule heating. The upper and lower surfaces of the cavity are adiabatic, while the left and right sides are assumed to be cold. The temperature distribution and flow field within the cavity are visualized through the plotting of isotherms, streamlines, and total entropy generation. The graphical representations are performed to examine the impacts of physical parameters on temperature, flow field, and total entropy generation, including the Reynolds number , Weissenberg number ( 1.0 ≤ We ≤ 6.0 ) , Power law index ( 0.1 ≤ n ≤ 5.0 ) , Hartman number ( 1.0 ≤ M ≤ 50 ) , and Richardson number ( 0.2 ≤ Ri ≤ 5.0 ) and ( 0.2 ≤ Rc ≤ 0.8 ) . It is found that the Joule heating and Richardson number are the main reasons to increase the thermal performance within the cavity. The r ( 100 ≤ Re ≤ 500 ) ising values of power law index and Weissenberg numbers caused to reduce the velocity field, irreversible effects and circulation within cavity. While the entropy is augmented by elevating the Eckert and Hartmann numbers. Furthermore, heat transfer improved with the increased magnetic field and Eckert number.

Investigation of internal thermal distribution in inverted perovskite solar cell and improvement of its heat dissipation performance

Abstract COMSOL Multiphysics software was used to construct a numerical opto-electro-thermal coupling model to investigate the mechanisms of internal heat generation, conduction, and dissipation in inverted (p-i-n architecture) perovskite solar cells (PSCs). The research results indicate that Joule heating and Shockley-Read-Hall (SRH) recombination are the primary sources of heat, leading to significant accumulation of heat at the interfaces between the perovskite and the electron transport layer (ETL), as well as between the ETL and the electrode. This concentration of heat not only affects the performance of the device but also poses challenges for overall thermal management. Therefore, we compared four different top electrode materials (Ag, Cu, Al, and reduced graphene oxide) to assess their performance in terms of heat dissipation efficiency. The results showed that reduced graphene oxide (RGO) performed exceptionally well in heat dissipation efficiency, primarily due to its high thermal conductivity, which enables it to effectively reduce heat accumulation at the interfaces, thereby improving performance of PSCs. This finding provides important material selection criteria for optimizing the thermal management of PSCs.

Effect of soft magnetic particles content on multi-physics field of magnetorheological composite gel clutch with complex flow channel excited by Halbach array arrangement

In this paper we investigate the influence of soft magnetic particle mass fraction on the multi-physics field of proposed magnetorheological (MR) gel clutch. A novel MR gel clutch with complex cup-shaped gap is described, whose performance is based on the relative placement between Halbach array arrangement and MR gel. Smart MR gel with 40 wt%, 50 wt%, 60 wt%, 70 wt% and 80 wt% of soft magnetic particles used as the transfer medium and Halbach array is adopted as magnetic excitation system. Its magnetic/mechanical analysis is carried out based on Bingham-plastics model using COMSOL Multiphysics software, which takes into account the variation of dynamic viscosity with magnetic flux density. The distribution of the magnetic flux density, shear yield stress, post-yield viscosity, shear stress and torque in the four flow channels during the transition from engagement state to disengagement state are obtained and analyzed in detail. Multi-physics field characteristics of proposed MR gel clutch with five kinds of MR gels are studied and compared in order to give some useful suggestions in the design phase. Finally, the dynamic torque of the MR clutch with different MR gel is experimentally evaluated.

Exploration of micro-miniature venturi pump for continuous glucose monitoring through interstitial fluid (ISF) extraction

Interstitial Fluid (ISF) is a naturally occurring body fluid that surrounds cells and tissues and contains certain markers useful for continuously monitoring glucose levels in the blood without constantly extracting it. Miniaturizing can reduce its cost and timescale. In this study, for continuous glucose monitoring the numerical analysis of a Polymethyl methacrylate (PMMA) based Micro-miniature venturi pump was carried out to extract Interstitial Fluid (ISF). The numerical study was carried out using COMSOL Multiphysics software by employing the k-ε turbulence model to obtain output as a vacuum at the throat of the pump. The micro-miniature venturi pump was numerically analyzed for several key parameters, including the dimensions of the throat, the convergent angle, and the divergent angle. This study obtained a minimum pressure of 49.41 kPa numerically and 70 kPa experimentally when an input pressure of 252 kPa of compressed air was given to the system. Further influence of turbulence kinetic energy and turbulence dissipation rate on throat pressure output have been notified. Notably, the pressure output reaches its minimum value of 49.41 kPa when the turbulent kinetic energy and the turbulent dissipation rate reach maximum values of 5.0 × 109 m2/s3 and 7000 m2/s2 respectively. As per the previous research findings, for continuous glucose monitoring through ISF extraction, a pressure level of 95 kPa is deemed adequate hence the designed micro-miniature venturi pump is suitable for this application.

Spiral shaped highly sensitive rectangular PCF-based cancer cells detector in terahertz regime

Abstract The increasing fatality rate of cancer in the modern world has become a global concern. To address this issue, early detection of cancer cells in human blood is essential. In this article, a rectangular spiral shaped photonic crystal fiber has been proposed for sensing cancer cells in terahertz regime. The proposed PCF based sensor features a square shaped core surrounded by rectangular air holes of different sizes in a spiral manner. Finite element method based COMSOL Multiphysics software has been used for the design and analysis of the model. The model is sensitive to refractive index variation for the identification of cancer cells. The values of effective refractive index, effective mode area, relative sensitivity, non-linearity, and confinement loss are obtained for MCF-7, MDA-MB-231, Jurkat, HeLa, Basal and PC-12. The PCF based sensor resulted in excellent results of relative sensitivity and Confinement loss as good as 99.377% and 8.59 × 10−12 dBm−1 respectively. Furthermore, the present PCF based sensor model is fabrication feasible, and we hope that it can be utilized at industrial level for biosensing applications in terahertz waveguide regime.

Numerical simulation of an innovative design of solar air heater using spiral duct inserts for enhancing thermal efficiency

ABSTRACT An innovative design of solar air heater (SAH) with spiral shape is numerically analyzed in this paper. The spiral-shaped SAH hydro-thermal performance is studied for different air mass flow rates using CFD COMSOL Multiphysics software. The outlet air temperature difference, thermal efficiency, flow structure, thermal field and pressure drop are the main parameters which are considered to estimate the fluid flow and heat transfer behaviors of the SAH. The simulation is realized for 1000 W/m2 of solar heat flux and 25°C of ambient temperature as typical climatic conditions in a middle day. The results show that the maximum outlet temperature difference of the SAH reached 93.08°C with 0.0083 kg/s of air mass flow rate. For the range of 0.0083–0.0332 kg/s, the thermal efficiency is varied from 91.41% to 99.01%. The ambitious percentage of the thermal efficiency of 99.01% was observed for 0.012 kg/s. The proposed design provides high thermal efficiency of SAH compared to other configurations of SAH.

Evaluating the Accuracy of the COMSOL-Based Finite-Element Method for Simulating Plasmon-Modified Fluorescence.

Accurately modeling plasmon-modified fluorescence is important for understanding and guiding the design of experimental nanostructures that reliably enhance fluorescence. They are of particular interest due to their potential to allow localized "hot spots" of high fluorescence enhancement in a reproducible manner. Given the increasingly prevalent use of the COMSOL Multiphysics software package for simulating these phenomena, we investigate its accuracy using an analytically tractable model consisting of a gold nanosphere interacting with either a plane wave or a radiating point dipole. COMSOL simulation results were compared with a formally exact analytical theory. It was found that simulation parameters commonly used for plane-wave scattering do not necessarily produce accurate results for the nanoparticle-plasmon-coupled dipole emission case. Instead, user-input adaptive meshing parameters were found to be helpful in achieving quantitative agreements between COMSOL and analytical theory results for plasmon-modified fluorescence. Our studies suggest convergence to analytically calculated values when a minimum of two additional user-input mesh elements separate the point-dipole position and the nanoparticle surface. This practical insight is expected to aid in the application of COMSOL simulations to planning and interpreting fluorescence modification experiments.

Heat and mass transfer during conductive drying of washing phosphate sludge

This study investigates the influence of drying temperature on the conductive drying behavior of washing phosphate sludge. A mathematical model, based on Fourier's Fick diffusion law and the heat transfer equation, was developed to simulate coupled heat and mass transfer during the drying process, incorporating the shrinkage effect. The sludge samples, obtained from the Khouribga phosphate plant in Morocco, were subjected to drying at temperatures of 70°C and 95°C. The experimental setup included continuous monitoring of temperature, humidity, and mass loss throughout the process. COMSOL Multiphysics software was used to solve the nonlinear equations governing heat and mass transfer. Results indicate that increasing the drying temperature significantly accelerates the drying rate. The model successfully predicts the moisture content and temperature distribution over time, aligning well with experimental data. Statistical validation, using the determination coefficient (R²) and root mean square error (RMSE), confirms that the model accurately describes the drying behavior, with R² values of 0.9957 and 0.9967 at 70°C and 95°C, respectively. The inclusion of the shrinkage effect enhances the model's precision in describing the drying kinetics under varying conditions.

Tri-objective and multi-parameter geometric optimization of two-stage radioisotope thermoelectric generator based on NSGA-II

Radioisotope thermoelectric generator (RTG) is one of ideal power sources for deep space exploration due to the long service life and high stability. To further extend service life and improve electrical output, this study is first to simultaneously optimize electrical output and mechanical stability induced by thermal stress of the two-stage RTG using non-dominated sorting genetic algorithm-II. Skutterudite and bismuth telluride are used as the thermoelectric materials in hot-stage and cold-stage thermoelectric generators, respectively. The maximum von Mises stress of the skutterudite-stage (VonSKD) and bismuth telluride-stage (VonBT) thermoelectric generator and the conversion efficiency of the two-stage RTG (ηtwo) are set to the three objectives, and twelve geometrical parameters of two-stage RTG are used as variables. Electrical output and mechanical stability are simulated with COMSOL 6.0 Multiphysics software based on the finite element method. Results show that the ranges of ηtwo, VonSKD, and VonBT in the Pareto front are 4.81–7.57 %, 112.2–228.46 and 19.24–87.33 MPa, respectively, when the temperature difference is 300 K. TOPSIS decision-making method is used to select the specific optimal solution on the Pareto front according to the weight of the different objectives. Comparisons between tri- and bi-objective as well as single-stage optimization are also carried. Results show that there are strong influences among these three objectives. The performance of one objective would improve at the expense of the performance of the remaining objectives decrease if a larger weight was given. The VonSKD has the greatest importance among the three objectives, followed by ηtwo, and the smallest is VonBT. Skutterudite-stage thermoelectric generator should be particularly considered in the design of two-stage RTGs. The two-stage RTGs with different weight factors exhibit different electrical performances due to different geometrical parameters. This research can provide references for two-stage thermoelectric generators in more scenarios to improve electrical performances and extend service life.

Design optimization for enhancing performances of integrated planar inductor for power electronics applications

Goal. In this work, the performance of an integrated planar inductor with a square geometric shape using different materials for the substrate: Ni-Fe, Mn-Zn and Ni-Zn have been analyzed and investigated in order to assess the impact of the substrate material on the performance of the integrated planar inductor and to determine the optimal material in the various applications of power modules in power electronics. Methods. To this end, we carried out an in-depth analysis of the geometric dimensions of the integrated planar inductor, by calculating all the geometric parameters of the proposed structure, to establish an equivalent physical model of the integrated planar inductor in order to evaluate its different electrical specifications. The numerical simulation, based on the three-dimensional mathematical model of the system using Maxwell’s equations, was realized by COMSOL Multiphysics software. Results show the importance of the substrate material for the performance of the integrated planar inductor, and specify that the use of Ni-Fe ferrite as a substrate of the integrated planar inductor gives very interesting performance compared to other materials studied. The presented results provide valuable information on the influence of substrate material on the performance of embedded integrated planar inductor and can help to design and optimize these components for use in power electronic systems. Practical value. These results are significant for a wide range of applications, where the integrated planar inductor performance and efficiency can have a significant impact on the overall performance and cost-effectiveness of the power electronic device. References 37, tables 2, figures 7.

Method for prediction magnetic silencing of uncertain energy-saturated extended technical objects in prolate spheroidal coordinate system

Aim. Development of method for prediction by energy-saturated extended technical objects magnetic silencing based on magnetostatics geometric inverse problems solution and magnetic field spatial spheroidal harmonics calculated in prolate spheroidal coordinate system taking into account of technical objects magnetic characteristics uncertainties. Methodology. Spatial prolate spheroidal harmonics of extended technical objects magnetic field model calculated as magnetostatics geometric inverse problems solution in the form of nonlinear minimax optimization problem based on near field measurements for prediction far extended technical objects magnetic field magnitude. Nonlinear objective function calculated as the weighted sum of squared residuals between the measured and predicted magnetic field COMSOL Multiphysics software package used. Nonlinear minimax optimization problems solutions calculated based on particle swarm nonlinear optimization algorithms. Results. Results of prediction extended technical objects far magnetic field magnitude based on extended technical objects spatial prolate spheroidal harmonics of the magnetic field model in the prolate spheroidal coordinate system using near field measurements with consideration of extended technical objects magnetic characteristics uncertainty. Originality. The method for prediction by extended technical objects magnetic cleanliness based on spatial prolate spheroidal harmonics of the magnetic field model in the prolate spheroidal coordinate system with consideration of magnetic characteristics uncertainty is developed. Practical value. The important practical problem of prediction extended technical objects magnetic silencing based on the spatial prolate spheroidal harmonics of the magnetic field model in the prolate spheroidal coordinate system with consideration of extended technical objects magnetic characteristics uncertainty solved. References 48, figures 2.

Performance evaluation of different new channel box photovoltaic thermal systems

Current research has highlighted three major gaps in the literature: (a) a rarity of studies on photovoltaic thermal systems with channel boxes, (b) a lack of data on pressure drop and surface temperature distribution in these systems, and (c) absence of analysis regarding the impact of materials (metals and polymers) used in photovoltaic thermal systems on performance, cost, payback time and weight. The novelty of this paper lies in addressing these gaps by studying a new channel-box photovoltaic thermal system, varying the number of channels and the height of the collectors that distribute water to these channels, in order to identify the most efficient design in terms of electrical, and thermal efficiency as well as pressure drop. In addition, the impact of different heat exchanger materials (metals and polymers) on the performance, cost, payback time, and weight of photovoltaic thermal systems was investigated. 3D numerical simulations were carried out using COMSOL Multiphysics software, based on the finite element method, and validated by outdoor experimental research. The results show that the photovoltaic thermal system 8 design performs best, with a total efficiency of 89.42%. In addition, the aluminum photovoltaic thermal system stands out as the most cost-effective option, with a payback time of 0.823 y (300 d) and annual savings of $459.262, while being 38% lighter than the copper system. Based on this study, it is recommended that the photovoltaic thermal system 8 be implemented because of its exceptional performance and ease of implementation.

An Algorithm for Creating a Synaptic Cleft Digital Phantom Suitable for Further Numerical Modeling

One of the most significant applications of mathematical numerical methods in biology is the theoretical description of the convectional reaction–diffusion of chemical compounds. Initial biological objects must be appropriately mimicked by digital domains that are suitable for further use in computational modeling. In the present study, an algorithm for the creation of a digital phantom describing a local part of nervous tissue—namely, a synaptic contact—is established. All essential elements of the synapse are determined using a set of consistent Boolean operations within the COMSOL Multiphysics software 6.1. The formalization of the algorithm involves a sequence of procedures and logical operations applied to a combination of 3D Voronoi diagrams, an experimentally defined inner synapse area, and a simple ellipsoid under different sets of biological parameters. The obtained digital phantom is universal and may be applied to different types of neuronal synapses. The clear separation of the designed domains reveals that the boundary’s conditions and internal flux dysconnectivity functions can be set up explicitly. Digital domains corresponding to the parts of a synapse are appropriate for further application of the derived numeric meshes, with various capacities of the included elements. Thus, the obtained digital phantom can be effectively used for further modeling of the convectional reaction–diffusion of chemical compounds in nervous tissue.

Numerical simulation analysis of effect of borehole gas extraction under fluid‐structure interaction

AbstractTo study the effect of coal pressure on gas extraction under fluid‐solid coupling, the effect of gas extraction radius under different pore sizes and treatment methods was simulated by using COMSOL Multiphysics software, the time evolution law of gas pressure on coal surface was obtained, and the numerical simulation results were verified by the observation of flow rate and concentration in experimental mine combined with fluid‐solid interaction. The results show that the distance of extraction radius increases with the increase of borehole diameter, and the relation between extraction time and extraction radius is a power function, but the increase gets smaller and smaller until it becomes zero. For Φ98 mm borehole and Φ120 mm borehole, hydraulic treatment can increase the efficiency of gas extraction by 31.3% and 22.7%, respectively. For hydraulic treatment and conventional treatment, the ratio of gas drainage effect by enlarging hole size is 6.3% and 13.8% respectively. Compared with the areas without gas extraction under the four conditions, the descending speed of gas pressure from fast to slow is Φ120 mm hydraulic flushing treatment, Φ98 mm hydraulic flushing treatment, Φ120 mm conventional treatment, Φ98 mm conventional treatment. Compared with four different conditions, after 180 days of extraction the coal gas pressure decreased by 75.3% within reasonable hole spacing. At the same time, in multi‐hole pumping, the influence area of adjacent borehole is larger than that of single‐hole pumping, and the spacing of borehole should be less than twice the radius of pumping.

RF NEMS switches based on graphene for low pull-in voltage and excellent RF performance

This study introduces a novel graphene RF NEMS capacitive switch and conducts an extensive analysis of its RF performance using CST and COMSOL Multiphysics software. The switch’s characteristic impedance matching is accomplished through a coplanar waveguide (CPW) structure, incorporating a dielectric layer topped with conductive graphene beam electrodes. Simulation results reveal that the monolayer graphene RF NEMS switch maintains an insertion loss of 0.003 ~ 0.015 dB and an isolation greater than 40 dB across an ultra-wideband (UWB) frequency range of DC ~ 140 GHz. Moreover, the switch demonstrates a switching time of 4.03 ps at a pull-in voltage of 1 V. The performance of multilayer graphene RF NEMS switches was also evaluated. The study further investigates the operational mechanism of the graphene RF NEMS switch and performs finite element analysis on the proposed design. The graphene RF NEMS switch discussed herein offers significant advantages, including minimal size, lightweight, cost-effectiveness, low pull-in voltage, rapid response time, excellent RF performance, and reduced space consumption in circuitry. Its potential applications span a range from L to F bands, encompassing data storage, communication systems, sensors, remote sensing and military uses.

Temperature Analysis of Single busbar With Variation of Total Harmonic Distortion in Current

Typically, busbars are used to transmit and distribute currents in a bus duct system. The alarming use of nonlinear loads in the industrial sector or at residentials, such as arc welding, computers, ballast lighting, variable speed drives, and so on, has resulted in the generation of harmonics in current distortion, which are uncontrolled and thus increase heat generation within the system. The research conducted in this paper focuses on the prediction of the heat distribution as well as the analysis on operating temperature of a single busbar with compliance to the British National and International Standard (BS 159: 2014) using the Finite Element Method (FEM) in COMSOL Multiphysics software. The copper busbar dimension used for this research was 20mm x 6mm x 300mm, and the fundamental Root Mean Square (RMS) current was 419.1 A. The size of this busbars can withstand the maximum current of 430 A at a maximum operating temperature of 90°C, which complies with the standard requirement. The fundamental current is injected with variation of total harmonic distortion in current up to 55% with an interval of 5%. According to the findings, the operating temperature increases in direct proportion to the increase in total harmonic distortion with the current injections. With the presence of 55% of total harmonics in the current, the current was increased up to 57.73 A from the fundamental current, while the operating temperature was increased up to 140C from the fundamental temperature. The total harmonics in current produced by the nonlinear loads could affect the operating temperature of the busbars, and this continuous operation of current flow will affect the busbars' lifespan due to the occurrence of overheating.

Electrical properties of acetone imprinted hematite nanomaterials doped with Pd & Ag for gas sensing and simulation of their wireless devices

This article presents a novel technique for wireless hydrothermally improved gas sensor devices enhanced by molecular imprinting technique (MIT) for hematite nanomaterials with the existence of acetone using different ratios (10 %, 15 %, and 25 %) as well as doping with palladium/ silver to detect the imprinted gas. α-Fe2O3 are characterized utilizing field emission scanning electron microscope (FESEM), photoluminescence, Fourier transforms infrared spectroscopy (FTIR), X-ray diffraction (XRD), X-ray photo-spectroscopy (XPS), and High-Resolution Transmission Electron Microscopy (HRTEM), Thermogravimetric Analysis (TGA) and utilizes Brunauer-Emmett-Teller (BET) to ascertain their shape, optical properties, and values for surface area and thermal characteristics, respectively. Depending on the preparation conditions, non-regular cubes and rods have typical particle sizes ranging from 25 to 120 nm, acetone imprinted and doped with palladium sample (Ap) offers smaller particles than acetone imprinted sample (A15) making adsorption easier. Where XRD showed all diffraction peaks for α-Fe2O3 as well as XPS authorized the oxidation state, also FTIR showed the characteristic peaks of the stretching mode Fe-O which indicates the formation of α-Fe2O3. Lower defects for (Ap) as it has the highest intensity of 176.15 a.u. in the photoluminescent spectrum, for TGA analysis results explain that pure sample has more thermal stability and lower weight loss. Six hematite films are fabricated using the spin coating technique where acetone molecules are presented during the synthesis process of the nanomaterials to imprint their shapes on the surface of hematite. By measuring the response of the gas sensors and their electrical properties, the I-V curve for Ap showed a rectifying behavior and its shunt resistance (Rsh) is higher than series resistance (Rs) which ensures a high response of (115 %). To develop a wireless gas sensor COMSOL Multiphysics software 5.3(a) software is used for the simulation of three devices model by depositing a layer by specified dimensions of hematites (A15, As, and Ap) with their different electrical conductivities on the surface of rectangular patch antenna by showing geometry, the microstrip patch antenna bandwidth, resonant frequency, scattering parameters, and radiation patterns in the E-plane and 3-D far-field pattern and directivity is calculated also, Sp gives improved reflection coefficient, gain quality factor and directivity due to its high electrical conductivity.

ВЕРИФИКАЦИЯ РАСПРЕДЕЛЕНИЯ ТЕМПЕРАТУРНЫХ ПОЛЕЙ ПРИ МОДЕЛИРОВАНИИ ПРОЦЕССА ЛАЗЕРНОЙ ОБРАБОТКИ БИМЕТАЛЛА ВТ1-0+МН45

The results of modeling the distribution of temperature fields during laser treatment of bimetal VT1-0+MN45 are presented. It is shown that the use of a finite element model in the COMSOL Multiphysics software package makes it possible to estimate the distribution of the temperature field in an explosion-welded bimetal VT1-0 + + MN45 when a laser beam passes over its surface.

Citric acid functionalized chitosan hydrogel beads for enhanced Cr(VI) removal: Experimental assessment and COMSOL modelling of adsorption in a fixed-bed system

Chitosan hydrogel beads (CHB) and cross-linked CHB functionalized with tricarboxylic citric acid (CA-GLA-CHB) were synthesized and used for the removal of Cr(VI) from aqueous solutions in batch and dynamic adsorption systems. Applicability of COMSOL Multiphysics software as a tool for predictive modelling of column dynamics using a minimum set of experimental parameters was tested and compared with experimental results. Batch experiments demonstrated that CA-GLA-CHB showed almost 20% higher adsorption capacity and wider operational pH range compared to non-functionalized CHB. In a dynamic system, CA-GLA-CHB was used as column packing material and the column dynamics was investigated at different conditions. The breakthrough and exhaustion time of the column increased with an increase in bed height and a decrease in feed flow rate and initial Cr(VI) concentration. Computational modeling by COMSOL software with implemented Advection-Dispersion-Reaction (ADR) equation was used to simulate breakthrough curves for the applied dynamic system in order to validate the possibility of its application for process design in the scaled-up fixed-bed column systems for wastewater treatment. The simulated breakthrough curves of Cr(VI) adsorption matched well with the experimental, thus proving that COMSOL Multiphysics software can be used for scaling-up fixed-bed column systems packed with low-cost and effective CA-GLA-CHB bioadsorbent.

Process intensification in the direct contact membrane distillation (DCMD) desalination by patterning membrane surface: A CFD study

In this research, the impact of pattern geometry on the intensification of the performance of surface patterned membranes for saline water desalination by direct contact membrane distillation (DCMD) is investigated by computational fluid dynamics (CFD) simulation. The Comsol Multiphysics software was applied to solve the governing transport equations for heat, momentum, and mass transfer. The result of the model was validated by the published experimental data, and the maximum deviation was less than 10%. The target was to maximize the permeate flux by altering surface pattern geometry and dimensions. Based on the previous studies, a prism pattern was chosen in this work, and the influences of pattern type (3 types), pattern dimension (25-150 µm valley depth), and the distance between the valleys (0-400 µm) were studied on the temperature polarization coefficient (TPC) and DCMD permeate flux. The results showed that the pattern with a valley depth of 25 µm and a distance between the valleys of 300 µm had the best performance in DCMD operation. In this situation and feed temperature of 80°C, a TPC of 0.78 and a water flux of 49.3 kg/m2.h were attained. The characteristics of flow close to the patterned membrane surface were also investigated, and it was observed that there are weak shear stresses in the lower zone of the valleys, while stronger shear stresses are created in the upper regions that are responsible for improving the TPC and water flux in the patterned membranes.