Frequency monitoring network (FNET) is an Internet based GPS synchronized wide-area frequency monitoring network deployed at distribution level. At first part of this thesis, FNET structure and characteristics are introduced. After analysis and smoothing FDR signals, the algorithm of event trigger is present with Visual C++ DLL programming. Estimation of disturbance location method is discussed based on the time delay of arriving (TDOA) in the second part of this work. In this section, author shows the multiply method to calculate event time, which is important when deal with pre-disturbance frequency in TDOA part. Two event kinds are classified by the change of frequency and the linear relationship between change of frequency and imbalance of generation and load power is presented. Prove that Time Delay of Arrival (TDOA) is a good algorithm for estimation event location proved by real cases. At last, the interface of DLL module and the key word to import and export DLL variables and function is described. At last, PSS compensation optimization with a set of nonlinear differential algebraic equations (DAE) is introduced in detail. With combining the bifurcation theory of nonlinear system and the optimization theory, the optimal control of small-signal stability of power electric systems are solved. From the perspective of stability margin, global coordination of controller parameters is studied to ensure the stable operation of power grids. The main contents of this thesis include: ï¼ 1ï¼ Models of power systems and test power electric systems. Tht5e dynamic and static models of the elements of power systems, such as generatorbbs, AVRs, PSSs, loads and FACTS controllers are presented. Method of power system linearization modeling is introduced. Three test power systems, WSCC 9-bus system, 2-area system, New England 39-bus system, are used in thesis. ï¼ 2ï¼ Multi-objective optimizations based on bifurcation theory. The optimization models, damping control-Hopf bifurcation control, voltage control-damping control, are presented. Pareto combined with evolutionary strategy (ES) are used to solve multi-objective optimizations. Based on traditional PSS parameters optimizations, it can be formulated as a multi-objective problem, in which, two objectives should be taken into account. The minimum damping torque should be identified. / Master of Science

<p dir="ltr">The ability to predict and manipulate how a particulate material will flow in a process is challenging for industry and researchers alike. This dissertation presents the results of a model-directed, experimental approach using a concentric cylinder rheometer titled along an axis to enable continuous chute flow of granular media. Experiments were performed using draining flows for constant and oscillatory applied shear rates. Multiple flow and stress sensors were used to investigate the interaction of mass holdup, shear rate, specific torque, particle velocity, discharge mass flow rate, and wall pressure. Depending on the flow configuration, linear ranges were observed wherein the specific torque remained steady during draining. This finding enabled systematic testing of flow behavior as a function of dimensionless shear rates. Results suggest changes in the specific torque, wall slip, and outflow variance occur with the transition from the quasi-static to dense-inertial flow regimes. A pump-curve analogy was also identified for the relationship between the outlet mass flow rate and the specific power relationship for the constant shear rate experiments. Oscillatory shear rate experiments show a significant influence of the phase shift between the applied shear rate and the specific torque. Adding an asperity to the rotor revealed rate-dependent patterns in bulk flow and force chain dynamics. Overall, the study offers valuable insights into the effects of shear rate and boundary conditions on dense granular flows. The effects of particle characteristics (e.g., size and shape distributions, friction, cohesivity) and material properties (e.g., density, modulus) remain topics for future work. </p>

Les vagues de bateau ajoutent une pression supplémentaire sur les berges de rivières et doivent être considérées dans les modèles de prédiction des taux de recul des berges. L’objectif de cette étude est d’examiner le rôle des vagues de bateau sur l’écoulement et le transport en suspension le long des berges en milieu fluvial. Pour atteindre cet objectif, nous utilisons un transect perpendiculaire à la berge de quatre courantomètres électromagnétiques (ECMs) mesurant deux dimensions de l’écoulement et deux turbidimètres (OBSs) placés dos à dos, orientés vers la berge et le large pour mesurer les conditions moyennes et turbulentes de l’écoulement longitudinal et vertical ainsi que les flux de sédiments en suspension provoqués par les vagues. Une chaloupe à moteur de 16 pieds, équipée d’un moteur 40 hp, a été utilisée afin de générer des vagues. Nous avons mesuré l’effet de trois distances à partir de la berge (5, 10, 15 m) et trois vitesses de bateau (5, 15 et 25 km/h) et cinq répliques de chaque combinaison de distance et de vitesse ont été réalisées, totalisant 45 passages. Nous avons caractérisé la variabilité des conditions d’écoulement, de vagues et de transport de sédiments et nous avons réalisé des analyses spectrales afin de séparer les portions oscillatoire et turbulente de l’écoulement généré par les vagues de bateau. L’effet de la distance et de la vitesse du bateau sur le transport de sédiments est non-linéaire et la réponse sédimentaire induite par les passages de bateau montre une variabilité importante entre les répliques et les deux sondes OBS, ce qui suggère un changement morphologique induit par les vagues de bateau. Les corrélations entre les variables d’écoulement et de transport montrent l’importance des relations entre le cisaillement et la puissance de la portion turbulente de l’écoulement avec le transport de sédiments. Cette étude a permis de quantifier les relations entre la dynamique des vagues et les flux de concentrations de sédiments en suspension, ce qui représente une contribution importante au développement de mesures de mitigation dans les environnements fluviaux où les berges sont fragilisées par le trafic plaisancier. / Boat induced waves generate additional stress on banks and should be included in bank erosion predictive models. The objective of this study is to investigate the role of boat generated waves on near bank flow in a fluvial environment. We used a longitudinal array of four bi-directional eletromagnetic current meters (ECMs) perpendicular to the bank and two optical backcattering sensors (OBSs) set back to back, one facing the river bank (inner looking) and the other towards the river (outer looking) to measure mean and turbulent properties of the horizontal and vertical flow velocities and suspended sediment fluxes generated by the waves. The waves were generated using a 16 ft long boat equipped with a 40hp motor. The effect of three different distances between the bank and sailing line (5, 10, 15 m) and of three different speeds of the vessel (5, 15 and 25 km/h) was evaluated. Five replicates of each combination of distance and speed were realized for a total of 45 passages. Variability of flows conditions, wave characteristics and turbidity properties were characterized and spectral analyses were performed to separate oscillatory and turbulent flow induced by boat generated waves. An important variability in suspended sediment response to the passages of boat wave is observed between replications and the effect of speed and distance of sailing line is nonlinear. It was often unbalanced between the inner and outer looking OBS suggesting that the bank changed its morphology. Cross-correlation between hydraulic and turbidity variables show the significance of the relation between shear stress and suspended sediment transport variables. This study provides quantitative relations between wave dynamics and plumes of suspended sediment that could help to develop mitigation measures in fluvial environments where vessel traffic represents a major issue for bank erosion and retreat.

The goals of this thesis are the characterization of polymer melts using mainly non-linear shear and extensional rheological techniques. The fabrication of scaffolds with excellent physical and mechanical properties using solution electrospinning technology for tissue engineering applications and the development of melt electrospinning equipment to facilitate the fabrication of solvent free scaffolds. To achieve the first goal, we focused on the characterization of entangled polymer melts in the linear and nonlinear viscoelastic regimes. The influence of molecular weight, Mw, molecular weight distribution (MWD), long-chain branching (LCB) and addition of particles to the polymer matrix on polymer melt properties were investigated using shear and extensional rheological techniques. The resulting structure-property relationships were established using newly introduced mechanical parameters under large amplitude oscillatory shear (LAOS) flow.

Les vagues de bateau ajoutent une pression supplémentaire sur les berges de rivières et doivent être considérées dans les modèles de prédiction des taux de recul des berges. L’objectif de cette étude est d’examiner le rôle des vagues de bateau sur l’écoulement et le transport en suspension le long des berges en milieu fluvial. Pour atteindre cet objectif, nous utilisons un transect perpendiculaire à la berge de quatre courantomètres électromagnétiques (ECMs) mesurant deux dimensions de l’écoulement et deux turbidimètres (OBSs) placés dos à dos, orientés vers la berge et le large pour mesurer les conditions moyennes et turbulentes de l’écoulement longitudinal et vertical ainsi que les flux de sédiments en suspension provoqués par les vagues. Une chaloupe à moteur de 16 pieds, équipée d’un moteur 40 hp, a été utilisée afin de générer des vagues. Nous avons mesuré l’effet de trois distances à partir de la berge (5, 10, 15 m) et trois vitesses de bateau (5, 15 et 25 km/h) et cinq répliques de chaque combinaison de distance et de vitesse ont été réalisées, totalisant 45 passages. Nous avons caractérisé la variabilité des conditions d’écoulement, de vagues et de transport de sédiments et nous avons réalisé des analyses spectrales afin de séparer les portions oscillatoire et turbulente de l’écoulement généré par les vagues de bateau. L’effet de la distance et de la vitesse du bateau sur le transport de sédiments est non-linéaire et la réponse sédimentaire induite par les passages de bateau montre une variabilité importante entre les répliques et les deux sondes OBS, ce qui suggère un changement morphologique induit par les vagues de bateau. Les corrélations entre les variables d’écoulement et de transport montrent l’importance des relations entre le cisaillement et la puissance de la portion turbulente de l’écoulement avec le transport de sédiments. Cette étude a permis de quantifier les relations entre la dynamique des vagues et les flux de concentrations de sédiments en suspension, ce qui représente une contribution importante au développement de mesures de mitigation dans les environnements fluviaux où les berges sont fragilisées par le trafic plaisancier. / Boat induced waves generate additional stress on banks and should be included in bank erosion predictive models. The objective of this study is to investigate the role of boat generated waves on near bank flow in a fluvial environment. We used a longitudinal array of four bi-directional eletromagnetic current meters (ECMs) perpendicular to the bank and two optical backcattering sensors (OBSs) set back to back, one facing the river bank (inner looking) and the other towards the river (outer looking) to measure mean and turbulent properties of the horizontal and vertical flow velocities and suspended sediment fluxes generated by the waves. The waves were generated using a 16 ft long boat equipped with a 40hp motor. The effect of three different distances between the bank and sailing line (5, 10, 15 m) and of three different speeds of the vessel (5, 15 and 25 km/h) was evaluated. Five replicates of each combination of distance and speed were realized for a total of 45 passages. Variability of flows conditions, wave characteristics and turbidity properties were characterized and spectral analyses were performed to separate oscillatory and turbulent flow induced by boat generated waves. An important variability in suspended sediment response to the passages of boat wave is observed between replications and the effect of speed and distance of sailing line is nonlinear. It was often unbalanced between the inner and outer looking OBS suggesting that the bank changed its morphology. Cross-correlation between hydraulic and turbidity variables show the significance of the relation between shear stress and suspended sediment transport variables. This study provides quantitative relations between wave dynamics and plumes of suspended sediment that could help to develop mitigation measures in fluvial environments where vessel traffic represents a major issue for bank erosion and retreat.

L’utilisation de réseaux de neurones artificiels pour divers types de traitements d’information bio-inspirés est une technique de plus en plus répandue dans le domaine de l’intelligence artificielle. Leur fonctionnement diffère avantageusement de celui des ordinateurs conventionnels en permettant une plus grande parallélisation des calculs, ce qui explique pourquoi autant d’efforts sont déployés afin de réaliser une plate-forme matérielle dédiée à leur simulation. Pour ce projet, une architecture matérielle flexible simulant efficacement un réseau de neurones à décharges est présentée. Celle-ci se distingue des architectures existantes notamment parce qu’elle utilise une approche de simulation à événements discrets et parce qu’elle permet une détection efficace des événements simultanés. Ces caractéristiques en font une plate-forme de choix pour la simulation de réseaux de neurones à décharges de plus de 100 000 neurones où un niveau important de synchronie des décharges neuronales est atteint. Afin d’en démontrer les performances, une application en traitement d’images utilisant cette architecture a été réalisée sur FPGA. Cette application a permis de démontrer que la structure proposée pouvait simuler jusqu’à 2[indice supérieur 17] neurones et traiter des dizaines de millions d’événements par secondes lorsque cadencé à 100 MHz.

Objetivos: avaliar as alterações morfofuncionais da retina de pacientes com esclerose múltipla (EM) e espectro da neuromielite óptica (ENMO), com ou sem histórico de neurite óptica (NO), por meio de eletrorretinografia de campo total e multifocal (ERGct / ERGmf) e a tomomografia de coerência óptica (TCO). Avaliar as correlações de tais achados entre si e com o potencial evocado visual (PEV), o campo visual (CV) e a sensibilidade ao contraste (SC). Métodos: Pacientes com EM (n = 30), ENMO (n = 30) e controles (n = 29) foram submetidos a avaliação oftalmológica completa incluindo CV, TCO, ERGct, ERGmf, PEV e medida da SC. Os olhos foram distribuídos em 5 grupos: EM com ou sem história de NO (EM + NO e EM - NO), ENMO com ou sem NO (ENMO + NO e ENMO - NO) e controles. Com a TCO foram medidas as espessuras das camadas de fibras nervosas da retina na região macular (CFNRm), camada de células ganglionares (CCG), camada plexiforme interna (CPI), camada nuclear interna (CNI), camada plexiforme externa (CPE), camada nuclear externa (CNE), complexo camada de fotorreceptores + epitélio pigmentário da retina(CFR+EPR) e da camada de fibras nervosas peripapilar (CFNRp). Dados das ondas a, b e potenciais oscilatórios (POs) do ERGct, medidas de amplitude e latência de N1 e P1 do ERGmf, ondas do PEV e medidas do desvio da normalidade do CV foram analisados. Os grupos foram comparados usando equações estimadas generalizadas. Correlações entre as medidas foram avaliadas. Resultados: Redução da sensibilidade do desvioda normalidade do CV, da SC e aumento da latência das ondas do PEV foram identificadas para ambos os grupos. Nos pacientes ENMO+NO apresentaram redução de amplitude ao PEV em relação aos controles, diferentemente aos demais grupos. As medidas de espessura da CCG e da CPI foram significativamente menores nos grupos de olhos dos pacientes em relação aos controles. A RNFLm foi menor em todos os grupos de olhos de pacientes, exceto o grupo ENMO-NO. Não foi observada diferença significativa entre os grupos de olhos estudados nas comparações referenetes às demais camadas da retina. Comparado aos controles, as amplitudes do POs foram maiores nos olhos de pacientes com ENMO, enquanto as latências de N1 e P1 do ERGmf foram menor nos olhos de pacientes com EM. Essas anormalidades foram fortemente correlacionadas com a espessura da camada retiniana intermediária e externa. Os achados do PEV, assim como do CV e SC, se correlacionaram fortemente com as camadas retinianas. Conclusões: as camadas retinianas internas se mostraram reduzidas à TCO tanto nos olhos de pacientes com EM quanto no naqueles com ENMO, mas os achados POs e ERGmf sugerem também envolvimento das outras camadas retinianas nessas afecções. O PEV apresentou alterações distintas para cada doença. O uso combinado da TCO, do ERG e PEV podem ajudar a entender como as duas condições diferem em relação aos danos retinianos / Objectives: To evaluate the morphofunctional alterations of the retina of patients with multiple sclerosis (MS) and spectrum of neuromyelitis optica spectrum disorder (NMOSD), with or without a history of optic neuritis (ON), using full field electroretinogram (ERG) and multifocal electroretinography (mf- ERG) and optical coherence tomography (OCT). To evaluate the correlations of such findings among themselves and with visual evoked potential (VEP), visual field (VF) and contrast sensitivity (CS). Methods: Patients with MS (n = 30), NMOSD (n = 30) and healthy controls (n = 29) were submitted to a complete ophthalmologic evaluation including VF, OCT, ERG, mf-ERG, VEP and SC measurement. The eyes were distributed in 5 groups: MS with or without history of ON (MS + ON and MS - ON), NMOSD with or without ON (NMOSD + ON and NMOSD - ON) and controls. With the OCT were measured the thickness of the retinal nerve fiber layers in the macular region (mRNFL), ganglion cell layer (GCL), inner plexiform layer (IPL), inner nuclear layer (INL), outer plexiform layer (OPL), outer nuclear layer (ONL), photoreceptor layer (PHOT) and peripapillary retinal nerve fiber layer (pRNFL). The data of the a- and b-waves and oscillatory potentials (OPs) of the ERG, amplitude and peak time of N1 and P1 of mf-ERG, waves of the VEP and VF deviation were analyzed. The groups were compared using generalized estimating equations. Correlations between measurements were evaluated. Results: Reduction of the oVF deviation, CS and increased VEP wave\'s peak times were identified for both groups. In the NMOSD + ON patients, an amplitude redeuction was found in the VEP in relation to the controls, unlike the other groups. The thickness of GCL and IPL was significantly lower in patients\' eyes than controls. mRNFL was lower in all patients, but NMOSD-ON. No significant difference was observed for the remaining layers. Compared to controls, the amplitudes of the POs were higher in the NMOSD group, whereas the mf-ERG\'s N1 and P1 peak time was lower in the MS patients. These abnormalities were strongly correlated with the thickness of the intermediate and outer retinal layers. The findings of the VEP, as well as the VF and SC, correlated strongly with the retinal layers. Conclusions: the inner retinal layers were reduced in both the MS and the NMOSD, but the findings OPs and mf-ERG suggest involvement of the other retinal layers. The VEP presented different alterations for each disease. The combination of OCT, ERG and VEP may help to understand how the two conditions differ in relation to retinal damage

Les mécanismes volontaires (V) et involontaires (I) de l’attention reposent sur les réseaux dorsal et ventral, convergeant dans le cortex préfrontal latéral (lPFC). La distractibilité accrue liée au vieillissement ou à une lésion frontale pourrait être due à une altération de l’équilibre entre ces mécanismes V et I, essentiel mais rarement étudié. Notre objectif est de tester, dans la modalité auditive, si (1) les oscillations alpha coordonnent l'activité du réseau dorsal, (2) les oscillations gamma celle du réseau ventral, (3) le couplage oscillatoire dans le lPFC maintient l’équilibre entre les deux réseaux. Ce travail vise également à étudier les corrélats oscillatoires de la distractibilité accrue liée au vieillissement ou à une atteinte frontale. Des données MEEG ont été enregistrées alors que des participants réalisaient le Competitive Attention Test, qui permet d’étudier simultanément les mécanismes V et I de l’attention. Nous avons montré que les oscillations alpha reflètent l’activation des mécanismes facilitateurs et suppresseurs de l’attention V, et la communication au sein du réseau dorsal ; alors que les oscillations gamma indexent l’activation du réseau ventral. De plus, le lPFC serait impliqué dans la communication au sein des deux réseaux, et le PFC médian dans l’équilibre attentionnel V/I. Nous avons également montré que la distractibilité accrue était liée à un déficit d’attention V au cours du vieillissement, et à une altération des processus V et I après lésion frontale. Ce travail de thèse offre donc une meilleure compréhension de la dynamique cérébrale oscillatoire sur laquelle repose l'équilibre attentionnel V/I, et donc la distractibilité / Top-down (TD) and bottom-up (BU) mechanisms of attention are supported by dorsal and ventral networks that mainly overlap in the lateral prefrontal cortex (lPFC). A balance between these mechanisms is essential, yet rarely investigated. Increased distractibility observed during ageing or after frontal damage could result from jeopardizing this balance. It has been proposed that distinct oscillatory frequencies support the activation of these two attention networks. Our main aim was to test, in the auditory modality, whether (1) alpha oscillations would coordinate activity within the dorsal TD network, (2) gamma activity would index the activation of the ventral BU network, (3) the lPFC would support the balance between these networks through oscillatory coupling. We also aimed to investigate the oscillatory correlates of the increased distractibility associated with ageing or frontal damage. MEEG data were recorded while participants performed the Competitive Attention Test, which enables simultaneous investigation of BU and TD attention mechanisms. We showed that alpha oscillations indexed facilitatory and suppressive mechanisms of TD attention, and communication within the dorsal network; while gamma oscillations indexed the ventral network activation. Moreover, the lPFC subtended communication in the two networks; with the TD/BU interaction occurring in the medial PFC. We also showed that ageing-related distractibility was of TD deficit origin. Finally, preliminary results suggest that lPFC damage can impact both TD and BU attention. This thesis provides novel insights into the brain oscillatory dynamics of the TD/BU attentional balance supporting distractibility

The structure and dynamics of the smectic-A liquid crystal 8CB (4 cyano-4 octylbiphenyl) when sheared and confined to mesoscale gaps (with crossed cylindrical geometry and mica confining surfaces) were studied using a Surface Forces Apparatus (SFA). Triangular shear patterns with frequencies of 0.01, 0.1, 1.0 and 10 Hz, and amplitudes of 62.5 nm, 625 nm and 6.25 m were applied to samples at gap sizes of 0.5 and 5.0 m. The study was performed at room temperature (20.5C) and at two higher temperatures (22C and 27C). In order to minimize the thermal fluctuations within the test chamber and hence to allow for the rapid re-initialization of test runs, the SFA was modified to allow for quick, precise and remote control of the confining surfaces. The procedure maximized the number of tests that could be undertaken with a single pair of surfaces so that a single gap geometry could be maintained for the duration of the test run. In order to run the SFA remotely, scripts written with a commercial software package, LabVIEW, were used to control of the SFA components, its FECO-monitoring camera and all its peripheral electronic equipment as well. Samples were agitated to disrupt any shear-induced liquid crystal domain alignment from previous testing following each shear test, and methodologies were developed to ascertain the extent of confinement quickly and remotely following agitation. Separate methods were developed for gap sizes at each extreme of the mesoscale regime, where the transition from bulklike structure and dynamics to nano-confinement occurs (between 1 and 10 microns for smectic-A 8CB). The results revealed that the greater amplitude-gap aspect ratio and surface-to-domain contact associated with smaller gaps facilitated reorientation of the domains in the shear direction. Evidence was also presented of domains at the higher end or outside of the mesoscale regime that, while straining and accreting, were unable to reorient and thereby led to an overall increase of viscoelastic response. The effective viscosity was found to obey a simple power law with respect to shear rate, , and the flow behaviour indices, n, slightly in excess of unity indicate shear thickening occurs with large enough shear amplitude, and that the viscosity reached a plateau near unity over shear rates of 0.005 to 500 s-1 within the mesoscale regime. Different K and n values were observed depending on the shear amplitude used. Unlike bulk smectic 8CB, whose domains do not align well in the shear direction with large shear-strain amplitude, at mesoscale levels of confinement large amplitude shearing (up to 12.5 shear strain amplitude) was found to be very effective at aligning domains. In general domain reorientation is found to be much more rapid within the mesoscale regime than has been reported in bulk. Aggressive shearing was found to result in a complete drop in viscoelastic response within seconds, while gentler shearing is found to produce a very gradual increase that persists for more than six hours, with individual shear periods exhibiting frequent and significant deviations from the expected smooth shear path that may be a product of discrete domain reorientations. From these findings, certain traits of the smectic 8CB domain structures under mesoscale confinement were deduced, including how they respond to shear depending on the level of confinement, and how their reorientation due to shear varies not only with shear rate but also independently with shear amplitude. An equation describing the viscosity change as a function of both shear rate and shear amplitude is proposed. The shear amplitude dependence introduces the notion of shearing beyond the proposed smectic 8CB “viscoelastic limit”, which was shown to exhibit behaviour in accordance with Large Amplitude Oscillatory Shear (LAOS) techniques developed for Fourier Transform rheology. The findings provided an understanding of the behavioural changes that occur as one reduces the level of confinement of smectic materials from bulk to nanoconfinement.

[1] Four specimens of zoned tourmaline from granitic pegmatites are characterised in detail, each having unusual compositional and/or morphologic features: (1) a crystal from Black Rapids Glacier, Alaska, showing a central pink zone of elbaite mantled by a thin rim of green liddicoatite; (2) a large (~25 cm) slab of Madagascar liddicoatite cut along (001) showing complex patterns of oscillatory zoning; and (3) a wheatsheaf and (4) a mushroom elbaite from Mogok, Myanmar, both showing extensive bifurcation of fibrous crystals originating from a central core crystal, and showing pronounced discontinuous colour zoning. Crystal chemistry and crystal structure of these samples are characterised by SREF, EMPA, and 11B and 27Al MAS NMR and Mössbauer spectroscopies. For each sample, compositional change, as a function of crystal growth, is characterised by EMPA traverses, and the total chemical variation is reduced to a series of linear substitution mechanisms. Of particular interest are substitutions accommodating the variation in [4]B: (1) TB + YAl ↔ TSi + Y(Fe, Mn)2+, where transition metals are present, and (2) TB2 + YAl ↔ TSi2 + YLi, where transition metals are absent. Integration of all data sets delineates constraints on melt evolution and crystal growth mechanisms. [2] Uncertainty has surrounded the occurrence of [4]Al and [4]B at the T-site in tourmaline, because B is difficult to quantify by EMPA and Al is typically assigned to the octahedral Y- and Z-sites. Although both [4]Al and [4]B have been shown to occur in natural tourmalines, it is not currently known how common these substituents are. Using 11B and 27Al MAS NMR spectroscopy, the presence of [4]B and [4]Al is determined in fifty inclusion-free tourmalines of low transition-metal content with compositions corresponding to five different species. Chemical shifts of [4]B and [3]B in 11B spectra, and [4]Al and [6]Al in 27Al spectra, are well-resolved, allowing detection of very small (< ~0.1 apfu) amounts of T-site constituents. Results show that contents of 0.0 < [4]B, [4]Al < 0.5 apfu are common in tourmalines containing low amounts of paramagnetic species, and that all combinations of Si, Al and B occur in natural tourmalines.