共查询到20条相似文献,搜索用时 15 毫秒
1.
A. E. J. Bogaers S. Kok A. G. Malan 《International journal for numerical methods in engineering》2011,86(8):935-952
In this paper, we implement the method of proper orthogonal decomposition (POD) to generate a reduced order model (ROM) of an optimization‐based mesh movement scheme. In this study it is shown that POD can be used effectively to generate an ROM, that accurately reproduces the full order mesh movement algorithm, with a decrease in computational time of over 99%. We further introduce a novel training procedure whereby the POD models are generated in a fully automated fashion. The technology is applicable to any mesh movement method and enables potential reductions of up to four orders of magnitude in mesh movement related costs. The proposed model can be implemented without having to pre‐train the POD model, to any fluid–structure interaction code with an existing mesh movement scheme. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
2.
战庆亮;刘鑫;白春锦;葛耀君 《工程力学》2025,42(6):1-10, 19
工程流动问题中,在结构转角附近一般会发生流动分离,对下游流场的发展有较大影响。因此,对转角处及其附近流场的准确模拟比较重要,尤其是明确转角处的网格划分策略对空间离散精度的影响有较大的研究意义。该文对非结构化网格中的控制方程离散方法进行了分析,结合高雷诺数方柱绕流算例,对非结构化网格离散时的边界层划分策略、网格质量以及结果准确性等问题进行了研究。结果表明:保持转角附近边界层网格的正交性、偏斜度和长宽比对整体流场的模拟结果具有较大影响,并通过流场特征比较给出了转角附近非结构化边界层网格的建议划分方式,可为类似形状工程问题的网格划分与流场模拟提供一定的参考。 相似文献
3.
A novel Lagrangian gradient smoothing method (L‐GSM) is developed to solve “solid‐flow” (flow media with material strength) problems governed by Lagrangian form of Navier‐Stokes equations. It is a particle‐like method, similar to the smoothed particle hydrodynamics (SPH) method but without the so‐called tensile instability that exists in the SPH since its birth. The L‐GSM uses gradient smoothing technique to approximate the gradient of the field variables, based on the standard GSM that was found working well with Euler grids for general fluids. The Delaunay triangulation algorithm is adopted to update the connectivity of the particles, so that supporting neighboring particles can be determined for accurate gradient approximations. Special techniques are also devised for treatments of 3 types of boundaries: no‐slip solid boundary, free‐surface boundary, and periodical boundary. An advanced GSM operation for better consistency condition is then developed. Tensile stability condition of L‐GSM is investigated through the von Neumann stability analysis as well as numerical tests. The proposed L‐GSM is validated by using benchmarking examples of incompressible flows, including the Couette flow, Poiseuille flow, and 2D shear‐driven cavity. It is then applied to solve a practical problem of solid flows: the natural failure process of soil and the resultant soil flows. The numerical results are compared with theoretical solutions, experimental data, and other numerical results by SPH and FDM to evaluate further L‐GSM performance. It shows that the L‐GSM scheme can give a very accurate result for all these examples. Both the theoretical analysis and the numerical testing results demonstrate that the proposed L‐GSM approach restores first‐order accuracy unconditionally and does not suffer from the tensile instability. It is also shown that the L‐GSM is much more computational efficient compared with SPH, especially when a large number of particles are employed in simulation. 相似文献
4.
A. Rassineux P. Breitkopf P. Villon 《International journal for numerical methods in engineering》2003,57(3):371-389
We present a method to adapt a tetrahedron mesh together with a surface mesh with respect to a size criterion. The originality of our work lies in the fact that both surface and tetrahedron mesh adaptation are carried out simultaneously and that no CAD is required to adapt the surface mesh. The adaptation procedure consists of splitting or removing interior and surface edges which violate a given size criterion. The enrichment process is based on a bisection technique. In order to guarantee mesh conformity during the refinement process, all possible remeshing configurations of tetrahedra have been examined. Once the tetrahedron mesh has been adapted, surface nodes are projected on a geometrical model. The building of a surface model from discrete data has already been presented in this journal. The method is based on a mesh‐free technique called Hermite Diffuse Interpolation. Surface and volume mesh optimization procedures are carried out during the adaptation and at the end of the process to enhance the mesh. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
5.
Nestor A. Calvo Sergio R. Idelsohn 《International journal for numerical methods in engineering》2002,55(12):1439-1449
There has been some degree of success in all‐hexahedral meshing. Standard methods start with the object geometry defined by means of an all‐quadrilateral mesh, followed by the use of the combinatorial dual to the mesh in order to define the internal connectivities among elements. For all of the known methods using the dual concept, it is necessary to first prevent or eliminate self‐intersecting (SI) dual lines of the given quadrilateral mesh. The relevant features of SI lines are studied, giving a method to remove them, which avoids deforming the original geometry. Some examples of resulting meshes are shown where the current meshing method has been successfully applied. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
6.
Md. Nurtaj Hossain Debraj Ghosh 《International journal for numerical methods in engineering》2020,121(15):3417-3441
Multiquery problems such as uncertainty quantification (UQ), optimization of a dynamical system require solving a differential equation at multiple parameter values. Therefore, for large systems, the computational cost becomes prohibitive. This issue can be addressed by using a cheaper reduced order model (ROM) instead. However, the ROM entails error in the solution due to approximation in a lower dimensional subspace. Moreover, the ROM lacks robustness over a wide range of parameter values. To address these issues, first, an upper bound on the norm of the state transition matrix is derived. This bound, along with the residual in the governing equation, are then used to develop an error estimator for general nonlinear dynamical systems. Furthermore, this error estimator is used in conjunction with the modified greedy search algorithm proposed by Hossain and Ghosh (Int J Numer Methods Eng, 2018;116(12-13): 741-758) to adaptively construct a robust proper orthogonal decomposition-based ROM. This adaptive ROM is subsequently deployed for UQ by invoking it in a statistical simulation. Two numerical studies: (i) viscous Burgers' equation and (ii) beam on nonlinear Winkler foundation, showed an improved accuracy of the error estimator compared to the current literature. A significant computational speed-up in UQ is achieved. 相似文献
7.
Azam Moosavi Răzvan Ştefănescu Adrian Sandu 《International journal for numerical methods in engineering》2018,113(3):512-533
This paper introduces multivariate input‐output models to predict the errors and bases dimensions of local parametric Proper Orthogonal Decomposition reduced‐order models. We refer to these mappings as the multivariate predictions of local reduced‐order model characteristics (MP‐LROM) models. We use Gaussian processes and artificial neural networks to construct approximations of these multivariate mappings. Numerical results with a viscous Burgers model illustrate the performance and potential of the machine learning‐based regression MP‐LROM models to approximate the characteristics of parametric local reduced‐order models. The predicted reduced‐order models errors are compared against the multifidelity correction and reduced‐order model error surrogates methods predictions, whereas the predicted reduced‐order dimensions are tested against the standard method based on the spectrum of snapshots matrix. Since the MP‐LROM models incorporate more features and elements to construct the probabilistic mappings, they achieve more accurate results. However, for high‐dimensional parametric spaces, the MP‐LROM models might suffer from the curse of dimensionality. Scalability challenges of MP‐LROM models and the feasible ways of addressing them are also discussed in this study. 相似文献
8.
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. 相似文献
9.
Bangti Jin 《International journal for numerical methods in engineering》2008,76(2):230-252
This paper presents two techniques, i.e. the proper orthogonal decomposition (POD) and the stochastic collocation method (SCM), for constructing surrogate models to accelerate the Bayesian inference approach for parameter estimation problems associated with partial differential equations. POD is a model reduction technique that derives reduced‐order models using an optimal problem‐adapted basis to effect significant reduction of the problem size and hence computational cost. SCM is an uncertainty propagation technique that approximates the parameterized solution and reduces further forward solves to function evaluations. The utility of the techniques is assessed on the non‐linear inverse problem of probabilistically calibrating scalar Robin coefficients from boundary measurements arising in the quenching process and non‐destructive evaluation. A hierarchical Bayesian model that handles flexibly the regularization parameter and the noise level is employed, and the posterior state space is explored by the Markov chain Monte Carlo. The numerical results indicate that significant computational gains can be realized without sacrificing the accuracy. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
10.
Weimin Han Gregory J. Wagner Wing Kam Liu 《International journal for numerical methods in engineering》2002,53(6):1323-1336
Implementation of Dirichlet boundary conditions in mesh‐free methods is problematic. In Wagner and Liu (International Journal for Numerical Methods in Engineering 2001; 50 :507), a hierarchical enrichment technique is introduced that allows a simple implementation of the Dirichlet boundary conditions. In this paper, we provide some error analysis for the hierarchical enrichment mesh‐free technique. We derive optimal order error estimates for the hierarchical enrichment mesh‐free interpolants. For one‐dimensional elliptic boundary value problems, we can directly apply the interpolation error estimates to obtain error estimates for the mesh‐free solutions. For higher‐dimensional problems, derivation of error estimates for the mesh‐free solutions depends on the availability of an inverse inequality. Numerical examples in 1D and 2D are included showing the convergence behaviour of mesh‐free interpolants and mesh‐free solutions when the hierarchical enrichment mesh‐free technique is employed. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
11.
定日镜结构风振响应中有多个振型参与,利用完全二次型组合法(CQC)可以考虑各振型间的耦合作用,得到精确结果,但计算量很大,因此提出将随机风场本征正交分解(POD)技术运用到定日镜结构响应计算中.首先给出了POD法原理,讨论了如何将POD法与CQC法结合进行结构脉动响应的频域分析;然后利用POD法分解风洞试验中多通道电子式压力扫描阀系统同步采集到的定日镜表面风压时程,分析了其本征向量的分布特点,以及各阶模态的贡献情况;最后计算出定日镜风致位移响应的均方根,并且对比采用不同阶数模态缩减得到均方根响应的误差.结果表明将POD法引入到风振响应分析中可以大大减少计算工作量并且能保证足够的精度. 相似文献
12.
G. Kerschen J. C. Golinval 《International journal for numerical methods in engineering》2004,60(13):2147-2164
The objective of this paper is to present a model updating strategy of non‐linear vibrating structures. Because modal analysis is no longer helpful in non‐linear structural dynamics, a special attention is devoted to the features extracted from the proper orthogonal decomposition and one of its non‐linear generalizations based on auto‐associative neural networks. The efficiency of the proposed procedure is illustrated using simulated data from a three‐dimensional portal frame. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
13.
A new formulation of the element‐free Galerkin (EFG) method is presented in this paper. EFG has been extensively popularized in the literature in recent years due to its flexibility and high convergence rate in solving boundary value problems. However, accurate imposition of essential boundary conditions in the EFG method often presents difficulties because the Kronecker delta property, which is satisfied by finite element shape functions, does not necessarily hold for the EFG shape function. The proposed new formulation of EFG eliminates this shortcoming through the moving kriging (MK) interpolation. Two major properties of the MK interpolation: the Kronecker delta property (?I( s J)=δIJ) and the consistency property (∑In?I( x )=1 and ∑In?I( x )xIi=xi) are proved. Some preliminary numerical results are given. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
14.
Annika Radermacher Stefanie Reese 《International journal for numerical methods in engineering》2016,107(6):477-495
The constantly rising demands on finite element simulations yield numerical models with increasing number of degrees‐of‐freedom. Due to nonlinearity, be it in the material model or of geometrical nature, the computational effort increases even further. For these reasons, it is today still not possible to run such complex simulations in real time parallel to, for example, an experiment or an application. Model reduction techniques such as the proper orthogonal decomposition method have been developed to reduce the computational effort while maintaining high accuracy. Nonetheless, this approach shows a limited reduction in computational time for nonlinear problems. Therefore, the aim of this paper is to overcome this limitation by using an additional empirical interpolation. The concept of the so‐called discrete empirical interpolation method is translated to problems of solid mechanics with soft nonlinear elasticity and large deformations. The key point of the presented method is a further reduction of the nonlinear term by an empirical interpolation based on a small number of interpolation indices. The method is implemented into the finite element method in two different ways, and it is extended by using different solution strategies including a numerical as well as a quasi‐Newton tangent. The new method is successfully applied to two numerical examples concerning hyperelastic as well as viscoelastic material behavior. Using the extended discrete empirical interpolation method combined with a quasi‐Newton tangent enables reductions in computational time of factor 10 with respect to the proper orthogonal decomposition method without empirical interpolation. Negligibly, orders of error can be reached. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
15.
将特征正交分解型谱表示法用于模拟汽车受路面激励。首先给出了路面不平度对汽车输入的位移随机激励的谱描述。基于路面激励的功率谱矩阵,结合特征正交分解(POD,Proper Orthogonal Decomposition,)型谱表示法的模拟表达式,给出了路面激励的显式POD分解,定义了汽车的“路面激励模态”,推导了路面对汽车输入激励随机模拟的简化计算公式。该方法可用FFT来减少计算量以提高计算速度。它由于完全消除掉了原型谱表示法的Cholesky分解过程而具有较高的计算效率和更明确的物理意义。最后,通过对一个四轮轿车在国标GB7031—87中的A级路面不平度下受到的位移随机激励进行模拟,说明了该方法的有效性。 相似文献
16.
Kevin Carlberg Charbel Farhat 《International journal for numerical methods in engineering》2011,86(3):381-402
A novel model reduction technique for static systems is presented. The method is developed using a goal‐oriented framework, and it extends the concept of snapshots for proper orthogonal decomposition (POD) to include (sensitivity) derivatives of the state with respect to system input parameters. The resulting reduced‐order model generates accurate approximations due to its goal‐oriented construction and the explicit ‘training’ of the model for parameter changes. The model is less computationally expensive to construct than typical POD approaches, since efficient multiple right‐hand side solvers can be used to compute the sensitivity derivatives. The effectiveness of the method is demonstrated on a parameterized aerospace structure problem. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
17.
Chung‐Hao Lee Jiun‐Shyan Chen 《International journal for numerical methods in engineering》2013,96(10):599-627
Model order reduction for molecular dynamics (MD) systems exhibits intrinsic complexities because of the highly nonlinear and nonlocal multi‐atomic interactions in high dimensions. In the present work, we introduce a proper orthogonal decomposition‐based method in conjunction with the radial basis function (RBF) approximation of the nonlinear and nonlocal potential energies and inter‐atomic forces for MD systems. This approach avoids coordinate transformation between the physical and reduced‐order coordinates, and allows the potentials and inter‐atomic forces to be calculated directly in the reduced‐order space. The RBF‐approximated potential energies and inter‐atomic forces in the reduced‐order space are discretized on the basis of the Smolyak sparse grid algorithm to further enhance the effectiveness of the proposed method. The good approximation properties of RBFs in interpolating scattered data make them ideal candidates for the reduced‐order approximation of MD inter‐atomic force calculations. The proposed approach is validated by performing the reduced‐order simulations of DNA molecules under various external loadings. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
18.
19.
推导了本征正交分解(Proper Orthogonal Decomposition,POD)型谱表示法模拟所得平稳正态脉动风场的偏度误差和随机误差.从POD型谱表示法的模拟公式出发,推导了Ⅳ变量风场模拟结果序列的样本均值、相关函数、功率谱函数和根方差等前二阶矩统计特征的时域估计表达式;并证明了时域估计相关函数是正态过程,功率谱函数为非正态随机过程.进一步,计算上述样本时域估计二阶矩特征的均值和根方差,即得到了POD型谱表示法模拟所得风场的各统计量时域估计的偏度误差和随机误差,并以此给出了误差计算的通式.算例中统计误差和理论误差值的对比验证了所推导的解析解. 相似文献
20.
Yan Liu Ted Belytschko 《International journal for numerical methods in engineering》2010,82(6):699-715
A new support integration technique is proposed, which is similar to those used in truly mesh‐free methods. The contribution of this paper is to exploit the divergence‐free condition for the support integrals to construct quadrature formulas that only require three integration points per particle in two dimensions. Numerical examples show that the proposed integration method can achieve results that agree with manufactured closed‐form solutions. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献