Issue 
J. Space Weather Space Clim.
Volume 12, 2022



Article Number  29  
Number of page(s)  13  
DOI  https://doi.org/10.1051/swsc/2022018  
Published online  15 August 2022 
Technical Article
Development of accelerated methods for calculating the pattern of current spreading over the surface of spacecraft
HSE University, 20 Myasnitskaya st., 101000 Moscow, Russian Federation
^{*} Corresponding author: avostrikov@hse.ru
Received:
30
June
2021
Accepted:
19
May
2022
Highenergy charged plasma particles pose a danger to space technology. The accumulation of charged particles on the body of the spacecraft generates discharges. Electrostatic discharge is a source of powerful electromagnetic interference that adversely affects the functioning of individual parts and entire systems. According to statistics, in about 30% of cases, the loss of satellites is the consequence of discharges. Before the operation of spacecraft, it is necessary to calculate the spreading of currents, which requires large machine and time costs. The article proposes original approaches for quickly constructing a picture of the spreading of currents over the surface of a spacecraft due to electrification. The key point of the first approach is the construction of a limited area for calculating the flow spreading. The calculation of transient currents will only take place in the electromagnetic compatibility area specified by the user without affecting the rest of it. The paper also developed new simplified computational schemes for a system of differential equations based on the Euler methods. With the help of new computational schemes, the time for calculating unknown quantities in a local area specified by the user has been reduced by several orders of magnitude compared to the calculation of unknown full models. The article presents conclusions on new computational schemes, indicating the complexity of their construction. The adequacy and accuracy of the new computational scheme are confirmed by a practical example.
Key words: Electrification / Electrostatic discharge / Spacecraft / Highlighted area method / Macromodel / Euler methods
© A.V. Vostrikov & E.N. Prokofeva, Published by EDP Sciences 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Work to solve the problem of electrification of spacecraft (SC) began in the early 70s of the 20th century during the development of the geostationary orbit (De Forest & Mc. Jlwain, 1971). By charging a spacecraft, we mean the accumulation of charge on the surfaces of the spacecraft (Garrett & Whittlesey, 2000). The spacecraft is charged as a result of its interaction with the surrounding plasma and solar energetic particles; therefore, charging effects are most pronounced during geomagnetic disturbances (substorms) in the Earth’s magnetosphere (Purvis et al., 1984; Lui, 2000; Lai, 2011). The penetration of the energetic electrons into the internal components or portions of the spacecraft will be stopped and accumulated, leading to high electric fields affecting a nearby victim circuit. This is what causes an electrostatic discharge (ESD). This occurs when the electric fields resulting from the charge exceed the breakdown strength of the outer surface materials of the spacecraft. ESD, a consequence of charging, adversely affects the operation of the onboard electronic equipment of the spacecraft. Electromagnetic interference generated by ESD can reach tens of volts and cause malfunctions in the operation of various onboard systems of a spacecraft, and intense discharge currents can lead to irreversible damage to equipment elements. The main receptors of impulse interference from ESD are fragments of the onboard cable network (OCN) laid along the outer surface of the spacecraft (Garrett, 1981). Table 1 shows some examples of the effects of electrostatic discharges.
Examples of the consequences of electrostatic discharges.
According to statistics, about 30% of the loss of satellites is due to the phenomenon of electrostatic discharge. According to statistics, in a third of the cases, the loss of a spacecraft occurs due to an electrostatic discharge that has arisen on it (Fig. 1). An electrical potential on the order of tens of volts can have detrimental effects on spacecraft components.
Figure 1 Types of spacecraft damage in geostationary orbit. 
Attempts to completely eliminate the possibility of ESD by selecting materials for the outer surface of the spacecraft or active protection of the spacecraft have not yet been very successful (Tyutnev et al., 2019). It is possible to reduce the frequency and power of ESD but not completely eliminate them (Saenko et al., 2015). Therefore, it is necessary to take additional measures for the troublefree operation of the spacecraft electronics when exposed to ESD. The spacecraft electronics must be manufactured in such a way as to guarantee its flawless operation in flight conditions. EMI levels should be assessed in advance during the spacecraft design sketch phase. Then the electronics engineers will have the necessary design information. Thus, the main goal of creating new solutions is aimed at recommendations for minimizing negative phenomena and hazards caused by space weather. To do this, before the operation of the spacecraft, it is necessary to calculate the pickups in the OCN of the spacecraft. Let us present the procedure for determining interference in the OCN of the spacecraft in Figure 2 (the stage considered in the article is marked with a red marker).
Figure 2 The procedure for determining interference in the OCN SC. 
First, identify the most likely locations for electrostatic discharge. The National Aeronautics and Space Administration uses a wellknown computer program for this purpose. NASA Charging Analyzer Program (NASCAP).
Secondly, it is necessary to estimate the levels of impulse noise in the electronic input circuits connected with the outer surface of the spacecraft through the fragments of the cable network. To solve this problem, we use the software “SatelliteMIEM,” (developed by Novikov et al.) is used in Russia. For the software “SattelliteMIEM” the place of the discharge is a priori information obtained earlier. The tasks of the software “SattelliteMIEM” are to build a picture of the spreading of currents over the surface of the spacecraft and to calculate pickups in fragments of the onboard cable network of the spacecraft. This computer system uses Windows 7 or higher. With the help of the SatelliteMIEM software, we calculated the pickups in the cables of the SpektrR space radio telescope, which was launched on July 18, 2011, and served a warranty period in space. A preliminary result of this design phase may be the current distribution resulting from an electrostatic discharge localized at a specified point.
To do this, we used the SatelliteMIEM software (Marchenkov et al., 2008), adapted for modeling and calculating noise signals generated in cable networks and electronic inputs of a spacecraft. This approach is based on the formal representation of the model in the form of a system of differential equations describing thermal, oscillatory, and structural processes in various devices. SatelliteMIEM software for a given 3D model generates a lumped element model (LEM) of the spacecraft, which was originally used to assess the emerging potentials on the spacecraft surface due to cosmic particle fluxes (De Forest & Mc. Jlwain, 1971; Rosen, 1976a). One of the first implementations can be seen in the NASCAP software (H.B. Garret, 1981). NASCAP Software Update Confirms Concerns From NASA Spacecraft Charging (V.A. Davis et al., 2004).
The structural electrophysical model (SEM) of a spacecraft is an equivalent electrical circuit of the spacecraft surface made of lumped R, L, and C elements (Saenko et al., 2015). To construct an SEM, the spacecraft body with attached elements is divided into geometric bodies such as a cylinder, sphere, torus, plane, rod, etc.
SatelliteMIEM software synthesizes SEM based on a given polygonal 3D model of the spacecraft (Fig. 3). A polygonal model consisting of a collection of rectangles or triangles can be easily and quickly built in the 3dsMax package. Then the model is exported to the SatelliteMIEM software, where it is transformed into a surface mesh: a collection of connected nodes. Each connection (branch) is represented in the form of elements of an electric circuit as a whole, forming an equivalent electric circuit of the spacecraft surface.
Figure 3 Converting a collection of quadrangles to an SEM. 
If there is information about possible places of ESD occurrence, it is possible to calculate the spreading pattern of transient currents along the body and attached elements of the spacecraft. In this case, ESDs are presented in the form of a pulsed current source with predetermined characteristics corresponding to the discharge parameters. These sources are connected to points between which ESD may occur. Next, the transient currents are calculated.
To calculate the magnitude of the interference signal in the OCN fragment, the software is used to lay the trace of this fragment along the outer surface of the spacecraft, and the value of the transformation ratio of the current flowing along the spacecraft surface is introduced into the pickup voltage in the OCN (Agapov et al., 2010). The calculated data on the levels of interference signals at the inputs of electronic units, obtained at the stage of the preliminary design of the spacecraft, are included in terms of reference (TOR) for manufacturing these units. In accordance with such TOR, the electronic unit must work without failures at this calculated level of the interference signal.
The cutoff frequency in the ESD spectrum is defined as the rise time of an electrostatic discharge pulse and is 200 MHz for a rise time of 5 ns (equivalent wavelength λ = 1.5 m). For acceptable computational accuracy, the grid cell size is equivalent to the length of the branch on the outer side of the CS. Surface should not exceed λ/10, i.e., 0.15 m. In this case, the use of LEM is permissible. With a typical satellite size of 10 m, the number of grid nodes will approach (1–2)10^{5.}
To calculate the spreading pattern of transient currents along the SEM of the spacecraft, the SatelliteMIEM software uses the wellknown program for calculating electrical circuits PSpice from the OrCAD package. It should be noted that SPICE programs such as OrCad, Pspice, MicroCap, and LTspice compute R, L, and C circuits with this many nodes, using conventional multicore personal computers very slowly. Analysis of an EMS consisting of 150,000 nodes or more on modern PCs in such a program is impossible due to the enormous complexity of the computer spent on solving a system of ordinary differential equations (ODE) by the Runge–Kutta method (Rosen, 1976b). Studies have shown that the time spent solving a system of 2000 or more ODEs in PSpice software is unacceptable for calculations. We have carried out a test analysis of equivalent electrical network (EEN) models of spacecraft of various dimensions on a modern machine. The signal had a duration of 150 ns, and the sampling step was 1 ns. Figure 4 shows the dependence of the computer time consumption for experimental calculations of the spreading of transient currents from ESD to the EEN of the spacecraft on the number of nodes. It can be seen that with an increase in the dimension of the model by n times, the time for calculating the transient currents increases by approximately n^{2} times.
Figure 4 The dependence of the calculation time on the number of nodes in the EEN model. 
The following conclusion was made: the result of spreading currents on large EEN (the number of circuit nodes is 150,000) on modern personal computers takes a lot of computer time (about 80 h). Therefore, it becomes necessary to reduce the time of calculating the spreading of currents over the surface of the spacecraft.
2 Development of the method of highlighted areas
In this paper, we propose methods for the accelerated calculation of the spreading of currents from an electrostatic discharge over the surface of a spacecraft. A feature of the first method is the calculation of transient currents only in the local zone of the spacecraft circuit. In the course of the experiments carried out in the SatelliteMIEM program, an area of significant transient currents of 400 nodes around a given ESR site was identified. Insignificant transient currents can be considered currents that do not exceed 1–2% of the ESD value. Minor currents are considered due to the fact that the electromagnetic interference from such transient currents has too little effect on the onboard electronic equipment of the spacecraft (Garret H.B., 1981; Gaines E.E. et al., 1981).
For an accelerated calculation of the spreading of currents from an electrostatic discharge over the surface of a spacecraft, it is rational to use the method of selected areas developed by us.
The idea is to apply Kirchhoff’s law to every node of the equivalent circuit. The node gives five branches, four of which are connected to other nodes, and the fifth is connected to the ground through a capacitor. Consequently, the current drops sharply after several cell sizes from the ESD of node (node), for example, drops from 100 to 1 A or less. At this point, the cable voltage caused by electrostatic discharge can be disregarded.
The essence of the method consists in constructing a narrow area for calculating the spreading of currents. The calculation of transient currents will only occur in the userlimited EMS area without affecting the rest. To do this, it is proposed to convert the file received as input to the PSpice program. The idea of the method is as follows: based on the found node in which ESD occurred, select the area adjacent to the place of discharge of nodes (Fig. 5). The size of the selected area depends on the value of the error in calculating the spreading of transient currents over the spacecraft body, entered by the user. So, for example, specifying an error equal to 9%, an area equal to 441 nodes around the node where ESD occurred will be selected, and with an error equal to 1%, the selected area will consist of 1681 nodes.
Figure 5 Block diagram of the algorithm for selecting areas. 
When the local area is selected, which means that when the dimension of the model is reduced by n times, the time for calculating the transient currents decreases by n^{3} times. The application of the method in selected areas is justified in the presence of a large spacecraft scheme. The algorithm of the method for calculating the interference in the selected area in the fragments of the spacecraft OCN is shown in Figure 6.
Figure 6 Block diagram of the algorithm of the method for calculating the interference in the selected area in the fragments of the OCN. 
The result of using this method is a significant reduction in the time for calculating transient currents over the surface of the spacecraft in the presence of an error. The error of the currents was calculated on a plane in the region, the values of the transient currents in the branches of which are no more than 1–2% of the ESR value. It was found that with the expansion of the selected area, and hence, with an increase in the number of nodes around a given ESD, the error decreases (Fig. 7). When 1681 nodes are selected around the discharge point, the error becomes insignificant (equal to 1%) and is quite acceptable for real calculations. The main advantage of this approach is the ability to read largescale RLC circuits on modern computers in a short time.
Figure 7 The decrease in the error with an increase in the number of nodes in the selected area. 
This method is based on a heuristic approach and is used when the accuracy of calculations is 3–5%. The main idea is to apply Kirchhoff’s law to every node of an equivalent electrical circuit. The node gives five branches, four of which are connected to other nodes, and the fifth is connected to the ground through a capacitor. Consequently, the current drops sharply after a few cell sizes from the ESD node, for example, drops from 100 to 1 A or less. At this point, the cable voltage caused by electrostatic discharge can be disregarded.
The main disadvantage of this method is that thirdparty software (PSpice) is required for its implementation. In this work, we propose another method that does not have this drawback.
2.1 Development of a macromodel of a spacecraft based on explicit and implicit Euler methods
Let us present the problem before us in the language of mathematics. SEM of the spacecraft can be formed in an extended homogeneous coordinate basis and written in the form of a system of linear ordinary differential equations
where are C, G – numeric order matrices (n × n), moreover C – capacitance and inductance matrix, G – resistance matrix, – the vector of the sought phase variables (voltages at all nodes of the circuit and currents flowing through inductive elements), – vector of inputs.
It is necessary to carry out the solution of the system of equations at the moment of time with minimal time expenditures t*, i.e., calculate the numerical vector .
To calculate the transient process of the system of equations, it is necessary to solve the system (1). It is possible to use various methods for solving the problem, numerical, numericalanalytical, based on macromodeling, etc.
Explicit Euler’s Method for solving a system of equations
can be obtained by expanding the solution into a truncated Taylor series in the neighborhood h_{i} any point t_{i}
From where
The implicit Euler method is also derived from the expansion of the solution into a truncated Taylor series
From where
Let us show the idea of constructing a macromodel using Euler’s formulas and assuming that the SEM of the spacecraft is devoid of any specificity; that is, the model will be presented in a general form. Let us represent the mathematical model of the spacecraft SEM in block form
where Q are variable parameters, matrices G and C are nondegenerate. After specifying the discharge points, the vector will have only two nonzero elements in the subvector .
The macromodeling algorithm consists of the preliminary transformation of problems (2) and (3) into some reduced spacecraft model, the number of equations in which will be equal to m ≪ n, which reduces the complexity of the calculations.
Substituting expressions for derivatives (2) and (3) into the model (4), we obtain a system of 4 equations
Let us carry out a series of substitutions and transformations. The aim of the transformation is to reduce the system (4). Substitute (4.2) into (4.4), then substitute (4.2) into (4.1), then substitute the result of the second substitution into the result of the first substitution. The result will be a new macro model or computational scheme (5):
Based on the fact that the subvector is zero, then after multiplying all matrices in (5), the final expression will have the form
where are A_{1}, A_{2}, A_{3}, A_{4} – numeric (m × m) – matrices.
But the model (1) considered by the authors has a number of peculiarities. Let us note the specifics of the scheme model.
The matrix is diagonal and degenerate. The diagonal consists of two groups of coefficients, and the coefficients within each group are the same.
Matrix  nondegenerate, symmetric, and sparse matrix.
Vector contains only one or two nonzero coefficients of the form y_{i}t, where the subscript “i” can be any number from 1 to n.
We will assume that the original problem has the following form:
Taking into account the specific form of the system of equations (7), we obtain the following formula of the Euler method
Taking into account the specific form of the system of equations (6), we obtain the following formula of the Euler method
The algorithm of macromodeling developed by us consists of the preliminary transformation of problems (8) and (9) into some reduced spacecraft model, the number of equations which will be equal, which reduces the complexity of calculations. A feature of the method is the calculation of transient currents only in the local zone of the spacecraft circuit. As previously indicated, this action is correct due to too little influence on the spacecraft radioelectronic equipment of the calculated values of electromagnetic interference outside the selected area equal to 400 nodes. Therefore, for the accelerated calculation of the spreading of currents from ESD over the surface of the spacecraft, it is rational to use the method proposed in this work. To do this, we write both types of computational processes in block form:
where , . The specificity is that block matrix A is sparse and has a high order. In this regard, it is proposed, before the second stage of model reduction, to transform the matrix A to a triangular form with a border using the method of determining quantities SangiovanniVincentelli & Bickart (1979). After transforming the matrix A we get the blocks: A_{11} is a nondefective lower triangular matrix of order (n*n), A_{12} – matrix (n*m), A_{21} – matrix (m*n), A_{22} – matrix (m*m).
We will assume that the subvectors and contain on m ≪ n the desired solution coefficients. It is necessary, on the basis of expressions (10) and (11), to formulate a new computational scheme from which the subvectors and .
For this, we write (10) and (11) in the form of the following subsystems of equations
Expressions (10.1), (10.2), (11.1), and (11.2) will be considered as an analog of a system of four equations with four unknowns , , , , based, generally speaking, on computational processes of different types (explicit and implicit).
Now we substitute the obtained expression (12) into (10.2).
After multiplying all matrices in (12), the final expression will have the form
where B_{1}, B_{2}, B_{3} – numeric (m × m) – matrices, B_{4} – numeric (m × n) – matrix. In practical problems, the vector contains only a few nonzero coefficients. Therefore, in practice, expression (14) will have a simpler form
By construction, the matrix B_{1} nondegenerate and independent of the parameter h. Therefore, the inverse matrix can be calculated from it, after multiplying by which expressions (15), we obtain a computational scheme in the form
As a result, a macromodel was built, and a new computational scheme was obtained, in which only m ≪ n vector coefficients .
The complexity T of the process of constructing a computational scheme (13):
The complexity of calculations according to the scheme (16) with a changed step will be 3m^{2} real multiplicative operations (RMO), and with a constant step – m^{2} RMO. If n = 1,50,000 and for dense matrices, the computation speed will increase 7,031,250,000 times as compared to the implicit Euler method and 46,875 times as compared to the explicit Euler’s method.
The application of this approach to the macromodeling of the EES of the spacecraft is justified in the presence of a spacecraft of a large dimension, with a significant time interval t^{*}, that is when it is required to calculate the numerical vector , having taken a large number of steps.
According to the proposed method, test calculations were carried out and compared with calculations performed by traditional methods. A comparison of the calculation results showed their complete adequacy.
2.2 Comparison of time costs
The practice has shown that calculating large mathematical models according to the computational scheme is 2–3 orders of magnitude faster than the methods used in specialized software. The time for analyzing an EPS, consisting of 1000 nodes, by different methods with a step h = 0.001 s and the number of steps 1000 is given in Table 2, and the time for constructing a reduced computational scheme is given in Table 2. The experiments were carried out on a PC with a dualcore processor (clock frequency 1.8 GHz per core) and 2 GB of RAM. In this case, the results of the calculation by the reduced computational scheme differ from the results of the implicit method ≪ 1%. Note that the time of constructing the macromodel is comparable to the time of a single analysis of the circuit by the implicit Euler method (Table 3). Thus, the use of the proposed reduced scheme is justified in case of its repeated use.
Time for calculating EES by different methods.
Time of building a reduced computational scheme containing a different number of nodes in a subvector .
2.3 Practical example
The design scheme was tested on large problems, but it does not seem rational to illustrate these examples due to a large amount of data. To confirm the results of the operability of the given computational circuit, we will conduct an experiment without any impact on the investigated electrical circuit in the laboratory (room conditions). To do this, on the model, we will build a test equivalent electrical circuit consisting of resistors with a nominal value of 100 Ω, capacitors of 100 pF, and intactivities of 1 μH each, as a result, we get a circuit consisting of 23 nodes. We will connect a signal generator to the 15th node of the assembled circuits, to which we will supply a pulsed periodic signal (meander) as a disturbing effect; we will also connect an oscilloscope to fix the potentials in the circuit nodes. The meander has the following characteristics (Fig. 8): Von [V] = 5, T_{rise} [s] = 0.3e–6, T_{fall} [s] = 0.3e–6, T_{on} [s] = 0.3e–6, Tperiod [s] = 1e–6. In Figure 8, the oscilloscope captures signals at nodes 15 and 23 of the circuit. Figure 9 shows the applied signal (meander), while the abscissa shows the time with a grid step of 0.2e–6 s, and the yaxis shows the voltage with a grid step of 0.5 V.
Figure 8 Signal captured at node 15. 
Figure 9 EES containing 23 nodes. 
In the future, we will obtain a similar result with the same electrical circuit in the LTSpice program, the adequacy of which in the field of computer modeling of digital and analog and digital electrical circuits has been confirmed by world experience (Fig. 9). Figure 9 shows the applied signal (meander), while the abscissa shows the time with a grid step of 0.2e–6 s, and the yaxis shows the voltage with a grid step of 0.5 V. Figure 9 below shows the diagram assembled in the LTSpice software, corresponding to the layout. Applying the tester to each of the 23 nodes, we get 23 marked signals. The signals obtained in the LTSpice software completely coincide with the experiment.
Next, we will carry out a similar calculation with the same electrical circuit according to the proposed reduction method. Let us evaluate the results obtained after analyzing the scheme with various computational schemes: the explicit Euler method, the implicit Euler method, and the proposed reduced method. Let the perturbation be associated with the 15th node, and the desired local area for calculation according to the proposed calculation scheme is located at nodes: 15, 16, and 17. Let us estimate the values of the potentials at the nodes at one point in time. Moment 400 ns by various calculation methods. (Table 4).
Calculation with h = 0.001 s and number of steps = 1000, moment of time = 400 ns.
The result of calculating the circuit using the LTSpice software coincides with the result of calculating the implicit Euler method. The analysis results obtained by the implicit Euler method and the reduced computational scheme differ by 0.0006%. At the stage of the preliminary design of the pattern of current spreading through the body of the spacecraft, such an error is more than acceptable (Table 4).
3 Conclusion

The research questions or problem defined. We considered the process of charging spacecraft and the electrostatic discharge associated with it. Then we examined the consequences of an electrostatic discharge on the onboard radioelectronic equipment of a spacecraft and presented statistics according to which up to a third of spacecraft fail due to charging. Therefore, before launching a spacecraft into outer space, it is necessary to carry out calculations of pickups delivered to the input of electronic units. Next, we discussed the existing software for calculating the current spreading pattern. In this case, the calculation takes too much time on a conventional personal computer. The urgent problem is to develop a new accelerated method for constructing a pattern of current spreading over the spacecraft body.

Data and assumptions. To solve the problem posed, it is proposed to reduce the system of equations describing the process of current spreading based on the wellknown explicit Euler method and implicit Euler method. The computer time of system transformation should be no more than the time of analysis by traditional methods. The calculation according to the proposed scheme should be faster than traditional methods.

Rationale of the proposed method. To justify the proposed method, we write the electrical circuit model in the extended homogeneous coordinate basis (1). Next, we write the computational process in the form of explicit and implicit Euler methods (4). After the reduction of computational processes, a macromodel or a new computational scheme (6) was obtained. Taking into account the peculiarities, the new computational scheme takes the form (13) or (16).

Results. To confirm the performance of the new computational scheme (16), we conducted a laboratory experiment. To do this, we assembled a working circuit from resistors, inductors, capacitors, and a source of disturbing influences. The scheme consisted of 23 nodes. A meander was applied as a disturbing influence. The result of the signal passing through the circuit is shown in Figure 8. Further, to check the adequacy of the result, the same circuit was built in the LTSpice simulator. We obtained the values of the currents in the nodes of the circuit at a time instant of 400 ns, then the same result with an error of 0.0006% was obtained using a new computational scheme.

Assessment based on author’s validation. The calculation time for the new computational scheme is shown in Table 2. For the 23node scheme, the implicit Euler method took 11.4 s; the explicit Euler method took 1057.01 s, and the new computational scheme 0.063 s. The time for calculating unknown quantities in a userspecified local area has been reduced by several orders of magnitude compared to the calculation of unknown full models. It took 851.4 s to build a new computational circuit. The calculation error according to the proposed computational scheme was 0.0006%; such an error is more than acceptable. Thus, the use of a new computational scheme is justified if it is used repeatedly.
Acknowledgments
This article is an output of a research project implemented as part of the Basic Research Program at the National Research University Higher School of Economics (HSE University). The editor thanks two anonymous reviewers for their assistance in evaluating this paper.
References
 Agapov VV, Marchenkov KV, Saenko VS, Sokolov AB. 2010. Method for determining the transformation coefficient converting surface structure currents into voltages in fragments of cable harnesses. RU Patent 2 378 657 C2. [Google Scholar]
 Davis VA, Mandell MJ, Gardner BM, Mikellides IG, Neergaard Jacobs LF, Cooke DL, Minor J. 2004. Validation of NASCAP2K Spacecraftenvironment interactions calculations. [Google Scholar]
 De Forest SE, Mc Jlwain GE. 1971. Plasma clouds in the magnetosphere. J. Geophys. Res. 76: 3587–3611. https://doi.org/10.1029/JA076i016p03587. [CrossRef] [Google Scholar]
 Gaines EE, Nightingale RW, Imhof WL, Reagan JB. 1981. Enhanced Radiation Doses to HighAltitude Spacecraft during June 1980. IEEE Trans Nuclear Sci 28(6): 4502–4504. [CrossRef] [Google Scholar]
 Garrett HB. 1981. The Charging of Spacecraft Surfaces. Rev Geophys 19(4): 577–616. ISSN: 87551209. https://doi.org/10.1029/RG019i004p00577. [CrossRef] [Google Scholar]
 Garrett HB, Whittlesey AC. 2000. Spacecraft charging, an update. IEEE Trans Plasma Sci 28(6): 2017–2028. https://doi.org/10.1109/27.902229. [CrossRef] [Google Scholar]
 Lai ST. 2011. Fundamentals of spacecraft charging: spacecraft interactions with space plasmas. Princeton University Press. pp. 272. ISBN 1400839092, 9781400839094. [Google Scholar]
 Lui ATY. 2000. Tutorial on geomagnetic storms and substorms. IEEE Trans Plasma Sci 28(6): 1854–1866. https://doi.org/10.1109/27.902214. [CrossRef] [Google Scholar]
 Marchenkov KV, Sokolov AB, Saenko VS, Pozhidaev ED. 2008. A newgeneration ‘SatelliteMIEM’ software for calculation of interference in fragments of the cable network laid on the external surface of the space vehicle (in Russian). Proc Technol Electromagn Compat, Moscow, Russia 1: 39–44. [Google Scholar]
 Purvis CK, Garrett HB, Whittlesey AC, Stevens NJ. 1984. Design guidelines for assessing and controlling spacecraft charging effects. Tech. Paper 2361, NASA, Washington, DC, USA, pp. 1–44. [Google Scholar]
 Robinson PA, Holman Jr AB. 1977. Pioneer Venus spacecraft charging model. In: Proceedings of the Spacecraft Charging Technology Conference, February 1977, USA, pp. 297–308. [Google Scholar]
 Rosen A. 1976a. Spacecraft charging by magnetospheric plasmas. IEEE Trans Nuclear Sci 23(6): Electronic ISSN: 15581578. https://doi.org/10.1109/TNS.1976.4328575. [Google Scholar]
 Rosen A. 1976b. Spacecraft charging problems. physics of solar planetary environments. In: Proceedings of the International Symposium on SolarTerrestrial Physics, June 7–18, 1976, Boulder, Colorado, Vol. II, pp. 1024–1038. https://doi.org/10.1029/SP008p1024. [Google Scholar]
 Saenko VS, Tyutnev AP, Nikolski EV, Bakutov AE. 2015. Protection of the SpectrR spacecraft against ESD effects using satelliteMIEM computer code IEEE Trans Plasma Sci 43(9): 2828–2831. https://doi.org/10.1109/TPS.2015.2403955. [CrossRef] [Google Scholar]
 SangiovanniVincentelli A, Bickart T. October, 1979. Bipartite graphs and an optimal bordered triangular form of a matrix. IEEE Trans Circuits Syst 26(10): 880–889. [CrossRef] [Google Scholar]
 Tyutnev A, Saenko V, Ikhsanov R, Krouk E. 2019. Radiationinduced conductivity in polymers under pulsed and longtime smallsignal irradiations combined to determine their stepfunction response. J Appl Phys 126(9): 095501. ISSN: 00218979. https://doi.org/10.1063/1.5109768. [CrossRef] [Google Scholar]
Cite this article as: Vostrikov AV & Prokofeva EN. 2022. Development of accelerated methods for calculating the pattern of current spreading over the surface of spacecraft. J. Space Weather Space Clim. 12, 29. https://doi.org/10.1051/swsc/2022018.
All Tables
Time of building a reduced computational scheme containing a different number of nodes in a subvector .
Calculation with h = 0.001 s and number of steps = 1000, moment of time = 400 ns.
All Figures
Figure 1 Types of spacecraft damage in geostationary orbit. 

In the text 
Figure 2 The procedure for determining interference in the OCN SC. 

In the text 
Figure 3 Converting a collection of quadrangles to an SEM. 

In the text 
Figure 4 The dependence of the calculation time on the number of nodes in the EEN model. 

In the text 
Figure 5 Block diagram of the algorithm for selecting areas. 

In the text 
Figure 6 Block diagram of the algorithm of the method for calculating the interference in the selected area in the fragments of the OCN. 

In the text 
Figure 7 The decrease in the error with an increase in the number of nodes in the selected area. 

In the text 
Figure 8 Signal captured at node 15. 

In the text 
Figure 9 EES containing 23 nodes. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.