基于非均匀参数化的自由终端时间最优控制问题求解

针对自由终端时间最优控制问题,提出了一种基于非均匀控制向量参数化的数值解法.将控制时域离散化为不同长度的时间段,各时间段长度作为新的控制变量.通过引入标准化的时间变量,原问题转化为均匀参数化的固定终端时间最优控制问题.建立目标和约束函数的Hamilton函数,通过求解伴随方程获得目标和约束函数的梯度,采用序列二次规划(SQP)获得数值解.针对两个经典的化工过程自由终端时间最优控制问题进行仿真研究,验证了所提出算法的可行性和有效性.     

Subspace Design of Compactly Supported Orthonormal Wavelets

We derive a new matrix parameterization of compactly supported orthonormal wavelets where the coefficients of the wavelet filter are the solution of a linear system of equations that is parameterized by an arbitrary vector. The parameterization shows that the vector of the wavelet filter coefficients is the kernel of a subspace of the condition matrix row-space. This property is exploited to develop a new design procedure for orthonormal wavelets of compact support. The proposed parameterization also describes the class of two-channel orthogonal filter banks where in this case we have two extra degrees of freedom in the design. The effectiveness of the proposed procedure is illustrated by design examples of common orthonormal wavelets.     

A model of a general elastic curved rod

We indicate a new approach to the deformation of three‐dimensional curved rods with variable cross‐section. The model consists of a system of nine ordinary differential equations for which we prove existence and uniqueness via the coercivity of the association bilinear form. From the geometrical point of view, we are using the Darboux frame or a new local frame requiring just a C¹‐parameterization of the curve. Our model also describes the deformation occurring in the cross‐sections of the rod. Copyright © 2002 John Wiley & Sons, Ltd.     

The Runge‐Kutta DG finite element method and the KFVS scheme for compressible flow simulations

This article presents a new type of second‐order scheme for solving the system of Euler equations, which combines the Runge‐Kutta discontinuous Galerkin (DG) finite element method and the kinetic flux vector splitting (KFVS) scheme. We first discretize the Euler equations in space with the DG method and then the resulting system from the method‐of‐lines approach will be discretized using a Runge‐Kutta method. Finally, a second‐order KFVS method is used to construct the numerical flux. The proposed scheme preserves the main advantages of the DG finite element method including its flexibility in handling irregular solution domains and in parallelization. The efficiency and effectiveness of the proposed method are illustrated by several numerical examples in one‐ and two‐dimensional spaces. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A finite‐volume scheme for a spinorial matrix drift‐diffusion model for semiconductors

An implicit Euler finite‐volume scheme for a spinorial matrix drift‐diffusion model for semiconductors is analyzed. The model consists of strongly coupled parabolic equations for the electron density matrix or, alternatively, of weakly coupled equations for the charge and spin‐vector densities, coupled to the Poisson equation for the electric potential. The equations are solved in a bounded domain with mixed Dirichlet–Neumann boundary conditions. The charge and spin‐vector fluxes are approximated by a Scharfetter–Gummel discretization. The main features of the numerical scheme are the preservation of nonnegativity and bounds of the densities and the dissipation of the discrete free energy. The existence of a bounded discrete solution and the monotonicity of the discrete free energy are proved. For undoped semiconductor materials, the numerical scheme is unconditionally stable. The fundamental ideas are reformulations using spin‐up and spin‐down densities and certain projections of the spin‐vector density, free energy estimates, and a discrete Moser iteration. Furthermore, numerical simulations of a simple ferromagnetic‐layer field‐effect transistor in two space dimensions are presented. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 819–846, 2016     

Comparison of the Solution of the Riemann Problem for Inviscid/Viscous Flows

The aim of this paper is to present the influence of the method used to solve the convective fluxes on the transonic internal flow field. The governing equations are discretized using an upwind method based on different solutions of the Riemann problem, flux vector splitting (FVS) or flux difference splitting schemes (FDS). Turbulence effects are simulated by means of the low‐Reynolds‐number k – ϵ and the SST (Shear‐Stress‐Transport) turbulence models.     

Low frequency tangential filtering decomposition

For block‐tridiagonal systems of linear equations arising from the discretization of partial differential equations, a composite preconditioner is proposed and tested. It combines a classical ILU0 factorization for high frequencies with a tangential filtering preconditioner. The choice of the filtering vector is important: the test‐vector is the Ritz eigenvector corresponding to the approximate lowest eigenvalue, obtained after a limited number of iterations of a ILU0 preconditioned Krylov method. Numerical tests are carried out for this method. Copyright © 2006 John Wiley & Sons, Ltd.     

A numerical solution of problems in calculus of variation using direct method and nonclassical parameterization

In this paper a direct method for solving variational problems using nonclassical parameterization is presented. A nonclassical parameterization based on nonclassical orthogonal polynomials is first introduced to reduce a variational problem to a nonlinear mathematical programming problem. Then, using the Lagrange multiplier technique the problem is converted to that of solving a system of algebraic equations. Illustrative examples are given to demonstrate the validity and applicability of the technique.     

Steger–Warming flux vector splitting method for special relativistic hydrodynamics

This paper discusses the properties of the rotational invariance and hyperbolicity in time of the governing equations of the ideal special relativistic hydrodynamics and proves for the first time that the ideal relativistic hydrodynamical equations satisfy the homogeneity property, which is the footstone of the Steger–Warming flux vector splitting method [J. L. Steger and R. F. Warming, J. Comput. Phys., 40(1981), 263–293]. On the basis of this remarkable property, the Steger–Warming flux vector splitting (SW‐FVS) is given. Two high‐resolution SW‐FVS schemes are also given on the basis of the initial reconstructions of the solutions and the fluxes, respectively. Several numerical experiments are conducted to validate the performance of the SW‐FVS method. Copyright © 2013 John Wiley & Sons, Ltd.     

Some predictor–corrector‐type iterative schemes for solving nonsymmetric algebraic Riccati equations arising in transport theory

It is as well known that nonsymmetric algebraic Riccati equations arising in transport theory can be translated to vector equations. In this paper, we propose six predictor–corrector‐type iterative schemes to solve the vector equations. And we give the convergence of these schemes. Unlike the previous work, we prove that all of them converge to the minimal positive solution of the vector equations by the initial vector (e,e), where e = (1,1, ? ,1)^T. Moreover, we prove that all the sequences generated by the iterative schemes are strictly and monotonically increasing and bounded above. In addition, some numerical results are also reported in the paper, which confirm the good theoretical properties of our approach. Copyright © 2014 John Wiley & Sons, Ltd.     

Approximating Rough Stochastic PDEs

We study approximations to a class of vector‐valued equations of Burgers type driven by a multiplicative space‐time white noise. A solution theory for this class of equations has been developed recently in Probability Theory Related Fields by Hairer and Weber. The key idea was to use the theory of controlled rough paths to give definitions of weak/mild solutions and to set up a Picard iteration argument. In this article the limiting behavior of a rather large class of (spatial) approximations to these equations is studied. These approximations are shown to converge and convergence rates are given, but the limit may depend on the particular choice of approximation. This effect is a spatial analogue to the Itô‐Stratonovich correction in the theory of stochastic ordinary differential equations, where it is well known that different approximation schemes may converge to different solutions.© 2014 Wiley Periodicals, Inc.     

Weak solutions for fractional differential equations in nonreflexive Banach spaces via Riemann‐Pettis integrals

The aim of this paper is to develop fractional calculus for vector‐valued functions using the weak Riemann integral. Also, we establish the existence of weak solutions for a class of fractional differential equations with fractional weak derivatives.     

多体系统动力学微分／代数方程组的一类新的数值分析方法

本文讨论了多体系统动力学微分／代数混合方程组的数值离散问题．首先把参数t并入广义坐标讨论，简化了方程组及其隐含条件的结构，并将其化为指标1的方程组．然后利用方程组的特殊结构，引入一种局部离散技巧并构造了相应的算法．算法结构紧凑，易于编程，具有较高的计算效率和良好的数值性态，且其形式适合于各种数值积分方法的的实施．文末给出了具体算例．     

Continuum Limit of the Triple Tau-Function Model

We present a system of integrable second-order differential equations for three fields in the three-dimensional space–time. The system is obtained as the continuum limit of discrete equations for a triplet of tau-functions. We give a parameterization of the soliton solutions of equations of motion, describe the linear problem, and establish the integrability of the corresponding classical field theory.     

Numerical solutions of distributed order fractional differential equations in the time domain using the Müntz–Legendre wavelets approach

In this paper, a numerical method is presented to obtain and analyze the behavior of numerical solutions of distributed order fractional differential equations of the general form in the time domain with the Caputo fractional derivative. The suggested method is based on the Müntz–Legendre wavelet approximation. We derive a new operational vector for the Riemann–Liouville fractional integral of the Müntz–Legendre wavelets by using the Laplace transform method. Applying this operational vector and collocation method in our approach, the problem can be reduced to a system of linear and nonlinear algebraic equations. The arising system can be solved by the Newton method. Discussion on the error bound and convergence analysis for the proposed method is presented. Finally, seven test problems are considered to compare our results with other well‐known methods used for solving these problems. The results in the tabulated tables highlighted that the proposed method is an efficient mathematical tool for analyzing distributed order fractional differential equations of the general form.     

Intersection of a ruled surface with a free-form surface

This paper presents a simple method for computing the intersection curve of a ruled surface and a free-form surface. The basic idea is to reduce the problem of surface intersection to the one of projecting an appropriate curve such as a directrix of the ruled surface, along its indicatrix curve (direction vector field of its generating lines), onto the free-form surface; the projection curve is just the intersection curve. With techniques in classical differential geometry, we derive the differential equations of the intersection curve in the cases of parametrically and implicitly defined free-form surfaces. The intersection curve naturally inherits the parameter of the chosen directrix. Moreover, it is independent of the base surface geometry and its parameterization, and is obtained by numerically solving the initial-value problem for a system of first-order ordinary differential equations in the parametric domain associated to the surface representation for parametric case or in 3D space for implicit case. Some experimental examples are also given to demonstrate that the presented method is effective and potentially useful in computer aided design and computer graphics. An erratum to this article can be found at     

Simulation of Multibody Systems with Nonholonomic Constraints

The paper deals with the nonholonomic multibody system dynamics from a point of view resulting from some present applications in high‐tec areas like high‐speed train technology or biomechanics of some disciplines in high‐performance sports. A formulation of nonholonomic constraints which are linear related to generalized velocities is based on a derivative‐free approach for generating Lagrangian motion equations of multibody systems with kinematical tree structure as well as for constrained multibody systems. This has been done by using di.erential‐geometric concepts in a Riemannian space. The ideas are illustrated by the classical edge condition on double‐curved surfaces. The surfaces are described by C²‐vector functions, for example by NURBS‐approximation. As an example a bobsleigh is regarded moving on a double‐curved surface.     

Finitely smooth solutions of nonlinear singular partial differential equations

We solve a Fuchsian system of singular nonlinear partial differential equations with resonances. These equations have no smooth solutions in general. We show the solvability in a class of finitely smooth functions. Typical examples are a homology equation for a vector field and a degenerate Monge–Ampère equation. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Two‐grid methods for mixed finite‐element solution of coupled reaction‐diffusion systems

We develop 2‐grid schemes for solving nonlinear reaction‐diffusion systems: where p = (p, q) is an unknown vector‐valued function. The schemes use discretizations based on a mixed finite‐element method. The 2‐grid approach yields iterative procedures for solving the nonlinear discrete equations. The idea is to relegate all the Newton‐like iterations to grids much coarser than the final one, with no loss in order of accuracy. The iterative algorithms examined here extend a method developed earlier for single reaction‐diffusion equations. An application to prepattern formation in mathematical biology illustrates the method's effectiveness. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 589–604, 1999