<p>Autoinjectors are pen-like devices that typically deliver drug products of 2 mL or less. They shield the needle before and after use, reducing patient anxiety from needle phobia and mitigating the risk of needlestick injuries and accidental contamination. Additionally, automatic delivery ensures more consistent needle penetration depth and injection force than manual injection methods. </p> <p><br></p> <p>To optimize autoinjector design, this thesis presents experimentally validated computational models that describe the processes of needle insertion, drug delivery, and transport of subcutaneously administered therapeutic proteins in the body. A multi-objective optimization framework is also proposed to guide the design of autoinjectors.</p> <p><br></p> <p>This thesis focuses on spring-driven autoinjectors, the most common type of autoinjector. It begins with an overview of the interactions between the spring-driven autoinjector, tissue, and therapeutic proteins. Moving on to Chapter 2, a computational model is presented to accurately predict the kinematics of the syringe barrel and plunger during the needle insertion process.</p> <p><br></p> <p>In Chapter 3, we present a quasi-steady model for the drug delivery process, which considers the rheology of therapeutic proteins. The Carreau model is adopted to describe protein viscosity, and explicit relationships between flow rate and pressure drop in the needle are derived. Furthermore, the applicable regime for the power-law model for protein viscosity is identified.</p> <p><br></p> <p>Chapter 4 quantifies the impact of sloshing and cavitation on therapeutic proteins in the syringe. Additionally, a workflow is presented to integrate available simulation tools to predict the performance of spring-driven autoinjectors. The influence of each design parameter of spring-driven autoinjectors on their performance is also discussed. </p> <p><br></p> <p>The spring-driven autoinjector delivers therapeutic proteins through subcutaneous administration. To gain insights into the transport process of therapeutic proteins, Chapter 5 presents a physiologically-based pharmacokinetic model that has been validated against experimental data for humans and rats. The lymph flow rate significantly affects the bioavailability of therapeutic proteins. This finding highlights the importance of studying the transport of therapeutic proteins in the lymphatic system in future research.</p> <p><br></p> <p>Chapter 6 provides a multi-objective design optimization framework for the spring-driven autoinjector. The computational model is replaced with an accurate deep neural network surrogate to improve the computational efficiency.  Using this surrogate model, we conduct a sensitivity analysis to identify essential design parameters. After that, we perform multi-objective optimization to find promising design candidates.</p> <p><br></p> <p>Chapter 7 presents a model for bubble dynamics in a protein solution. An explicit expression for the bubble dissolution rate is derived, enabling extraction of the interfacial properties of the protein-coated interface from the measured bubble radii. Moreover, analytical solutions for the response of a protein-coated bubble to an imposed acoustic pressure are derived. This work provides insight into protein-coated bubbles, which are used as vehicles to deliver drugs, as active miniature tracers to probe the rheology of soft and biological materials, or as contrast agents to enhance the ultrasound backscatter in ultrasonic imaging.</p> <p><br></p> <p>At last, in Chapter 8, we introduce a model for laser-induced cavitation that considers several key factors, such as liquid compressibility, heat transfer, and non-equilibrium evaporation and condensation. Our model's predictions for the time-course of bubble radii have been validated with experimental data. Moreover, our model reveals that the reduction of the bubble's oscillation amplitude is primarily due to a decrease in the number of vapor molecules inside the bubble, highlighting the crucial role of phase change in laser-induced cavitation bubbles.</p> <p><br></p> <p>The developed computational models and framework provide crucial insights into the development of spring-driven autoinjectors and cavitation bubbles. These studies can also enhance the efficacy and safety of the delivery of therapeutic proteins, ultimately improving patient outcomes.</p>

[pt] A adsorção de moléculas de asfalteno na interface óleo-água induz a formação de uma microestrutura complexa, que estabiliza as emulsões e prejudica a eficiência dos processes de refino de petróleo. Neste trabalho, desenvolvemos um conjunto de novos protocolos de reologia de cisalhamento para avaliar o efeito de solventes polares e apolares na adsorção de genuínos asfaltenos brasileiros. Além disso, a morfologia do asfalteno, após a adição de solventes com aromaticidades distintas, é investigada por microscopia de varredura (MEV). Os resultados indicam que os asfaltenos estão organizados em uma estrutura do tipo ilha com unidades aromáticas e policondensadas, que formam filmes interfaciais reversíveis com a adição de solventes polares. O estudo interfacial também revela que solventes apolares podem "prender" os nanoagregados de asfalteno na mistura. Esses agregados, na presença de solventes fracamente polares, podem se consolidar em um padrão mais compactado, sugerindo que o crescimento do filme e o autoarranjo do asfalteno estão diretamente relacionados ao conteúdo aromático. Explora-se as diferenças na estruturação do asfalteno e como afetam a extensão da emulsificação espontânea. É proposto que a taxa de emulsificação está diretamente relacionada à configuração química dos asfaltenos. Finalmente, estuda-se a adição de ácido esteárico (AE) a soluções de asfalteno em conteúdo de água deionizada (AD) e água sintética (AS) para melhor compreender como as propriedades reológicas e superficiais são afetadas pela competição das coespécies. Verifica-se que interfaces formadas puramente por AEs originam filmes mais viscosos do que elásticos na interface ar-água. A atividade interfacial dos asfaltenos brasileiros é evidente e significativa na presença de eletrólitos e dependente da aromaticidade do solvente. / [en] Adsorption of asphaltene molecules at the oil-water interface induces the formation of a complex microstructure, which stabilizes emulsions and impairs the efficiency of crude oil refining. In this work, we design a set of new shear rheology protocols to assess the effect of polar and non-polar solvents on indigenous Brazilian (BR) asphaltene adsorption. Moreover, the asphaltene morphology upon addition of solvents with distinct aromaticities is investigated by SEM microscopy. Our findings indicate that asphaltenes are a polycondensate aromatic island-type structure that forms reversible films when polar solvents are placed on top of the adsorbed film. The interfacial study also reveals that non-polar solvents may lock up asphaltene nanoaggregates in mixture. These aggregates, upon the presence of weakly polar solvents, can consolidate into a more close-packed pattern, suggesting that network growth and asphaltene self-arrangement are directly related to the aromatic content. We explore the differences in asphaltene structuring and how it affects the extent of spontaneous emulsification. We find that the rate of emulsification is directly related to the chemical configuration of asphaltenes. Finally, we study the addition of stearic acid (SA) to asphaltene solutions in deionized water (DW) and synthetic water (SW) to better understand how surface and rheological properties are affected by competitive adsorption. We find that single SA are more prone to form liquid-like rather than solid-like films at the air-water interface. Furthermore, we show that the interfacial activity of our asphaltenes is enhanced in the presence of electrolytes and is dependent of the solvent aromaticity.

[pt] Inúmeros estudos têm sido realizados para melhor entender a formação e estabilidade de emulsões. Em algumas situações, é desejável ter emulsões estáveis; em outros, a separação de fases por coalescência de gotas é benéfica. Em ambos os casos, é importante entender os mecanismos associados ao processo de coalescência. O presente trabalho investiga a relação entre as propriedades reológicas das interfaces óleo-água e o tempo de drenagem de um filme fino de óleo entre duas gotas aquosas. A tensão interfacial e a reologia dilatacional interfacial foram medidas usando a análise axissimétrica da forma da gota. Foram avaliadas diferentes concentrações de um tensoativo não iônico (Span 80) dissolvido em óleo mineral (Primol 352). Os resultados indicam uma relação direta entre as propriedades da estrutura formada na interface óleo-água e a ausência de coalescência das gotas. Para concentrações de surfactante abaixo da concentração micelar crítica (CMC), a interface é fracamente elástica (fluid-like) e o processo de coalescência sempre ocorre; o tempo de drenagem não está relacionado ao tempo de envelhecimento da interface. Para concentrações de surfactante acima da CMC, os módulos elástico e viscoso mostraram mudanças significativas com o envelhecimento, levando à formação de um filme sólido na interface, impedindo a coalescência entre as gotas. Usamos experimentos de coalescência gota/gota para avaliar o efeito da reologia interfacial na dinâmica de coalescência. Para entender melhor o fenômeno da não coalescência, estudamos microscopicamente a estrutura do filme interfacial e observamos o aparecimento de pequenas gotas de água formadas na interface através de emulsificação espontânea. Descobrimos que a taxa de surgimento dessas microgotículas está diretamente relacionada à concentração de surfactante. À medida que a concentração de surfactante aumenta, mais rápido ocorre o processo de emulsificação espontânea, o que confirma os resultados obtidos com a reologia interfacial. Finalmente, um novo método para promover a desestabilização da emulsão impondo uma perturbação do filme interfacial pelo escoamento das gotas através de capilares constritos é proposto e testado. / [en] Several studies have been conducted to understand emulsions formation and stability. In some situations, it is desirable to have stable emulsions; in others, phase separation through drop coalescence is beneficial. In both cases, it is important to understand the mechanisms associated to the coalescence process. The present work investigates the relationship between rheological properties of oil-water interfaces and the drainage time of a thin oil film between two aqueous drops. Interfacial tension and dilatational rheology were measured using the axisymmetric drop shape analysis. We evaluated different concentrations of a nonionic surfactant (Span 80) dissolved in mineral oil (Primol 352) phase. The results indicate a direct relationship between the properties of the structure formed at the oil-water interface and the absence of droplet coalescence. For low surfactant concentrations, below the critical micelle concentration (CMC), the interface is weakly elastic (fluid-like) and the coalescence process always occurs; the draining time is not to related to the aging time of the interface. For surfactant concentrations above CMC, the elastic and viscous moduli showed significant changes with aging leading to the formation of a solid-like film at the interface preventing further coalescence. We used a drop/drop coalescence experiments to evaluate the effect of interfacial rheology on the coalescence dynamics. To better understand the phenomenon of non-coalescence, we study the structure of interfacial film microscopically and observe the appearance of small water droplets formed at the interface by spontaneous emulsification. We found that the emergence rate of these microdroplets is directly related to the surfactant concentration. As the surfactant concentration increases, faster the spontaneous emulsification process occurs, which confirms the results obtained with the interfacial rheology. Finally, a new method to promote emulsion destabilization by imposing a perturbation of the interfacial film by flowing the drops through constricted capillaries is proposed and tested.

L’objectif de cette étude a été de formuler et caractériser des émulsions simples huile/eau d’intérêt pharmaceutique stabilisées par de la β-lactoglobuline (β-lg), de la gomme arabique (GA), de la gomme xanthane (GX) et des mélanges β-lg:GA et β-lg:GX. Les concentrations massiques totales des dispersions de biopolymères étaient de 1 % et ont été augmentées à 2,5 % si les émulsions formulées n’étaient pas stables. Le mélange β-lg:GA a été réalisé à pH 4,2 afin de permettre la formation de complexes par interactions électrostatiques attractives entre la β-lg et la GA. Deux ratios β-lg:GA ont été étudiés : 2:1 et 1:2. Enfin, le mélange β-lg:GX a été effectué à pH 7, où les deux biopolymères étant chargés négativement ne se complexent pas et à un ratio de 1:1. Une étude de stabilité des émulsions a été menée sur 6 mois. Les stabilités obtenues ont pu être classées par ordre croissant : GA 2,5 % < β-lg:GA 2,5 % < β-lg 2,5 % < GX 1 % = β-lg:GX 1 %. Plusieurs mécanismes de stabilisation ont été mis en évidence grâce à l’étude des propriétés interfaciales des biopolymères, à l’étude des propriétés rhéologiques des émulsions et à des observations au microscope confocal à balayage laser des émulsions après marquage des biopolymères à la fluorescence. La β-lg et la GA sont toutes deux capables de s’adsorber à l’interface des globules huileux alors que la GX augmente la viscosité de la phase continue. L’association β-lg:GA conduit à la formation d’une double couche interfaciale stabilisante. Enfin, l’association β-lg:GX combine les mécanismes de stabilisation de la protéine, par adsorption interfaciale et de la gomme, par augmentation de la viscosité de la phase continue. / The main objective of this study was to formulate and characterize oil-in-water simple emulsions of pharmaceutical interest stabilized by β-lactoglobulin (β-lg), gum arabic (GA), xanthan gum (XG), and mixtures of β-lg:GA and β-lg:XG. The total biopolymer final concentration in the dispersions was 1 (w/w) % and could be raised to 2.5 (w/w) % if the formulated emulsions were not stable. β-lg:GA mixing was performed at pH 4.2 to allow attractive electrostatic interactions between the two biopolymers and thus the formation of complexes. Two protein:polysaccharide ratios were investigated: 2:1 and 1:2. Conversely, β8lg:XG mixing was performed at pH 7, where both biopolymers are negatively charged, in order to avoid the complex formation, and with a 1:1 ratio. A stability study was conducted for emulsions over a 6-month period. The obtained stabilities could be classified increasingly: GA 2.5 % < β-lg:GA 2.5 % < β-lg 2.5 % < XG 1 % = β-lg:XG 1 %. Several stabilization mechanisms were evidenced by the study of the biopolymer interfacial properties, the study of emulsion rheology and by confocal laser scanning microscopy observations with labeled fluorescent biopolymers. β-lg and GA were both able to adsorb at the interface of oil globule. XG enhanced the continuous phase viscosity. β-lg:GA mixing led to the formation of a stabilizing interfacial double layer. Finally, β-lg:XG association combined the stabilization mechanisms of both biopolymers, respectively: interfacial adsorption and enhancement of the continuous phase viscosity.

Les travaux de cette thèse concernent des études fondamentales liées à la rhéologie interfaciale des systèmes polymères multicouches. Les matériaux choisis sont à base de deux polymères compatibles,PVDF et PMMA de différentes masses molaires. Ces systèmes ont été étudiés sous sollicitations en cisaillement et en élongation suivant les deux régimes en viscoélasticité linéaire (VEL) et non-linéaire (VENL). Les études en VEL ont permis d’étudier la cinétique de développement de l’interphase. Quant aux études en VENL, elles ont permis d’étudier les propriétés intrinsèques de l’interphase simulant ainsi les conditions de mise en œuvre proches de celles des procédés usuels. On démontre ainsi que la rhéologie joue le rôle d’une sonde très fine pour explorer les propriétés aux interfaces des matériaux multicouches. Des modélisations ont été établies en se basant sur les mécanismes physiques mis en jeu. Dans un premier temps, le comportement rhéologique à l’état fondu des multicouches a été étudié par spectrométrie mécanique dynamique en VEL. Les cinétiques d'interdiffusion ainsi que le développement de l’interphase générée aux interfaces de bicouches symétriques et asymétriques ont été étudiés. Les résultats obtenus ont été analysés et modélisés en se basant sur les concepts de la dynamique moléculaire en l’occurrence le modèle de Doi et Edwards. De plus, un nouveau modèle rhéologique a été développé. Il a permis de quantifier les coefficients d'interdiffusion. Les coefficients de friction des chaines et les propriétés rhéologiques de l’interphase ont été modélisés à leur tour. Les résultats obtenus corroborent ceux de la littérature, obtenus par des méthodes spectroscopiques sophistiquées. Le modèle a permis de quantifier les grandeurs viscoélastiques et l’épaisseur de l’interphase. Dans un second temps, des expérimentations en VENL ont été réalisées. Un modèle original a été également proposé pour décrire le comportement relatif à la relaxation des multicouches et de l'interphase. De plus, la sensibilité de la densité d’enchevêtrement a été étudiée pendant et après sollicitations. On démontre que sa présence retarde l'écoulement interfacial surtout sous hautes déformations et vitesses de déformation. En outre, les études des structures multicouches sous sollicitation élongationnelle ont montrées que les propriétés dépendent du rapport de viscosité des couches et les propriétés de l’interphase diffuse. Les travaux de cette thèse mettent en lumière la compétition entre l’effet négatif de l'orientation des chaînes et l'effet favorable de l'écoulement sur les cinétiques de la diffusion. Ensuite, des cartes de stabilités des écoulements stratifiés ont été établies.. La présence de l'interphase diffuse a contribué à une élimination des instabilités. On montre ainsi qu’outre la cinématique de l’écoulement en cisaillement et en élongation, les propriétés de l’interphase ont un rôle important dans la stabilité des écoulements stratifiés en coextrusion. / Fundamental studies have been devoted in this work to probe and modelize the interfacial phenomena at multilayered polymer systems based on two model compatible polymers of PVDF and PMMA with varying molar masses. Linear and nonlinear rheology have been demonstrated to be sensitive to the presence of diffuse interphase triggered from interdiffusion at polymer/polymer interface. Firstly, the interdiffusion kinetics as well as the development of the interphase decoupled to flow as generated at a symmetrical (self diffusion) and an asymmetrical (mutual diffusion) bilayer have been investigated using small-amplitude oscillatory shear measurements. Results were analyzed according to Doi-Edwards theory (tube model) and the effects of annealing factors as well as structural properties on the diffusion kinetics have been studied. The PMMA/PVDF mixtures have been examined to be a couple of weak thermorheological complexity, owning close monomeric friction coefficients of each species in the present experimental conditions. Based on this physics, a new rheological model was developed to quantify the interdiffusion coefficients by taking into account the component dynamics in mixed state and the concept of interfacial rheology. Rheological and geometrical properties of the interphase have been able to be quantified through this model, as validated by scanning electron microscopy coupled with energy dispersive X-ray analysis (SEM-EDX) and transmission electron microscopy (TEM). Secondly, experiments of step strain, startup in simple shear and in uni-axial extension have been carried out on the PMMA/PVDF multilayer structures. An original model was proposed to fit the stress relaxation behavior of multilayer structures and to estimate the relaxation behavior of the interphase. Lack of entanglement at the interface and weak entanglement intensity at the diffuse interphase make them to be subsequently readily to suffer from interfacial yielding even interfacial failure during and after continuous large deformations. Interphase delays the interfacial yielding to a larger external deformation or a higher deformation rate. Besides, elongational properties of the multilayer structures have been shown to be a function of composition as controlled by layer number(interfacial area) and interphase properties (rheology related to entanglement intensity). Finally, the diffuse interphase development coupled to flow in practical coextrusion process has been considered. The compromising result between negative effect of chain orientation and favorable effect of flow on diffusion kinetics gives rise to a broadening interphase after coextrusion. Presence of the diffuse interphase was demonstrated to significantly weaken (or even eliminate) the viscous and elastic instabilities despite of the high rheological contrast. Hence, this work gives guidelines on the key role of the interphase plays in structure-property-processing relationships.