Thinsat Charging and Array Spreading Forces

Thinsats intercept electrons and protons in the van Allen belt; they can be discharged by Fowler-Nordheim tunnelling through carefully aligned arrays of carbon nanotubes at the corners, spaced far enough apart on those corners to maximize the electric field. The charging rate is about 180 nA/m² for protons and 3400 nA/m² for electrons. Thinsats stop low energy protons more effectively, but the low mass electrons move faster and consequently the flux is higher. The thinsat will charge negatively.

Mass-sorted, calibrated CNT in solution will be used, these are commercially available and used for CNT-enabled silicon substrate integrated circuits. Presumably some kind of integrated circuit tool-chip will be used to pluck the CNT macromolecules out of solution, align them, and prepare them to be bonded to a carrier chip. Unlike an electronic circuit, the yield of this bonding process does not need to be perfect or even high; what is most important is that most of the CNT are bonded radially for charge emission into the concentrated field near the corners of the thinsat. Per-CNT currents can be fairly high, and is crudely exponential with voltage; what is most important for reducing voltage is maximizing the electric field gradient and lowering the work function of the CNT tips.

Hydrogen (proton) "contamination" of the CNT tips lowers their work function and lowers thinsat voltage further.

What is the field gradient at the tip of a CNT on the corner of a thinsat? Watcherone (2008) suggests a field gradient enhancement β over the "free space" field at the rounded corner of a graphene square peaks at around 400. The graphene sheet modelled in that paper is 1 nm thick and 1 μm square; This is less sharp than a 70 μm thick and 160 mm square thinsat, so we can presume that the corner β is better for the thinsat.

The CNTs off the corner provide additional field concentration factor; Figure 2 of Eletskii (2010) suggests a β factor of 700 for a CNT with a height to radius factor of 700. If the CNT has a diameter of 2 nm, the height would be about 1.4 μm. If the CNT on each corner are spaced 5 μm apart, that might reduce the β factor to 500. The multiplied β factor might be 2e5. Assuming 50% fill and an array 70 by 70 μm, that is 100 CNT per corner or 400 CNT per thinsat.

A 250 cm² thinsat will collect a net charge of about 80 nA, or about 200 pA per CNT.

Lets assume the free space field is proportional to the thinsat voltage divided by thinsat edge length. If the voltage is 100 V, then the field would be around 500 volts per meter. Times the multiplied beta factor, that is 1e7 volts per meter.

So, what's the voltage?

MoreLater

The capacitance of a square (pointy corner!) thinsat with side S is C = 0.367 × 4πε₀ × S; for a 16 cm thinsat, C₁ = 0.367 × 111.265 pF/m × 0.16 m = 6.530 pF. If the voltage on the thinsat is 200V, the charge is Q_1 = C_1 V = 5223 picoColomb per thinsat.

Assume N thinsats are charged to voltage V_th and are spaced L units apart in an approximately triangular geodesic grid mapped on a sphere. The area of the unit cell is √3/2 L², the area of the sphere is √3/2 L² N = 4 π R², so

R ~=~ \sqrt{ { \Large { \sqrt{ 3 } \over { 8 \pi } } } N } ~ L ~\approx~ 0.26252 \sqrt{ N } ~ L

If the thinsat sphere is 7842 thinsats (icosahedral geodesic sphere with tiling ratio V = 28) and L = 5 meters, then R = 116 meters .

The charge on the whole sphere of thinsats creates a Colomb force on a single thinsat. This resembles the math that shows that the gravitational force at (or above) a spherical shell is like the force from a concentrated mass at the center; we can move the charge to the center, and the effect on a charge at the edge will be the same. This simplifies the math; the force on a thinsat is approximately the simple Colomb force between the charge of N thinsats in the center and the charge of one thinsat distance R away:

If V=800 volts, S = 0.16 meters, and L = 5 meters, then F = 143 nanoNewton. The light pressure on a thinsat (ignoring albedo and thermal re-radiation) is P/c = 1368W / 2.998e8 m/s = 4.56 nN/m² times the S^2 area, or 0.117 nanoNewton, much smaller than the electrostatic pressure. This doesn't look good ...

Suppose we make S 4 times larger, increasing thinsat area by a factor of 16, reducing the tiling ratio V from 28 to 7 and N from 7842 to 492. Thinsat area increases to 0.41 m² and light pressure to 1.87 nanoNewton. That will increase turn and stop time by the square root of the thinsat length, from 15 minutes to about 30 minutes for a 45 degree turn. Further, let's increase the array radius by a factor of 4 (increasing cross-array maneuvering time by 2) and spacing L by a factor of 16. The electrostatic force is reduced by a factor of ( 4/16 )² to 8.9 nanoNewtons. That is still almost 5 times larger than the light pressure ...

MoreLater - how does this actually affect array rotation? what are the tidal forces via the Clohessy-Hill equation, and how does this force affect the array rotation rate?

MoreLater - where does the 800V come from? Can we lower the voltage? this may actually be differential voltage, which we can prevent by design.

MoreLater - a positive voltage on a thinsat implies it is either missing electrons or has stopped some extra protons from the radiation belt that it hasn't boiled off. How many protons is that?


Capacitance of a Square


Possible Mitigations


Guesstimated Thinsat Charging at 6411 km altitude

chargeflux.ods

70.00

μm thinsat thickness

2018/01/06

2.76

mil thinsat thickness

250.00

cm² thinsat area

Proton

Electron

1.602E-10

nA/electron/s

total nA

4.47

83.45

Differential Flux/cm2/s

Depth, mil

Captured Flux/cm2/s

MeV

Proton

Electron

Proton

Electron

Proton

Electron

0.1

1.885E+07

1.734E+09

3.6E-2

3.0E+0

1.9E+7

1.0E+9

0.15

1.780E+07

1.463E+09

5.0E-2

4.0E+0

1.8E+7

7.3E+8

0.2

1.623E+07

8.860E+08

7.5E-2

9.0E+0

1.6E+7

2.3E+8

0.3

1.344E+07

3.787E+08

1.2E-1

1.7E+1

1.3E+7

5.7E+7

0.4

1.105E+07

1.562E+08

1.7E-1

2.6E+1

1.1E+7

1.6E+7

0.5

8.969E+06

6.367E+07

2.2E-1

3.4E+1

9.0E+6

5.0E+6

0.6

7.705E+06

2.586E+07

2.9E-1

4.4E+1

7.7E+6

1.6E+6

0.7

6.529E+06

8.673E+06

3.5E-1

5.5E+1

6.5E+6

4.2E+5

1.0

4.883E+06

3.121E+06

6.0E-1

8.5E+1

4.8E+6

1.0E+5

1.5

3.186E+06

6.610E+05

1.0E+0

1.3E+2

3.0E+6

1.4E+4

2

2.349E+06

9.823E+04

1.8E+0

1.9E+2

1.8E+6

1.4E+3

3

1.285E+06

3.911E+04

3.6E+0

2.9E+2

6.9E+5

3.7E+2

4

8.026E+05

1.546E+04

5.5E+0

3.8E+2

3.2E+5

1.1E+2

5

5.040E+05

7.439E+03

8.0E+0

4.8E+2

1.5E+5

4.3E+1

6

3.101E+05

3.579E+03

1.1E+1

5.7E+2

6.9E+4

1.7E+1

7

1.959E+05

1.889E+03

1.4E+1

6.5E+2

3.5E+4

8.0E+0

10

7.315E+04

9.742E+02

2.7E+1

9.0E+2

7.1E+3

3.0E+0

15

1.514E+04

3.494E+02

4.8E+1

1.2E+3

8.4E+2

8.0E-1

20

3.778E+03

1.251E+02

8.5E+1

1.7E+3

1.2E+2

2.0E-1

30

6.655E+02

5.032E+01

1.8E+2

2.2E+3

1.0E+1

6.3E-2

40

1.601E+02

2.036E+01

3.0E+2

2.7E+3

1.5E+0

2.1E-2

50

7.601E+01

9.270E+00

4.5E+2

3.1E+3

4.6E-1

8.2E-3

60

3.151E+01

5.311E+00

6.3E+2

3.5E+3

1.4E-1

4.2E-3

70

1.957E+01

2.282E+00

8.0E+2

4.0E+3

6.7E-2

1.6E-3

100

8.669E+00

0.000E+00

1.6E+3

4.8E+3

1.5E-2

0.0E+0

150

2.543E+00

0.000E+00

2.8E+3

6.0E+3

2.5E-3

0.0E+0

200

1.045E+00

0.000E+00

5.3E+3

7.0E+3

5.4E-4

0.0E+0

300

2.052E-01

0.000E+00

1.0E+4

8.0E+3

5.7E-5

0.0E+0

400

0.000E+00

0.000E+00

1.6E+4

9.0E+3

0.0E+0

0.0E+0

Flux estimates from SPENVIS 4.6.9.3223 . . Project: RE2

6-Jan 2018

Penetration depth from NASA-HDBK-4002 w/CHANGE-1

graph on pg 28 of 181

ArrayCharging (last edited 2018-02-15 04:21:22 by KeithLofstrom)