## Abstract

Coupled multidisciplinary systems involve different disciplines/subsystems with feedback-coupled interactions, illustrating the complex interdependencies inherent in real-world engineering systems. Effective monitoring of a coupled multidisciplinary system is crucial for real-time assessment of the interactions between various disciplines within the system. This monitoring provides the data necessary for detecting and addressing issues in a timely manner and facilitates adaptive decision-making for taking reliable design or control actions. However, processing and analyzing data in real time is computationally intensive, and limited resources, such as computational power, sensor capabilities, and budget, may constrain the extent to which a system can be monitored comprehensively. To address this, this article develops a particle-based approach that dynamically selects a subset of sensors that provides the highest information about the state of the system in real time. The proposed approach first predicts the amount of uncertainty in the estimation of the state of the system given noisy measurements from different subsets of available sensors. Then, it selects the sensors that reduce this uncertainty the most, enhancing the precision and efficiency of the monitoring process. The efficacy of the proposed framework is demonstrated via two coupled multidisciplinary systems in the numerical experiments.

## 1 Introduction

Many real-world systems consist of several disciplines interacting with each other dynamically and typically with a great deal of uncertainty associated with each discipline. The interactions between disciplines can be feedforward coupled (unidirectional), in which the output of the upstream discipline becomes the input to the downstream discipline, or feedback coupled (bidirectional), in which the output from one discipline is the input to the other discipline and vice versa. The integration of disciplines and the interactions between them play a key role in the overall behavior of coupled multidisciplinary systems. These systems are common in many applications such as in aerospace and cyber-physical systems [1–9], fluid-thermal–structural modeling and analysis of hypersonic structures [10,11], turbine engine cycle analysis [12–15], satellite performance analysis [16,17], topology optimization [18–20], and more.

The complex nature of coupled multidisciplinary systems, characterized by feedback-coupled interactions, uncertainties, and dynamic behaviors, necessitates monitoring for proactive assessment, optimization, and the continuous improvement of system performance and reliability. The need for effective monitoring of these systems is highlighted by its pivotal role in providing valuable insights into the dynamic interplay between disciplines, offering a comprehensive understanding of the system’s behavior. The information obtained from such monitoring facilitates a deeper comprehension of the complex relationships within the system, which is critical for the early detection of anomalies, deviations, and faults, thereby enabling preemptive intervention. This insight plays a crucial role in conducting reliable control strategies and informing the design of these complex, interconnected engineering systems, ensuring their adaptability, resilience, and optimal performance in dynamic and uncertain operational environments. However, challenges arise in balancing the trade-off between comprehensive monitoring and resource constraints, emphasizing the need for innovative approaches and technologies to overcome these complexities and enhance the effectiveness of monitoring in coupled multidisciplinary systems.

Coupled multidisciplinary systems are often partially observable through sensors that provide noisy measurements from different parts of the system with different levels of uncertainty. While all these sensors provide information about the system and collecting all these data could highly benefit the monitoring of coupled multidisciplinary systems, the limits on computational and/or economic resources prevent the analysis of a large amount of data collected from a large number of sensors for real-time monitoring of the system. Processing and analyzing this diverse and often high-dimensional data in real time is computationally intensive, and delays in data acquisition or processing can impact the system’s responsiveness. Therefore, there is a need to select the optimal subset of sensors in real time that provides the highest information about the system by taking all sources of uncertainty into account. This requires dynamic sensor selection, where the goal is to maximize the accuracy and efficiency of the system monitoring while minimizing the cost and power consumption.

The sensor subset selection problem can be broadly categorized into two main categories developed for static and dynamic systems. In static systems, where the underlying conditions are time independent, the focus is on optimizing the placement and configuration of sensors to capture essential information about the system state. Classical methods for sensor subset selection in static systems include greedy algorithms which iteratively select sensors to maximize information gain or reduce uncertainty [21–24], information theory-based methods utilizing metrics like mutual information and entropy [25–27], optimization techniques aiming to maximize system performance subject to constraints [28–31], heuristic approaches employing intuitive strategies such as thresholding or random selection [32,33], and statistical methods modeling sensor–system relationships through regression analysis or principal component analysis [34,35].

Sensor selection in dynamic systems presents additional challenges due to their time-varying nature, requiring adaptation of selection to changing conditions to capture systems’ transient behavior. Several methods have been developed for sensor scheduling in dynamical systems [36–39]. Classical approaches such as those based on Kalman filtering aim to minimize error covariance and optimize sensor placement for linear dynamic systems [40,41]. Extensions of these techniques to nonlinear systems include methods like extended and unscented Kalman filtering, which provide robust estimation in the presence of nonlinearities [42–44]. These methodologies are primarily applicable to systems that can be linearized, and the stochasticity in state and measurement can be represented by Gaussian noise.

The complex and nonlinear dynamics of coupled multidisciplinary systems with nondifferentiable state and measurement processes and non-Gaussian noises hinder the applicability of existing sensor selection methods. More specifically, the uncertainty inherent in both system states and sensor measurements introduces complexity into the selection process. Also, the need for real-time operation in dynamical multidisciplinary systems necessitates selection methods that are computationally efficient and capable of operating in real time. Addressing these limitations requires the development of advanced sensor selection techniques that can adapt to the complexities of dynamical multidisciplinary systems while accounting for uncertainty, resource constraints, and real-time operation.

In this article, we consider the problem of dynamically selecting an informative subset of sensors for collecting measurement data in coupled multidisciplinary systems to estimate the state of the system in real time. To achieve this, we present a general nonlinear hidden Markov model (HMM) representation of coupled multidisciplinary systems, accounting for different sources of uncertainty in the system’s state transition and sensor measurement noises. Given this representation, we develop a particle-based sensor selection approach that minimizes the predicted expected uncertainty of state estimation in a Bayesian filtering framework. The proposed framework quantifies the expected posterior uncertainty in a particle-based approach, measured by conditional entropy, and chooses the next subset of sensors in a sequential fashion. The key contributions of the proposed framework, differentiating it from existing methods, are its applicability to a general class of hidden Markov models with nonlinearity of processes and non-Gaussianity of noise characteristics, and offering generalizable monitoring solutions for a wide range of arbitrary objectives that can be expressed through the posterior distribution of states. These features make the proposed method suitable for adaptable monitoring in complex multidisciplinary systems.

This article is organized as follows. In Sec. 2, a detailed description of the proposed framework, including a description of coupled multidisciplinary systems, HMM representation of coupled multidisciplinary systems, dynamic sensor selection, and the proposed particle-based sensor selection approach, is provided. This is followed by Sec. 3, which presents the numerical experiments demonstrating the performance of the proposed framework applied to a coupled aerodynamics–structures system and a coupled heart dipole problem. Finally, Sec. 4 contains the concluding remarks.

## 2 Proposed Framework

### 2.1 Coupled Multidisciplinary Systems.

### 2.2 Hidden Markov Model Representation of Coupled Multidisciplinary Systems.

In multidisciplinary systems, the coupling variables are the variables that are shared among multiple disciplines, indicating the disciplines’ internal behavior, affecting other disciplines’ behavior as they get fed into other disciplines, and overall affecting the behavior of the whole system over time as a result of all couplings. Therefore, coupling variables are parts of the system that determine the interactions between disciplines and provide insights into the system-level behavior. In coupled multidisciplinary systems, precisely estimating the values of coupling variables at each time is crucial in order to achieve effective control, design, and analysis of the system. In fact, estimating the correct state of the coupling variables allows one to have accurate knowledge about the current state of the multidisciplinary system as a whole and, consequently, to make informed decisions when choosing control or design inputs to the system or disciplines. Having the functional model of coupled multidisciplinary systems in Eq. (1), we present a general nonlinear/non-Gaussian HMM representation of coupled multidisciplinary systems, consisting of the state and measurement models described below.

*State and Measurement Models:*Let $xk\u2208X$ represent the state vector of the system at time-step $k$ ($k=1,2,\cdots $), containing all coupling variables of the multidisciplinary system, defined as follows:

It should be noted that variations in the type and the number of sensors employed to monitor the system lead to differences in the size of the measurement vector, which may be smaller or larger than the size of the state vector. In real-world large multidisciplinary systems, the limitation of resources, including computational power, energy, and budget, imposes restrictions on the feasibility of comprehensive monitoring. Thus, the strategic identification of a subset of sensors becomes not only a necessity dictated by resource constraints but also a means to enhance the efficiency and effectiveness of the monitoring process. In the following subsection, we describe our proposed approach designed to address this imperative by systematically selecting an optimized subset of sensors to monitor a coupled multidisciplinary system.

### 2.3 Dynamic Sensor Selection.

Real-world complex systems need to be monitored properly for timely detection and reaction to possible changes in the system or for prescribing an appropriate intervention for a particular control or design goal. However, the monitoring process is often resource intensive, time consuming, and costly, as it requires a large number of sensors and a large amount of computing power for processing sensor measurements in real time. Hence, to alleviate this computational burden, we propose a framework to select a small subset of sensors to monitor a coupled multidisciplinary system at any given time.

### 2.4 Proposed Particle-Based Sensor Selection Approach.

#### Proposed Particle-Based Sensor Selection Framework

1: $x0i\u223c\pi 0,\omega 0i=1N,$ for $i=1,\u2026,N$

2: **for**$k=1,2,\u2026$**do**

3: $xki=f(xk\u22121i,vki),$ for $i=1,\u2026,N$⊳ $vki\u223cpv$

4: **for**$\u2200s\u2208S$**do**

5: $zki=h(xki,sk\u22121=s,wki)$⊳ $wki\u223cpw$

6: $\omega ~ki=p(zki|xki,sk\u22121=s)\omega k\u22121i,$ for $i=1,\u2026,N$

7: $\omega ki=\omega ~ki/\u2211j=1N\omega ~kj,$ for $i=1,\u2026,N$

8: $H(p(xk|sk\u22121=s,zk,s0:k\u22122,z1:k\u22121))$

$=\u2212\u2211i=1Nlog(p(zki|xki,sk\u22121=s)\omega k\u22121i)\omega ki$

9: **end for**

10: $sk\u22121*=argmins\u2208SH(p(xk|sk\u22121=s,zk,s0:k\u22122,z1:k\u22121))$

11: Collect observation $zk$ associated with the selected sensors $sk\u22121*$

12: $\omega ~ki=p(zk|xki,sk\u22121=sk\u22121*)\omega k\u22121i,$ for $i=1,\u2026,N$.

13: $\omega ki=\omega ~ki/\u2211j=1N\omega ~kj,$ for $i=1,\u2026,N$.

14: **end for**

## 3 Numerical Experiments

In this section, we demonstrate the performance of the proposed particle-based sensor selection framework for dynamically selecting the most informative sensors to monitor coupled multidisciplinary systems. The proposed framework is implemented on a coupled aerodynamics–structures system and a heart dipole problem.

### 3.1 Coupled Aerodynamics–Structures System.

In the first part of the numerical experiments, we investigate the effectiveness of the proposed framework in the estimation of the state distribution after monitoring the selected optimal subsets with “one” and “two” sensors. The measurement noise variances are considered to be $\sigma wL2=1,\sigma w\varphi 2=2\xd710\u22127$, and $\sigma wL+\varphi 2=1$. Figure 3 shows the particles and the $95%$ confidence intervals associated with the true state distributions and the estimated distributions at time-steps ($a$) $2$, ($b$) $10$, and ($c$) $20$, after monitoring the sensors selected by the proposed framework and “All” and “Random” sensor selection policies. At all time-steps, it can be seen that observing “All” available sensors results in the most accurate estimation of the true state distribution. Nevertheless, this is achieved at the expense of increased processing time and resource utilization due to the computational/economic demands involved in observing all available sensors. On the other hand, the “Random” policy results in the least accurate estimation as it does not take into account the informativeness of the sensor in the selection process and acquiring the observation. However, it can be seen that the proposed framework in both cases of selecting one and two sensors outperforms the random selection strategy and performs close to the case of monitoring all available sensors while offering substantial potential in reducing computational costs.

Figure 4(a) represents $DKL$ as a measure of distance between the true and the estimated state distributions derived from observing different sets of sensors selected by each method over time. The smaller values of $DKL$ indicate the closeness of the true and the estimated state distributions, therefore reflecting the accuracy of the estimation process. As expected, monitoring all available sensors and monitoring one sensor randomly have respectively achieved the smallest and the largest $DKL$ values. It can be seen that the $DKL$ values obtained by the proposed framework in both cases of selecting one and two sensors are considerably small, which demonstrates the effectiveness of the proposed framework in monitoring coupled multidisciplinary systems with lower computational cost.

Figure 4(b) presents the entropy of the state particles at each time-step, averaged over 300 independent runs obtained by the proposed framework after monitoring one and two sensors and the two scenarios of “All” and “Random” sensor selection. The entropy serves as a measure of the amount of uncertainty associated with the estimated state and quantifies the amount of information needed to describe the probability distribution of the state. High entropy indicates greater uncertainty, while low entropy implies more certainty in the estimation process. It can be seen that observing all three sensors has achieved the minimum entropy values at all time-steps, and randomly selecting one sensor has resulted in the maximum entropy value. However, the proposed framework, by selecting only one informative sensor among all available sensors, achieves a lower entropy in state estimation compared to the random scenario. Consistent with expectations, observing two sensors by the proposed framework yields smaller entropy values and more desirable results, lying between observing “All” and one sensor by our method, with a closer trajectory to the “All” sensor selection scenario. This demonstrates the superiority of the proposed framework compared to other scenarios, in which the entropy of the state particles is reduced by observing only the most informative sensor. As explained earlier, higher entropy represents a higher uncertainty in the state estimation. It is shown that the proposed framework achieves a remarkable level of accuracy with high robustness and low computational cost. It should be noted that the proposed framework can be generalized to the selection of any number of sensors among the available ones.

Figures 4(c) and 4(d) illustrate the average absolute errors between the true and the estimated state of the system over 300 independent runs. In accordance with the previous results, observing “All” sensors yields the smallest error values in the estimation of both $L$ and $\varphi $, while “Random” observation has the highest error values among all methods. The estimation errors after observing the sensors selected by the proposed method are within these two extremes, showing a satisfactory level of accuracy and tending to the lower error boundaries in estimating the state of the system.

Figure 5 presents the results obtained by the proposed framework selecting one sensor and the “All” and “Random” sensor selection strategies for different values of measurement noises (noise 1–9) over $100$ independent runs of the coupled aerodynamics-structures system in $20$ time-steps. Figures 5(a), 5(d), and 5(g) show the ratio of each sensor with different measurement noise selected by the proposed framework in $20$ time-steps, averaged over $100$ independent runs. As noted earlier, in rows 1, 2, and 3, the noise associated with sensors 1, 2, and 3 increases, respectively, while others remain constant. As can be seen, the increase in the noise variance of each sensor results in a decrease in the selection rate associated with that sensor. For instance, in Fig. 5(a), as sensor 1 has the smallest variance ($\sigma wL2$) in Noise 1, it has the highest selection ratio compared to noise 2 and 3 as it provides more informative observation, and as its uncertainty in the corresponding measurement increases, the selection ratio for this sensor decreases and for other sensors increases.

Furthermore, as a sensor is more frequently chosen, there is a corresponding decrease in the error of estimating the state associated with the sensor. This is reflected in the second and third columns of Fig. 5, where the mean average errors (MAEs) between the true and the estimated state values of $L$ and $\varphi $ obtained by the proposed framework and two scenarios of “All” and “Random” are compared for different measurement noises. It can be seen that the proposed framework achieves reasonably high accuracy, i.e., low MAE, in estimating the system states by selecting only one, that is the most informative, sensor with estimation errors that are comparable to the “All” scenario. It should be emphasized that in the “All” scenario, all three sensors are observed for the state estimation, while in the proposed framework, only one sensor is selected for acquiring a measurement. Consistent with the previous results and as expected, the “Random” scenario has the highest estimation error unless a sensor is not frequently selected in the proposed framework because of the corresponding high measurement noise. This can be seen, for instance, in noise 6 in Fig. 5($f$), where the MAE of $\varphi $ obtained by the proposed framework is higher than the random selection scenario as shown in noise 6 in Fig. 5($d$), the selection ratio of sensor 2 is noticeably small.

It is important to emphasize that the proposed framework prioritizes minimizing system-level uncertainty, quantified in terms of entropy, through the intelligent selection of sensors. The observed higher mean average error in some cases is attributed to the inherent characteristics of the sensors, particularly the noise levels of sensors and their informativeness. In scenarios where certain sensors are not selected for measurement due to their low informativeness or high noise levels, the corresponding state variables exhibit higher errors compared to those where sensors have been selected. However, it is crucial to emphasize that the proposed method aims to minimize the overall entropy of the posterior distribution of the state variables, thereby enhancing the overall estimation performance of the system. It is also worth noting that the proposed method can be adapted to minimize the entropy of the marginalized state distributions, accounting for the importance of states to be estimated accurately instead of the entropy of the joint distribution. By incorporating this refinement, the framework can prioritize reducing uncertainty in critical state variables while ensuring that less critical variables are adequately estimated, allowing for a more targeted sensor selection strategy tailored to the specific requirements and objectives of the system under consideration.

### 3.2 Heart Dipole Problem.

The heart dipole problem represents a well-established test case in the field of multidisciplinary design optimization, originally introduced in the literature by NASA [47]. This problem involves the experiment of measuring a synthetic dipole moment generated by a human heart through the use of two artificial dipoles contained within an electrolyte disk. The human heart dipole has been converted into a test problem for multidisciplinary system analysis methods by defining two disciplines, shown in Fig. 6. In this problem, the two disciplines are structured to address different aspects of the measurement process. The first discipline focuses on optimizing the placement of artificial electrodes within an electrolyte disk to simulate the measurement of the dipole moment generated by a human heart. This discipline, known as electrode placement, focuses on determining the optimal locations and orientations of the electrodes to maximize the accuracy and reliability of the measured signal while minimizing interference and noise. The second discipline involves signal processing techniques applied to the recorded data from the electrodes. This discipline, called signal processing, encompasses various methods, such as filtering, feature extraction, and signal reconstruction, aimed at extracting meaningful information about the dipole moment parameters from the recorded signals. Optimizing electrode placement and signal processing disciplines enables a comprehensive analysis and enhancement of the measurement system’s performance in capturing the underlying physiological activity of the heart. This problem can be stated as a set of four nonlinear equations and eight variables ($x1,\u2026,x8$ as coordinates and magnitudes of each dipole on $x$ and $y$ axis), where four variables ($x1,x4,x6andx7$) and four pieces of measured data (parameters $dmx,dmy,dAanddB$) are the inputs to the system, and four variables ($x2,x3,x5,$ and $x8$) are the coupling variables [48], which define the system state, represented as follows:

We implement the proposed sensor selection approach on the heart dipole problem to demonstrate its performance in the state estimation of this coupled multidisciplinary system given noisy measurement from a single selected sensor at each time. Similar to previous experiments, the proposed framework is compared to two scenarios of “All” and “Random” sensor selection. In the “All” scenario, all four sensors are observed simultaneously at each time-step, and their measurements are used for state estimation. In the “Random” scenario, one sensor is selected randomly to acquire an observation, regardless of the amount of information the sensor provides.

Figure 7 presents the results obtained over $500$ independent runs of the heart dipole problem in $20$ time-steps. Figure 7(a) shows the ratio of each sensor selected by the proposed framework in $20$ time-steps, averaged over $500$ independent runs. It can be seen that sensors 1 and 2 have been chosen more frequently, indicating that they possess the highest information compared to the other two sensors. This serves as the advantage of the proposed approach, as it eliminates the necessity of allocating computational resources toward noninformative sensors. It should be noted that the computational cost associated with acquiring measurements from all sensors is four times of the proposed framework and the random selection strategy, as in the latter cases, only one measurement from a single sensor is acquired.

Figure 7(b) illustrates the entropy of the state particles at each time-step averaged over $500$ independent runs of sensor selection by the proposed framework as well as the two scenarios of “All” and “Random” sensor selection. It is seen that the entropy values at each time-step obtained by the proposed method fall between the entropy values obtained by the “All” and “Random” sensor selection methods. It is essential to acknowledge that while the proposed method has not achieved the same level of entropy reduction as scenario “All” and it is closer to the “Random” entropy levels, it has achieved a high efficiency in estimating the state distribution by recognizing the most informative sensor at each time. This efficiency is demonstrated in Fig. 7(c) through the Kullback–Leibler divergence ($DKL$) between the true and the estimated state distributions at each time-step over $500$ independent runs. It can be seen that the scenario of observing all available sensors achieves the minimum $DKL$ values, while the “Random” sensor selection yields significantly higher $DKL$ values. Notably, in this result, one can observe that the proposed framework has achieved $DKL$ values close to “All” scenario. This result illustrates that the proposed method has been effective in state estimation with a high accuracy level by observing only one informative sensor instead of all four available sensors.

In this part of numerical experiments, we evaluate the performance of the proposed framework in the presence of different combinations of sensors monitoring the heart dipole problem. The experiment considers various configurations of available sensors to monitor the coupled multidisciplinary system. The proposed framework is implemented to select the most informative sensor among all available ones. Table 1 presents the average entropy of state particles per step and mean absolute error of state per step over $500$ independent runs for various combinations of available sensors in the heart dipole problem. As with previous experiments, the proposed framework is compared to two other scenarios, “All” and “Random” sensor selection. The average entropy values are computed to assess the performance of the proposed framework in minimizing the uncertainty in the state estimation process, and the mean absolute errors quantify the accuracy of the state estimation. The results show that the proposed framework, which selects only one informative sensor, achieves both low uncertainty and high accuracy in state estimation. This performance is comparable to the scenario where all available sensors are selected and surpasses the performance of the random selection method. It can be seen that regardless of the sensor combination or the presence/absence of any specific sensor, the proposed framework consistently identifies the most informative sensor among the available ones. These results demonstrate the adaptability and generalizability of the proposed framework across different monitoring scenarios and highlight its capability to maintain estimation accuracy while optimizing resource efficiency.

Average entropy per step | Mean absolute error of state per step | |||||
---|---|---|---|---|---|---|

Available sensors | Proposed framework | All | Random | Proposed framework | All | Random |

$S1,S2,S3,S4$ | $\u22124.4911$ | $\u22124.8507$ | $\u22124.2226$ | $0.0717$ | $0.0617$ | $0.0925$ |

$S1,S2,S3$ | $\u22124.3976$ | $\u22124.7481$ | $\u22124.2342$ | $0.0751$ | $0.0647$ | $0.0991$ |

$S1,S2,S4$ | $\u22124.4805$ | $\u22124.8482$ | $\u22124.2972$ | $0.0730$ | $0.0617$ | $0.0872$ |

$S1,S3,S4$ | $\u22124.4910$ | $\u22124.5032$ | $\u22124.1861$ | $0.0648$ | $0.0641$ | $0.0903$ |

$S2,S3,S4$ | $\u22124.4910$ | $\u22124.5318$ | $\u22124.1682$ | $0.0822$ | $0.0801$ | $0.1037$ |

$S1,S2$ | $\u22124.3978$ | $\u22124.7454$ | $\u22124.3624$ | $0.0749$ | $0.0649$ | $0.0847$ |

$S1,S3$ | $\u22124.3979$ | $\u22124.4006$ | $\u22124.1854$ | $0.0673$ | $0.0672$ | $0.0919$ |

$S1,S4$ | $\u22124.4704$ | $\u22124.5015$ | $\u22124.2901$ | $0.064$ | $0.0634$ | $0.0752$ |

$S2,S3$ | $\u22124.3184$ | $\u22124.3329$ | $\u22124.1488$ | $0.1042$ | $0.1101$ | $0.1083$ |

$S2,S4$ | $\u22124.4911$ | $\u22124.5271$ | $\u22124.2562$ | $0.1030$ | $0.0810$ | $0.0941$ |

$S3,S4$ | $\u22124.1792$ | $\u22124.1838$ | $\u22124.0858$ | $0.0833$ | $0.0833$ | $0.0970$ |

Average entropy per step | Mean absolute error of state per step | |||||
---|---|---|---|---|---|---|

Available sensors | Proposed framework | All | Random | Proposed framework | All | Random |

$S1,S2,S3,S4$ | $\u22124.4911$ | $\u22124.8507$ | $\u22124.2226$ | $0.0717$ | $0.0617$ | $0.0925$ |

$S1,S2,S3$ | $\u22124.3976$ | $\u22124.7481$ | $\u22124.2342$ | $0.0751$ | $0.0647$ | $0.0991$ |

$S1,S2,S4$ | $\u22124.4805$ | $\u22124.8482$ | $\u22124.2972$ | $0.0730$ | $0.0617$ | $0.0872$ |

$S1,S3,S4$ | $\u22124.4910$ | $\u22124.5032$ | $\u22124.1861$ | $0.0648$ | $0.0641$ | $0.0903$ |

$S2,S3,S4$ | $\u22124.4910$ | $\u22124.5318$ | $\u22124.1682$ | $0.0822$ | $0.0801$ | $0.1037$ |

$S1,S2$ | $\u22124.3978$ | $\u22124.7454$ | $\u22124.3624$ | $0.0749$ | $0.0649$ | $0.0847$ |

$S1,S3$ | $\u22124.3979$ | $\u22124.4006$ | $\u22124.1854$ | $0.0673$ | $0.0672$ | $0.0919$ |

$S1,S4$ | $\u22124.4704$ | $\u22124.5015$ | $\u22124.2901$ | $0.064$ | $0.0634$ | $0.0752$ |

$S2,S3$ | $\u22124.3184$ | $\u22124.3329$ | $\u22124.1488$ | $0.1042$ | $0.1101$ | $0.1083$ |

$S2,S4$ | $\u22124.4911$ | $\u22124.5271$ | $\u22124.2562$ | $0.1030$ | $0.0810$ | $0.0941$ |

$S3,S4$ | $\u22124.1792$ | $\u22124.1838$ | $\u22124.0858$ | $0.0833$ | $0.0833$ | $0.0970$ |

## 4 Conclusion

This article presented a dynamic sensor selection approach for efficient monitoring of partially observable coupled multidisciplinary systems. In many multidisciplinary systems, a vast array of sensors is employed to capture various aspects of the system’s behavior. However, analyzing data from all sensors can be computationally intensive and time consuming. To address this challenge, we proposed a sensor selection framework that identifies the most informative sensors for state estimation in real time, making a balance between estimation accuracy and resource efficiency. This is achieved through a particle-based Bayesian filtering technique that probabilistically predicts the uncertainty in the state estimation after monitoring any set of sensors. The efficiency, accuracy, and robustness of state estimation by the proposed framework are demonstrated in the numerical experiments. Future work will include enhancing the scalability of the proposed framework to accommodate larger-scale coupled multidisciplinary systems by integrating adaptive sampling strategies alongside approximate inference methods and optimized particle resampling techniques.

## Funding Data

• This work was supported by the 10.13039/100006754 ARMY Research Laboratory (Award numbers W911NF2410098 and W911NF2320179).

## Conflict of Interest

There are no conflicts of interest.

## Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.