Optimal control of shallow water flow using time-delay function

This paper presents an optimal control system that includes a time-delay function for application to flood control setups with a retardation area. This system consists of the present and past controls that express flow behaviour in the retardation area. Optimal control theory is used to obtain a control discharge that satisfies the state equation including the time-delay function and minimizes the performance function. The optimal control and the delayed control discharges are obtained by the solution of an adjoint equation. The weighted gradient method is employed as a minimization algorithm. The Galerkin finite element procedure is employed to discretize the state and adjoint equations in the spatial direction. The bubble function interpolation, originated by the authors' group, using a stabilized term, is employed for the discretization in space. The flood flow in the Tsurumi river is presented as a numerical model. We show in this paper that floods can be controlled by means of a time-delay function.     

Computational water purification controls using second‐order adjoint equations with stability confirmation

This paper presents a computational method for water purification using second‐order adjoint equations. In Japan, the waters of polluted rivers are purified by conveying the waters from other rivers into the main rivers or by using outflows from sewage plants. The shallow water flow equation based on the water velocity and elevation and the advection diffusion equation of COD concentration are governing equations. The control problem involves finding a flow velocity into the main river that can reduce the COD concentration as close to the target value as possible. In other words, the problem is to find a water velocity to minimize the performance function, which is the square sum of the discrepancy between the computed and the observed COD concentrations. The present research was motivated by the need to apply water purification controls to practical projects. We have found that the controls occasionally tend to be unstable, and the stability of control must be ensured. By expanding the extended performance function into the Taylor series, the necessary condition for the stationary state is derived. Based on this condition, the first‐and second‐order adjoint equations can be obtained. The backward solution of the adjoint equation leads to the gradient and the Hessian product; these serve as the basis of the quasi‐Newton method. From the condition that the performance function must be minimum, the stability confirmation index can be determined. Using this index, we have derived the trust region method, the computation of which confirms the stability of control. Verification was carried out using a simple channel model. By varying the peak value of the inflow velocity, the outlet velocity has been determined such that the water elevation at the target point is zero. Depending on the peak value of the inflow, unstable control arises; this is determined by the stability confirmation index presented in this paper. The trust region method with the stability confirmation index is shown to be adaptable to judge the stability of control. The present method was applied to the water purification of Teganuma river in Japan. The steady fundamental state was computed with the inflow, outflow, and COD concentration at the inlet being specified. The control velocity at the control point can be determined for a fixed control duration with and without the stability confirmation index. The inflow, outflow, and COD concentration are specified as functions of time. It is shown that this method is suitable for practical use because control stability can be ensured. Moreover, it is also noted that the maximum flow velocity for stable control depending on the given control duration can be obtained. Copyright © 2011 John Wiley & Sons, Ltd.     

研究柔性体撞击问题的子结构离散方法

本文用子结构方法研究刚性小球和均质柔性杆的纵向撞击以及和均质柔性梁的横向撞击问题,导出了用模态坐标表示的动力学方程,通过对刚性小球和柔性杆的纵向撞击的仿真教育处发现,在总单元数相同的情况下,取8个子结构比较合理,在此基础上对刚性小球和均质柔性梁的横向撞击的特性进行研究,发现撞击力在变化过程中会产生上下波动,当梁的弹性模量增加时,撞击力增大,撞击时间缩短。     

求解接触问题的一种新的实验误差法

提出了一种带松弛因子的UZAW算法求解实验误差法中给定状态下的位移和接触力满足的等式方程，并证明了该算法是R超线性收敛的。整个区域被划分为多个子区域，不同子区域位移场的求解是独立的。还提出了一种带参数的以不完全因子分解为基础的预条件子共轭梯度法求解不同子区域位移场，该算法在块体规模较大时更加有效。     

Moving immersed boundary method

A novel implicit immersed boundary method of high accuracy and efficiency is presented for the simulation of incompressible viscous flow over complex stationary or moving solid boundaries. A boundary force is often introduced in many immersed boundary methods to mimic the presence of solid boundary, such that the overall simulation can be performed on a simple Cartesian grid. The current method inherits this idea and considers the boundary force as a Lagrange multiplier to enforce the no‐slip constraint at the solid boundary, instead of applying constitutional relations for rigid bodies. Hence excessive constraint on the time step is circumvented, and the time step only depends on the discretization of fluid Navier‐Stokes equations, like the CFL condition in present work. To determine the boundary force, an additional moving force equation is derived. The dimension of this derived system is proportional to the number of Lagrangian points describing the solid boundaries, which makes the method very suitable for moving boundary problems since the time for matrix update and system solving is not significant. The force coefficient matrix is made symmetric and positive definite so that the conjugate gradient method can solve the system quickly. The proposed immersed boundary method is incorporated into the fluid solver with a second‐order accurate projection method as a plug‐in. The overall scheme is handled under an efficient fractional step framework, namely, prediction, forcing, and projection. Various simulations are performed to validate current method, and the results compare well with previous experimental and numerical studies.     

广义Hamilton(控制)系统的离散梯度积分法

对于广义Hamilton系统及广义Hamilton控制系统，基于能量的Hamilton函数，用离散梯度方法给出了系统保持Hamilton函数特征的数值解法，证明了积分方法可有效地保持Hamilton函数随时间的变化率。通过算例说明了本文方法的有效性。     

Optimal control applied to water flow using second order adjoint method

This paper presents an optimal control applied to water flow using the first and second order adjoint equations. The gradient of the performance function with respect to control variables is analytically obtained by the first order adjoint equation. It is not necessary to compute the Hessian matrix directly using the second order adjoint equation. Two numerical studies have been performed to show the adaptability of the present method. The performance of the second order adjoint method is compared with that of the weighted gradient method, Broyden–Fletcher–Goldfarb–Shanno method and Lanczos method. The precise forms of the adjoint equations and the gradient to use for the minimisation algorithm are derived. The computation by the Lanczos method is shown as superior to those of the other methods discussed in this paper. The message passing interface library is used for the communication of parallel computing.     

Approximation of generalized Stokes problems using dual‐mixed finite elements without enrichment

In this work a finite element method for a dual‐mixed approximation of generalized Stokes problems in two or three space dimensions is studied. A variational formulation of the generalized Stokes problems is accomplished through the introduction of the pseudostress and the trace‐free velocity gradient as unknowns, yielding a twofold saddle point problem. The method avoids the explicit computation of the pressure, which can be recovered through a simple post‐processing technique. Compared with an existing approach for the same problem, the method presented here reduces the global number of degrees of freedom by up to one‐third in two space dimensions. The method presented here also represents a connection between existing dual‐mixed and pseudostress methods for Stokes problems. Existence, uniqueness, and error results for the generalized Stokes problems are given, and numerical experiments that illustrate the theoretical results are presented. Copyright © 2010 John Wiley & Sons, Ltd.     

A simple feedback control system: Bifurcations of periodic orbits and chaos

We consider a pendulum subjected to linear feedback control with periodic desired motions. The pendulum is assumed to be driven by a servo-motor with small time constant, so that the feedback control system can be approximated by a periodically forced oscillator. It was previously shown by Melnikov's method that transverse homoclinic and heteroclinic orbits exist and chaos may occur in certain parameter regions. Here we study local bifurcations of harmonics and subharmonics using the second-order averaging method and Melnikov's method. The Melnikov analysis was performed by numerically computing the Melnikov functions. Numerical simulations and experimental measurements are also given and are compared with the previous and present theoretical predictions. Sustained chaotic motions which result from homoclinic and heteroclinic tangles for not only single but also multiple hyperbolic periodic orbits are observed. Fairly good agreement is found between numerical simulation and experimental results.     

Shape Recovery of Viscoelastic Beams After Stowage

The deployment of viscoelastic structures that have been held stowed for a given time duration can be formulated as a viscoelastic boundary value problem in which the prescribed condition switches from constant displacement to constant traction. This paper presents closed-form expressions for the load relaxation and shape recovery of a linear viscoelastic beam subject to such time-varying constraints. It is shown that a viscoelastic beam recovers to its original shape asymptotically over time. The analytical solutions are employed to investigate the effect of temperature and stowage time on the time required to achieve recovery with a specified precision. Based on the time-temperature equivalence principle, the relationship between recovery time and holding duration is concisely presented on a single plot. It is found that recovery time increases with holding duration but with a diminishing effect.     

范式理论与平均法的等价性分析

分析了范式理论与平均法的等价性。得到的结论是：对含有一对纯虚根的二维非线性系统，使用两种方法得到的结果是等价的，并提供了两个算例来证实其结论的正确性。虽然分析是针对一类二维非线性系统，但其结论同样适合于高维非线性系统。     

Comparison of variants of the bi-conjugate gradient method for compressible Navier-Stokes solver with second-moment closure

Variants of the bi-conjugate gradient (Bi-CG) method are used to resolve the problem of slow convergence in CFD when it is applied to complex flow field simulation using higher-order turbulence models. In this study the Navier-Stokes and Reynolds stress transport equations are discretized with an implicit, total variation diminishing (TVD), finite volume formulation. The preconditioning technique of incomplete lower-upper (ILU) factorization is incorporated into the conjugate gradient square (CGS), bi-conjugate gradient stable (Bi-CGSTAB) and transpose-free quasi-minimal residual (TFQMR) algorithms to accelerate convergence of the overall itertive methods. Computations have been carried out for separated flow fields over transonic bumps, supersonic bases and supersonic compression corners. By comparisons of the convergence rate with each other and with the conventional approximate factorization (AF) method it is shown that the Bi-CGSTAB method gives the most efficient convergence rate among these methods and can speed up the CPU time by a factor of 2·4–6·5 as compared with the AF method. Moreover, the AF method may yield somewhat different results from variants of the Bi-CG method owing to the factorization error which introduces a higher level of convergence criterion.     

Dispersion error reduction for acoustic problems using the smoothed finite element method (SFEM)

The smoothed finite element method (SFEM), which was recently introduced for solving the mechanics and acoustic problems, uses the gradient smoothing technique to operate over the cell‐based smoothing domains. On the basis of the previous work, this paper reports a detailed analysis on the numerical dispersion error in solving two‐dimensional acoustic problems governed by the Helmholtz equation using the SFEM, in comparison with the standard finite element method. Owing to the proper softening effects provided naturally by the cell‐based gradient smoothing operations, the SFEM model behaves much softer than the standard finite element method model. Therefore, the SFEM can significantly reduce the dispersion error in the numerical solution. Results of both theoretical and numerical experiments will support these important findings. It is shown clearly that the SFEM suits ideally well for solving acoustic problems, because of the crucial effectiveness in reducing the dispersion error. Copyright © 2015 John Wiley & Sons, Ltd.     

A stabilized finite element method for the incompressible Navier–Stokes equations using a hierarchical basis

Stabilized finite element methods have been shown to yield robust, accurate numerical solutions to both the compressible and incompressible Navier–Stokes equations for laminar and turbulent flows. The present work focuses on the application of higher‐order, hierarchical basis functions to the incompressible Navier–Stokes equations using a stabilized finite element method. It is shown on a variety of problems that the most cost‐effective simulations (in terms of CPU time, memory, and disk storage) can be obtained using higher‐order basis functions when compared with the traditional linear basis. In addition, algorithms will be presented for the efficient implementation of these methods within the traditional finite element data structures. Copyright © 2001 John Wiley & Sons, Ltd.     

SIZE OPTIMIZATION FOR CANTILEVER I-BEAM BASED ON VARIATIONAL METHOD

以变分方法为工具,研究了分布载荷下工字型悬臂梁的尺寸优化问题,给出了在强度和位移约束下梁缘条沿展向的宽度函数,并简要介绍了优化准则法和多岛遗传算法. 以平直机翼为工程实例,将之简化为一个工字型悬臂梁,分别采用该3 种方法进行优化设计. 结果表明变分法不仅效率高,而且给出了尺寸分布的函数,这对飞行器概念设计阶段的机翼结构设计具有较大的指导意义.     

Optical measurement and scaling of blasts from gram-range explosive charges

Laboratory-scale experiments with gram-range explosive charges are presented. Optical shadowgraphy and high-speed digital imaging are used to measure the explosive-driven shock-wave position as a function of time. From this, shock Mach number-versus-distance from the explosion center can be found. These data then yield the peak overpressure and duration, which are the key parameters in determining the potential damage from an explosion as well as the TNT equivalent of the explosive. Piezoelectric pressure gage measurements of overpressure duration at various distances from the explosive charges compare well with theoretical calculations. A scaling analysis yields an approach to relate the gram-range blast to a large-scale blast from the same or different explosives. This approach is particularly suited to determining the properties and behavior of exotic explosives like triacetone triperoxide (TATP). Results agree with previous observations that the concept of a single TNT equivalence value is inadequate to fully describe an explosive yield, rather TNT equivalence factor and overpressure duration should be presented as functions of radius.      

A conservative local interface sharpening scheme for the constrained interpolation profile method

A conservative local interface sharpening scheme has been developed for the constrained interpolation profile method with the conservative semi‐Lagrangian scheme, because the conservative semi‐Lagrangian scheme does not feature a mechanism to control the interface thickness, thus causing an increase of numerical error with the advance of the time step. The proposed sharpening scheme is based on the conservative level set method proposed by Olsson and Kreiss. However, because their method can cause excessive deformation of the free‐surface in certain circumstances, we propose an improvement of the method by developing a local sharpening technique. Several advection tests are presented to assess the correctness of the advection and the improved interface sharpening scheme. This is followed by the validations of dam‐breaking flow and the rising bubble flows. The mass of the fluid is exactly conserved and the computed terminal velocity of the rising bubble agrees well with the experiments compared with other numerical methods such as the volume of fluid method (VOF), the front tracking method, and the level set method. Copyright © 2011 John Wiley & Sons, Ltd.     

High-efficiency improved symmetric successive over-relaxation preconditioned conjugate gradient method for solving large-scale finite element linear equations

Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering （CAE） and computer aided manufacturing （CAM）. This paper presents a high-efficiency improved symmetric successive over-relaxation （ISSOR） preconditioned conjugate gradient （PCG） method, which maintains lelism consistent with the original form. Ideally, the by 50% as compared with the original algorithm. the convergence and inherent paralcomputation can It is suitable for be reduced nearly high-performance computing with its inherent basic high-efficiency operations. By comparing with the numerical results, it is shown that the proposed method has the best performance.     

基于Wilson-θ和Newmark-β法的非线性动力学方程改进算法

Wilson-θ法和Newmark-β法是非线性动力学方程求解的常用方法。它们的一个基本步骤是,将方程改写为增量平衡的形式,在每一个积分步长内用状态参量修正平衡方程的系数矩阵,其本质是在单个步长内对系统的非线性环节进行了线性化处理。本文基于增量思想分别改进了Wilson-θ法和Newmark-β法,根据即时解给出下一步的猜测解,然后对猜测解进行迭代校正,最终得到收敛的近似解。算例表明,改进算法的精度更高,且收敛准则简单。更为重要的是,本文方法无须对非线性项进行线性化处理,因而计算效率更高,适应范围更广。     

时滞反馈控制对弹簧摆动力学行为的影响

利用解析和数值方法,以弹簧摆为对象讨论了线性的时滞位移反馈控制对一类平方非线性系统动力学行为的影响。根据多尺度法得到了1:2内共振情况下一次近似解的慢变方程,基于此讨论了反馈控制参数对零解的稳定性和周期解振幅的影响。结果表明:耦合的反馈项在平均方程中并不出现。根据罗斯-霍尔维茨判据发现,没有反馈控制时该系统的零解总是不稳定的,而通过调整反馈增益或反馈时滞就可以很容易地使零解稳定。反馈时滞对周期解振幅的影响呈现周期性,反馈增益或时滞发生变化时,周期解振幅的变化会表现出鞍结分岔现象;同时基于MATLAB软件的数值计算结果验证了该理论分析的正确性。