Study on Numerical Simulation of Mold Filling and Heat Transfer in Die Casting Process

A 3-D mathematical model considering turbulence phenomena has been established based on a computational fluid dynamics technique, so called 3-D SOLA-VOF (Solution Algorithm-Volume of Fluid), to simulate the fluid flow of mold filling process of die casting. In addition, the mathematical model for simulating the heat transfer in die casting process has also been established. The computation program has been developed by the authors with the finite difference method (FDM) recently. As verification, the mold filling process of a S-shaped die casting has been simulated and the simulation results coincide with that of the benchmark test. Finally, as a practical application, the gating design of a motorcycle component was modified by the mold filling simulation and the dies design of another motorcycle component was optimized by the heat transfer simulation. All the optimized designs were verified by the production practice.

Numerical Simulation on Subcooled Boiling Heat Transfer Characteristics of Water-Cooled W/Cu Divertors

In order to realize safe and stable operation of a water-cooled W/Cu divertor under high heating condition,the exact knowledge of its subcooled boiling heat transfer characteristics under different design parameters is crucial.In this paper,subcooled boiling heat transfer in a water-cooled W/Cu divertor was numerically investigated based on computational fluid dynamic（CFD）.The boiling heat transfer was simulated based on the Euler homogeneous phase model,and local differences of liquid physical properties were considered under one-sided high heating conditions.The calculated wall temperature was in good agreement with experimental results,with the maximum error of 5%only.On this basis,the void fraction distribution,flow field and heat transfer coefficient（HTC）distribution were obtained.The effects of heat flux,inlet velocity and inlet temperature on temperature distribution and pressure drop of a water-cooled W/Cu divertor were also investigated.These results provide a valuable reference for the thermal-hydraulic design of a water-cooled W/Cu divertor.

Numerical Study of Fluid Dynamics and Heat Transfer Induced by Plasma Discharges

A numerical investigation is conducted to explore the evolution of a plasma discharge and its interaction with the fluid flow based on a self-consistent fluid model which couples the discharge dynamics with the fluid dynamics.The effects of the applied voltage on the distribution of velocity and temperature in initially static air are parainetrically studied.Furthermore,the spatial structure of plasma discharge and the resulting force contours in streamwise and normal directions are discussed in detail.The result shows that the plasma actuator produces a net force that should always be directed away from the exposed electrode,which results in an ionic wind pushing particles into a jet downstream of the actuator.When the energy added by the plasma is taken into account,the ambient air temperature is increased slightly around the electrode,but the velocity is almost not affected.Therefore it is unlikely that the induced flow is buoyancy driven.For the operating voltages considered in this paper,the maximum induced velocity is found to follow a power law,i.e.,it is proportional to the applied voltage to the 3.5 power.This promises an efficient application in the flow control with plasma actuators.

Numerical and Experimental Investigation on the In-Flight Melting Behaviour of Granulated Powders in Induction Thermal Plasmas

An innovative in-flight glass melting technology with thermal plasmas was developed for the purpose of energy conservation and environment protection. In this study, modelling and experiments of argon-oxygen induction thermal plasmas were conducted to investigate the melting behaviour of granulated soda-lime glass powders injected into the plasma. A two-dimensional local thermodynamic equilibrium （LTE） model was performed to simulate the heat and momentum transfer between plasma and particle. Results showed that the particle temperature was strongly affected by the flow rate of carrier gas and the particle size of raw material. A higher flow rate of carrier gas led to lower particle temperature and less energy transferred to particles which resulted in lower vitrification. The incomplete melting of large particles was attributed to the lower central temperature of the particle caused by a larger heat capacity. The numerical analysis explained well the experimental results, which can provide valuable practical guidelines for the process control in the melting process for the glass industry.

Numerical simulation of flow and heat transfer in a random packed bed

Random packed beds have more complex interior structure than structured beds and are widely used in industry and engineering. CFD simulation was carried out to investigate and analyze the local flow and heat transfer in a 120-sphere random packed bed. 3D Navier-Stokes equation was solved with a finite volume formulation based on the Chimera meshing technique. Investigation was focused on low Reynolds number flow （Re=4.6-56.2）, which typically occurs in packed bed reactors in bio-chemical fields. Detailed temperature field information was obtained. Inhomogeneity of flow and heat transfer due to the non-uniform distribution of void fraction was discussed and analyzed.

Numerical Study and Performance Analyses of Mems-Based Particle Velocity Sensor with Combined Inclined Rib Pair

The combined inclined rib pair(CIRP)is the first time proposed to improve the sensor performance of particle velocity sensor(PVS)by using a three-dimensional numerical method.The method is verified by the experiment results in the literature.The optimal plain channel parameters are determined as the basic sensor structure.In comparison of plain channel,both heat transfer characteristics and sensor performance are enhanced effectively by arranging the CIRP.The reason is that the high flow rate region caused by the CIRP can maintain strongly in the whole fluid field if there are enough rib pairs.Furthermore,the produced longitudinal vortex pair can get a better fluid mix,which is more conductive to heat transfer.The increasing height and number of the CIRP can improve the heat transfer characteristics,but the flow resistance will increase as well.For the purpose of finding the best overall performance,the effects of the parameters including the geometric sizes and the position of the CIRP have been investigated.The results show that PVS will get the best sensitivity when the rib length and width are 0.2 mm and 0.03 mm respectively,and the distance between rib pair and between ribs in the same pair are 0.15 mm and 0.3 mm respectively.Besides,the most suitable crossing angle is 45°.Thus,the performance of PVS can be significantly improved by this novel structure.

A Numerical Methodology for Enforcing Maximum Principles and the Non-Negative Constraint for Transient Diffusion Equations

Transient diffusion equations arise in many branches of engineering and applied sciences(e.g.,heat transfer and mass transfer),and are parabolic partial differential equations.It is well-known that these equations satisfy important mathematical properties like maximum principles and the non-negative constraint,which have implications in mathematical modeling.However,existing numerical formulations for these types of equations do not,in general,satisfy maximum principles and the nonnegative constraint.In this paper,we present a methodology for enforcing maximum principles and the non-negative constraint for transient anisotropic diffusion equation.The proposed methodology is based on the method of horizontal lines in which the time is discretized first.This results in solving steady anisotropic diffusion equation with decay equation at every discrete time-level.We also present other plausible temporal discretizations,and illustrate their shortcomings in meeting maximum principles and the non-negative constraint.The proposedmethodology can handle general computational grids with no additional restrictions on the time-step.We illustrate the performance and accuracy of the proposed methodology using representative numerical examples.We also perform a numerical convergence analysis of the proposed methodology.For comparison,we also present the results from the standard singlefield semi-discrete formulation and the results froma popular software package,which all will violate maximum principles and the non-negative constraint.

Simulation of steel beam under ceiling jet based on a wind–fire–structure coupling model

For localized fires, it is necessary to consider the thermal and mechanical responses of building elements subject to uneven heating under the influence of wind. In this paper, the thermomechanical phenomena experienced by a ceiling jet and I-beam in a structural fire were simulated. Instead of applying the concept of adiabatic surface temperature (AST) to achieve fluid–structure coupling, this paper proposes a new computational fluid dynamics–finite element method numerical simulation that combines wind, fire, thermal, and structural analyses. First, to analyze the velocity and temperature distributions, the results of the numerical model and experiment were compared in windless conditions, showing good agreement. Vortices were found in the local area formed by the upper and lower flanges of the I-beam and the web, generating a local high-temperature zone and enhancing the heat transfer of convection. In an incoming-flow scenario, the flame was blown askew significantly;the wall temperature was bimodally distributed in the axial direction. The first temperature peak was mainly caused by radiative heat transfer, while the second resulted from convective heat transfer. In terms of mechanical response, the yield strength degradation in the highest-temperature region in windless conditions was found to be significant, thus explaining the stress distribution of steel beams in the fire field. The mechanical response of the overall elements considering the incoming flows was essentially elastic.

Computational fluid dynamics simulation of friction stir welding:A comparative study on different frictional boundary conditions

Numerical simulation based on computational fluid dynamics （CFD） is a useful approach for quantitatively investigating the underlying thermal-mechanical conditions during FSW, such as temperature field and material deformation field. One of the critical issues in CFD simulation of FSW is the use of the frictional boundary condition, which represents the friction between the welding tool and the workpiece in the numerical models. In this study, three-dimensional numerical simulation is conducted to analyze the heat transfer and plastic deformation behaviors during the FSW of AA2024. For comparison purposes, both the boundary velocity （BV） models and the boundary shear stress （BSS） models are employed in order to assess their performances in predicting the temperature and material deformation in FSW. It is interesting to note that different boundary conditions yield similar predictions on temperature, but quite different predictions on material deformation. The numerical predictions are compared with the experimental results. The predicted deformation zone geometry by the BSS model is consistent with the experimental results while there is large difference between the predictions by the BV models and the experimental measurements. The fact that the BSS model yields more reasonable predictions on the deformation zone geometry is attributed to its capacity to automatically adjust the contact state at the tool/workpiece interface. Based on the favorable predictions on both the temperature field and the material deformation field, the BSS model is suggested to have a better performance in numerical simulation of FSW than the BV model.

Multiphase Flow and Thermo-Mechanical Behaviors of Solidifying Shell in Continuous Casting Mold

The metallurgical phenomena occurring in the continuous casting mold have a significant influence on the performance and the quality of steel product.The multiphase flow phenomena of molten steel,steel/slag interface and gas bubbles in the slab continuous casting mold were described by numerical simulation,and the effect of electromagnetic brake（EMBR） and argon gas blowing on the process were investigated.The relationship between wavy fluctuation height near meniscus and the level fluctuation index F,which reflects the situation of mold flux entrapment,was clarified.Moreover,based on a microsegregation model of solute elements in mushy zone with δ/γ transformation and a thermo-mechanical coupling finite element model of shell solidification,the thermal and mechanical behaviors of solidifying shell including the dynamic distribution laws of air gap and mold flux,temperature and stress of shell in slab continuous casting mold were described.