Comparison of sequential data assimilation methods for the Kuramoto–Sivashinsky equation

The Kuramoto–Sivashinsky equation plays an important role as a low‐dimensional prototype for complicated fluid dynamics systems having been studied due to its chaotic pattern forming behavior. Up to now, efforts to carry out data assimilation with this 1‐D model were restricted to variational adjoint methods domain and only Chorin and Krause (Proc. Natl. Acad. Sci. 2004; 101 (42):15013–15017) tested it using a sequential Bayesian filter approach. In this work we compare three sequential data assimilation methods namely the Kalman filter approach, the sequential Monte Carlo particle filter approach and the maximum likelihood ensemble filter methods. This comparison is to the best of our knowledge novel. We compare in detail their relative performance for both linear and nonlinear observation operators. The results of these sequential data assimilation tests are discussed and conclusions are drawn as to the suitability of these data assimilation methods in the presence of linear and nonlinear observation operators. Copyright © 2009 John Wiley & Sons, Ltd.     

Some theoretical problems on variational data assimilation

Theoretical aspects of variational data assimilation (VDA) for a simple model with both global and local observational data are discussed. For the VDA problems with global observational data, the initial conditions and parameters for the model are revisited and the model itself is modified. The estimates of both error and convergence rate are theoretically made and the validity of the method is proved. For VDA problem with local observation data, the conventional VDA method are out of use due to the ill-posedness of the problem. In order to overcome the difficulties caused by the ill-posedness, the initial conditions and parameters of the model are modified by using the improved VDA method, and the estimates of both error and convergence rate are also made. Finally, the validity of the improved VDA method is proved through theoretical analysis and illustrated with an example, and a theoretical criterion of the regularization parameters is proposed.     

A variational ensemble Kalman filtering method for data assimilation using 2D and 3D version of COHERENS model

Decoupled implementation of data assimilation methods has been rarely studied. The variational ensemble Kalman filter has been implemented such that it needs not to communicate directly with the model, but only through input and output devices. In this work, an open multi‐functional three‐dimensional (3D) model, the coupled hydrodynamical‐ecological model for regional and shelf seas (COHERENS), has been used. Assimilation of the total suspended matter (TSM) is carried out in 154 km² lake Säkylän Pyhäjärvi. Observations of TSM were derived from high‐resolution satellite images of turbidity and chrolophyll‐a. For demonstrating the method, we have used a low‐resolution model grid of 1 km. The model was run for a period from May 16 to September 14. We have run the COHERENS model with two‐dimensional (2D) mode time steps and 3D mode time steps. This allows COHERENS to switch between 2D and 3D modes in a single run for computational efficiency. We have noticed that there is not much difference between these runs. This is because satellite images depict the derived TSM for the surface layer only. The use of additional 3D data might change this conclusion and improve the results. We have found that in this study, the use of a large ensemble size does not guarantee higher performance. The successful implementation of decoupled variational ensemble Kalman filter method opens the way for other methods and evolution models to enjoy the benefits without having to spend substantial effort in merging the model and assimilation codes together, which can be a difficult task. © 2016 The Authors. International Journal for Numerical Methods in Fluids Published by John Wiley & Sons Ltd.     

An efficient very large eddy simulation model for simulation of turbulent flow

Among the various hybrid methodologies, Speziale's very large eddy simulation (VLES) is one that was proposed very early. It is a unified simulation approach that can change seamlessly from Reynolds Averaged Navier–Stokes (RANS) to direct numerical simulation (DNS) depending on the numerical resolution. The present study proposes a new improved variant of the original VLES model. The advantages are achieved in two ways: (i) RANS simulation can be recovered near the wall which is similar to the detached eddy simulation concept; (ii) a LES subgrid scale model can be reached by the introduction of a third length scale, that is, the integral turbulence length scale. Thus, the new model can provide a proper LES mode between the RANS and DNS limits. This new methodology is implemented in the standard k ? ? model. Applications are conducted for the turbulent channel flow at Reynolds number of Re_τ = 395, periodic hill flow at Re = 10,595, and turbulent flow past a square cylinder at Re = 22,000. In comparison with the available experimental data, DNS or LES, the new VLES model produces better predictions than the original VLES model. Furthermore, it is demonstrated that the new method is quite efficient in resolving the large flow structures and can give satisfactory predictions on a coarse mesh. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical study of an oscillatory turbulent flow over a flat plate

Oscillatory turbulent flow over a flat plate is studied using large eddy simulation (LES) and Reynolds-average Navier-Stokes (RANS) methods. A dynamic subgrid-scale model is employed in LES and Saffman's turbulence model is used in RANS. The flow behaviors are discussed for the accelerating and decelerating phases during the oscillating cycle. The friction force on the wall and its phase shift from laminar to turbulent regime are also investigated for different Reynolds numbers. The project supported by the Youngster Funding of Academia Sinica and by the National Natural Science Foundation of China     

基于卡尔曼滤波的信息融合算法优化研究

通过比较采用联邦卡尔曼滤波的状态向量融合和量测信息融合,得出量测信息融合优于状态向量融合,因为只有当卡尔曼滤波一致时状态向量融合才有效.采用基于最小均方差估计的观测值加权融合法融合了多传感器数据,保持了观测向量的维数.这种方法具有高效性.为了提高该算法的速度和精度,对系统的量测空间进行了等价变换,而等价系统的状态空间却没有改变.给出了等价变换前后的系统误差方差阵和状态估计均一致性的证明.把矩阵分析中的L-D分解算法运用到该算法中以避免计算矩阵的逆,从而改善了算法的稳定性和精度.举例验证了所设计算法的这些优点,给出了采用联邦卡尔曼滤波和所优化滤波算法的状态估计和误差的仿真结果,并依次进行了分析.经过这种优化,算法的精度和速度得到很大提高,已经应用到实际工程中.     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

Input-output reduced-order modeling of unsteady flow over an airfoil at a high angle of attack based on dynamic mode decomposition with control

Due to the damage caused by stall flutter, the investigation and modeling of the flow over a wind turbine airfoil at high angles of attack are essential. Dynamic mode decomposition (DMD) and dynamic mode decomposition with control (DMDc) are used to analyze unsteady flow and identify the intrinsic dynamics. The DMDc algorithm is found to have an identification problem when the spatial dimension of the training data is larger than the number of snapshots. IDMDc, a variant algorithm based on reduced dimension data, is introduced to identify the precise intrinsic dynamics. DMD, DMDc and IDMDc are all used to decompose the data for unsteady flow over the S809 airfoil that are obtained by numerical simulations. The DMD results show that the dominant feature of a static airfoil is the adjacent shedding vortices in the wake. For an oscillating airfoil, the DMDc results may fail to consider the effect of the input and have an identification problem. IDMDc can alleviate this problem. The dominant IDMDc modes show that the intrinsic flow for the oscillating case is similar to the unsteady flow over the static airfoil. Moreover, the input–output model identified by IDMDc can give better predictions for different oscillating cases than the identified DMDc model.     

联合直接攻击弹药精确制导技术分析

从技术分析的角度对联合直接攻击弹药(JDAM)的武器制导系统进行了研究。安装JDAM制导控制装置的弹药武器,通过GPS/INS组合制导和尾翼调节,控制武器精确命中目标。在对JDAM系统组成和精度分析的基础上,详细讨论了JDAM系统的卡尔曼滤波器、捷联惯性导航技术、传递对准和GPS辅助惯性制导等技术。对设计中遇到主要的技术难点和相应的解决途径也进行了分析。     

Large eddy simulation of separation control over a backward-facing step flow by suction

Separation control over a backward-facing step (BFS) flow by continuous suction was numerically investigated using the turbulence model of large eddy simulation (LES). The effect of suction control on the flow fields was scrutinised by altering the suction flow coefficient, and the results indicate that suction is not only very effective in shortening the reattachment length but also very influential in reducing the tangential velocity gradient and turbulence fluctuations of the reattached flows. With increasing increments of the absolute suction flow coefficient, the effect of suction control is more significant. Furthermore, the detailed flow fields (including the time-averaged stream and velocity fields) and turbulence characteristics (including the time-averaged resolved kinetic energy and RMS velocity) for the BFS models with or without suction are presented to discuss the mechanism of suction control. Comparisons of the time-averaged statistics between the numerical simulations and corresponding experiments are conducted, and it shows that the LES based on the dynamic kinetic energy subgrid-scale model (DKEM) can acquire exact results. Therefore, feasibility of the numerical methods to simulate suction-controlled models is validated.     

Recent progress on the study of asymmetric vortex flow over slender bodies

The investigations of forebody vortex flow and its flow control have great importance in both academic field and engineering application areas. A large number of papers and many review papers have been published. However in this research field of forebody asymmetric vortices, three problems such as tip perturbation effect, Reynolds number effect and flow instability are less studied and thus not understood completely. So many researches are still working on the issues in recent years. The present paper attempts to provide a review of recent research progress on first two problems. The first problem is mainly concerned with how the vortex flow evolves after tip perturbation; how to solve the problem of repeatability and reproducibility of wind tunnel testing data; how to develop a conception of active flow control technique with tip perturbation based on the study of vortex flow response to tip perturbation. For the second problem one is mainly concerned that how the asymmetric vortices are developed with the increase of Reynolds number; how to classify the vortex flow patterns in different Reynolds number regimes; how to develop an appropriate boundary layer transition technique to simulate flows at high Reynolds number in the convention wind tunnels. Finally, some important ques- tions that deserve answers are proposed in the concluding remarks.     

An investigation of thermally stratified turbulent channel flow with temperature oscillation on the bottom wall by large eddy simulation

Thermally stratified shear turbulent channel flow with temperature oscillation on the bottom wall of the channel is calculated to investigate the behavior of turbulent flow and heat transfer by use of large eddy simulation (LES) approach coupled with dynamic subgrid-scale (SGS) models. The objective of this study is to deal with the effect of the temperature oscillation on turbulent behavior of thermally stratified turbulent channel flow and to examine the effectiveness of the LES technique for predicting statistically unsteady turbulent flow driven by time-varying buoyancy force. To validate the present calculation, thermally stratified shear turbulent channel flow is computed and compared with available data obtained by direct numerical simulation (DNS), which confirm that the present approach can be used to predict thermally stratified turbulent channel flow satisfactorily. Further, to illustrate the effect of the temperature oscillation with different Richardson numbers and periods of the oscillation on turbulence characteristics, the phase-averaged mean value and fluctuation of the resolved velocities and temperature, and instantaneous velocity fluctuation structures are analyzed.