Optimization of reservoirs in series on a river with a nonlinear storage curve for long-term regulation

We present in this paper a new approach to finding the monthly optimal operation of a multireservoir power system connected in series on a river. The hydroelectric power generation is a highly nonlinear function of the storage, and the conversion factor assigned to each power plant is also a nonlinear function of the storage. We use for both a quadratic function of the storage; the resulting problem has a highly nonlinear objective function and linear constraints. We propose a transformation such that the system equations are reduced to linear-quadratic form. Lagrange and Kuhn-Tucker multipliers are used to adjoin the equality and inequality constraints to the objective function. Numerical results are presented for a real system in operation consisting of two reservoirs in series on a river for widely different water conditions.This work was supported by the National Research Council of Canada, Grant No. A4146. The authors would like to acknowledge data obtained from B.C. Hydro.     

Mathematical modeling of a pressurizer in a pressurized water reactor for control design

A pressurizer is a key equipment to ensure the safe operation of a pressurized water reactor by maintaining the reactor coolant system pressure within allowed tolerances. In this paper, a nonequilibrium three-region pressurizer model for the pressurized water reactor was developed first base on the basic conservation laws of mass and energy for the steam and water in the pressurizer. Then the nonlinear model was linearized to introduce the transfer function models of the pressurizer during insurge and outsurge transients for the controller design of a small pressurized water reactor pressurizer. Based on the developed transfer function models, the closed-loop pressurizer pressure and water level control systems including conventional proportional-integral-derivative controllers were designed employing the control strategy that uses a spray valve and two electric heaters for pressure control and regulates the charging flow rate with the letdown flow rate keeping constant for water level control. Finally, a simulation platform for the small pressurized water reactor pressurizer was developed in MATLAB/Simulink with implementation of the proposed nonlinear and nonequilibrium three-region pressurizer model, a widely-used nonlinear and nonequilibrium three-region pressurizer model and the designed pressure and water level controllers. Two typical load change transients, including the 100% to 90% of full power step load decrease and load rejection from 100% to 25% of full power, were simulated based on the platform. Dynamic responses of the pressure and water level obtained using the two nonlinear pressurizer models were analyzed, and the controller performances were assessed. The analysis and assessment results have shown satisfactory control performance and good robustness of the designed controllers regardless of the pressurizer simulation model adopted, demonstrating the feasibility, effectiveness and accuracy of the developed nonlinear and transfer function models of the pressurizer for dynamic simulation and controller design.     

Finite element approximations to one‐phase nonlinear free boundary problem in groundwater contamination flow

In this article, we consider a single‐phase coupled nonlinear Stefan problem of the water‐head and concentration equations with nonlinear source and permeance terms and a Dirichlet boundary condition depending on the free‐boundary function. The problem is very important in subsurface contaminant transport and remediation, seawater intrusion and control, and many other applications. While a Landau type transformation is introduced to immobilize the free boundary, a transformation for the water‐head and concentration functions is defined to deal with the nonhomogeneous Dirichlet boundary condition, which depends on the free boundary function. An H¹‐finite element method for the problem is then proposed and analyzed. The existence of the approximation solution is established, and error estimates are obtained for both the semi‐discrete schemes and the fully discrete schemes. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Interaction behavior associated with a generalized (2+1)-dimensional Hirota bilinear equation for nonlinear waves

In this paper, we focus on the interaction behavior associated with a generalized (2+1)-dimensional Hirota bilinear equation. With symbolic computation, two types of interaction solutions including lump-kink and lump-soliton ones are derived through mixing two positive quadratic functions with an exponential function, or two positive quadratic functions with a hyperbolic cosine function in the bilinear equation. The completely non-elastic interaction between a lump and a stripe is presented, which shows the lump is drowned or shallowed by the stripe. The interaction between lump and soliton is also given, where the lump moves from one branch to the other branch of the soliton. These phenomena exhibit the dynamics of nonlinear waves and the solutions are useful for the study on interaction behavior of nonlinear waves in shallow water, plasma, nonlinear optics and Bose–Einstein condensates.     

An Analytical Model of Periodic Waves in Shallow Water

An explicit, analytical model is presented of finite-amplitude waves in shallow water. The waves in question have two independent spatial periods, in two independent horizontal directions. Both short-crested and long-crested waves are available from the model. Every wave pattern is an exact solution of the Kadomtsev-Petviashvili equation, and is based on a Riemann theta function of genus 2. These biperiodic waves are direct generalizations of the well-known (simply periodic) cnoidal waves. Just as cnoidal waves are often used as one-dimensional models of “typical” nonlinear, periodic waves in shallow water, these biperiodic waves may be considered to represent “typical” nonlinear, periodic waves in shallow water without the assumption of one-dimensionality.     

波浪上方飞行二维地效翼的非线性分析

针对二维波浪上方飞行的非定常二维地效翼进行了非线性分析．通过对二维奇点在规则波上方运动的Green函数的推导,利用离散涡方法解决了二维波浪上方飞行的非定常地效翼的升力问题．针对不同的几何参数和波浪参数对升力系数进行了研究．通过与定常情况的对比,验证了方法的有效性.     

两相流体非线性渗流模型及其应用

基于三参数非线性渗流运动定律、质量守恒定律及椭圆渗流的概念,建立了低渗透介质中两相流体椭圆非线性渗流数学模型,运用有限差分法与外推法求得了其解,导出了两相流体椭圆非线性渗流条件下油井见水前后开发指标的计算公式,进行了实例分析。结果表明:非线性渗流对含水饱和度分布影响较大;非线性渗流使得水驱油推进速度比线性渗流的快,使油井见水时间提前,使得石油开发指标变差;非线性渗流使得同一时刻的压差比线性渗流的大,使石油开发难度加大。这为低渗油藏垂直裂缝井开发工程提供了科学依据。     

On the interaction of deep water waves and exponential shear currents

A train of periodic deep-water waves propagating on a steady shear current with a vertical distribution of vorticity is investigated by an analytic method, namely the homotopy analysis method (HAM). The magnitude of the vorticity varies exponentially with the magnitude of the stream function, while remaining constant on a particular streamline. The so-called Dubreil–Jacotin transformation is used to transfer the original exponentially nonlinear boundary-value problem in an unknown domain into an algebraically nonlinear boundary-value problem in a known domain. Convergent series solutions are obtained not only for small amplitude water waves on a weak current but also for large amplitude waves on a strong current. The nonlinear wave-current interaction is studied in detail. It is found that an aiding shear current tends to enlarge the wave phase speed, sharpen the wave crest, but shorten the maximum wave height, while an opposing shear current has the opposite effect. Besides, the amplitude of waves and fluid velocity decay over the depth more quickly on an aiding shear current but more slowly on an opposing shear current than that of waves on still water. Furthermore, it is found that Stokes criteria of wave breaking is still valid for waves on a shear current: a train of propagating waves on a shear current breaks as the fiuid velocity at crest equals the wave phase speed. Especially, it is found that the highest waves on an opposing shear current are even higher and steeper than that of waves on still water. Mathematically, this analytic method is rather general in principle and can be employed to solve many types of nonlinear partial differential equations with variable coefficients in science, finance and engineering.     

Stochastic dual dynamic programming applied to nonconvex hydrothermal models

In this paper we apply stochastic dual dynamic programming decomposition to a nonconvex multistage stochastic hydrothermal model where the nonlinear water head effects on production and the nonlinear dependence between the reservoir head and the reservoir volume are modeled. The nonconvex constraints that represent the production function of a hydro plant are approximated by McCormick envelopes. These constraints are split into smaller regions and the McCormick envelopes are used for each region. We use binary variables for this disjunctive programming approach and solve the problem with a decomposition method. We resort to a variant of the L-shaped method for solving the MIP subproblem with binary variables at any stage inside the stochastic dual dynamic programming algorithm. A realistic large-scale case study is presented.     

On the interaction of deep water waves and exponential shear currents

A train of periodic deep-water waves propagating on a steady shear current with a vertical distribution of vorticity is investigated by an analytic method, namely the homotopy analysis method (HAM). The magnitude of the vorticity varies exponentially with the magnitude of the stream function, while remaining constant on a particular streamline. The so-called Dubreil–Jacotin transformation is used to transfer the original exponentially nonlinear boundary-value problem in an unknown domain into an algebraically nonlinear boundary-value problem in a known domain. Convergent series solutions are obtained not only for small amplitude water waves on a weak current but also for large amplitude waves on a strong current. The nonlinear wave-current interaction is studied in detail. It is found that an aiding shear current tends to enlarge the wave phase speed, sharpen the wave crest, but shorten the maximum wave height, while an opposing shear current has the opposite effect. Besides, the amplitude of waves and fluid velocity decay over the depth more quickly on an aiding shear current but more slowly on an opposing shear current than that of waves on still water. Furthermore, it is found that Stokes criteria of wave breaking is still valid for waves on a shear current: a train of propagating waves on a shear current breaks as the fiuid velocity at crest equals the wave phase speed. Especially, it is found that the highest waves on an opposing shear current are even higher and steeper than that of waves on still water. Mathematically, this analytic method is rather general in principle and can be employed to solve many types of nonlinear partial differential equations with variable coefficients in science, finance and engineering.     

Nonlinear Separation Approach to Constrained Extremum Problems

In this paper, by virtue of a nonlinear scalarization function, two nonlinear weak separation functions, a nonlinear regular weak separation function, and a nonlinear strong separation function are first introduced, respectively. Then, by the image space analysis, a global saddle-point condition for a nonlinear function is investigated. It is shown that the existence of a saddle point is equivalent to a nonlinear separation of two suitable subsets of the image space. Finally, some necessary and sufficient optimality conditions are obtained for constrained extremum problems.     

油田注水系统拓扑布局优化的混合遗传算法

以投资最小为目标函数,建立了注水系统拓扑布局优化数学模型.根据模型特点,将优化问题分为两层,分别采用遗传算法和非线性优化方法进行求解.并对遗传算法的操作过程进行了改进,调整了适应函数,改进了交叉和变异操作,结合了模拟退火算法,在操作过程中使约束条件得到满足,减少了不可行解的产生,使遗传算法的优化性能得到了提高.优化算例说明了该方法的有效性.     

Nonlinear separation approach to inverse variational inequalities

In this paper, we employ the image space analysis to investigate an inverse variational inequality (for short, IVI) with a cone constraint. By virtue of the nonlinear scalarization function commonly known as the Gerstewitz function, three nonlinear weak separation functions, two nonlinear regular weak separation functions and a nonlinear strong separation function are first introduced. Then, by these nonlinear separation functions, theorems of the weak and strong alternative and some optimality conditions for IVI with a cone constraint are derived without any convexity. In particular, a global saddle-point condition for a nonlinear function is investigated. It is shown that the existence of a saddle point is equivalent to a nonlinear separation of two suitable subsets of the image space. Finally, two gap functions and an error bound for IVI with a cone constraint are obtained.     

基于学习向量量化神经网络的水环境质量综合评价

针对水环境质量综合评价中指标权重确定方法的不足,利用学习向量量化(LVQ)神经网络具有的强大的非线性运算和相似特征聚类功能,提出一种基于学习向量量化(LVQ)神经网络的水质综合评价决策方法.将它应用于水质综合指标评价,为改进水质综合评价提供了一种简捷的分类评价方法.     

Periodic wave solutions and asymptotic analysis of the Hirota–Satsuma shallow water wave equation

In this paper, we devise a simple way to explicitly construct the Riemann theta function periodic wave solution of the nonlinear partial differential equation. The resulting theory is applied to the Hirota–Satsuma shallow water wave equation. Bilinear forms are presented to explicitly construct periodic wave solutions based on a multidimensional Riemann theta function. We obtain the one‐periodic and two‐periodic wave solutions of the equation. The relations between the periodic wave solutions and soliton solutions are rigorously established. Copyright © 2014 John Wiley & Sons, Ltd.     

Nonlinear low frequency water waves in a cylindrical shell subjected to high frequency excitations – Part I: Experimental study

This paper presents an experimental investigation on nonlinear low frequency gravity water waves in a partially filled cylindrical shell subjected to high frequency horizontal excitations. The characteristics of natural frequencies and mode shapes of the water–shell coupled system are discussed. The boundaries for onset of gravity waves are measured and plotted by curves of critical excitation force magnitude with respect to excitation frequency. For nonlinear water waves, the time history signals and their spectrums of motion on both water surface and shell are recorded. The shapes of water surface are also measured using scanning laser vibrometer. In particular, the phenomenon of transitions between different gravity wave patterns is observed and expressed by the waterfall graphs. These results exhibit pronounced nonlinear properties of shell–fluid coupled system.     

用非线性方程组求解等式约束非线性规划问题的降维算法

本文研究线性和非线性等式约束非线性规划问题的降维算法.首先,利用一般等式约束问题的降维方法,将线性等式约束非线性规划问题转换成一个非线性方程组,解非线性方程组即得其解;然后,对线性和非线性等式约束非线性规划问题用Lagrange乘子法,将非线性约束部分和目标函数构成增广的Lagrange函数,并保留线性等式约束,这样便得到一个线性等式约束非线性规划序列,从而,又将问题转化为求解只含线性等式约束的非线性规划问题.     

The nonlinear bilevel programming problem:formulations,regularity and optimality conditions

In this paper, we study regularity and optimality conditions for the BLPP by using a marginal function formulation, where the marginal function is defined by the optimal value function of the lower problem. We address the regularity issue by exploring the structure of the tangent cones of the feasible set of the BLPP. These regularity results indicate that the nonlinear/nonlinear BLPP is most likely degenerate and a class of nonlinear/linear BLPP is regular in the conventional sense. Existence of exact penalty function is proved for a class of nonlinear/linear BLPP. Fritz-John type optimality conditions are derived for nonlinear BLPP, while KKT type conditions are obtained for a class of nonlinear/linear BLPP in the framework of nonsmooth analysis. A typical example is examined for these conditions and some applications of these conditions are pointed out     

Relaxation approximation to bed-load sediment transport

In this work we propose and apply a numerical method based on finite volume relaxation approximation for computing the bed-load sediment transport in shallow water flows, in one and two space dimensions. The water flow is modeled by the well-known nonlinear shallow water equations which are coupled with a bed updating equation. Using a relaxation approximation, the nonlinear set of equations (and for two different formulations) is transformed to a semilinear diagonalizable problem with linear characteristic variables. A second order MUSCL-TVD method is used for the advection stage while an implicit–explicit Runge–Kutta scheme solves the relaxation stage. The main advantages of this approach are that neither Riemann problem solvers nor nonlinear iterations are required during the solution process. For the two different formulations, the applicability and effectiveness of the presented scheme is verified by comparing numerical results obtained for several benchmark test problems.     

Robust Stability and Stabilization of a Class of Nonlinear Itô-Type Stochastic Systems via Linear Matrix Inequalities

This article presents a new approach to robust quadratic stabilization of nonlinear stochastic systems. The linear rate vector of a stochastic system is perturbed by a nonlinear function, and this nonlinear function satisfies a quadratic constraint. Our objective is to show how linear constant feedback laws can be formulated to stabilize this type of stochastic systems and, at the same time maximize the bounds on this nonlinear perturbing function which the system can tolerate without becoming unstable. The new formulation provides a suitable setting for robust stabilization of nonlinear stochastic systems where the underlying deterministic systems satisfy the generalized matching conditions. Our sufficient conditions are written in matrix forms, which are determined by solving linear matrix inequalities (LMIs), which have significant computational advantage over any other existing techniques. Examples are given to demonstrate the results.