100% found this document useful (5 votes)
23 views

Fundamentals of RCS Prediction Methodology using Parallelized Numerical Electromagnetics Code NEC and Finite Element Pre processor Vineetha Joy pdf download

The document outlines a methodology for predicting radar cross section (RCS) using a parallelized version of the Numerical Electromagnetics Code (NEC) combined with a finite element pre-processor. It discusses the theoretical background, formulation, and modeling guidelines necessary for accurate RCS computation, particularly in the context of stealth technology. The work serves as a resource for researchers and engineers in the field of computational electromagnetics.

Uploaded by

zankovisco0u
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
100% found this document useful (5 votes)
23 views

Fundamentals of RCS Prediction Methodology using Parallelized Numerical Electromagnetics Code NEC and Finite Element Pre processor Vineetha Joy pdf download

The document outlines a methodology for predicting radar cross section (RCS) using a parallelized version of the Numerical Electromagnetics Code (NEC) combined with a finite element pre-processor. It discusses the theoretical background, formulation, and modeling guidelines necessary for accurate RCS computation, particularly in the context of stealth technology. The work serves as a resource for researchers and engineers in the field of computational electromagnetics.

Uploaded by

zankovisco0u
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 54

Fundamentals of RCS Prediction Methodology using

Parallelized Numerical Electromagnetics Code NEC


and Finite Element Pre processor Vineetha Joy
download
https://textbookfull.com/product/fundamentals-of-rcs-prediction-
methodology-using-parallelized-numerical-electromagnetics-code-
nec-and-finite-element-pre-processor-vineetha-joy/

Download more ebook from https://textbookfull.com


We believe these products will be a great fit for you. Click
the link to download now, or visit textbookfull.com
to discover even more!

Electromagnetics through the Finite Element Method: A


Simplified Approach Using Maxwell's Equations 1st
Edition José Roberto Cardoso

https://textbookfull.com/product/electromagnetics-through-the-
finite-element-method-a-simplified-approach-using-maxwells-
equations-1st-edition-jose-roberto-cardoso/

Numerical Solution of Differential Equations:


Introduction to Finite Difference and Finite Element
Methods 1st Edition Zhilin Li

https://textbookfull.com/product/numerical-solution-of-
differential-equations-introduction-to-finite-difference-and-
finite-element-methods-1st-edition-zhilin-li/

Harmonic Balance Finite Element Method: Applications in


Nonlinear Electromagnetics and Power Systems 1st
Edition Junwei Lu

https://textbookfull.com/product/harmonic-balance-finite-element-
method-applications-in-nonlinear-electromagnetics-and-power-
systems-1st-edition-junwei-lu/

Fundamentals of the Finite Element Method for Heat and


Mass Transfer 2nd Edition Perumal Nithiarasu

https://textbookfull.com/product/fundamentals-of-the-finite-
element-method-for-heat-and-mass-transfer-2nd-edition-perumal-
nithiarasu/
Finite Element Analysis of Weld Thermal Cycles Using
ANSYS 1st Edition G. Ravichandran (Author)

https://textbookfull.com/product/finite-element-analysis-of-weld-
thermal-cycles-using-ansys-1st-edition-g-ravichandran-author/

The finite element method: fundamentals and


applications in civil, hydraulic, mechanical and
aeronautical engineering Zhu B.

https://textbookfull.com/product/the-finite-element-method-
fundamentals-and-applications-in-civil-hydraulic-mechanical-and-
aeronautical-engineering-zhu-b/

Extended finite element and meshfree methods Timon


Rabczuk

https://textbookfull.com/product/extended-finite-element-and-
meshfree-methods-timon-rabczuk/

Finite Element Analysis A Primer 2nd Edition Musa

https://textbookfull.com/product/finite-element-analysis-a-
primer-2nd-edition-musa/

Numerical Methods for Engineering: An introduction


using MATLAB® and computational electromagnetics
examples (Electromagnetic Waves) 2nd Edition Karl F.
Warnick
https://textbookfull.com/product/numerical-methods-for-
engineering-an-introduction-using-matlab-and-computational-
electromagnetics-examples-electromagnetic-waves-2nd-edition-karl-
SPRINGER BRIEFS IN ELEC TRIC AL AND COMPUTER
ENGINEERING  COMPUTATIONAL ELECTROMAGNETICS

Vineetha Joy
G. L. Rajeshwari
Hema Singh
Raveendranath U. Nair

Fundamentals
of RCS Prediction
Methodology using
Parallelized Numerical
Electromagnetics
Code (NEC) and Finite
Element Pre-processor

123
SpringerBriefs in Electrical and Computer
Engineering

SpringerBriefs in Computational Electromagnetics

Series Editors
K. J. Vinoy, Indian Institute of Science, Professor, Electrical Communication
Engineering, Indian Institute of Science, Bangalore, Karnataka, India
Rakesh Mohan Jha (Late), Centre for Electromagnetics, CSIR-National Aerospace
Laboratories, Bangalore, Karnataka, India
The phenomenon of electric and magnetic field vector (wave) propagation through the
free-space, or any other medium is considered within the ambit of electromagnetics. The
media themselves, in general, could be of diverse type, such as linear/non-linear,
isotropic/non-isotropic, homogeneous/inhomogeneous, reciprocal/non-reciprocal, etc. Such
electromagnetic wave propagation problems are formulated with the set of Maxwell’s
equations. Computational Electromagnetics endeavors to provide the solution to the
Maxwell’s equations for a given formulation. It is often difficult to find closed forms
solutions to the Maxwell’s equation formulations. The advent of computers, and particularly
the initial developments of efficient coding for numerical analysis, encouraged the
development of numerical electromagnetics. A second motivation came from the interaction
of the electromagnetic wave with the matter. This could be visualized as scattering bodies,
which required incorporation of the phenomena of reflection, refraction, diffraction and
polarization. The finite/large nature of the scatterer required that problem of electromagnetics
is considered with respect to the operation wavelength leading to the classification of
low-frequency, high-frequency and resonance region problems. This also inspired various
asymptotic and grid-based finite-method techniques, for solving specific electromagnetic
problems. Surface modeling and ray tracing are also considered for such electromagnetic
problems. Further, design optimization towards hardware realization have led to the recourse
to various soft computing algorithms. Computational Electromagnetics is deemed to
encompass the numerical electromagnetics along with all other above developments. With
the wide availability of massively parallel high performance parallel computing platforms,
new possibilities have emerged for reducing the computation time and developing macro
models that can even be employed for several practical multi-physics scenarios. Both volume
and surface discretization methods have been given a new boost, and several acceleration
techniques including GPU based computation, learning based approaches, and model order
reduction have been attempted. Limitations of generating meshes and modifying these for
parametric estimation have been addressed by statistical approaches and smart solvers. Many
nature-inspired algorithms and other soft computing approaches have been employed for
electromagnetic synthesis problems. One of the recent additions is Game Theoretic
optimization.
Finally, the emergence of Computational Electromagnetics has been motivated by myriad
applications. These diverse application include but are not restricted to those in Electronics
and Communication, Wireless Propagation, Computer Hardware, Aerospace Engineering,
Biomedical Engineering, Radio-astronomy, Terahertz Technology, Photonics, etc for
modelling of devices, components, systems and even large structures. Some of the
well-known applications are in Analysis and design of radio frequency (RF) circuit, antennas
and systems, Analysis of antenna on structures, radar imaging, radar cross section
(RCS) computation and reduction, and analysis of electromagnetic wave-matter interactions
at discrete, random and periodic geometries including metamaterials. Authors are encouraged
to submit original research work in the area of Computational Electromagnetics. The content
could be either theoretical development, or specific to particular applications. This Series also
encourages state-of-the-art reviews and easy to comprehend tutorials.

More information about this subseries at http://www.springer.com/series/13885


Vineetha Joy G. L. Rajeshwari
• •

Hema Singh Raveendranath U. Nair


Fundamentals of RCS
Prediction Methodology
using Parallelized Numerical
Electromagnetics Code
(NEC) and Finite Element
Pre-processor

123
Vineetha Joy G. L. Rajeshwari
Centre for Electromagnetics (CEM) Centre for Electromagnetics (CEM)
CSIR-National Aerospace Laboratories CSIR-National Aerospace Laboratories
Bengaluru, Karnataka, India Bengaluru, Karnataka, India

Hema Singh Raveendranath U. Nair


Centre for Electromagnetics (CEM) Centre for Electromagnetics (CEM)
CSIR-National Aerospace Laboratories CSIR-National Aerospace Laboratories
Bengaluru, Karnataka, India Bengaluru, Karnataka, India

ISSN 2191-8112 ISSN 2191-8120 (electronic)


SpringerBriefs in Electrical and Computer Engineering
ISSN 2365-6239 ISSN 2365-6247 (electronic)
SpringerBriefs in Computational Electromagnetics
ISBN 978-981-15-7163-3 ISBN 978-981-15-7164-0 (eBook)
https://doi.org/10.1007/978-981-15-7164-0
© The Author(s), under exclusive license to Springer Nature Singapore Pte Ltd. 2021
This work is subject to copyright. All rights are solely and exclusively licensed by the Publisher, whether
the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of
illustrations, recitation, broadcasting, reproduction on microfilms or in any other physical way, and
transmission or information storage and retrieval, electronic adaptation, computer software, or by similar
or dissimilar methodology now known or hereafter developed.
The use of general descriptive names, registered names, trademarks, service marks, etc. in this
publication does not imply, even in the absence of a specific statement, that such names are exempt from
the relevant protective laws and regulations and therefore free for general use.
The publisher, the authors and the editors are safe to assume that the advice and information in this
book are believed to be true and accurate at the date of publication. Neither the publisher nor the
authors or the editors give a warranty, expressed or implied, with respect to the material contained
herein or for any errors or omissions that may have been made. The publisher remains neutral with regard
to jurisdictional claims in published maps and institutional affiliations.

This Springer imprint is published by the registered company Springer Nature Singapore Pte Ltd.
The registered company address is: 152 Beach Road, #21-01/04 Gateway East, Singapore 189721,
Singapore
To
Our Families
Preface

With the advent of stealth technology, precise computation of radar cross section
(RCS) has become an inevitable component in the efficient design and development
of military vehicles. However, the scattering characterization of stealth platforms is
undoubtedly the most demanding problem in modern applied electromagnetics.
When it comes to the computation of scattered fields from conducting structures,
the accuracy and efficacy of a method of moments (MoM)-based integral equation
formulation are difficult to be surpassed. Among MoM-based solvers, Numerical
Electromagnetics Code (NEC) is a versatile open-source computer program used
for the electromagnetic analysis of metallic structures in the presence of sources or
incident fields. Since its development, NEC has continued to be one of the most
widely used electromagnetic simulation codes even in the presence of commercial
MoM-based solvers. This can be attributed to the availability of the well-
documented computational engine. However, a serious drawback of NEC is the
absence of an appropriate module for the wire-grid-based meshing of geometries
and the consequent generation of segmentation data in an NEC compatible format.
This is one of the most intricate and time-consuming steps in the computation of
RCS using NEC. In this regard, this brief presents a detailed methodology for the
computation of RCS of metallic structures using a parallelized version of NEC in
conjunction with a finite element pre-processor, which has been strategically
incorporated for geometry modelling catering to NEC guidelines. It includes a
thorough overview of the theoretical background of NEC including all relevant
aspects of formulation and modelling. This book will serve as a valuable resource
for students, researchers, scientists, and engineers working in the field of RCS
predictions and measurements.

Bengaluru, India Vineetha Joy


G. L. Rajeshwari
Hema Singh
Raveendranath U. Nair

vii
Acknowledgements

First and foremost, we would like to thank the God Almighty for being the beacon
of hope during the course of this research work.
Further, we would like to thank Shri. Jitendra J. Jadhav, Director, CSIR-National
Aerospace Laboratories, Bengaluru, for the permission to write this SpringerBrief.
We would also like to acknowledge the valuable suggestions from our col-
leagues at the Centre for Electromagnetics during the course of writing this book.
But for the concerted support and encouragement of Ms. Swati Meherishi,
Editorial Director, Applied Sciences and Engineering, and Ms. Muskan Jaiswal of
Springer Nature India Private Limited, it would not have been possible to bring out
this book within such a short span of time.
We are also forever beholden to our family members for their incessant inspi-
ration which supported us to stay at rough tracks. Vineetha Joy would like to
specially thank her husband Johnu George for his constant support and encour-
agment during the course of this work.

ix
Contents

1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
2 Theoretical Background for the Computation of Radar
Cross-Section (RCS) Using NEC . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.1 Formulation of Integral Equations (IE) . . . . . . . . . . . . . . . . . . . . . 5
2.1.1 Electric-Field Integral Equation (EF-IE) . . . . . . . . . . . . . . 5
2.1.2 Magnetic-Field Integral Equation (MF-IE) . . . . . . . . . . . . . 6
2.1.3 EF-IE/MF-IE Combined Equation . . . . . . . . . . . . . . . . . . . 7
2.2 Computation of Scattered Fields . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.3 Numerical Solution of Integral Equations . . . . . . . . . . . . . . . . . . . 9
2.3.1 Selection of Basis Functions and Weight
Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.3.2 Solution of Matrix Equation . . . . . . . . . . . . . . . . . . . . . . . 11
2.4 Strategy for Parallelization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
3 Modelling Guidelines and Input/Output Formats in NEC . . . . . . . . 17
3.1 Modelling of Wires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.2 Modelling of Surfaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.3 Format of NEC Input File . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.4 Format of NEC Output File . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
Reference . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
4 Methodology for the Computation of RCS using Parallelized NEC
and Finite Element Pre-processor . . . . . . . . . . . . . . . . . . . . . . . . . .. 23
4.1 Overview of Finite Element Pre-processor . . . . . . . . . . . . . . . . .. 23
4.2 Revised Methodology for the Computation of RCS using
Parallelized NEC and Finite Element Pre-processor . . . . . . . . . .. 25
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 27

xi
xii Contents

5 Studies on Pre-selected Structures . . . . . . . . . . . . . . . . . . . . . . . . . . 29


5.1 Monostatic RCS of a PEC Dihedral . . . . . . . . . . . . . . . . . . . . . . . 29
5.2 Monostatic RCS of NASA Almond . . . . . . . . . . . . . . . . . . . . . . . 30
5.3 Monostatic RCS of Metallic Single Ogive . . . . . . . . . . . . . . . . . . 33
5.4 Monostatic RCS of Rocket-Shaped Model . . . . . . . . . . . . . . . . . . 35
5.5 Monostatic RCS of Micro Aerial Vehicle (MAV) . . . . . . . . . . . . . 37
5.6 Bi-Static RCS of Unmanned Aerial Vehicle (UAV) . . . . . . . . . . . 39
5.7 Monostatic RCS of Automobile . . . . . . . . . . . . . . . . . . . . . . . . . . 40
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
6 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53

Appendix A: Libraries used for Parallel Computations . . . . . . . . . . . . . . 55


Appendix B: Format of Command Lines in NEC . . . . . . . . . . . . . . . . . . 61
Appendix C: Format of Various Output Files . . . . . . . . . . . . . . . . . . . . . 67
Author Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
Subject Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
About the Authors

Mrs. Vineetha Joy is working as Scientist at Centre for Electromagnetics of


CSIR-National Aerospace Laboratories, Bangalore, India since March, 2016. She
obtained B. Tech in Electronics and Communication Engineering with Second
Rank from University of Calicut in 2011 and received M. Tech degree with First
Rank in RF & Microwave Engineering from Indian Institute of Technology (IIT),
Kharagpur in 2014. She is actively involved in various R&D programs spanning
domains like Computational Electromagnetics, Electromagnetic Design and
Performance Analysis of Radomes, Design and Development of broadband radar
absorbing structures (RAS), Hybrid Numerical Techniques for Scattering Analysis,
and Electromagnetic Characterization of Potential Materials for Airborne
Structures. She has authored/co-authored several research publications including
peer reviewed journal papers, symposium papers, technical documents and test
reports.

Ms. G. L. Rajeshwari obtained M. Tech. in Digital Electronics and Communi-


cation in 2019 from Department of Electronics & Communication Engineering,
Bangalore Institute of Technology, Bengaluru, India and B.E. in Electronics &
Communication Engineering, in 2017 from Nitte Meenakshi Institute of Tech-
nology, Bengaluru. Her research interest includes Computational Electromagnetics,
Hybrid Numerical Techniques for Scattering Analysis etc.

Dr. Hema Singh is working as Senior Principal Scientist in Centre for Electro-
magnetics, National Aerospace Laboratories (CSIR-NAL), Bangalore, India. She
has received Ph.D. degree in Electronics Engineering from IIT-BHU, Varanasi
India in Feb. 2000. For the period 1999-2001, she was Lecturer in Physics at P.G.
College, Kashipur, Uttaranchal, India. She was a Lecturer in EEE of Birla Institute
of Technology & Science (BITS), Pilani, Rajasthan, India, for the period
2001-2004. She joined CSIR-NAL as Scientist in January 2005. Her active areas of
research are Computational Electromagnetics for Aerospace Applications, EM
analysis of propagation in an indoor environment, Phased Arrays, Conformal
Antennas, Radar Cross Section (RCS) Studies including Active RCS Reduction.

xiii
xiv About the Authors

She has contributed in the projects not only sponsored by DRDO on low RCS
phased array, Active RCS reduction, RAS development, but also in the project
sponsored by Boeing USA on EM analysis of RF field build-up within Boeing 787
Dreamliner, She received Best Woman Scientist Award in CSIR-NAL, Bangalore
for period of 2007-2008 for her contribution in the area of active RCS reduction.
Dr. Singh co-authored 14 books, 2 book chapters, 7 software copyrights, 370
scientific research papers and technical reports.

Dr. Raveendranath U. Nair is working as Senior Principal Scientist & Head in


Centre for Electromagnetics, National Aerospace Laboratories (CSIR-NAL),
Bangalore, India. He holds a Ph.D. degree in Microwave Electronics from MG
University, India. He has over 20 years’ experience in the field of electromagnetic
design and analysis of radomes, design and development of FSS based structures
for airborne platforms, radar cross section (RCS) studies, design and development
of artificially engineered materials etc. He has contributed significantly to the
National Radome Programs including Doppler Weather Radar (DWR) radome
installed at SHAAR Sreeharikota, Fire Control Radar (FCR) radome for Jaguar
aircraft, Astra missile ceramic radome, nosecone radome for Saras aircraft, multi-
band radomes for TU-142M aircraft etc. He has authored/co-authored over 200
research publications including peer reviewed journal papers, symposium papers,
technical reports, and two books.
List of Figures

Fig. 2.1 a 1D column distribution, b 2D block-cyclic distribution . . . .. 14


Fig. 2.2 Variation in memory requirements of NEC at a particular node
with respect to number of segments . . . . . . . . . . . . . . . . . . . . .. 15
Fig. 3.1 Geometry of a wire segment . . . . . . . . . . . . . . . . . . . . . . . . . .. 18
Fig. 3.2 Surface patch a general geometry b rectangular patch
c triangular patch d quadrilateral patch . . . . . . . . . . . . . . . . . . .. 20
Fig. 4.1 Graphical user interface (GUI) of HyperMesh . . . . . . . . . . . . .. 24
Fig. 4.2 Workflow for the computation of RCS using NEC
in conjunction with HyperMesh . . . . . . . . . . . . . . . . . . . . . . . .. 25
Fig. 4.3 Procedure for computation of RCS using NEC
and HyperMesh . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 27
Fig. 5.1 Wire grid model of PEC dihedral in HyperMesh . . . . . . . . . . .. 30
Fig. 5.2 A portion of the NEC compatible segmentation data,
corresponding to PEC dihedral, extracted from Nastran file . . .. 31
Fig. 5.3 NEC input file for the computation of monostatic RCS
of PEC dihedral at 600 MHz . . . . . . . . . . . . . . . . . . . . . . . . . .. 31
Fig. 5.4 Monostatic RCS of PEC dihedral . . . . . . . . . . . . . . . . . . . . . . .. 32
Fig. 5.5 Wire grid model of NASA almond in HyperMesh . . . . . . . . . .. 33
Fig. 5.6 A portion of the NEC compatible segmentation data,
corresponding to NASA almond, extracted from Nastran file . .. 34
Fig. 5.7 NEC input file for the computation of monostatic RCS
of NASA almond (horizontal polarization) . . . . . . . . . . . . . . . .. 34
Fig. 5.8 NEC input file for the computation of monostatic RCS
of NASA almond (vertical polarization) . . . . . . . . . . . . . . . . . .. 35
Fig. 5.9 Monostatic RCS of metallic NASA almond. a Horizontal
polarization b vertical polarization . . . . . . . . . . . . . . . . . . . . . .. 36
Fig. 5.10 Wire grid model of metallic single ogive in HyperMesh . . . . .. 37
Fig. 5.11 A portion of the NEC compatible segmentation data,
corresponding to metallic single ogive, extracted
from Nastran file . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 38

xv
xvi List of Figures

Fig. 5.12 NEC input file for the computation of monostatic


RCS of metallic single ogive (horizontal polarization) . . . . . . . . 38
Fig. 5.13 NEC input file for the computation of monostatic
RCS of metallic single ogive (vertical polarization) . . . . . . . . . . 39
Fig. 5.14 Monostatic RCS of metallic single ogive. a Horizontal
polarization b vertical polarization . . . . . . . . . . . . . . . . . . . . . . . 40
Fig. 5.15 Cross-sectional view of the metallic rocket-shaped model . . . . . 41
Fig. 5.16 Wire grid model of metallic rocket shaped geometry
in HyperMesh . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
Fig. 5.17 A portion of the NEC compatible segmentation data,
corresponding to metallic rocket shaped model, extracted
from Nastran file . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
Fig. 5.18 NEC input file for the computation of monostatic
RCS of rocket-shaped model . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
Fig. 5.19 Monostatic RCS of metallic rocket shaped model. . . . . . . . . . . . 43
Fig. 5.20 CAD model of metallic MAV . . . . . . . . . . . . . . . . . . . . . . . . . . 44
Fig. 5.21 Wire grid model of Metallic MAV in HyperMesh . . . . . . . . . . . 44
Fig. 5.22 A portion of the NEC compatible segmentation data,
corresponding to metallic MAV, extracted from Nastran file . . . 45
Fig. 5.23 NEC input file for the computation of monostatic
RCS of metallic MAV . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
Fig. 5.24 Monostatic RCS of metallic MAV . . . . . . . . . . . . . . . . . . . . . . . 46
Fig. 5.25 CAD model of metallic UAV . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
Fig. 5.26 Wire grid model of metallic UAV in HyperMesh . . . . . . . . . . . 47
Fig. 5.27 A portion of the NEC compatible segmentation data,
corresponding to metallic UAV, extracted from Nastran file . . . 48
Fig. 5.28 NEC input file for the computation of bi-static
RCS of metallic UAV. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Fig. 5.29 Bi-static RCS of metallic UAV model . . . . . . . . . . . . . . . . . . . . 49
Fig. 5.30 CAD model of metallic car . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Fig. 5.31 Wire grid model of metallic car in HyperMesh . . . . . . . . . . . . . 50
Fig. 5.32 A portion of the NEC compatible segmentation data,
corresponding to metallic car, extracted from Nastran file . . . . . 51
Fig. 5.33 NEC input file for the computation of monostatic
RCS of metallic car . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
Fig. 5.34 Monostatic RCS of metallic car model . . . . . . . . . . . . . . . . . . . . 52
Fig. A.1 Structure of a general MPI program . . . . . . . . . . . . . . . . . . . . . . 58
Fig. A.2 Format of a typical PBS script . . . . . . . . . . . . . . . . . . . . . . . . . . 59
List of Tables

Table 3.1 List of geometry commands . . . . . . . . . . . . . . . . . . . . . . . . . .. 21


Table 3.2 List of program control commands . . . . . . . . . . . . . . . . . . . . .. 22
Table 5.1 Simulation and meshing parameters to be used for the
computation of monostatic RCS of PEC dihedral . . . . . . . . . .. 30
Table 5.2 Simulation and meshing parameters to be used for the
computation of monostatic RCS of NASA almond . . . . . . . . .. 33
Table 5.3 Simulation and meshing parameters to be used for the
computation of monostatic RCS of metallic single ogive . . . .. 37
Table 5.4 Simulation and meshing parameters to be used for the
computation of monostatic RCS of metallic rocket-shaped
geometry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 42
Table 5.5 Simulation and meshing parameters to be used for the
computation of monostatic RCS of metallic MAV . . . . . . . . .. 45
Table 5.6 Simulation and meshing parameters to be used for the
computation of bi-static RCS of metallic UAV . . . . . . . . . . . .. 48
Table 5.7 Simulation and meshing parameters to be used for the
computation of monostatic RCS of metallic car . . . . . . . . . . . . 50
Table B.1 Options available for Int-1 (EX command) . . . . . . . . . . . . . . . . 63
Table B.2 Options available for Int-1 (RP command) . . . . . . . . . . . . . . . . 64
Table B.3 Options available for P (RP command) . . . . . . . . . . . . . . . . . . 64
Table B.4 Options available for N (RP command) . . . . . . . . . . . . . . . . . . 65
Table B.5 Options available for G (RP command) . . . . . . . . . . . . . . . . . . 65
Table B.6 Options available for A (RP command) . . . . . . . . . . . . . . . . . . 65

xvii
Chapter 1
Introduction

The burgeoning research and developmental activities happening in various domains


like antenna design, wireless communication, stealth technologies, electromagnetic
interference/compatibility (EMI/EMC), etc., have made systems much more complex
than ever before. Considering the huge expenses related to fabrication, it is not
practically possible to modify a device if its performance is not within the stipulated
limits. Therefore, techniques capable of precise characterization of complex systems
are the need of the hour. Closed-form analytical solutions are available for only
particular problems, which can hardly be directly applied in real-world systems.
This gap has led to extensive research on computational electromagnetics (CEM)
where numerical methods are used for solving electromagnetic problems. Among
all the numerical techniques, method of moments (MoM) is highly accurate and
efficient especially for the analysis of scattering from conducting structures. In this
regard, Numerical Electromagnetics Code (NEC) is a versatile MoM-based open-
source code (written in Fortran) used for the electromagnetic analysis of metallic
structures in the presence of sources or incident fields (Burke and Poggio 1981;
Burke et al. 2004). The development of NEC has been spearheaded by Lawrence
Livermore Laboratory, USA, with the support of the Naval Ocean Systems Center
and the Air Force Weapons Laboratory.
Since its development, NEC has continued to be one of the most widely used
electromagnetic simulation codes even in the presence of commercial MoM solvers.
This attributes to the availability of the well-documented powerful computational
platform, which can be used as a foundation for further developments in the area of
electromagnetic analysis. NEC does not use many of the simplifying assumptions
employed by other high-frequency asymptotic techniques and thus provides a highly
accurate and resourceful tool for electromagnetic studies.
The numerical solution of integral equations, for the induced currents on metallic
structures by sources or impinging electromagnetic waves, is the essence of NEC.
The excitations can be defined in terms of voltage sources or plane waves of
linear/elliptical polarization. Various electromagnetic entities like non-radiating

© The Author(s), under exclusive license to Springer Nature Singapore Pte Ltd. 2021 1
V. Joy et al., Fundamentals of RCS Prediction Methodology using Parallelized Numerical
Electromagnetics Code (NEC) and Finite Element Pre-processor,
SpringerBriefs in Computational Electromagnetics,
https://doi.org/10.1007/978-981-15-7164-0_1
2 1 Introduction

networks/transmission lines, perfect/imperfect conductors, perfectly or imperfectly


conducting ground planes, etc., can be easily included in NEC. The output of NEC can
be tailored for specific requirements like induced current and charge densities, near
magnetic or electric fields, impedance or admittance, radiated fields, etc. Further-
more, many of the commonly encountered performance parameters such as gain and
directivity are also available. Therefore, the program is well suited for both antenna
analysis and radar cross section (RCS) predictions (Burke and Poggio 1981).
Although numerical methods like MoM can be used in a wide range of radiation
and scattering problems, the computational complexity increases with increase in the
electrical size of the structure under consideration. The memory requirements have
been found to exhibit a quadratic growth with increase in the number of segments. A
parallelized version of NEC has been introduced (Rubinstein et al. 2003, Rubinstein
2004) to tackle this problem. The two most time-consuming steps in evaluating the
response of a structure, viz. (i) computation of the interaction matrix elements and
(ii) solution of the resultant matrix equation, have been efficiently parallelized in
Rubinstein et al. (2003).
In addition to the memory constraints, a serious issue in the computation of RCS
using NEC is the wire-grid-based meshing of geometries and the generation of NEC
input file with the segmentation data included in accordance with NEC guidelines.
This is an extremely difficult and error-prone process, especially for complex struc-
tures (Hubing et al. 1994; Toit and Davidson 1995; Ross et al. 1999), and further-
more, the code will often run in the presence of such errors and produce results that
are plausible but incorrect. Therefore, an appropriate pre-processor is essential for
the accurate prediction of RCS using NEC.
The aim of this book is to provide a strong foundation for the computation of
radar cross section (RCS) of aerospace structures using a parallelized version of
NEC in conjunction with a finite element pre-processor. This book has five chapters.
It begins with a thorough overview of the theoretical background of NEC including
all relevant aspects of formulation of integral equation and its numerical solution.
The strategy used for parallelization of the code is also included in Chap. 2. The
geometrical modelling of structures and the input/output formats are described in
Chap. 3. The revised methodology for the computation of RCS incorporating all the
required steps and details is discussed elaborately in Chap. 4 followed by case studies
and validations in Chap. 5.

References

Burke, G.J., and A.J. Poggio. 1981 January. Numerical Electromagnetics Code (NEC)—Method
of Moments, Lawrence Livermore National Laboratory, Livermore, California, Technical
Document, Rep. UCID-18834, 719.
Burke, G.J., E.K. Miller, and A.J. Poggio. 2004 June. “The Numerical Electromagnetics Code
(NEC)—A Brief history.” In Proceedings of IEEE Antennas and Propagation Society Symposium,
vol. 3, 2871–2874, Monterey, California.
References 3

Hubing, T.H., C.H.H. Lim, and J. Drewniak. 1994 March.“A Geometry Description Language for
3D Electromagnetic Analysis Codes.” In Proceedings of the 10th Annual Review of Progress in
Applied Computational Electromagnetics, Monterey, CA, 417–422.
Rubinstein, A., F. Rachidi, M. Rubinstein, and B. Reusser. 2003. A Parallel Implementation of NEC
for the Analysis of Large Structures. IEEE Transactions on Electromagnetic Compatibility 45
(2): 177–188.
Ross, J.E., L.L. Nagy, and J. Szostka. 1999. “CAD tools for vehicular antennas.” In Proceedings of
the 1999 IX National Symposium of Radio Science, Poznan, Poland, March 16–17.
Rubinstein, A. 2004 January. Simulation of Electrically Large Structures in EMC Studies: Applica-
tion to Automotive EMC, Ph.D. dissertation, École polytechnique fédérale de Lausanne (EPFL),
https://doi.org/10.5075, 156.
Toit, C.D., and D.B. Davidson. 1995. Wiregrid: A NEC2 Pre-processor. Journal of Applied
Computational Electromagnetics Society 10: 31–39.
Chapter 2
Theoretical Background
for the Computation of Radar
Cross-Section (RCS) Using NEC

The entire theoretical framework of NEC starting from the definition of integral
equations to the solution of matrix equation is briefly introduced in this chapter. The
basis functions used for representing the unknown current densities and the weight
functions used for testing are also described in detail. Further, the strategy used for
parallelization of NEC (Rubinstein et al. 2003) is also included here.

2.1 Formulation of Integral Equations (IE)

NEC employs both an electric-field integral equation (EF-IE) and a magnetic-field


integral equation (MF-IE) for evaluating the electromagnetic behaviour of conducting
bodies under consideration. EF-IEs are mainly used in the case of thin-wire geome-
tries with small conductor volumes. On the other hand, MF-IE, which cannot be
applied for the thin-wire case, is more suitable for voluminous structures with smooth
surfaces. EF-IE and MF-IE are coupled together for a structure containing both wires
and surfaces. The theoretical formulation presented in the following sections is based
on the work reported in (Burke and Poggio 1981).

2.1.1 Electric-Field Integral Equation (EF-IE)

The integral representation for the electric field due to a surface current distribution


J S on a perfectly conducting body (S) is given by,
     −

→ →
 r ) = (− jη/4π k)
E( JS r  · Ḡ¯ r, r  dA (2.1)
S

© The Author(s), under exclusive license to Springer Nature Singapore Pte Ltd. 2021 5
V. Joy et al., Fundamentals of RCS Prediction Methodology using Parallelized Numerical
Electromagnetics Code (NEC) and Finite Element Pre-processor,
SpringerBriefs in Computational Electromagnetics,
https://doi.org/10.1007/978-981-15-7164-0_2
6 2 Theoretical Background for the Computation …

where,
     
G r, r = k 2 I + ∇∇ g r, r
 
 
  e− j k r −r 

g r, r = 
 
 
r − r 


√ μ0
k = ω μ0 ε0 ; η = ; I = x̂ x̂ + ŷ ŷ + ẑ ẑ
ε0

r and r represent observation point and source point, respectively.




The integral equation for the current induced on S by an incident field E Inc is
obtained by imposing the following boundary condition for r ∈ S:

→ −

r ) × E Sca (
n̂( r ) + E Inc (
r) = 0 (2.2)



r ) is the unit normal vector to the surface at r and E Sca is the scattered field
where n̂(

→ −
→Sca
due to the induced current J s . Substituting (2.1) for E in (2.2) yields the final
EF-IE as,


→ →   2
−   
r ) × E Inc (
−n̂( r ) = (− jη/4π k) × n̂(
r) × J s (r ) · k I + ∇∇ g r, r dA
S
(2.3)

Since NEC employs thin cylindrical wires for representing the conducting surface
(S), the vector integral equation in (2.3) will be transformed to a scalar integral
equation. The assumptions applicable for the thin-wire approximation used in most
of the cases are given below:
(i) Currents in the transverse direction can be ignored with respect to the currents
along axial direction of the wire.
(ii) Circumferential variation in the currents along axial direction can be ignored.
(iii) The unknown current may be expressed as a filament along the wire axis.
(iv) The boundary condition pertaining to the electric field is required to be imposed
only in the axial direction.

2.1.2 Magnetic-Field Integral Equation (MF-IE)

The integral representation for the magnetic field due to a surface current distribution


J s on a perfectly conducting body is given by,
2.1 Formulation of Integral Equations (IE) 7
   

→ −→ 
H (
r ) = (1/4π ) J s (r ) · ∇  g r, r dA (2.4)
S

where the differentiation is with respect to the integration variable r .



→ −

If an external incident field H Inc induces a current J s on S, then the total magnetic
field within the perfectly conducting body vanishes. Hence, the boundary condition
for r immediately within the surface S can be expressed as,

→Sca −

H ( r ) + H Inc (
r) = 0 (2.5)


→ −

where H Inc is the incident field without the structure and H Sca is the scattered field


due to the induced current J s given by (2.4).


Substituting H Sca in (2.5) yields the final MF-IE as,
     

→ −
→ 
r0 ) × H Inc (
−n̂( r0 ) = n̂(
r0 ) × (1/4π ) lim J s r × ∇  g r, r dA (2.6)
r→
r0
S

r0 ) is the normal vector directed outwards at the surface point r0 .


where n̂(
Similar to EF-IE, the vector integral Eq. (2.6) can also be resolved into two scalar
expressions in the directions of the orthogonal surface vectors tˆ1 and tˆ2 where

r0 ) × tˆ2 (
tˆ1 ( r0 ) = n̂(
r0 ) (2.7)

2.1.3 EF-IE/MF-IE Combined Equation

The EF-IE/MF-IE combined equation is used for a structure consisting of both wires
and surfaces. In such cases, r in (2.5) is restricted to wires and the integral for

→Sca
E ( r ) extends over the complete structure. Similarly, r0 in (2.6) is restricted to


surfaces and the integrals for H Sca (
r ) extends over the entire structure. The integral
can be simplified by using the thin-wire approximation on wires.
The final coupled integral equations (Burke and Poggio 1981) are,
    2   

→  ∂  dS 
−ŝ · E Inc (
r ) = (− jη/4π k) I (s  )
k 2 ŝ · s  − g 
r , r
∂s ∂s 
L
    
→     2 
− ∂ 
− ( jη/4π k) Js r k ŝ − ∇  g r, r dA ;
∂s
S1

for r on surfaces of wire (2.8)


8 2 Theoretical Background for the Computation …

  

→  
tˆ2 ( r ) = (−1/4π )tˆ2 (
r ) · H Inc ( r) · I s  sˆ × ∇  g r, r dS 
L

− (1/2)tˆ1 (
r ) · JS (
r ) − (1/4π ) tˆ2 (
r)
S1

→  
· J S (r ) × ∇  g r, r dA ;
for r on surfaces excluding wires (2.9)
  


−tˆ1 ( r ) = (1/4π )tˆ1 (
r ) · H Inc ( r) · I (s  ) sˆ × ∇  g r, r dS 
L

− (1/2)tˆ2 (
r ) · JS (
r ) + (1/4π ) tˆ1 (
r)
S1
 
· JS (r ) × ∇ g r, r dA ;


for r on surfaces excluding wires (2.10)

where s is the parameter corresponding to distance along


 the wire axis at r and ŝ
is the unit vector tangent to the wire axis at r. and represents integration over
L S1
wires and integration over surfaces excluding wires, respectively.

2.2 Computation of Scattered Fields

The scattered fields due to induced currents can be computed using a far-field approx-
imation. This is applicable when the observation points are far apart from the induced
current distribution by distances much larger compared to the wavelength as well
as the dimensions of the current distribution. The far-zone scattered fields from a
conducting body consisting of a wire region with contour L having linear current

→ −

distribution I (s) and a surface region S with surface current density J s ( r ) are
given by,


→ ⎨   − → 


r0 ) = ( jkη/4π ) × (e− jkro /ro ) ×
E ( (k · I (s))k − I(s) e j k·r

L

   ⎬
 −→  −
→ 
+ (k · J s (
r ))k − J s (
r ) e j k·r dA (2.11)

S
2.2 Computation of Scattered Fields 9


where k = |rroo | , k = 2π
λ
and ro is the position of the observation point.
The contour integral can be computed in closed form over each wire segment
for the cosine, sine, and constant components of the basis functions and this in turn
becomes a summation over the wire segments. Since the surface current on the patch
is represented by a delta function at the centre of the patch, the second integral
transforms to a summation over the patch elements.


For an incident plane wave with electric field E Inc , the NEC program constants are
defined such that the radar cross-section (RCS) σ/λ2 is displayed under the column
gain in the output file.
− 
→Sca 2
σ E 
= 4π 
→ 2
(2.12)
λ2 −
 E Inc 



where E Sca is given by (2.11).

2.3 Numerical Solution of Integral Equations

The integral equations described in the previous sections are solved numerically in
NEC by using the method of moments (MoM). Method of moments starts with a
conventional linear-operator equation as,

L( f ) = u (2.13)

where f is an unknown response, u is a known excitation and L is an integral operator.


The formulation proceeds by expressing the unknown function f as a sum of
known basis functions f j (defined in the domain of the integral operator L) as,


N
f = βj f j (2.14)
j=1

The unknown coefficients (β j ) are obtained by taking the inner product of (2.13)
with a set of weighting functions wi , i.e.

wi , L( f ) = wi , u ; i = 1, 2, . . . , N (2.15)

The inner product in the present case can be defined as,



 p, q = r ) q(
( p( r ))dA
S
10 2 Theoretical Background for the Computation …

where integration is over the surface of the structure.


Substituting (2.14) in (2.15) gives


N
  
β j wi , L f j = wi , u ; i = 1, 2, . . . , N (2.16)
j=1

The above equation can be expressed in matrix form as,

[G] [B] = [U ] (2.17)


  
where G i j = wi , L f j ; B j = β j ; Ui = wi , u .
The solution to (2.17) is

[B] = [G]−1 [U ] (2.18)

2.3.1 Selection of Basis Functions and Weight Functions

The critical factors determining the accuracy and convergence of MoM are the weight
functions and the basis functions. NEC employs different basis and weight functions
unlike Galerkin’s method. The weighing functions (wi ) are chosen as a set of delta
functions in accordance with the collocation method, where the integral equation is
enforced at specific sample points on the geometry, i.e.

wi (
r ) = δ(
r − ri )

where {ri } are the set of points on the conducting surface. In NEC, conducting wires
are sub-divided into short linear segments where a sample point is defined at the
centre of each segment. On the other hand, surfaces are represented by a set of flat
patches with a sample point defined at the patch centre.
The linear current density on wires and surface current density on surfaces are
the unknown quantities in NEC which are to be represented by a set of known basis
functions. NEC uses sub-domain basis functions to simplify the calculation of inner
product integrals and also to ensure that the matrix G is well conditioned.
In the case of wires, the current on each segment is approximated as a sum of
three terms, viz., a sine, a cosine, and a constant. The total current on ith segment
can thus be expressed as,

Ii (s) = Ai sin k(s − si ) + Bi cos k(s − si ) + Ci , |s − si | < i


(2.19)
2
2.3 Numerical Solution of Integral Equations 11

where si is the value of s at the centre of segment i, i is the length of segment i,


and k is the free space wave number.
The additional equations required to evaluate all the unknown coefficients in (2.19)
for each and every segment will be provided by the local boundary conditions. The
continuity of current and charge is one such condition to be enforced at a junction
of two segments with uniform radius. In the case of a junction with more than
two segments of different radii, the continuity of current may be generalized to
Kirchhoff’s current law.
The final solution requires the computation of the electric field at each segment
due to these currents. NEC uses three approximations of the integral equation kernel
for the computation of fields. These include (i) thin-wire approximation for most
of the cases, (ii) extended thin-wire approximation for thick wires, and (iii) current
element approximation for the case where segments are separated far apart by a
considerable distance. In all the three methods used for the evaluation of field at
a particular segment due to the induced current on a different segment, the field
computation is performed on the surface of the observation segment.
NEC models surfaces using small flat patches and the unknown surface current
on each patch are expressed as a summation of a set of pulse functions except where
a wire connection is present. When a wire connects to a surface, a special local
boundary condition has to be enforced as explained in (Burke and Poggio 1981).
The unknown surface current density over M p patches can be written as,


→  Mp
 
J s (
r) = (J1i t 1i + J2i t 2i )vi (
r) (2.20)
i=1

   
where t 1i = t 1 (
ri ), t 2i = t 2 (
ri ), ri is the position of the centre of patch i, vi (
r) = 1
for r on patch i and 0 otherwise. Thus, every patch has two unknown quantities and
the equation is enforced for each vector component at the sample point in a given
patch.

2.3.2 Solution of Matrix Equation

The matrix equation obtained by appropriately substituting the basis functions and
weight functions in (2.16) has to be finally solved for the unknown currents. The
equation can be alternately expressed as,

[G] [I ] = [E] (2.21)

where [G] is the interaction matrix, [I ] is the column vector with unknown currents,
and [E] is the excitation vector. NEC employs Gauss elimination method for solving
(2.21). Here, the matrix G is first factorized into the product of upper and lower
12 2 Theoretical Background for the Computation …

triangular matrices, i.e.

[G] = [L] [U ] (2.22)

The matrix equation then becomes,

[L] [U ] [I ] = [E] (2.23)

Let,

[L] [F] = [E] (2.24)

where,

[U ] [I ] = [F] (2.25)

Equation (2.24) is first solved for F by forward substitution. Then (2.25) is solved
for I by backward substitution. The order of the matrix Eq. (2.21) for a structure
having N w wires segments and N p patches is given by,

N = N p + Nw

The final matrix Eq. (2.21) can be written in terms of sub-matrices as,
    
P Q IW EW
= (2.26)
R S IP HP

where I w is the column vector of segment basis function amplitudes, I p is the column
vector of patch-current amplitudes, E w is the LHS of (2.8) evaluated at segment
centres, and H p is the LHS of (2.9) and (2.10) evaluated at patch centres.
The sub-matrix element Pij denotes the electric field at the centre of segment i
due to the basis function centred on segment j. On the other hand, the sub-matrix
element S ij denotes the tangential magnetic-field component at patch u due to the
surface current pulse on patch v where
 
i −1
u = 1 + Int ; Int[ ]indicates truncation
2
 
j −1
v = 1 + Int
2

The elements in sub-matrix Q denote the electric fields due to surface current
pulses. On the other hand, the elements in R denote magnetic fields due to segment
basis functions.
2.4 Strategy for Parallelization 13

2.4 Strategy for Parallelization

The structure of NEC can be broadly divided into two parts, namely the input section
and the calculation section. The input section reads in the geometrical data pertaining
to the structure under consideration as well as the program control commands. The
calculation section on the other hand computes the elements of the interaction matrix
[G] which holds the geometrical and electromagnetic information about the model.
This section also deals with the solution of (2.21) by Gauss elimination method. This
final step is carried out in two routines. The first one performs the computationally
intensive task of LU decomposition of [G]. Later, it is used along with the excitation
vector [E] in another routine to obtain the solution. The computation of the constituent
elements of [G] and evaluation of the final solution to the matrix equation accounts
for more than 90% of the computation time (Burke and Poggio 1981). These parts
of the code need to be parallelized for the same reason.
The parallelized version of NEC used in this book is in accordance with (Rubin-
stein et al. 2003) where the subroutines responsible for LU decomposition (FACTR)
and solution of matrix equation (SOLVE) have been revised to work in a parallel envi-
ronment (Rubinstein 2004). This version distributes the interaction matrix among the
available processors and the matrix elements are locally and individually calculated
by their respective processors. The final system of linear equations is also solved
using a parallel version of the Gauss-Doolittle algorithm. In this way, the combined
processing power and memory of several processors working together as a cluster can
be efficiently utilized. The code can automatically detect the number of processors
and hence exploit the available memory and computing resources. The parallelized
code (Rubinstein et al. 2003) is designed to run on distributed memory parallel super-
computers and it can work on any system/cluster supporting conventional message
passing parallel environments such as Message Passing Interface (MPI) and Parallel
Virtual Machine (PVM). Since it is completely based on the original NEC code, the
input files need not be changed and the output format remains the same too.
The efficient distribution of the information stored in the matrices among the
available processors is very much critical for the successful parallelization of LU
decomposition and the parallel solution of the system of equations. For this reason,
the interaction matrix and the excitation vector have to be divided into smaller sub-
matrices which will be local to the respective processors. In (Rubinstein et al. 2003),
two-dimensional block-cyclic distribution has been used to uniformly distribute the
computation effort among all available processors and it is discussed in the following
paragraphs.
To illustrate two-dimensional block-cyclic distribution, consider a parallel
computer with four processors numbered as 1, 2, 3, and 4, as shown in Fig. 2.1. A
matrix is to be distributed among them. In a conventional way, as shown in Fig. 2.1a,
the matrix is distributed column wise to each processor. For ease of understanding,
the contents of processor 1 are shaded blue. Here, once all the elements that corre-
spond to processor 1 are assigned, it will become idle for the rest of the time. The
14 2 Theoretical Background for the Computation …

Fig. 2.1 a 1D column


distribution, b 2D
block-cyclic distribution

1 2 3 4 (a)

1 2 1 2 1 2 1 2
3 4 3 4 3 4 3 4
1 2 1 2 1 2 1 2
3 4 3 4 3 4 3 4
(b)
1 2 1 2 1 2 1 2
3 4 3 4 3 4 3 4
1 2 1 2 1 2 1 2
3 4 3 4 3 4 3 4

same situation repeats for the rest of the processors too. On the other hand, two-
dimensional block-cyclic distribution as shown in Fig. 2.1b guarantees that, on an
average, all the processors remain active from the start to the end of the algorithm.
A square processor grid is a desirable attribute for better efficiency as it maintains
equitable distribution of the computational load. But, in a practical scenario, the
number of processors varies from one machine to another. Therefore, the squarest
possible processor grid on the running platform is automatically calculated at the
beginning of the algorithm. ScaLAPACK library has been used in (Rubinstein et al.
2003) for implementing the block-cyclic decomposition in NEC. A brief description
of ScaLAPACK and MPI is given in Appendix A.
It is noted that the formation of the interaction matrix on a single processor
and then the distribution among other processors is not an efficient solution for
memory problem. The interaction matrix computation is therefore parallelized in
such a manner so as to produce an already distributed matrix as it is being computed.
Each processor on the grid will have adequate information required to determine
the elements of [G] once the input section is completed. An additional function
is employed to notify each processor whether the element (i, j) being calculated
corresponds to its local sub-matrix. In that case, another sub-routine will compute
the local values of (i, j) for the same element so as to be assigned to the corresponding
local matrix. Once the matrix generation is completed, each processor will hold its
local version of the interaction matrix and there would not be any further need for
2.4 Strategy for Parallelization 15

communication among the processors. Once the matrix equation is solved, data is
compiled by a block-cyclic composition routine and all the processors become aware
of the solution.
An approximation of the memory requirements in bytes at a local node for the
parallelized model of NEC is given by (Rubinstein et al. 2003),
 
16 × Nseg
2

Memory =
P
where P is the number of processors being used during the execution of the program
and Nseg is the number of segments. On the other hand, the original version of NEC
which is not parallelized requires (16 × Nseg
2
) bytes of memory. Figure 2.2 presents
the difference in memory requirements at a particular node between the original
version of NEC and its parallelized counterpart.
Considering a structure meshed with a particular number of segments, the vari-
ation in the run time of NEC with respect to the number of processors has been
studied in (Rubinstein et al. 2003), keeping all other parameters fixed. When the
number of processors was increased from one to four, it has been found that the
total run time reduced by almost 98%, thereby emphasizing the efficiency of paral-
lelization. However, the number of processors has to be judiciously chosen based on
the electrical dimensions of the structure under consideration, as employing large
number of processors may also increase the run time due to communication delay.

700

600 Original NEC


Memory requirements (MB)

Parallel NEC (P=2)


Parallel NEC (P=4)
500 Parallel NEC (P=6)
Parallel NEC (P=8)
400

300

200

100

0
1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000

Number of segments

Fig. 2.2 Variation in memory requirements of NEC at a particular node with respect to number of
segments
Another Random Document on
Scribd Without Any Related Topics
Wirthshaus tritt. Zu Hause, auf Reisen, wo es auch sei!
Schreiben Sie ihnen, und Sie haben den ersten Posttag
Antwort. Die Gräfin Augusta — vom Morgen bis im Abend
laufen die Depeschen bei ihr ein, wie bei einem
Staatsminister, und werden sorgfältiger abgefertigt als in
einer Kanzlei.
Letzthin allegorisirte ich darüber mit Tellow (der
Liebesname seines und meines Freundes Cramer).
Wo ist nun die Gräfin wieder? fragte ich (Klopstock).
Cramer: Oben; schreibt Briefe.
Klopstock: Das ist wahr! Die Stolbergs! Sie liegen am
Briefschreiben recht krank danieder.
Cramer: Freilich, es ist eine Krankheit zum Tode.
Klopstock: O! sie sind schon gestorben.
Cramer: Und begraben dazu.
Klopstock: Was? Sie sind schon auferstanden.
Cramer: Ei! sie sind schon selig.
Klopstock: Ja, nun — kann ich nicht weiter.
Hierüber kommt die Gräfin herunter.
Wir sprachen, sagt ihr Klopstock, eben zusammen von
Ihrer Krankheit, Ihrem Begräbniß, Ihrer Auferstehung,
Ihrer Seligkeit!
Wie so?
Ja, gestehen Sie es nur, schöne Gräfin, Ihr Briefschreiben
ist doch eine wahre Krankheit, eine Schwachheit, eine
Seuche!
Sie mögen aber doch wol selbst gern Briefe haben?
Das mag ich wohl; — o, das Briefel e s e n ist eine ganz
vortreffliche Sache; aber das S c h r e i b e n! Es ist eine
Schwachheit, ein Fehler, sage ich, aber eine nicht eben
unliebenswürdige Schwachheit! Wenn sich die Briefe, die
Antworten wenigstens, nur selbst schrieben!
Meine Anekdote ist zu Ende. Die mußte man freilich von
Klopstock selbst erzählen hören!

In einem spätern Briefe, vom 30. März, an Ludwig


findet sich eine hübsche Stelle, die unter Weglassung
anderer nicht hierher gehöriger Bemerkungen hier
noch folgen möge:
Bald, vielleicht wenige Stunden später, als Sie diese Zeilen
erhalten, drücke ich Sie an meine Brust und sage Ihnen
mündlich, wie sehr ich Sie liebe und verehre. Mehr wie je.
Es ist mit der Freundschaft wie mit der Liebe. Die
Entfernung tödtet schwache, sie stärkt die echte und
wahre. Den Frauen Ihres Hauses küsse ich mit Verehrung
die schönen Hände.

Aus Amsterdam und dann von der leipziger Messe


nach Altenburg wieder zurückgekehrt, schreibt
Brockhaus in dem schon erwähnten Briefe vom 21.
Mai an Bornträger:
Meine freundschaftlichen Verhältnisse mit Ludwigs,
Hempels und Andern dauern ununterbrochen fort und
consolidiren sich selbst immer mehr. Seit einer Reihe von
Jahren ist dieser Sommer der erste, wo ich meines Lebens
wieder froh bin. Die Pfingstfeiertage werde ich mit meinen
Freunden eine Tour nach Dresden machen. Fritz und Lina
sind in demselben Hause, wo ich wohne, in Kost und unter
Aufsicht. Wahrscheinlich werde ich Fritz, um ihm mehr
Reibung zu geben, hier in der Nähe in ein sehr gutes
Institut thun. Lina behalte ich aber bei mir. Was aus
Sophiechen werden soll? Ich habe von Ihnen beständig
Nachricht erwartet über ihre Unterbringung in dem Dorfe
bei Muiden. Wo das arme Kind gut ist, da ist es mir recht
bis dahin, daß ich es von dort zu mir nehmen kann. Wäre
es einmal hier, so wäre es gut aufgehoben. Aber wie
hierher bringen?

Die für die Pfingstfeiertage 1811 beabsichtigte Reise


mit Ludwigs nach Dresden fand erst Mitte Juli statt.
Brockhaus hatte von derselben vielen Genuß,
besonders von dem Aufenthalte in Dresden selbst,
das er wol zum ersten male sah, »dieser an
Kunstschätzen und Naturschönheiten einzigen
Stadt«, wie er schreibt, ebenso von dem
Zusammentreffen mit interessanten Leuten. Unter
diesen nennt er einen Baron von Heinse aus Lübeck,
mit dem zusammen er Ludwigs, die sich in Dresden
von ihm getrennt hatten, um einen längern
Aufenthalt in dem Bade Teplitz zu nehmen, dort
besuchte und dann nach Dresden zurückkehrte.
Anfang August war er wieder in Altenburg, während
Ludwigs erst am 4. September wieder dort eintrafen.
Die Schwester von Frau Ludwig, Jeannette von
Zschock, hatte nebst ihrer Freundin Karoline Hempel
ebenfalls an der Reise nach Dresden und Teplitz
theilgenommen, Beide waren aber, wie es scheint, in
Brockhaus' Begleitung gleich mit nach Altenburg
zurückgereist.
Diese gemeinschaftliche Reise und die unmittelbar
darauffolgende Zeit brachten in Brockhaus einen
Entschluß zur Reife, den er schon lange mit sich
herumtrug und der auch oft zwischen den Zeilen
seiner von uns mitgetheilten Briefe an Herrn und
Frau Ludwig durchschimmert: er verlobte sich mit
Jeannette von Zschock, und der Verlobung, die nach
der Rückkehr ihrer Schwester und ihres Schwagers
zuerst nur im Stillen gefeiert, bald darauf aber auch
öffentlich erklärt wurde, folgte gegen Ende des
nächsten Jahres die Verheirathung.
Jeannette von Zschock (mit ihren vollen Vornamen
Johanne Charlotte Luise Rosine, aber gewöhnlich nur
die französische Form des erstern führend) war am
7. September 1775 in Offenbach geboren und lebte
seit dem Tode ihrer Mutter und ihres Vaters, der
Rittmeister in schwäbischen Diensten gewesen war,
bei ihrer Schwester in Altenburg. Sie stand bei ihrer
Verlobung im siebenunddreißigsten Lebensjahre,
Brockhaus im vierzigsten.
Die einzigen Mittheilungen über die mit der
Verlobung zusammenhängenden Umstände finden
sich wieder in einem Briefe von Brockhaus an
Bornträger. Er schreibt diesem aus Altenburg vom 30.
August 1811:
Heute endlich die Beantwortung Ihrer mehrmals
geäußerten Wünsche, mein jetziges inneres Leben zu
kennen, meine Verhältnisse hier zur Welt, zu meinen
Freunden. Je wichtigere Nachrichten ich Ihnen über das
Höchste im Leben mitzutheilen habe, je mehr haben Sie
Recht, darüber etwas zu wissen, da Sie mit seltener
Freundschaft mein Schicksal theilen. Es hat sich in diesen
Tagen viel entschieden.

Hier folgt die früher schon mitgetheilte Stelle über


die in dieser Zeit von der Hofräthin Spazier
gemachten neuen Anknüpfungsversuche, und daran
schließen sich folgende von uns dort absichtlich noch
weggelassenen Worte:
Da aber der Verstand, beleidigte Ehre, Pflichtgefühl und
auch zarte und edle Neigung für ein anderes weibliches
Wesen mich stärken und schützen, so werde ich der
Sirenenstimme, die von der Spree her zu mir
herüberschallt, nicht folgen.

Er fährt dann fort:


Mein Verhältniß hier zur Welt ist im ganzen noch dasselbe,
wie ich es Ihnen geschildert. Innige Freundschaft mit allen
Gliedern des Ludwig'schen Hauses ist jedoch das, was
mich allein sehr anzieht. Sie sind es auch allein, die mich
ganz verstehen und würdigen. Ich habe hier nämlich
wieder das Schicksal, daß viele Menschen gegen mich
sind, daß mich diese für stolz, üppig, eingebildet und Gott
weiß wofür Alles halten, wozu ich freilich durch mein
schneidendes, auch oft sonst nie vorsichtiges Betragen
Veranlassung gegeben habe. Ich bin über die Ursachen
und die einzelnen Gravamina lange in Unsicherheit
gewesen, da ich nur die Spuren in den Folgen entdeckte,
ohne die Ursachen errathen zu können, da die Winke, die
ich von einer Seite erhielt, nicht hinreichten, mir die
nöthige Aufklärung zu geben. Jetzt kenne ich aber alle
Fäden der geheimsten Verhandlungen darüber und auch
alle Intriguen, die dabei stattgefunden und -finden.
Mein Genius, der mir jene Winke und jetzt alle
Offenbarungen gegeben hat; der mein Interesse vom
ersten Augenblicke, daß ich hier vor einem Jahre
aufgetreten bin, zum eigensten gemacht hat; dem ich und
die Hofräthin alles Gute und Liebe verdanken, das wir hier
genossen; der mich und sie mit gleicher Energie
verfochten und vertreten; der durch einen wunderbar
sympathetischen Zug sich zu mir wie ich mich zu ihm
hinneigte, als auch noch nicht die allerentfernteste
Möglichkeit da war, daß je ein näheres Verhältniß eintreten
könnte — dieser Genius ist jenes herrliche Mädchen,
Ludwig's Schwägerin, Fräulein Jeannette von Zschock —
seit einer Woche meine still Verlobte! Sie wird mir fürs
Leben angehören, wenn ich es vermag, mein bürgerliches
Schicksal ganz zu ordnen, die Einwilligung Ludwig's und
ihrer Schwester, die noch nicht von Teplitz zurück sind, zu
erhalten und die Welt ganz mit mir zu versöhnen. Wir
werden aber Vieles zu kämpfen haben, ehe wir ans Ziel
kommen.
Unsere Wahlverwandtschaft hat um so weniger
unbeobachtet bleiben können, da durch die Eifersucht der
Hofräthin, die zu einer Zeit, als der Gedanke daran zu den
Märchen aus dem Monde gehörte, mich und die arme
Jeannette aufs Blut damit verfolgte, dies unser Verhältniß
die allgemeinste Aufmerksamkeit auf sich zog, da es
psychologisch allerdings höchst interessant war und jetzt
von neuem die Behauptung Schubert's und Anderer
gewissermaßen bestätigt, wie diese Art Nervenkranker die
Gabe der Voraussehung und Voraussagung haben.
Außerordentlich ist's, daß sie im Wahnsinn ihres Fiebers
prophetisch Alles ausgesprochen hat. »Ich bin«, sagte sie,
indem sie unsere Hände zusammenlegte, »Donna Elvira,
mein Fräulein; ich werde nun gehen«....
Meinerseits bin ich überzeugt, daß meine Kinder eine
vortreffliche Mutter und Erzieherin, ich eine edle Freundin
und treue Genossin fürs Leben errungen habe, wenn es
mir gelingt, unsere Verbindung zu vollenden. Sie wissen,
wie verarmt mein Leben war und wie es das
außerordentlichste Glück ist, wenn ich es auf diese Weise
neu und schön ordnen kann. Ich gedeihe nur in einem
edeln Familienkreise, und ohne solchen bin ich nichts. Und
was wird und kann aus meinen Kindern werden, wenn sie
nicht wieder eine edle Mutter finden? Es ist der
lebhafteste Wunsch meiner Freundin, die kleine Sophie
von Amsterdam herüberzuhaben, und es werden daher
ernste Ueberlegungen stattfinden müssen, wenn Sie
herüberkommen, wie dies zu bewerkstelligen.

Der Schluß des Briefs enthält eine anziehende


Schilderung des altenburger Vogelschießens, eines
damals mit weit mehr Glanz als jetzt gefeierten
Volksfestes:
Ich würde Ihnen diesen Brief schon vor acht Tagen
geschrieben haben, wo er freilich noch nicht so klar und
bestimmt hätte melden können, was er jetzt enthält, wenn
nicht in dieser Zeit gerade das hiesige große
Vogelschießen stattgefunden, das jede geregelte Arbeit
beinahe unmöglich macht. Sie können sich keinen Begriff
davon bilden, mit welchem Pompe, mit welchen
Feierlichkeiten es begleitet ist, und wie sich Geschmack
und alle schönen Künste vereinigen, die Belustigungen
dabei zu veredeln und zu verschönern. Dies Jahr war noch
eine neue Loge erbaut worden, in der sich nun die Elite
der Gesellschaft versammelte, wo des Morgens Déjeuner
dansant, dann Dîner, Abend Bal paré, Spiel und Souper
unter den Colonnaden des Saals und in den
Nebenzimmern war. An einem Tage in den Zwischenzeiten
war noch Lotterie für Damen, an einem andern Tage
Concert. Jeannette gewann auf zwei von mir geschenkte
Lose zwei Ringe! Sie können denken, wie glücklich uns
dieser Zufall oder diese Schicksalsdeutung machte. Blos in
dieser Loge speisten gewöhnlich 4-500 Personen. Ich
glaube nicht, daß es irgendwo brillantere oder
angenehmere Bälle und Partien geben könnte, als es die
hier waren. Aus der ganzen Gegend bis von Dresden her
hatten sich lebenslustige Fremde in außerordentlicher Zahl
eingefunden, die täglich ab- und zuwogten. Dazu die
zahlreichen Buden auf der an einer sanften Anhöhe
gelegenen Vogelwiese, von der man eine wunderschöne
Aussicht hat; die herrliche Witterung, die schönen
mondhellen Nächte, der Jubel der Volksmengen, denen
diese Woche das ist, was den Römern ihr Carneval; das
Werfen der Schwärmer, der Raketen, womit sich Jung und
Alt amusirt; das ewige Musiciren von zwanzig Orten her,
das Trommeln bei jedem Schusse, der den Vogel
verwundet; die militärische Haltung aller Freunde und
Bekannten, die sämmtlich in ihren ebenso
geschmackvollen als wohlkleidenden Uniformen
(dunkelgrün aufs brillanteste mit Silber gestickt,
französischer Offiziersschnitt) mit großen russischen Hüten
und rothen Schwungfedern erscheinen; die geputzten
Weiber und Mädchen, von denen es wimmelt.
Sie wissen, wie arm man in Holland an allem ist, was
Vergnügen heißt, und Sie können daher denken, wie sehr
es auf mich einwirken mußte.
Ich war dazu doppelt glücklich, aber auch doppelt mäßig
in jedem Genusse, da am Vorabend des Festes meine edle
Freundin mir mit ihrem Herzen auch ihre Hand zugesagt
hatte. Das geheimnißvolle Glück, das eine ausgesprochene
edle Liebe begleitet, deren Höhe von keinem Aber
geahndet wurde, goß einen besondern Reiz über unsere
beiderseitige Haltung und Wesen, die von unsern nähern
Freunden nicht übersehen wurde.
Ich komme nochmal auf unsere Widersacher. Niemand
ahndet zwar, daß zwischen uns eine Erklärung
stattgefunden und wir uns Beide vollkommen verstehen;
allein Jeder bemerkt leicht unsere gegenseitige Neigung,
und da ich in allen öffentlichen Orten ihr den Arm gebe,
bei Tisch ihr immer zur Seite bin, jeden ersten Tanz mit ihr
tanze, mich par préférence mit ihr unterhalte, sie
beständig nach Hause führe, so hat man natürlich unsere
gegenseitige Neigung nicht übersehen können, davon
abgesehen, daß da, wo ich nicht bin, sie mich hebt oder
nöthigenfalls vertheidigt, wie ich schon oben gedacht
habe.

Die weitere günstige Entwickelung der


Verlobungsangelegenheit schildert Brockhaus in
einem fernern Briefe an Bornträger aus Leipzig vom
21. September:
Ich bin seit meinem vorigen Briefe ein paar mal in Leipzig
gewesen, wo ich auch jetzt mich schon wieder seit acht
Tagen befinde, um mehrere Expeditionen zu
beschleunigen und vieles Andere zu reguliren, da die
Niederlage muß geräumt werden und hundert andere
Dinge zu thun sind.
Seit jenem Briefe hat sich in den dort geschilderten
Verhältnissen viel geändert und zum Guten, sodaß ich
hoffen darf, es werde sich Alles schön und edel lösen.
Ludwigs kamen den 4. September zurück. Ich war in
Leipzig und kam erst den 8. wieder nach Altenburg. Meine
Freundin hatte sich ihnen gleich erklärt und mit der
entschiedensten Energie sich ausgesprochen, daß nichts
sie zurückhalten würde, ihr Leben mit dem meinigen zu
vereinigen, wenn meine bürgerlichen Verhältnisse sich
ordnen ließen. Ludwigs hatten es gut aufgenommen und
ihr allen Beistand zugesagt.

Zu den »Widersachern«, von denen er mehrfach


spricht, gehörte besonders der mit dem
Ludwig'schen Hause eng befreundete Bankier August
Reichenbach. Indessen gelang es Brockhaus und
Frau Ludwig, auch ihn zu gewinnen, ja er wurde ihm
bald ein treuer Freund, der ihn auch materiell durch
Credit bei seinen Verlagsunternehmungen
unterstützte.
Wie sehr Brockhaus seine künftige Schwägerin Frau
Ludwig verehrte, zeigt folgendes
Gratulationsschreiben, das er zu ihrem Geburtstage,
27. December 1811, an sie richtete:
Als ich vor einem Jahre der »Schönen und Guten« am
heutigen Tage ein Zeichen meiner Verehrung brachte, wie
wenig kannte ich da noch den Umfang Ihres herrlichen
Geistes, den Adel Ihrer Seele, die Tiefe Ihres Gemüths, die
Wärme Ihres Herzens, die, zusammen vereint, Sie zum
Stolze und zur Ersten Ihres Geschlechtes machen, und
Allen, die Ihnen nahen und die Ihnen angehören, der
sicherste Leitstern sind fürs eigene Streben. Sie werden
heute vielfach begrüßt werden, liebe Ludwig, und gewiß
von Vielen in Liebe und Treue und Wahrheit. Ich geselle
mich zu den Vielen, und in kunstloser Rede sage ich Ihnen
denn auch, Keinem wenigstens an Wahrheit, Treue und
Freundschaft nachstehend, wie sehr ich Sie verehre und
wie meine heißesten Wünsche für Ihr Glück, für Ihren
Seelenfrieden, für Ihr Wohlsein sich mit denen Ihrer ältern
und besten Freunde vereinigen! Seien Sie so glücklich, als
Sie verdienen es zu sein!
Wie fern stand ich Ihnen vor einem Jahre! Wie unglücklich
war ich damals! Vieles, wie Vieles hat sich in den schnell
verflossenen Monden geändert! Ich sehe für mich die
Morgenröthe eines neuen Glücks aufgehen, das um so
größern Reiz für mich haben wird, je näher ich Ihnen,
Verehrte, dadurch zu stehen komme!
Möge ich Sie zur nächsten Feier des heutigen Tags mit
einem Namen begrüßen dürfen, der für mich, außer dem
Herrlichen, was er an sich in sich faßt, die schönste
Lebensmusik sein wird.

Im Laufe dieses und des folgenden Jahres hatten


sich Brockhaus' geschäftliche Verhältnisse immer
mehr befestigt. Die Verlobung wurde jetzt
veröffentlicht und den Verwandten und Freunden
mitgetheilt. Von allen Seiten kamen herzliche
Glückwünsche; der kurze, aber treffende
Glückwunsch eines dortmunder Jugendfreundes,
Johannes Rappe, an Brockhaus lautete:
Dein Genie hat Dich durch so mancherlei Labyrinthe des
bürgerlichen Lebens gejagt und geführt, daß Du meinen
Glückwunsch zu Deinem frohen Lebensgenuß in ruhiger
Wirksamkeit für aufrichtig anerkennen und Deiner
praktischen Vernunft zur Ausführung anvertrauen wirst.
Auch sein Bruder Gottlieb schrieb sehr herzlich, und
die in Dortmund noch weilenden drei Kinder feierten
dort die Hochzeit ihres Vaters wol deshalb besonders
freudig, weil sie ihnen die Aussicht bot, wieder eine
Mutter zu bekommen und nunmehr bald in das
älterliche Haus zurückkehren zu können.
Die Hochzeit fand in Altenburg am 26. November
1812 statt, unter regster Theilnahme der neuen und
alten Freunde des Bräutigams, die sich in zahlreichen
ernsten und humoristischen Gedichten kundgab.
Mit Bedauern vermißte Brockhaus unter seinen
anwesenden Freunden den Professor Ersch aus
Halle. Derselbe war im September bei ihm zu Besuch
gewesen und hatte ihm dann geschrieben:
Immer wird die Erinnerung meines Aufenthalts in
Altenburg an die erfreulichsten meines Lebens sich
anreihen; immer werde ich mit frohem Gefühle der
Stunden denken, in welchen ich Bekanntschaften mit
guten Menschen erneuerte und stiftete.

Jetzt durch Krankheit abgehalten, an der Hochzeit


theilzunehmen, schrieb er an Brockhaus aus Halle
vom 21. December:
Wahrlich, Sie hätten nicht nöthig gehabt, durch Ihre
Nachrichten von Ihrer frohen Hochzeit und den
Feierlichkeiten, mit welchen Ihre Freunde sie ausstatteten,
meine Trauer über die Entbehrung dieser Freuden zu
schärfen, und doch waren sie mir ungemein lieb und
interessant, vorzüglich erfreuend aber die Bemerkung, daß
Sie und Ihre gute Jeannette in Altenburg so viele Freunde
haben. Wer, wie ich, den höchsten Lebensgenuß in dem
Besitz von Freunden findet, weiß dies Gut zu würdigen.

Die nächsten vier Jahre, 1813-1816, verbrachte


Brockhaus meist in Altenburg, im ruhigen Genusse
seiner neuen Häuslichkeit, aber auch in
angestrengter Thätigkeit für den Wiederaufbau
seines Geschäfts und unter lebhafter Theilnahme an
den großen Ereignissen dieser Zeit. Außer häufigen
Fahrten nach Leipzig machte er nur im Sommer 1814
in Erbschaftsangelegenheiten seiner Frau eine
dreimonatliche Reise nach Stuttgart, Augsburg und
München, von wo er über Strasburg, Frankfurt a. M.
und Braunschweig zurückkehrte, und kleinere
Ausflüge nach Dresden, Weimar, Dessau, Wittenberg,
Berlin.
Von seinen Kindern hatte er Auguste und Hermann
im April 1814 von Dortmund nach Altenburg kommen
lassen, während Heinrich erst im Mai 1816 folgte und
die jüngste Tochter, Sophie, endlich im August 1817
von ihrem ältesten Bruder Friedrich aus Amsterdam
abgeholt und nach Altenburg gebracht wurde.
Friedrich war im Herbst 1813 zu dem Pastor
Schlosser in Großzschocher bei Leipzig gekommen,
wo er mit andern Knaben zusammen erzogen und
unterrichtet wurde; zu Neujahr 1816 nahm ihn auf
Wunsch seines Vaters der mit diesem befreundete
und schon seit der amsterdamer Zeit einen großen
Theil seiner Verlagswerke druckende Buchhändler
und Buchdrucker Hans Friedrich Vieweg in
Braunschweig zu sich in die Lehre; er sollte hier
gleichzeitig mit Vieweg's fast gleichaltrigem Sohne
Eduard die Buchdruckerkunst erlernen, weil sein
Vater die Absicht hatte, mit dem immer größere
Ausdehnung erlangenden Verlagsgeschäfte eine
Druckerei zu errichten. Der jüngste Sohn Hermann
erhielt mit den Ludwig'schen Kindern zusammen
Privatunterricht und kam später gleich seinem ältern
Bruder Heinrich, der an diesem Unterrichte auch mit
theilnahm, in die Erziehungsanstalt zu
Wackerbarthsruhe bei Dresden Die älteste Tochter
Auguste wurde im Januar 1815 in eine Pension nach
Dresden gebracht und war dort bis zur
Uebersiedelung ihres Vaters nach Leipzig; die zweite
Tochter, Karoline, blieb in Altenburg.
Von seiner Frau wurden ihm in dieser Zeit zwei
Kinder geboren, Alexander und Luise, die aber bald
wieder starben, ersterer am 20. August 1814,
letztere am 4. August 1818. Später wurden ihm noch
zwei Töchter geboren: Johanne Wilhelmine am 29.
December 1817 noch in Altenburg und Marie Ottilie
am 18. Mai 1821 in Leipzig.
Brockhaus' langjähriger vertrauter Gehülfe und treuer
Freund Bornträger war nach dem Verkaufe des
amsterdamer Geschäfts noch bis zum Frühjahre 1812
in Amsterdam geblieben, um die von dem Käufer,
Johannes Müller, nicht mit übernommenen
Außenstände einzuziehen und alle sonstigen
Verhältnisse daselbst zu regeln. Als ihm dies
gelungen war und er am 4. März 1812 den schon
früher erwähnten Vertrag mit dem amsterdamer
Buchhändler Sülpke abgeschlossen, schrieb ihm
Brockhaus offen: er könne ihm augenblicklich keine
feste Stellung in Altenburg anbieten, da seine
Verhältnisse noch zu wenig consolidirt seien, und
rathe ihm seiner selbst wegen eine andere Condition
anzunehmen, zu deren Erlangung er ihm gern
behülflich sein werde; für alle Fälle sei ihm in seinem
Hause ein Asyl gesichert. Durch diese Mittheilung
und wol auch durch manche Vorwürfe verletzt, die
ihm während der allerdings sehr schwierigen Zeit
seiner Geschäftsführung gemacht worden waren,
kündigte Bornträger und nahm eine untergeordnete
Stellung bei dem Buchhändler Tasché in Gießen an.
Er schrieb aber bald darauf selbst an Brockhaus, daß
er seinen Entschluß bereue, und in spätern Jahren,
bei Bornträger's regelmäßigem Besuche der leipziger
Messe, glichen sich alle Differenzen zwischen ihnen
vollständig aus. Bornträger rühmt selbst in einem
spätern Briefe, er habe sich der Freundschaft seines
frühern Principals bis zu dessen Tode zu erfreuen
gehabt.
In Gießen blieb Bornträger bis Anfang 1815, ging
dann nach Berlin zu Amelang und errichtete 1818 in
Gemeinschaft mit seinem jüngern Bruder Ludwig in
Königsberg unter der Firma Gebrüder Bornträger
eine Sortimentsbuchhandlung, mit der bald auch
Verlagsbuchhandel vereinigt wurde. Diese
Buchhandlung leitete er erst mit seinem Bruder, dann
nach dessen Tode (1843) allein bis zu seinem am 6.
März 1866 in hohem Alter (er war am 17. September
1787 zu Osterode am Harz geboren) erfolgten Tode
und wußte seiner (noch jetzt unter einem andern
Besitzer in Berlin fortblühenden)
Verlagsbuchhandlung Ansehen zu verschaffen; sein
Sortimentsgeschäft war schon 1842 an Tag & Koch
verkauft worden. Auch persönlich genoß er hohe
Achtung bei seinen Mitbürgern, die ihn 1843 zum
Stadtrath wählten.
Die Verdienste, die sich Bornträger um Brockhaus als
treuer Freund und Berather in schwieriger Zeit
erworben, werden auch von dessen Nachkommen
vollkommen gewürdigt und sein Andenken wird bei
ihnen stets in Ehren gehalten werden.
Bornträger's Nachfolger als Brockhaus' vertrauter
Gehülfe und bald in noch höherm Grade wie dieser
als Freund des Hauses wurde Karl Ferdinand
Bochmann, der am 10. Juli 1813 in das Geschäft
eintrat. Am 11. Februar 1788 zu Thurm bei Glauchau
geboren, hatte er in der Buchhandlung des Magister
Sommer in Leipzig sechs Jahre lang den Buchhandel
erlernt und dann, von Wanderlust getrieben, im
August 1808 eine Gehülfenstelle in Amsterdam bei
dem Buchhändler Hesse angenommen. Bezeichnend
für die damaligen Verhältnisse ist es, daß Hesse mit
seinem neuen Gehülfen einen förmlichen Vertrag
abschloß, in dem sich dieser verpflichten mußte, sich
niemals in Amsterdam zu etabliren, ja selbst »nie mit
den Principalen der andern zwei dortigen deutschen
Buchhandlungen und mit deren Leuten sich
einzulassen und allen Umgang mit denselben zu
vermeiden, ansonsten Er augenblickliche Entlassung
seiner Condition zu erwarten hat«. Trotzdem war er
in Amsterdam mit Brockhaus bekannt geworden. Als
Hesse im Sommer 1813 seine amsterdamer
Buchhandlung aufgab und nach Paris zog, nahm
Bochmann die ihm jetzt durch Vermittelung seiner an
Dr. Bernhardi in Altenburg verheiratheten Schwester
angebotene Stelle bei Brockhaus um so lieber an, als
es ihm bei den aufgeregten politischen Verhältnissen
Hollands in Amsterdam nicht mehr gefiel und er sich
nach der Heimat sehnte. War er doch 1809 sogar
gezwungen worden, in die amsterdamer Bürgerwehr
(Schutterij) einzutreten. So ergriff er am 11. Juni
1813 den Wanderstab und legte die Reise nach
Altenburg, wo er am 26. Juni eintraf, zu Fuße zurück.
Ueber seine Wanderung wie über die nächste so
ereignißreiche kriegerische Periode führte er ein
Tagebuch, das manches Interessante enthält. Er
gewann bald Brockhaus' vollständiges Vertrauen und
war schon während der altenburger Zeit dessen
Hauptstütze im Geschäft.

Außer Bochmann hatte Brockhaus noch zwei Männer


an sich gezogen, die ihn bei seiner literarischen und
redactionellen Thätigkeit unterstützten, während
Bochmann das rein Buchhändlerische besorgte: Dr.
Ludwig Hain, der im August 1812 eintrat, um ihn
zunächst bei der Redaction des »Conversations-
Lexikon«, später auch bei der Herausgabe der
»Deutschen Blätter« zu unterstützen, und bis 1820
bei ihm blieb, und Dr. Sievers, der im Herbst 1813 zu
Hain's Unterstützung kam, seine Stellung aber schon
1815 wieder aufgab.
Während dieser Zeit vollzog sich auch die
Umänderung der bisherigen Firma des Geschäfts
»Kunst- und Industrie-Comptoir« in die seitdem
beibehaltene Firma: »F. A. Brockhaus.« Und zwar
erfolgte diese Umänderung in ganz formloser Weise,
da man überhaupt auf alle solche Dinge damals
wenig Gewicht legte.
Wie schon früher erwähnt, gebrauchte Brockhaus
seit Aufgabe des amsterdamer Geschäfts die Firma
desselben auch in Altenburg noch fort, nur mit einem
Zusatz, indem er »Kunst- und Industrie-Comptoir von
Amsterdam« firmirte und auf den Büchertiteln bald
Altenburg, bald Leipzig, bald beide Städte als
Verlagsort nannte. In einem vom 15. Januar 1814,
und noch dazu nicht aus Altenburg, sondern aus
Leipzig (wo sich Brockhaus damals zufällig befand),
datirten Circulare über Rechnungsverhältnisse finden
sich am Schluß ganz beiläufig folgende Zeilen:
Noch bemerken wir Ihnen, daß wir von jetzt an blos nach
dem Eigenthümer unserer Handlung mit F. A.
B r o c k h a u s firmiren werden.

Diese Firmenzeichnung findet sich seitdem auf allen


seinen Verlagsartikeln, im Anfang noch abwechselnd
mit »Altenburg« oder »Leipzig« oder beiden Städten
als Verlagsorten, seit 1817 meist und seit 1819
ausschließlich mit dem Verlagsort »Leipzig«.
2.
Neue Verlagsthätigkeit.

In Altenburg entfaltete Brockhaus, sobald er die


Verwickelungen aus der amsterdamer Periode zum
Abschluß gebracht, gleich eine überaus rege und
umfassende Thätigkeit. Mit neuer Kraft und mit
gewohnter Energie gelang es ihm, von dem rasch
wiederkehrenden Vertrauen der Buchhändlerwelt
gehoben und von seinen neugewonnenen Freunden
in Altenburg moralisch und materiell unterstützt, sein
Verlagsgeschäft bald zu größerer Bedeutung zu
bringen, als es in Amsterdam gehabt, und dadurch
auch seine äußere Lage wieder zu einer günstigen zu
gestalten.

Seine Verlagsthätigkeit in dieser altenburger Periode


erstreckte sich besonders nach drei Richtungen hin.
Die eine umfaßt seine Thätigkeit auf politisch-
publicistischem Gebiete während der ereignißreichen
und für Deutschland so hochbedeutsamen Jahre
1813-1815. Die zweite Hauptthätigkeit betrifft das
»Conversations-Lexikon«, das er wesentlich in diesen
Jahren zu dem gestaltete, was es für ihn und für die
deutsche Literatur geworden ist. Die dritte Seite
Welcome to our website – the ideal destination for book lovers and
knowledge seekers. With a mission to inspire endlessly, we offer a
vast collection of books, ranging from classic literary works to
specialized publications, self-development books, and children's
literature. Each book is a new journey of discovery, expanding
knowledge and enriching the soul of the reade

Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.

Let us accompany you on the journey of exploring knowledge and


personal growth!

textbookfull.com

You might also like