Solar power mirror arrays for radio astronomy Olaf Wucknitz, Alan - - PowerPoint PPT Presentation

solar power mirror arrays for radio astronomy
SMART_READER_LITE
LIVE PREVIEW

Solar power mirror arrays for radio astronomy Olaf Wucknitz, Alan - - PowerPoint PPT Presentation

Solar power mirror arrays for radio astronomy Olaf Wucknitz, Alan Roy wucknitz@mpifr-bonn.mpg.de aroy@mpifr-bonn.mpg.de Scintillometry Conference, Bonn, 7th November 2019 Solar power mirror arrays for radio astronomy Solar Power Mirror


slide-1
SLIDE 1

Solar power mirror arrays for radio astronomy

Olaf Wucknitz, Alan Roy

wucknitz@mpifr-bonn.mpg.de aroy@mpifr-bonn.mpg.de

Scintillometry Conference, Bonn, 7th November 2019

slide-2
SLIDE 2

Solar power mirror arrays for radio astronomy

  • Solar Power Mirror Arrays
  • Phased array feeds
  • Test case J¨

ulich

  • Simulations of test observations
  • Issues, plans
  • Prospects
  • O. Wucknitz 2019

2/33

slide-3
SLIDE 3

Gemasolar as SKA? (Alan Roy, Ivan Camara, Olaf Wucknitz, …)

slide-4
SLIDE 4
slide-5
SLIDE 5

Gemasolar Basics

Solar field: 2650 heliostats, each 120 m2, total 304 750 m2, equivalent to 620 m diameter single dish Tower height: 140 m Heat-Transfer Fluid: Molten salts (sodium + potassium nitrate) Receiver inlet temp: 290 °C Receiver outlet temp: 565 °C Turbine capacity: 19.9 MW Construction cost: 230 M€ (5 M€ from EU FP5, 80 M€ loan EIB) Timeline: 2007 begin, 2011 online Electricity sales: 110 000 MWh/yr = 30 M€/yr Ownership: Torresol Energy, subsidiary of consortium: 60 % SENER Grupo de Ingeniería (private company, Spain) 40 % MASDAR (alternative energy company of Abu Dhabi)

slide-6
SLIDE 6

The Solar power array problem

  • many mirrors, different delays
  • signal spread over larger area
  • cannot catch the signal with one big feed
  • PAF

⋆ sample focal area ⋆ re-align phases ⋆ scale with signal strength

  • ptimal weights from speckle pattern
  • Need to test concept!
  • O. Wucknitz 2019

6/33

slide-7
SLIDE 7

Some large facilities

[ https://en.wikipedia.org/wiki/Solar_power_tower ] collecting areas Ivanpah 2.6 km2 [ https://solarpaces.nrel.gov/ ] Ashalim 1.1 km2 Crescent Dunes 1.2 km2 Gemasolar 0.3 km2 J¨ ulich 0.018 km2 (> 2× Effelsberg)

  • O. Wucknitz 2019

7/33

slide-8
SLIDE 8

Solar Tower J¨ ulich: 150-m equivalent

  • O. Wucknitz 2019

8/33

slide-9
SLIDE 9

Solar Tower J¨ ulich: Experimental platform

  • O. Wucknitz 2019

9/33

slide-10
SLIDE 10

Practical issues

  • heat (use dedicated tower?)
  • RFI
  • mirrors: do they reflect radio waves?

⋆ expectation: must be thicker than skin depth ⋆ ca. 0.5 – 2 µm for 10 – 1 GHz ⋆ metal mirrors (Gemasolar) okay ⋆ J¨ ulich: 0.2 µm ⋆ actually seems to work!

  • can we predict speckle pattern?
  • beamforming techniques
  • O. Wucknitz 2019

10/33

slide-11
SLIDE 11

Technical setup for test in J¨ ulich

  • tripole antennas (Uppsala, Onsala)
  • simple uncooled receivers, mostly OTS parts
  • Rubidium/GPS clock
  • DBBC3 for sampling (512 MHz bandwidth,
  • max. 6 channels)
  • Mark 6 or dedicated server for recording
  • correlate with Effelsberg for calibration (VLBI)
  • record, correlate, beamform, analyse
  • O. Wucknitz 2019

11/33

slide-12
SLIDE 12

Representative mirror positions for J¨ ulich

200 150 100 50 50 100 150 200 east [m] 50 100 150 200 250 300 north [m]

  • O. Wucknitz 2019

12/33

slide-13
SLIDE 13

Parameters

  • field size ca. D = 300m
  • 2150 mirrors, each 3.2×2.5m2
  • distance to receiver ca. L = 150m (ca. 50 m height)
  • assume 1.5 GHz (20 cm)
  • ‘focus’ size ca. 10 m
  • approximate speckle size

⋆ 0.2 m size ⋆ 1 MHz in frequency ⋆ 10 sec in time ⋆ 0.04 deg on sky

  • O. Wucknitz 2019

13/33

slide-14
SLIDE 14

Speckle image for J¨ ulich experimental platform

4 2 2 4 x [m] 3 2 1 1 2 3 z [m]

  • O. Wucknitz 2019

14/33

slide-15
SLIDE 15

Speckle dynamic spectrum for single feed

100 200 300 400 500 600 700 time [sec] 1500 1510 1520 1530 1540 1550 freq [MHz]

  • O. Wucknitz 2019

15/33

slide-16
SLIDE 16

Instrument response

  • mirror rm, focus/PAF element f p
  • time t and frequency ν generally omitted
  • delays

c τpm = |rm −f p|−θ·rm

  • voltage response for signal E(t)

⋆ time domain Vpm(t) = gpmE(t −τpm) ⋆ freq. domain Vpm = E gpme2πiντpm

  • total voltage response

Vp = E Bp Bp = ∑

m

gpme2πiντpm

  • O. Wucknitz 2019

16/33

slide-17
SLIDE 17

Beamforming theory

  • fit of field E per t,ν

Efit =

p

Bp Vp

p

|Bp|2

  • power estimate from all t,ν

Sfit =

p

Bp Vp

  • 2

p

|Bp|22

  • O. Wucknitz 2019

17/33

slide-18
SLIDE 18

Beamforming result: one antenna element, one sample

20 10 10 20 x [arcmin] 20 10 10 20 y [arcmin]

beam map

0.0 0.5 1.0 1.5 2.0 2.5 3.0 400 200 200 400 U [m] 400 200 200 400 V [m]

FT of beam

10 20 30 40 50

no position from one sample

  • O. Wucknitz 2019

18/33

slide-19
SLIDE 19

Beamforming result: 5 antenna elements, one sample

20 10 10 20 x [arcmin] 20 10 10 20 y [arcmin]

beam map

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 400 200 200 400 U [m] 400 200 200 400 V [m]

FT of beam

10 20 30 40 50

5 elements spread over ∼ 5m

  • O. Wucknitz 2019

19/33

slide-20
SLIDE 20

Beamforming result: 20 antenna elements, one sample

20 10 10 20 x [arcmin] 20 10 10 20 y [arcmin]

beam map

0.0 0.2 0.4 0.6 0.8 1.0 400 200 200 400 U [m] 400 200 200 400 V [m]

FT of beam

10 20 30 40 50

20 elements spread over ∼ 5m

  • O. Wucknitz 2019

20/33

slide-21
SLIDE 21

Beamforming result: 100 antenna elements, one sample

20 10 10 20 x [arcmin] 20 10 10 20 y [arcmin]

beam map

0.0 0.2 0.4 0.6 0.8 1.0 400 200 200 400 U [m] 400 200 200 400 V [m]

FT of beam

10 20 30 40 50

100 elements spread over ∼ 5m

  • O. Wucknitz 2019

21/33

slide-22
SLIDE 22

Beamforming result: one antenna element, 100 samples

20 10 10 20 x [arcmin] 20 10 10 20 y [arcmin]

beam map

0.0 0.2 0.4 0.6 0.8 1.0 400 200 200 400 U [m] 400 200 200 400 V [m]

FT of beam

10 20 30 40 50

10 times (∆t = 10sec), 10 frequencies (∆ν = 1MHz)

  • O. Wucknitz 2019

22/33

slide-23
SLIDE 23

Beamforming result: one antenna element, 2500 samples

20 10 10 20 x [arcmin] 20 10 10 20 y [arcmin]

beam map

0.0 0.2 0.4 0.6 0.8 1.0 400 200 200 400 U [m] 400 200 200 400 V [m]

FT of beam

10 20 30 40 50

50 times (∆t = 10sec), 50 frequencies (∆ν = 1MHz)

  • O. Wucknitz 2019

23/33

slide-24
SLIDE 24

Beamforming result: one antenna element, 10000 samples

20 10 10 20 x [arcmin] 20 10 10 20 y [arcmin]

beam map

0.0 0.2 0.4 0.6 0.8 1.0 400 200 200 400 U [m] 400 200 200 400 V [m]

FT of beam

10 20 30 40 50

100 times (∆t = 10sec), 100 frequencies (∆ν = 1MHz)

  • O. Wucknitz 2019

24/33

slide-25
SLIDE 25

Alternative beamforming

  • formal result

Efit ∝ ∑

p

Bp Vp Bp = ∑

m

gpme2πiντpm

  • reorder

Efit ∝ ∑

m

p

gpme−2πiντpm Vp

  • c τpm = |rm −f p|−θ·rm
  • split delay

Efit ∝ ∑

m

e2πiν θ·rm/c ∑

p

e−2πiν |rm−f p|/c Vp two-stage delay beamformer (optical/analog?)

  • O. Wucknitz 2019

25/33

slide-26
SLIDE 26

Summary

  • solar power array radio telescope may actually work
  • PAF is essential, big PAFs not trivial

ulich: tests in preparation

  • many practical issues to consider
  • multi-beaming provides huge field of view
  • true Square Kilometre Array within reach !?
  • advanced beamforming, only cross-corr?, polarisation
  • build dedicated optimised array?
  • synergy with interstellar scattering/scintillation

see backup slides from Alan Roy (scintillometry 2016)

  • O. Wucknitz 2019

26/33

slide-27
SLIDE 27

Model Mirror Locations

slide-28
SLIDE 28

Dynamic Spectrum: Amplitude

80 MHz

(1400 MHz to 1480 MHz)

12 min

slide-29
SLIDE 29

Dynamic Spectrum: Phase

80 MHz

(1400 MHz to 1480 MHz)

12 min

slide-30
SLIDE 30

Secondary Spectrum

10 μs 1.1 Hz

slide-31
SLIDE 31

Ivanpah Solar Power Facility (USA)

  • O. Wucknitz 2019

31/33

slide-32
SLIDE 32

Cerro Dominador (Chile)

[ https://cerrodominador.com/ ]

  • O. Wucknitz 2019

32/33

slide-33
SLIDE 33

Some Google maps links

https://en.wikipedia.org/wiki/Solar_power_tower

ulich

https://maps.google.de/maps?ll=50.915,6.387778&t=h&z=15

  • Gemasolar

https://maps.google.de/maps?ll=37.558,-5.329&t=h&z=15

  • Crescent Dunes

https://maps.google.de/maps?ll=38.233,-117.366&t=h&z=15

  • Ivanpah

https://maps.google.de/maps?ll=35.57,-115.47&t=h&z=13

  • Cerro Dominador

https://www.google.de/maps?ll=-22.771,-69.485&t=h&z=15

  • O. Wucknitz 2019

33/33