scholarly journals APPLYING THE METHOD OF MAXIMUM CONTRIBUTIONS TO THE MAGNETOGRAM INVERSION TECHNIQUE

2020 ◽  
Vol 6 (4) ◽  
pp. 57-65
Author(s):  
Yury Penskikh
Keyword(s):  
Geomagnetic Field ◽  
Linear Equations ◽  
Original Method ◽  
Spherical Harmonic Analysis ◽  
System Of Linear Equations ◽  
Inversion Technique ◽  
Comprehensive Method ◽  
Maximum Contribution ◽  
Inherent Error ◽  
Ground Magnetic

Fundamentals of the spherical harmonic analysis (SHA) of the geomagnetic field were created by Gauss. They acquired the classical Chapman — Schmidt form in the first half of the XXth century. The SHA method was actively developed for domestic geomagnetology by IZMIRAN, and then, since the start of the space age, by ISTP SB RAS, where SHA became the basis for a comprehensive method of MIT (magnetogram inversion technique). SHA solves the inverse problem of potential theory and calculates sources of geomagnetic field variations (GFV) - internal and external electric currents. The SHA algorithm forms a system of linear equations (SLE), which consists of 3K equations (three components of the geomagnetic field, K is the number of ground magnetic stations). Small changes in the left and (or) right side of such SLE can lead to a significant change in unknown variables. As a result, two consecutive instants of time with almost identical GFV are approximated by significantly different SHA coefficients. This contradicts both logic and real observations of the geomagnetic field. The inherent error of magnetometers, as well as the method for determining GFV, also entails the instability of SLE solution. To solve such SLEs optimally, the method of maximum contribution (MMC) was developed at ISTP SB RAS half a century ago. This paper presents basics of the original method and proposes a number of its modifications that increase the accuracy and (or) speed of solving the SLEs. The advantage of MMC over other popular methods is shown, especially for the Southern Hemisphere of Earth.

scholarly journals APPLYING THE METHOD OF MAXIMUM CONTRIBUTIONS TO THE MAGNETOGRAM INVERSION TECHNIQUE

2020 ◽  
Vol 6 (4) ◽  
pp. 67-76
Author(s):  
Yury Penskikh
Keyword(s):  
Geomagnetic Field ◽  
Linear Equations ◽  
Original Method ◽  
Spherical Harmonic Analysis ◽  
System Of Linear Equations ◽  
Inversion Technique ◽  
Comprehensive Method ◽  
Maximum Contribution ◽  
Inherent Error ◽  
Ground Magnetic

Fundamentals of the spherical harmonic analysis (SHA) of the geomagnetic field were created by Gauss. They acquired the classical Chapman — Schmidt form in the first half of the XXth century. The SHA method was actively developed for domestic geomagnetology by IZMIRAN, and then, since the start of the space age, by ISTP SB RAS, where SHA became the basis for a comprehensive method of MIT (magnetogram inversion technique). SHA solves the inverse problem of potential theory and calculates sources of geomagnetic field variations (GFV) - internal and external electric currents. The SHA algorithm forms a system of linear equations (SLE), which consists of 3K equations (three components of the geomagnetic field, K is the number of ground magnetic stations). Small changes in the left and (or) right side of such SLE can lead to a significant change in unknown variables. As a result, two consecutive instants of time with almost identical GFV are approximated by significantly different SHA coefficients. This contradicts both logic and real observations of the geomagnetic field. The inherent error of magnetometers, as well as the method for determining GFV, also entails the instability of SLE solution. To solve such SLEs optimally, the method of maximum contribution (MMC) was developed at ISTP SB RAS half a century ago. This paper presents basics of the original method and proposes a number of its modifications that increase the accuracy and (or) speed of solving the SLEs. The advantage of MMC over other popular methods is shown, especially for the Southern Hemisphere of Earth.

scholarly journals Approximate solution of system of equations arising in interior-point methods for bound-constrained optimization

2021 ◽  
Author(s):  
David Ek ◽  
Anders Forsgren
Keyword(s):  
Approximate Solution ◽  
Interior Point ◽  
Interior Point Methods ◽  
Linear Equations ◽  
Approximate Solutions ◽  
System Of Nonlinear Equations ◽  
Intermediate Step ◽  
System Of Linear Equations ◽  
Constrained Nonlinear Optimization ◽  
Asymptotic Error

AbstractThe focus in this paper is interior-point methods for bound-constrained nonlinear optimization, where the system of nonlinear equations that arise are solved with Newton’s method. There is a trade-off between solving Newton systems directly, which give high quality solutions, and solving many approximate Newton systems which are computationally less expensive but give lower quality solutions. We propose partial and full approximate solutions to the Newton systems. The specific approximate solution depends on estimates of the active and inactive constraints at the solution. These sets are at each iteration estimated by basic heuristics. The partial approximate solutions are computationally inexpensive, whereas a system of linear equations needs to be solved for the full approximate solution. The size of the system is determined by the estimate of the inactive constraints at the solution. In addition, we motivate and suggest two Newton-like approaches which are based on an intermediate step that consists of the partial approximate solutions. The theoretical setting is introduced and asymptotic error bounds are given. We also give numerical results to investigate the performance of the approximate solutions within and beyond the theoretical framework.

scholarly journals Numerical approximation for the solution of linear sixth order boundary value problems by cubic B-spline

2019 ◽  
Vol 2019 (1) ◽  
Author(s):  
A. Khalid ◽  
M. N. Naeem ◽  
P. Agarwal ◽  
A. Ghaffar ◽  
Z. Ullah ◽  
...  
Keyword(s):  
Numerical Approximation ◽  
Analytic Solution ◽  
Linear Equations ◽  
Computational Cost ◽  
Spline Method ◽  
System Of Linear Equations ◽  
Approximation Solution ◽  
B Spline ◽  
Novel Technique ◽  
Derivatives Of

AbstractIn the current paper, authors proposed a computational model based on the cubic B-spline method to solve linear 6th order BVPs arising in astrophysics. The prescribed method transforms the boundary problem to a system of linear equations. The algorithm we are going to develop in this paper is not only simply the approximation solution of the 6th order BVPs using cubic B-spline, but it also describes the estimated derivatives of 1st order to 6th order of the analytic solution at the same time. This novel technique has lesser computational cost than numerous other techniques and is second order convergent. To show the efficiency of the proposed method, four numerical examples have been tested. The results are described using error tables and graphs and are compared with the results existing in the literature.

scholarly journals Geomagnetic Virtual Observatories: monitoring geomagnetic secular variation with the Swarm satellites

2021 ◽  
Vol 73 (1) ◽  
Author(s):  
Magnus D. Hammer ◽  
Grace A. Cox ◽  
William J. Brown ◽  
Ciarán D. Beggan ◽  
Christopher C. Finlay
Keyword(s):  
Time Series ◽  
Secular Variation ◽  
Geomagnetic Field ◽  
Time Windows ◽  
Low Earth Orbit ◽  
Data Selection ◽  
Spherical Harmonic Analysis ◽  
Core Field ◽  
The Core ◽  
Swarm Satellites

AbstractWe present geomagnetic main field and secular variation time series, at 300 equal-area distributed locations and at 490 km altitude, derived from magnetic field measurements collected by the three Swarm satellites. These Geomagnetic Virtual Observatory (GVO) series provide a convenient means to globally monitor and analyze long-term variations of the geomagnetic field from low-Earth orbit. The series are obtained by robust fits of local Cartesian potential field models to along-track and East–West sums and differences of Swarm satellite data collected within a radius of 700 km of the GVO locations during either 1-monthly or 4-monthly time windows. We describe two GVO data products: (1) ‘Observed Field’ GVO time series, where all observed sources contribute to the estimated values, without any data selection or correction, and (2) ‘Core Field’ GVO time series, where additional data selection is carried out, then de-noising schemes and epoch-by-epoch spherical harmonic analysis are applied to reduce contamination by magnetospheric and ionospheric signals. Secular variation series are provided as annual differences of the Core Field GVOs. We present examples of the resulting Swarm GVO series, assessing their quality through comparisons with ground observatories and geomagnetic field models. In benchmark comparisons with six high-quality mid-to-low latitude ground observatories we find the secular variation of the Core Field GVO field intensities, calculated using annual differences, agrees to an rms of 1.8 nT/yr and 1.2 nT/yr for the 1-monthly and 4-monthly versions, respectively. Regular sampling in space and time, and the availability of data error estimates, makes the GVO series well suited for users wishing to perform data assimilation studies of core dynamics, or to study long-period magnetospheric and ionospheric signals and their induced counterparts. The Swarm GVO time series will be regularly updated, approximately every four months, allowing ready access to the latest secular variation data from the Swarm satellites.

scholarly journals Hybrid Bernstein Block-Pulse Functions Method for Second Kind Integral Equations with Convergence Analysis

2014 ◽  
Vol 2014 ◽  
pp. 1-8 ◽  
Author(s):  
Mohsen Alipour ◽  
Dumitru Baleanu ◽  
Fereshteh Babaei
Keyword(s):  
Integral Equation ◽  
Approximate Solution ◽  
Convergence Analysis ◽  
Bernstein Polynomials ◽  
Linear Equations ◽  
The Other ◽  
New Combination ◽  
System Of Linear Equations ◽  
Numerical Examples ◽  
Block Pulse Functions

We introduce a new combination of Bernstein polynomials (BPs) and Block-Pulse functions (BPFs) on the interval [0, 1]. These functions are suitable for finding an approximate solution of the second kind integral equation. We call this method Hybrid Bernstein Block-Pulse Functions Method (HBBPFM). This method is very simple such that an integral equation is reduced to a system of linear equations. On the other hand, convergence analysis for this method is discussed. The method is computationally very simple and attractive so that numerical examples illustrate the efficiency and accuracy of this method.

Sign in / Sign up

Export Citation Format

Share Document