(Ebook PDF) ELECTRIMACS 2019 Selected Papers Volume 1 1st Edition by Walter Zamboni, Giovanni Petrone 3030371611 9783030371616 Full Chapters Download
(Ebook PDF) ELECTRIMACS 2019 Selected Papers Volume 1 1st Edition by Walter Zamboni, Giovanni Petrone 3030371611 9783030371616 Full Chapters Download
https://ebookball.com/product/ebook-pdf-consciousness-based-
leadership-and-management-volume-1-1st-edition-by-anil-
maheshwari-3031062345-9783031062346-full-chapters-22016/
https://ebookball.com/product/ebook-pdf-gate-2019-electrical-
engineering-2019-1st-edition-by-trishna-knowledge-
systems-9352868439-978-9352868438-full-chapters-22770/
https://ebookball.com/product/ebook-pdf-strategic-human-resource-
management-volume-1-text-and-cases-1st-edition-by-feza-tabassum-
azmi-1108642209-9781108642200-full-chapters-20610/
Seminars in Orthodontics March 2019 Volume 25 Issue 1 1st edition by
Peter Ngan ISSN 1073-8746
https://ebookball.com/product/seminars-in-orthodontics-
march-2019-volume-25-issue-1-1st-edition-by-peter-ngan-
issn-1073-8746-4056/
https://ebookball.com/product/journal-of-endodontics-
january-2019-volume-45-issue-1-1st-edition-by-the-journal-of-
endodontics-issn-00992399-3888/
https://ebookball.com/product/ebook-pdf-managing-healthcare-
ethically-volume-1-leadership-roles-and-responsibilities-3rd-
edition-by-william-nelson-paul-
hofmann-1640552472-9781640552470-full-chapters-21952/
https://ebookball.com/product/ebook-pdf-biomedical-engineering-
technologies-volume-1-1st-edition-by-miguel-ossandon-houston-
baker-avraham-rasooly-1071618024aeurz-978-1071618028-full-
chapters-22416/
Literature and Society Selected Papers from the English Institute 1978
1st edition by Professor Edward Said ISBN 0801822947 9780801822940
https://ebookball.com/product/literature-and-society-selected-
papers-from-the-english-institute-1978-1st-edition-by-professor-
edward-said-isbn-0801822947-9780801822940-19098/
Lecture Notes in Electrical Engineering 604
Walter Zamboni
Giovanni Petrone Editors
ELECTRIMACS
2019
Selected Papers - Volume 1
Lecture Notes in Electrical Engineering
Volume 604
Series Editors
Leopoldo Angrisani, Department of Electrical and Information Technologies Engineering, University of Napoli
Federico II, Naples, Italy
Marco Arteaga, Departament de Control y Robótica, Universidad Nacional Autónoma de México, Coyoacán,
Mexico
Bijaya Ketan Panigrahi, Electrical Engineering, Indian Institute of Technology Delhi, New Delhi, Delhi, India
Samarjit Chakraborty, Fakultät für Elektrotechnik und Informationstechnik, TU München, Munich, Germany
Jiming Chen, Zhejiang University, Hangzhou, Zhejiang, China
Shanben Chen, Materials Science and Engineering, Shanghai Jiao Tong University, Shanghai, China
Tan Kay Chen, Department of Electrical and Computer Engineering, National University of Singapore,
Singapore, Singapore
Rüdiger Dillmann, Humanoids and Intelligent Systems Lab, Karlsruhe Institute for Technology, Karlsruhe,
Baden-Württemberg, Germany
Haibin Duan, Beijing University of Aeronautics and Astronautics, Beijing, China
Gianluigi Ferrari, Università di Parma, Parma, Italy
Manuel Ferre, Centre for Automation and Robotics CAR (UPM-CSIC), Universidad Politécnica de Madrid,
Madrid, Spain
Sandra Hirche, Department of Electrical Engineering and Information Science, Technische Universität
München, Munich, Germany
Faryar Jabbari, Department of Mechanical and Aerospace Engineering, University of California, Irvine, CA,
USA
Limin Jia, State Key Laboratory of Rail Traffic Control and Safety, Beijing Jiaotong University, Beijing, China
Janusz Kacprzyk, Systems Research Institute, Polish Academy of Sciences, Warsaw, Poland
Alaa Khamis, German University in Egypt El Tagamoa El Khames, New Cairo City, Egypt
Torsten Kroeger, Stanford University, Stanford, CA, USA
Qilian Liang, Department of Electrical Engineering, University of Texas at Arlington, Arlington, TX, USA
Ferran Martin, Departament d’Enginyeria Electrònica, Universitat Autònoma de Barcelona, Bellaterra,
Barcelona, Spain
Tan Cher Ming, College of Engineering, Nanyang Technological University, Singapore, Singapore
Wolfgang Minker, Institute of Information Technology, University of Ulm, Ulm, Germany
Pradeep Misra, Department of Electrical Engineering, Wright State University, Dayton, OH, USA
Sebastian Möller, Quality and Usability Lab, TU Berlin, Berlin, Germany
Subhas Mukhopadhyay, School of Engineering & Advanced Technology, Massey University,
Palmerston North, Manawatu-Wanganui, New Zealand
Cun-Zheng Ning, Electrical Engineering, Arizona State University, Tempe, AZ, USA
Toyoaki Nishida, Graduate School of Informatics, Kyoto University, Kyoto, Japan
Federica Pascucci, Dipartimento di Ingegneria, Università degli Studi “Roma Tre”, Rome, Italy
Yong Qin, State Key Laboratory of Rail Traffic Control and Safety, Beijing Jiaotong University, Beijing, China
Gan Woon Seng, School of Electrical & Electronic Engineering, Nanyang Technological University,
Singapore, Singapore
Joachim Speidel, Institute of Telecommunications, Universität Stuttgart, Stuttgart, Baden-Württemberg,
Germany
Germano Veiga, Campus da FEUP, INESC Porto, Porto, Portugal
Haitao Wu, Academy of Opto-electronics, Chinese Academy of Sciences, Beijing, China
Junjie James Zhang, Charlotte, NC, USA
The book series Lecture Notes in Electrical Engineering (LNEE) publishes the latest developments
in Electrical Engineering - quickly, informally and in high quality. While original research reported
in proceedings and monographs has traditionally formed the core of LNEE, we also encourage
authors to submit books devoted to supporting student education and professional training in
the various fields and applications areas of electrical engineering. The series cover classical and
emerging topics concerning:
• Communication Engineering, Information Theory and Networks
• Electronics Engineering and Microelectronics
• Signal, Image and Speech Processing
• Wireless and Mobile Communication
• Circuits and Systems
• Energy Systems, Power Electronics and Electrical Machines
• Electro-optical Engineering
• Instrumentation Engineering
• Avionics Engineering
• Control Systems
• Internet-of-Things and Cybersecurity
• Biomedical Devices, MEMS and NEMS
For general information about this book series, comments or suggestions, please contact leontina.
[email protected].
To submit a proposal or request further information, please contact the Publishing Editor in
your country:
China
Jasmine Dou, Associate Editor ([email protected])
India
Swati Meherishi, Executive Editor ([email protected])
Aninda Bose, Senior Editor ([email protected])
Japan
Takeyuki Yonezawa, Editorial Director ([email protected])
South Korea
Smith (Ahram) Chae, Editor ([email protected])
Southeast Asia
Ramesh Nath Premnath, Editor ([email protected])
USA, Canada:
Michael Luby, Senior Editor ([email protected])
** Indexing: The books of this series are submitted to ISI Proceedings, EI-Compendex,
SCOPUS, MetaPress, Web of Science and Springerlink **
ELECTRIMACS 2019
Selected Papers - Volume 1
Editors
Walter Zamboni Giovanni Petrone
DIEM DIEM
Università degli Studi di Salerno Università degli Studi di Salerno
Fisciano, Italy Fisciano, Italy
This Springer imprint is published by the registered company Springer Nature Switzerland AG.
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Preface
v
vi Preface
based devices, and safety issues in electrical engineering contain papers focussed
on new technologies for renewable energy, storage systems and electronic devices,
mainly selected from the special session papers.
Organising Committee
General Chairs
Giovanni Petrone, Università degli studi di Salerno, Italy
Walter Zamboni, Università degli studi di Salerno, Italy
Technical Program Chairs
Giovanni Spagnuolo, Università degli studi di Salerno, Italy
Éric Monmasson, Université de Cergy-Pontoise, France
Benoît Robyns, HEI Lille, France
Ramon Blasco-Gimenez, Universitat Politécnica de Valéncia, Spain
Track Chairs
Federico Baronti, Università di Pisa, Italy
Efstratios Batzelis, Imperial College, London, United Kingdom
Pavle Boškoski, Jožef Stefan Institute, Ljubljana, Slovenia
Mario Cacciato, Università degli studi di Catania, Italy
Maria Carmela Di Piazza, ISSIA-CNR Palermo, Italy
Seiichiro Katsura, Keio University, Japan
Eftichios Koutroulis, Technical university of Crete, Greece
Marie-Cecile Pera, Université de Franche-Comté, France
Carlos Andres Ramos Paja, Universidad National de Colombia, Colombia
Bruno Sareni, ENSEEIHT, France
João Pedro Trovão, University of Sherbrooke, Canada
Dmitri Vinnikov, Tallinn University of Technology Estonia
Special Session Chair
Ilhem Slama-Belkodja, Université de Tunis El Manar, Tunisia
Local Organising Committee
Gina Scorziello
Raffaele Raimo
Andrea Contrada
Antonio Guarino
Rudy Alexis Guejia Burbano
Vittorio Mattei
Luigi Mattia
Brian Ospina Agudelo
Carmine Russomando
Preface vii
Scientific Committee
Yacine Amara, France
Federico Baronti, Italy
Efstratios Batzelis, United Kingdom
Ramon Blasco-Gimenez, Spain
Pavle Boskoski, Slovenia
Alain Bouscayrol, France
Mario Cacciato, Italy
Carlo Cecati, Italy
Chandan Chakraborty, India
Gerard Champenois, France
Bruno Dehez, Belgium
Louis-A. Dessaint, Canada
Maria Carmela Di Piazza, Italy
Christian Dufour, Canada
Maurice Fadel, France
Gabriel Garcerá, Spain
Leopoldo García Franquelo, Spain
Guillaume Gateau, France
Luis Gomes, Portugal
Gabriele Grandi, Italy
Francesco Grasso, Italy
Hamid Gualous, France
Lennart Harnefors, Sweden
Sergio Junco, Argentina
Ðani Juričić, Slovenia
Hadi Y. Kanaan, Lebanon
Seiichiro Katsura, Japan
Samir Kouro, Chile
Eftichios Koutroulis, Greece
Marco Liserre, Germany
Luc Loron, France
Óscar Lucía, Spain
Massimiliano Luna, Italy
Chengbin Ma, China
Jean Mahseredjian, Canada
Mariusz Malinowski, Poland
Patrizio Manganiello, Belgium
Sébastien Mariethoz, France
Fabrizio Marignetti, Italy
Bogdan Marinescu, France
Pascal Maussion, France
Tuomas Messo, Finland
Rosario Miceli, Italy
viii Preface
Guest Editors
Fisciano, Italy Giovanni Petrone
Walter Zamboni
Contents
ix
x Contents
Walter Zamboni received the PhD degree in Electrical Engineering from the Uni-
versità degli studi di Napoli Federico II, Italy, in 2004. From 2016 he is Associate
Professor of Electrical Engineering at the Università degli studi di Salerno, Italy,
where he served as assistant professor from 2008 to 2016. He is a member of IEEE
and IMACS. At present, his main scientific interests include battery modelling and
SoC–SoH estimation, identification and diagnostics of batteries and fuel cells, and
energy and power management in systems with energy storage. He co-authored
more than 80 papers published in international journals and conference proceedings,
and one patent.
Giovanni Petrone is Associate Professor of Circuit Theory at the University of
Salerno, Italy. His research areas include analysis and design of power electronics
and controls for photovoltaic, fuel cell and wind systems, tolerance analysis
of electronic circuits, non-linear control techniques, wireless power transfer and
identification and diagnosis of renewable sources and storage systems. He is
Associate Editor of IEEE Journal of Photovoltaics and of IET Power Electronics.
He is member of the editorial board of MPDI Applied Sciences journal. Since 2017,
he is a senior member of IEEE. Prof. Petrone is co-author of more than 150 papers
published in international journals and conference proceedings. He is also co-author
of two books, two IEEE e-Learning library courses and five patents in the field of
power electronics for photovoltaic applications.
xv
Part I
Modelling and Simulation
Modelling of Photovoltaic Systems for
Real-Time Hardware Simulation
Abstract The real-time simulation is a valid help to test electrical systems when
a physical device is not available. This is significantly evident when used in
hardware and software co-simulation environment, where it is possible to connect
the emulator to a real subsystem to test or validate it. In this paper, a model of
the photovoltaic system is presented that can be implemented within a hardware
simulator to be able to interface it with a real circuit, the hardware simulator used is
the National Instruments RIO system.
1 Introduction
Real-time simulation (RT) is one of the frontiers of electrical system analysis. Its
task is to control or react to events that occur in the real world, with the same speed
as the real world clock. With the aim of having systems that can be interfaced with
other software or real systems in order to simulate increasingly complex systems
and on a large scale, also taking communication infrastructure into account, as in
the case of smart grids [1–4].
It is necessary to be aware of the execution time of the simulation system
to ensure that it is shorter or equivalent to the selected time step especially in
simulating scenarios that cannot normally be performed, for example, to model
the high switching characteristics of the converters, a fast processor is required to
reach the minimum step size of up to microseconds [5, 6]. To simulate very quick
phenomena in the order of microseconds, it is necessary to think of architectures
with more than one processor, and based on FPGAs, that allow high sampling rates,
a faster and parallel computing capacity [7]. As the penetration of PV generation
increases, its impact on stability and security of the power system will become
more and more significant, due to the characteristic of randomness and volatility
[8]. Modeling and simulation are the basic technologies to study the impact on the
power grid in which, large-scale PV generation systems are integrated. Simulations
in the usual platform may give good results, but they are not able to deliver the
results for dynamic change in input as present in real world in run time, the model
may not respond for such a change. When we try to simulate to know the long-
term performance of a system, the normal simulation requires a very long time to
deliver the results and the accuracy of the results may also get compromised. While
simulating the complex model like PV system interaction with the grid in real-time,
we may encounter many problems. The important ones to mention are,
1. The need for a PV mathematical model that can deliver results faster to keep the
real-time simulation properties during execution. Solving the algebraic loop in
the PV model is an important task, as algebraic loops are not supported in the
real-time hardware.
2. There is a need of cost-effective test bench/platform for simulating PV systems
in real-time that can be used for control validation, studies of the storage system,
and integration of PV system to the power drive train or grid.
3. Model-based design of process and systems is very popular, there are tools
available, for automatic code generation for the developed model, it is required
to use these tools, that can deliver C code from the model, which can be used for
the cost-effective target hardware.
4. The memory of the real-time digital simulator (RTDS) is the main constraint
while simulating a complex model, like the grid, this memory is used for storing
and executing the compiled C code in real-time. It may be necessary to split the
model into two or more separate systems and bridge them using an appropriate
interface.
5. Interfacing the two models using respective interfacing algorithm introduces
some errors in the execution, that results in, instability of the system during
run time and also the accuracy of the results varies according to the interfacing
algorithm used.
The main objective of this work is to conduct the real-time simulation of the
PV considering as hardware under test (HUT), with the MPPT (Maximum Power
Tracking Point) system, storage and the load connected to DC bus acting as a rest of
the system (ROS) of hardware-in-the-loop (HIL) simulation methodology. The idea
is that the photovoltaic system and the rest of the system are executed in parallel in
real-time and the maximum power is monitored for irradiation variation in the real
world and the extracted power is supplied to the DC bus to which storage element
and the fixed load resistance are connected. The modeling of these systems is done
in Simulink [9]. The modeled systems are compiled and deployed in NI myRIO’s
FPGA [10] for real-time simulation.
The rest of this article is organized as follows: Sect. 2 deals with the modeling
of photovoltaic and MPPT systems, Sect. 3 explains the real-time simulation
architecture used in this work, Sect. 4 tells the real-time simulation with the results
and conclusions are drawn in Sect. 5.
Modelling of Photovoltaic Systems for Real-Time Hardware Simulation 5
2 PV System Modelling
The mathematical model of the PV system is drawn accordingly with [11, 12]. The
equivalent described in Fig. 1 consists of a current source in parallel with the diode
and Rsh and Rs are the intrinsic shunt and series resistances of the cell, respectively.
The I –V characteristic of the cell is determined by the diode. PV cells are grouped
into larger units to form PV modules, these modules are connected in series or
parallel to create a PV array to obtain required electrical power.
Iph = Isc + Ki (Tcell − Tref ) λ (1)
Isc
Irs = (2)
exp NcellqV oc
KnTcell − 1
3
Tcell qEg 1 1
Io = Irs exp − (3)
Tref nK Tref Tcell
KTcell
Vt = (4)
q
N
Vpv NcellpNs + Ipv Rs
Ish = (5)
Rsh
⎡ ⎛ Vpv Ipv Rs
⎞ ⎤
Ncell Ns + Np
Ipv = Np Iph − Np Io ⎣exp ⎝ ⎠ − 1⎦ − Ish (6)
nVt
where
Ipv Cell current
Vpv Array terminal voltage
Iph Photo current
Isc Short-circuit current at 25 ◦ C and 1 kW/m2
Ki Short-circuit current temperature coefficient
Tcell , Tref PV cell and reference temperature in K
450
400
350
PV array current [A]
300
250
200
150
100
Ir = 1kW/sq.m
Ir = 800W/sq.m
50 Ir = 500W/sq.m
Ir = 300W/sq.m
0
0 50 100 150 200 250 300 350
Voltage [V]
In this section the test benches created for conducting HIL experiment will be
presented. The mathematical model will be simulated in myRIO in real-time.
FPGA IO’s are used to integrate the model in-ports and out-ports to the real world
measurement systems (Fig. 4).
The model will be implemented in MATLAB-Simulink [9] and translated in into
C code with its internal C code generating tool using the NI VeriStand compilers
which adds on to the MATLAB and it can generate the file “.SO” which defines our
model with the libraries. NI VeriStand Engine has to be deployed in the RT target
to make it compatible to run with the host computer. .lvbitx file defines the custom
FPGA personality used for the simulation in the RT target. The main operation
is shown in Fig. 5. National Instruments myRIO-1900 is a Linux based device, it
8 H. Palahalli et al.
comes with built-in Xilinx Z-7010 FPGA [10] and it also allows us to build custom
FPGA applications with real-time capabilities. In this study, the real-time simulation
of PV system is done in myRIO using NI VeriStand [10], it is a configuration-based
testing software, specifically designed to run the real-time testing application, thus
allowing us to develop and test control system using hardware I/O and simulation
models.
This test bench consists of a 100 kW PV array model as HUT, whose interaction
with the ROS for the change in temperature and the irradiation is studied. The
ROS consists of MPPT, DC–DC converter, DC bus is created where the load
and the storage are connected with the breaker and the bi-directional intelligent
converter, respectively. The temperature and the irradiation data are measured using
a temperature sensor and a small 5 W PV cell, respectively. The model used to
simulate using this test bench is shown in Fig. 6.
The power rating and other parameters of the proposed RT simulated system
are given in Table 2, the voltage of DC bus is maintained at 250 V all the time
with the help of storage and the bi-directional DC-DC converter, which allows the
power flow in both the directions. PID control is used for the voltage control of the
converter. When the power required by the load is more than the generation, DC bus
voltage drops. Power will be compensated by the energy stored and the voltage is
brought regulated. When the power generated is more than the load requirement, the
DC bus voltage increases, and the controller takes action to bring back the voltage
to the set reference of 250 V. In this experiment, the fixed resistive load of 1 is
connected to the DC bus at 5th s, the outputs are driven to the real world using FPGA
IO and read using the digital oscilloscope, the test bench setup is shown in Fig. 7.
The irradiation and the temperature measured in real-time are given in Fig. 8.
They both are measured by myRIO and fed to the model during run time.
The load is connected to the DC bus at 5th s, the DC bus voltage and its variations
during run time is given in Fig. 9. It can be observed that the controller is able to
maintain the bus voltage constant even after connecting load at 5th s.
To ensure the simulation running in real-time the bus voltage is also measured
using a digital oscilloscope as shown in Fig. 10. The real-time simulation started
at −9 s, this corresponds to 0th s of the model time, the load is configured to get
connected to the DC bus at 5th s of the model time, and in the oscilloscope the
voltage variation due to load application is recorded at −4th s, this proves that the
simulation is running in real-time, just the voltage is scaled down by 100 V as 1 V
to drive output through DAC, finally the simulation is stopped at 1.6 s as seen in
oscilloscope.
The maximum power delivered by PV and the power consumed by the load after
connecting it to the DC bus is shown in Fig. 11. The power delivered by PV before
connecting the load is used to charge the storage system of 10 Ah. The state of
10 H. Palahalli et al.
MPPT
L1 Boost Converter
+ DC bus
–
Vpv Ipv
pulse
C1
Rload
100kW PV array
Bi-directional DC-DC Converter
Storage Rint L2
C2
ic
Not
gate
Vref
+ ∫ kp +
-
kd – Pulse
Generator
ki –
Vbus
charge measured in the storage system and the current flowing to the batteries are
given in Fig. 12. Initial SOC is about 50% and raised rapidly till connecting the load
to the DC bus, later the slope of SOC decreases since the flow of current is decreased
to the storage system. The voltage seen across the storage system throughout the
simulation is given in Fig. 13.
Modelling of Photovoltaic Systems for Real-Time Hardware Simulation 11
798
797
796
0 2 4 6 8 10 12 14 16 18
time (s)
24.945
Tempperatue (C)
24.94
24.935
24.93
24.925
0 2 4 6 8 10 12 14 16 18
time (s)
300
250
200
Voltage (V)
150
100
50
0
0 2 4 6 8 10 12 14 16 18
time (s)
5 Conclusion
104
8
6
Power (W)
4
0
0 2 4 6 8 10 12 14 16 18
time (s)
104
8
6
Power (W)
0
0 2 4 6 8 10 12 14 16 18
time (s)
55
SOC (%)
50
45
0 2 4 6 8 10 12 14 16
time (s)
400
200
Current (A)
-200
-400
0 2 4 6 8 10 12 14 16
time (s)
250
200
150
Voltage (V)
100
50
0
0 2 4 6 8 10 12 14 16 18
time (s)
References
1. Y. Huo, G. Gruosso, L. Piegari, Power hardware in the loop simulator of photovoltaic plant
for smart grid interation analysis, in 2017 IEEE International Conference on Environment
and Electrical Engineering and 2017 IEEE Industrial and Commercial Power Systems Europe
(EEEIC / I CPS Europe), June, pp. 1–5
2. M.D. Omar Faruque, T. Strasser, G. Lauss, V. Jalili-Marandi, P. Forsyth, C. Dufour, V.
Dinavahi, A. Monti, P. Kotsampopoulos, J.A. Martinez et al. Real-time simulation technologies
for power systems design, testing, and analysis. IEEE Power Energy Technol. Syst. J. 2(2), 63–
73 (2015)
3. A. Ingalalli, H. Satheesh, M. Kande, Platform for hardware in loop simulation, in 2016
International Symposium on Power Electronics, Electrical Drives, Automation and Motion
(SPEEDAM) (IEEE, New York, 2016), pp. 41–46
4. W. Li, G. Joós, J. Bélanger, Real-time simulation of a wind turbine generator coupled with
a battery supercapacitor energy storage system. IEEE Trans. Ind. Electron. 57(4), 1137–1145
(2010)
5. P.M. Menghal, A. Jaya Laxmi, Real time control of electrical machine drives: a review, in
2010 International Conference on Power, Control and Embedded Systems (ICPCES) (IEEE,
New York, 2010), pp. 1–6
6. B.L. Schenkman, D.G. Wilson, R.D. Robinett, Photovoltaic distributed generation for lanai
power grid real-time simulation and control integration scenario, in 2010 International
Symposium on Power Electronics Electrical Drives Automation and Motion (SPEEDAM)
(IEEE, New York, 2010), pp. 154–157
7. M. Dagbagi, L. Idkhajine, E. Monmasson, I. Slama-Belkhodja, FPGA implementation of
power electronic converter real-time model, in 2012 International Symposium on Power
Electronics, Electrical Drives, Automation and Motion (SPEEDAM) (IEEE, New York, 2012),
pp. 658–663
Modelling of Photovoltaic Systems for Real-Time Hardware Simulation 15
8. J.V. Paatero, P.D. Lund, Effects of large-scale photovoltaic power integration on electricity
distribution networks. Renew. Energy 32(2), 216–234 (2007)
9. The Mathworks. http://www.mathworks.com/
10. National Instruments. http://www.ni.com
11. F. Bizzarri, A. Brambilla, G. Gruosso, C. Guardiani, A. Sangiovanni Vincentelli, G. Storti
Gajani, Modeling and estimating yield and efficiency of photovoltaic solar parks, in 2013
IEEE International Conference on Industrial Technology (ICIT) (2013), pp. 734–739
12. F. Bizzarri, M. Bongiorno, A. Brambilla, G. Gruosso, G.S. Gajani, Model of photovoltaic
power plants for performance analysis and production forecast. IEEE Trans. Sustain. Energy
4(2), 278–285 (2013)
13. S. Sumathi, L. Ashok Kumar, P. Surekha, Application of MATLAB/SIMULINK in solar PV
systems, in Solar PV and Wind Energy Conversion Systems (Springer, New York, 2015),
pp. 59–143
14. M. Pisano, F. Bizzarri, A. Brambilla, G. Gruosso, G.S. Gajani, Micro-inverter for solar power
generation, in International Symposium on Power Electronics Power Electronics, Electrical
Drives, Automation and Motion, June (2012), pp. 109–113
Hardware-in-the-Loop Simulation
of High-Power Modular Converters
and Drives
Abstract This paper explains several industrial cases involving the HIL simulation
of MW-range drives and inverters using CPU cores with FPGAs to compute model
equations. The use of HIL simulators is common today in the industry to accelerate
design cycles, mitigate financial and human risks and support software updates
throughout the product life cycle.
The first case presented is a 2-level inverter scheme in which increasing power
specifications are met by adding parallel IGBT-modules. The second case is a multi-
level motor drive with low harmonic injection on the AC-side. The third case is a
modular multi-level converter in a grid application. We also discuss a new T-type
inverter model that uses an industry PV-to-grid power converter.
In each case, all power system modelling was done using Simulink and Sim-
PowerSystems in conjunction with the SSN solver from the ARTEMiS blockset in
addition to code generation for CPU execution at time steps in the 20–50 μm range,
with an exception for MMC models on FPGA. In all cases the firing accuracy of
the IGBTs remains in the nanosecond range using time-stamping techniques and
an FPGA board. In the case of the parallel 2-level inverters, there is significant
difficulty regarding the small firing delays (typically <500 ns) between modules
that create circulating currents. These circulating currents are rendered correctly on
the HIL bench.
Also discussed in the paper are the various optimisations, solvers and methods
that enable these performances.
C. Dufour ()
Opal-RT Technologies, Montréal, QC, Canada
e-mail: [email protected]
K. Palaniappan · B. J. Seibel
Rockwell Automation Inc., Milwaukee, WI, USA
e-mail: [email protected]; [email protected]
1 Introduction
The main objective of the HIL test system is to verify the functionality of the real
system in controlled laboratory conditions and at low power levels, before actual
release, including software upgrades and commissioning.
In particular, the HIL test system must:
– Check the drive start-up and shutdown sequences. Models must be able to
output truly null currents before the starting sequence for this to be checked
correctly in HIL mode. Otherwise an error occurs and the controller goes into safe
mode.
– Check the circulating currents between IGBT inverter modules in parallel
configurations. Small circulating currents are normal in the real system and
are due to variations in firing caused by wiring and element tolerances. The
drive closely monitors these currents and shuts them down if they exceed a pre-
determined level, which is an indication of system malfunction.
– The HIL system must be able to adequately reproduce the PWM inverter
characteristics on a CPU-based simulator running the model within a 25–50 μs
range. This is achieved using the time-stamping technique [1] and special inverter
models called Time-Stamped Bridges or TSB. Direct connection of controller
PWM pulses to the HIL simulator is mandatory (i.e., use of averaged models is
not possible).
– These TSBs must be able to work correctly in natural rectifying mode; this is
especially important for drives with Active-Front-End rectifiers.
Hardware-in-the-Loop Simulation of High-Power Modular Converters and Drives 19
This topology comes from the Rockwell PowerFlex 750-Series products for the
Low-Voltage market (from 160 to 6000 kW) and is depicted in Fig. 1 This topology
provides a scalable power level by using parallel 2-level IGBT modules up to 6 MW.
More details can be found in [2].1
An OPAL-RT Hardware-In-the-loop (HIL) system was chosen to perform a
wide variety of product software and hardware verification and validation during
the product design phase and will be used for regression testing over the life of
the product. As with any simulation, fidelity and accuracy of the simulation must
adequately match the product itself to guarantee usefulness and confidence in any
testing and verification. This section shows the PowerFlex 755TM Common Bus
Inverter. The inverter PWM frequency is 1.33 kHz and the complete model runs at
50 μs in HIL mode. In HIL tests, close to 2000 I/Os were required in the simulator
to interface and test the actual controllers in the real-time simulated drive system,
with up to ten converters at any given time during operation.
IGBT modules are connected together through interphase reactance to smooth
out any small difference of voltage output between the parallel inverters.
The drive controller closely monitors this current difference between the same
phases of parallel inverters and puts the system in fault mode if it rises above a
certain threshold.
1 The authors acknowledge the re-use of some of their own verbatim of this reference.
20 C. Dufour et al.
This inverter topology poses several difficulties, such as the accurate computation
of inter module currents and the need for accurate current levels at low power.
Furthermore, the model must be decoupled in small parts, running on different cores
of the CPU to maintain the simulation time step within an acceptable range.
The system has a large DC-link, which makes the separation of the model into
separate parts easier. DC-bus voltages and currents can be transmitted between
separated parts with an inserted delay, with negligible errors.
The model uses modified interpolating inverter models (Time-Stamped Bridge,
TSB); these inverter models are able to incorporate I/O gating signal time-stamps
and to interpolate voltage outputs, resulting in very accurate simulation, mainly
limited by the simulator FPGA sampling rate (in this instance, 5 ns).
Of special interest is the treatment of cases in which dead-time occurs near zero-
crossing. In this case, the standard TSB model may actually output opposite voltages
on different modules of the same phase.
For example, for a DC voltage of 640 V and an interphase inductance of 25 μH,
the circulating current rises at the rate of 640/(2 × 25e−6) × τ (for two parallel
IGBT inverters), where τ is the firing delay between parallel IGBTs. For a 500 ns
delay, the current would rise to about 7 A, within the controller’s tolerance level.
As a simple test, we ran the model offline at 50 μs, in a simple motor start-up
sequence with two parallel IGBT modules, with a high-resolution PWM modulator
and variable inserted delay between the parallel IGBTs. The result is depicted in
Fig. 2 and is consistent with the theory.
The requirement for null current during blocked modes, as well as accurate
natural rectification modes, led us to model the diode as binary switches. This in
turn caused the total number of connected diodes to be very high (60 diodes for
ten parallel modules), dramatically increasing the total number of possible electric
modes and permutations of state-space equations (possibly 260 ). By using SSN, each
TSB module, containing six binary diodes each, was set into a different SSN group,
solving this potential issue for real-time calculation.
Finally, the induction motor model used was an SSN, meaning that no minimum
load or stabilisation snubber was required to run the model stably, even at very low
power levels [3].
The multi-level motor drive depicted in Fig. 3 is designed to provide low harmonic
AC-voltage to the load, as well as low input current harmonic to the feeding grid.
This is achieved by arranging the multiple diode rectifiers with zig-zag transformers
with varying phase shifts. The system runs in HIL with attached industrial controller
under test at a time step of 25 μs.
Hardware-in-the-Loop Simulation of High-Power Modular Converters and Drives 21
Fig. 2 Motor start-up with variable firing delay between two IGBT modules. Top: motor current;
bottom: current difference between the parallel inverters
This inverter topology presents several difficulties for HIL simulation. One of these
is the common inductive connection of all zig-zag transformer primary windings;
the SSN solver makes it possible to break down this huge state-space system of
equations into several smaller ones, without delays. In one version of the model,
where the 3 × n 6-pulse rectifier (here n is the number of inverter stages) was
separated using SSN nodes, SSN was also required to decouple the 3 × n 6
pulse rectifiers. In the final implementation, stublines were used to provide a better
decoupling at the expense of slightly lower precision.
Also, the presence of a DC-bus at each 6-pulse rectifier allowed us to simulate
the rectifiers and the inverters in different cores of the simulator.
The SSN node location allows for the creation of three large SSN groups, one for
each phase of zig-zag transformers, which allowed us to compute these groups on
different cores.
Stublines are artificial one-time-step-delay transmission lines and are used to
completely decouple a circuit from one end of the stubline to the other. As they add
a little capacitance to the circuit, they are less precise than SSN, which does not
approximate or add delays to the simulated equations.
To achieve HIL simulation of this MMC system with 3240 IGBTs and 1620
capacitors, the modelling approach is different from the previous cases: all IGBT
switching and capacitor voltage calculations are done directly on the FPGA, and a
Hardware-in-the-Loop Simulation of High-Power Modular Converters and Drives 23
Thevenin equivalent is sent to the CPU for the SSN solver updates with iterative
surge arresters. This model runs in HIL at a time step of 25 μs [6].
Among recent developments is the creation of a TSB-type model for the T-type 3-
level inverter in ARTEMiS-SSN. The topology is tested in this section for a 25 kW
industrial PV converter for grid, depicted in Fig. 5.
The converter’s PWM switching frequency is 20 kHz. This is typical for PV
inverters rated below 30 kW for residential and commercial PV applications where
audible noise is not acceptable.
24 C. Dufour et al.
S1 400V
600 W/m2 S4 S3
50Hz
400μH
The model was simulated in a multi-rate offline model as follows: the grid and
converter were simulated at 15 μs while the controller and PWM modulator were
simulated at 1 μs and interfaced with the electric system using time-stamping.
This time-stamping approach very closely imitates the HIL time-stamping method,
described in Sect. 4.1.
This test case is one where the converter delivers 25 kW of power to the grid.
Then at 0.25 s, the DC-link voltage reference is changed from 800 to 850 V. This
results in a small gap in the currents; while the solar panel charges up the DC link, it
does not deliver power to the grid; then at 0.26 s, the DC-link meets the commanded
values and power flow to the grid is restored. The simulation result is compared with
a reference model in which the entire model is run at 1 μs, depicted in Fig. 6.
The figure is zoomed in on the currents to show that the PWM component is
still present in the simulation but somehow attenuated. This is normal considering
the low ratio of simulation sampling frequency to PWM frequency. For Ts = 15 μs
and FPWM = 20 kHz, this ratio is only 3.33! Nevertheless, the current amplitudes
are very accurate and free of amplitude jitter, which is sufficient to fully test the
controller.
This also clears up some common confusion about TSBs: they are NOT averaged
models, because the PWM component is included in the simulation. Rather, TSBs
could be called ‘Per time step averaging models’, while standard averaged models
are, by definition, ‘Per fundamental period averaging models’.
4 Miscellaneous Topics
Ts=15 us
50
Currents (A)
-50
Ts=1 us
50
Currents (A)
-50
0.2 0.21 0.22 0.23 0.24 0.25 0.26 0.27 0.28 0.29 0.3
Time (s)
HIL simulator, as shown in Fig. 7 [7, 8]. A similar approach is done without I/Os by
calculating the intra-step zero-crossing of PWM modulators.
In certain cases, such as MMC, some parts of the model can be computed on the
FPGA.
All the machine models used in these models were SSN machines; that is, machine
models developed using the nodal admittance method of SSN [3, 9, 10]. The
main reason for this choice was to ensure stability without any parasitic loads or
stabilisation snubbers at the machine terminals.
This was especially important for the Two-Level Parallel IGBT Modular Motor
Drive because of the inductive connection of the motor and the interface inductor
of the different modules. Standard SPS machine models use the current injection
technique, which is prone to numerical instability, and is actually unstable when
injected into such inductive node.
26 C. Dufour et al.
FPGA chips can be powerful enough to compute drive equations in HIL simulators.
Special advantages of FGPA include:
• Smaller time steps, which are typically under 1 μs.
• Lower input–output latency in HIL systems (considering that I/Os and models
are together on the FPGA).
• A natural way to sample high-frequency PWM drives.
Recent advances in the field of electric circuit simulation on FPGA [11, 12]
now enable the simulation of variable topology, variable parameter circuits, all
without the need for bitstream generation, which until recently has been a powerful
disincentive in FPGA technology since it can take hours to place and route (the
equivalent of compiling with a CPU). Today, FPGAs are still considered as a
specialist’s domain with several limitations, in particular the coding of complex
and flexible algorithms such as SSN. This is the main reason that CPU-based HIL
simulation is still preferred today.
The best approach here is to use both technologies in a pragmatic approach:
unless very low-latency (<5 μs) or PWM components are absolutely required in the
simulation (e.g., resonant converters), CPUs are preferable. When the model is run
on a CPU, there is a typical two-time step latency involved, caused by data transfer
to and from I/Os at each time step.
The case of the MMC is particularly interesting in its pragmatism: the MMC
topology is indeed better suited to FPGA simulation because of its extremely
high modularity and number of devices. FPGA makes it possible to simulate this
topology efficiently, using pipelining methods. It also allows for direct connectivity
of thousands of I/O points. The more complex model parts of the complete MMC
system, such as iterative surge arresters, are then computed on the simulator’s CPU
using the SSN solver.
28 C. Dufour et al.
It is worth noting that high PWM frequency is NOT always a good reason to run
a model on FPGA. With the time-stamping techniques described in this paper,
accuracy is excellent even at large simulation time steps. Even when the sampling
frequency is close to the PWM frequency, this technique still provides good
accuracy by filtering out the PWM component in the simulation. It happens that
most converters are designed so that PWM components are filtered out (with some
exceptions such as resonant converters); therefore this PWM filtering effect often
has no impact on tests. Also, the time-stamping method has an obvious advantage
over average models: in HIL mode, power electronic devices can be directly driven
by the controller firing pulses.
In this paper, we have shown a good example using a 3-level T-type inverter in
which the model sampling frequency was only 3.33 times higher than the PWM
frequency where accuracy remains high. PWM modulation is simply a way to
convert DC to AC and the PWM component is not normally tested (again, there
are exceptions, such as resonant converters). So, for control testing purposes, it is
often acceptable to have this PWM component damped.
The inverter used in the models of this paper are of a novel generation
(called TSB-RD), a kind of hybrid between interpolated switching-functions, for
the active mode (similar to older generation of TSB) and binary switches (in
SimPowerSystems, these are 2-state switches with very low resistance if ON or
very high resistance if OFF) for the natural rectification modes. It is worth noting
that older generations of TSB, in which the rectifying modes were simulated using
current-nulling feedback loops, can still be used because they are sometimes more
computationally efficient than using TSB-RD with SSN. That was actually the case
in the Multi-level Motor Drive with Low Input Harmonic (Sect. 3.2).
5 Conclusions
This paper discussed the different modelling techniques used to simulate several
large motor drives and converters. All these topologies have been successfully
commissioned at industrial client sites.
Here the motto of HIL simulation is ‘Make it work!’, and this motto is
additionally the reason for providing all these various models, options, decoupling
techniques and solvers in ARTEMiS: to meet all the demanding client objectives.
References
1. P. Terwiesch, T. Keller, E. Scheiben, Rail vehicle control system integration testing using digital
hardware-in-the-loop simulation. IEEE Trans. Cont. Syst. Technol. 7(3), 352–362 (1999)
2. K. Palaniappan, B.J. Seibel, C. Dufour, Real time hardware in the loop validation of common
bus inverter low voltage drives, in IEEE Applied Power Electronics Conference & Exposition
(APEC-2019), Anaheim, California, 17–21 Mar 2019
Hardware-in-the-Loop Simulation of High-Power Modular Converters and Drives 29
3. C. Dufour, Highly stable rotating machine models using the state-space-nodal real-time solver,
in COMPENG-2018 conference, Florence, Italia, 10–12 Oct 2018
4. G. Li, Y. Dong, J. Tian, W. Wang, W. Li, J. Bélanger, Factory acceptance test of a five-terminal
MMC control and protection system using hardware-in-the-loop method, in Proc. of 2015 IEEE
Power & Energy Society General Meeting, Denver, USA, 26–30 July 2015
5. Z. Zhu, X. Li, H. Rao, W. Wang, W. Li, Testing a complete control and protection system for
multi-terminal MMC HVDC links using hardware-in-the-loop simulation, in IECON 2014 -
40th Annual Conference of the IEEE Industrial Electronics Society, pp. 4402–4408
6. C. Dufour, W. Li, X. Xiao, J.-N. Paquin, J. Bélanger, Fault studies of MMC-HVDC links using
FPGA and CPU on a real-time simulator with iteration capability, in 11th International Con-
ference on Compatibility, Power Electronics and Power Engineering (IEEE CPE-POWERENG
2017, Cadiz, Spain, 4–6 Apr 2017
7. C. Dufour, S. Abourida, J. Belanger, Real-time simulation of electrical vehicle motor drives on
a PC cluster, in 10th European Conference on Power Electronics EPE, Sep 2003
8. M. Harakawa, H. Yamasaki, T. Nagano, S. Abourida, C. Dufour, J. Bélanger, Real-time
simulation of a complete PMSM drive at 10 μs time step, in Proceedings of the 2005
International Power Electronics Conference - Niigata (IPEC-Niigata 2005), 2005
9. C. Dufour, J. Mahseredjian, J. Bélanger, A combined state-space nodal method for the
simulation of power system transients. IEEE Trans. Power Deliv. 26(2), 928–935 (2011). (ISSN
0885–8977)
10. C. Dufour, J. Mahseredjian, J. Bélanger, J.L. Naredo, An advanced real-time electro-magnetic
simulator for power systems with a simultaneous state-space nodal solver, in IEEE/PES T&D
2010 - Latin America, São Paulo, Brazil, 8–10 Nov 2010
11. C. Dufour, S. Cense, J. Bélanger, An induction machine and power electronic test system
on a field-programmable gate array. Math. Comput. Simul., 112–123 (2016). https://doi.org/
10.1016/j.matcom.2016.03.014
12. T. Ould Bachir, C. Dufour, J. Bélanger, J. Mahseredjian, J.P. David, A fully automated
reconfigurable calculation engine dedicated to the real-time simulation of high switching
frequency power electronic circuits. Math. Comput. Simul. 91, 167–177 (2013)
Concept and Implementation of a Rotor
Position Detection Method for Permanent
Magnet Synchronous Machines Based
on Linear Halls
Yuchen Wang, Kai Liu, Wei Hua, Xiaofeng Zhu, and Baoan Wang
Abstract The purpose of this paper is to propose and implement a novel rotor
position detection method for permanent magnet synchronous machines (PMSMs)
based on linear Halls, which are embedded inside of stator of PMSMs. A three-
phase 9-slots/8-poles PMSM is exampled to verify the method. Firstly, a special
point located in stator yoke (back-iron) is found by two-dimensional finite element
analysis (2D-FEA), where the open-circuit flux-density due to permanent magnets
versus rotor position (BPM ) shows a high amplitude and good linearity, while the
armature-reaction flux-density (Barmature ) due to armature currents exhibits a low
amplitude and good linearity versus armature currents. Then, an analytical model
is built and the analytical relationship between armature currents and the Barmature
is derived. Based on the analytical mode, BPM can be obtained by separating the
Barmature from the synthetic magnetic field (BSynthetic ). Thereafter, the resultant BPM
can be used to detect the rotor position information with differential-type piecewise-
linear analytical method. The feasibility of the proposed detection method is verified
by co-simulations and experiments. The simulation results show that the novel
linear Hall-based angle sensor can achieve the accuracy equivalent to 3000-line. The
experimental results indicate that compared with an encoder, the maximum error of
electric angle position at different speeds is less than 0.3%.
1 Introduction
Ideally, the magnetic field detected by the linear Hall should be in the form of an
analytical expression versus rotor position, which can enhance the accuracy of angle
decoding. Firstly, two-dimensional finite element analysis (2D-FEA) is carried out
by Ansoft Maxwell to find suitable mounting locations for the linear Hall.
From the perspective of magnetic field amplitude, the amplitude of the BPM
should be as large as possible, while the amplitude of the Barmature should be as
small as possible. From the perspective of solving accuracy, both the BPM versus
rotor position and Barmature versus armature currents should be linear.
A 9-slot/8-pole PMSM is exampled as shown in Fig. 1. Considering the stator is
symmetrical in the circumference, nine special points are selected as magnetic field
detection points within a tooth slot range.
Secondly, produced separately by PMs and armature currents, both the tangential
and radial components of the flux-densities of the nine detection points are analysed.
Among them, the tangential component of the BPM versus rotor position is shown in
Concept and Implementation of a Rotor Position Detection Method for. . . 33
Fig. 1 Distributions of
particular points in a
9-slot/8-pole PMSM
0
0 60 120 180 240 300 360
-1
-2
-3
Angle(°)
Fig. 2. The simulation results are shown in Table 1. In summary, point 2 is selected
as the magnetic detection point due to its high linearity for both BPM and Barmature .
For convenience, ‘point D’ represents ‘point 2’ in the following.
Moreover, considering the geometrical periodic of the 9-slot/8-pole PMSM, there
are other eight ‘point Ds’ on the stator, i.e., totally nine ‘point Ds’ with the detailed
locations listed in Table 2. Moreover, the nine ‘point Ds’ are divided into three
groups, e.g., D1 , D4 , and D7 for Group 1. The three Barmature versus time waveforms
in same group are distributed symmetrically by 120◦ in electrical degrees, and the
34 Y. Wang et al.
BT = Bm sin (θ + ϕ)
= Bm (sin
θ cos ϕ + cos θ sin ϕ) (2)
= Bm I√b −Ic cos ϕ − IbI+Ic
sin ϕ
3Im m
= EI b + F I c + G
◦ ◦
where, E = − √2Bm sin ϕ − 30 , F = − √2Bm sin ϕ + 30 , and G represents the
3Im 3Im
error term.
Concept and Implementation of a Rotor Position Detection Method for. . . 35
0
0 10 20 30 40
-0.3
-0.6
Time(ms)
The Barmature with armature currents of point D1 can be fit by using the curve
fitting toolbox in Matlab and the following result is obtained:
⎧
⎨ Ua = Um sin (θ
)
+ UOQ + ΔU
◦
U = Um sin θ − 120 + UOQ + ΔU (4)
⎩ b ◦
Uc = Um sin θ + 120 + UOQ + ΔU
⎧ √ ◦
⎨ Uab = Ua − Ub = √3Um sin θ + 30
◦
U = Ub − Uc =√ 3Um sin θ − 90 (5)
⎩ bc ◦
Uca = Uc − Ua = 3Um sin θ + 150
where Ua , Ub , and Uc are signals from three linear Halls located at points D1 , D4 ,
and D7 , as shown in Fig. 4; Uab , Ubc , and Uca are signals after the approximate
differential process, as shown in Fig. 5.
As for Eq. (5), a six-quadrant symmetric partition is performed according to
Fig. 5. In each quadrant region, the signal with highest linearity (marked with red
line) will be utilized by linear analytical solving, which treats the waveform as a line
segment, accelerating the dynamic response.
This novel method eliminates the DC offset without √ adapting addition differ-
ential linear Hall and increases the signal amplitude to 3 times, improving the
Angle(°)
Angle(°)
Concept and Implementation of a Rotor Position Detection Method for. . . 37
The rotor position can be obtained, according to the flow chart (Fig. 6), where two
sets of raw data are needed in angle solving system
1. Three-phase armature currents are obtained by high-precision current sensor
detection.
2. The voltage signals of Hall group are derived from three linear Halls embedded
in the stator.
Then,
Three-phase
Voltage signal
armature
from linear Hall
currents
Voltage-magnetic field
conversion
Current-magnetic field
conversion
Synthetic magnetic
field signal
Armature magnetic
Subtraction
field signal
Permanent magnet
magnetic field
Differential piecewise
linear analytic method
Rotor Position
End
Other documents randomly have
different content
preceding side, and, therefore, had it been a star, must have passed
over it, whereas it preserved the same apparent distance all the
while from the edge of the ring. (N.B. Saturn not very far from the
zenith on merid.)
Next night, Friday, May 5, Saturn most gloriously seen: quite as
sharp as any copper-plate engraving, with power 240 and full
aperture. All the five old satellites seen and measured, being now on
the opposite side. Now considerably short of its greatest following
elongation; distance just as before; and, as on Thursday, it was kept
in view long enough for Saturn to have left it behind by its own
motion had it been a star. The change of situation agrees perfectly
with the period 1d 9h, which is also the reason why it was not seen
May 5th, being on that night near its inferior conjunction. So this is
at last a thing made out. As for No. Seven I have no hope of ever
seeing it.
Fig. 2.
My dear Aunt,—
* * * * *
I need hardly say how much we are rejoiced to see your
handwriting once more, though that joy is damped by your
complaints of winter indisposition. And such a winter! by all
accounts. May this prove a better! and may we hope to find you in
no worse health and spirits when we come to see you next summer
in Hanover. For so, if it please God to lead us safe home, according
to our present altered plans, we most assuredly propose to do.
I say our altered plans, for you know our intention was to have
embarked next March for Rio Janeiro, and there to have spent two
or three months, after which to have taken passage in the Brazilian
packet for England, which would have probably detained us till
October, and have rendered a visit to Hanover that season
impracticable. But by striking off this Brazilian trip, and taking our
course directly homewards, so much time will be saved, and all the
rest of our domestic arrangements become so much simplified that it
seems like finding a treasure, as a fund of time will thereby be
placed at our disposal, the first fruits of which, as in all love and duty
bound, we have determined to devote to you; or rather, I should say,
that, when in talking over with Margaret all the pro’s and con’s of the
question, whether to return home direct, or viâ Brazil?—this
consideration at once decided it in favour of the direct course, her
desire to see you outweighing every consideration of amusement or
temporary gratification which a visit to Rio could offer. So now be
sure, dear aunty, and keep yourself well, and let us find you in your
best looks and spirits; and, although what you say respecting our
good Mrs. Beckedorff’s health is somewhat deplorable, yet I will
indulge the hope that she too will perform a part in the dramatis
personæ of that happy meeting. Meanwhile, as the time of our
departure hence approaches, we shall take care and apprise you of
all our movements, respecting which it is impossible at present to
speak more precisely.
“Madam,—
“In transmitting to you the accompanying Diploma from the Royal
Irish Academy, I wish to be allowed to add, as I thus do, the
expression of my own high sense of your services to Astronomy, and
of the eminent degree in which you have deserved the present
testimonial.
“I have the honour to be, Madam, &c., &c.,
“William Rowan Hamilton,
“P.R.I.A.”
1839. Life in
Hanover.
MISS HERSCHEL TO SIR JOHN HERSCHEL.
Hanover, Dec. 1, 1839.
Dear Nephew,—
Do not you think I have been very good to let the most dismal
month in the year pass without troubling you for accounts of the
progress my dear niece is making in her recovery?
My dear niece said once, I should write often, and in few lines
inform her how I go on, so I must say—I get up as usual every day,
change my clothing, eat, drink, and go to sleep again on the sofa,
except I am roused by visitors; then I talk till I can no more—
nineteen to the dozen! N.B. I don’t tell fibs, though they may not
always like what I say.
I have been twice at the concert, and each time been honoured
with a wie gehts?[51] by His Majesty, and the notice of many
acquaintances whom I have no opportunity of seeing elsewhere, the
public concerts being the only place where I can go with the least
trouble to myself or others. You say when I talk of the Gelehrten
then all goes well, but I know nothing about them....
But one piece of news I must tell you, which is, that a fortnight
after Dr. Mädler had been the conductor of Mde. Witte (the Moon)
and her daughter to the meeting at Pyrmont, I received two cards,
the one, “Professor Dr. Mädler,” under it, “Minna Witte-Verlobt.[52]”
The reason Madame Witte gives for this hasty courtship is, that it is
Dr. M.’s first love, and that he would not wait, so the lady said yes!
As you have seen this lady, I would give you this piece of news.
* * * * *
I beg you will give a true account of my dear niece’s and the
children’s health, not forgetting the babe and how she will be
named, that I may enter the same in my biographical account.
I remain, my dear nephew,
Your most affectionate aunt,
Car. Herschel.
Aug. 3rd.—I went to buy some clothing for wearing at home, and
went to my mantua-maker to give directions. I had to climb up to
the third story, and I was of course quite knocked up when I came
home, but it is my intention to continue to take some exercise as
long as the weather and the length of the afternoon will permit.
1839-1840. Her Aug. 26th.—My niece Knipping came in
Day-Book. the afternoon to assist me in some
needlework—we did not do much!
Sept. 25th.—To-day I was made happy by a visit of Alexander
Humboldt; which, though it was extended to the utmost limit of the
time which this interesting man could spare me, was too short for all
I wished to hear and had to say, which, as the theme of our
conversation was my nephew, may be easily imagined.
Oct. 5th & 6th.—Mr. Hohenbaum and the carpenter were with me
to pack up the seven-foot telescope. I assisted as well as I could,
being very ill all the while.
Oct. 7th.—Dr. G. called for a moment, but nobody else!
Dec. 10th.—I went in the evening to the concert, where I exposed
myself most sadly by falling a-crying when the King most kindly
came to me to inquire after my health. I do not think I shall have
the courage to show myself there again in a hurry.
Jan. 27th.—This is the first day since the 30th December that the
ice is detached from my sitting-room window.
Jan. 31st.—Mr. Hausmann brought me some Journals, and talked
for an hour of old times with me, as he ever does, good man!
Feb. 7th.—A letter from my niece came this morning by the
Hamburger post, which will make me happy for some time, and
make me bear my painful solitude more patiently.
March 17th.—Thank God the 7th and 16th March are got over, and
I begin to recollect that I have much else to do than bewail myself at
being obliged to spend such days severed from all that are, or were,
so dear!... I found my poor friend [Mrs. Beckedorff] very much
altered, but before I left her I thought she looked a twelvemonth
younger for our two hours’ chat. But we both were obliged to part,
for we could no more. Yesterday she sent me some fine flowers, as
usual on my birthday. Dr. Mühry left a card; two of my nieces called,
and Hofräthin Ubelode brought me some flowers. They left me
fatigued to death, to spend the long evening in solitude.
June 18th.—Yesterday Mr. Hausmann came to see me, and
brought the Philosophical Magazine for June, in which I had the
pleasure to see that Dr. Lamont has observed three of the Georgium
Sidus satellites.
July 3rd.—Dr. G. brought me an extract from The Sun that my
nephew has been created a baronet on the occasion of the
coronation.
July 9th.—My nephew arrived in Hanover in the evening.
July 10th.—In the afternoon I saw him and my little grand-nephew
for a few hours.
July 25th.—My nephew and his son took tea with me, and we
soon parted, without taking leave, and next morning I am told they
left Hanover at four in the morning. More I cannot say!
Oct. 24th.—Mr. Hausmann came in the forenoon and took the box
with the mirror of my sweeper with him, and in the evening he came
to receive the stand. I am glad my poor sweeper is now in good
hands!
Oct. 29th.—Mrs. Knipping spent an hour with me in the dusk of
the evening, and read an act of a play.
Dec. 30th.—In the afternoon Fraulein S. came to see me, but she
is deaf. I talked with her for a couple of hours without either of us
being the wiser.
Jan. 5th.—Went in the evening to the concert; had some talk with
the Levies, who delighted the company with their performance,
especially the youngest son, eight years of age, who gave several
pieces on the French horn. Conversed with several persons besides
the Prince Solms.
Jan. 20th.—I have been to the concert last night to hear the
wonderful violinist, Ole Bull. It was very crowded for the confined
room, though the largest in Hanover next the play-house. By the
help of Miss B. and the M.’s I got safely through the crowd to my
chair. But I was somewhat disappointed, for, by the report of those
who had heard Ole Bull before, I expected to hear a virtuoso on the
violin who would have given us an idea of the manner of
performance of a Jordine, Kramer, Jacob Herschel, and Dietrich too;
but it is more like conjuration than playing on a violin.
Feb. 12th.—Dr. Lissing paid me a visit. He wished me to subscribe
to a work on Magnetism, but I think it would look only like
affectation to let my name appear among the learned subscribers on
a subject of which I know so little.
March 16th.—Mrs. Beckedorff sent me two beautiful flowers,
accompanied by her good wishes, which she never forgets to do on
my birthday. Mde. Knipping, and others, came to wish me to live
many more years,—but what can I say?
March 23rd.—I was at the last subscription concert. His Majesty
was there, and asked me how I did? I said, tolerably! This was all
our conversation.
July 16th.—The whole of yesterday I had no other prospect but
that it would have been the last of the days of sorrow, trouble, and
disappointment I have spent from the moment I had any recollection
of my existence, which is from between my third and fourth year....
In the night I fell out of one fainting fit into another, and when I
came to my recollection, between six and seven in the morning, I
found Dr. G. sitting before me talking loud in his usual nonsensical
way. Him had Betty called in her fright, for his wife (who is of use to
nobody) is gone to spend the summer months in the country. Mde.
Knipping also is away.
1839. End of Day- July 25th.—Mr. Hausmann, junior, and Mr.
Book. Hohenbaum called to look at the
photographical drawing. I am told it is the only specimen of the kind
in Hanover.
This Day-book, No. 2, is now full, and I shall not be easy till it is
deposited in a portfolio, in which will also be found the Mem.-book
9.... It often enables me to contradict erroneous impertinent notions
concerning my brother William’s disinterested character.
I am now not able even to look over, much less to correct, what I
have scribbled, but it must go as it is. Perhaps my dear niece may
look into them at some leisure moment, and she will see what a
solitary and useless life I have led these seventeen years, all owing
to not finding Hanover, nor anyone in it, like what I left, when the
best of brothers took me with him to England in August, 1772!
1840. Anecdote of
the Old Telescope.
MISS HERSCHEL TO LADY HERSCHEL.
Jan. 10, 1840.
My dearest Niece,—
* * * * *
Perhaps you may have heard that in the early part of its [the forty-
foot telescope’s] existence, “God save the King” was sung in it by
the whole company, who got up from dinner and went into the tube,
among the rest two Misses Stows, the one a famous pianoforte
player, some of the Griesbachs, who accompanied on the oboe, or
any instrument they could get hold of, and I, you will easily imagine,
was one of the nimblest and foremost to get in and out of the tube.
But now!—lack-a-day!—I can hardly cross the room without help.
But what of that? Dorcas, in the Beggar’s Opera, says, “One cannot
eat one’s cake and have it too!”... I will only thank you once more for
your charming letter, and beg to be kindly remembered to all who
are dear to you, and to give an embrace extraordinary to the dear
little ones around you, and not forgetting to include my dear nephew
in the general hug! and believe me,
My dearest niece,
Yours and his most affectionate aunt,
Car. Herschel.
P.S.—One anecdote of the old tube (if you have not heard it) I
must give you. Before the optical parts were finished, many visitors
had the curiosity to walk through it, among the rest King George III.,
and the Archbishop of Canterbury, following the King, and finding it
difficult to proceed, the King turned to give him the hand, saying,
“Come, my Lord Bishop, I will show you the way to Heaven!”
This was in the year 1787, August 17th, when the King and
Queen, the Duke of York and some of the Princesses were of the
company.
I hope the book where the visitors were noted, has been
preserved? Some time after it was kept by other hands; but before I
parted with it, I copied some pages which put me sometimes in
mind of persons who were interesting to me.
These scribblings will come to you among the rest of my scraps.
Good-bye!
ebookball.com