Element Size Research Articles

Simulation-free determination of microstructure representative volume element size via Fisher scores

A representative volume element (RVE) is a reasonably small unit of microstructure that can be simulated to obtain the same effective properties as the entire microstructure sample. Finite element (FE) simulation of RVEs, as opposed to much larger samples, saves computational expenses, especially in multiscale modeling. Therefore, it is desirable to have a framework that determines the RVE size prior to FE simulations. Existing methods select the RVE size based on when the FE-simulated properties of samples of increasing sizes converge with insignificant statistical variations, with the drawback being that many samples must be simulated. We propose a simulation-free alternative that determines the RVE size based only on a micrograph. The approach utilizes a machine learning model trained to implicitly characterize the stochastic nature of the input micrograph. The underlying rationale is to view RVE size as the smallest moving window size for which the stochastic nature of the microstructure within the window is stationary as the window moves across a large micrograph. For this purpose, we adapt a recently developed Fisher score-based framework for microstructure nonstationarity monitoring. Because the resulting RVE size is based solely on the micrograph and does not involve any FE simulation of specific properties, it constitutes an RVE for any property of interest that solely depends on the microstructure characteristics. Through numerical experiments of simple and complex microstructures, we validate our approach and show that our selected RVE sizes are consistent with when the chosen FE-simulated properties converge.

Direct maximum power injection control of grid-connected PV micro-inverter systems connected to the grid

The increasing use of fossil fuels such as gas and oil has caused many environmental problems, including air pollution, depletion of underground resources, the destructive phenomenon of global warming, destruction of the ozone layer, and many other problems for the environment. Photovoltaic systems are one of the most promising solutions to overcome these problems. Due to the low voltage of solar panels and the higher controllability of photovoltaic systems against changes in sunlight intensity during the day, DC-DC power converters are required. In this project, the aim is to design and build a DC-DC converter with the ability to provide high voltage gain at a suitable efficiency for a 1-string photovoltaic system with the ability to track the maximum power point 2. The designed converter is designed and analyzed with the ability to provide a voltage gain of 5/6 at the critical operating point, a minimum efficiency of 94%, and a power of 500 watts. This converter uses classical soft switching, so the switching losses are reduced, which in turn increases the efficiency of the converter. This converter has a lower voltage stress on the switch (4/0) than the existing techniques. On the other hand, while new converters use a large number of switches and magnetic cores in the cell to provide high gain, this converter uses the least number of switches and magnetic cores to achieve high gain. Also, the passive elements of this converter, such as inductor and capacitor, are in the order of a few microhenries and farads, and therefore the size of these components is also reduced. In addition, due to the high gain provided in lower duty cycles, the efficiency of the switch is increased, and therefore cheaper switches can be used. Also, due to the reduction in the size of passive elements, the size of the magnetic core is also reduced, which in turn optimizes the converter economically.

Small planar antenna array design using length and spacing through Matlab‐HFSS interfacing

SummaryThis study presents a novel method for reducing return loss and achieving lower side lobe levels in small planar antenna arrays (PAAs). Quantum particle swarm optimisation (QPSO) findings are utilised for antenna design modelling and performance assessment utilising the MATLAB‐based moment's code approach and MATLAB‐HFSS interface. The mutual coupling between antenna components is considered and computed using the method of moments (MoM). Changing the array elements' length and spacing will accomplish the abovementioned objectives. The Matlab‐based moments code calculates the antenna components' length and spacing, and the antennas' length and spacing are used to produce the radiation pattern via the Matlab‐HFSS interface. Four examples of small PAA designs are given to illustrate how well the applied methodology works. The method discussed in this paper is generic and can be employed in the future for different sizes of the array elements and the different antenna array structure designs.

Assembly collapsing versus heterozygosity oversizing: detection of homokaryotic and heterokaryotic Laccaria trichodermophora strains by hybrid genome assembly.

Genome assembly and annotation using short-paired reads is challenging for eukaryotic organisms due to their large size, variable ploidy and large number of repetitive elements. However, the use of single-molecule long reads improves assembly quality (completeness and contiguity), but haplotype duplications still pose assembly challenges. To address the effect of read length on genome assembly quality, gene prediction and annotation, we compared genome assemblers and sequencing technologies with four strains of the ectomycorrhizal fungus Laccaria trichodermophora. By analysing the predicted repertoire of carbohydrate enzymes, we investigated the effects of assembly quality on functional inferences. Libraries were generated using three different sequencing platforms (Illumina Next-Seq, Mi-Seq and PacBio Sequel), and genomes were assembled using single and hybrid assemblies/libraries. Long reads or hybrid assemby resolved the collapsing of repeated regions, but the nuclear heterozygous versions remained unresolved. In dikaryotic fungi, each cell includes two nuclei and each nucleus has differences not only in allelic gene version but also in gene composition and synteny. These heterokaryotic cells produce fragmentation and size overestimation of the genome assembly of each nucleus. Hybrid assembly revealed a wider functional diversity of genomes. Here, several predicted oxidizing activities on glycosyl residues of oligosaccharides and several chitooligosaccharide acetylase activities would have passed unnoticed in short-read assemblies. Also, the size and fragmentation of the genome assembly, in combination with heterozygosity analysis, allowed us to distinguish homokaryotic and heterokaryotic strains isolated from L. trichodermophora fruit bodies.

A Structured Mesh Generation Based on Improved Ray-Tracing Method for Finite Difference Time Domain Simulation

A structured mesh generation method based on an improved ray-tracing technique is proposed for finite time difference domain (FDTD) simulation in this paper. This method converts the triangular faces provided by the STL file into the structured mesh required by the FDTD method. The size of the generated structured mesh is determined by the size of the triangular element, so it is non-uniform and can adapt to the shape of the target. In addition, this method proposes an improved ray-tracing technique to avoid the singularities of the object, so it has much better accuracy than the conventional method. Two representative electromagnetic objects, an Su-30 aircraft and a complex bow-tie antenna, representing a typical electrically large model and fine structure, respectively, are used to validate the accuracy of the mesh generation method. The results show that the presented method has high accuracy, regardless of whether it is used for large targets or fine structures such as a small hole, strip, thin layer, etc.

An evaluation of large diameter through-thickness metallic pins in composites

There is increasing demand for functional through-thickness reinforcement (TTR) in composites using elements whose geometry exceeds limitations of existing TTR methods like tufting, stitching, and z-pinning. Recently, static insertion of large diameter TTR pins into heated prepreg stacks has proven a feasible and robust reinforcement process capable of providing accurate TTR element placement with low insertion forces and lower tow damage compared with existing methods for similar element sizes (>1mm diameter) like post-cure drilling. Local mechanical performance and failure mechanics of these pinned laminates are reported here. Laminates with a single statically inserted pins (1.2, 1.5, and 2.0 mm) can mostly retain their in-plane integrity alongside a local improvement in mode I delamination toughness in carbon fibre-benzoxazine laminates. Tensile strength is mostly unaffected by the pins resulting from delamination suppression, whereas there is up to a doubling of Young’s modulus. Compressive strength is significantly diminished (up to 42 %) in pinned laminates. Interlaminar toughness is improved, and peak toughness is pushed ahead of the crack as pin diameter increases. The lack of significant deterioration in in-plane tensile properties in pinned laminates produced using static insertion can expand the range and forms of materials that can be inserted compared to existing TTR.

Research and fabrication of bimetallic plate transducer based on finite element analysis

By enhancing the type 1–3 piezoelectric composite material, a highly sensitive hydroacoustic transducer sensitive element structure is suggested for uses in sonar systems, such as small target identification and fine imaging. A lower metal plate, a piezoelectric ceramic column array, and an upper metal plate make up the piezoelectric material structure. The resonance frequency of the bimetal plate type piezoelectric ceramic array material was theoretically analyzed, the sensitive element geometry of the bimetal plate piezoelectric material was studied using finite element simulation, the material choice for the metal plate affected the material’s performance, and the parameters of each size of the sensitive element were determined. These findings served as a guide for the production of a highly sensitive hydroacoustic transducer, as well as the structure of the piezoelectric material transducer in the water performance test. Results indicate that the structure transducer has a maximum transmit voltage response of 178 dB and a −3 dB bandwidth of approximately 25 kHz; the transducer may be efficiently tuned for a maximum receive sensitivity of −172 dB.

A combined cyclic viscoplasticity and entropy generation approach for modelling fatigue crack growth behavior of a nickel-based superalloy at high temperature

This paper developed an entropy-based approach to estimate the fatigue crack growth (FCG) behavior of GH4169 at high temperature. Firstly, systematic FCG experiments of GH4169 at 650 °C were conducted to obtain the experimental data. Then, a finite element analysis combined with Chaboche viscoplasticity and node release technology was developed to obtain cyclic stress–strain responses and entropy generation fields. Based on thermodynamic analysis, numerous discrete material elements as sub-systems were divided along the crack growth direction in simplified representative 2D middle section. Fracture process zone (FPZ) was defined equal to the size of discrete material elements (ρ) which can be determined by distribution of normal stress perpendicular to the crack plane. Subsequently, an effective cyclic entropy generation related to ρ was defined as the crack driving force. Effective cumulative entropy generation was calculated when crack increment was ρ and the results indicated that effective cumulative entropy generation float within a certain range. The average of all effective cumulative entropy generation values was defined as crack growth entropy. Finally, an entropy-based FCG rate model was established to estimate the FCG behavior of GH4169, and a good correlation between experimental results and model predictions are achieved for all the high temperature FCG tests.

Microstructural and Corrosion Behavior of Thin Sheet of Stainless Steel-Grade Super Duplex 2507 by Gas Tungsten Arc Welding

<div>Super duplex stainless steel (SDSS) is a type of stainless steel made of chromium (Cr), nickel (Ni), and iron (Fe). In the present work, a 1.6 mm wide thin sheet of SDSS is joined using gas tungsten arc welding (GTAW). The ideal parameter for a bead-on-plate trial is found, and 0.216 kJ/mm of heat input is used for welding. As an outcome of the welding heating cycle and subsequent cooling, a microstructural study revealed coarse microstructure in the heat-affected zone and weld zone. The corrosion rate for welded joints is 9.3% higher than the base metal rate. Following the corrosion test, scanning electron microscope (SEM) analysis revealed that the welded joint’s oxide development generated a larger corrosive attack on the weld surface than the base metal surface. The percentages of chromium (12.5%) and molybdenum (24%) in the welded joints are less than those in the base metal of SDSS, as per energy dispersive X-ray (EDX) analysis. Corrosion modeling is done using the COMSOL Multiphysics software. Electrochemical corrosion modeling is used to determine the electrolyte potential (i.e., 0.09 V) and current density (i.e., 0.2 A/m<sup>2</sup> to 1.8 A/m<sup>2</sup>). An entire mesh model contains 6240 elements. The largest and smallest element sizes are 4 mm and 0.1 mm, respectively. The maximum element rate of growth is 1.2.</div>

A New Technique for Urban and Rural Settlement Boundary Extraction Based on Spectral–Topographic–Radar Polarization Features and Its Application in Xining, China

Highly accurate data on urban and rural settlement (URS) are essential for urban planning and decision-making in response to climate and environmental changes. This study developed an optimal random forest classification model for URSs based on spectral–topographic–radar polarization features using Landsat 8, NASA DEM, and Sentinel-1 SAR as the remote-sensing data sources. An optimal urban and rural settlement boundary (URSB) extraction technique based on morphological and pixel-level statistical methods was established to link discontinuous URSs and improve the accuracy of URSB extraction. An optimal random forest classification model for URSs was developed, as well as a technique to optimize URSB, using the Google Earth Engine (GEE) platform. The URSB of Xining, China, in 2020 was then extracted at a spatial resolution of 30 m, achieving an overall accuracy and Kappa coefficient of 96.21% and 0.92, respectively. Compared to using a single spectral feature, these corresponding metrics improved by 16.21% and 0.35, respectively. This research also demonstrated that the newly constructed Blue Roof Index (BRI), with enhanced blue roof features, is highly indicative of URSs and that the URSB was best extracted when the window size of the structural elements was 13 × 13. These results can be used to provide technical support for obtaining highly accurate information on URSs.

Crack nucleation in heterogeneous bars: h- and p-FEM of a phase field model

Failure initiation and subsequent crack trajectory in heterogeneous materials, such as functionally graded materials and bones, are yet insufficiently addressed. The AT1 phase field model (PFM) is investigated herein in a 1D setting, imposing challenges and opportunities when discretized by h- and p-finite element (FE) methods. We derive explicit PFM solutions to a heterogeneous bar in tension considering heterogeneous E(x) and GIc(x)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$G_{Ic}(x)$$\\end{document}, necessary for verification of the FE approximations. GIc(x)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$G_{Ic}(x)$$\\end{document} corrections accounting for the element size at the damage zone in h-FEMs are suggested to account for the peak stress underestimation. p-FEMs do not require any such corrections. We also derive and validate penalty coefficient for heterogeneous domains to enforce damage positivity and irreversibility via penalization. Numerical examples are provided, demonstrating that p-FEMs exhibit faster convergence rates comparing to classical h-FEMs. The new insights are encouraging towards p-FEM discretization in a 3D setting to allow an accurate prediction of failure initiation in human bones.

A metric-based adaptive mesh refinement criterion under constrain for solving elliptic problems on quad/octree grids

In this work we propose and investigate the performance of a metric-based refinement criteria for adaptive meshing used for improving the numerical solution of an elliptic problem. We show that in general, when solving elliptic equations such as the Poisson-Helmholtz equation, the minimization of the interpolation error often used as local refinement criteria does not always guarantee the minimization of the total numerical error. Numerical and theoretical arguments are given to unveil the critical role of the mesh compression – the size aspect ratio between the finest cell size and the mean cell size of an adapted mesh – to determine whether the estimated error is purely local meaning that the interpolation error is a good enough error model for the total error or if other, non-local, sources of error need to be accounted for. We show through particular examples that a slightly sub-optimal mesh in terms of interpolation error may significantly reduce the total error of a numerical solution, depending on the value of the compression ratio and not on the number of grid points. Based on this observation, we propose a new method to exclude the grids where non-local errors can control the accuracy of the solution. This is achieved by an automatic estimation of the optimal compression ratio, which is imposed as an additional constraint in the minimal element size during the mesh adaptation process. The method is tested on quadtree and octree grids, showing very satisfactory performances in reducing the total numerical error despite the additional constrain imposed.

3-D visual analysis and quantitative characterization on coal matrix based on Micro-CT

Coal matrix, the storage space for adsorbed gas, is highly related to gas recovery and gas outburst disasters. However, existing research predominantly relies on theoretical assumptions, lacking experimental validation. In this study, we employ X-ray computed tomography (CT) to obtain high-resolution (13.45–15.157 μm) internal structure images of various coal samples. The three-dimensional model of coal matrix elements is reconstructed using entropy threshold segmentation and watershed segmentation. A systematic quantification of the size, shape, and spatial arrangement of these matrix elements is then performed. The results reveal that coal matrix element equivalent diameters span 148.45–2800.35 μm, with volume-weighted average equivalent diameters ranging from 1320.82 to 1720.87 μm. Additionally, we explore particle size distribution (PSD) models to characterize size-quantity relationships and identify the Jaky model based on the logarithmic-exponential equation as the optimal PSD model. Observation of the matrix element models shows that large matrix elements typically exhibit flattened morphologies with rough surfaces, while smaller ones tend to be blocky with smooth surfaces. The dimensionless shape factor results suggest that, when scale effects are not considered, the shape of larger matrix elements is more conducive to matrix-fracture gas mass transfer. Based on this, a functional relationship between the matrix shape factor and size, considering changes in matrix morphology, is established. Moreover, we algorithmically determine that the average number of adjacent matrix elements per individual matrix element falls between 13 and 15. For matrix clusters with large volume fractions, the size ratio of the central matrix element to adjacent matrix elements primarily falls within the range of 0.4–0.6, approximately following a Weibull distribution. According to this feature, a method for building a matrix cluster model used for simulations is proposed.

Melting and Heat Transfer Characteristics of Phase Change Material: A Comparative Study on Wire Mesh Finned Structure and Other Fin Configurations

Abstract In the present paper, a two-dimensional transient numerical study has been performed to investigate the influence of different fin designs on the melting and heat transfer characteristics of a phase change material (PCM), i.e., Paraffin wax, filled in square enclosures equipped with fin structures. Five distinct fin designs were examined: single rectangular, double rectangular, double triangular, double angled, and wire mesh. It is worth noting that all these fin designs have the equal heat transfer area. An isothermal heat source of temperature 350 K is provided at the left wall of the square enclosure and the remaining walls are assumed to be adiabatic. Six parameters were evaluated to determine the best fin configurations: melting time, enhancement ratio (ER), time savings, energy stored, mean power, and Nusselt number. The results show that all the fin designs outperformed as compared to model 1 (no fin configuration). Among the finned configurations, model 2 had the poorest performance, taking 1314 s to complete the melting, while model 6 had the most efficient fin design, with a melting time reduced by 67.53% compared to model 1. Model 6 also had the highest ER and mean power, i.e., 70.43% and 199.51%, respectively and as the melting process continued, the Nusselt number decreased. In addition to the above, we optimized the element size of the wire-mesh fin design using RSM methodology. This optimized design decreases the melting period by 70.04%. Overall, present study provides a comprehensive analysis of different finned configurations for improving the melting performance of the PCM in square enclosures and found wire-mesh fin design most appropriate and promising.

Estimation of the in-plane ultimate stress of lamellar tissue as a function of bone mineral density and osteocyte lacunae porosity

Background and ObjectiveDetailed finite element models based on medical images (μ-CT) are commonly used to analyze the mechanical behavior of bone at microscale. In order to simulate the tissue failure onset, isotropic failure criteria of lamellar tissue are often used, despite its non-isotropic and heterogeneous nature. The main goal of the present work is to estimate the in-plane ultimate stress of lamellar bone, considering the influence of mineral content and the porosity due to the osteocyte lacunae concentration. MethodsTo this aim, a representative volume cell of lamellar tissue is modeled numerically, including: (1) non-isotropic elastic properties of tissue as a function of the bone mineral density and (2) explicit modeling of the osteocyte lacunae, considering the range of porosity content, size and orientation of ellipsoid-shaped lacunae. Firstly, the element size for the finite element models have been defined from a preliminary convergence analysis. Bounds on the ultimate stress of non-porous lamellar tissue are estimated for two values of bone mineral density, considering the results of tensile and compressive tests of wet osteons from the literature. Subsequently, the ultimate stress of lamellar tissue considering several values of micro-porosity is addressed. ResultsResults obtained in this work show that the strength of lamellar bone decreases exponentially with the increase of lacunae porosity concentration. Ultimate stress of non-porous tissue (p=0%) increases with high mineral content, reaching a value of S¯transc=355.40±39.80 MPa for compression in the transversal direction of the fiber bundles, being BMD=1.246g/cm3. The mean value for the longitudinal to transverse strength ratio evaluated for porosity p=0%,1% and 5% and a mineral content BMD=1.2g/cm3, is 2.47:1 for tension and 1.55:1 for compression. These values are in agreement with literature. ConclusionsOsteocyte lacunae act as stress concentrators, acting as potential stimulus for the bone regeneration process. A novel micromechanical model for the in-plane ultimate stress of lamellar tissue as a function of mineral content and lacunae concentration is presented. Additional considerations about the intralamellar shear stress evolution are also given.

Formation and significance of topologically close-packed Laves phases in refractory high-entropy alloys

Refractory high-entropy alloys (RHEAs) inherently have a high potential to form topologically close-packed (TCP) Laves phases. Phase engineering to control the amount and type of TCP-Laves phases plays an important role in physical and mechanical properties of RHEAs. The present investigation addressed key parameters to TCP-Laves phase formation in RHEAs by considering thermodynamic calculations and empirical relations. Two novel TiVCrZrCo and TiVCrZrFe RHEAs were designed and fabricated as alloy models with high melting points (>2000 K) and relatively low density (∼6.5 g.cm−3) with multiphase microstructure including BCC matrix together with C14 and C15 Laves phases. Existence of atoms such as V, Co, Cr and Fe together with Zr, Ti, Nb and Hf refractory elements promote TCP-Laves phase formation. Very negative values for mixing enthalpy of atom pairs (ΔHmix) and remarkable differences between atomic sizes of the constitutive elements (δ) are key factors in the Laves phase formation in multicomponent alloys. The parameters including ΔHmix and δ must be in the defined range of −20 ≤ ΔHmix ≤ −3 kJ.mol−1 and 4 <δ < 10 to appear Laves phase in the microstructures of HEAs. TCP-Laves phase formation led to increase and decrease in hardness and fracture toughness, respectively. The understanding of key factors in RHEAs design can lead to phase engineering and the fabrication alloys with a desirable amount of TCP-Laves phase, which are suited for hydrogen storage applications or for applications in harsh operating environments exposed to high temperatures, irradiations, wear and erosion.

Analysis method useful for calculating various interface stress intensity factors efficiently by using a proportional stress field of a single reference solution modeling

This paper has proposed an efficient analysis method to calculate interface stress intensity factors (SIFs) based on a proportional stress field of a reference problem whose exact solution is available. In the previous proportional methods, the same crack length and the same element size were applied to both reference and unknown problems so that the same FEM error can be produced. Therefore, when analyzing many unknown problems, the conventional method needs to analyze many reference problems at the same time. Since this approach is time-consuming, this paper considers how to calculate many crack lengths efficiently by using only one single reference solution modeling. For this purpose, several general relations of SIFs are derived for the unknown and the reference problems when both crack length and element size are different. To analyze many unknown problems accurately by using a single reference solution modeling, how to choose the most suitable element dimension of the reference model is clarified. The proposed method is especially useful for crack propagation analysis.

DEVELOPMENT OF A VORTEX CHAMBER MODEL TO ESTIMATE THE VALUES OF THE LIQUID AND GASEOUS PHASES

When preparing well products from gas fields at complex gas treatment plants, dry stripped gas and unstable condensate are separated. During the exploitation of a field, development indicators (flow rate, reservoir pressure, composition of extracted products, etc.) change. At the same time, the efficiency of the apparatus for separating the gas and liquid phases decreases, which is accompanied by an increase in the entrainment of the liquid phase from the gaseous phase and can reach about 30–40 % of the total volume. This affects economic efficiency indicators, since the volume of unstable condensate decreases, and, consequently, the profit from its sale.The article is devoted to studying the effectiveness of using vortex chambers to reduce the entrainment of the liquid phase from the gaseous phase during the preparation of well products. The results of the research analysis are presented, methods for assessing phase transitions in vortex chambers and the mathematical models used are described.The work defines the basic equations that describe the process of entrainment of the liquid phase from the gaseous phase. Thus, to describe the phase boundaries, the volume of fluid model is used, and the standard k – ε is used as a flow turbulence model. A model of a vortex chamber is being developed in ANSYS FLUENT, the mesh model is built using elements – tetrahedrons. The size of the mesh elements was determined from the condition of stability of the resulting solution and amounted to 10 mm.The influence of chamber geometry, liquid mass velocity and other parameters on the process of phase entrainment is considered.A comparison of the results of experimental studies and modeling results shows the possibility of reducing the entrainment of the liquid phase as a result of the use of vortex chambers and justify the selected parameters for developing the model.

Comparison of textured nylon surfaces manufactured by CNC micromachining and 3D printing

Patterned surfaces can be designed and engineered to control friction and wear resistance in various applications. In the present work, a comparative analysis of two contrasting manufacturing processes for surface texturing of polyamide, namely 3D printing by fused filament fabrication (FFF) and CNC micromachining (micro CNC) as additive and subtractive manufacturing techniques respectively was carried out. The analysis included the evaluation of the topography, roughness, mechanical properties, and tribological performance of the texturized surfaces. Both manufacturing routes were suitable for fabricating deterministic surfaces, being the micro CNC method the one that provided better repeatability and surface finishing as well as higher hardness. These factors influenced the tribological behavior of the polyamide when in contact with AISI 304 stainless steel. The effects of size, distribution and height of the texture elements and the manufacturing technique on the Coefficient of Friction (COF) were also discussed.

HEAT TRANSFER ANALYSIS ON BASE OF HEAT SINK WITH UNIFORM HEAT GENERATION

The present work describes the important results of problem of combined conduction, convection and radiation from a rectangular heat sink base subjected to uniform internal heat generation in its base. The governing equations for temperature distribution within the base are obtained by a relevant energy balance between rate of heat generation, heat conducted, convected and radiated. The derived non-linear partial differential equations are converted into algebraic form using finite difference method with 2nd order accuracy. The obtained algebraic equations have been solved simultaneously by Gauss-siedel iterative method with stringent convergence criteria and a C code has been written for above purpose. The computational domain has been discretized uniformly in all directions. A grid convergence test has been carried out to optimize the size of element. A complete parametric study has been carried out to study the effect of material of heat sink, regime of convection and role of radiation on local temperature distribution and maximum heat sink base temperature. It has been observed that material with high conductivity decreases the local temperature gradient and thus the maximum temperature. The regime of convection has a great influence on maximum temperature of heat sink with 15K of drop in Tmax. for a change in flow regime from free convection to forced convection. A special study on exclusive role of convection on maximum temperature reveals that one should not ignore the role forced convection in heat dissipation. Key Words: Heat generation, Multimode heat transfer, Radiation