Speckle suppression is extremely important for understanding and utilising synthetic aperture radar (SAR) images, while the emphasis of the traditional methods for speckle suppression is usually focused on removing th...
详细信息
Speckle suppression is extremely important for understanding and utilising synthetic aperture radar (SAR) images, while the emphasis of the traditional methods for speckle suppression is usually focused on removing the noise instead of keeping the scattering character of imaging objects, which has caused serious interference to the subsequent applications of SAR images. The authors aim to develop sparse partial differential equation (PDE) for speckle suppression of SAR image, where the PDE model is associated with the sparse prior of objects and the statistical property of speckle. The PDE model has been proved to have the ability of denoising and edge-preserving by proper design, and the sparse prior of the point- and line-like objects on SAR images has been also illustrated, both of which help for keeping the scattering characters. To solve the proposed sparse PDE model, a numerical algorithm is designed and the sparse constraint is realised in each step of the diffusion process. In experiments, several real SAR images are utilised to validate the performance of the proposed method.
This study presents a new analytic approach to optimal coordination of directional overcurrent relays. This approach is based on the selection of optimum pickup current and time dial setting, in order to obtain minimu...
详细信息
This study presents a new analytic approach to optimal coordination of directional overcurrent relays. This approach is based on the selection of optimum pickup current and time dial setting, in order to obtain minimum operating time for the relays, while satisfying various coordination and boundary constraints. Based on the new optimal relay setting procedure, an iterative numerical solution is proposed. The proposed numerical algorithm converges to the global optimum values, which are independent of initial values and the order of relay setting. The proposed method is applied to three different test systems. The new method is compared with some previously proposed analytic and evolutionary approaches. The results demonstrate the advantages of the proposed method over the previous works.
In this paper, a numerical method for computing the relaxation modulus of a linearly viscoelastic material is presented. The method is valid for relaxation tests where a constant strain rate is followed by a constant ...
详细信息
In this paper, a numerical method for computing the relaxation modulus of a linearly viscoelastic material is presented. The method is valid for relaxation tests where a constant strain rate is followed by a constant strain. The method is similar to the procedure suggested by Zapas and Phillips. Unlike Zapas-Phillips approach, this new method can be also applied for times shorter than t (1)/2, where t (1) denotes time when the maximum strain is achieved. Therefore this method is very suitable for materials that experiences fast relaxation. The method is verified with numerical simulations. Results from the simulations are compared with analytical solution and Zapas-Phillips method. Results indicate that the presented approach is suitable for estimating the relaxation modulus.
A novel, simple iterative algorithm is used to calculate the temperature distribution in a finite medium for the case of non-Fourier (hyperbolic) heat conduction. In this algorithm the temperature is calculated explic...
详细信息
A novel, simple iterative algorithm is used to calculate the temperature distribution in a finite medium for the case of non-Fourier (hyperbolic) heat conduction. In this algorithm the temperature is calculated explicitly in one simple calculation that is repeated for each time step as the heat wave propagates through the medium with constant speed. When the wave reaches a boundary of the medium, it bounces back and moves in the opposite direction. All simple initial and boundary conditions can be modelled. An example of using the algorithm for the case of a finite, thermally insulated medium is given, and the results are compared with an exact analytical solution.
Aiming at the shortcomings of the analytical algorithm and numerical algorithm for calculating the eddy current losses in the solid rotor of high speed squirrel cage induction motor, this study proposes a semi-analyti...
详细信息
Aiming at the shortcomings of the analytical algorithm and numerical algorithm for calculating the eddy current losses in the solid rotor of high speed squirrel cage induction motor, this study proposes a semi-analytical algorithm. Based on the magnetic field of the rotor surface and Maxwell's equations, the analytical formula of this algorithm is derived. According to the law of energy conservation, the losses in rotor and rotor bars are calculated, respectively, so that the eddy current losses in the solid rotor are calculated indirectly. Taking a high speed squirrel cage induction motor as an example, the eddy current losses in the solid rotor are calculated by the semi-analytical method. The calculation result is verified by the two-dimensional transient finite-element method, and the calculation result meets the basic requirements of engineering practice.
The mathematical model of the Baltic Sea dynamics developed at the Institute of numerical Mathematics of RAS is considered. The problem of variational assimilation of average daily data for the sea surface temperature...
详细信息
The mathematical model of the Baltic Sea dynamics developed at the Institute of numerical Mathematics of RAS is considered. The problem of variational assimilation of average daily data for the sea surface temperature (SST) is formulated and studied with the use of covariance matrices of observation data errors. Based on variational assimilation of satellite observation data, we propose an algorithm for solving the inverse problem of the heat flux reconstruction on the sea surface. The results of numerical experiments on reconstruction of the heat flux function are presented for the problem of variational assimilation of observation SST data.
Gauss-Newton method is mainly designed and exploited to solve the nonlinear least squares problems of static overdetermined systems with order-2 precision. When the system comes to the future time-variant situation, t...
详细信息
Gauss-Newton method is mainly designed and exploited to solve the nonlinear least squares problems of static overdetermined systems with order-2 precision. When the system comes to the future time-variant situation, the disadvantages of the Gauss-Newton method in effectiveness and precision need to be overcome. On the basis of Zhang neurodynamics, the first Zhang neurodynamics model and the second Zhang neurodynamics model are constructed and discussed to obtain the nonlinear least squares solution of the time-variant overdetermined system in the future time-variant situation. Then, we propose a novel eleven-instant Zhang time discretization formula whose convergenceness and truncation errors of order-6 are proved through rigorous mathematical derivation. The proposed eleven-instant Zhang time discretization formula is thus utilized to discretize the constructed first Zhang neurodynamics model and second Zhang neurodynamics model, respectively. Consequently, a novel discrete-time eleven-instant first Zhang neurodynamics-based algorithm and a novel discrete-time eleven-instant second Zhang neurodynamics-based algorithm are proposed and generalized to obtain the predictive nonlinear least squares solution of the time-variant overdetermined system in the future time-variant situation. Furthermore, the precisions of the proposed discrete-time eleven-instant first Zhang neurodynamics-based and second Zhang neurodynamics-based algorithms are proved to be of order-6 through mathematical derivation. Moreover, this work reveals the facts that the Gauss-Newton method is actually the simplified version of the discrete-time two-instant first Zhang neurodynamics-based algorithm and the discrete-time two-instant second Zhang neurodynamics-based algorithm, respectively. Finally, five challenging numerical experiments, including the problem-solving of the critical force of a crane in a complicated environment, are carried out. The experimental results further substantiate the appli
In this article, our damped-heat wave (DHW) algorithm for the calculation of temperature distribution in a homogeneous finite medium is compared with Fourier's algorithm for transmission of heat between discrete b...
详细信息
In this article, our damped-heat wave (DHW) algorithm for the calculation of temperature distribution in a homogeneous finite medium is compared with Fourier's algorithm for transmission of heat between discrete bodies.
We have extended our previous 2D method [Gerya, T.V., Yuen, D.A., 2003. Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable tra...
详细信息
We have extended our previous 2D method [Gerya, T.V., Yuen, D.A., 2003. Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable transport properties. Phys. Earth Planet. Interiors 140, 295-320], which is a combination of conservative finite-differences with marker-in-cell techniques to include the effects of visco-elasto-plastic rheology, self-gravitation and a self-consistently derived evolving curvilinear planetary surface. This code is called 12ELVIS and can solve a new class of computationally challenging problems in geodynamics, such as shear localization with large strains, crustal intrusion emplacement of magmas, bending of realistic visco-elasto-plastic plates and core-formation by vigorous shell tectonics activities related to a global Rayleigh-Taylor instability of a metal layer formed around silicate-rich lower density (primordial) core during planetary accretion. We discuss in detail the computational strategy required the rheological constraints to be satisfied at each time step and spatial location. We show analytical benchmarks and examples drawn from comparing between numerical and analogue experiments in structural geology, subducting slab bending with a visco-elasto-plastic rheology and equilibrium spherical configurations from self-gravitation. We have also tested possibilities of future applications by addressing 3D geometries based on multigrid method and including inertial effects in the momentum equation with tracers in order to simulate meteoritic impact events and eventually earthquake instabilities. (c) 2007 Elsevier B.V. All rights reserved.
We present a hybrid imaging method for a challenging traveltime tomography problem which includes both unknown medium and unknown scatterers in a bounded domain. The goal is to recover both the medium and the boundary...
详细信息
We present a hybrid imaging method for a challenging traveltime tomography problem which includes both unknown medium and unknown scatterers in a bounded domain. The goal is to recover both the medium and the boundary of the scatterers from the scattering relation data on the domain boundary. Our method is composed of three steps: (1) preprocess the data to classify them into three different categories of measurements corresponding to nonbroken rays, broken-once rays, and others, respectively, (2) use the the nonbroken ray data and an effective data-driven layer stripping strategy-an optimization based iterative imaging method-to recover the medium velocity outside the convex hull of the scatterers, and (3) use selected broken-once ray data to recover the boundary of the scatterers-a direct imaging method. By numerical tests, we show that our hybrid method can recover both the unknown medium and the not-too-concave scatterers efficiently and robustly.
暂无评论