Fundamentals of RCS Prediction Methodology Using Parallelized Numerical Electromagnetics Code NEC and Finite Element Pre Processor Vineetha Joy
Fundamentals of RCS Prediction Methodology Using Parallelized Numerical Electromagnetics Code NEC and Finite Element Pre Processor Vineetha Joy
com
https://textbookfull.com/product/reimagining-philosophy-and-
technology-reinventing-ihde-glen-miller/
textbookfull.com
The Empire of Disgust: Prejudice, Discrimination, and
Policy in India and the US 1st Edition Zoya Hasan
https://textbookfull.com/product/the-empire-of-disgust-prejudice-
discrimination-and-policy-in-india-and-the-us-1st-edition-zoya-hasan/
textbookfull.com
https://textbookfull.com/product/the-coconut-palm-cocos-nucifera-l-
research-and-development-perspectives-v-krishnakumar/
textbookfull.com
https://textbookfull.com/product/handbook-of-sustainability-in-
additive-manufacturing-volume-1-1st-edition-subramanian-senthilkannan-
muthu/
textbookfull.com
https://textbookfull.com/product/neurorheumatology-a-comprehenisve-
guide-to-immune-mediated-disorders-of-the-nervous-system-tracey-a-cho/
textbookfull.com
https://textbookfull.com/product/the-making-of-islamic-heritage-
muslim-pasts-and-heritage-presents-1st-edition-trinidad-rico-eds/
textbookfull.com
Puppy Love (Her Secret Menagerie #1) 1st Edition Katelyn
Beckett
https://textbookfull.com/product/puppy-love-her-secret-
menagerie-1-1st-edition-katelyn-beckett/
textbookfull.com
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
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.
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
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.
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
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.
xv
xvi List of Figures
xvii
Chapter 1
Introduction
© 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
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.
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 ) 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.
The integral representation for the magnetic field due to a surface current distribution
−
→
J s on a perfectly conducting body is given by,
Visit https://textbookfull.com
now to explore a rich
collection of eBooks, textbook
and enjoy exciting offers!
2.1 Formulation of Integral Equations (IE) 7
−
→ −→
H (
r ) = (1/4π ) J s (r ) · ∇ g r, r dA (2.4)
S
−
→ −
→
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 ) × tˆ2 (
tˆ1 ( r0 ) = n̂(
r0 ) (2.7)
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
−
→
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 ;
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).
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)
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.
N
β j wi , L f j = wi , u ; i = 1, 2, . . . , N (2.16)
j=1
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,
−
→ 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.
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,
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
Random documents with unrelated
content Scribd suggests to you:
Motolinia, Torquemada, Ixtlilxochitl, and Gomara agree in placing
the two oratories or shrines on the extreme eastern edge of the
platform, so that there was only just room for a man to pass round
them on the east side. The two oratories were separate one from the
other, each being enclosed within its own walls with a doorway
towards the west. The oratory of Huitzilopochtli was the larger of the
two and stood to the south. The oratory of Tlaloc stood to the north.
No measurements are given of the area covered by these two
oratories, but there is no suggestion that they were large buildings[7]
except in height. The roof and probably the upper stages were made
of wood (Torquemada), and we know that they were burnt during the
siege.
Height:—
Recent Excavations.
They build a square tower one hundred and fifty paces, or rather
more, in length, and one hundred and fifteen or one hundred in
breadth. The foundation of this building is solid; when it reaches the
height of two men, a passage is left two paces wide on three sides,
and on one of the long sides steps are made until the height of two
more men is reached, and the edifice is throughout solidly built of
masonry. Here, again, on three sides they leave the passage two
paces wide, and on the other side they build the steps, and in this
way it rises to such a height that the steps total one hundred and
twenty or one hundred and thirty.
There is a fair-sized plaza on the top and from the middle [of it]
arise two other towers which reach the height of ten or twelve men’s
bodies and these have windows above. Within these tall towers stand
the Idols in regular order and well adorned, and the whole house
highly decorated. No one but their high priest was allowed to enter
where the principal God was kept, and this god had distinct names in
different provinces; for in the great city of Mexico he was called
Horchilobos (Huitzilopochtli), and in another city named Chuennila
(Cholula) he was called Quecadquaal (Quetzalcoatl), and so on in the
others.
Whenever they celebrated the festivals of their Idols, they
sacrificed many men and women, boys and girls; and when they
suffered some privation, such as drought or excess of rain, or found
themselves hard pressed by their enemies, or suffered any other
calamity, then they made these sacrifices in the following manner....
The court of this Temple was very large, almost two hundred
fathoms square; it was all paved, and had within it many buildings
and towers. Some of these were more lofty than others, and each one
of them was dedicated to a God.
The principal tower of all was in the middle and was higher than
the others, and was dedicated to the God Vitzilopuchtli
Tlacavepancuexcotzin.
This tower was divided in the upper part, so that it looked like two,
and had two chapels or altars on the top, each one covered by its
dome (chapitel) and each one of them had on the summit its
particular badges or devices. In the principal one of them was the
statue of Vitzilopuchtli, also called Ilhuicatlxoxouhqui, and in the
other the image of the God Tlaloc. Before each one of these was a
round stone like a chopping-block, which they call Texcatl, where
they killed those whom they sacrificed in honour of that God, and
from the stone towards the ground below was a pool of blood from
those killed on it; and so it was on all the other towers; these faced
the West, and one ascended by very narrow straight steps to all these
towers.
(Sahagun mentions seventy-eight edifices in connection with the
Great Temple, but it is almost certain that these were not all within
the Temple enclosure.)
Sahagun, Hist. de la Conquista, Book 12, Ch. XXII.
They [the Mexicans] ascended a Cu, the one that was nearest to the
royal houses [i. e. of Axayacatl], and they carried up there two stout
beams so as to hurl them from that place on to the royal houses and
beat them down so as to force an entry. When the Spaniards
observed this they promptly ascended the Cu in regular formation,
carrying their muskets and crossbows, and they began the ascent
very slowly, and shot with their crossbows and muskets at those
above them, a musketeer accompanying each file and then a soldier
with sword and shield, and then a halberdier: in this order they
continued to ascend the Cu, and those above hurled the timbers
down the steps, but they did no damage to the Spaniards, who
reached the summit of the Cu and began to wound and kill those who
were stationed on the top, and many of them flung themselves down
from the Cu: finally, all those [Mexicans] who had ascended the Cu
perished.
Hernando Cortés, 2nd Letter. (The attack on the Great Teocalli.)
We fought from morning until noon, when we returned with the
utmost sadness to our fortress. On account of this they [the enemy]
gained such courage that they came almost up to the doors, and they
took possession of the great Mosque[15], and about five hundred
Indians, who appeared to me to be persons of distinction, ascended
the principal and most lofty tower, and took up there a great store of
bread and water and other things to eat, and nearly all of them had
very long lances with flint heads, broader than ours and no less
sharp.
From that position they did much damage to the people in the fort,
for it was very close to it. The Spaniards attacked this said tower two
or three times and endeavoured to ascend it, but it was very lofty,
and the ascent was steep, for it had more than one hundred steps,
and as those on the top were well supplied with stones and other
arms, and were protected because we were unable to occupy the
other terraces, every time the Spaniards began the ascent they were
rolled back again and many were wounded. When those of the enemy
who held other positions saw this, they were so greatly encouraged
they came after us up to the fort without any fear.
Then I (seeing that if they could continue to hold the tower, in
addition to the great damage they could do us from it, that it would
encourage them to attack us) set out from the fort, although maimed
in my left hand by a wound that was given me on the first day, and
lashing the shield to my arm, I went to the tower accompanied by
some Spaniards and had the base of it surrounded, for this was easily
done, although those surrounding it had no easy time, for on all sides
they were fighting with the enemy who came in great numbers to the
assistance of their comrades. I then began to ascend the stairway of
the said tower with some Spaniards supporting me, and although the
enemy resisted our ascent very stubbornly, so much so that they
flung down three or four Spaniards, with the aid of God and his
Glorious Mother (for whose habitation that tower had been chosen
and her image placed in it), we ascended the said tower and reaching
the summit we fought them so resolutely that they were forced to
jump down to some terraces about a pace in width which ran round
the tower. Of these the said tower had three or four, thrice a man’s
height from one [terrace] to the other.
Some fell down the whole distance [to the ground], and in addition
to the hurt they received from the fall, the Spaniards below who
surrounded the tower put them to death. Those who remained on the
terraces fought thence very stoutly, and it took us more than three
hours to kill them all, so that all died and none escaped ... and I set
fire to the tower and to the others which there were in the Mosque.
Juan de Torquemada, Monarchia Indiana, Vol. II. Book 8, Ch. XI.
p. 144. [Giving a description of the Great Temple.]
This Temple was rebuilt and added to a second time; and was so
large and of such great extent, that it was more than a crossbow-shot
square.
It was all enclosed in masonry of well squared stone.
There were in the square four gateways which opened to the four
principal streets, three of them by which the city was approached
along the causeways from the land, [the fourth] on the east in the
direction of the lake whence the City was entered by water.
In the middle of this enormous square was the Temple which was
like a quadrangular tower (as we have already stated) built of
masonry, large and massive.
This Temple (not counting the square within which it was built)
measured three hundred and sixty feet from corner to corner, and
was pyramidal in form and make, for the higher one ascended the
narrower became the edifice, the contractions being made at
intervals so as to embellish it.
On the top, where there was a pavement and small plaza rather
more than seventy feet wide, two very large altars had been built, one
apart from the other, set almost at the edge or border of the tower on
the east side, so that there was only just sufficient ground and space
for a man to walk [on the east side] without danger of falling down
from the building.
These altars were five palms in height with their walls inlaid with
stone, all painted with figures according to the whim and taste of him
who ordered the painting to be done. Above the altars were the
chapels roofed with very well dressed and carved wood.
Each of these chapels had three stories one above the other, and
each story or stage was of great height, so that each one of them [of
the chapels] if set on the ground (not on that tower, but on the
ground level whence the edifice sprang) would have made a very lofty
and sumptuous building, and for this reason the whole fabric of the
Temple was so lofty that its height compelled admiration. To behold,
from the summit of this temple, the city and its surroundings, with
the lakes and all the towns and cities that were built in it and on its
banks, was a matter of great pleasure and contentment.
On the West side this building had no stages [contractions], but
steps by which one ascended to the level of the chapels, and the said
steps had a rise of one foot or more. The steps, or stairs, of this
famous temple numbered one hundred and thirteen, and all were of
very well dressed stone.
From the last step at the summit of this Temple to the Altars and
entrance to the Chapels was a considerable space of ground, so that
the priests and ministers of the Idols could carry out their functions
unimpeded and thoroughly.
On each of the two altars stood an Idol of great bulk, each one
representing the greatest God they possessed, which was
Huitzilupuchtli or by his other name Mexitli.
Near and around this Great Temple there were more than forty
lesser ones, each one of them dedicated and erected to a God, and its
tower and shape narrowed up to the floor on which the Chapel and
altar began to arise, and it was not as large as the Great Temple, nor
did it approach it by far in size, and all these lesser Temples and
towers were associated with the Great Temple and tower which there
was in this City.
The difference between the Great Temple and the lesser ones was
not in the form and structure, for all were the same, but they differed
in site and position [orientation], for the Great Temple had its back
to the East, which is the practice the large temples ought to follow, as
we have noticed that the ancients assert, and their steps and
entrance to the West (as we are accustomed to place many of our
Christian Churches), so that they paid reverence in the direction of
the sun as it rose, the smaller temples looked in the other direction
towards the East and to other parts of the heaven [that is to] the
North and South.
PLATE A.
Part of the City of Mexico from a modern Map.
TRACING A1.
After J. García Icazbalceta.
TRACING A2.
Suggested site of the Great Teocalli and enclosure.
PLATE B.
Suggested plan and section of the Great Teocalli.
PLATE C.
Plan by Padre Sahagun, after Dr. E. Seler.
PLATE D.
View of the Great Teocalli, after Ixtlilxochitl (‘Codice Goupil’).
PLATE E.
View of the Great Teocalli and enclosure, from ‘The Chronicle of
Mexico,’ 1576. (Manuscript in British Museum, No. 31219.
Additional.)
PRINTED BY TAYLOR AND FRANCIS,
RED LION COURT, FLEET STREET.
Tracing A1