Multiphysics simulation of microwave curing in micro-electronics packaging applications

Purpose – This paper aims to present an open-ended microwave curing system for microelectronics components and a numerical analysis framework for virtual testing and prototyping of the system, enabling design of physical prototypes to be optimized, expediting the development process. Design/methodology/approach – An open-ended microwave oven system able to enhance the cure process for thermosetting polymer materials utilised in microelectronics applications is presented. The system is designed to be mounted on a precision placement machine enabling curing of individual components on a circuit board. The design of the system allows the heating pattern and heating rate to be carefully controlled optimising cure rate and cure quality. A multi-physics analysis approach has been adopted to form a numerical model capable of capturing the complex coupling that exists between physical processes. Electromagnetic analysis has been performed using a Yee finite-difference time-domain scheme, while an unstructured finite volume method has been utilized to perform thermophysical analysis. The two solvers are coupled using a sampling-based cross-mapping algorithm. Findings – The numerical results obtained demonstrate that the numerical model is able to obtain solutions for distribution of temperature, rate of cure, degree of cure and thermally induced stresses within an idealised polymer load heated by the proposed microwave system. Research limitations/implications – The work is limited by the absence of experimentally derived material property data and comparative experimental results. However, the model demonstrates that the proposed microwave system would seem to be a feasible method of expediting the cure rate of polymer materials. Originality/value – The findings of this paper will help to provide an understanding of the behaviour of thermosetting polymer materials during microwave cure processing.

Coupled FDTD-FVM simulation of microwave heating of polymer materials for micro-systems packaging applications

Heating in an idealised polymer load in a novel open-ended variable frequency microwave oven is numerically simulated using a couple solver approach. The frequency-agile microwave oven bonding system (FAMOBS)is developed to meet rapid polymer curing requirements in microelectronics and optoelectronics manufacturing. The heating of and idealised polymer load has been investigated through numerical modelling. Assessment of the system comprises of simulation of electromagnetic fields and of temperature distribution within the load. Initial simulation results are presented and contrasted with experimental analysis of field distribution

Numerical simulation of microwave thawing using a coupled solver approach

Thawing of a frozen food product in a domestic microwave oven is numerically simulated using a coupled solver approach. The approach consists of a dedicated electromagnetic FDTD solver and a closely coupled UFVM multi-physics package. Two overlapping numerical meshes are defined; the food material and container were meshed for heat transfer and phase change solution, whilst the microwave oven cavity and waveguide were meshed for the microwave irradiation. The two solution domains were linked using a cross-mapping routine. This approach allowed the rotation of the food load to be captured. Power densities obtained on the structured FDTD mesh were interpolated onto the UFVM mesh for each timestep/turntable position. The UFVM solver utilised the power density data to advance the temperature and phase distribution solution. The temperature-dependant dielectric and thermo-physical properties of the food load were updated prior to revising the electromagnetic solution. Changes in thermal/electric properties associated with the phase transition were fully accounted for as well as heat losses from product to cavity. Two scenarios were investigated: a centric and eccentric placement on the turntable. Developing temperature fields predicted by the numerical solution are validated against experimentally obtained data. Presented results indicate the feasibility of fully coupled simulations of the microwave heating of a frozen product. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Numerical simulation of encapsulant curing within a variable frequency microwave processing system

Curing of encapsulant material in a simplified microelectronics package using an open oven Variable Frequency Microwave (VFM) system is numerically simulated using a coupled solver approach. A numerical framework capable of simulating electromagnetic field distribution within the oven system, plus heat transfer, cure rate, degree of cure and thermally induced stresses within the encapsulant material is presented. The discrete physical processes have been integrated into a fully coupled solution, enabling usefully accurate results to be generated. Numerical results showing the heating and curing of the encapsulant material have been obtained and are presented in this contribution. The requirement to capture inter-process coupling and the variation in dielectric and thermophysical material properties is discussed and illustrated with simulation results.

Multiphysics simulation of microwave processing using a multidomain coupled solver approach

Microwave processing of materials is numerically simulated using a coupled solver approach. Microwave heating is a complex coupled process due to the variation in dielectric properties during heating. The effects of heating an object in a electromagnetic field directly influence the manner in which it interacts with the field. Simplifying assumptions and empirical solutions do not capture the fundamental physics involved and, in general, do not provide usefully accurate solutions in a number of practical problems. In order to capture the underlying processes involved in microwave heating, the problem must be looked at in a holistic manner rather than a number of discrete processes. This contribution outlines a coupled-solver multiphysics analysis approach to the solution of practical microwave heating problems.

Computer simulation of crack propagation in power electronics module solder joints

A numerical modelling method for the analysis of solder joint damage and crack propagation has been described in this paper. The method is based on the disturbed state concept. Under cyclic thermal-mechanical loading conditions, the level of damage that occurs in solder joints is assumed to be a simple monotonic scalar function of the accumulated equivalent plastic strain. The increase of damage leads to crack initiation and propagation. By tracking the evolution of the damage level in solder joints, crack propagation path and rate can be simulated using Finite Element Analysis method. The discussions are focused on issues in the implementation of the method. The technique of speeding up the simulation and the mesh dependency issues are analysed. As an example of the application of this method, crack propagation in solder joints of power electronics modules under cyclic thermal-mechanical loading conditions has been analyzed and the predicted cracked area size after 3000 loading cycles is consistent with experimental results.

Computer simulation of aluminium wirebonds with globtop in power electronics modules

In this paper, computer modelling techniques are used to analyse the effects of globtops on the reliability of aluminium wirebonds in power electronics modules under cyclic thermal-mechanical loading conditions. The sensitivity of the wirehond reliability to the changes of the geometric and the material property parameters of wirebond globtop are evaluated and the optimal combination of the Young's modulus and the coefficient of thermal expansion have been predicted.

Numerical simulation of heat transfer in rectangular microchannel

Numerical simulation of heat transfer in a high aspect ratio rectangular microchannel with heat sinks has been conducted, similar to an experimental study. Three channel heights measuring 0.3 mm, 0.6mmand 1mmare considered and the Reynolds number varies from 300 to 2360, based on the hydraulic diameter. Simulation starts with the validation study on the Nusselt number and the Poiseuille number variations along the channel streamwise direction. It is found that the predicted Nusselt number has shown very good agreement with the theoretical estimation, but some discrepancies are noted in the Poiseuille number comparison. This observation however is in consistent with conclusions made by other researchers for the same flow problem. Simulation continues on the evaluation of heat transfer characteristics, namely the friction factor and the thermal resistance. It is found that noticeable scaling effect happens at small channel height of 0.3 mm and the predicted friction factor agrees fairly well with an experimental based correlation. Present simulation further reveals that the thermal resistance is low at small channel height, indicating that the heat transfer performance can be enhanced with the decrease of the channel height.

Three-dimensional simulation of the influence of convection on dendritic solidification

Bulk and interdendritic flow during solidification alters the microstructure development, potentially leading to the formation of defects. In this paper, a 3D numerical model is presented for the simulation of dendritic growth in the presence of fluid flow in both liquid and semi-solid zones during solidification. The dendritic growth was solved by the combination of a stochastic nucleation approach with a finite difference solution of the solute diffusion equation and. a projection method solution of the Navier-Stokes equations. The technique was applied first to simulate the growth of a single dendrite in 2D and 3D in an isothermal environment with forced fluid flow. Significant differences were found in the evolution of dendritic morphology when comparing the 2D and 3D results. In 3D the upstream arm has a faster growth velocity due to easier flow around the perpendicular arms. This also promotes secondary arm formation on the upstream arm. The effect of fluid flow on columnar dendritic growth and micro-segregation in constrained solidification conditions is then simulated. For constrained growth, 2D simulations lead to even greater inaccuracies as compared to 3D.