Thermal Field Simulation Research Articles

Multiphysics Modeling of Low-Intensity Pulsed Ultrasound Induced Chemotherapeutic Drug Release from the Surface of Gold Nanoparticles.

Currently, no numerical model for low-intensity pulsed ultrasound (LIPUS)-triggered anticancer drug release from gold nanoparticle (GNP) drug carriers exists in the literature. In this work, LIPUS-induced doxorubicin (DOX) release from GNPs was achieved in an ex vivo tissue model. Transmission electronic microscopy (TEM) imaging was performed before and after LIPUS exposure, and significant aggregation of the GNPs was observed upon DOX release. Subsequently, GNP surface potential was determined before and after LIPUS-induced DOX release, using a Zetasizer. A numerical model was then created to predict GNP aggregation, and the subsequent DOX release, via combining a thermal field simulation by solving the bioheat transfer equation (in COMSOL) and the Derjaguin, Landau, Verwey, and Overbeek (DLVO) total interaction potential (in MATLAB). The DLVO model was applied to the colloidal DOX-loaded GNPs by summing the attractive van der Waals and electrostatic repulsion interaction potentials for any given GNP pair. DLVO total interaction potential was found before and after LIPUS exposure, and an energy barrier for aggregation was determined. The DLVO interaction potential peak amplitude was found to drop from 1.36 kBT to 0.24 kBT after LIPUS exposure, translating to an 82.4% decrease in peak amplitude value. It was concluded that the interaction potential energy threshold for GNP aggregation (and, as a result, DOX release) was equal to 0.24 kBT.

Thermal Field Simulation and Analysis of Tulip Contact Based on Virtual Material Method

In order to analyze the temperature distribution of circuit breaker tulip contact, a thermal field simulation model for tulip contact is built combined with the electromagnetic-stress-thermal coupled finite element simulation and virtual material method which models the contact resistance characteristics of tulip contact. The temperature rise characteristic of tulip contacts under rated current conditions and short-circuit conditions is analyzed with the model. Results show that the virtual material method can describe the contraction resistance and additional resistance at the end of the contact finger simultaneously, so as to simulate the partial high temperature of the contact interface. Then, the influence of spring pressure and pressure uniformity on the thermal field of tulip contact is analyzed. Results show that the decrease of spring pressure and the imbalance of contact finger pressure may lead to an increase in the temperature rise of the tulip contact. The proposed method can be used as a reference for the design and the overheating analysis of tulip contact.

A theoretical approach to the residual stress assessment based on thermal field evaluation in laser beam welding

Residual stresses are one of the major issues in welded parts, since they could be detrimental to the integrity of components and structure. Their determination is rather complex and could be an arduous task, both when it is based on experimental methods and on numerical simulations. The proposed work presents a theoretical approach to the prediction of the longitudinal residual stress distribution, based on a parameterized multi-source model for thermal field simulation in laser welding previously introduced. Reference is made to the case of “keyhole” full penetration welding mode obtained by CO2 laser beam single pass on butt-positioned AISI 304L plates. The resolution of the thermal field allows the analytical calculation of the distribution of the longitudinal residual stresses in two ways: one makes use of a simplified formulation of the distribution well-known in the literature; a second modality makes use of a procedure for residual stress generation, which is based on a combined processing of thermal profiles and the corresponding heating–cooling cycles calculated in single points as their distance from the welding axis varies, and provides a complete characterization of the distribution of longitudinal residual tensile stresses. After the introduction of thermal field modeling, both the proposed residual stress calculation procedures are detailed, applied to the analyzed case, and validated, highlighting the differences in the approaches and results.

Control of flow, thermal and pollutant concentration fields by entrainer air streams to improve fresh air quality intake into a semiconductor manufacture/ processing plant

Simulations of turbulent flow, thermal and concentration fields around a semiconductor manufacture/processing plant are reported. Monitored pollutants are ammonia (NH 3 ), hydrofluoric acid (HF), total volatile organic compounds (TVOC) and PM 2.5 that are released from chimneys on the roofs of various buildings. Some of the fresh air for the plant is supplied by the make-up air units (MAU) located on the roofs of the buildings. Returned pollutants through the MAUs can end up in the cleanroom thus decreasing the quality of wafer products. A 3-D transient large eddy simulation model was developed to predict the environment near the plant and transport of pollutants during the summer and winter seasons. Concentrations of pollutants at the inlet planes of the MAUs exhibited marked variations during the day. Pollutants with concentrations higher than industry thresholds were generally re-directed back into the plant through the six MAUs except two that were free of HF and TVOC. Modified flow, thermal and pollutant concentration fields in the vicinity of the plant brought about by introducing 12 vertically-oriented entrainer air streams (EAS) into the environment were also studied in order to quantify their effects on reducing the concentrations of pollutants at the inlet planes of the MAUs below accepted industry thresholds. It is shown that dispersion of the pollutants was markedly affected by the introduction of these novel entrainer air streams. Ingestion of pollutants into the plant was completely eliminated during (i) summer and (ii) winter except for MAU5 (NH 3 ) and MAU6 (NH 3 and TVOC) that exhibited noticeable reduction of concentrations by at least 73.3%. The proposed system can lead to cleaner production practices through eliminating product rejection, lower utilization of raw materials, reduced introduction of pollutants in the environment, etc.

Mechanical behavior of thermally-induced shape memory polymer composites thin-wall tube in non-uniform thermal field: Experiments and simulations

Thermally-induced shape memory polymer composites (SMPCs) have the capability to self-deploy when stimulated by a specific temperature. As a basic component of deployable truss structure, SMPCs thin-wall circular tubes will undergo complex non-linear deformation and changing thermal boundary conditions during application. This article aims to research the mechanical behavior of thermally-induced SMPCs thin-wall circular tubes in a non-uniform thermal field. Thermal properties of SMPCs including thermal conductivities and coefficients of thermal expansion are tested as fundamental temperature-dependent parameters. A 500 mm-long SMPCs thin-wall tube is conducted in the flattening-folding-cooling-reheating experiments through a specially manufactured loading apparatus. Combining phase transition theory and Fourier heat conduction, a coupling algorithm is proposed to describe the thermo-mechanical behavior of SMPCs, which is implemented by the user subroutines VUMAT by commercial finite element analysis packages ABAQUS. By comparing the results from the simulation of the non-uniform thermal field, the simulation of the uniform thermal field and experiments, the coupling algorithm is verified, and the effect on SMPCs from the non-uniform thermal field is given.

Simulation of Electro-Thermal Runaway and Thermal Limits of a Loaded HVDC Cable

This paper presents simultaneous simulation of interdependent electric and thermal fields inside a dc cable under loaded conditions using analogous distributed circuit models for thermal and electrical phenomena. Certain nonlinearities are integrated in such circuit models, leading to simulation of electro-thermal runaway, for the first time, in a dc cable. Using the suggested, nonlinear circuit models, complete dynamics of the temperature profile and electric field inside the cable dielectric can also be obtained under various test and operating conditions until electro-thermal runaway limits in time domain. The results are useful for design and development of dc cables with safety factors, apart from understanding the electro-thermal limits. The models are validated with experimental results, found to be in excellent agreement.

Thermal Field and Stress Analysis of Induction Motor with Stator Inter-Turn Fault

Inter-turn fault (ITF), a typical motor fault, results in significant variations in the thermal characteristics of a motor. For fault, temperature rise (TR) experiments and thermal field-stress field simulations of an induction motor are carried out to reveal the fault characteristics related to ITF. First, based on the actual structure and the cooling type of the motor, a whole-domain simulation model of the fault thermal field was established. The reasonable equivalence of the motor and the calculation of the heat transfer boundaries were conducted during the modeling process. Then, the three-dimensional transient thermal field under a rated load before and after the fault was obtained, and the accuracy of the simulation could be validated through the comparison of the measured TR at several temperature-measuring points. The heat-transfer law and the notable thermal characteristics of the fault can be presented by analyzing the simulated and measured temperature data. In addition, a fault feature is proposed to provide a reference for diagnosis using the temperature difference of winding at different positions at different moments. Finally, the rotor thermal stress distribution of the normal and faulty motor is obtained by thermal-stress-coupled calculation, which can be used to evaluate the possibility of rotor fault caused by ITF.

Simulation of Thermal and Electric Fields in Electrodes with Allowance for Melting and Fritting of Films

Simulation of Thermal and Electric Fields in Electrodes with Allowance for Melting and Fritting of Films

Simulation of Thermal Field in Eutectic Microwave Bonding for Electrical Connection of Photovoltaic Cells

Sustainable development requires green energy and low carbon footprint in manufacturing sector of photovoltaic systems. The electrical connections of photovoltaic cells need to have low electrical resistance in order to reduce the electrical losses and therefore to improve the performance of the photovoltaic panels. This paper aims to present researches related to bonding of wires that connect solar cells by using microwave technology. The microwave bonding has the main advantage that offers fast bonding but, in the same time, this technology does not offer stability of the thermal heating. Two different unwanted phenomena like thermal runaway and plasma arc discharge often lead to the damaging of copper and aluminum wires used in electrical connection. The study presented in this paper is focused on simulation of the thermal field developed in copper wires in order to optimize the bonding process and increase the quality of products. The simulation of the thermal field has been done using Fourier equations for conducting heating in copper materials and eutectic alloys. The simulation model has been validated through experimental heating using a 6 kW water-cooled microwave generator controlled by a matching load auto-tuner for best transfer of the power from generator to copper wires. The temperature has been measured in real time using an infrared pyrometer for metals with 2.3 μm spectral range and measurement range between 0o C and 7000 C. The study is finalized with elaboration of mathematical model for microwave-injected power as function for temperature developed in copper wires that can be applied with success in further microwave bonding applications of copper wires. In addition, the electrical resistance of bonded wires was measured in order to collect feedback for improving the microwave bonding process.

CAD Approach to Control High-Temperature Processes in SiC Technology

Manufacturing silicon carbide semiconductor devices may require high-temperature treatment in closed graphite reactors. This makes temperature control of processed SiC substrates difficult, since their temperature cannot be measured directly. As the monitoring of the SiC sample temperature is critically important for proper process flow, an indirect method involving the use of the CAD approach has been developed. A numerical model of a furnace reactor was created on the basis of the commercial ANSYS package, allowing for the simulation of thermal fields under given heat-dissipation conditions in the modeled area and in the presence of gaseous and liquid media participating in heat exchange and transport. Obtained simulation results remain very consistent with the reference temperature measurements of selected areas of the reactor. The model acts as an accurate tool for temperature distribution verification during the high-temperature annealing of and diffusion of dopants for silicon carbide.

Growth of large size near-stoichiometric lithium niobate single crystals with low coercive field for manufacturing high quality periodically poled lithium niobate

Periodically poled lithium niobate is an important component in mid-infrared lasers due to its perfect nonlinear optical property. However, fabrication of high-power laser based on periodically poled lithium niobate is still a great challenge, because lithium niobate crystal wafer with large thickness for large clear aperture periodically poled lithium niobate is hard to be completely poled due to the high coercive field of commercial congruent lithium niobate (LiNbO3). In this work, we have successfully grown 3-inch near-stoichiometric LiNbO3 by efficiently regulation of temperature gradient with the aid of thermal field simulation. Characteristics on the as-grown large size near-stoichiometric LiNbO3 proved that the [Li]/[Nb] of crystal is around 1:1, and the crystal possesses a typical single domain with high optical homogeneity. Compared with congruent and Mg-doped LiNbO3, the near-stoichiometric LiNbO3 crystal has better piezoelectric performance. The forward and backward coercive fields of the near-stoichiometric LiNbO3 crystal were 2.90 kV/mm and 2.43 kV/mm, respectively, which were approximately 8 times lower than that for CLN (21 kV/mm). The periodically polarization experiments proved that only 1.36 kV can realize the complete polarization for a 1 mm thick near-stoichiometric LiNbO3 wafer. The mass-productive near-stoichiometric LiNbO3 crystals grown in this work will have great applications in high-power mid-infrared solid-state laser.

Multigrid method for numerical modelling of high temperature superconductors

An approach to numerical simulation of three-dimensional electrical and thermal fields in high-temperature superconductors is described. In such a semiconductor, the phenomena of superconductivity are observed at high temperatures above the temperature of liquid nitrogen. The absence of a generally accepted theory of superconductivity leads to the need to study physical processes in semiconductor structures using mathematical simulations. The main attention is paid to the calculation of temperature and electric current distributions in large-size mesas with a self-heating effect. An efficient algorithm for solving the equations describing these distributions is constructed. The basis of the algorithm is an adaptive multigrid method on structured Cartesian grids. The adaptability is based on the Chebyshev iterative method for constructing the smoothing procedures at each grid level and for solving the coarsest grid equations. The adaptive technique allows us to realistically simulate the anisotropic phenomena. The functionality of the algorithm is demonstrated along with an example of solving an anisotropic model problem with discontinuous coefficients.

Directly fabricated Al2O3/GdAlO3 eutectic ceramic with large smooth surface by selective laser melting: Rapid solidification behavior and thermal field simulation

Selective laser melting (SLM), a novel approach for one-step melting and solidifying ceramic powder beds layer by layer without post-process of degreasing and sintering, has been developed to directly prepare highly dense (>95 %) Al2O3/GdAlO3(GAP) eutectic composite ceramics with large smooth surfaces. Compact net-shaped plates with the maximum size of 73 × 24 × 5 mm3 are obtained by different strategies of laser pre-heating and multi-tracks’ deposition without any binders. Combined with the finite element thermodynamic coupling simulation results, it is proved that the stress between the substrate and depositions during SLM can be greatly reduced by the step-up preheating, and thus effectively improving the ceramic forming quality. The macro-morphology, microstructure evolution, rapid solidification behavior and mechanical properties of the SLM-ed eutectic ceramics are systematically investigated at different laser processing parameters. The microstructure transforms from ultra-fine irregular eutectic to complex regular eutectic with the increase of the scanning rate. The average eutectic spacing, and solidification rate has an approximately linear relationship consistent with the Jackson-Hunt (JH) model. The microhardness and fracture toughness can reach 17.1 ± 0.2 GPa and 4.5 ± 0.1 MPa·m1/2, respectively. The results indicate that SLM method is a highly effective technique for fabricating high-performance net-shaped structural composite ceramics.

Investigation on Laser Beam Weldability of AISI 304L Plates Based on Thermal Field Simulation by Experimentally-Fitted Analytical Modeling

The objective of the work is to present a new approach to the simulation of the thermal field in laser beam welding, based on an experimentally-fitted analytical model, applied to investigate the weldability of AISI 304L austenitic steel. Reference is made to the welding trial in a single pass of two 10 mm thick butt-positioned plates. Welding was performed under the keyhole full penetration mode, which is characteristic of high-power laser beam, and simulated by an analytical model based on a multipoint-line thermal source system and fitted on the experimental fusion zone profile. The model was applied to simulate the effects of welding speed changes on thermal fields and cooling rates, in order to determine how they can affect the weld composition, the solidification mode and the possible formation of a sensitized zone in the heat affected zone. A limit value of welding speed, which allows the weld formation without lack of fusion, was identified. For all the welding speeds considered, the formation of a sensitized zone can be excluded. The contribution of welding speed on cooling rate, not significant near the welding axis, results to be determinant at the boundary of the fused zone with base metal. The combined choice of the filler material and the welding speed, which in all cases gives rise to primary ferrite solidification modes, affects the content of residual ferrite, which must be balanced to enhance the resistance to solidification cracking, avoiding the adverse effects due to too high contents. As a conclusion, the model proves to be a valid support in investigating the thermal effects, which result from the setup of welding parameters, on the weldability of the base metal-filler system.Graphic

Thermal field distribution investigation and simulation of silver paste heating fabric by screen printing based on Joule heating effect

Thermal field distribution investigation and simulation of silver paste heating fabric by screen printing based on Joule heating effect

Повышение эффективности поверхностно-термического упрочнения деталей машин в условиях совмещения обрабатывающих технологий, интегрируемых на единой станочной базе

Introduction. In the manufacturing industry, there is a particular interest in the development of a new type of technological equipment, which makes it possible to implement methods for modifying the parts surface layers by processing it with concentrated energy sources. The combination of two processing technologies (mechanical and surface-thermal operations) in the conditions of integrated equipment makes it possible to neutralize the disadvantages of monotechnologies and obtain new effects that are unattainable when using technologies separately. The use of hybrid machine tools in conjunction with the developed technological recommendations will allow achieving a multiple increase in the technical and economic efficiency of production, resource and energy saving, which in turn will contribute to an increase in the competitiveness of products and the renewal of the technological paradigm. Purpose of work: increasing productivity and reducing energy consumption during surface-thermal hardening of machine parts by exposure to concentrated energy sources under conditions of integrated processing. Theory and methods: studies of the possible structural composition and layout of hybrid equipment during the integration of mechanical and surface-thermal processes are carried out taking into account the main provisions of structural synthesis and the components of metalworking systems. Theoretical studies are carried out using the basic provisions of system analysis, geometric theory of surface formation, design of metalworking machines, finite-element method, mathematical and computer simulation. Mathematical simulation of thermal fields and structural-phase transformations in the case of HEH HFC is carried out in the ANSYS and SYSWELD software packages, using numerical methods for solving the differential equations of unsteady thermal conductivity (Fourier's equation), carbon diffusion (Fick's second law), and elastoplastic behavior of the material. The verification of the simulation results is carried out by conducting field experiments using: optical and scanning microscopy; mechanical and X-ray methods for determining residual stresses. In the study, Uone JD520 and Form Talysurf Series 2 profilograph-profilometers are used to simultaneously measure shape deviations, waviness and surface roughness. Surface topography is assessed using a Zygo New View 7300 laser profilograph-profilometer. The microhardness of the hardened surface layer of parts is evaluated on a Wolpert Group 402MVD device. Results and discussion. An original method of structural-kinematic analysis for pre-design research of hybrid metalworking equipment is presented. Methodological recommendations are developed for the modernization of metal-cutting machine tools, the implementation of which will allow the implementation of high-energy heating by high-frequency currents (HEH HFC) on a standard machine-tool system and ensure the formation of high-tech technological equipment with expanded functionality. A unified integral parameter of the temperature-time effect on a structural material is proposed when the modes of hardening by concentrated heating sources are assigned, which guarantee the required set of quality indicators of the surface layer of machine parts, while ensuring energy efficiency and processing productivity in general. It is experimentally confirmed that the introduction into production of the proposed hybrid machine tool in conjunction with the developed recommendations for the purpose of the HEH HFC modes in the conditions of integral processing of a “Plunger bushing” type part in relation to the factory technology allows increasing the productivity of surface hardening by 3.5 ... 4.1 times, and reduce energy consumption by 9.5 ... 11.3 times.

Numerical Simulation of Thermal Field in Mass Concrete With Pipe Water Cooling

Water pipe cooling is mainly used to control temperature in the construction of mass concrete structures. It is important to reveal how to accurately stimulate the temperature field of mass concrete under action of this water pipe cooling. This paper presents a new method for this purpose. In this method, the contact surface of the water pipe and the concrete is used as the heat dissipation surface into the control equation and the composite Multiquadrics radial basis function (MQ-RBF) and low-order linear polynomial combination are used to discrete the spatial domain. The heat dissipation surface of the water pipe is included in the boundary conditions so that there is no need to build the refined water pipe modeling. This new method not only reduces the calculation cost but also ensures calculation accuracy. Through four calculation examples, this paper show that the algorithm has advantages in the numerical simulation of the concrete temperature field with water pipe cooling.

Simulation and experiment investigaton of electromagnetic and thermal field during induction cladding

A fully electrical–magnetic–thermal coupled model is proposed to investigate the heating mechanism of induction cladding. The simulation results reveal that the coating and the thin layer close to the substrate surface were heated layer by layer. Furthermore, the non-uniform electromagnetic field generated an uneven temperature distribution, which had a larger impact in the longitudinal direction than in the cross section. The maximum relative errors between the simulated and measured temperatures indiacate that the numerical calculation fit the experimental results well. The impact of the non-uniform temperature distribution on the coating properties was experimentally studied. The specimen in the middle position had a more complicated phase composition, lower porosity, and higher microhardness than those at the two ends.

An LBM-based investigation of thermal buoyancy and arrangement angle on flow characteristics and heat transfer over four heated square cylinders

In this work, the effects of thermal buoyancy and rotation angle on the hydrodynamics and mixed convection heat transfer characteristics over four heated identical square cylinders placed in a parallel horizontal channel are numerically studied. For the flow field, the Lattice Boltzmann method (LBM) with multiple relaxation time collision operator (MRT) is implemented and the Bhatnagar-Gross-Ktook collision operator (BGK) is adopted for the thermal field simulation. The distance between two adjacent square cylinders is fixed for a parametric study and a uniform velocity is imposed at the inlet region. Ranges of the influence parameters are defined as: α = 45°, 90°, 120°; 50 ≤ Re ≤ 150; 0 ≤ Ri ≤ 0.5. Several flow and heat transfer parameters including the time-averaged lift/drag coefficients, the global time-averaged Nusselt number and the distributions of local Nusselt number are analyzed, respectively. The vorticity and isotherm patterns are also plotted to illustrate the visual disturbance of the thermal buoyancy and the rotation angle on the flow and thermal fields. Numerical results show that the flow and thermal patterns for all arrangements become disordered and nonsystematic with the introduction of thermal buoyancy, and the thermal pattern for the downstream cylinders has an upward drift in the wake region. The time-averaged lift coefficient for each square cylinder is relatively sensitive to the influence parameters and the arrangement, and it is always less than zero when the Richardson number is not equal to zero. The global time-averaged Nusselt number is directly proportional to the Reynolds number and Richardson number, and the local Nusselt numbers on the front and rear faces for each square cylinder reach a maximum and minimum, respectively. It is also observed that the distribution of local Nusselt number for all rotation angles is intensively affected by the thermal buoyancy.

Characterization and simulation of the nonlinear thermal field of the aramid/bismaleimide composites caused by the dielectric heating effects of the microwave radiations

AbstractThe aramid/bismaleimide composites are severely affected by the dielectric heating effects of the microwave radiations during use as the microwave functional materials. This study aims to characterize and predict the nonlinear thermal fields and thermal damages of the aramid/bismaleimide composites caused by the microwave radiations. The microwave radiation experiments were performed using a tailor‐made system, and the equilibrium temperature and thermal damage threshold of the aramid/bismaleimide composites were experimentally determined. The temperature variation, thermal damage mechanism, and effect of the composite thickness associated with the dielectric heating were analyzed. Further, a simulation approach was developed to predict the nonlinear thermal fields of the composites caused by the microwave radiations. The dielectric heating process was considered an internal heat source, and a loop program was subsequently designed to implement the nonlinear thermal load. The developed simulation method efficiently predicted the equilibrium temperature and thermal damage of the composite materials due to the dielectric heating effect. The microwave radiation power was suggested to remain below 550 W for 3 mm thick aramid composites during service as a wave‐transmitting material.