Sponsor: U.S. National Science Foundation
Project No.: CMMI-2317172
Duration: September 2023 – August 2026
Principal Investigator: Professor Sharif Rahman
Graduate Students: Md. Rashel Talukdar
SUMMARY
The principal goal of this research is to conduct foundational research on inventing efficient computational algorithms and creating practical computational tools for robust and reliability-based design optimization (RDO and RBDO) of high-dimensional complex systems subject to random input following an arbitrary dependent probability distribution.
The technical activities to meet the project goal, as outlined in the original proposal, comprise the following three tasks:
- Task 1: Novel mathematical developments of a generalized analysis-of-variance expansion, leading to a generalized spline dimensional decomposition (GSDD) for tackling dependent random variables head on;
- Task 2: New scalable algorithms of the GSDD method for calculating relevant probabilistic response characteristics and design sensitivities of a high-dimensional, complex mechanical system; and
- Task 3: Innovative GSDD-driven optimization algorithms for efficiently solving high-dimensional RDO and RBDO problems, including stochastic shape and topology design.
Research
Year 1:
Higher-Order Moment Analysis of Spline Chaos Expansion:
Year 2:
A Generalized Polynomial Chaos Expansion (GPCE) for Global Sensitivity Analysis with Dependent Random Variables
Time-dependent Uncertainty Quantification Analysis of Complex Dynamical Systems:
Year 3:
Year 1:
Higher-Order Moment Analysis of Spline Chaos Expansion: A UQ analysis entailing high-order moments calculated from SCE approximations of a real-valued, general output function of input random variable was conducted. The approximation quality of SCE was assessed in terms of the modulus of smoothness of the function. When the largest element of the mesh from SCE approaches zero, the modulus of smoothness vanishes, resulting in the higher-order moment convergence of SCE to the correct limit. Therefore, the moment of SCE of an arbitrary order converges to the exact moment for any degree of splines as the element size decreases. Moreover, the weaker modes of convergence, such as those in probability and in distribution, transpire naturally.
Numerical computations of moments by SCE and polynomial chaos expansion (PCE), conducted for a collection of simple yet relevant examples, indicate the following:
- When the output function is smooth or nonsmooth and the random domain is bounded or unbounded, SCE and PCE both provide convergent estimates of the variance. However, their relative convergent rates may differ, depending on the smoothness of the function.
- When the function is globally smooth, PCE is likely to provide more efficient estimates of variance than SCE for the same computational effort. However, for a nonsmooth function, the trend reverses, and the convergence properties of PCE in estimating variance may degrade appreciably.
- Higher-order moments, such as skewness and kurtosis, calculated using SCE converge for all examples considered in this study. In contrast, moments of PCE of orders larger than two may or may not converge, depending on the regularity of the output function or the probability measure of input random variables.
Year 2:
A Generalized Polynomial Chaos Expansion (GPCE) for Global Sensitivity Analysis with Dependent Random Variables: A GPCE method was developed for variance-based global sensitivity analysis of uncertain systems with statistically dependent input random variables following a general, non-product-type probability distribution. The method is built upon three key components:(1) a general three-step process for constructing multivariate orthonormal polynomials consistent with the probability measure of input variables; (2) a Fourier series expansion of a square-integrable output function with respect to these multivariate polynomials; and (3) novel analytical formulae for global sensitivity indices in terms of the expansion coefficients and expectations of the products between relevant orthonormal polynomials. For statistically independent input variables, the proposed GPCE simplifies to the classical PCE, recovering established expressions of sensitivity analysis. By tackling dependent random variables head-on, GPCE circumvents the need for a measure transformation between dependent and independent variables, thereby retaining the original output function and hence calculating the sensitivity indices accurately and efficiently. Once the GPCE approximation is established — typically as part of uncertainty quantification analysis — the sensitivity indices can be generated with minimal extra effort. Therefore, both uncertainty quantification and sensitivity analyses can be performed simultaneously.
Results from four illustrative examples involving elementary mathematical functions show that GPCE yields exact or rapidly converging sensitivity indices for both Gaussian and non-Gaussian input distributions. Furthermore, application to an industrial-scale problem demonstrates GPCE’s practical utility, integrating seamlessly with legacy simulation code to identify influential and non-influential variables during random eigenvalue analysis of a complex mechanical system.
Time-dependent Uncertainty Quantification Analysis of Complex Dynamical Systems: A novel computational methodology, supported by robust numerical algorithms, was developed for time-dependent uncertainty quantification (UQ) analysis on complex dynamical systems. The proposed approach is based on (1) a new stochastic adaptation of the nonlinear autoregressive with exogenous input (NARX) model, utilizing dimension-wise tensor product expansion to effectively capture the behavior of dynamical systems, (2) a polynomial dimensional decomposition (PDD) technique to propagate uncertainty in input random variables to the NARX coefficients, and (3) a unique integration between NARX and PDD, resulting in the PDD-NARX approximation for time-dependent UQ analysis. The PDD-NARX method distinguishes itself from conventional deterministic system identification tools by considering uncertainties originating from both the system’s dynamic properties (e.g., mass, stiffness, and damping) and external forces (e.g., amplitude and frequency content of excitation time series). Unlike traditional methods, which rely on an intuitive selection of NARX basis functions, this approach employs dimensional decomposition and importance factors to systematically construct the NARX model function. Furthermore, PDD, due to its hierarchical, dimension-wise expansion, is better equipped to manage high-dimensional UQ problems than many existing methods, including the widely recognized polynomial chaos expansion.
Numerical results demonstrate that low-order PDD-NARX approximations provide accurate and computationally efficient estimates of the probabilistic characteristics of simple dynamical systems. Moreover, the probabilistic vehicle dynamic analysis of a pick-up truck traversing road bumps underscores the effectiveness of the PDD-NARX method in addressing industrial-scale complex problems.
Year 3:
Education
Supervision of Graduate Students: A new PhD student, Mr. M. R. Talukdar, joined the PI’s research group on August 15, 2024. The student, appointed as a half-time graduate research assistant, is currently supported by this NSF-funded project. The PI serves as the student’s academic advisor.
Insertion of Research Results into Undergraduate and Graduate Courses: In 2022, the PI developed a new intermediate-level course entitled Uncertainty Quantification and Design Optimization, designed to introduce undergraduate and graduate students to contemporary methods of probabilistic analysis and design optimization under uncertainty. During its second offering in Fall 2024 (August–December), the course was enhanced with novel uncertainty quantification and robust/reliability-based design optimization methods derived from the current and previous NSF-funded projects.
Invited Seminars and Lectures: The PI’s research group received invitations to present five lectures on uncertainty quantification and design optimization at the following prominent conferences:
14th International Conference on Structural Safety and Reliability (ICOSSAR25), Los Angeles, CA, June 1-6, 2025.
Engineering Mechanics Institute Conference and Probabilistic Mechanics & Reliability Conference (EMI/PMC 2024), Chicago, IL, May 28-31, 2024.
2024 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference (IDETC-CIE2024), Washington, DC, August 25-28, 2024.
Publications
Journal Articles
- Rahman, S., “A Generalized Polynomial Chaos Expansion for Global Sensitivity Analysis with Dependent Random Variables,” submitted to SIAM Journal on Uncertainty Quantification, 2025.
- Ebadollahi, M. and Rahman, S., “Time-dependent Uncertainty Quantification Analysis of Complex Dynamical Systems,” Probabilistic Engineering Mechanics, Vol. 80, Article 103776, pp. 1-19, 2025.
- Rahman, S., “Higher-Order Moments of Spline Chaos Expansion,” Probabilistic Engineering Mechanics, Vol. 77, Article No. 103666, 2024.
- Lee, D. and Rahman, S., “High-Dimensional Stochastic Design Optimization under Dependent Random Variables by a Dimensionally Decomposed Generalized Polynomial Chaos Expansion,” International Journal for Uncertainty Quantification, Volume 13, No. 4, pp. 23-59, 2023.
Papers in Conferences, Symposiums, and Congresses
- Talukdar, M. R. and Rahman, S., “Practical Uncertainty Quantification Analysis in the Presence of Statistically Dependent Random Variables,” Proceedings of the 14th International Conference on Structural Safety and Reliability (ICOSSAR25), Los Angeles, CA, June 1-6, 2025.
- Ebadollahi, M. and Rahman, S., “Time-dependent Uncertainty Quantification Analysis of Complex Dynamical Systems,” Proceedings of the 14th International Conference on Structural Safety and Reliability, Los Angeles, CA, June 1-6, 2025.
- Ebadollahi, M. and Rahman, S., “A PDD-NARX Approximation for Uncertainty Quantification in Complex Dynamical Systems,” Proceedings of the 2024 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference, Washington, DC, August 25-28, 2024.
- Rahman, S., “A Spline Chaos Expansion for Uncertainty Quantification in Linear Dynamical Systems,” Proceedings of 2024 Engineering Mechanics Institute Conference and Probabilistic Mechanics & Reliability Conference (EMI/PMC 2024), Chicago, IL, May 28-31, 2024.
- Lee, D. and Rahman, S., “A Generalized Polynomial Chaos Expansion for High-Dimensional Design Optimization under Dependent Random Variables,” Proceedings of 2024 Engineering Mechanics Institute Conference and Probabilistic Mechanics & Reliability Conference (EMI/PMC 2024), Chicago, IL, May 28-31, 2024.
Others
Software
Several computer codes have been developed to perform the various numerical simulations required for this project. These codes were created exclusively for research purposes. For further information, please contact the PI (S. Rahman) at sharif-rahman@uiowa.edu.