Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correspondence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Analysis of Voyager Observed High-Energy Electron Fluxes in the Heliosheath Using MHD Simulations

The Voyager spacecraft (V1 and V2) observed electrons of 6-14 MeV in the heliosheath which showed several incidences of flux variation relative to a background of gradually increasing flux with distance from the Sun. The increasing flux of background electrons is thought to result from inward radial diffusion. We compare the temporal electron flux variation with dynamical phenomena in the heliosheath that are obtained from our MHD simulations. Because our simulation is based on V2 observed plasma data before V2 crossed the termination shock, this analysis is effective up to late 2008, i.e., about a year after the V2-crossing, during which disturbances, driven prior to the crossing time, survived in the heliosheath. Several electron flux variations correspond to times directly associated with interplanetary shock events. One noteworthy example corresponds to various times associated with the March 2006 interplanetary shock, these being the collision with the termination shock, the passage past the V1 spacecraft, and the collision with the region near the heliopause, as identified by W.R. Webber et al. for proton/helium of 7-200 MeV. Our simulations indicate that all other electron flux variations, except one, correspond well to the times when a shock-driven magneto-sonic pulse and its reflection in the heliosheath either passed across V1/V2, or collided with the termination shock or with the plasma sheet near the heliopause. This result suggests that variation in the electron flux should be due to either direct or indirect effects of magnetosonic pulses in the heliosheath driven by interplanetary shocks

Washimi, Haruichi↗

Propagation of Convective Complexes Observed by TRMM in the Eastern Tropical Atlantic

Precipitation maxima during the West African summer monsoon propagate generally westward in tandem with African easterly waves. A heretofore unreported, repeating pattern of northward drift of precipitation maxima is detected on Tropical Rainfall Measurement Mission (TRMM satellite) time-latitude distributions of daily accumulations over the eastern tropical Atlantic. Corresponding 3-hourly TRMM accumulations show that the northward drifting envelopes of precipitation during August 2006 are often comprised of individual swaths propagating towards the southwest, presumably as mesoscale squall lines. The implied northward drift on the time-latitude distribution is a component of a resultant northwestward movement. The study examines the entire available record of TRMM precipitation observations, 1998-2010, to summarize TRMM maxima propagation over the eastern tropical Atlantic. Meridional displacements of precipitation maxima are most prevalent in June-September 2006, occurring less frequently during other summers. An investigation of geopotential and circulation fields, limited to two case studies, suggests mechanisms to explain some of the observed propagation of TRMM maxima. In one event, northward drift of the precipitation envelope is consistent with the corresponding displacement of the intertropical convergence zone trough, although the southwest propagation of individual mesoscale convection maxima does not correspond to any synoptic feature on reanalysis circulation or reanalysis downscaled by a regional model. One speculation is that southwestward propagation of precipitation maxima could be caused by regeneration of convection at outflow boundaries of mature thunderstorms.

Druyan, Leonard M.↗

In Situ Detection of Strong Langmuir Turbulence Processes in Solar Type III Radio Bursts

The high time resolution observations obtained by the WAVES experiment of the STEREO spacecraft in solar type III radio bursts show that Langmuir waves often occur as intense localized wave packets. These wave packets are characterized by short durations of only a few ms and peak intensities, which well exceed the supersonic modulational instability (MI) thresholds. These timescales and peak intensities satisfy the criterion of the solitons collapsed to spatial scales of a few hundred Debye lengths. The spectra of these wave packets consist of primary spectral peaks corresponding to beam-resonant Langmuir waves, two or more sidebands corresponding to down-shifted and up-shifted daughter Langmuir waves, and low frequency enhancements below a few hundred Hz corresponding to daughter ion sound waves. The frequencies and wave numbers of these spectral components satisfy the resonance conditions of the modulational instability (MI). Moreover, the tricoherences, computed using trispectral analysis techniques show that these spectral components are coupled to each other with a high degree of coherency as expected of the MI type of four wave interactions. The high intensities, short scale lengths, sideband spectral structures and low frequency spectral enhancements and, high levels of tricoherences amongst the spectral components of these wave packets provide unambiguous evidence for the supersonic MI and related strong turbulence processes in type III radio bursts. The implication of these observations include: (1) the MI and related strong turbulence processes often occur in type III source regions, (2) the strong turbulence processes probably play very important roles in beam stabilization as well as conversion of Langmuir waves into escaping radiation at the fundamental and second harmonic of the electron plasma frequency, fpe, and (3) the Langmuir collapse probably follows the route of MI in type III radio bursts.

Langmuir↗

Field Programmable Gate Array Apparatus, Method, and Computer Program

An apparatus is provided that includes a plurality of modules, a plurality of memory banks, and a multiplexor. Each module includes at least one agent that interfaces between a module and a memory bank. Each memory bank includes an arbiter that interfaces between the at least one agent of each module and the memory bank. The multiplexor is configured to assign data paths between the at least one agent of each module and a corresponding arbiter of each memory bank based on the assigned data path. The at least one agent of each module is configured to read data from the corresponding arbiter of the memory bank or write modified data to the corresponding arbiter of the memory bank.

Morfopoulos, Arin C.↗

Supplementary Spin Tests of a 1/20-Scale Model of the Bell XP-39E Airplane

As a result of the loss of the first XP-39E airplane due to an uncontrolled spin, and at the request of the War Department received in their letter to the NACA, March 31, 1942, additional tests have been made of the corresponding 1/20-scale spin-tunnel model. Tests were made for the normal full load condition and also for the loading condition corresponding to that of the airplane at the time of the crash, that Is, with guns, ammunition, and half of the fuel removed from the wings. The model tests were made with the tail arrangement corresponding to that of the airplane at the time of the crash, as well as with numerous modified tail and contour arrangements In an attempt to evolve satisfactory recovery characteristics. The tests were started on April 2, 1942, and completed in approximately two weeks, with the development of several alternate tail arrangements giving reasonably satisfactory recoveries for normal spinning control configuration. The tests were witnessed by Mr. Melvin Shorr, of the Materiel Center at Wright Field, and Mr. Paul C. Emmons, of the Bell Aircraft Corporation, to whom records of the test data were furnished. Colonel H. P. Swofford, of the Materiel Center, Wright Field, also witnessed part of the tests. A preliminary report of the results of the tests has been forwarded to the Materiel Command at Wright Field. A brief check of the parachute installation for the airplane was also made.

A I Neihouse↗

A Metric to Quantify Shared Visual Attention in Two-Person Teams

Introduction: Critical tasks in high-risk environments are often performed by teams, the members of which must work together efficiently. In some situations, the team members may have to work together to solve a particular problem, while in others it may be better for them to divide the work into separate tasks that can be completed in parallel. We hypothesize that these two team strategies can be differentiated on the basis of shared visual attention, measured by gaze tracking. 2) Methods: Gaze recordings were obtained for two-person flight crews flying a high-fidelity simulator (Gontar, Hoermann, 2014). Gaze was categorized with respect to 12 areas of interest (AOIs). We used these data to construct time series of 12 dimensional vectors, with each vector component representing one of the AOIs. At each time step, each vector component was set to 0, except for the one corresponding to the currently fixated AOI, which was set to 1. This time series could then be averaged in time, with the averaging window time (t) as a variable parameter. For example, when we average with a t of one minute, each vector component represents the proportion of time that the corresponding AOI was fixated within the corresponding one minute interval. We then computed the Pearson product-moment correlation coefficient between the gaze proportion vectors for each of the two crew members, at each point in time, resulting in a signal representing the time-varying correlation between gaze behaviors. We determined criteria for concluding correlated gaze behavior using two methods: first, a permutation test was applied to the subjects' data. When one crew member's gaze proportion vector is correlated with a random time sample from the other crewmember's data, a distribution of correlation values is obtained that differs markedly from the distribution obtained from temporally aligned samples. In addition to validating that the gaze tracker was functioning reasonably well, this also allows us to compute probabilities of coordinated behavior for each value of the correlation. As an alternative, we also tabulated distributions of correlation coefficients for synthetic data sets, in which the behavior was modeled as a first-order Markov process, and compared correlation distributions for identical processes with those for disparate processes, allowing us to choose criteria and estimate error rates. 3) Discussion: Our method of gaze correlation is able to measure shared visual attention, and can distinguish between activities involving different instruments. We plan to analyze whether pilots strategies of sharing visual attention can predict performance. Possible measurements of performance include expert ratings from instructors, fuel consumption, total task time, and failure rate. While developed for two-person crews, our approach can be applied to larger groups, using intra-class correlation coefficients instead of the Pearson product-moment correlation.

Gontar, Patrick↗

Propellant Feed System for Swirl-Coaxial Injection

A propellant feed system for swirl-coaxial injection of a liquid propellant includes a reservoir having a bottom plate and at least one tube originating in the bottom plate and extending therefrom. The tube has rectangular slits defined in and distributed tangentially and evenly about a portion of the tube that is disposed in the bottom plate. Drain holes are provided in the bottom plate and tunnels are defined in the bottom plate. Each tunnel fluidly couples one of the drain holes to a corresponding one of the rectangular slits. Each tunnel includes (i) a bend of at least 90.degree., and (ii) a straight portion leading to its corresponding rectangular slit wherein the straight portion is at least five times as long as a hydraulic diameter of the corresponding rectangular slit.

Reynolds, David Christopher↗

Motion Estimation Utilizing Range Detection-Enhanced Visual Odometry

A motion determination system is disclosed. The system may receive a first and a second camera image from a camera, the first camera image received earlier than the second camera image. The system may identify corresponding features in the first and second camera images. The system may receive range data comprising at least one of a first and a second range data from a range detection unit, corresponding to the first and second camera images, respectively. The system may determine first positions and the second positions of the corresponding features using the first camera image and the second camera image. The first positions or the second positions may be determined by also using the range data. The system may determine a change in position of the machine based on differences between the first and second positions, and a VO-based velocity of the machine based on the determined change in position.

Friend, Paul Russell↗

Small vs. Large Convective Cloud Objects from CERES Aqua Observations: Where are the Intraseasonal Variation Signals?

During inactive phases of Madden-Julian oscillation (MJO), there are plenty of deep but small convective systems and far fewer deep and large ones. During active phases of MJO, a manifestation of an increase in the occurrence of large and deep cloud clusters results from an amplification of large-scale motions by stronger convective heating. This study is designed to quantitatively examine the roles of small and large cloud clusters during the MJO life cycle. We analyze the cloud object data from Aqua CERES observations for tropical deep convective (DC) and cirrostratus (CS) cloud object types according to the real-time multivariate MJO index. The cloud object is a contiguous region of the earth with a single dominant cloud-system type. The size distributions, defined as the footprint numbers as a function of cloud object diameters, for particular MJO phases depart greatly from the combined (8-phase) distribution at large cloud-object diameters due to the reduced/increased numbers of cloud objects related to changes in the large-scale environments. The medium diameter corresponding to the combined distribution is determined and used to partition all cloud objects into "small" and "large" groups of a particular phase. The two groups corresponding to the combined distribution have nearly equal numbers of footprints. The medium diameters are 502 km for DC and 310 km for cirrostratus. The range of the variation between two extreme phases (typically, the most active and depressed phases) for the small group is 6-11% in terms of the numbers of cloud objects and the total footprint numbers. The corresponding range for the large group is 19-44%. In terms of the probability density functions of radiative and cloud physical properties, there are virtually no differences between the MJO phases for the small group, but there are significant differences for the large groups for both DC and CS types. These results suggest that the intreseasonal variation signals reside at the large cloud clusters while the small cloud clusters represent the background noises resulting from various types of the tropical waves with different wavenumbers and propagation directions/speeds.

Xu, Kuan-Man↗

Background Noises Versus Intraseasonal Variation Signals: Small vs. Large Convective Cloud Objects From CERES Aqua Observations

During inactive phases of Madden-Julian Oscillation (MJO), there are plenty of deep but small convective systems and far fewer deep and large ones. During active phases of MJO, a manifestation of an increase in the occurrence of large and deep cloud clusters results from an amplification of large-scale motions by stronger convective heating. This study is designed to quantitatively examine the roles of small and large cloud clusters during the MJO life cycle. We analyze the cloud object data from Aqua CERES (Clouds and the Earth's Radiant Energy System) observations between July 2006 and June 2010 for tropical deep convective (DC) and cirrostratus (CS) cloud object types according to the real-time multivariate MJO index, which assigns the tropics to one of the eight MJO phases each day. The cloud object is a contiguous region of the earth with a single dominant cloud-system type. The criteria for defining these cloud types are overcast footprints and cloud top pressures less than 400 hPa, but DC has higher cloud optical depths (=10) than those of CS (<10). The size distributions, defined as the footprint numbers as a function of cloud object diameters, for particular MJO phases depart greatly from the combined (8-phase) distribution at large cloud-object diameters due to the reduced/increased numbers of cloud objects related to changes in the large-scale environments. The medium diameter corresponding to the combined distribution is determined and used to partition all cloud objects into "small" and "large" groups of a particular phase. The two groups corresponding to the combined distribution have nearly equal numbers of footprints. The medium diameters are 502 km for DC and 310 km for cirrostratus. The range of the variation between two extreme phases (typically, the most active and depressed phases) for the small group is 6-11% in terms of the numbers of cloud objects and the total footprint numbers. The corresponding range for the large group is 19-44%. In terms of the probability density functions of radiative and cloud physical properties, there are virtually no differences between the MJO phases for the small group, but there are significant differences for the large groups for both DC and CS types. These results suggest that the intreseasonal variation signals reside at the large cloud clusters while the small cloud clusters represent the background noises resulting from various types of the tropical waves with different wavenumbers and propagation speeds/directions.

Xu, Kuan-Man↗

Modular Battery Controller

Some embodiments of the present invention describe a battery including a plurality of master-less controllers. Each controller is operatively connected to a corresponding cell in a string of cells, and each controller is configured to bypass a fraction of current around the corresponding cell when the corresponding cell has a greater charge than one or more other cells in the string of cells.

Button, Robert M↗

Where are the r-modes? Chandra Observations of Millisecond Pulsars

We present the results of Chandra observations of two non-accreting millisecond pulsars, PSRs J1640+2224(J1640) and J1709+2313 (J1709), with low inferred magnetic fields and spin-down rates in order to constrain their surface temperatures, obtain limits on the amplitude of unstable r-modes in them, and make comparisons with similar limits obtained for a sample of accreting low-mass X-ray binary (LMXB) neutron stars. We detect both pulsars in the X-ray band for the first time. They are faint, with inferred soft X-ray fluxes(0.3-3 keV) of approx. 6 x 10(exp -15) and 3 x 10( exp -15) erg/sq cm for J1640 and J1709, respectively. Spectral analysis assuming hydrogen atmosphere emission gives global effective temperature upper limits (90% confidence) of 3.3-4.3 x 10(exp 5) K for J1640 and 3.6-4.7 x 10(exp 5) K for J1709, where the low end of the range corresponds to canonical neutron stars (M = 1.4 Stellar Mass), and the upper end corresponds to higher-mass stars (M = 2.21 Stellar Mass). Under the assumption that r-mode heating provides the thermal support, we obtain dimensionless r-mode amplitude upper limits of 3.2-4.8 x 10(exp -8) and 1.8-2.8 x 10(exp -7) for J1640 and J1709, respectively, where again the low end of the range corresponds to lower-mass, canonical neutron stars (M =1.4 Stellar Mass). These limits are about an order of magnitude lower than those we derived previously for a sample of LMXBs, except for the accreting millisecond X-ray pulsar SAX J1808.43658, which has a comparable amplitude limit to J1640 and J1709.

oscillations↗

Asymmetric Dianhyrdrides and Methods of Making the Same

A method for preparing asymmetric 3,4'-(hexafluroisopropylidene) dipthalic anhydride (a-6FDA) and asymmetric 3,4'-(methylene)diphthalic anhydride (a-MDPA) from asymmetric 2,3,3',4'-benzophenone dianhydride (a-BTDA). First, a-BTDA is converted to its corresponding N-alkyl or N-phenyl bisimides by reacting with either methylamine or aniline. The sequential addition of two trifluoromethyl (CF.sub.3) groups to the carbonyl unit of a-BTDA produce the corresponding bisimides of asymmetric 6FDA, which is hydrolyzed to the asymmetric 6F-tetracarboxylic acid followed by cyclodehydration to obtain asymmetric 6F-dianhydride (a-6FDA). Separately, the carbonyl unit of bisimides of a-BTDA is reduced to a methylene unit to afford the corresponding bisimides of asymmetric 3,4'-(methylene)diphthalic anhydride, which is hydrolyzed to the asymmetric methylene tetracarboxylic acid followed by cyclodehydration to obtain 3,4'-(methylene)diphthalic anhydride (a-MDPA).

Chuang, Chun-Hua↗

Effect of Distributed Patch of Smooth Roughness Elements on Transition in a High-Speed Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of hypersonic vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, smooth-shaped roughness elements on the laminar boundary layer over a 7-degree half-angle, circular cone for flow conditions corresponding to a selected trajectory point from the ascent phase of the HIFiRE-1 flight experiment. For peak-to-valley roughness heights of 50 percent or less in comparison with the thickness of the unperturbed boundary layer, the computations converge to a stationary flow, suggesting that the flow is globally stable. Analysis of convective instabilities in the wake of the roughness patch indicates two dominant families of unstable disturbances, namely, a high frequency mode that corresponds to Mack mode waves modified by the wake and a lower frequency mode that corresponds to shear layer instabilities associated with the streaks in the roughness wake. Even though the peak growth rate of the later mode is more than 35 percent greater than the peak growth rates of the Mack modes, the latter modes achieve higher amplification ratios, and hence, are likely to dominate the onset of transition, which is estimated to occur slightly later than that in the unperturbed, i.e., smooth surface boundary layer. Additional computations are performed to investigate the effects of various roughness patch configurations on a Mach 3.5 flat plate boundary layer, to help guide an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at NASA Langley Research Center. In this case, the cumulative reinforcement of basic state distortion over the length of the roughness patch is predicted to yield a significantly earlier transition than that over a smooth plate or a plate with a shorter length roughness patch.

Meelan Choudhari↗

Multi-point Jig Twist Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft and High-Fidelity CFD Validation

This paper presents a multi-point jig twist optimization study of the Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft using an in-house developed aero-structural analysis solver VSPAERO coupled to BEAM3D. A vortex-lattice model of the TTBW model is developed, and a transonic small disturbance /integral boundary layer correction method is implemented in the VSPAERO model to account for transonic and viscous flow effects. A correction method for the wing-strut interference aerodynamics is developed and applied to the VSPAERO solver. Also, a structural dynamic finite-element model of the TTBW aircraft is developed. This finite-element model includes the geometric nonlinear effect due to the tension in the struts which causes a deflection-dependent nonlinear stiffness. The VSPAERO model is coupled to the corresponding finite-element model to provide a rapid aero-structural analysis. The VSPAERO model has been validated against wind tunnel test data and FUN3D CFD simulation data. Three different flight conditions corresponding to Mach 0.745 at 40,000 ft are selected for the TTBW aircraft jig twist optimization to reduce the drag coefficients at the design and off-design cruise lift coefficients. The preliminary optimization results show that the TTBW aircraft with the optimized jig twist achieves a drag reduction of about 9 counts, 7 counts, and8 counts corresponding to the lift coefficients 0.68, 0.73, and 0.78, respectively. A high-fidelity CFD solver FUN3D is used to validate the jig twist optimization.

Multi-point↗

Summer temperatures during the last glaciation (MIS 5c to MIS 3) inferred from a 50,000-year chironomid record from Füramoos, southern Germany

There is a sparsity of long, continuous palaeotemperature records for the last glacial period in central Europe, particularly for the interval corresponding to Marine Isotope Stages (MIS) 4 and 3. Here we present a new, ca. 50-thousand year (ka)-long chironomid record from Füramoos, southern Germany, covering the interval from MIS 5a to MIS 3 that we use to examine lake development and then to quantitatively reconstruct mean July air temperatures. Chironomid assemblages with high abundances of taxa such as Polypedilum nubeculosum-type, Microtendipes pedellus-type, Cladopelma lateralis-type and Dicrotendipes nervosus-type imply a shallow-lake setting for the majority of the examined interval, which is corroborated by other aquatic remains such as oribatid mites, Sialidae and Ceratopogonidae. Assemblages from the interval ca. 99 to 80 ka (in the region corresponding to the Brörup Interstadial, Stadial B and early Odderade Interstadial) are dominated by taxa such as Tanytarsus glabrescens-type and Tanytarsus mendax-type and indicate relatively warm temperatures. Assemblages from the interval covering ca. 80 to 54 ka (corresponding to the late Odderade, Stadial C, Dürnten Interstadial and Stadial D) are dominated by taxa such as Sergentia coracina-type and Tanytarsus lugens-type and are typical for cooler conditions. Reconstructed July temperatures for the early Würmian (Brörup to early Odderade; ca. 99–80 ka) are 13–14 °C. Values decline to <10 °C during the late Odderade and Stadial C (ca. 80–77 ka) around the MIS 5a/4 transition. This decrease is coeval with a pronounced decrease in Northern Hemisphere summer insolation. Values stay in the range of 9–11 °C during the Dürnten and Stadial D (ca. 54–74.5 ka) and increase again to 12.5 °C during the Bellamont 1 interstadial (ca. 54–46 ka). Reconstructed July temperatures track changes in arboreal pollen percentages at Füramoos and agree with a summer-temperature decrease during the early to mid-Würmian as reported by other palaeotemperature records from Europe and the North Atlantic. Our chironomid record from Füramoos provides valuable new insights into Würmian climate dynamics in Central Europe, and corroborates other temperature reconstructions from the early to mid-Würmian glacial period.

Chironomids↗

The Atacama Cosmology Telescope: A search for Planet 9

We use Atacama Cosmology Telescope (ACT) observations at 98 GHz (2015–2019), 150 GHz (2013–2019), and 229 GHz (2017–2019) to perform a blind shift-and-stack search for Planet 9. The search explores distances from 300 au to 2000 au and velocities up to 6\farcm3 per year, depending on the distance (r). For a 5 Earth-mass Planet 9 the detection limit varies from 325 au to 625 au, depending on the sky location. For a 10 Earth-mass planet the corresponding range is 425 au to 775 au. The predicted aphelion and most likely location of the planet corresponds to the shallower end of these ranges. The search covers the whole 18,000 square degrees of the ACT survey. No significant detections are found, which is used to place limits on the millimeter-wave flux density of Planet 9 over much of its orbit. Overall we eliminate roughly 17% and 9% of the parameter space for a 5 and 10 Earth-mass Planet 9, respectively. These bounds approach those of a recent INPOP19a ephemeris-based analysis, but do not exceed it. We also provide a list of the 10 strongest candidates from the search for possible follow-up. More generally, we exclude (at 95% confidence) the presence of an unknown solar system object within our survey area brighter than 4–12 mJy (depending on position) at 150 GHz with current distance 300 au < r < 600 au and heliocentric angular velocity 1\farcm5 per yr < v x (500au/r) < 2\farcs3 per yr, corresponding to low-to-moderate eccentricities. These limits worsen gradually beyond 600 au, reaching 5–15 mJy by 1500 au.

Sigurd Naess↗

A Living Catalog of Stream Interaction Regions in the Parker Solar Probe Era

Stream interaction regions (SIRs) and corotating interaction regions (CIRs) are important phenomena in heliospheric physics. These large-scale structures vary temporally and spatially, both in latitude and with radial distance. The additions of Parker Solar Probe(PSP) and Solar Orbiter have allowed for investigations into the radial evolution of these structures over a wide range of heliocentric distances for the first time since the Helios era. To better enable investigations of SIRs and CIRs within the inner heliosphere, we have developed a living catalog of SIR and CIR observations by Parker Solar Probe with corresponding observations by STEREO-A as well as ACE and Wind at 1 au. The methodology used for the identification of events and the generation of this catalog, the initial catalog of PSP observations spanning orbits one through five along with corresponding 1 au observations, and information on accessing the living catalog for future studies is described. This list of SIR and CIR events from PSP and corresponding observations from other heliophysics missions will enable case studies utilizing unique orbital arrangements, as well as aid in future statistical studies to further understand the properties and evolution of these structures.

Catalogs↗