-
Multiobjective Optimization of Non-Smooth PDE-Constrained Problems
Authors:
Marco Bernreuther,
Michael Dellnitz,
Bennet Gebken,
Georg Müller,
Sebastian Peitz,
Konstantin Sonntag,
Stefan Volkwein
Abstract:
Multiobjective optimization plays an increasingly important role in modern applications, where several criteria are often of equal importance. The task in multiobjective optimization and multiobjective optimal control is therefore to compute the set of optimal compromises (the Pareto set) between the conflicting objectives. The advances in algorithms and the increasing interest in Pareto-optimal s…
▽ More
Multiobjective optimization plays an increasingly important role in modern applications, where several criteria are often of equal importance. The task in multiobjective optimization and multiobjective optimal control is therefore to compute the set of optimal compromises (the Pareto set) between the conflicting objectives. The advances in algorithms and the increasing interest in Pareto-optimal solutions have led to a wide range of new applications related to optimal and feedback control - potentially with non-smoothness both on the level of the objectives or in the system dynamics. This results in new challenges such as dealing with expensive models (e.g., governed by partial differential equations (PDEs)) and developing dedicated algorithms handling the non-smoothness. Since in contrast to single-objective optimization, the Pareto set generally consists of an infinite number of solutions, the computational effort can quickly become challenging, which is particularly problematic when the objectives are costly to evaluate or when a solution has to be presented very quickly. This article gives an overview of recent developments in the field of multiobjective optimization of non-smooth PDE-constrained problems. In particular we report on the advances achieved within Project 2 "Multiobjective Optimization of Non-Smooth PDE-Constrained Problems - Switches, State Constraints and Model Order Reduction" of the DFG Priority Programm 1962 "Non-smooth and Complementarity-based Distributed Parameter Systems: Simulation and Hierarchical Optimization".
△ Less
Submitted 2 August, 2023;
originally announced August 2023.
-
On the nonmonotone linesearch for a class of infinite-dimensional nonsmooth problems
Authors:
Behzad Azmi,
Marco Bernreuther
Abstract:
This paper provides a comprehensive study of the nonmonotone forward-backward splitting (FBS) method for solving a class of nonsmooth composite problems in Hilbert spaces. The objective function is the sum of a Fréchet differentiable (not necessarily convex) function and a proper lower semicontinuous convex (not necessarily smooth) function. These problems appear, for example, frequently in the co…
▽ More
This paper provides a comprehensive study of the nonmonotone forward-backward splitting (FBS) method for solving a class of nonsmooth composite problems in Hilbert spaces. The objective function is the sum of a Fréchet differentiable (not necessarily convex) function and a proper lower semicontinuous convex (not necessarily smooth) function. These problems appear, for example, frequently in the context of optimal control of nonlinear partial differential equations (PDEs) with nonsmooth sparsity promoting cost functionals. We discuss the convergence and complexity of FBS equipped with the nonmonotone linesearch under different conditions. In particular, R-linear convergence will be derived under quadratic growth-type conditions. We also investigate the applicability of the algorithm to problems governed by PDEs. Numerical experiments are also given that justify our theoretical findings.
△ Less
Submitted 3 March, 2023;
originally announced March 2023.
-
An adaptive certified space-time reduced basis method for nonsmooth parabolic partial differential equations
Authors:
Marco Bernreuther,
Stefan Volkwein
Abstract:
In this paper, a nonsmooth semilinear parabolic partial differential equation (PDE) is considered. For a reduced basis (RB) approach, a space-time formulation is used to develop a certified a-posteriori error estimator. This error estimator is adopted to the presence of the discrete empirical interpolation method (DEIM) as approximation technique for the nonsmoothness. The separability of the esti…
▽ More
In this paper, a nonsmooth semilinear parabolic partial differential equation (PDE) is considered. For a reduced basis (RB) approach, a space-time formulation is used to develop a certified a-posteriori error estimator. This error estimator is adopted to the presence of the discrete empirical interpolation method (DEIM) as approximation technique for the nonsmoothness. The separability of the estimated error into an RB and a DEIM part then guides the development of an adaptive RB-DEIM algorithm, combining both offline phases into one. Numerical experiments show the capabilities of this novel approach in comparison with classical RB and RB-DEIM approaches.
△ Less
Submitted 28 December, 2022;
originally announced December 2022.
-
ms2: A molecular simulation tool for thermodynamic properties, new version release
Authors:
Colin W. Glass,
Steffen Reiser,
Gábor Rutkai,
Stephan Deublein,
Andreas Köster,
Gabriela Guevara Carrión,
Amer Wafai,
Martin Horsch,
Martin F. Bernreuther,
Thorsten Windmann,
Hans Hasse,
Jadran Vrabec
Abstract:
A new version release (2.0) of the molecular simulation tool ms2 [S. Deublein et al., Comput. Phys. Commun. 182 (2011) 2350] is presented. Version 2.0 of ms2 features a hybrid parallelization based on MPI and OpenMP for molecular dynamics simulation to achieve higher scalability. Furthermore, the formalism by Lustig [R. Lustig, Mol. Phys. 110 (2012) 3041] is implemented, allowing for a systematic…
▽ More
A new version release (2.0) of the molecular simulation tool ms2 [S. Deublein et al., Comput. Phys. Commun. 182 (2011) 2350] is presented. Version 2.0 of ms2 features a hybrid parallelization based on MPI and OpenMP for molecular dynamics simulation to achieve higher scalability. Furthermore, the formalism by Lustig [R. Lustig, Mol. Phys. 110 (2012) 3041] is implemented, allowing for a systematic sampling of Massieu potential derivatives in a single simulation run. Moreover, the Green-Kubo formalism is extended for the sampling of the electric conductivity and the residence time. To remove the restriction of the preceding version to electro-neutral molecules, Ewald summation is implemented to consider ionic long range interactions. Finally, the sampling of the radial distribution function is added.
△ Less
Submitted 25 July, 2015;
originally announced July 2015.
-
ls1 mardyn: The massively parallel molecular dynamics code for large systems
Authors:
Christoph Niethammer,
Stefan Becker,
Martin Bernreuther,
Martin Buchholz,
Wolfgang Eckhardt,
Alexander Heinecke,
Stephan Werth,
Hans-Joachim Bungartz,
Colin W. Glass,
Hans Hasse,
Jadran Vrabec,
Martin Horsch
Abstract:
The molecular dynamics simulation code ls1 mardyn is presented. It is a highly scalable code, optimized for massively parallel execution on supercomputing architectures, and currently holds the world record for the largest molecular simulation with over four trillion particles. It enables the application of pair potentials to length and time scales which were previously out of scope for molecular…
▽ More
The molecular dynamics simulation code ls1 mardyn is presented. It is a highly scalable code, optimized for massively parallel execution on supercomputing architectures, and currently holds the world record for the largest molecular simulation with over four trillion particles. It enables the application of pair potentials to length and time scales which were previously out of scope for molecular dynamics simulation. With an efficient dynamic load balancing scheme, it delivers high scalability even for challenging heterogeneous configurations. Presently, multi-center rigid potential models based on Lennard-Jones sites, point charges and higher-order polarities are supported. Due to its modular design, ls1 mardyn can be extended to new physical models, methods, and algorithms, allowing future users to tailor it to suit their respective needs. Possible applications include scenarios with complex geometries, e.g. for fluids at interfaces, as well as non-equilibrium molecular dynamics simulation of heat and mass transfer.
△ Less
Submitted 20 August, 2014;
originally announced August 2014.
-
Static and dynamic properties of curved vapour-liquid interfaces by massively parallel molecular dynamics simulation
Authors:
Martin T. Horsch,
Svetlana K. Miroshnichenko,
Jadran Vrabec,
Colin W. Glass,
Christoph Niethammer,
Martin F. Bernreuther,
Erich A. Müller,
George Jackson
Abstract:
Curved fluid interfaces are investigated on the nanometre length scale by molecular dynamics simulation. Thereby, droplets surrounded by a metastable vapour phase are stabilized in the canonical ensemble. Analogous simulations are conducted for cylindrical menisci separating vapour and liquid phases under confinement in planar nanopores. Regarding the emergence of nanodroplets during nucleation, a…
▽ More
Curved fluid interfaces are investigated on the nanometre length scale by molecular dynamics simulation. Thereby, droplets surrounded by a metastable vapour phase are stabilized in the canonical ensemble. Analogous simulations are conducted for cylindrical menisci separating vapour and liquid phases under confinement in planar nanopores. Regarding the emergence of nanodroplets during nucleation, a non-equilibrium phenomenon, both the non-steady dynamics of condensation processes and stationary quantities related to supersaturated vapours are considered. Results for the truncated and shifted Lennard-Jones fluid and for mixtures of quadrupolar fluids confirm the applicability of the capillarity approximation and the classical nucleation theory.
△ Less
Submitted 20 October, 2011;
originally announced October 2011.
-
Poiseuille flow of liquid methane in nanoscopic graphite channels by molecular dynamics simulation
Authors:
Martin Horsch,
Jadran Vrabec,
Martin Bernreuther,
Hans Hasse
Abstract:
MD simulations of methane confined between graphite walls with up to 4,800,000 interaction sites, i.e., carbon atoms and methane molecules, are conducted, where the channel width is varied to include both the boundary-dominated regime and the transition to the continuum regime. This proves that MD can be used today to cover the entire range of characteristic lengths for which continuum methods f…
▽ More
MD simulations of methane confined between graphite walls with up to 4,800,000 interaction sites, i.e., carbon atoms and methane molecules, are conducted, where the channel width is varied to include both the boundary-dominated regime and the transition to the continuum regime. This proves that MD can be used today to cover the entire range of characteristic lengths for which continuum methods fail.
△ Less
Submitted 29 November, 2009;
originally announced November 2009.
-
Homogeneous nucleation in supersaturated vapors of methane, ethane, and carbon dioxide predicted by brute force molecular dynamics
Authors:
Martin Horsch,
Jadran Vrabec,
Martin Bernreuther,
Sebastian Grottel,
Guido Reina,
Andrea Wix,
Karlheinz Schaber,
Hans Hasse
Abstract:
Molecular dynamics (MD) simulation is applied to the condensation process of supersaturated vapors of methane, ethane, and carbon dioxide. Simulations of systems with up to a million particles were conducted with a massively parallel MD program. This leads to reliable statistics and makes nucleation rates down to the order of 10^30/(m^3 s) accessible to the direct simulation approach. Simulation…
▽ More
Molecular dynamics (MD) simulation is applied to the condensation process of supersaturated vapors of methane, ethane, and carbon dioxide. Simulations of systems with up to a million particles were conducted with a massively parallel MD program. This leads to reliable statistics and makes nucleation rates down to the order of 10^30/(m^3 s) accessible to the direct simulation approach. Simulation results are compared to the classical nucleation theory (CNT) as well as the theory of Laaksonen, Ford, and Kulmala (LFK) which introduces a size dependence of the specific surface energy. CNT describes the nucleation of ethane and carbon dioxide excellently over the entire studied temperature range, whereas LFK provides a better approach to methane at low temperatures.
△ Less
Submitted 6 April, 2009;
originally announced April 2009.