首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
    
This paper introduces a sparse matrix discrete interpolation method to effectively compute matrix approximations in the reduced order modeling framework. The sparse algorithm developed herein relies on the discrete empirical interpolation method and uses only samples of the nonzero entries of the matrix series. The proposed approach can approximate very large matrices, unlike the current matrix discrete empirical interpolation method, which is limited by its large computational memory requirements. The empirical interpolation indices obtained by the sparse algorithm slightly differ from the ones computed by the matrix discrete empirical interpolation method as a consequence of the singular vectors round‐off errors introduced by the economy or full singular value decomposition (SVD) algorithms when applied to the full matrix snapshots. When appropriately padded with zeros, the economy SVD factorization of the nonzero elements of the snapshots matrix is a valid economy SVD for the full snapshots matrix. Numerical experiments are performed with the 1D Burgers and 2D shallow water equations test problems where the quadratic reduced nonlinearities are computed via tensorial calculus. The sparse matrix approximation strategy is compared against five existing methods for computing reduced Jacobians: (i) matrix discrete empirical interpolation method, (ii) discrete empirical interpolation method, (iii) tensorial calculus, (iv) full Jacobian projection onto the reduced basis subspace, and (v) directional derivatives of the model along the reduced basis functions. The sparse matrix method outperforms all other algorithms. The use of traditional matrix discrete empirical interpolation method is not possible for very large dimensions because of its excessive memory requirements. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

2.
GPS/INS 组合导航系统降阶滤波器设计   总被引:1,自引:0,他引:1  
本文研究了某型号飞机上GPS/平台惯性组合导航系统中降阶卡尔曼滤波器的设计方法;提出一种降阶滤器结构,并讨论了模型误差的伪随机噪声补偿方法。采用蒙特—卡洛分析法对降阶滤波器的实际误差进行了数字仿真,结果表明,本文提出的这种降阶滤波器能保证组合系统的导航精度,而且具有实际使用价值  相似文献   

3.
    
In this article, a reduced‐order modeling approach, suitable for active control of fluid dynamical systems, based on proper orthogonal decomposition (POD) is presented. The rationale behind the reduced‐order modeling is that numerical simulation of Navier–Stokes equations is still too costly for the purpose of optimization and control of unsteady flows. The possibility of obtaining reduced‐order models that reduce the computational complexity associated with the Navier–Stokes equations is examined while capturing the essential dynamics by using the POD. The POD allows the extraction of a reduced set of basis functions, perhaps just a few, from a computational or experimental database through an eigenvalue analysis. The solution is then obtained as a linear combination of this reduced set of basis functions by means of Galerkin projection. This makes it attractive for optimal control and estimation of systems governed by partial differential equations (PDEs). It is used here in active control of fluid flows governed by the Navier–Stokes equations. In particular, flow over a backward‐facing step is considered. Reduced‐order models/low‐dimensional dynamical models for this system are obtained using POD basis functions (global) from the finite element discretizations of the Navier–Stokes equations. Their effectiveness in flow control applications is shown on a recirculation control problem using blowing on the channel boundary. Implementational issues are discussed and numerical experiments are presented. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

4.
    
Model order reduction of the two‐dimensional Burgers equation is investigated. The mathematical formulation of POD/discrete empirical interpolation method (DEIM)‐reduced order model (ROM) is derived based on the Galerkin projection and DEIM from the existing high fidelity‐implicit finite‐difference full model. For validation, we numerically compared the POD ROM, POD/DEIM, and the full model in two cases of Re = 100 and Re = 1000, respectively. We found that the POD/DEIM ROM leads to a speed‐up of CPU time by a factor of O(10). The computational stability of POD/DEIM ROM is maintained by means of a careful selection of POD modes and the DEIM interpolation points. The solution of POD/DEIM in the case of Re = 1000 has an accuracy with error O(10?3) versus O(10?4) in the case of Re = 100 when compared with the high fidelity model. For this turbulent flow, a closure model consisting of a Tikhonov regularization is carried out in order to recover the missing information and is developed to account for the small‐scale dissipation effect of the truncated POD modes. It is shown that the computational results of this calibrated ROM exhibit considerable agreement with the high fidelity model, which implies the efficiency of the closure model used. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

5.
6.
    
A novel proper orthogonal decomposition (POD) model has been developed for use with an advanced unstructured mesh finite‐element ocean model, the Imperial College Ocean Model (ICOM, described in detail below), which includes many recent developments in ocean modelling and numerical analysis. The advantages of the POD model developed here over existing POD approaches are the ability:
  • 1. To increase accuracy when representing geostrophic balance (the balance between the Coriolis terms and the pressure gradient). This is achieved through the use of two sets of geostrophic basis functions where each one is calculated by basis functions for velocities u and v.
  • 2. To speed up the POD simulation. To achieve this a new numerical technique is introduced, whereby a time‐dependent matrix in the discretized equation is rapidly constructed from a series of time‐independent matrices. This development imparts considerable efficiency gains over the often‐used alternative of calculating each finite element over the computational domain at each time level.
  • 3. To use dynamically adaptive meshes in the above POD model.
Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

7.
V. I. Sakharov 《Fluid Dynamics》2007,42(6):1007-1016
The results of numerical simulation are presented for thermally and chemically nonequilibrium air plasma flows in a plasmatron discharge channel and underexpanded dissociated and partially ionized air jets flowing past a cylindrical model with a blunt leading edge and cooled copper surface under the experimental conditions realized in a VGU-4 100 kW induction plasmatron (Institute for Problems in Mechanics of the Russian Academy of Sciences) (see, for example, [1, 2]). The nonequilibrium excitation of the vibrational degrees of freedom of the molecules in the modal approximation and the difference between the electron and translational heavy-particle temperatures are taken into account in the calculations. The calculated data on the heat transfer and pressure at the stagnation point are compared with the results obtained within the framework of the thermally equilibrium model. Comparison with the experimental data obtained in the Institute for Problem in Mechanics of the Russian Academy of Sciences (Laboratory for interaction between plasma and radiation and materials) and kindly provided for comparison purposes gives satisfactory agreement.  相似文献   

8.
    
A novel acceleration technique using a reduced‐order model is presented to speed up convergence of continuous adjoint solvers. The acceleration is achieved by projecting to an improved solution within an iterative process solely using early solution results. This is achieved by forming basis vectors from early iteration adjoint solutions to perform model order reduction of the adjoint equations. The reduced‐order model of the adjoint equations is then substituted into the full‐order discretized governing equations to determine weighting coefficients for each basis vector. With these coefficients, a linear combination of the basis vectors is used to project to an improved solution. The method is applied to 3 inviscid quasi‐1D nozzle flow cases including fully subsonic flow, subsonic inlet to supersonic outlet flow, and transonic flow with a shock. Significant cost reductions are achieved for a single application as well as repeated applications of the convergence acceleration technique.  相似文献   

9.
    
Intrusive reduced order modeling techniques require access to the solver's discretization and solution algorithm, which are not available for most computational fluid dynamics codes. Therefore, a nonintrusive reduction method that identifies the system matrix of linear fluid dynamical problems with a least-squares technique is presented. The methodology is applied to the linear scalar transport convection-diffusion equation for a 2D square cavity problem with a heated lid. The (time-dependent) boundary conditions are enforced in the obtained reduced order model (ROM) with a penalty method. The results are compared and the accuracy of the ROMs is assessed against the full order solutions and it is shown that the ROM can be used for sensitivity analysis by controlling the nonhomogeneous Dirichlet boundary conditions.  相似文献   

10.
    
Inflow and outflow boundary conditions are essential for the application of computational fluid dynamics to many engineering scenarios. In this paper we present a new boundary condition implementation that enables the simulation of flow through permeable boundaries in the Lagrangian mesh‐free method, smoothed particle hydrodynamics (SPH). Each permeable boundary is associated with an inflow or outflow zone outside the domain, in which particles are created or removed as required. The analytic boundary condition is applied by prescribing the appropriate variables for particles in an inflow or outflow zone, and extrapolating other variables from within the domain. Characteristic‐based non‐reflecting boundary conditions, described in the literature for mesh‐based methods, can be implemented within this framework. Results are presented for simple one‐dimensional flows, quasi‐one‐dimensional compressible nozzle flow, and two‐dimensional flow around a cylinder at Reynolds numbers of 40 and 100 and a Mach number of 0.1. These results establish the capability of SPH to model flows through open domains, opening a broad new class of applications. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

11.
A reduced order cyclic method was developed to compute limit-cycle oscillations for large, nonlinear, multidisciplinary systems of equations. Method efficacy was demonstrated for two simplified models: a typical-section airfoil with nonlinear structural coupling and a nonlinear panel in high-speed flow. The cyclic method was verified to maintain second-order temporal accuracy, yield converged limit cycles in about 10 Newton iterates, and provide precise estimates of cycle frequency. This method was projected onto a low-order space using a set of variables governing the amplitudes of empirically derived modes, which were computed with the proper orthogonal decomposition. In this reduced order form, the cyclic Jacobian was greatly compressed, allowing accurate limit cycle solutions to be very efficiently computed.The U.S. Governments right to retain a non-exclusive, royalty-free licence in and to any copyright is acknowledged.  相似文献   

12.
This paper discusses modal filtering of experimental data and the corresponding identification of linear and nonlinear parameters in reduced order space. Specifically, several experimental configurations will be discussed in order to provide insight into such identification issues as spatial discretization, observability, and the linear independence of the assumed filter or basis. The two experiments considered herein represent different measurement configurations of the same clamped–clamped beam. First, asymmetric inertial loading via asymmetric sensor location was considered, while the second scenario presents a symmetric sensor configuration. Several important conclusions can be drawn from the two experimental scenarios. First, by asymmetrically loading the beam, a corresponding asymmetric beam mode was excited yet not observable. In the second scenario, the symmetric distribution of sensors minimized the impact of the respective asymmetric mode. The resulting spatial information allowed for the proper filtering of the remnants of the asymmetric mode. Nonlinear parameters in modal space as well as the underlying linear parameters were successfully identified simultaneously in both experimental scenarios, although the usefulness of the asymmetrically loaded beam was limited. Finally, successful comparisons were made between the identified reduced order model and experimental response at the beam quarter point using the symmetric case and the beam midpoint using both experimental scenarios.  相似文献   

13.
Fluid–structure interaction (FSI) phenomena are of significant importance in several engineering fields. Recently developed active flow control devices regulate the FSI in order to control the dynamic response of the structure that is involved. As a first step to use active control, computationally efficient reduced-order models are required. The reduced-order models must be able to predict the nonlinear structural dynamic response given an incoming flow condition. This paper presents a computationally efficient method for the construction of a hybrid reduced-order model for FSI problems based on data obtained through high-fidelity numerical simulations. The model splits the force and the structural dynamic response into two separate blocks and uses model reduction techniques to account for the flow field information. The current model is tested on a vibrating rigid cylinder submerged in a flow at low Reynolds number regime.  相似文献   

14.
A class of semilinear nonautonomous parabolic equations subjected to additive white noise is considered. The existence of a family of randomN-dimensional approximate inertial manifolds (AIMs) whose neighborhoods of thickness of order exp(- N+1 ) attract exponentially in the mean all the trajectories is proved forN large enough. Here N+1 is the (N+1)th eigenvalue of the corresponding linear problem, and and are positive constants. We also construct a sequence of AIMs which converges to the exact inertial manifold, when a spectral gap condition is satisfied. These results remain true for deterministic autonomous and nonautonomous cases.  相似文献   

15.
16.
    
Reduced‐Order Models (ROMs) have been the focus of research in various engineering situations, but it is only relatively recently that such techniques have begun to be introduced into the CFD field. The purpose of generating such models is to capture the dominant dynamics of the full set of CFD equations, but at much lower cost. One method that has been successfully implemented in the field of fluid flows is based on the calculation of the linear pulse responses of the CFD scheme coupled with an Eigensystem Realization algorithm (ERA), resulting in a compact aerodynamic model. The key to the models is the identification of the linear responses of the non‐linear CFD code. Two different methods have been developed and reported in literature for linear response identification; the first method linearizes the CFD code and the second method uses Volterra theory and the non‐linear code. As these methods were developed independently they have not previously been brought together and compared. This paper first explains the subtle, but fundamental differences between the two methods. In addition, a series of test cases are shown to demonstrate the performance and drawbacks of the ROMs derived from the different linear responses. The conclusions of this study provide useful guidance for the implementation of either of the two approaches to obtain the linear responses of an existing CFD code. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

17.
    
In this paper, a new vector‐filtering criterion for dynamic modes selection is proposed that is able to extract dynamically relevant flow features from dynamic mode decomposition of time‐resolved experimental or numerical data. We employ a novel modes selection criterion in parallel with the classic selection based on modes amplitudes, in order to analyze which of these procedures better highlight the coherent structures of the flow dynamics. Numerical tests are performed on two distinct problems. The efficiency of the proposed criterion is proved in retaining the most influential modes and reducing the size of the dynamic mode decomposition model. By applying the proposed filtering mode technique, the flow reconstruction error is shown to be significantly reduced. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

18.
安效民 《计算力学学报》2014,31(2):273-276,284
传统气动弹性的时域计算耗费了大量时间,为了提高计算效率,本文发展了基于边界元方法的降阶模型技术。首先基于边界元方法建立非定常流场的求解模型,结合特征值分析技术建立了非定常气动力的低阶模型;然后,利用边界元方法建立了气动网格和结构网格之间的信息转换矩阵;最后将非定常气动力降阶模型和结构动力学方程联合,建立了气动弹性系统的低阶状态空间模型。将所发展的降阶模型方法应用于NACA0012翼型的非定常气动力求解中,结果表明降阶模型可以在保证原系统计算精度的同时提高了计算效率;将降阶模型技术应用到三维机翼的气动弹性响应计算中,在系统阶数仅为12阶的情况下可以得到与原系统一致的极限环响应,说明降阶模型技术在求解气动弹性问题中的巨大优势。  相似文献   

19.
    
ABSTRACT

The purpose of this paper is the identification of a reduced order model (ROM) from numerical code output by non-intrusive techniques (i.e. not requiring projecting of the governing equations onto the reduced basis modes). In this paper, we perform a comparison between two methods of model order reduction based on dynamic mode decomposition (DMD). The first method is a deterministic (classic) DMD technique endowed with a dynamic filtering criterion of selection of modes used in the ROM model. The second method is an adaptive randomised DMD algorithm (ARDMD) based on a randomised singular value decomposition. This produced an accelerating algorithm, which is endowed with a few additional advantages. In addition, the reduced order model is guaranteed to satisfy the boundary conditions of the full model, which is crucial for surrogate modelling. For numerical illustration, we use the shallow water equations model.  相似文献   

20.
基于气动力降阶模型的跨音速气动弹性稳定性分析   总被引:6,自引:0,他引:6       下载免费PDF全文
基于离散型输入输出差分模型,运用非定常CFD方法训练信号,然后运用最小二乘方法进行参数辨识,得到降阶的非定常气动力模型,再将该离散差分模型转换为连续时间域内的状态方程.耦合气动状态方程和结构状态方程,得到耦合系统的气动弹性状态方程.求解不同动压下状态矩阵的特征值,根据根轨迹图分析系统的稳定性特性.分析结果与直接耦合CFD/CSD方法结果相吻合,可以计算跨音速非线性气动弹性问题.其计算效率比直接耦合CFD/CSD方法提高1~2个数量级.针对Isogai wing在跨音速出现的S型颤振边界进行了较为细致的分析,阐述了该现象是由于系统诱发颤振的分支随着速度(来流动压)的提高而发生转移所导致的.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号