Das Be- und Entnetzen von Oberflächen durch Fluide ist in unserer Natur allgegenwärtig, spielt aber auch in vielen technischen Prozessen und Anwendungen eine entscheidende Rolle. Beispiele dafür sind das Beschichten und Drucken, Mikrofluidik und Lab-on-a-Chip-Technologien oder auch das Kühlen in bestimmten Reaktorgeometrien oder anderen Anlagen der Prozessindustrie. Im Allgemeinen können dynamische Benetzungsvorgänge als Mehrphasenströmungen dargestellt werden, die sich mathematisch mithilfe der Navier-Stokes-Gleichungen beschreiben lassen. Zusätzlich werden Sprungbedingungen benötigt, die die Strömung der verschiedenen Fluide oder Phasen über ihre Grenzflächen hinweg verbinden, an welchen zusätzlich eine Oberflächenspannung anliegt. Kommt neben den Flüssigkeits-Gas-Grenzflächen eine dynamische Kontaktlinie vor, an der Flüssigkeit und Gas die Festkörperoberfläche berühren, wird ein etwas anderer Modellierungsansatz verfolgt. Hierfür werden die Navier-Stokes-Gleichungen um eine zusätzliche Transportgleichung ergänzt, welche die Advektion der Strömung beschreibt. Die Transportgleichung wird für einen algebraischen Volume-of-Fluid-Ansatz benötigt, der zu einer Ein-Feld-Formulierung des Problems führt. Das Modell wird durch geeignete Anfangs- und Randbedingungen vervollständigt, wobei der dynamische Kontaktwinkel als Randbedingung eingeht. Es entsteht ein System partieller Differentialgleichungen, welches in die simulationsbasierte Optimierung von Benetzungsvorgängen einfließt. Die betrachteten Optimierungsprobleme gehören zu der Klasse der Optimalsteuerungsprobleme, in denen eine Zielfunktion bezüglich eines Zustandes und einer Steuerung bzw. Kontrolle optimiert wird. Für Existenz- und Eindeutigkeitsaussagen wird die Differenzierbarkeit der zugehörigen Steuerungs-Zustands-Abbildung benötigt, wobei L^p-maximale Regularität des zugrundeliegenden linearen Zweiphasenproblems erarbeitet wird. Die Differenzierbarkeit zu zeigen ist ein zentraler Bestandteil dieser Arbeit, wobei sich bei den theoretischen Betrachtungen die Kontrolle aus einem initialen Geschwindigkeitsfeld und einer über dem gesamtem Gebiet verteilten Kontrolle auf der rechten Seite der Impulsgleichung zusammensetzt. Damit wird eine Basis geschaffen, um Optimalsteuerungsprobleme von Benetzungsprozessen mit Methoden der ableitungsbasierten Optimierung zu lösen. Jedoch gestaltet sich die Lösung eines restringierten Optimierungsproblems mit partiellen Differentialgleichungen als Nebenbedingungen im Kontext von Benetzungsproblemen als schwierig und kann nicht wie üblich mithilfe der Lagrange-Funktion gelöst werden. Daher verfolgen wir einen Sensitivitätsansatz und leiten formal Sensitivitätsgleichungen für das kontinuierliche Strömungsproblem her. Zustandsgleichungen und Sensitivitätsgleichungen werden nun numerisch mit geeigneten Diskretisierungsverfahren gelöst, da eine analytische Lösung für solche Probleme bislang nicht bekannt ist. Dafür wird der aus der OpenFOAM-Bibliothek stammende zwei-Phasen Löser interFoam dahingehend angepasst, dass simultan zu den Zustandsgleichungen auch die zugehörigen Sensitivitätsgleichungen gelöst werden. Das entwickelte Verfahren wird an einem Demonstratorbespiel getestet, welches durch den Tiefdruck motiviert ist. Für ein gutes Druckergebnis ist es essenziell, überschüssige Farbe bis auf einen dünnen Film von der Druckplatte zu entfernen. Zu diesem Zweck wird ein Stahlband über die Oberfläche gezogen, welches auch als Rakel bekannt ist. Numerische Ergebnisse werden beispielhaft für ein Parameteridentifikationsproblem und die Optimierung hinsichtlich geometrischer Aspekte für den genannten Benetzungsprozess dargestellt.