<p dir="ltr">Distributed energy resources (DERs), such as rooftop solar generation and energy storage systems, are becoming more prevalent in distribution systems. DERs are connected to the distribution system via power electronic converters, introducing faster dynamics in the system. Understanding the system dynamics under a high penetration of inverter-based DERs is critical for power system researchers and practitioners, driving the development of modeling techniques and simulation software. Aiming to reduce computational complexity, existing tools and techniques often employ various approximations. Meanwhile, modern advancements in computational hardware capabilities provide opportunities to include the faster time-scale dynamics. To address this, the primary objective of this thesis is to develop an open-source Python simulation package, Dynamic Simulation using Sparse Tableau Approach in Python, DynaSTPy (pronounced “dynasty”), capable of capturing the dynamics of all components in a distribution feeder. The distribution feeder is modeled as a system of Differential-Algebraic Equations (DAEs). Further, each component in the feeder is modeled based on the Sparse Tableau Approach (STA), which involves the representation of component model equations using sparse matrices, facilitating a systematic procedure to model the components and construct the system DAEs. In sinusoidal steady state, the DAEs can be represented in phasor form, extending the approach to perform power flow analysis of distribution feeders.</p>

Nous développons dans ce travail de thèse un outil de simulation numérique pour les matériaux à changement de phase (MCP), en tenant compte du phénomène de convection naturelle dans la phase liquide, pour des configurations en deux et trois dimensions. Les équations de Navier-Stokes incompressible avec le modèle de Boussinesq pour la prise en compte des forces de flottabilité liées aux effets thermiques, couplées avec une formulation de l’équation d’énergie suivant la méthode d’enthalpie, sont résolues par une méthode d’éléments finis adaptatifs. Une approche mono-domaine, consistant à résoudre les mêmes systèmes d’équations dans les phases solide et liquide, est utilisée. La vitesse est ramenée à zéro dans la phase solide, en introduisant un terme de pénalisation dans l’équation de quantité de mouvement, suivant le modèle de Carman-Kozeny, consistant à freiner la vitesse à travers un milieu poreux. Une discrétisation spatiale des équations utilisant des éléments finis de Taylor-Hood, éléments finis P2 pour la vitesse et éléments finis P1 pour la pression, est appliquée, avec un schéma d’intégration en temps implicite d’ordre deux (GEAR). Le système d’équations non-linéaires est résolu par un algorithme de Newton. Les méthodes numériques sont implémentées avec le logiciel libre FreeFem++ (www.freefem.org), disponible pour tout système d’exploitation. Les programmes sont distribués sous forme de logiciel libre, sous la forme d’une forme de toolbox simple d’utilisation, permettant à l’utilisateur de rajouter d’autres configurations numériques pour des problèmes avecchangement de phase. Nous présentons dans ce manuscrit des cas de validation du code de calcul, en simulant des cas tests bien connus, présentés par ordre de difficulté croissant : convection naturelle de l’air, fusion d’un MCP, le cycle complet fusion-solidification, chauffage par le bas d’un MCP, et enfin, la solidification de l’eau. / In this thesis we develop a numerical simulation tool for computing two and three-dimensional liquid-solid phase-change systems involving natural convection. It consists of solving the incompressible Navier-Stokes equations with Boussinesq approximation for thermal effects combined with an enthalpy-porosity method for the phase-change modeling, using a finite elements method with mesh adaptivity. A single-domain approach is applied by solving the same set of equations over the whole domain. A Carman-Kozeny-type penalty term is added to the momentum equation to bring to zero the velocity in the solid phase through an artificial mushy region. Model equations are discretized using Galerkin triangular finite elements. Piecewise quadratic (P2) finite-elements are used for the velocity and piecewise linear (P1) for the pressure. The coupled system of equations is integrated in time using a second-order Gear scheme. Non-linearities are treated implicitly and the resulting discrete equations are solved using a Newton algorithm. The numerical method is implemented with the finite elements software FreeFem++ (www.freefem.org), available for all existing operating systems. The programs are written and distributed as an easy-to-use open-source toolbox, allowing the user to code new numerical algorithms for similar problems with phase-change. We present several validations, by simulating classical benchmark cases of increasing difficulty: natural convection of air, melting of a phase-change material, a melting-solidification cycle, a basal melting of a phase-change material, and finally, a water freezing case.