Carreau-Yasuda Model
$$\eta = \eta_0\left[1+(\lambda\dot{\gamma})^a\right]^{\frac{n-1}{a}}$$
Power Law: $\eta = m\,\dot{\gamma}^{n-1}$
Zero-shear viscosity vs melt temperature ±40°C (Arrhenius approximation)
Compute viscosity curves using Carreau-Yasuda and power-law models for common polymers. Calculate gate pressure drop and cooling time to support early-stage injection molding CAE studies.
$$\eta = \eta_0\left[1+(\lambda\dot{\gamma})^a\right]^{\frac{n-1}{a}}$$
Power Law: $\eta = m\,\dot{\gamma}^{n-1}$
Zero-shear viscosity vs melt temperature ±40°C (Arrhenius approximation)
The Carreau-Yasuda model describes how polymer viscosity changes over a wide range of shear rates, capturing the zero-shear plateau, the shear-thinning region, and the transition between them.
$$\eta = \eta_0\left[1+(\lambda\dot{\gamma})^a\right]^{\frac{n-1}{a}}$$η: Viscosity (Pa·s)
η₀: Zero-shear viscosity – the viscosity at very low flow rates.
λ: Relaxation time – dictates when shear-thinning begins.
ṙ: Shear rate (1/s) – the speed gradient within the flow.
n: Flow index – power-law exponent; lower n means more shear-thinning.
a: Yasuda parameter – controls the sharpness of the transition zone.
For pressure drop in a cylindrical gate (a common simplification in injection molding), the power-law form is often used for calculation, where the consistency coefficient (m) is derived from the Carreau-Yasuda model at the relevant shear rate.
$$\Delta P = \frac{2(3n+1)}{n \pi}\left( \frac{Q}{R^3} \right)^n m L$$ΔP: Pressure drop across the gate (Pa).
Q: Volumetric flow rate (m³/s).
R: Gate radius (m).
L: Gate length (m).
m: Consistency coefficient (Pa·sⁿ).
This equation shows how pressure is highly sensitive to gate radius (R³ in the denominator) and the non-Newtonian nature of the polymer (through n and m).
Injection Molding Process Window Definition: Engineers use these curves to define the allowable injection speeds and pressures for a given material and part. Setting the injection speed too low can cause the melt to freeze prematurely; too high can cause degradation from excessive shear heat.
Gate and Runner System Design: The calculated pressure drop is critical for designing balanced runner systems in multi-cavity molds. It ensures all cavities fill at the same time and pressure, preventing defects like short shots or overpacking in specific cavities.
Material Selection and Substitution: When switching polymer grades or suppliers, engineers overlay viscosity curves to ensure the new material will behave similarly in the existing mold. A significant shift in the curve could require retooling the mold.
Predicting Wear on Machine Components: High viscosity at processing shear rates leads to higher pressure and stress on injection screws, check rings, and mold gates. This analysis helps predict maintenance schedules and component life.
When starting with this tool, there are several pitfalls that newcomers to CAE often encounter. The first is understanding that "the parameters for preset materials are not absolute constants." For example, even a material like "PC" can have vastly different viscosity characteristics depending on the manufacturer and grade (high-flow, flame-retardant, etc.). The tool's preset values are merely representative. It's crucial to compare them with the viscosity curve from your actual molding material's datasheet and be prepared to fine-tune λ and n if necessary.
The second point is overlooking the assumptions behind the gate pressure loss calculation. The formula assumes "fully developed laminar flow in a circular channel." This means it does not account for entrance effects at the actual gate constriction or the influence of a frozen layer (frozen skin) at the wall due to mold temperature. You should consider the calculated value as a theoretical minimum pressure loss; practical wisdom suggests applying a safety margin of 1.5 to 2 times this value.
The third point is not over-relying on the cooling time calculation. The tool's cooling time is a theoretical calculation based solely on heat conduction to a mold wall at constant temperature. However, in actual molding, complex factors like cooling channel layout and flow rate, or resin crystallization heat, come into play. Use this calculated value as a reference "to grasp the order of magnitude for cooling." For instance, if the calculation yields 10 seconds, the actual cycle time will often be around 12-15 seconds—understand it in terms of such correlations.
The non-Newtonian fluid viscosity models and flow analysis handled by this tool form the basis for a wide range of engineering fields beyond injection molding. For example, the rheology of paints and cosmetics. The difference in feel when slowly squeezing hair gel from a tube versus quickly spreading it on your hand is precisely shear-rate dependence. In quality control, models similar to the Carreau model used in this tool are employed to evaluate these characteristics.
Another field is food engineering. Ketchup and mayonnaise don't flow when tilted slowly (high zero-shear viscosity) but do flow when shaken or a strong force is applied (shear thinning). This behavior also involves thixotropy, a time-dependent property, but the fundamental understanding of shear-thinning behavior is the same. Furthermore, in hemorheology (bio-rheology), the phenomenon where blood viscosity changes with shear rate inside blood vessels is mathematically analogous to the analysis of polymer melts.
Thus, the "way of thinking about non-Newtonian fluids" you learn with one tool is a universal engineering language that remains applicable even when the material changes.
Once you're comfortable with this tool's calculations and start wondering "why?", it's time to take the next step. First, we recommend getting an overview of the field of "Rheology." In particular, the first Newtonian plateau that appears in the low shear-rate region of a viscosity curve is deeply related to relaxation phenomena in a state of entangled molecular chains. To understand the physical meaning of the parameter "λ (relaxation time)", learning basic viscoelastic models like the Maxwell model can be truly enlightening.
Regarding the mathematical background, try following the derivation of the power law in non-Newtonian fluid mechanics. Starting from an assumed velocity profile for flow in a circular pipe and a force balance on shear stress, you arrive at this foundational equation: $$ \Delta P \cdot \pi r^2 = \tau \cdot 2\pi r L $$ From here, by relating shear stress τ and shear rate $\dot{\gamma}$ via the power law ($\tau = m \dot{\gamma}^n$), you can derive the pressure loss formula used in this tool. Understanding this derivation process will help the meaning of coefficients like "3n+1" in the formula click into place.
For a more practical next topic, challenge yourself with "extending the analysis to various channel shapes in molds (rectangular, thin plate)." There are methods to approximate the formula for circular channels using the concept of hydraulic radius. Mastering this brings you a step closer to predicting more realistic pressure losses in actual mold cavities and runners.