<p>The heat transfer in the filling phase of injection moulding is studied, based on Gunnar Aronsson’s distance model for flow expansion ([Aronsson], 1996).</p><p>The choice of a thermoplastic materials model is motivated by general physical properties, admitting temperature and pressure dependence. Two-phase, per-phase-incompressible, power-law fluids are considered. The shear rate expression takes into account pseudo-radial flow from a point inlet.</p><p>Instead of using a finite element (FEM) solver for the momentum equations a general analytical viscosity expression is used, adjusted to current axial temperature profiles and yielding expressions for axial velocity profile, pressure distribution, frozen layer expansion and special front convection.</p><p>The nonlinear energy partial differential equation is transformed into its conservative form, expressed by the internal energy, and is solved differently in the regions of streaming and stagnant flow, respectively. A finite difference (FD) scheme is chosen using control volume discretization to keep truncation errors small in the presence of non-uniform axial node spacing. Time and pseudo-radial marching is used. A local system of nonlinear FD equations is solved. In an outer iterative procedure the position of the boundary between the “solid” and “liquid” fluid cavity parts is determined. The uniqueness of the solution is claimed. In an inner iterative procedure the axial node temperatures are found. For all physically realistic material properties the convergence is proved. In particular the assumptions needed for the Newton-Mysovskii theorem are secured. The metal mould PDE is locally solved by a series expansion. For particular material properties the same technique can be applied to the “solid” fluid.</p><p>In the circular plate application, comparisons with the commercial FEM-FD program Moldflow (Mfl) are made, on two Mfl-database materials, for which model parameters are estimated/adjusted. The resulting time evolutions of pressures and temperatures are analysed, as well as the radial and axial profiles of temperature and frozen layer. The greatest differences occur at the flow front, where Mfl neglects axial heat convection. The effects of using more and more complex material models are also investigated. Our method performance is reported.</p><p>In the polygonal star-shaped plate application a geometric cavity model is developed. Comparison runs with the commercial FEM-FD program Cadmould (Cmd) are performed, on two Cmd-database materials, in an equilateral triangular mould cavity, and materials model parameters are estimated/adjusted. The resulting average temperatures at the end of filling are compared, on rays of different angular deviation from the closest corner ray and on different concentric circles, using angular and axial (cavity-halves) symmetry. The greatest differences occur in narrow flow sectors, fatal for our 2D model for a material with non-realistic viscosity model. We present some colour plots, e.g. for the residence time.</p><p>The classical square-root increase by time of the frozen layer is used for extrapolation. It may also be part of the front model in the initial collision with the cold metal mould. An extension of the model is found which describes the radial profile of the frozen layer in the circular plate application accurately also close to the inlet.</p><p>The well-posedness of the corresponding linearized problem is studied, as well as the stability of the linearized FD-scheme.</p> / Report code: LiU-TEK-LIC-2002:66.

The heat transfer in the filling phase of injection moulding is studied, based on Gunnar Aronsson’s distance model for flow expansion ([Aronsson], 1996). The choice of a thermoplastic materials model is motivated by general physical properties, admitting temperature and pressure dependence. Two-phase, per-phase-incompressible, power-law fluids are considered. The shear rate expression takes into account pseudo-radial flow from a point inlet. Instead of using a finite element (FEM) solver for the momentum equations a general analytical viscosity expression is used, adjusted to current axial temperature profiles and yielding expressions for axial velocity profile, pressure distribution, frozen layer expansion and special front convection. The nonlinear energy partial differential equation is transformed into its conservative form, expressed by the internal energy, and is solved differently in the regions of streaming and stagnant flow, respectively. A finite difference (FD) scheme is chosen using control volume discretization to keep truncation errors small in the presence of non-uniform axial node spacing. Time and pseudo-radial marching is used. A local system of nonlinear FD equations is solved. In an outer iterative procedure the position of the boundary between the “solid” and “liquid” fluid cavity parts is determined. The uniqueness of the solution is claimed. In an inner iterative procedure the axial node temperatures are found. For all physically realistic material properties the convergence is proved. In particular the assumptions needed for the Newton-Mysovskii theorem are secured. The metal mould PDE is locally solved by a series expansion. For particular material properties the same technique can be applied to the “solid” fluid. In the circular plate application, comparisons with the commercial FEM-FD program Moldflow (Mfl) are made, on two Mfl-database materials, for which model parameters are estimated/adjusted. The resulting time evolutions of pressures and temperatures are analysed, as well as the radial and axial profiles of temperature and frozen layer. The greatest differences occur at the flow front, where Mfl neglects axial heat convection. The effects of using more and more complex material models are also investigated. Our method performance is reported. In the polygonal star-shaped plate application a geometric cavity model is developed. Comparison runs with the commercial FEM-FD program Cadmould (Cmd) are performed, on two Cmd-database materials, in an equilateral triangular mould cavity, and materials model parameters are estimated/adjusted. The resulting average temperatures at the end of filling are compared, on rays of different angular deviation from the closest corner ray and on different concentric circles, using angular and axial (cavity-halves) symmetry. The greatest differences occur in narrow flow sectors, fatal for our 2D model for a material with non-realistic viscosity model. We present some colour plots, e.g. for the residence time. The classical square-root increase by time of the frozen layer is used for extrapolation. It may also be part of the front model in the initial collision with the cold metal mould. An extension of the model is found which describes the radial profile of the frozen layer in the circular plate application accurately also close to the inlet. The well-posedness of the corresponding linearized problem is studied, as well as the stability of the linearized FD-scheme. / <p>Report code: LiU-TEK-LIC-2002:66.</p>

Some normally consolidated soft soils manifest strength sensitivity, ie these soil manifest strain softening when shear in an undrained mode. These soils, referred to as sensitive soft soils, have the typical features of strain hardening in drained shearing and strain softening in undrained shearing. The consolidation lines of these soils are also curved (concave upwards) in the semi-log space. However, under high consolidation stress or upon large shearing, these soils re-gain the features of re-constituted soil. Ground improvement methods like stone columns were reported as not effective when installed in the sensitive soft clays. But mechanism of the un-effectiveness of the stone columns remains unknown because of lack of a suitable and simple model for simulating the stress-strain behaviours of sensitive soft soils. Although these soils have a meta-stable micro-structure, models that developed for simulating structured firm soils are not suitable for simulating sensitive soft soil features. Thus, a new model was formulated. The new model can degenerate back to a Modified Cam Clay model. The ability of new model in simulating a range of behaviour was verified by using the finite difference (FD) method in solving the partial differential equations of the soil model for a range of tri-axial test conditions. The model was further implemented in coupled analysis formulation and coded into FEM program AFENA. Various cases with different soil parameters were then simulated and compared with the FD solutions for various triaxial tests so as to check the stability of the FEM code. The coupled FEA was then used to simulate the performance of geosynthetic-encased stone columns. A new stone column element and a geo-encasement element were developed and coded into AFENA. The stone column simulations were then done for both non-sensitive soils (represented by Modified Cam Clay model) and sensitive soft soil (represented by the new model). Parametric study was conducted to examine the performance of the geo-encased stone columns in both types of soils. Furthermore, two different installation methods: wished-in installation and full displacement installation were studied numerically. Cross comparison was done to investigate how the sensitive soft soil features interact with the installation method in affecting the performance of the geo-encased stone columns. A range of factors that influence the geosynthetic-encased stone columns performance installed in soft soils were also made clear.