Contingent derivatives and regularization for noncoercive inverse problems

Abstract

We study the inverse problem of parameter identification in noncoercive variational problems that commonly appear in applied models. We examine the differentiability of the set-valued parameter-to-solution map using the first-order and the second-order contingent derivatives. We explore the inverse problem using the output least-squares and the modified output least-squares objectives. By regularizing the noncoercive variational problem, we obtain a single-valued regularized parameter-to-solution map and investigate its smoothness and boundedness. We also consider optimization problems using the output least-squares and the modified output least-squares objectives for the regularized variational problem. We give a complete convergence analysis showing that for the output least-squares and the modified output least-squares, the regularized minimization problems approximate the original optimization problems suitably. We also provide the first-order and the second-order adjoint method for the computation of the first-order and the second-order derivatives of the output least-squares objective. We provide discrete formulas for the gradient and the Hessian calculation and present numerical results.     

An improved depth-averaged nonhydrostatic shallow water model with quadratic pressure approximation

Phase-resolved information is necessary for many coastal wave problems, for example, for the wave conditions in the vicinity of harbor structures. Two-dimensional (2D) depth-averaging shallow water models are commonly used to obtain a phase-resolved solution near the coast. These models are in general more computationally effective compared with computational fluid dynamics software and will be even more capable if equipped with a parallelized code. In the current article, a 2D wave model solving the depth-averaged continuity equation and the Euler equations is implemented in the open-source hydrodynamic code REEF3D. The model is based on a nonhydrostatic extension and a quadratic vertical pressure profile assumption, which provides a better approximation of the frequency dispersion. It is the first model of its kind to employ high-order discretization schemes and to be fully parallelized following the domain decomposition strategy. Wave generation and absorption are achieved with a relaxation method. The simulations of nonlinear long wave propagations and transformations over nonconstant bathymetries are presented. The results are compared with benchmark wave propagation cases. A large-scale wave propagation simulation over realistic irregular topography is shown to demonstrate the model's capability of solving operational large-scale problems.     

基于正弦控制因子的Lateral变异鲨鱼优化算法

针对传统鲨鱼优化算法在求解高维目标函数时,易早熟收敛,陷入局部最优的缺陷.提出一种基于正弦控制因子的Lateral变异鲨鱼优化算法.通过正弦曲线的特性和自适应惯性权重,改善了传统鲨鱼优化算法中由于随机选取控制因子数值大小可能导致算法在迭代后期全局搜索能力降低的问题,提高了算法在迭代后期的全局收敛能力,并对最佳鲨鱼位置引入Lateral变异策略,加强了算法跳出局部最优的可能性.改进后的算法对多个shifted单峰,多峰以及固定维测试函数进行求解,实验结果表明,对比多种不同优化算法而言,本文所提LSSO算法具有更高的收敛精度和搜索速度.     

Finite deformation swelling of a temperature-sensitive hydrogel cylinder under combined extension-torsion

The swelling behavior of a temperature-sensitive poly-N-isopropylacrylamide(PNIPAM) hydrogel circular cylinder is studied subjected to combined extension-torsion and varied temperature. In this regard, a semi-analytical solution is proposed for general combined loading. A finite element(FE) analysis is conducted, subjecting a hydrogel cylinder to the combined extension-torsion and the varied temperature to evaluate the validity and accuracy of the solution. A user-defined UHYPER subroutine is developed and verified under free and constrained swelling conditions. The FE results illustrate excellent agreement with the semi-analytical solution. Due to the complexity of the problem, some compositions and applied loading factors are analyzed. It is revealed that for larger cross-linked density and larger ending temperature, the cylinder yields higher stresses and smaller radial swelling deformation. Besides, the radial and hoop stresses increase by applying larger twist and axial stretch. The hoop stresses intersect at approximately R/Rout = 0.58, where the hoop stress vanishes. Besides, the axial force has direct and inverse relationships with the axial stretch and the twist, respectively. However, the resultant torsional moment behavior is complex, and the position of the maximum point varies significantly by altering the axial stretch and the twist.     

Limit theorems for a class of critical superprocesses with stable branching

We consider a critical superprocess $\{X;{\mathbf{P}}_{\mu }\}$ with general spatial motion and spatially dependent stable branching mechanism with lowest stable index ${\gamma }_0>1$. We first show that, under some conditions, ${\mathbf{P}}_{\mu }(\left|X_t\right|\ne 0)$ converges to 0 as $t\to \infty$ and is regularly varying with index ${({\gamma }_0-1)}^{-1}$. Then we show that, for a large class of non-negative testing functions $f$, the distribution of $\{X_t\left(f\right);{\mathbf{P}}_{\mu }\left(\cdot \right|6X_{t}6\ne 0)\}$, after appropriate rescaling, converges weakly to a positive random variable ${\mathbf{z}}^{({\gamma }_0-1)}$ with Laplace transform $E\left[e^{-u{\mathbf{z}}^{({\gamma }_0-1)}}\right]=1-{(1+u^{-({\gamma }_0-1)})}^{-1∕({\gamma }_0-1)}.$     

关于不定方程x3-1=749y2

利用同余式、平方剩余、Pell方程的解的性质、递归序列证明了:不定方程x³-1=749y²仅有整数解(x,y)=(1,0).