SEARCH · Engineering Papers

Results for “Non-linear Dynamics”

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.

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRS2003

System Identification of X-33 Neural Network

Modern flight control research has improved spacecraft survivability as its goal. To this end we need to have a failure detection system on board. In case the spacecraft is performing imperfectly, reconfiguration of control is needed. For that purpose we need to have parameter identification of spacecraft dynamics. Parameter identification of a system is called system identification. We treat the system as a black box which receives some inputs that lead to some outputs. The question is: what kind of parameters for a particular black box can correlate the observed inputs and outputs? Can these parameters help us to predict the outputs for a new given set of inputs? This is the basic problem of system identification. The X33 was supposed to have the onboard capability of evaluating the current performance and if needed to take the corrective measures to adapt to desired performance. The X33 is comprised of both rocket and aircraft vehicle design characteristics and requires, in general, analytical methods for evaluating its flight performance. Its flight consists of four phases: ascent, transition, entry and TAEM (Terminal Area Energy Management). It spends about 200 seconds in ascent phase, reaching an altitude of about 180,000 feet and a speed of about 10 to 15 Mach. During the transition phase which lasts only about 30 seconds, its altitude may increase to about 190,000 feet but its speed is reduced to about 9 Mach. At the beginning of this phase, the Main Engine is Cut Off (MECO) and the control is reconfigured with the help of aerosurfaces (four elevons, two flaps and two rudders) and reaction control system (RCS). The entry phase brings down the altitude of X33 to about 90,000 feet and its speed to about Mach 3. It spends about 250 seconds in this phase. Main engine is still cut off and the vehicle is controlled by complex maneuvers of aerosurfaces. The last phase TAEM lasts for about 450 seconds and the altitude and speed, both are reduced to zero. The present attempt, as a start, focuses only on the entry phase. Since the main engine remains cut off in this phase, there is no thrust acting on the system. This considerably simplifies the equations of motion. We introduce another simplification by assuming the system to be linear after some non-linearities are removed analytically from our consideration. Under these assumptions, the problem could be solved by Classical Statistics by employing the least sum of squares approach. Instead we chose to use the Neural Network method. This method has many advantages. It is modern, more efficient, can be adapted to work even when the assumptions are diluted. In fact, Neural Networks try to model the human brain and are capable of pattern recognition.

Aggarwal, Shiv↗

NASA NTRS2003

Mesospheric Non-Migrating Tides Generated With Planetary Waves: Characteristics - 1

We discuss results from a modeling study with our Numerical Spectral Model (NSM) that specifically deals with the non-migrating tides generated in the mesosphere. The NSM extends from the ground to the thermosphere, incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GWs), and it describes the major dynamical features of the atmosphere including the wave driven equatorial oscillations (QBO and SAO), and the seasonal variations of tides and planetary waves. Accounting solely for the excitation sources of the solar migrating tides, the NSM generates through dynamical interactions also non-migrating tides in the mesosphere that are comparable in magnitude to those observed. Large non-migrating tides are produced in the diurnal and semi-diurnal oscillations for the zonal mean (m = 0) and in the semidiurnal oscillation for m = 1. In general, significant eastward and westward propagating tides are generated for all the zonal wave numbers m = 1 to 4. To identify the cause, the NSM is run without the solar heating for the zonal mean (m = 0), and the amplitudes of the resulting non-migrating tides are then negligibly small. In this case, the planetary waves are artificially suppressed, which are generated in the NSM through instabilities. This leads to the conclusion that the non-migrating tides are generated through non-linear interactions between planetary waves and migrating tides, as Forbes et al. and Talaat and Liberman had proposed. In an accompanying paper, we present results from numerical experiments, which indicate that gravity wave filtering contributes significantly to produce the non-linear coupling that is involved.

Mayr, H. G.↗

NASA NTRS2000

Kinetics of DSB rejoining and formation of simple chromosome exchange aberrations

PURPOSE: To investigate the role of kinetics in the processing of DNA double strand breaks (DSB), and the formation of simple chromosome exchange aberrations following X-ray exposures to mammalian cells based on an enzymatic approach. METHODS: Using computer simulations based on a biochemical approach, rate-equations that describe the processing of DSB through the formation of a DNA-enzyme complex were formulated. A second model that allows for competition between two processing pathways was also formulated. The formation of simple exchange aberrations was modelled as misrepair during the recombination of single DSB with undamaged DNA. Non-linear coupled differential equations corresponding to biochemical pathways were solved numerically by fitting to experimental data. RESULTS: When mediated by a DSB repair enzyme complex, the processing of single DSB showed a complex behaviour that gives the appearance of fast and slow components of rejoining. This is due to the time-delay caused by the action time of enzymes in biomolecular reactions. It is shown that the kinetic- and dose-responses of simple chromosome exchange aberrations are well described by a recombination model of DSB interacting with undamaged DNA when aberration formation increases with linear dose-dependence. Competition between two or more recombination processes is shown to lead to the formation of simple exchange aberrations with a dose-dependence similar to that of a linear quadratic model. CONCLUSIONS: Using a minimal number of assumptions, the kinetics and dose response observed experimentally for DSB rejoining and the formation of simple chromosome exchange aberrations are shown to be consistent with kinetic models based on enzymatic reaction approaches. A non-linear dose response for simple exchange aberrations is possible in a model of recombination of DNA containing a DSB with undamaged DNA when two or more pathways compete for DSB repair.

DNA Damage↗

NASA NTRS2000

Vascular mechanics of the coronary artery

This paper describes our research into the vascular mechanics of the coronary artery and plaque. The three sections describe the determination of arterial mechanical properties using intravascular ultrasound (IVUS), a constitutive relation for the arterial wall, and finite element method (FEM) models of the arterial wall and atheroma. METHODS: Inflation testing of porcine left anterior descending coronary arteries was conducted. The changes in the vessel geometry were monitored using IVUS, and intracoronary pressure was recorded using a pressure transducer. The creep and quasistatic stress/strain responses were determined. A Standard Linear Solid (SLS) was modified to reproduce the non-linear elastic behavior of the arterial wall. This Standard Non-linear Solid (SNS) was implemented into an axisymetric thick-walled cylinder numerical model. Finite element analysis models were created for five age groups and four levels of stenosis using the Pathobiological Determinants of Atherosclerosis Youth (PDAY) database. RESULTS: The arteries exhibited non-linear elastic behavior. The total tissue creep strain was epsilon creep = 0.082 +/- 0.018 mm/mm. The numerical model could reproduce both the non-linearity of the porcine data and time dependent behavior of the arterial wall found in the literature with a correlation coefficient of 0.985. Increasing age had a strong positive correlation with the shoulder stress level, (r = 0.95). The 30% stenosis had the highest shoulder stress due to the combination of a fully formed lipid pool and a thin cap. CONCLUSIONS: Studying the solid mechanics of the arterial wall and the atheroma provide important insights into the mechanisms involved in plaque rupture.

NASA Program Biomedical Research and Countermeasur↗

NASA NTRS2011

Assessment of CFD-based Response Surface Model for Ares I Supersonic Ascent Aerodynamics

The Ascent Force and Moment Aerodynamic (AFMA) Databases (DBs) for the Ares I Crew Launch Vehicle (CLV) were typically based on wind tunnel (WT) data, with increments provided by computational fluid dynamics (CFD) simulations for aspects of the vehicle that could not be tested in the WT tests. During the Design Analysis Cycle 3 analysis for the outer mold line (OML) geometry designated A106, a major tunnel mishap delayed the WT test for supersonic Mach numbers (M) greater than 1.6 in the Unitary Plan Wind Tunnel at NASA Langley Research Center, and the test delay pushed the final delivery of the A106 AFMA DB back by several months. The aero team developed an interim database based entirely on the already completed CFD simulations to mitigate the impact of the delay. This CFD-based database used a response surface methodology based on radial basis functions to predict the aerodynamic coefficients for M > 1.6 based on only the CFD data from both WT and flight Reynolds number conditions. The aero team used extensive knowledge of the previous AFMA DB for the A103 OML to guide the development of the CFD-based A106 AFMA DB. This report details the development of the CFD-based A106 Supersonic AFMA DB, constructs a prediction of the database uncertainty using data available at the time of development, and assesses the overall quality of the CFD-based DB both qualitatively and quantitatively. This assessment confirms that a reasonable aerodynamic database can be constructed for launch vehicles at supersonic conditions using only CFD data if sufficient knowledge of the physics and expected behavior is available. This report also demonstrates the applicability of non-parametric response surface modeling using radial basis functions for development of aerodynamic databases that exhibit both linear and non-linear behavior throughout a large data space.

Hanke, Jeremy L.↗

NASA NTRS2012

Digital Architecture for a Trace Gas Sensor Platform

A digital architecture has been implemented for a trace gas sensor platform, as a companion to standard analog control electronics, which accommodates optical absorption whose fractional absorbance equivalent would result in excess error if assumed to be linear. In cases where the absorption (1-transmission) is not equivalent to the fractional absorbance within a few percent error, it is necessary to accommodate the actual measured absorption while reporting the measured concentration of a target analyte with reasonable accuracy. This requires incorporation of programmable intelligence into the sensor platform so that flexible interpretation of the acquired data may be accomplished. Several different digital component architectures were tested and implemented. Commercial off-the-shelf digital electronics including data acquisition cards (DAQs), complex programmable logic devices (CPLDs), field-programmable gate arrays (FPGAs), and microcontrollers have been used to achieve the desired outcome. The most completely integrated architecture achieved during the project used the CPLD along with a microcontroller. The CPLD provides the initial digital demodulation of the raw sensor signal, and then communicates over a parallel communications interface with a microcontroller. The microcontroller analyzes the digital signal from the CPLD, and applies a non-linear correction obtained through extensive data analysis at the various relevant EVA operating pressures. The microcontroller then presents the quantitatively accurate carbon dioxide partial pressure regardless of optical density. This technique could extend the linear dynamic range of typical absorption spectrometers, particularly those whose low end noise equivalent absorbance is below one-part-in-100,000. In the EVA application, it allows introduction of a path-length-enhancing architecture whose optical interference effects are well understood and quantified without sacrificing the dynamic range that allows quantitative detection at the higher carbon dioxide partial pressures. The digital components are compact and allow reasonably complete integration with separately developed analog control electronics without sacrificing size, mass, or power draw.

Gonzales, Paula↗

NASA NTRS2012

Comparison of Aircraft Models and Integration Schemes for Interval Management in the TRACON

Reusable models of common elements for communication, computation, decision and control in air traffic management are necessary in order to enable simulation, analysis and assurance of emergent properties, such as safety and stability, for a given operational concept. Uncertainties due to faults, such as dropped messages, along with non-linearities and sensor noise are an integral part of these models, and impact emergent system behavior. Flight control algorithms designed using a linearized version of the flight mechanics will exhibit error due to model uncertainty, and may not be stable outside a neighborhood of the given point of linearization. Moreover, the communication mechanism by which the sensed state of an aircraft is fed back to a flight control system (such as an ADS-B message) impacts the overall system behavior; both due to sensor noise as well as dropped messages (vacant samples). Additionally simulation of the flight controller system can exhibit further numerical instability, due to selection of the integration scheme and approximations made in the flight dynamics. We examine the theoretical and numerical stability of a speed controller under the Euler and Runge-Kutta schemes of integration, for the Maintain phase for a Mid-Term (2035-2045) Interval Management (IM) Operational Concept for descent and landing operations. We model uncertainties in communication due to missed ADS-B messages by vacant samples in the integration schemes, and compare the emergent behavior of the system, in terms of stability, via the boundedness of the final system state. Any bound on the errors incurred by these uncertainties will play an essential part in a composable assurance argument required for real-time, flight-deck guidance and control systems,. Thus, we believe that the creation of reusable models, which possess property guarantees, such as safety and stability, is an innovative and essential requirement to assessing the emergent properties of novel airspace concepts of operation.

Neogi, Natasha↗

NASA NTRS2013

Trajectory Reconstruction and Uncertainty Analysis Using Mars Science Laboratory Pre-Flight Scale Model Aeroballistic Testing

As part of the Mars Science Laboratory (MSL) trajectory reconstruction effort at NASA Langley Research Center, free-flight aeroballistic experiments of instrumented MSL scale models was conducted at Aberdeen Proving Ground in Maryland. The models carried an inertial measurement unit (IMU) and a flush air data system (FADS) similar to the MSL Entry Atmospheric Data System (MEADS) that provided data types similar to those from the MSL entry. Multiple sources of redundant data were available, including tracking radar and on-board magnetometers. These experimental data enabled the testing and validation of the various tools and methodologies that will be used for MSL trajectory reconstruction. The aerodynamic parameters Mach number, angle of attack, and sideslip angle were estimated using minimum variance with a priori to combine the pressure data and pre-flight computational fluid dynamics (CFD) data. Both linear and non-linear pressure model terms were also estimated for each pressure transducer as a measure of the errors introduced by CFD and transducer calibration. Parameter uncertainties were estimated using a "consider parameters" approach.

Lugo, Rafael A.↗

NASA NTRS2013

3D Radiative Transfer Effects in Multi-Angle/Multi-Spectral Radio-Polarimetric Signals from a Mixture of Clouds and Aerosols Viewed by a Non-Imaging Sensor

When observing a spatially complex mix of aerosols and clouds in a single relatively large field-of-view, nature entangles their signals non-linearly through polarized radiation transport processes that unfold in the 3D position and direction spaces. In contrast, any practical forward model in a retrieval algorithm will use only 1D vector radiative transfer (vRT) in a linear mixing technique. We assess the difference between the observed and predicted signals using synthetic data from a high-fidelity 3D vRT model with clouds generated using a Large Eddy Simulation model and an aerosol climatology. We find that this difference is signal--not noise--for the Aerosol Polarimetry Sensor (APS), an instrument developed by NASA. Moreover, the worst case scenario is also the most interesting case, namely, when the aerosol burden is large, hence hase the most impact on the cloud microphysics and dynamics. Based on our findings, we formulate a mitigation strategy for these unresolved cloud adjacency effects assuming that some spatial information is available about the structure of the clouds at higher resolution from "context" cameras, as was planned for NASA's ill-fated Glory mission that was to carry the APS but failed to reach orbit. Application to POLDER (POLarization and Directionality of Earth Reflectances) data from the period when PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar) was in the A-train is briefly discussed.

APS↗

NASA NTRS2015

Regional Rates of Young US Forest Growth Estimated From Annual Landsat Disturbance History and IKONOS Stereo Imagery

Forests of the Contiguous United States (CONUS) have been found to be a large contributor to the global atmospheric carbon sink. The magnitude and nature of this sink is still uncertain and recent studies have sought to define the dynamics that control its strength and longevity. The Landsat series of satellites has been a vital resource to understand the long-term changes in land cover that can impact ecosystem function and terrestrial carbonstock. We combine annual Landsat forest disturbance history from 1985 to 2011 with single date IKONOS stereoimagery to estimate the change in young forest canopy height and above ground live dry biomass accumulation for selected sites in the CONUS. Our approach follows an approximately linear growth rate following clearing over short intervals and does not estimate the distinct non-linear growth rate over longer intervals.We produced canopy height models by differencing digital surface models estimated from IKONOS stereo pairs with national elevation data (NED). Correlations between height and biomass were established independently using airborne LiDAR, and then applied to the IKONOS-estimated canopy height models. Graphing current biomass against time since disturbance provided biomass accumulation rates. For 20 study sites distributed across five regions of the CONUS, 19 showed statistically significant recovery trends (p is less than 0.001) with canopy growth from 0.26 m yr-1to 0.73 m yr-1. Above ground live dry biomass (AGB) density accumulation ranged from 1.31 t/ha yr-1 to 12.47 t/ha yr-1. Mean forest AGB accumulationwas 6.31 t/ha yr-1 among all sites with significant growth trends. We evaluated the accuracy of our estimates by comparing to field estimated site index curves of growth, airborne LiDAR data, and independent model predictions of C accumulation. Growth estimates found with this approach are consistent with site index curves and total biomass estimates fall within the range of field estimates. This is aviable approach to estimate forest biomass accumulation in regions with clear-cut harvest disturbances.

Neigh, Christopher S. R.↗

NASA NTRS2022

Present-Day and Future PM 2.5 and O 3 -Related Global and Regional Premature Mortality in the EVAv6.0 Health Impact Assessment Model

We used the EVAv6.0 system to estimate the present (2015) and future (2015–2050) global PM 2.5 and O 3 -related premature mortalities, using simulated surface concentrations from the GISS-E2.1-G Earth system model. The PM 2.5 -related global premature mortality is estimated to be 4.3 and 4.4 million by the non-linear and linear models, respectively. Ischemic heart diseases are found to be the leading cause of PM 2.5 -related premature deaths, contributing by 35% globally. Both long-term and short-term O 3 -related premature deaths are estimated to be around 1 million, globally. Overall, PM 2.5 and O 3 -related premature mortality leads to 5.3–5.4 million premature deaths, globally. The global burden of premature deaths is mainly driven by the Asian region, which in 2015 contributes by 75% of the total global premature deaths. An increase from 6.2% to 8% in the PM 2.5 relative risk as recommended by the WHO leads to an increase of PM 2.5 -related premature mortality by 28%, to 5.7 million. Finally, bias correcting the simulated PM 2.5 concentrations in 2015 leads to an increase of up to 73% in the global PM 2.5 -related premature mortality, leading to a total number of global premature deaths of up to 7.7 million, implying the necessity of bias correction to get more robust health burden estimates. PM 2.5 and O 3 -related premature mortality in 2050 decreases by up to 57% and 18%, respectively, due to emission reductions alone. However, the projected increase and aging of the population leads to increases of premature mortality by up to a factor of 2, showing that the population exposed to air pollution is more important than the level of air pollutants, highlighting that the population dynamics should be considered when setting up health assessment systems.

Premature mortality↗

NASA NTRSDate not supplied

Stagnation Point Convective Heating Correlations for Entry Into H2/He Atmospheres

New convective heating correlations have been developed over a broad range of entry conditions into H2/He atmospheres. The correlations are based on calculations of stagnation point convective heat fluxes over a large range of velocities, densities, and effective nose radii. The heat fluxes and corresponding flowfields are computed using NASA’s hypersonic computational fluid dynamics code, DPLR. Cases in this study have been computed using both a temperaturelimited, radiative equilibrium wall and a cold wall temperature boundary condition. The results for the scaled heat flux, as a function of the difference in freestream and wall enthalpies, show a clear deviation from the linear trend predicted by the method of Sutton and Graves. Therefore, updated non-linear correlations are proposed for the prediction of stagnation point convective heat fluxes. Generally, these updated correlations predict lower heat loads and therefore smaller required TPS thicknesses when compared to previous correlations for a range of test cases in the Neptune, Saturn, and Uranus atmospheres.

Alex T. Carroll↗

NASA NTRSDate not supplied

Stagnation Point Convective Heating Correlations for Entry Into H2/He Atmospheres

New convective heating correlations have been developed over a broad range of entry conditions into H2/He atmospheres. The correlations are based on calculations of stagnation point convective heat fluxes over a large range of velocities, densities, and effective nose radii. The heat fluxes and corresponding flow fields are computed using NASA’s hypersonic computational fluid dynamics code, DPLR. Cases in this study have been computed using both a temperature limited, radiative equilibrium wall and a cold wall temperature boundary condition. The results for the scaled heat flux, as a function of the difference in freestream and wall enthalpies, show a clear deviation from the linear trend predicted by the method of Sutton and Graves. Therefore, updated non-linear correlations are proposed for the prediction of stagnation point convective heat fluxes. Generally, these updated correlations predict lower heat loads and therefore smaller required TPS thicknesses when compared to previous correlations for a range of test cases in the Neptune, Saturn, and Uranus atmospheres.

Alex T. Carroll↗

DOE OSTI2018

The Effects of Market Concentration on Residential Solar PV Prices: Competition, Installer Scale, and Soft Costs

Competition among residential solar photovoltaic (PV) installers may reduce PV price markups and yield lower prices. At the same time, competition may reduce installer experience and opportunities for cost reductions through learning-by-doing and economies of scale. These dynamics suggest that PV non-hardware or 'soft' costs and installed prices depend on the distribution of market shares among installers, also known as market structure. This study leverages a rich data set of 226,769 residential PV systems to examine the relationship between market structure, soft costs, and PV prices. The results show that PV prices are lower, on average, in more concentrated markets, i.e., markets with fewer installers or where a subset of installers holds a large market share. The study provides evidence that this relationship is non-linear, such that prices are minimized in markets with an optimal balance of the benefits of market concentration and the benefits of competition.

14 SOLAR ENERGY↗

DOE OSTI2019

Complete the Design of Benchtop Fiber-Optic Based Photo-Thermal Radiometry System

Thermal transport properties of nuclear fuels, i.e., thermal conductivity and diffusivity, are critical parameters in the thermal design of a nuclear reactor. It is well known that thermal transport properties change temporally and spatially due to dynamic microstructural evolution as a result of neutron irradiation. A real-time instrument to monitor the local thermal property change, which can survive the harsh in-pile environment, will give researchers unprecedented access to fuel behavior during irradiation. Laser-based photothermal techniques are an ideal candidate to accomplish this specific task, because of the remote, non-contact, and non-destructive features. In this report, we present the progress in the development of a bench-top photothermal radiometry (PTR) system to measure the thermal diffusivity of solid samples, and the preliminary design of a fiber-based PTR instrument to measure the real-time, in-pile thermal diffusivity of nuclear fuels. The PTR technique measures the local temperature response from a modulated heat source by detecting the sample radiation heat flux. Initially a well-developed thermal wave analytical model was fit to the data to obtain the thermal diffusivity. However, systematic errors from the diffraction effect and non-linearity issues resulted in thermal diffusivity measurements higher than literature values. These effects have been studied and an improved analytical model has been developed, and the physics behind the phenomena are still under investigation. A procedure for optimizing measurements from the bench-top system has been provided. The experiences from the bench top measurement were then used to guide the fiber-based system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Midfidelity Model Verification for a Point-Absorbing Wave Energy Converter with Linear Power Takeoff: Preprint

In the preliminary design stage of a wave energy converter (WEC), fast and reliable simulation tools are required. High-fidelity numerical models are usually employed to study the wave-structure interaction, yet the computational cost is demanding. Instead, mid-fidelity models provide simulations in the order of real time. In this study, Uppsala University WEC operates in a relatively mild sea state and is modeled using WEC-Sim. The model is verified based on OpenFOAM simulations. To analyze the ability of mid-fidelity model to capture the WEC dynamics, the system is investigated separately in 1, 2 and 3 DoF. The contribution of viscous phenomena is examined. Moreover, linear and weakly non-linear solutions provided by WEC-Sim are studied. The results obtained indicate that the viscous effects in heave and surge motion can be neglected but not for the pitch. The weakly non-linear WEC-Sim solution successfully agrees with the CFD, whereas the linear solution could drive to misleading results.

CFD↗

DOE OSTI2021

Midfidelity Model Verification for a Point-Absorbing Wave Energy Converter with Linear Power Take-Off

CFD↗

DOE OSTI2022

Stability analysis of a one-dimensional multiphysics model of a molten salt fast reactor

Reactor designs with a fast neutron spectrum have recently moved into the research focus for possible Generation IV reactor concepts - a representative of which is the molten salt fast reactor (MSFR). In order to assess the safety and reliability of this reactor concept, it is necessary to obtain knowledge of the system behaviour. An important part of this research concerns the dynamical stability of the MSFR. Research regarding the MSFR stability so far has either been using linearized models for stability criteria or only been examining certain transients at fixed parameter values. This work delivers a comprehensive stability analysis for a non-linearized MSFR model, and a wide range of values for all relevant parameters. For this purpose, a one-dimensional MSFR model was set up, taking neutron kinetics and thermal hydraulics in the form of a system of differential equations into account. The stability of this model was investigated by means of the numerical tool MATCONT, which was used to monitor the evolution of a so-called fixed-point solution, here referring to the steady state at operating conditions of the system, while system parameters got varied. The MATCONT analyses showed no loss of stability for any of the considered parameter variations and no solutions that might co-exist in parallel with the steady-state fixed-point solution. Therefore these results indicate a stable fixed point to which all solution transients converge, making the considered MSFR model insusceptible to deviations from the equilibrium. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗