1 Introduction

Multiobjective optimization (MO) deals with optimization problems that have at least two (conflicting) objective functions to be optimized simultaneously. The conflict between the objectives means that their optimal per objective values cannot be obtained by a single solution and, thus, some balancing (trading off) between the objectives is needed. Mathematically, an MO problem can be formulated as:

$$\begin{aligned} \begin{aligned} \min \;&{\mathbf {f}}({\mathbf {x}}) = \left\{ f_1({\mathbf {x}}),\dots ,f_k({\mathbf {x}})\right\} \\ \text {s.t.}\;&{\mathbf {x}}\in S, \end{aligned} \end{aligned}$$
(1)

where the objective vector \({\mathbf {f}}\) consists of k \((\ge 2)\) objective functions \(f_i:{\mathbb {R}}^n\rightarrow {\mathbb {R}}, i=1,\dots ,k\). The vector \({\mathbf {x}}=(x_1,\dots ,x_n)^T\), consisting of n decision variables, is constrained by the feasible region \(S \subset {\mathbb {R}}^n\).

The optimal solutions for problem (1) are often defined by using Pareto optimality. A feasible solution \({\mathbf {x}}^*\in S\) is called Pareto optimal if there does not exist another feasible solution \({\mathbf {x}}\in S\) \(({\mathbf {x}}\ne {\mathbf {x}}^*)\) such that \(f_i({\mathbf {x}})\le f_i({\mathbf {x}}^*)\) for all \(i=1,\dots ,k\) and \(f_j({\mathbf {x}})< f_j({\mathbf {x}}^*)\) for at least one j. The set of all Pareto optimal solutions is said to form the Pareto set, whose image in the objective space is known as the Pareto front. All Pareto optimal solutions are mathematically incomparable without any additional information.

When solving real-world MO problems, the idea is often to find a Pareto optimal solution or a small number of Pareto optimal solutions that are most preferred in the context of the problem considered. Then, a decision maker (DM) is a person who has the knowledge and understanding of the problem domain and can provide preference information regarding the conflicting objectives. There exist many methods for solving MO problems and they can be divided into categories based on the role of the DM in helping the methods to find the most preferred solution (Miettinen 1999): (1) In a priori methods, the DM is first asked to provide preferences that are then taken account in optimization to find a Pareto optimal solution that is closest to the preferences. (2) In a posteriori methods, an approximation of the whole Pareto front is first obtained and shown to the DM, who then selects the most preferred solution amongst them. (3) In interactive methods, the DM provides preference information iteratively and thus provides feedback to the method during the optimization (Miettinen et al. 2016).

Both mathematical and evolutionary optimization approaches have been used to solve MO problems. Mathematical optimization methods (Miettinen 1999) are efficient when functions are differentiable and gradient information is available or the structure of the problem can be exploited, while evolutionary optimization methods (Deb 2005; Emmerich and Deutz 2018) perform well on problems where not much or no information of the problem structure is available (e.g. black-box problems) or the problem features non-standard properties such as noisy/uncertain functions and multiple objectives (as they allow for the discovery of multiple Pareto optimal solutions in a single algorithmic run). Further, mathematical optimization methods offer guarantees of converge, especially for convex problems, while evolutionary optimization methods do not guarantee convergence but are well suited for non-convex problems and can provide a set of solutions for a single optimization run. Good results have been obtained by hybridizing these two approaches (Branke et al. 2008).

When solving real-world MO problems, a common challenge is the non-availability of closed-form expressions of the objective and/or constraint functions, also known as black-box functions. Consequently, the evaluation of these functions may involve the execution of computer simulations and/or real (physical, biological, chemical) experiments. Such evaluations are typically time-consuming, costly and/or resource-intense, forcing the optimizer to arrive at a good solution within a limited number of evaluations. Hence, a lot of thought goes into deciding which solutions to evaluate at each iteration of the optimization. A common approach to cope with expensive objective/constraint function(s) is to replace them with a cheap to evaluate surrogate model (Allmendinger et al. 2017; Chugh et al. 2019). By optimizing over the surrogate model using an acquisition function, such as expected improvement, one can identify one or multiple solutions that are likely to lead to an improvement of the current best solution; these promising solutions can then be evaluated on the original expensive objective function.

In the presence of multiple objective functions, there is no reason to believe that the objectives share the same properties. In fact, the objectives may be heterogeneous varying in aspects such as evaluation time (slow vs fast to evaluate objectives) (Allmendinger et al. 2015), the induced landscape (rugged vs smooth), safety risks (causing serious consequences when evaluated below some safety threshold vs no consequences) (Kim et al. 2021), level of noise (deterministic vs noisy) and there may also be varying levels of correlations between the objectives throughout the search space; for a review of different types of heterogeneity please refer to Allmendinger and Knowles (2021). The use of surrogate models serves as a good starting point to cope with heterogeneous objectives. Surrogate models can also be used when the problem formulation has to be based only on some existing data (e.g. measurement data from experiments or sensors, business data, health data), and obtaining new data is either impossible or out of the hands of the optimizer; this scenario is also referred to as offline data-driven optimization as opposed to online data-driven optimization where the optimizer can decide when to evaluate which solution (Jin et al. 2019).

2 Overview of the special issue articles

The presence of multiple objectives in a problem increases the complexity of the decision making process as it requires both the optimizer and DM to trade-off multiple conflicting factors. An increase in the number of objectives poses also an additional layer of complexity when it comes to theoretical research in this domain. The eight papers in this special issue address both experimental and theoretical challenges arising in engineering sciences.

In Korondi et al. (2021) propose a novel optimization method for MO (and single-objective) problems featuring both expensive and uncertain evaluations of candidate solutions. The problems under investigation have highly limited computational budgets. A multi-fidelity surrogate-assisted approach is proposed, and has shown, on various test problems, to outperform standard surrogate techniques if cheap low-fidelity evaluations are available. The approach can be applied to problems where new data can be obtained by evaluating solutions, or where new data cannot be obtained but there is a surrogate model available of the problem under investigation. These are common setups in engineering applications.

An expensive MO problem was also considered by Prinz et al. (2021). More precisely, the authors consider the optimization of a mixing process in electrically conducting fluids. To tackle this problem, a trust-region method was modified to tackle a problem with multiple expensive objective functions. Moreover, the proposed method guarantees an improvement of the objective function values from a starting guess. The study allowed the authors not only to optimize the mixing process but also gain insight into the relationship between the decision variables and the objectives.

A new approach to solve multi-attribute decision making problems with subjective preferences and non-preferences of the DM is proposed by Mehlawat et al. (2021). Triangular fuzzy preference relations and triangular fuzzy non-preference relations are used to represent the subjective preferences and non-preferences, respectively. The proposed approach is validated using a multi-attribute decision making problem from the literature, and compared against similar existing approaches with good results being reported.

In Sekizaki et al. (2021) propose a novel methodology to tackle a many-objective engineering problem in the area of network reconfiguration of an electric distribution system. The challenge in the application considered is to determine how best to downsize equipment (distribution lines) such that costs are reduced while feasibility of the network reconfigurations (due to maintenance and/or replacement of distribution lines) is guaranteed. The proposed methodology utilises a modified version of the NSGA-III algorithm (Deb and Jain 2014) to identify a set of trade-off solutions to be presented to a DM (distribution system operator in this case), who then makes a final decision on which distribution lines to downsize. In collaboration with a power grid company, the authors then validate the methodology on a large-scale distribution system model and identify, within reasonable computational time, distribution lines suitable for downsizing. This very practical approach can assist the distribution system operator in long-term planning.

Yet another real-world application is addressed in the manuscript of Ortiz-Martin et al. (2021). Here the application is concerned with the design of a wastewater treatment process such that two objectives, operational costs and effluent quality, are optimized. The authors formulate this problem as a dynamic MO problem with the aim being to find the time varying control profiles (aeration factor in the last aerated reactor and the internal recycle flow rate) in order to optimize the objectives. An NSGA-II algorithm (Deb et al. 2002) coupled with scenario-based simulations is used to solve the problem. This methodology is of particular use in situations where the energy costs are fluctuating hourly, supplying the decision maker a set of possible wastewater treatment processes to implement depending on the actual economical, technical or legal circumstances.

The design of a fin to optimize its effectiveness and efficiency using MO is considered by Weerasena et al. (2021). The authors use a piecewise linear approximation of the design of a fin. The problem itself is then optimized using \(\epsilon\)-constraint method. It is demonstrated that the approach is able to find trade-off solutions of similar quality regardless of the initial design of the fin and the Biot number (a dimensionless index used in heat transfer calculations). This research supports the design of optimized fins and, as noted by the authors, the discretization method proposed can also find application in other areas of heat transfer and mechanics.

In Yano and Nishizaki (2021) focus on the application of MO concepts to two-level programming problems, also known as Stackelberg games. In these games, a leader (e.g. natural gas shipping company) chooses a strategy to start off, subsequently prompting a follower (e.g. pipeline operating company) to choose a strategy having knowledge of the leader’s decision. The authors augment such problems with the notion of simple resources, which are shortages and excesses arising from the violation of constraints with discrete random variables. The simple resources are accounted for in the objective functions as penalty functions. By combining the concept of Pareto Stackelberg solutions with the Kuhn–Tucker approach, the authors are able to formulate the problem as a mixed-integer programming problem, which they solve using an interactive algorithm. Numerical experiments show that the proposed methodology can react efficiently to simple resource changes.

Finally, in Schultes et al. (2021) tackle a bi-objective shape optimization problem. The goal is to maximize the mechanical stability of a ceramic component under tensile load while minimizing its volume. Hypervolume scalarization is suggested to tackle the (discretized) optimization problem as it does not require an explicit handling of additional constraints and, compared to the weighted sum approach, allows for the discovery of unsupported non-dominated solutions. When augmenting hypervolume maximization on an iterative ascent algorithm, the authors show that for convex problems a non-trivial reference set generating the set of properly non-dominated solutions can be determined without prior knowledge on the Pareto front. The approach is also validated on a real ceramic component design problem allowing the authors to gain valuable insights into the relationship between the component shape and its reliability and cost.

3 Conclusions

The articles in the special issue cover a wide range of topics related to MO, including both mathematical and evolutionary optimization approaches, preference elicitation techniques, approaches for solving expensive optimization problems, and approaches drawing from game theory. A strong engineering aspect is also present since many of the papers are dealing with a concrete real-world case study, which emphasizes a clear need for MO approaches in engineering sciences. It is also great to see that in many of the papers, the decision making aspect was present since the authors were not just happy to only produce Pareto optimal solutions but did also consider the decision making aspect of identifying practically relevant solutions. To summarize, MO has become an essential tool in many fields not least in engineering sciences, and it will continue to provide a platform to rely on in the future.

This special issue was motivated by the MCDM2019 conference held in Istanbul, Turkey, which belongs to the trademark bi-annual conference series of the International Society on Multiple Criteria Decision Making. We highly appreciate the support from the conference organizers. Altogether, 27 manuscripts were submitted for the special issue of which 8 were finally accepted. This gives an acceptance rate of about 30%. We would like to thank all the reviewers that were spending their precious time and provided high-quality reviews helping us to identify the manuscripts that were accepted for this special issue.