A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

双掺(Tm3+,Tb3+)LiYF4激光器1.5 μm波长激光阈值分析

由速率方程推出了双掺（Tm^3 ,Tb^3 ）离子准四能级系统的激光阈值解析式，讨论了Tm^3 和Tb^3 离子之间的相互作用。分析了1.5μm波长附近的激光阈值和Tm^3 、Tb^3 离子的掺杂原子数分数及晶体长度的关系。结果表明，对于对应Tm^3 离子^3H4→^3F4跃迁的约1.5μm波长的激光，激活离子Tm^3 的掺杂原子数分数过大时，交叉弛豫作用将使系统阈值迅速增加。Tb^3 离子的加入，一方面能抽空激光下能级，起到降低阈值的作用；另一方面亦减少了激光上能级的寿命，使阈值升高。故Tb^3 离子有最佳掺杂原子数分数。对于Tm原子数分数为y=0.01的Tm:LiYF4晶体，Tb^3 离子的最佳掺杂原子数分数为0.002左右，同时表明，激光阈值与晶体长度有关。最佳晶体长度与Tm^3 、Tb^3 离子的掺杂原子数分数以及晶体的衍射损耗和吸收损耗有关。     

Approximate inverse computation using Frobenius inner product

Parallel preconditioners are presented for the solution of general linear systems of equations. The computation of these preconditioners is achieved by orthogonal projections related to the Frobenius inner product. So, min_M∈??∥AM?I∥ _F and matrix M₀∈?? corresponding to this minimum (?? being any vectorial subspace of ??_n(?)) are explicitly computed using accumulative formulae in order to reduce computational cost when subspace ?? is extended to another one containing it. Every step, the computation is carried out taking advantage of the previous one, what considerably reduces the amount of work. These general results are illustrated with the subspace of matrices M such that AM is symmetric. The main application is developed for the subspace of matrices with a given sparsity pattern which may be constructed iteratively by augmenting the set of non‐zero entries in each column. Finally, the effectiveness of the sparse preconditioners is illustrated with some numerical experiments. Copyright © 2002 John Wiley & Sons, Ltd.     

一阶最优性条件研究

本对由Botsko的关于多变量函数取极值的一阶导数检验条件定理^[1]进行了分析研究，给出了更实用而简捷的差别条件。最后，举出若干例子予以说明。     

Optimal Planar Interception with Fixed End Conditions: Approximate Solutions

In a previous paper (Ref. 1), an exact solution of the optimal planar interception with fixed end conditions was derived in closed form. The optimal control was expressed as an explicit function of the state variables and two fixed parameters, obtained by solving a set of nonlinear algebraic equations involving elliptic integrals. In order to facilitate the optimal control implementation, the present paper derives a highly accurate simplified solution assuming that the ratio of the pursuer turning radius to the initial range is small. An asymptotic expansion further reduces the computational workload. Construction of a near-optimal open-loop control, based on the approximations, completes the present paper.     

Polarized Spacetime Foam

An approximate model of a spacetime foam is presented. It is supposed that in the spacetime foam each quantum handle is like to an electric dipole and therefore the spacetime foam is similar to a dielectric. If we neglect of linear sizes of the quantum handle then it can be described with an operator containing a Grassman number and either a scalar or a spinor field. For both fields the Lagrangian is presented. For the scalar field it is the dilaton gravity + electrodynamics and the dilaton field is a dielectric permeability. The spherically symmetric solution in this case give us the screening of a bare electric charge surrounded by a polarized spacetime foam and the energy of the electric field becomes finite one. In the case of the spinor field the spherically symmetric solution give us a ball of the polarized spacetime foam filled with the confined electric field. It is shown that the full energy of the electric field in the ball can be very big.     

一类机器人系统的最优控制

本文通过把结构阻尼系数当作控制变量来讨论一类弹性机器人系统的最优控制问题 ,并利用Banach空间几何性质证明了最优控制元的存在唯一性     

Sensitivity analysis for discrete optimal control problems

We present local sensitivity analysis for discrete optimal control problems with varying endpoints in the case when the customary regularity of boundary conditions can be violated. We study the behavior of the optimal solutions subject to parametric perturbations of the problem.     

线性分式规划最优解集的求法

本文使用多面集的表示定理，导出了线性分式规划最优解集的结构，并给出确定全部最优解的计算步骤。     

OPTIMAL INVESTMENT IN THE PROTECTION OF A VULNERABLE BIOLOGICAL RESOURCE

It is assumed that the probability of destruction of a biological asset by natural hazards can be reduced through investment in protection. Specifically a model, in which the hazard rate depends on both the age of the asset and the accumulated invested protection capital, is assumed. The protection capital depreciates through time and its effectiveness in reducing the hazard rate is subject to diminishing returns. It is shown how the investment schedule to maximize the expected net present value of the asset can be determined using the methods of deterministic optimal control, with the survival probability regarded as a state variable. The optimal investment pattern involves “bang-bang-singular” control. A numerical scheme for determining jointly the optimal investment policy and the optimal harvest (or replacement) age is outlined and a numerical example involving forest fire protection is given.