The goal of this paper is to study a mixed finite element approximation of the general convex optimal control problems governed by quasilinear elliptic partial differential equations. The state and co-state are approx...
详细信息
The goal of this paper is to study a mixed finite element approximation of the general convex optimal control problems governed by quasilinear elliptic partial differential equations. The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise constant functions. We derive a priori error estimates both for the state variables and the control variable. Finally, some numerical examples are given to demonstrate the theoretical results.
This paper investigates an optimal control problem governed by an elliptic equation with integral control and state *** control problem is approxi-mated by the hp spectral element method with high accuracy and geometr...
详细信息
This paper investigates an optimal control problem governed by an elliptic equation with integral control and state *** control problem is approxi-mated by the hp spectral element method with high accuracy and geometricfl*** conditions of the continuous and discrete optimal control problems are presented,*** a posteriori error estimates both for the control and state variables are established in *** addition,illustrative numerical examples are carried out to demonstrate the accuracy of theoretical results and the validity of the proposed method.
Based on W -transformation, some parametric symplectic partitioned Runge–Kutta (PRK) methods depending on a real parameter α are developed. For α = 0 , the corresponding methods become the usual PRK methods, includ...
Based on W -transformation, some parametric symplectic partitioned Runge–Kutta (PRK) methods depending on a real parameter α are developed. For α = 0 , the corresponding methods become the usual PRK methods, including Radau IA – I A ¯ and Lobatto IIIA – IIIB methods as examples. For any α ≠ 0 , the corresponding methods are symplectic and there exists a value α ∗ such that energy is preserved in the numerical solution at each step. The existence of the parameter and the order of the numerical methods are discussed. Some numerical examples are presented to illustrate these results.
In this paper, we present a novel indirect convergent Jacobi spectral collocation method for fractional optimal control problems governed by a dynamical system including both classical derivative and Caputo fractional...
Numerical simulation of the stress-strain state of a composite material may be difficult due to large computational complexity associated with a grid resolution of a large number of inclusions. To overcome the problem...
详细信息
In this paper, the fractional variational integrators developed by Wang and Xiao (2012) [28] are extended to the fractional Euler–Lagrange (E–L) equations with holonomic constraints. The corresponding fractional dis...
详细信息
In this paper, the fractional variational integrators developed by Wang and Xiao (2012) [28] are extended to the fractional Euler–Lagrange (E–L) equations with holonomic constraints. The corresponding fractional discrete E–L equations are derived, and their local convergence is discussed. Some fractional variational integrators are presented. The suggested methods are shown to be efficient by some numerical examples.
The Gross–Pitaevskii equation is the model equation of the single-particle wave function in a Bose–Einstein condensation. A computation difficulty of the Gross–Pitaevskii equation comes from the semiclassical probl...
详细信息
The Gross–Pitaevskii equation is the model equation of the single-particle wave function in a Bose–Einstein condensation. A computation difficulty of the Gross–Pitaevskii equation comes from the semiclassical problem in supercritical case. In this paper, we apply a diffeomorphism to transform the original one-dimensional Gross–Pitaevskii equation into a modified equation. The adaptive grids are constructed through the interpolating wavelet method. Then, we use the time-splitting finite difference method with the wavelet-adaptive grids to solve the modified Gross–Pitaevskii equation, where the approximation to the second-order derivative is given by the Lagrange interpolation method. At last, the numerical results are given. It is shown that the obtained time-splitting finite difference method with the wavelet-adaptive grids is very efficient for solving the one-dimensional semiclassical Gross–Pitaevskii equation in supercritical case and it is suitable to deal with the local high oscillation of the solution.
Particulate flows are commonly encountered in both engineering and environmental applications. The discrete element method (DEM) has attracted plentiful attentions since it can predict the whole motion of the particul...
详细信息
It is important to calculate the reachable domain(RD)of the manned lunar mission to evaluate whether a lunar landing site could be reached by the spacecraft. In this paper, the RD of free return orbits is quickly eval...
详细信息
It is important to calculate the reachable domain(RD)of the manned lunar mission to evaluate whether a lunar landing site could be reached by the spacecraft. In this paper, the RD of free return orbits is quickly evaluated and calculated via the classification and regression neural networks. An efficient databasegeneration method is developed for obtaining eight types of free return orbits and then the RD is defined by the orbit’s inclination and right ascension of ascending node(RAAN) at the perilune. A classify neural network and a regression network are trained respectively. The former is built for classifying the type of the RD, and the latter is built for calculating the inclination and RAAN of the RD. The simulation results show that two neural networks are well trained. The classification model has an accuracy of more than 99% and the mean square error of the regression model is less than 0.01°on the test set. Moreover, a serial strategy is proposed to combine the two surrogate models and a recognition tool is built to evaluate whether a lunar site could be reached. The proposed deep learning method shows the superiority in computation efficiency compared with the traditional double two-body model.
Anisotropic meshes are known to be well-suited for problems which exhibit anisotropic solution features. Defining an appropriate metric tensor and designing an efficient algorithm for anisotropic mesh generation are t...
详细信息
Anisotropic meshes are known to be well-suited for problems which exhibit anisotropic solution features. Defining an appropriate metric tensor and designing an efficient algorithm for anisotropic mesh generation are two important aspects of the anisotropic mesh methodology. In this paper, we are concerned with the natural metric tensor for use in anisotropic mesh generation for anisotropic elliptic problems. We provide an algorithm to generate anisotropic meshes under the given metric tensor. We show that the inverse of the anisotropic diffusion matrix of the anisotropic elliptic problem is a natural metric tensor for the anisotropic mesh generation in three aspects: better discrete algebraic systems, more accurate finite element solution and superconvergence on the mesh nodes. Various numerical examples demonstrating the effectiveness are presented.
暂无评论