A&A Supplement series, Vol. 122, April II 1997, 381-396
Received March 14; accepted July 22, 1996
S. Cuperman - C. Bruma - D. Heristchi
Send offprint request: S. Cuperman
DASOP, Observatoire de Paris, Section Meudon,
F-92195 Meudon Principal Cedex,
France
The three-dimensional (3D) reconstruction of magnetic configurations
above the photosphere is considered within the framework of the
nonlinear force-free-field (FFF) model. The physical- computational algorithm
proposed and tested incorporates, for the first time, the following
basic features:
1) Both photospheric vector field,
and chromospheric line of sight
field component,
data are utilized; this reduces
significantly the degree of ill-posedness characterizing the Cauchy
problem corresponding to the case when only
- values
are used as boundary conditions.
2) A high-order, very efficient computational algorithm is developed and
used: horizontal derivatives are evaluated by 14 - terms
formulas in 14 different forms, selected such as to provide
optimal computational accuracy; the vertical integration
is achieved by the use of ``moving" 10 - term formulas expressed in terms
of 10 derivatives and the first
values (i=x,y,z).
3) At neutral points, where inherent computational singularities in
the values of the FFF-function
arise, rather than using
smoothing techniques based on four-neighbouring- values averages,
suitable procedures ensuring continuity are developed and used.
The overall result of the incorporation of these novel features is an
improvement by orders of magnitude of the accuracy with which the
chromospheric fields are reconstructed in the case in which one uses
(i) only
- values as boundary conditions and
(ii) relative simple computational formulas and smoothing techniques;
at
,
!
The elimination/minimization of measurement errors as well as the
fitting of the corrected date to FFF-model-states is also discussed.
keywords: MHD -- methods: numerical -- Sun: magnetic fields