NUMERICAL SIMULATION AND MOLDABILITY INVESTIGATION OF MICRO-FEATURES

Similar documents
Dynamic Behavior of Core-material Penetration in Multi-Cavity Co-Injection Molding

OPTICAL PROPERTY INVESTIGATION IN SEQUENTIAL MULTI-COMPONENT MOLDING

A Numerical Study of an Injection-Compression Molding Process by using a Moving Grid Bambang Arip Dwiyantoro

INVESTIGATION ON WARPAGE AND ITS BEHAVIOR IN SEQUENTIAL OVERMOLDING

TEMPERATURE INVESTIGATION IN REAL HOT RUNNER SYSTEM USING A TRUE 3D NUMERICAL METHOD

Multiscale simulations of the injection moulding of a diffractive optical element

International Communications in Heat and Mass Transfer 36 (2009) Contents lists available at ScienceDirect

reported that the available simple contact conductance model was expressed as [5][6]: h sum = h solid + h fluid (1) Where h sum, h solid and h fluid a

The effect of Entry Region on Thermal Field

Nicholas Cox, Pawel Drapala, and Bruce F. Finlayson Department of Chemical Engineering, University of Washington, Seattle, WA, USA.

3D Compression Molding

On Clean Cooling Systems for Wind Turbine Nacelle operating in Hot Climate

Computer Simulation of the Filling Process in Gas- Assisted Injection Molding Based on Gas-Penetration Modeling

CASE MF Report of pipe connector

2009 Best Paper Understanding and Quantification of Die Drool Phenomenon During Polypropylene Extrusion Process

Drag reduction in a class 8 truck - scaled down model

Lattice Boltzmann Simulation of One Particle Migrating in a Pulsating Flow in Microvessel

CFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate

COMPUTATIONAL MODELING OF POLYMER FLOW IN MICROCAVITIES THROUGH A MICROSCREEN

Analysis of Mixing Chambers for the Processing of Two-Component Adhesives for Transport Applications

Natural Convection from Horizontal Rectangular Fin Arrays within Perforated Chassis

Differential relations for fluid flow

Computational Fluid Dynamics 2

Rheometry. II.1 Introduction

Viscosity and Polymer Melt Flow. Rheology-Processing / Chapter 2 1

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER

Chapter 1 Introduction

Micro Cooling of SQUID Sensor

Getting started: CFD notation

MHD and Thermal Dispersion-Radiation Effects on Non-Newtonian Fluid Saturated Non-Darcy Mixed Convective Flow with Melting Effect

Direct numerical simulation database for supercritical carbon dioxide

Impacts of Electroosmosis Forces on Surface-Tension- Driven Micro-Pumps

Numerical modeling of magnetic induction and heating in injection molding tools

Lecture 3. Properties of Fluids 11/01/2017. There are thermodynamic properties of fluids like:

Study on the flow and thermal characteristics of a heat storage system

A multiscale framework for lubrication analysis of bearings with textured surface

Measurement of the Transverse and Longitudinal Viscosities of Continuous Fibre Reinforced Composites

Computational Analysis for Mixing of Fluids Flowing through Micro- Channels of Different Geometries

A NEW DISPERSIVE AND DISTRIBUTIVE STATIC MIXER FOR THE COMPOUNDING OF HIGHLY VISCOUS MATERIALS

Component & Board Level Cooling.

ABSTRACT 1. INTRODUCTION

CE MECHANICS OF FLUIDS UNIT I

Laminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid

Parash Moni Thakur. Gopal Ch. Hazarika

NUMERIC SIMULATION OF ULTRA-HIGH PRESSURE ROTARY ATOMIZING WATERJET FLOW FIELD

MAGNETOHYDRODYNAMICS

Nonlinear Analysis: Modelling and Control, 2008, Vol. 13, No. 4,

Final Polymer Processing

NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE

Wall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes

Computational Methods and Experimental Measurements XII 909

A Finite Element Model for Numerical Analysis of Sintering

Impact of Magnetic Field Strength on Magnetic Fluid Flow through a Channel

THE 3D VISCOELASTIC SIMULATION OF MULTI-LAYER FLOW INSIDE FILM AND SHEET EXTRUSION DIES

Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel

Magnetic Assisted Laser Micromachining for Highly Reflective Metals

Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders

Homework #4 Solution. μ 1. μ 2

NUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP REGIME

COMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE

AMME2261: Fluid Mechanics 1 Course Notes

G. C. Hazarika 2 Department of Mathematics Dibrugarh University, Dibrugarh

Numerical study of entropy generation and melting heat transfer on MHD generalised non-newtonian fluid (GNF): Application to optimal energy

Constructal multi-scale design of compact micro-tube heat sinks and heat exchangers

Polymer Engineering (MM3POE)

Supplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer

Phan-Thien-Tanner Modeling of a Viscoelastic Fluid in the Stick-Slip Scenario. Andrew Griffith

Polymerization Technology Laboratory Course

Chapter 3 Rheological Studies of Polymer Gel Systems

An Essential Requirement in CV Based Industrial Appliances.

PERMEABILITY PREDICTION IN POROUS MEDIA WITH COMPLEX 3D ARCHITECTURES IN A TRI- PERIODIC COMPUTATIONAL DOMAIN

Authors: Correspondence: ABSTRACT: Keywords:

Excerpt from the Proceedings of the COMSOL Users Conference 2006 Boston

MODENA. Deliverable 3.2. WP s leader: TU/e. Simulations for foams, dispersion and mixing and developed SW. Principal investigator:

Numerical simulation of polyurethane foaming processes on bubble scale

Morphology Evolution in PS/LDPE Blends in a Twin Screw Extruder: Effects of Compatibilizer

Levitation by Using Electromagnetic Force

Design optimization of first wall and breeder unit module size for the Indian HCCB blanket module

Liquefaction Resistance and Internal Erosion Potential of Non-Plastic Silty Sand

Defense Technical Information Center Compilation Part Notice

Zhang Yue 1, Guo Wei 2, Wang Xin 3, Li Jiawu 4 1 School of High way, Chang an University, Xi an, Shanxi, China,

Outlines. simple relations of fluid dynamics Boundary layer analysis. Important for basic understanding of convection heat transfer

Physics 54 Lecture March 1, Micro-quiz problems (magnetic fields and forces) Magnetic dipoles and their interaction with magnetic fields

IJSRD - International Journal for Scientific Research & Development Vol. 4, Issue 05, 2016 ISSN (online):

ENMA490 Capstone: Design of Microfluidics Mixer

Flow Structure Investigations in a "Tornado" Combustor

Chapter 5. The Differential Forms of the Fundamental Laws

JJMIE Jordan Journal of Mechanical and Industrial Engineering

Department of Engineering and System Science, National Tsing Hua University,

Studies on flow through and around a porous permeable sphere: II. Heat Transfer

Excerpt from the Proceedings of the COMSOL Conference 2009 Bangalore

CFD Simulation Of Hydrodynamics And Heat Transfer In Flow Of Liquids On Inclined Plates

Study of Aero-thermal Effects with Heat Radiation in Optical Side Window

Dimension measurement. By Mr.Vuttichai Sittiarttakorn

The Meaning and Significance of Heat Transfer Coefficient. Alan Mueller, Chief Technology Officer

Investigation on relation of viscosity and inner structure of suspension under a magnetic field

Transport by convection. Coupling convection-diffusion

ENGR Heat Transfer II

THERMAL HYDRAULIC ANALYSIS IN REACTOR VESSEL INTERNALS CONSIDERING IRRADIATION HEAT

STUDY OF HEAT TRANSFER MECHANISMS DURING THE LENS TM PROCESS

Transcription:

NUMERICAL SIMULATION AND MOLDABILITY INVESTIGATION OF MICRO-FEATURES Shia-Chung Chen¹, ², Li-Chi Su¹, ², Cheng-Yi Chang³, Hsin-Wei Hung³ and Wen-Hsien Yang³ 1. Department of Mechanical Engineering, Chung Yuan Christian University Taiwan. 2. R&D Center for Membrane Technology, Chung Yuan Christian University Taiwam. 3. CoreTech System Co., Ltd., Hsinchu, Taiwan. Abstract In micro-injection molding, the preservation of precise micro-feature is one of the most important indications to ensure proper functionality and quality. A new technique, Induction Heating, which is advanced in heating up the mold quickly and accurately, is adopted to control mold temperature during filling phase. This paper aims to analyze the technique specifically for a part with micro-feature of a high aspect ratio. Meanwhile, it probes into the result of numerical simulation and actual experimental investigation. The result shows that some critical factors have a dominant effect on the molding mechanism, and this result will be beneficial to the development of micro-injection molding technology. Introduction Micro-injection molding can be classified into two categories. One is a micro part which is entirely small; the other is a part with tiny features on a relatively large body. Nowadays, with the advancement of manufacturing technique and economics, products with micro-features, such as light guide plate of LCD, DVD and so on, have been widely prevailed all over the world. However, it is still a challenge for engineers to ensure proper functionality and quality by correctly preserving the micro-feature. To preserve the micro-feature as complete as possible, variotherm technique is often introduced into the molding process. With the technique, the mold surface with the micro-feature is allowed to be heated up to a temperature higher than general recommended mold temperature to keep the polymer from solidification in micro-features during injection. The extra benefit of the technology is that molders can heat up the mold surface efficiently and accurately to achieve the required temperature [3-6]. This paper aims to analyze the moldability of a micro-feature vs. its surface temperature, meanwhile comparing the experiment with the 3D simulated results, in order to identify a proper heat transfer coefficient (HTC) that would be suitable for simulations of micro-injection with similar micro-feature. The true 3D numerical approach developed in the past few years [1-2] is applied to the simulation. Hele-Shaw approach method is not suitable for the simulation because it lacks transverse flow information and the fountain flow effect around the flow advancement. Under the circumstances, the thermal history of a polymer particle could be not fully characterized; it might lead to a significant influence on the temperature estimation. Experiment An experiment of micro-feature molding whose purpose is to get a good reference data regarding the subject had been conducted. The schematic of the experimental part is shown as Fig. 1. The thickness of the part is 1.2 mm. The picture of the mold is shown as Fig. 2. To control the temperature of mold surface which containing the micro-feature, a variotherm technique, "Induction Heating", is introduced to the project. The device is composed of a rectifier, a frequency converter, a transformer and a coil. Following Faraday s and Lenz s Law, having an AC current passing through a heating coil, an alternating magnetic field will be created. Simply by placing the coil near by the mold surface, the difference in magnetic flux through the mold surface will result in Eddie current flow through it. The Eddie current flow in the mold will cause the object to heat up due to the electrical resistance, thus achieving the goal of heating mold surface. The detailed molding conditions are shown in Table 1. Fig. 3 shows the schematics of the mold surface with the micro-feature and induction-heating coil. Before filling, the mold temperature of side surface had always been kept as 60 C, meanwhile the temperature of the core side containing the micro-feature varied from 70 C ~140 C. The experiment is performed with 2 conditions, with and without vacuum. Numerical modeling The polymer melt is assumed to behave as Generalized Newtonian Fluid (GNF). Hence the non-isothermal 3D flow motion can be mathematically described by the followings: (1) Conservation of Mass ρ + ρu = 0 (2) Conservation of Momentum

( ρ u) + ( ρuu = ρg T where = pi + η ( u + u ) (3) Conservation of Energy T ρ C P + u T = 2 ( k T ) + ηγ& where u is the velocity vector, T the temperature, t the time, p the pressure, the total stress tensor, the density, the viscosity, k the thermal conductivity, C p the specific heat and γ& is the shear rate. In this work, the modified-cross model with Arrhenius temperature dependence is employed to describe the viscosity of polymer melt: (4) Modified-Cross model with η ( T, γ& ) η o ( ) η T o = 1 1+ * ( η γ& τ ) n b ( T ) = BExp o T T where n is the power law index, ηo the zero shear viscosity, and * τ is the parameter that describes the transition region between zero shear rate and the power law region of the viscosity curve. A volume fraction function f is introduced to track the evolution of the melt front. Here, f=0 is defined as the air phase, f=1 as the polymer melt phase, and then the melt front is located within cells with 0<f<1. The advancement of f over time is governed by the following transport equation: (5) Transport Equation f + ( u f ) = 0 The flow rate or injection pressure is prescribed at mold inlet. No slip is assumed at mold wall. Note that only inlet boundary condition is necessary for the hyperbolic transport equation of volume fraction function. Simulation Utilizing the simulation of injection molding to evaluate the moldability of micro-feature of part, the most difficult issue is to determine the heat transfer coefficient that can properly describe the non-continuous temperature distribution at both sides of the interface between the melt and the mold. Besides, because of the huge thickness difference between the micro-feature and the plate, the mesh modeling for sufficient resolution with a suitable element count for calculation becomes a challenge as well. To solve the problem, different element topologies were adopted to build the mesh model. Fig. 4 shows the cross section of the mesh in the joint of the plate and the micro-feature. The material used in this study is PMMA-CHIMEI CM205. Fig. 5 (a) and (b) are the viscosity and pvt curves for the material. The process condition of simulation is the same with the experiment, besides, various HTC (from 25000 to 500) were taken into simulation. Results and discussions To precisely measure the height of the micro-feature, a 3D laser digital microscope is applied to scan and measure it. Fig. 6 shows one of the scanning images of a molded micro-feature. Fig. 7 shows the location to measure the altitude of the micro-feature. The experimental result is shown in Fig. 8, which shows that the vacuum condition doesn t have a significant influence on the moldability. The most important phenomenon observed in the result is that the height of the microstructure increases quickly when the temperature of the induction-heated surface exceeds 100 ºC. Also, when the surface is heated up to 110 ºC, the altitude of the structure is drastically increased to more than 90 % of the height. Refer to [3] for the details of the experiments. Simulations with the same process condition but different HTC were conducted for comparisons. The temperature boundary condition of the induction-heating surface is programmed to be the same to meet the experimental condition. Fig. 9 shows the surface with micro-feature where the boundary temperature is controlled to simulate the effect of induction heating. Fig. 10 shows the melt advancement of one of the simulations. With a series of simulations regarding different HTC, a set of simulations where HTC = 1250 is observed to have similar moldability tendency at the same temperature domain. Fig. 11 shows the filling status of the micro-feature with different induction-heating temperature, where the deep blue color indicates the unfilled area. Fig.12, the X-Y plot of the filling altitude with HTC =1250, shows a very good agreement with the experimental result. Fig.13 shows the simulation results with HTC = 25000, 5000 and 750. The correlations with experiment are not as good as HTC=1250 at all. Conclusion Heat transfer coefficient is used to describe the convection heat transfer between fluid and solid. In the study, the exact injection experiments have been conducted as well as the injection simulations. By comparing the results between simulations and experiments, a proper heat transfer coefficient is found in the simulation that shows the same physical tendency with the experiment. The coefficient found in the study will be a good reference for engineers who want to evaluate the moldability of a micro-feature part.

References [1]. R.Y.Chang and W.H.Yang, Numerical Simulation of Mold Filling in Injection Molding Using a Three-Dimensional Finite Volume Approach, Int. J. Numer. Methods Fluids, 37, 125-148 (2001) [2]. L.Liu, W.H. Yang, D.Hsu and V.Yang, ANTEC2002, 621-626 (2002) [3]. S. C. Chen, W. R. Jong, Y. J. Chang, J. A. Chan g, and J. C. Cin, 2006, Rapid Mold Temperature Variation for Assisting Micro Injection of High Aspect Ratio Micro-Feature Parts Using Induction Heating Technology, Journal of Micromechanics and Microengineering, Vol. 16, pp 1783-1791. [4]. Shia-Chung Chen,Wen-Ren Jong, Yaw-Jen Chang, Jen-An Chang and Jin-Chuan Cin, Rapid mold temperature variation for assisting the micro injection of high aspect ratio micro-feature parts using induction heating technology, Journal Of Micromechanics And Microengineering, 1783 1791(2006) [5]. Shia-Chung Chen, Wen-Ren Jong, Jen-An Chang, Dynamic Mold Surface Temperature Control Using Induction Heating and Its Effects on the Surface Appearance of Weld Line, Journal of Applied Polymer Science, Vol. 101, 1174 1180 (2006) [6]. J. A. Chang, S. C. Chen and J. C. Cin, Rapid Mold Temperature Controlon Micro Injection Molded Partswith Highaspect Ratio Micro-Features, ANTEC2006, 1295-1297 (2006) Fig. 1. Schematics of micro-feature part Fig. 2. Micro-feature mold for experiment Keywords micro injection, induction heating, heat transfer coefficient, 3D simulation Acknowledgement Fig. 3. Schematics of micro-feature part and the induction coil This work was partially supported by National Science Council NSC grant 94-2745-E-033-001-URD. Table 1. Experimental molding conditions Experimental molding conditions Material PMMA-CHIMEI:CM205 Inject velocity 100 mm/s Inject ratio 25.434 cm 3 /s Packing pressure 140 MPa Packing time 3 seconds Melt temperature 250 Mold temperature 60.80.90.100.110.120.130.140( )

Fig. 4. Schematics of 3D mesh for micro-structure Fig. 8. The altitude of micro-feature vs. various induction-heating temperature (experiment) Fig. 5. Material properties used in this analysis (a)viscosity (b) pvt Fig. 9. The surface boundary used to simulate induction-heating effect. Fig. 6. Scanning image of the 3D laser digital microscope Fig. 10. The simulated melt advancement Fig. 7. The location to measure the altitude of the micro-feature Fig. 11. Shows the filling status of the micro-feature with different induction-heating temperature with HTC=1250

Fig. 12. The simulated X-Y plot of the filling altitude with HTC =1250 Fig. 13. The simulated X-Y plot of the filling altitude with the other HTC value