Supervised machine learning (ML) projects require data for model training, validation, and testing. However, the confidential nature of field and well production data often hinders the progress of ML projects. To address this issue, we developed a well simulator that generates realistic well production data based on physical, governing differential equations. The simulation models the reservoir, wellbore, flowline, and choke coupled using transient nodal analysis to solve for transient flow rate, pressure, and temperature as a function of variable choke opening over time in addition to a wide range of static parameters for each component. The simulator’s output is then perturbed using the gauge transfer function to introduce systematic and random errors, creating a dataset for ML projects without the need for confidential production data. We then generated a simulated dataset to train a recurrent neural network (RNN) on the task of classifying well on/off times. This task typically requires a significant number of manhours to manually filter and verify data for hundreds or thousands of wells. Our RNN model achieves high accuracy in classifying the correct on/off labels, representing a promising step towards a fully-automated rate allocation process. Our simulator for well production data can be used for other ML projects, circumventing the need for confidential data, and enabling the study and development of different ML models to streamline and automate various oil and gas work processes. Overall, the success of our RNN model demonstrates the potential of ML to improve the operational efficiency of various oil and gas work processes.

Power distribution networks (PDNs) are conducting structures employed in semiconductor systems with the aim of providing circuits with reliable and constant operating voltage. This network has non-neglible electrical parasitics. Consequently, when digital circuits inside the chip switch, the supply voltage delivered to them does not remain ideal and exhibits spatial and temporal voltage fluctuations. These fluctuations in the supply voltage, known as the power-supply noise (PSN), can affect the functionality and the performance of modern microprocessors. The design of this PDN in the chip is an important part in ensuring power integrity. Modeling and simulation of the PSN in on-chip PDNs is important to reduce the cost of processors. These PDNs have irregular geometries, which affect the PSN. As a result, they have to be modeled. The problem sizes encountered in this simulation are usually large (on the order of millions), necessitating computationally efficient simulation approaches. Existing approaches for this simulation do not guarantee at least one of the following three required properties: computationally efficiency, accuracy, and numerically robustness. Therefore, there is a need to develop accurate, numerically robust, and efficient algorithms for this simulation. For many interconnects (e.g., transmission lines, board connectors, package PDNs), only their frequency responses and SPICE circuits (e.g., nonlinear switching drivers, equivalent circuits of interconnects) terminating them are known. These frequency responses are usually available only up to a certain maximum frequency. Simulating the electrical behavior of these systems is important for the reliable design of microprocessors and for their faster time-to-market. Because terminations can be nonlinear, a transient simulation is required. There is a need for a transient simulation of band-limited frequency-domain data characterizing a multiport passive system with SPICE circuits. The number of ports can be large (greater than or equal to 100 ports). In this simulation, unlike in traditional circuit simulators, normal properties like stability and causality of transient results are not automatically met and have to be ensured. Existing techniques for this simulation do not guarantee at least one of the following three required properties: computationally efficiency for a large number of ports, causality, and accuracy. Therefore, there is a need to develop accurate and efficient time-domain techniques for this simulation that also ensure causality. The objectives of this Ph.D. research are twofold: 1) To develop accurate, numerically robust, and computationally efficient time-domain algorithms to compute PSN in on-chip PDNs with irregular geometries. 2) To develop accurate and computationally efficient time-domain algorithms for the causal cosimulation of band-limited frequency-domain data with SPICE circuits.

Zur numerischen L\"osung von Algebro-Differentialgleichungen (ADGln) m\"ussen konsistente Anfangswerte berechnet werden. Diese Arbeit befasst sich mit einem Ansatz zur Behandlung dieses Problems f\"ur Index-2 DAEs unter Verwendung von Projektoren auf die zur DAE zugeh\"origen Unterr\"aume. Die Arbeit hat zwei Schwerpunkte.\\ Zum einen werden neue Struktureigenschaften aus schwachen Voraussetzungen hergeleitet. Anschlie{\ss}end wird eine Vorgehensweise zur Auswahl von geeigneten Gleichungen einer Index-2 ADGln vorgeschlagen, deren Differentiation zu einer Indexreduktion f\"uhrt. Diese Indexreduktion liefert neue Existenz- und Eindeutigkeitsergebnisse f\"ur L\"osungen von Index-2 ADGln. Die Ergebnisse umfassen eine allgemeinere Aufgabenklasse als die bisherigen Resultate. Beruhend auf dieser Vorgehensweise wird ein stufenweiser Ansatz zur Berechnung konsistenter Anfangswerte hergeleitet. Auf diese Weise werden neue Einsichten hinsichtlich der Ausnutzung von Struktureigenschaften von Index-2 ADGln gewonnen. Insbesondere stellt sich heraus, dass im Vergleich zu Index-1 ADGln der zus\"atzliche Schritt oft in der L\"osung eines linearen Systems besteht. Die sich hieraus ergebenden numerischen Folgen werden f\"ur zwei in der Schaltungssimulation h\"aufig verwendete Verfahren, das implizite Eulerverfahren und die Trapezregel, erl\"autert. \\ Zum anderen wird die Anwendung der erhaltenen Ergebnisse auf die Gleichungen, die bei der Schaltungssimulation mittels modifizierter Knotenanalyse entstehen, ausgearbeitet. Abschlie{\ss}end wird eine kurze \"Ubersicht der durchgef\"uhrten Umsetzung gegeben.\\ / For solving DAEs numerically, consistent initial values have to be calculated. This thesis deals with an approach for handling this problem for index-2 DAEs by considering projectors onto the spaces related to the DAE. There are two major aspects in this work.\\ On the one hand, new structural properties are deduced from weak assumptions. Subsequently, a method is proposed to choose suitable equations of an index-2 DAE, whose differentiation leads to an index reduction. This index reduction yields new theoretical results for the existence and uniqueness of solutions of index-2 DAEs which apply to a wider class of applications than previous results. Based on this method, a step-by-step approach to compute consistent initial values is developed. In this way, we gain new insights about how to deal with structural properties of index-2 DAEs. In particular, it turns out that, in comparison to index-1 DAEs, the additional step that has to be undertaken in practice often consists in solving a linear system. The numerical consequences of this fact are exemplified for two methods commonly used in circuit simulation, the implicit Euler method and the trapezoidal rule.\\ On the other hand, the application of the obtained results to the equations arising in circuit simulation by means of the modified nodal analysis (MNA) is worked out. Finally, a short overview of the specifics of their realization is given.

Avec l’accroissement des exigences des avionneurs, les problématiques thermique interviennent de plus en plus tôt dans les phases de conception d’un équipement roue et frein. La première phase de conception, appelée pré-étude, doit permettre de définir les grands choix technologiques pour répondre aux spécifications. Cette phase se déroulant sur une période de l’ordre du mois, la prévision du comportement thermique induit par les différentes solutions doit être de l’ordre de quelques jours. Aussi, les outils de simulation (de type CFD) ne peuvent pas être utilisé, car le temps d’implémentation (maillage) et de calcul ne sont pas adaptés. En outre, compte tenu de la faible maturité du design, aucune CAO n’est disponible. En revanche, le degré de précision et de résolution exigé sur les résultats est plus faible. Dans cette thèse nous proposons un outil de modélisation simplifiée adapté aux contraintes de cette phase de conception, à savoir, un temps d’implémentation et de calcul de l’ordre de quelques heures, et, un très fort niveau de paramétrisation. Pour cela, nous avons tout d’abord établi une représentation géométrique simplifié de l’équipement dans une grille en coordonnées R,T,Z. La simplification géométrique se traduit par une adaptation de la forme des pièces à la grille. L’obtention de cette représentation devient alors automatisable pour réduire le temps d’implémentation. Le modèle thermique utilisé est basé sur l’association d’une approche nodale par l’utilisation de formulation RC et un maillage spatiale adaptable sur chaque pièce qui se rapproche du principe d’une méthode de type volume finis. La complexité du problème thermique (liée à la géométrie, aux transferts multi-modes instationnaires, aux hautes températures et aux non-linéarités) va nous imposer, à des fins d’analyses, la mise au point d’un outil de représentation spatiale et temporelle des grandeurs thermiques. Une décomposition géométrique, associée à une décomposition temporelle du refroidissement par rapport à des temps caractéristiques, a permis d’établir une représentation intelligible des informations sous formes de cartographies. Dans une première étape, des cartographies de températures ont été établies à partir de données d’essais. Ces cartographies ont permis d’énoncer quelques hypothèses sur le sens et l’intensité des transferts, en particulier, conductifs. Dans une seconde étape, les résultats d’un modèle thermo-aéraulique complet de l’équipement, recalés sur les températures d’essais, ont offert la possibilité d’obtenir des cartographies énergétiques par l’intermédiaire des flux et des décompositions spatiale et temporelle. Ces dernières mettent en évidence les lieux et l’intensité des transferts énergétiques, ainsi que, les pièces stockant la chaleur. Cet outil de représentation a pu être utilisé sur d’autres équipements roue et frein, et sur plusieurs énergies de freinage (16 à 64MJ). / For several years, expectations of aircraft manufacturers and airlines have risen as far as reduction of aircraft weight and turn-around time are concerned. Mass reduced wheel and brake equipments tend to reach higher maximal temperature and need to be cooled down more efficiently. Therefore, equipments need to be better designed in terms of heat transfer, especially from the early design phases (i.e. right from the request for proposal). Brake cooling is insured by transient and multi-modal heat transfer, involving particularly, surface-to-surface thermal radiation, natural convection and heat conduction. Furthermore, because of the geometric complexity of the system, three-dimensional heat transfer should be considered. The first aim of this work is to develop a simplified representation of the heat transfer, inside the equipment and with the outdoor, leading to a better acknowledgment of its thermal comportment. This would eventually provide to the engineers some guidelines to a way to help the decision forecast and improve the thermal design of wheel and brake equipments. Secondly, a simplified thermal model should be developed, in order to fill the lack of dedicated tools ,in the preliminary design stage, to predict the thermal comportment. To predict thermal fields, in such a complex system, modelisation approaches, such as finite volume or finite-element methods (the use of CFD in the case of coupled convection/conduction problems) are commonly used in industry. Despite their ability to provide suitable and accurate results, that kind of numerical methods don't fit with the preliminary design stage requirements, especially in term of implementation and calculation time. On the opposite, model based on RC thermal network seems to meet the requirements of fast calculation time and implementation time thanks to fully adapted parameterization.

Carbon-based conductors like carbon nanotubes (CNTs) and graphene nanoribbons (GNRs) have many properties, which make them relevant for potential electronic applications. Among them are high conductances and tunable band gap sizes. These properties make CNTs and GNRs useful in many circumstances, e.g. as channel material in transistors or transparent electrodes in solar cells. Plenty of literature can be found on the topic of single linear CNTs/GNRs. Some applications however require a large network of these conductors. In addition, a single conductor has only a small impact on the network conductance, which reduces the need to control the properties of each individual nanotube/-ribbon. This leads to networks being easier to apply. In this work, the conductance of large networks of GNRs is calculated using the quantum-transport formalism (QT). This has not been done before in literature. In order to apply QT to such a large amount of atoms, the recursive Green's function formalism is used. For this the networks are devided into subcells, which are represented by tight-binding matrices. Similar networks are also examined using two different nodal analysis (NA) approaches, where the nanoribbons are treated as ohmic conductors. For NA with one-dimensional conductors, major discrepancies are found in regards to the QT model. However, networks consisting of two-dimensional conductors (NA-2D) have many properties similar to the QT networks. A recipe to approximate the QT results with NA-2D is presented.:1. Introduction 2. Theoretical Principles 2.1 Carbon-based Conductors 2.1.1 Structure and Properties 2.1.2 Networks 2.2 Tight-Binding Model 2.3 Quantum Transport 2.3.1 Introduction 2.3.2 Level Broadening 2.3.3 Current Flow 2.3.4 Transmission 2.4 Nodal Analysis 3. Implementation 3.1 Quantum Tranport 3.1.1 Network Generation 3.1.2 Density-Functional based Tight-Binding Method 3.1.3 Recursive Green's Function Algorithm 3.1.4 Conductance 3.2 Nodal Analysis 3.2.1 One-dimensional Conductors 3.2.2 Two-dimensional Conductors 4. Results 4.1 Quantum Transport 4.1.1 Band Structures and Fermi Energies 4.1.2 Ideal Transmission and Consistency Tests 4.1.3 Percolation 4.1.4 Transmission 4.1.5 Conductance 4.1.6 Power Law Scaling 4.1.7 Size Dependence and Confinement Effects 4.1.8 Calculation Time 4.2 Nodal Analysis 4.2.1 One-dimensional Conductors 4.2.2 Two-dimensional Conductors 4.2.3 Calculation Time 4.3 Approximating QT with NA 4.3.1 Optimal Parameters 4.3.2 Percolation 4.3.3 Conductance 4.3.4 Power Law Scaling 5. Conclusions / Graphenbasierte Leiter wie Kohlenstoff-Nanoröhrchen (engl. 'carbon nanotubes', CNTs) oder Graphen-Nanobänder (engl. 'graphene nanoribbons', GNRs) haben viele Eigenschaften, die sie für potenzielle elektronische Anwendungen interessant machen. Darunter sind hohe Leitfähigkeiten und einstellbare Bandlückengrößen. Dadurch sind CNTs und GNRs in vielen Bereichen nützlich, z.B. als Kanalmaterial in Transistoren oder als transparente Elektroden in Solarzellen. Es gibt viel Literatur über einzelne, lineare CNTs/GNRs. Einige Anwendungen benötigen jedoch ein großes Netzwerk dieser Leiter. Zusätzlich hat ein einzelner Leiter wenig Einfluss auf die Leitfähigkeit des Netzwerks, wodurch die Eigenschaften der einzelnen Nanoröhrchen/-streifen weniger streng kontrolliert werden müssen. Dies führt dazu, dass es einfacher ist Netzwerke zu nutzen. In dieser Arbeit wird die Leitfähigkeit von großen GNR-Netzwerken mittels Quantentransport (QT) berechnet. Dies wurde in der Literatur noch nicht getan. Um QT auf eine so große Menge an Atomen anzuwenden wird der rekursive Greenfunktions-Formalismus benutzt. Dazu werden die Netzwerke in Unterzellen unterteilt, die durch Tight-Binding-Matrizen dargestellt werden. Ähnliche Netzwerke werden auch mit zwei Versionen der Knotenanalyse (engl. 'nodal analysis', NA) untersucht, welche die Nanobänder wie ohmische Leiter behandelt. Die Ergebnisse der NA mit eindimensionalen Leitern weisen deutliche Unterschiede zu den mit QT erzielten Ergebnissen auf. Wenn jedoch zweidimensionale Leiter in NA verwendet werden (NA-2D) gibt es viele parallelen zu den QT Ergebnissen. Zuletzt wird ein Vorgehen präsentiert, mit dem QT Resultate durch NA-2D Rechnungen genähert werden können.:1. Introduction 2. Theoretical Principles 2.1 Carbon-based Conductors 2.1.1 Structure and Properties 2.1.2 Networks 2.2 Tight-Binding Model 2.3 Quantum Transport 2.3.1 Introduction 2.3.2 Level Broadening 2.3.3 Current Flow 2.3.4 Transmission 2.4 Nodal Analysis 3. Implementation 3.1 Quantum Tranport 3.1.1 Network Generation 3.1.2 Density-Functional based Tight-Binding Method 3.1.3 Recursive Green's Function Algorithm 3.1.4 Conductance 3.2 Nodal Analysis 3.2.1 One-dimensional Conductors 3.2.2 Two-dimensional Conductors 4. Results 4.1 Quantum Transport 4.1.1 Band Structures and Fermi Energies 4.1.2 Ideal Transmission and Consistency Tests 4.1.3 Percolation 4.1.4 Transmission 4.1.5 Conductance 4.1.6 Power Law Scaling 4.1.7 Size Dependence and Confinement Effects 4.1.8 Calculation Time 4.2 Nodal Analysis 4.2.1 One-dimensional Conductors 4.2.2 Two-dimensional Conductors 4.2.3 Calculation Time 4.3 Approximating QT with NA 4.3.1 Optimal Parameters 4.3.2 Percolation 4.3.3 Conductance 4.3.4 Power Law Scaling 5. Conclusions

No description available.

Diese Arbeit widmet sich der Simulation von elektrischen/elektronischen Schaltungen welche um elektromagnetische Bauelemente erweitert werden. Im Fokus stehen unterschiedliche Kopplungen der Schaltungsgleichungen, modelliert mit der modifizierten Knotenanalyse, und den elektromagnetischen Bauelementen mit deren verfeinerten Modell basierend auf den vollen Maxwell-Gleichungen in der Lorenz-geeichten A-V Formulierung welche durch Finite-Integrations-Technik räumlich diskretisiert werden. Eine numerische Analyse erweitert die topologischen Kriterien für den Index der resultierenden differential-algebraischen Gleichungen, wie sie bereits in anderen Arbeiten mit ähnlichen Feld/Schaltkreis-Kopplungen hergeleitet wurden. Für die Simulation werden sowohl ein monolithischer Ansatz als auch Waveform-Relaxationsmethoden untersucht. Im Mittelpunkt stehen dabei Zeitintegration, Skalierungsmethoden, strukturelle Eigenschaften und ein hybride Ansatz zur Lösung der zugrundeliegenden linearen Gleichungssysteme welcher den Einsatz spezialisierter Löser für die jeweiligen Teilsysteme erlaubt. Da die vollen Maxwell-Gleichungen zusätzliche Ableitungen in der Kopplungsstruktur verursachen, sind bisher existierende Konvergenzaussagen für die Waveform-Relaxation von gekoppelten differential-algebraischen Gleichungen nicht anwendbar und motivieren eine neue Konvergenzanalyse. Auf dieser Analyse aufbauend werden hinreichende topologische Kriterien entwickelt, welche eine Konvergenz von Gauß-Seidel- und Jacobi-artigen Waveform-Relaxationen für die gekoppelten Systeme garantieren. Schließlich werden numerische Benchmarks zur Verfügung gestellt, um die eingeführten Methoden und Theoreme dieser Abhandlung zu unterstützen. / This work is devoted to the simulation of electrical/electronic circuits incorporating electromagnetic devices. The focus is on different couplings of the circuit equations, modeled with the modified nodal analysis, and the electromagnetic devices with their refined model based on full-wave Maxwell's equations in Lorenz gauged A-V formulation which are spatially discretized by the finite integration technique. A numerical analysis extends the topological criteria for the index of the resulting differential-algebraic equations, as already derived in other works with similar field/circuit couplings. For the simulation, both a monolithic approach and waveform relaxation methods are investigated. The focus is on time integration, scaling methods, structural properties and a hybrid approach to solve the underlying linear systems of equations with the use of specialized solvers for the respective subsystems. Since the full-Maxwell approach causes additional derivatives in the coupling structure, previously existing convergence statements for the waveform relaxation of coupled differential-algebraic equations are not applicable and motivate a new convergence analysis. Based on this analysis, sufficient topological criteria are developed which guarantee convergence of Gauss-Seidel and Jacobi type waveform relaxation schemes for introduced coupled systems. Finally, numerical benchmarks are provided to support the introduced methods and theorems of this treatise.

Die Hauptthemen dieser Arbeit sind einerseits eine tiefgehende Analyse von nichtlinearen differential-algebraischen Gleichungen (DAEs) vom Index 2, die aus der modifizierten Knotenanalyse (MNA) von elektrischen Schaltkreisen hervorgehen, und andererseits die Entwicklung von Konvergenzkriterien für Waveform Relaxationsmethoden zum Lösen gekoppelter Probleme. Ein Schwerpunkt in beiden genannten Themen ist die Beziehung zwischen der Topologie eines Schaltkreises und mathematischen Eigenschaften der zugehörigen DAE. Der Analyse-Teil umfasst eine detaillierte Beschreibung einer Normalform für Schaltkreis DAEs vom Index 2 und Abschätzungen, die für die Sensitivität des Schaltkreises bezüglich seiner Input-Quellen folgen. Es wird gezeigt, wie diese Abschätzungen wesentlich von der topologischen Position der Input-Quellen im Schaltkreis abhängen. Die zunehmend komplexen Schaltkreise in technologischen Geräten erfordern oftmals eine Modellierung als gekoppeltes System. Waveform relaxation (WR) empfiehlt sich zur Lösung solch gekoppelter Probleme, da sie auf die Subprobleme angepasste Lösungsmethoden und Schrittweiten ermöglicht. Es ist bekannt, dass WR zwar bei Anwendung auf gewöhnliche Differentialgleichungen konvergiert, falls diese eine Lipschitz-Bedingung erfüllen, selbiges jedoch bei DAEs nicht ohne Hinzunahme eines Kontraktivitätskriteriums sichergestellt werden kann. Wir beschreiben allgemeine Konvergenzkriterien für WR auf DAEs vom Index 2. Für den Fall von Schaltkreisen, die entweder mit anderen Schaltkreisen oder mit elektromagnetischen Feldern verkoppelt sind, leiten wir außerdem hinreichende topologische Konvergenzkriterien her, die anhand von Beispielen veranschaulicht werden. Weiterhin werden die Konvergenzraten des Jacobi WR Verfahrens und des Gauss-Seidel WR Verfahrens verglichen. Simulationen von einfachen Beispielsystemen zeigen drastische Unterschiede des WR-Konvergenzverhaltens, abhängig davon, ob die Konvergenzbedingungen erfüllt sind oder nicht. / The main topics of this thesis are firstly a thorough analysis of nonlinear differential-algebraic equations (DAEs) of index 2 which arise from the modified nodal analysis (MNA) for electrical circuits and secondly the derivation of convergence criteria for waveform relaxation (WR) methods on coupled problems. In both topics, a particular focus is put on the relations between a circuit's topology and the mathematical properties of the corresponding DAE. The analysis encompasses a detailed description of a normal form for circuit DAEs of index 2 and consequences for the sensitivity of the circuit with respect to its input source terms. More precisely, we provide bounds which describe how strongly changes in the input sources of the circuit affect its behaviour. Crucial constants in these bounds are determined in terms of the topological position of the input sources in the circuit. The increasingly complex electrical circuits in technological devices often call for coupled systems modelling. Allowing for each subsystem to be solved by dedicated numerical solvers and time scales, WR is an adequate method in this setting. It is well-known that while WR converges on ordinary differential equations if a Lipschitz condition is satisfied, an additional convergence criterion is required to guarantee convergence on DAEs. We present general convergence criteria for WR on higher index DAEs. Furthermore, based on our results of the analysis part, we derive topological convergence criteria for coupled circuit/circuit problems and field/circuit problems. Examples illustrate how to practically check if the criteria are satisfied. If a sufficient convergence criterion holds, we specify at which rate of convergence the Jacobi and Gauss-Seidel WR methods converge. Simulations of simple benchmark systems illustrate the drastically different convergence behaviour of WR depending on whether or not the circuit topological convergence conditions are satisfied.