{"id":4803,"date":"2024-06-24T14:21:41","date_gmt":"2024-06-24T13:21:41","guid":{"rendered":"https:\/\/www.springboard.pro\/?p=4803"},"modified":"2026-09-29T09:08:09","modified_gmt":"2026-09-29T09:08:09","slug":"through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations","status":"publish","type":"blog","link":"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/","title":{"rendered":"Through Thick and Thin: Modeling the Delivery of Non-Newtonian Formulations"},"content":{"rendered":"\n<section id=\"armstrong-block-editor-block_f0530a45299ede5f738091738cfb091b\" class=\"armstrong-block armstrong-block-editor\" >\n    <div class=\"container-fluid\">\n        <div class=\"row\">\n            <div class=\"col-lg-12\">\n                                <p>The pharmaceutical industry is currently undergoing a period of rapid innovation, propeled by several factors. These include advancements in the development of biologics, which are becoming increasingly diverse<sup>1<\/sup>, decentralized healthcare practices, and a shift towards administering injections less frequently in home settings<sup>2<\/sup>. Additionally, small molecules are being designed through AI drug discovery techniques<sup>3<\/sup>. While these therapeutic innovations hold promise, they can also bring forth unforeseen challenges in formulation and delivery.<\/p>\n<p>These can include ensuring the stability of large or fragile biologics at higher therapeutic concentrations or improving the solubility of poorly soluble APIs. Overcoming these hurdles can require the implementation of sophisticated solutions such as microcarriers, adjuvants, excipients, and stabilizing matrices<sup>4<\/sup>. It&rsquo;s crucial to recognize that these formulations profoundly impact the physical properties of drugs and can present challenges for parenteral drug delivery, especially in rapid, automated systems like autoinjectors.<\/p>\n<p>One example of this are high-concentration monoclonal antibodies which can experience aggregation under agitation-induced stress. Stabilizing matrices composed of polymer, lipid crystal, or polysaccharide gels, offer improved stability and efficacy but can exhibit non-Newtonian flow behavior. This is also true for long-acting injections composed of particulate suspensions where particle-particle interactions and alignment and result in the same.<\/p>\n<p>The modelling of syringeability and delivery of parenteral formulations has often focused on Newtonian models and the Haggen-Poiseuille equation. While this generally performs well, it is not suitable for predicting the behavior of non-Newtonian fluids. This is because non-Newtonian fluids do not uniformly thicken or thin as they pass through a needle; rather, this effect occurs near the needle wall where the fluid experiences the highest shear.<\/p>\n<p>For shear-thinning fluids, this phenomenon results in plug-like flow, which can be seen in the dispensation of fluids such as toothpaste. This results in a very different volumetric flowrate at a given pressure. In autoinjectors, this issue is further complicated by the decay of delivery force over the stroke, which is characteristic of most devices powered by springs or gas.<\/p>\n<p>To tackle this challenge, this paper presents an overview of solutions to non-Newtonian fluids described by the power-law, Cross and Carreau models. The aim is to provide a toolbox to assist the integration of formulation development and device design. Specifically, we examine model solutions in the context of an autoinjector where the delivery force varies throughout the delivery process. Adopting a standardized approach to predictive modeling can enable the industry to make well-informed decisions early on and provide valuable insights into potential device parameters during formulation development.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Autoinjector-scaled-1.png\" sizes=\"auto, (max-width: 2560px) 100vw, 2560px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Autoinjector-scaled-1.png 2560w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Autoinjector-scaled-1-300x190.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Autoinjector-scaled-1-1024x649.png 1024w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Autoinjector-scaled-1-768x487.png 768w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Autoinjector-scaled-1-1536x974.png 1536w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Autoinjector-scaled-1-2048x1298.png 2048w\" alt=\"\" width=\"2560\" height=\"1623\" \/><\/p>\n<p><strong>Figure 1: Illustration of autoinjector parameters for the definition of injection time, omitting losses in available force to bubble compression and friction.<\/strong><\/p>\n<p>&nbsp;<\/p>\n<h3>Defining the parameters of an autoinjector<\/h3>\n<p>&nbsp;<\/p>\n<p>Consider an autoinjector (Figure 1) comprising a primary container with a liquid volume <em>V<\/em> and an internal diameter <em>D<\/em>. For illustration it is assumed that any bubble volume in the container vanishes once the container is pressurized, although this can be accounted for if needed. The liquid occupies a fraction of the container length L equivalent to the delivery stroke:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-1.png\" alt=\"\" width=\"297\" height=\"65\" \/><\/p>\n<p>In this example the work required to deliver the drug is performed by the force of a compression spring which follows Hooke&rsquo;s law, pushing on a rubber plunger stopper within the container:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-2.png\" sizes=\"auto, (max-width: 303px) 100vw, 303px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-2.png 303w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-2-300x55.png 300w\" alt=\"\" width=\"303\" height=\"56\" \/><\/p>\n<p>Here <em>F<\/em> is the force at a particular instance during delivery, <em>F<\/em><sub>0<\/sub>\u00a0\u00a0is the force at the start of delivery, <em>k <\/em>is the spring constant and x represents the spring extension, which for simplicity is assumed to correlate to the plunger stopper displacement. The system generates pressure within the container, which due to the spring constant <em>k<\/em> diminishes as delivery progresses. This case is true for most autoinjectors powered by helical compression springs or gas springs. In the event of a system providing constant force such as an electronic driver or subliming substance, the delivery occurs at near steady state.<\/p>\n<p>The force <em>F<\/em> exerted by the delivery system after accounting for frictional losses (glide force <em>F<sub>f<\/sub><\/em>) creates a pressure<em> P<\/em> within the container proportional to its cross-sectional area <em>A. <\/em>This produces a flow rate <em>Q<\/em> through the needle proportional to the hydrodynamic resistance of the needle (length and diameter) as well as the effective viscosity of the fluid. We assume here that no back pressure is present and therefore that P within the container is equivalent to <em>\u0394P<\/em> across the needle.<\/p>\n<h3>Modelling non-Newtonian fluids<\/h3>\n<p>The viscosity and shear relationships for non-Newtonian fluids<sup>5,6<\/sup> described by the power-law, Carreau and Cross fluid models are shown in Table 1.<\/p>\n<h4>Defining volumetric flowrate<\/h4>\n<p>The local velocity of the liquid in a needle depends on the distance r from the center axis. Here a definition of volumetric flow rate adapted from a method attributed to Weissenberg, Rabinowitsch, Mooney and Schofield (WRMS)<sup>5,6<\/sup> is used, and combined it with definitions for an autoinjector mechanism:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-3.png\" sizes=\"auto, (max-width: 1772px) 100vw, 1772px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-3.png 1772w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-3-300x27.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-3-1024x93.png 1024w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-3-768x70.png 768w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-3-1536x140.png 1536w\" alt=\"\" width=\"1772\" height=\"161\" \/><\/p>\n<p>From this, an integral <em>I<\/em> is derived<sup>6<\/sup>. It is this integral that requires solving to determine the volumetric flow rate for a specific fluid model:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-4.png\" sizes=\"auto, (max-width: 339px) 100vw, 339px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-4.png 339w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-4-300x127.png 300w\" alt=\"\" width=\"339\" height=\"143\" \/><\/p>\n<p>Substituting specific terms for \u03c4 and d\u03c4 followed by integration provides an analytical solution for a model of interest. In the case of a Newtonian fluid the results in the familiar Hagge-Poiseuille equation:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-5.png\" sizes=\"auto, (max-width: 863px) 100vw, 863px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-5.png 863w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-5-300x127.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-5-768x326.png 768w\" alt=\"\" width=\"863\" height=\"366\" \/><\/p>\n<h3>Non-Newtonian behavior<\/h3>\n<p>It is useful to initially define power-law fluid behavior, as it underpins both the Carreau and Cross fluid models. A power-law fluid flow profile is not parabolic (plug flow) and affects its volumetric flow rate at a given differential pressure. The magnitude of shear dependency is governed by a flow behavior index (n), while the overall viscosity is regulated by the flow consistency index (K). Unlike the Cross and Carreau models, the power-law fluid equation for injection time can be analytically solved, as previously reported<sup>7<\/sup>. The normalized transition from parabolic to plug flow, as the value of <em>n<\/em> decreases, is illustrated in Figure 2.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/figure-2.jpg\" sizes=\"auto, (max-width: 1000px) 100vw, 1000px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/figure-2.jpg 1000w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/figure-2-300x241.jpg 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/figure-2-768x616.jpg 768w\" alt=\"\" width=\"1000\" height=\"802\" \/><\/p>\n<p><strong>Figure 2 \u2013 Illustration of the change from a parabolic flow profile to plug flow for a power-law fluid with decreasing flow behavior index n.<\/strong><\/p>\n<p>The analytical solution for injection time is found from the definition of the Power Law fluid volumetric flow rate:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-6.png\" sizes=\"auto, (max-width: 1104px) 100vw, 1104px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-6.png 1104w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-6-300x101.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-6-1024x343.png 1024w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-6-768x257.png 768w\" alt=\"\" width=\"1104\" height=\"370\" \/><\/p>\n<p>Expanding for <em>I<\/em> and substituting container Diameter D for container radius R gives:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-7.png\" sizes=\"auto, (max-width: 810px) 100vw, 810px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-7.png 810w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-7-300x82.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-7-768x210.png 768w\" alt=\"\" width=\"810\" height=\"221\" \/><\/p>\n<p>Where <em>F <\/em>is the applied force, <em>D<\/em> is the container diameter, <em>d<\/em> is needle diameter,<em> l<\/em> is needle length, <em>K <\/em>flow consistency index and <em>n<\/em> the flow behavior index.<\/p>\n<p>This is combined with the two equations describing the autoinjector to obtain a relation for plunger motion over time<sup>7<\/sup> (note k for spring index is similar in appearance to K the flow consistency index):<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-8.png\" sizes=\"auto, (max-width: 1003px) 100vw, 1003px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-8.png 1003w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-8-300x106.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-8-768x272.png 768w\" alt=\"\" width=\"1003\" height=\"355\" \/><\/p>\n<p>For the purpose of integration this is simplified by introducing a constant C:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-9.png\" sizes=\"auto, (max-width: 1340px) 100vw, 1340px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-9.png 1340w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-9-300x32.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-9-1024x110.png 1024w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-9-768x83.png 768w\" alt=\"\" width=\"1340\" height=\"144\" \/><\/p>\n<p>&nbsp;<\/p>\n<p>Integration over the stroke results in an analytical solution for injection time for a Power Law fluid from an autoinjector with a specific needle geometry and spring design:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-10.png\" sizes=\"auto, (max-width: 1737px) 100vw, 1737px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-10.png 1737w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-10-300x74.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-10-1024x253.png 1024w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-10-768x190.png 768w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-10-1536x380.png 1536w\" alt=\"\" width=\"1737\" height=\"430\" \/><\/p>\n<h4>Formulations with time-dependent non-Newtonian behavior<\/h4>\n<p>Cross and Carreau are more complex, behaving as Power law fluid at intermediate rates and Newtonian fluids at both high and low shear rates (Figure 3). They require the definition of four-parameters including low-shear and high-shear viscosities , a characteristic time <em>\u03bb<\/em> and a flow behavior index <em>n<\/em>.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Figure-3.jpg\" sizes=\"auto, (max-width: 1000px) 100vw, 1000px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Figure-3.jpg 1000w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Figure-3-300x220.jpg 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Figure-3-768x562.jpg 768w\" alt=\"\" width=\"1000\" height=\"732\" \/><\/p>\n<p><strong>Figure 3 \u2013 Illustration of the Power law, Cross and Carreau apparent viscosity as a function of shear rate.<\/strong><\/p>\n<p>Finding an analytical solution as the one shown for a power-law fluid is not possible because the relationship between flowrate and pressure drop needs to be defined implicitly, requiring a numerical solver.<\/p>\n<p>The two solutions start with the same definition from the WRMS method<sup>4,5<\/sup> as before, where:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-11.png\" sizes=\"auto, (max-width: 880px) 100vw, 880px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-11.png 880w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-11-300x50.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-11-768x127.png 768w\" alt=\"\" width=\"880\" height=\"146\" \/><\/p>\n<h3>Carreau Fluid:<\/h3>\n<p>To solve the integral and calculate volumetric flowrate the equation defining viscosity for the Carreau fluid is introduced, substituting <em>\u03b4<\/em>= (<em>\u03bc<\/em><sub>0<\/sub>\u2013<em>\u03bc<\/em><sub>\u221e<\/sub>) and\u00a0<em>n<\/em>\u00a0-1 =\u00a0<em>n<\/em>\u2019\u00a0 :<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-12.png\" sizes=\"auto, (max-width: 1993px) 100vw, 1993px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-12.png 1993w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-12-300x46.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-12-1024x158.png 1024w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-12-768x119.png 768w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-12-1536x237.png 1536w\" alt=\"\" width=\"1993\" height=\"308\" \/><\/p>\n<p>&nbsp;<\/p>\n<p>Substituting \u03c4 and\u00a0\u00a0<em>d<\/em>\u03c4\u00a0\u00a0into the definition of integral <em>I<\/em>, changing the integration limits, and solving analytically results in a large equation (See reference [6] equation 30).<\/p>\n<p>The only thing needed to calculate \/ is \u00a0 <em>\u03b3 \u0307<sub>w<\/sub><\/em> \u00a0which can be obtained numerically using:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-13.png\" sizes=\"auto, (max-width: 761px) 100vw, 761px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-13.png 761w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-13-300x51.png 300w\" alt=\"\" width=\"761\" height=\"129\" \/><\/p>\n<h3>Cross Fluid:<\/h3>\n<p>The solution of a cross fluid is similar:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-14.png\" sizes=\"auto, (max-width: 1727px) 100vw, 1727px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-14.png 1727w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-14-300x47.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-14-1024x162.png 1024w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-14-768x121.png 768w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-14-1536x243.png 1536w\" alt=\"\" width=\"1727\" height=\"273\" \/><\/p>\n<p>These are substituted into <em>I<\/em> \u00a0and solved (See reference [6] equation 39 and 40). As with the Carreau model, only\u00a0 <em>\u03b3 \u0307<sub>w<\/sub><\/em>\u00a0\u00a0is required to do this numerically using the equation below:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-15.png\" sizes=\"auto, (max-width: 670px) 100vw, 670px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-15.png 670w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-15-300x63.png 300w\" alt=\"\" width=\"670\" height=\"140\" \/><\/p>\n<h3>Application to an autoinjector problem<\/h3>\n<p>In an autoinjector <em>Q<\/em> can be converted to the rate of change of the plunger position over time <em>dx\/dt<\/em>, which correlates to an incremental reduction in available delivery force. The computation of injection time for Cross and Carreau fluids requires an incremental, numerical approach.<\/p>\n<p>Initially the autoinjector system is considered at <em>x = 0<\/em> and <em>t = 0<\/em>. The procedure is to starts with the determination of <em>\u0394P<\/em>, <em>\u00a0I <\/em>\u00a0and finally flowrate <em>Q <\/em>at <em>x(t=0) = 0<\/em> and <em>F(x=0) = F<sub>0<\/sub><\/em>. <em>Q <\/em>is subsequently converted into a change of plunger stopper position over an increment. The new level of force at this spring extension is used to again calculate <em>\u0394P<\/em>, ,<em> I <\/em>\u00a0and Q. Injection time is the sum of time intervals until the point where x = L, with smaller increments being preferable.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-16.png\" sizes=\"auto, (max-width: 1162px) 100vw, 1162px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-16.png 1162w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-16-300x28.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-16-1024x94.png 1024w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Equation-16-768x71.png 768w\" alt=\"\" width=\"1162\" height=\"107\" \/><\/p>\n<p>If the autoinjector supplies a constant force, the system can be treated as steady state over the interval from x=0 to x=L.<\/p>\n<p>Determination of constants<\/p>\n<p>The determination of the best fitting model as well as the relevant model parameters is achieved by fitting the models presented in Table 1 to data of log dynamic viscosity against the log of the shear rate. Solvers in programming languages such as MATLAB or Python as well as other commercially available rheology software can be used for this. Before analyzing the data, it is important to consider if the experimental conditions cover a sufficiently large span of shear rates. This is to ensure that any Newtonian to non-Newtonian transitions are identified, and that the shear rates are representative of those experienced within an autoinjector.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"full-width-image\" src=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Table-1-1.png\" sizes=\"auto, (max-width: 2107px) 100vw, 2107px\" srcset=\"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Table-1-1.png 2107w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Table-1-1-300x162.png 300w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Table-1-1-1024x552.png 1024w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Table-1-1-768x414.png 768w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Table-1-1-1536x827.png 1536w, https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/Table-1-1-2048x1103.png 2048w\" alt=\"\" width=\"2107\" height=\"1135\" \/><\/p>\n<h3>Conclusion:<\/h3>\n<p>The models provided in this article are a step towards a common toolbox for drug delivery device design and process development. Understanding the non-Newtonian behavior of formulations is essential for ensuring the efficacy and safety of advanced therapies. By embracing standardized modelling approaches and considering the complexities of non-Newtonian behavior early in the development process, researchers and developers can de-risk future device development. This becomes more relevant as the industry focuses on established platform devices which may one day have to deliver increasingly more sophisticated formulations in the future.<\/p>\n<p>Written by Alex Vasiev, Head of Biointerface, Sanner UK and featured in <a href=\"https:\/\/ondrugdelivery.com\/through-thick-and-thin-modelling-the-delivery-of-non-newtonian-formulations\/\">ONdrugDelivery<\/a>; edited to appear here with American writing conventions<\/p>\n<p>&nbsp;<\/p>\n<p>References:<\/p>\n<ol>\n<li>Whitepaper: \u201cThe Rise of Biologics: Emerging Trends and Opportunities\u201d, American Chemical Society, (2023).<\/li>\n<li>Baryakova, T.H., Pogostin, B.H., Langer, R. et al., \u201cOvercoming barriers to patient adherence: the case for developing innovative drug delivery systems\u201d. Nat Rev Drug Discov 22, 387\u2013409 (2023).<\/li>\n<li>Debleena P., Gaurav S., Snehal S., Dnyaneshwar K., Kiran K., Tekade R.K., \u201cArtificial intelligence in drug discovery and development.\u201d Drug discovery today vol. 26,1 (2021): 80-93.<\/li>\n<li>Jindal A. B., Bhide A. R., Salave S., Rana D., Benival D., \u201cLong-acting parenteral drug delivery systems for the treatment of chronic diseases\u201d, Advanced Drug Delivery Reviews, (2023), Volume 198,<\/li>\n<li>A.H.P., \u201cNon-Newtonian Flow and Heat Transfer\u201d. John Wiley and Sons Inc., (1967).<\/li>\n<li>Sochi, T., \u201cAnalytical solutions for the flow of Carreau and Cross fluids in circular pipes and thin slits.\u201d Rheologica Acta. (2015). 10.1007\/s00397-015-0863-x.<\/li>\n<li>Vilaplana M., Vasiev A., White S., \u201cFunction-Based Primary Packaging Design for Injecting Viscous Non-Newtonian Formulations\u201d. ONdrugDelivery Magazine, Issue 101 (Oct 2019), pp 20-25.<\/li>\n<\/ol>\n            <\/div>\n        <\/div>\n    <\/div>\n<\/section>","protected":false},"featured_media":7681,"template":"","blog_category":[75],"class_list":["post-4803","blog","type-blog","status-publish","has-post-thumbnail","hentry"],"acf":[],"yoast_head":"<!-- This site is optimized with the Yoast SEO plugin v28.6 - https:\/\/yoast.com\/product\/yoast-seo-wordpress\/ -->\n<title>Through Thick and Thin: Modeling the Delivery of Non-Newtonian Formulations - Sanner Group<\/title>\n<meta name=\"description\" content=\"How can you ensure the stability of biologic non-newtonian drug formulations? Read the full article to find out.\" \/>\n<meta name=\"robots\" content=\"index, follow, max-snippet:-1, max-image-preview:large, max-video-preview:-1\" \/>\n<link rel=\"canonical\" href=\"https:\/\/www.sanner-group.com\/fr\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/\" \/>\n<meta property=\"og:locale\" content=\"fr_FR\" \/>\n<meta property=\"og:type\" content=\"article\" \/>\n<meta property=\"og:title\" content=\"Through Thick and Thin: Modeling the Delivery of Non-Newtonian Formulations - Sanner Group\" \/>\n<meta property=\"og:description\" content=\"How can you ensure the stability of biologic non-newtonian drug formulations? Read the full article to find out.\" \/>\n<meta property=\"og:url\" content=\"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/\" \/>\n<meta property=\"og:site_name\" content=\"Sanner Group\" \/>\n<meta property=\"article:modified_time\" content=\"2026-09-29T09:08:09+00:00\" \/>\n<meta name=\"twitter:card\" content=\"summary_large_image\" \/>\n<meta name=\"twitter:label1\" content=\"Dur\u00e9e de lecture estim\u00e9e\" \/>\n\t<meta name=\"twitter:data1\" content=\"1 minute\" \/>\n<script type=\"application\/ld+json\" class=\"yoast-schema-graph\">{\"@context\":\"https:\\\/\\\/schema.org\",\"@graph\":[{\"@type\":\"WebPage\",\"@id\":\"https:\\\/\\\/www.sanner-group.com\\\/insights\\\/blog\\\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\\\/\",\"url\":\"https:\\\/\\\/www.sanner-group.com\\\/insights\\\/blog\\\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\\\/\",\"name\":\"Through Thick and Thin: Modeling the Delivery of Non-Newtonian Formulations - Sanner Group\",\"isPartOf\":{\"@id\":\"https:\\\/\\\/www.sanner-group.com\\\/fr\\\/#website\"},\"primaryImageOfPage\":{\"@id\":\"https:\\\/\\\/www.sanner-group.com\\\/insights\\\/blog\\\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\\\/#primaryimage\"},\"image\":{\"@id\":\"https:\\\/\\\/www.sanner-group.com\\\/insights\\\/blog\\\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\\\/#primaryimage\"},\"thumbnailUrl\":\"https:\\\/\\\/www.sanner-group.com\\\/wp-content\\\/uploads\\\/2024\\\/06\\\/2022-02-21-CMF-scaled-1.avif\",\"datePublished\":\"2024-06-24T13:21:41+00:00\",\"dateModified\":\"2026-09-29T09:08:09+00:00\",\"description\":\"How can you ensure the stability of biologic non-newtonian drug formulations? Read the full article to find out.\",\"breadcrumb\":{\"@id\":\"https:\\\/\\\/www.sanner-group.com\\\/insights\\\/blog\\\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\\\/#breadcrumb\"},\"inLanguage\":\"fr-FR\",\"potentialAction\":[{\"@type\":\"ReadAction\",\"target\":[\"https:\\\/\\\/www.sanner-group.com\\\/insights\\\/blog\\\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\\\/\"]}]},{\"@type\":\"ImageObject\",\"inLanguage\":\"fr-FR\",\"@id\":\"https:\\\/\\\/www.sanner-group.com\\\/insights\\\/blog\\\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\\\/#primaryimage\",\"url\":\"https:\\\/\\\/www.sanner-group.com\\\/wp-content\\\/uploads\\\/2024\\\/06\\\/2022-02-21-CMF-scaled-1.avif\",\"contentUrl\":\"https:\\\/\\\/www.sanner-group.com\\\/wp-content\\\/uploads\\\/2024\\\/06\\\/2022-02-21-CMF-scaled-1.avif\",\"width\":2560,\"height\":1919},{\"@type\":\"BreadcrumbList\",\"@id\":\"https:\\\/\\\/www.sanner-group.com\\\/insights\\\/blog\\\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\\\/#breadcrumb\",\"itemListElement\":[{\"@type\":\"ListItem\",\"position\":1,\"name\":\"Home\",\"item\":\"https:\\\/\\\/www.sanner-group.com\\\/fr\\\/\"},{\"@type\":\"ListItem\",\"position\":2,\"name\":\"Through Thick and Thin: Modeling the Delivery of Non-Newtonian Formulations\"}]},{\"@type\":\"WebSite\",\"@id\":\"https:\\\/\\\/www.sanner-group.com\\\/fr\\\/#website\",\"url\":\"https:\\\/\\\/www.sanner-group.com\\\/fr\\\/\",\"name\":\"Sanner Group\",\"description\":\"\",\"potentialAction\":[{\"@type\":\"SearchAction\",\"target\":{\"@type\":\"EntryPoint\",\"urlTemplate\":\"https:\\\/\\\/www.sanner-group.com\\\/fr\\\/?s={search_term_string}\"},\"query-input\":{\"@type\":\"PropertyValueSpecification\",\"valueRequired\":true,\"valueName\":\"search_term_string\"}}],\"inLanguage\":\"fr-FR\"}]}<\/script>\n<!-- \/ Yoast SEO plugin. -->","yoast_head_json":{"title":"Through Thick and Thin: Modeling the Delivery of Non-Newtonian Formulations - Sanner Group","description":"How can you ensure the stability of biologic non-newtonian drug formulations? Read the full article to find out.","robots":{"index":"index","follow":"follow","max-snippet":"max-snippet:-1","max-image-preview":"max-image-preview:large","max-video-preview":"max-video-preview:-1"},"canonical":"https:\/\/www.sanner-group.com\/fr\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/","og_locale":"fr_FR","og_type":"article","og_title":"Through Thick and Thin: Modeling the Delivery of Non-Newtonian Formulations - Sanner Group","og_description":"How can you ensure the stability of biologic non-newtonian drug formulations? Read the full article to find out.","og_url":"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/","og_site_name":"Sanner Group","article_modified_time":"2026-09-29T09:08:09+00:00","twitter_card":"summary_large_image","twitter_misc":{"Dur\u00e9e de lecture estim\u00e9e":"1 minute"},"schema":{"@context":"https:\/\/schema.org","@graph":[{"@type":"WebPage","@id":"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/","url":"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/","name":"Through Thick and Thin: Modeling the Delivery of Non-Newtonian Formulations - Sanner Group","isPartOf":{"@id":"https:\/\/www.sanner-group.com\/fr\/#website"},"primaryImageOfPage":{"@id":"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/#primaryimage"},"image":{"@id":"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/#primaryimage"},"thumbnailUrl":"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/2022-02-21-CMF-scaled-1.avif","datePublished":"2024-06-24T13:21:41+00:00","dateModified":"2026-09-29T09:08:09+00:00","description":"How can you ensure the stability of biologic non-newtonian drug formulations? Read the full article to find out.","breadcrumb":{"@id":"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/#breadcrumb"},"inLanguage":"fr-FR","potentialAction":[{"@type":"ReadAction","target":["https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/"]}]},{"@type":"ImageObject","inLanguage":"fr-FR","@id":"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/#primaryimage","url":"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/2022-02-21-CMF-scaled-1.avif","contentUrl":"https:\/\/www.sanner-group.com\/wp-content\/uploads\/2024\/06\/2022-02-21-CMF-scaled-1.avif","width":2560,"height":1919},{"@type":"BreadcrumbList","@id":"https:\/\/www.sanner-group.com\/insights\/blog\/through-thick-and-thin-modeling-the-delivery-of-non-newtonian-formulations\/#breadcrumb","itemListElement":[{"@type":"ListItem","position":1,"name":"Home","item":"https:\/\/www.sanner-group.com\/fr\/"},{"@type":"ListItem","position":2,"name":"Through Thick and Thin: Modeling the Delivery of Non-Newtonian Formulations"}]},{"@type":"WebSite","@id":"https:\/\/www.sanner-group.com\/fr\/#website","url":"https:\/\/www.sanner-group.com\/fr\/","name":"Sanner Group","description":"","potentialAction":[{"@type":"SearchAction","target":{"@type":"EntryPoint","urlTemplate":"https:\/\/www.sanner-group.com\/fr\/?s={search_term_string}"},"query-input":{"@type":"PropertyValueSpecification","valueRequired":true,"valueName":"search_term_string"}}],"inLanguage":"fr-FR"}]}},"permalink_manager":null,"_links":{"self":[{"href":"https:\/\/www.sanner-group.com\/fr\/wp-json\/wp\/v2\/blog\/4803","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sanner-group.com\/fr\/wp-json\/wp\/v2\/blog"}],"about":[{"href":"https:\/\/www.sanner-group.com\/fr\/wp-json\/wp\/v2\/types\/blog"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.sanner-group.com\/fr\/wp-json\/wp\/v2\/media\/7681"}],"wp:attachment":[{"href":"https:\/\/www.sanner-group.com\/fr\/wp-json\/wp\/v2\/media?parent=4803"}],"wp:term":[{"taxonomy":"blog_category","embeddable":true,"href":"https:\/\/www.sanner-group.com\/fr\/wp-json\/wp\/v2\/blog_category?post=4803"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}