共查询到20条相似文献，搜索用时 15 毫秒
1.
M. Arnst R. Ghanem E. Phipps J. Red‐Horse 《International journal for numerical methods in engineering》2014,97(5):352-376
We address the curse of dimensionality in methods for solving stochastic coupled problems with an emphasis on stochastic expansion methods such as those involving polynomial chaos expansions. The proposed method entails a partitioned iterative solution algorithm that relies on a reduced‐dimensional representation of information exchanged between subproblems to allow each subproblem to be solved within its own stochastic dimension while interacting with a reduced projection of the other subproblems. The proposed method extends previous work by the authors by introducing a reduced chaos expansion with random coefficients. The representation of the exchanged information by using this reduced chaos expansion with random coefficients enables an expeditious construction of doubly stochastic polynomial chaos expansions that separate the effect of uncertainty local to a subproblem from the effect of statistically independent uncertainty coming from other subproblems through the coupling. After laying out the theoretical framework, we apply the proposed method to a multiphysics problem from nuclear engineering. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
2.
M. Arnst R. Ghanem E. Phipps J. Red‐Horse 《International journal for numerical methods in engineering》2012,92(12):1044-1080
Coupled problems with various combinations of multiple physics, scales, and domains are found in numerous areas of science and engineering. A key challenge in the formulation and implementation of corresponding coupled numerical models is to facilitate the communication of information across physics, scale, and domain interfaces, as well as between the iterations of solvers used for response computations. In a probabilistic context, any information that is to be communicated between subproblems or iterations should be characterized by an appropriate probabilistic representation. Although the number of sources of uncertainty can be expected to be large in most coupled problems, our contention is that exchanged probabilistic information often resides in a considerably lower‐dimensional space than the sources themselves. In this work, we thus propose to use a dimension reduction technique for obtaining the representation of the exchanged information, and we propose a measure transformation technique that allows subproblem implementations to exploit this dimension reduction to achieve computational gains. The effectiveness of the proposed dimension reduction and measure transformation methodology is demonstrated through a multiphysics problem relevant to nuclear engineering. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
3.
M. Arnst R. Ghanem E. Phipps J. Red‐Horse 《International journal for numerical methods in engineering》2012,92(11):940-968
Coupled problems with various combinations of multiple physics, scales, and domains are found in numerous areas of science and engineering. A key challenge in the formulation and implementation of corresponding coupled numerical models is to facilitate the communication of information across physics, scale, and domain interfaces, as well as between the iterations of solvers used for response computations. In a probabilistic context, any information that is to be communicated between subproblems or iterations should be characterized by an appropriate probabilistic representation. Although the number of sources of uncertainty can be expected to be large in most coupled problems, our contention is that exchanged probabilistic information often resides in a considerably lower dimensional space than the sources themselves. This work thus presents an investigation into the characterization of the exchanged information by a reduced‐dimensional representation and in particular by an adaptation of the Karhunen‐Loève decomposition. The effectiveness of the proposed dimension–reduction methodology is analyzed and demonstrated through a multiphysics problem relevant to nuclear engineering. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
4.
P. Surya Mohan Prasanth B. Nair Andy J. Keane 《International journal for numerical methods in engineering》2011,85(7):874-895
We present stochastic projection schemes for approximating the solution of a class of deterministic linear elliptic partial differential equations defined on random domains. The key idea is to carry out spatial discretization using a combination of finite element methods and stochastic mesh representations. We prove a result to establish the conditions that the input uncertainty model must satisfy to ensure the validity of the stochastic mesh representation and hence the well posedness of the problem. Finite element spatial discretization of the governing equations using a stochastic mesh representation results in a linear random algebraic system of equations in a polynomial chaos basis whose coefficients of expansion can be non‐intrusively computed either at the element or the global level. The resulting randomly parametrized algebraic equations are solved using stochastic projection schemes to approximate the response statistics. The proposed approach is demonstrated for modeling diffusion in a square domain with a rough wall and heat transfer analysis of a three‐dimensional gas turbine blade model with uncertainty in the cooling core geometry. The numerical results are compared against Monte–Carlo simulations, and it is shown that the proposed approach provides high‐quality approximations for the first two statistical moments at modest computational effort. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
5.
M. Raisee D. Kumar C. Lacor 《International journal for numerical methods in engineering》2015,103(4):293-312
In this paper, a non‐intrusive stochastic model reduction scheme is developed for polynomial chaos representation using proper orthogonal decomposition. The main idea is to extract the optimal orthogonal basis via inexpensive calculations on a coarse mesh and then use them for the fine‐scale analysis. To validate the developed reduced‐order model, the method is implemented to: (1) the stochastic steady‐state heat diffusion in a square slab; (2) the incompressible, two‐dimensional laminar boundary‐layer over a flat plate with uncertainties in free‐stream velocity and physical properties; and (3) the highly nonlinear Ackley function with uncertain coefficients. For the heat diffusion problem, the thermal conductivity of the slab is assumed to be a stochastic field with known exponential covariance function and approximated via the Karhunen–Loève expansion. In all three test cases, the input random parameters are assumed to be uniformly distributed, and a polynomial chaos expansion is found using the regression method. The Sobol's quasi‐random sequence is used to generate the sample points. The numerical results of the three test cases show that the non‐intrusive model reduction scheme is able to produce satisfactory results for the statistical quantities of interest. It is found that the developed non‐intrusive model reduction scheme is computationally more efficient than the classical polynomial chaos expansion for uncertainty quantification of stochastic problems. The performance of the developed scheme becomes more apparent for the problems with larger stochastic dimensions and those requiring higher polynomial order for the stochastic discretization. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
6.
D. Lucor C.‐H. Su G. E. Karniadakis 《International journal for numerical methods in engineering》2004,60(3):571-596
We present a new approach to obtain solutions for general random oscillators using a broad class of polynomial chaos expansions, which are more efficient than the classical Wiener–Hermite expansions. The approach is general but here we present results for linear oscillators only with random forcing or random coefficients. In this context, we are able to obtain relatively sharp error estimates in the representation of the stochastic input as well as the solution. We have also performed computational comparisons with Monte Carlo simulations which show that the new approach can be orders of magnitude faster, especially for compact distributions. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
7.
É. Savin B. Faverjon 《International journal for numerical methods in engineering》2017,111(12):1192-1200
Because of the complexity of fluid flow solvers, non‐intrusive uncertainty quantification techniques have been developed in aerodynamic simulations in order to compute the quantities of interest required in an optimization process, for example. The objective function is commonly expressed in terms of moments of these quantities, such as the mean, standard deviation, or even higher‐order moments. Polynomial surrogate models based on polynomial chaos expansions have often been implemented in this respect. The original approach of uncertainty quantification using polynomial chaos is however intrusive. It is based on a Galerkin‐type formulation of the model equations to derive the governing equations for the polynomial expansion coefficients. Third‐order, indeed fourth‐order moments of the polynomials are needed in this analysis. Besides, both intrusive and non‐intrusive approaches call for their computation provided that higher‐order moments of the quantities of interest need to be post‐processed. In most applications, they are evaluated by Gauss quadratures and eventually stored for use throughout the computations. In this paper, analytical formulas are rather considered for the moments of the continuous polynomials of the Askey scheme, so that they can be evaluated by quadrature‐free procedures instead. Matlab© codes have been developed for this purpose and tested by comparisons with Gauss quadratures. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
8.
Christophe Audouze Prasanth B. Nair 《International journal for numerical methods in engineering》2012,92(4):370-398
In this paper, we consider the problem of constructing reduced‐order models of a class of time‐dependent randomly parametrized linear partial differential equations. Our objective is to efficiently construct a reduced basis approximation of the solution as a function of the spatial coordinates, parameter space, and time. The proposed approach involves decomposing the solution in terms of undetermined spatial and parametrized temporal basis functions. The unknown basis functions in the decomposition are estimated using an alternating iterative Galerkin projection scheme. Numerical studies on the time‐dependent randomly parametrized diffusion equation are presented to demonstrate that the proposed approach provides good accuracy at significantly lower computational cost compared with polynomial chaos‐based Galerkin projection schemes. Comparison studies are also made against Nouy's generalized spectral decomposition scheme to demonstrate that the proposed approach provides a number of computational advantages. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
9.
Tittu Varghese Mathew Jayamanideep Rebbagondla Sundararajan Natarajan 《International journal for numerical methods in engineering》2020,121(19):4329-4346
Nouy and Clement introduced the stochastic extended finite element method to solve linear elasticity problem defined on random domain. The material properties and boundary conditions were assumed to be deterministic. In this work, we extend this framework to account for multiple independent input uncertainties, namely, material, geometry, and external force uncertainties. The stochastic field is represented using the polynomial chaos expansion. The challenge in numerical integration over multidimensional probabilistic space is addressed using the pseudo-spectral Galerkin method. Thereafter, a sensitivity analysis based on Sobol indices using the derived stochastic extended Finite Element Method solution is presented. The efficiency and accuracy of the proposed novel framework against conventional Monte Carlo methods is elucidated in detail for a few one and two dimensional problems. 相似文献
10.
11.
Debraj Ghosh Charbel Farhat 《International journal for numerical methods in engineering》2008,74(8):1219-1239
This paper focuses on the computation of statistical moments of strains and stresses in a random system model where uncertainty is modeled by a stochastic finite element method based on the polynomial chaos expansion. It identifies the cases where this objective can be achieved by analytical means using the orthogonality property of the chaos polynomials and those where it requires a numerical integration technique. To this effect, the applicability and efficiency of several numerical integration schemes are considered. These include the Gauss–Hermite quadrature with the direct tensor product—also known as the Kronecker product—Smolyak's approximation of such a tensor product, Monte Carlo sampling, and the Latin Hypercube sampling method. An algorithm for reducing the dimensionality of integration under a direct tensor product is also explored for optimizing the computational cost and complexity. The convergence rate and algorithmic complexity of all of these methods are discussed and illustrated with the non‐deterministic linear stress analysis of a plate. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
12.
Wael Slika George Saad 《International journal for numerical methods in engineering》2017,112(12):1869-1885
The polynomial chaos Kalman filter (PCKF) has been gaining popularity as a computationally efficient and robust alternative to sampling methods in sequential data assimilation settings. The PCKF's sampling free scheme and attractive structure to represent non‐Gaussian uncertainties makes it a promising approach for data filtering techniques in nonlinear and non‐Gaussian frameworks. However, the accuracy of PCKF is dependent on the dimension and order of the polynomial chaos expansion used to represent all sources of uncertainty in the system. Thus, when independent sources of errors, like process noise and time independent sensors' errors are incorporated in the system, the curse of dimensionality hinders the efficiency and the applicability of PCKF. This study sheds light on this issue and presents a practical framework to maintain an acceptable accuracy of PCKF without scarifying the computational efficiency of the filter. The robustness and efficiency of the presented implementation is demonstrated on 3 typical numerical examples to illustrate its ability to achieve considerable accuracy at a low computational tax. 相似文献
13.
This article aims to present a combination of stochastic finite element and spectral finite element methods as a new numerical tool for uncertainty quantification. One of the well-established numerical methods for reliability analysis of engineering systems is the stochastic finite element method. In this article, a commonly used version of the stochastic finite element method is combined with the spectral finite element method. Furthermore, the spectral finite element method is a numerical method employing special orthogonal polynomials (e.g., Lobatto) and quadrature schemes (e.g., Gauss-Lobatto-Legendre), leading to suitable accuracy, and much less domain discretization with excellent convergence as well. The proposed method of this article is a hybrid method utilizing efficiencies of both methods for analysis of stochastically linear elastostatic problems. Moreover, a spectral finite element method is proposed for numerical solution of a Fredholm integral equation followed by the present method, to provide further efficiencies to accelerate stochastic computations. Numerical examples indicate the efficiency and accuracy of the proposed method. 相似文献
14.
Hongguan Zhang Tadahiro Shibutani 《International journal for numerical methods in engineering》2019,118(1):18-37
In this paper, a new method is proposed that extend the classical deterministic isogeometric analysis (IGA) into a probabilistic analytical framework in order to evaluate the uncertainty in shape and aim to investigate a possible extension of IGA in the field of computational stochastic mechanics. Stochastic IGA (SIGA) method for uncertainty in shape is developed by employing the geometric characteristics of the non-uniform rational basis spline and the probability characteristics of polynomial chaos expansions (PCE). The proposed method can accurately and freely evaluate problems of uncertainty in shape caused by deformation of the structural model. Additionally, we use the intrusive formulation approach to incorporate PCE into the IGA framework, and the C++ programming language to implement this analysis procedure. To verify the validity and applicability of the proposed method, two numerical examples are presented. The validity and accuracy of the results are assessed by comparing them to the results obtained by Monte Carlo simulation based on the IGA algorithm. 相似文献
15.
研究包装件参数不确定性对振动可靠性变化的影响,并分析振动可靠性指标对各不确定参数的灵敏度.采用Karhunen-Loeve展开将具有一定谱特征的平稳随机振动表示在标准正态随机变量空间中,应用一阶可靠性方法分析线性包装件振动可靠性指标.考虑缓冲材料弹性特性、阻尼特性、产品主体和脆弱部件之间的弹性特性、阻尼特性四个随机参数... 相似文献
16.
Xiaodong Wang Charles Hirsch Zhiyi Liu Shun Kang Chris Lacor 《International journal for numerical methods in engineering》2013,94(2):111-127
The aerodynamic performance of a compressor is highly sensitive to uncertain working conditions. This paper presents an efficient robust aerodynamic optimization method on the basis of nondeterministic computational fluid dynamic (CFD) simulation and multi‐objective genetic algorithm (MOGA). A nonintrusive polynomial chaos method is used in conjunction with an existing well‐verified CFD module to quantify the uncertainty propagation in the flow field. This method is validated by comparing with a Monte Carlo method through full 3D CFD simulations on an axial compressor (National Aeronautics and Space Administration rotor 37). On the basis of the validation, the nondeterministic CFD is coupled with a surrogate‐based MOGA to search for the Pareto front. A practical engineering application is implemented to the robust aerodynamic optimization of rotor 37 under random outlet static pressure. Two curve angles and two sweep angles at tip and hub are used as design variables. Convergence analysis shows that the surrogate‐based MOGA can obtain the Pareto front properly. Significant improvements of both mean and variance of the efficiency are achieved by the robust optimization. The comparison of the robust optimization results with that of the initial design, and a deterministic optimization demonstrate that the proposed method can be applied to turbomachinery successfully. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
17.
由于加工、制造等原因,实际结构系统往往所具有很多不确定性,准确评估随机系统的动力学行为不仅具有实际意义,而且是近年来结构动力学理论的一个研究热点。本文研究了同时考虑结构模型参数与所受外激励载荷具有不确定性的复合随机振动问题。结构模型参数的不确定性采用随机变量模拟,外激励载荷的不确定性采用随机过程模拟,提出了结构随机振动响应评估的混合混沌多项式-虚拟激励(PC-PEM)方法。数值算例研究了参数不确定性在21杆桁架中的传播,讨论了响应的一阶、二阶统计矩,并同蒙特卡洛方法进行对比表明提出方法的正确性和有效性。本文的工作对于考虑不确定的复杂装备与结构系统的随机振动分析具有很好的借鉴意义。 相似文献
18.
Kai Cheng Zhenzhou Lu 《International journal for numerical methods in engineering》2020,121(9):2068-2085
In this article, hierarchical surrogate model combined with dimensionality reduction technique is investigated for uncertainty propagation of high-dimensional problems. In the proposed method, a low-fidelity sparse polynomial chaos expansion model is first constructed to capture the global trend of model response and exploit a low-dimensional active subspace (AS). Then a high-fidelity (HF) stochastic Kriging model is built on the reduced space by mapping the original high-dimensional input onto the identified AS. The effective dimensionality of the AS is estimated by maximum likelihood estimation technique. Finally, an accurate HF surrogate model is obtained for uncertainty propagation of high-dimensional stochastic problems. The proposed method is validated by two challenging high-dimensional stochastic examples, and the results demonstrate that our method is effective for high-dimensional uncertainty propagation. 相似文献
19.
Velamur Asokan Badri Narayanan Nicholas Zabaras 《International journal for numerical methods in engineering》2004,60(9):1569-1593
An adjoint‐based functional optimization technique in conjunction with the spectral stochastic finite element method is proposed for the solution of an inverse heat conduction problem in the presence of uncertainties in material data, process conditions and measurement noise. The ill‐posed stochastic inverse problem is restated as a conditionally well‐posed L2 optimization problem. The gradient of the objective function is obtained in a distributional sense by defining an appropriate stochastic adjoint field. The L2 optimization problem is solved using a conjugate‐gradient approach. Accuracy and effectiveness of the proposed approach is appraised with the solution of several stochastic inverse heat conduction problems. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
20.
In the existing probabilistic hierarchical optimization approaches, such as probabilistic analytical target cascading (PATC), all the stochastic interrelated responses are characterized only by the first two statistical moments. However, due to the high nonlinear relation between the inputs and outputs, the interrelated responses are not necessarily normally distributed. The existing approaches, therefore, may not accurately quantify the probabilistic characteristics of the interrelated responses, and would further prevent achieving the real optimal solution. To overcome this deficiency, a novel PATC approach, named PATC-PCE is developed. By employing the polynomial chaos expansion (PCE) technique, the entire distribution of interrelated response can be characterized by a PCE coefficients vector, and then matched and propagated in the hierarchy. Comparative studies show that PATC-PCE outperforms PATC in terms of yielding more accurate optimal solutions and fewer design cycles when the interrelated response are random non-normal quantities, while at a sacrifice of extra function evaluations. 相似文献