An experimental/numerical study of the turbulent boundary layer development along a surface with a sudden change in roughness

A theory for the description of turbulent boundary layer flows over surfaces with a sudden change in roughness is considered. The theory resorts to the concept of displacement in origin to specify a wall function boundary condition for a kappa-epsilon model. An approximate algebraic expression for the displacement in origin is obtained from the experimental data by using the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963). This expression is subsequently included in the near wall logarithmic velocity profile, which is then adopted as a boundary condition for a kappa-epsilon modelling of the external flow. The results are compared with the lower atmospheric observations made by Bradley(Q. J. Roy. Meteo. Soc., vol. 94, pp. 361-379, 1968) as well as with velocity profiles extracted from a set of wind tunnel experiments carried out by Avelino et al.( 7th ENCIT, 1998). The measurements are found to be in good agreement with the theoretical computations. The skin-friction coefficient was calculated according to the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963) and to a balance of the integral momentum equation. In particular, the growth of the internal boundary layer thickness obtained from the numerical simulation is compared with predictions of the experimental data calculated by two methods, the "knee" point method and the "merge" point method.

Application of a non isotropic turbulence model to stable atmospheric flows and dispersion over 3D topography

A non isotropic turbulence model is extended and applied to three dimensional stably stratified flows and dispersion calculations. The model is derived from the algebraic stress model (including wall proximity effects), but it retains the simplicity of the "eddy viscosity" concept of first order models. The "modified k-epsilon" is implemented in a three dimensional numerical code. Once the flow is resolved, the predicted velocity and turbulence fields are interpolated into a second grid and used to solve the concentration equation. To evaluate the model, various steady state numerical solutions are compared with small scale dispersion experiments which were conducted at the wind tunnel of Mitsubishi Heavy Industries, in Japan. Stably stratified flows and plume dispersion over three distinct idealized complex topographies (flat and hilly terrain) are studied. Vertical profiles of velocity and pollutant concentration are shown and discussed. Also, comparisons are made against the results obtained with the standard k-epsilon model.

Veden desinfiointi UV-LED-valolla

Veden riittämätön puhdistus aiheuttaa riskin veden käyttäjille. Miljoonia kuolemia vuosittain aiheuttavien vesiteitse leviävien sairauksien ehkäisemiseksi vaaditaan tehokkaita juomaveden desinfiointimenetelmiä. Kuivuuden ja väestönkasvun myötä veden tarve on lisääntynyt ja vedenkulutus tulee yhä kasvamaan. Tästä syystä mahdollisuus kierrättää vettä hyödyntäen sitä esimerkiksi kasteluun on saanut yhä enemmän huomiota. Kierrätettävä vesi on kuitenkin käsiteltävä huolellisesti sen sisältämän mikrobiologisen kontaminaatioriskin vuoksi. Ultraviolettisäteily luokitellaan fysikaaliseksi desinfiointimenetelmäksi. Sen tehokkuus perustuu mikro-organismien absorboimaan UV-säteilyyn, jonka aiheuttamien DNA:ssa tai RNA:ssa tapahtuvien muutoksien seurauksena mikro-organismi inaktivoituu ja estyy lisääntymästä. UV-desinfioinnissa on tyypillisesti käytetty elohopeahöyrylamppuja. Vaihtoehtoinen UV-säteilyn lähde ovat LEDit eli valoa emittoivat diodit. Matalapaine-elohopeahöyrylamppujen emittoima säteily on aallonpituudella 254 nm ja keskipaine-elohopeahöyrylamppujen emittoima säteily on laajakaistaista säteilyä. Energiatehokkuuden lisäksi LEDien etuna on, että niillä voidaan tuottaa kapeakaista säteilyä aallonpituudella, joka parhaiten absorboituu DNA:han. Tämän diplomityön tarkoituksena oli tutkia, onko UVC-alueen aallonpituuksien yhdistelmillä synergistisiä etuja LEDien desinfiointitehokkuuteen, kun desinfioidaan virtaavaa vettä useilla säteilyannoksilla ja indikaattorimikrobina käytetään kolibakteeria. Tavoitteena oli myös tutkia tällä hetkellä saatavissa olevien LEDien desinfiointitehokkuutta energiatehokkuuden näkökulmasta. Yksittäisistä aallonpituuksista desinfiointitehokkuudeltaan parhaimmaksi osoittautui 260 nm, aallonpituuksien yhdistelmistä tehokkain oli 265 nm:n ja 260 nm:n yhdistelmä. Muilla aallonpituuksien yhdistelmillä ei saavutettu odotettua parempaa desinfiointitehokkuutta. Optiselta teholtaan parhaimmat LEDit, 265 nm, 270 nm ja 275 nm olivat kokeiden perusteella myös energiatehokkuuden kannalta tarkasteltuina parhaimmat sekä yksittäin että yhdistelminä. UVC-aallonpituuksia emittoivien LEDien optisen tehokkuuden paraneminen on edellytys LEDien hyödyntämiselle desinfioinnissa.

Power semiconductor nonlinearities in active du/dt output filtering

This doctoral thesis introduces an improved control principle for active du/dt output filtering in variable-speed AC drives, together with performance comparisons with previous filtering methods. The effects of power semiconductor nonlinearities on the output filtering performance are investigated. The nonlinearities include the timing deviation and the voltage pulse waveform distortion in the variable-speed AC drive output bridge. Active du/dt output filtering (ADUDT) is a method to mitigate motor overvoltages in variable-speed AC drives with long motor cables. It is a quite recent addition to the du/dt reduction methods available. This thesis improves on the existing control method for the filter, and concentrates on the lowvoltage (below 1 kV AC) two-level voltage-source inverter implementation of the method. The ADUDT uses narrow voltage pulses having a duration in the order of a microsecond from an IGBT (insulated gate bipolar transistor) inverter to control the output voltage of a tuned LC filter circuit. The filter output voltage has thus increased slope transition times at the rising and falling edges, with an opportunity of no overshoot. The effect of the longer slope transition times is a reduction in the du/dt of the voltage fed to the motor cable. Lower du/dt values result in a reduction in the overvoltage effects on the motor terminals. Compared with traditional output filtering methods to accomplish this task, the active du/dt filtering provides lower inductance values and a smaller physical size of the filter itself. The filter circuit weight can also be reduced. However, the power semiconductor nonlinearities skew the filter control pulse pattern, resulting in control deviation. This deviation introduces unwanted overshoot and resonance in the filter. The controlmethod proposed in this thesis is able to directly compensate for the dead time-induced zero-current clamping (ZCC) effect in the pulse pattern. It gives more flexibility to the pattern structure, which could help in the timing deviation compensation design. Previous studies have shown that when a motor load current flows in the filter circuit and the inverter, the phase leg blanking times distort the voltage pulse sequence fed to the filter input. These blanking times are caused by excessively large dead time values between the IGBT control pulses. Moreover, the various switching timing distortions, present in realworld electronics when operating with a microsecond timescale, bring additional skew to the control. Left uncompensated, this results in distortion of the filter input voltage and a filter self-induced overvoltage in the form of an overshoot. This overshoot adds to the voltage appearing at the motor terminals, thus increasing the transient voltage amplitude at the motor. This doctoral thesis investigates the magnitude of such timing deviation effects. If the motor load current is left uncompensated in the control, the filter output voltage can overshoot up to double the input voltage amplitude. IGBT nonlinearities were observed to cause a smaller overshoot, in the order of 30%. This thesis introduces an improved ADUDT control method that is able to compensate for phase leg blanking times, giving flexibility to the pulse pattern structure and dead times. The control method is still sensitive to timing deviations, and their effect is investigated. A simple approach of using a fixed delay compensation value was tried in the test setup measurements. The ADUDT method with the new control algorithm was found to work in an actual motor drive application. Judging by the simulation results, with the delay compensation, the method should ultimately enable an output voltage performance and a du/dt reduction that are free from residual overshoot effects. The proposed control algorithm is not strictly required for successful ADUDT operation: It is possible to precalculate the pulse patterns by iteration and then for instance store them into a look-up table inside the control electronics. Rather, the newly developed control method is a mathematical tool for solving the ADUDT control pulses. It does not contain the timing deviation compensation (from the logic-level command to the phase leg output voltage), and as such is not able to remove the timing deviation effects that cause error and overshoot in the filter. When the timing deviation compensation has to be tuned-in in the control pattern, the precalculated iteration method could prove simpler and equally good (or even better) compared with the mathematical solution with a separate timing compensation module. One of the key findings in this thesis is the conclusion that the correctness of the pulse pattern structure, in the sense of ZCC and predicted pulse timings, cannot be separated from the timing deviations. The usefulness of the correctly calculated pattern is reduced by the voltage edge timing errors. The doctoral thesis provides an introductory background chapter on variable-speed AC drives and the problem of motor overvoltages and takes a look at traditional solutions for overvoltage mitigation. Previous results related to the active du/dt filtering are discussed. The basic operation principle and design of the filter have been studied previously. The effect of load current in the filter and the basic idea of compensation have been presented in the past. However, there was no direct way of including the dead time in the control (except for solving the pulse pattern manually by iteration), and the magnitude of nonlinearity effects had not been investigated. The enhanced control principle with the dead time handling capability and a case study of the test setup timing deviations are the main contributions of this doctoral thesis. The simulation and experimental setup results show that the proposed control method can be used in an actual drive. Loss measurements and a comparison of active du/dt output filtering with traditional output filtering methods are also presented in the work. Two different ADUDT filter designs are included, with ferrite core and air core inductors. Other filters included in the tests were a passive du/dtfilter and a passive sine filter. The loss measurements incorporated a silicon carbide diode-equipped IGBT module, and the results show lower losses with these new device technologies. The new control principle was measured in a 43 A load current motor drive system and was able to bring the filter output peak voltage from 980 V (the previous control principle) down to 680 V in a 540 V average DC link voltage variable-speed drive. A 200 m motor cable was used, and the filter losses for the active du/dt methods were 111W–126 W versus 184 W for the passive du/dt. In terms of inverter and filter losses, the active du/dt filtering method had a 1.82-fold increase in losses compared with an all-passive traditional du/dt output filter. The filter mass with the active du/dt method was 17% (2.4 kg, air-core inductors) compared with 14 kg of the passive du/dt method filter. Silicon carbide freewheeling diodes were found to reduce the inverter losses in the active du/dt filtering by 18% compared with the same IGBT module with silicon diodes. For a 200 m cable length, the average peak voltage at the motor terminals was 1050 V with no filter, 960 V for the all-passive du/dt filter, and 700 V for the active du/dt filtering applying the new control principle.

Resonant tunneling effects in semiconductor heterostructures

The thesis is devoted to a theoretical study of resonant tunneling phenomena in semiconductor heterostructures and nanostructures. It considers several problems relevant to modern solid state physics. Namely these are tunneling between 2D electron layers with spin-orbit interaction, tunnel injection into molecular solid material, resonant tunnel coupling of a bound state with continuum and resonant indirect exchange interaction mediated by a remote conducting channel. A manifestation of spin-orbit interaction in the tunneling between two 2D electron layers is considered. General expression is obtained for the tunneling current with account of Rashba and Dresselhaus types of spin-orbit interaction and elastic scattering. It is demonstrated that the tunneling conductance is very sensitive to relation between Rashba and Dresselhaus contributions and opens possibility to determine the spin-orbit interaction parameters and electron quantum lifetime in direct tunneling experiments with no external magnetic field applied. A microscopic mechanism of hole injection from metallic electrode into organic molecular solid (OMS) in high electric field is proposed for the case when the molecules ionization energy exceeds work function of the metal. It is shown that the main contribution to the injection current comes from direct isoenergetic transitions from localized states in OMS to empty states in the metal. Strong dependence of the injection current on applied voltage originates from variation of the number of empty states available in the metal rather than from distortion of the interface barrier. A theory of tunnel coupling between an impurity bound state and the 2D delocalized states in the quantum well (QW) is developed. The problem is formulated in terms of Anderson-Fano model as configuration interaction between the carrier bound state at the impurity and the continuum of delocalized states in the QW. An effect of this interaction on the interband optical transitions in the QW is analyzed. The results are discussed regarding the series of experiments on the GaAs structures with a -Mn layer. A new mechanism of ferromagnetism in diluted magnetic semiconductor heterosructures is considered, namely the resonant enhancement of indirect exchange interaction between paramagnetic centers via a spatially separated conducting channel. The underlying physical model is similar to the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction; however, an important difference relevant to the low-dimensional structures is a resonant hybridization of a bound state at the paramagnetic ion with the continuum of delocalized states in the conducting channel. An approach is developed, which unlike RKKY is not based on the perturbation theory and demonstrates that the resonant hybridization leads to a strong enhancement of the indirect exchange. This finding is discussed in the context of the known experimental data supporting the phenomenon.

Large-eddy simulation of wind flows over complex terrains for wind energy applications

Wind energy has obtained outstanding expectations due to risks of global warming and nuclear energy production plant accidents. Nowadays, wind farms are often constructed in areas of complex terrain. A potential wind farm location must have the site thoroughly surveyed and the wind climatology analyzed before installing any hardware. Therefore, modeling of Atmospheric Boundary Layer (ABL) flows over complex terrains containing, e.g. hills, forest, and lakes is of great interest in wind energy applications, as it can help in locating and optimizing the wind farms. Numerical modeling of wind flows using Computational Fluid Dynamics (CFD) has become a popular technique during the last few decades. Due to the inherent flow variability and large-scale unsteadiness typical in ABL flows in general and especially over complex terrains, the flow can be difficult to be predicted accurately enough by using the Reynolds-Averaged Navier-Stokes equations (RANS). Large- Eddy Simulation (LES) resolves the largest and thus most important turbulent eddies and models only the small-scale motions which are more universal than the large eddies and thus easier to model. Therefore, LES is expected to be more suitable for this kind of simulations although it is computationally more expensive than the RANS approach. With the fast development of computers and open-source CFD software during the recent years, the application of LES toward atmospheric flow is becoming increasingly common nowadays. The aim of the work is to simulate atmospheric flows over realistic and complex terrains by means of LES. Evaluation of potential in-land wind park locations will be the main application for these simulations. Development of the LES methodology to simulate the atmospheric flows over realistic terrains is reported in the thesis. The work also aims at validating the LES methodology at a real scale. In the thesis, LES are carried out for flow problems ranging from basic channel flows to real atmospheric flows over one of the most recent real-life complex terrain problems, the Bolund hill. All the simulations reported in the thesis are carried out using a new OpenFOAM® -based LES solver. The solver uses the 4th order time-accurate Runge-Kutta scheme and a fractional step method. Moreover, development of the LES methodology includes special attention to two boundary conditions: the upstream (inflow) and wall boundary conditions. The upstream boundary condition is generated by using the so-called recycling technique, in which the instantaneous flow properties are sampled on aplane downstream of the inlet and mapped back to the inlet at each time step. This technique develops the upstream boundary-layer flow together with the inflow turbulence without using any precursor simulation and thus within a single computational domain. The roughness of the terrain surface is modeled by implementing a new wall function into OpenFOAM® during the thesis work. Both, the recycling method and the newly implemented wall function, are validated for the channel flows at relatively high Reynolds number before applying them to the atmospheric flow applications. After validating the LES model over simple flows, the simulations are carried out for atmospheric boundary-layer flows over two types of hills: first, two-dimensional wind-tunnel hill profiles and second, the Bolund hill located in Roskilde Fjord, Denmark. For the twodimensional wind-tunnel hills, the study focuses on the overall flow behavior as a function of the hill slope. Moreover, the simulations are repeated using another wall function suitable for smooth surfaces, which already existed in OpenFOAM® , in order to study the sensitivity of the flow to the surface roughness in ABL flows. The simulated results obtained using the two wall functions are compared against the wind-tunnel measurements. It is shown that LES using the implemented wall function produces overall satisfactory results on the turbulent flow over the two-dimensional hills. The prediction of the flow separation and reattachment-length for the steeper hill is closer to the measurements than the other numerical studies reported in the past for the same hill geometry. The field measurement campaign performed over the Bolund hill provides the most recent field-experiment dataset for the mean flow and the turbulence properties. A number of research groups have simulated the wind flows over the Bolund hill. Due to the challenging features of the hill such as the almost vertical hill slope, it is considered as an ideal experimental test case for validating micro-scale CFD models for wind energy applications. In this work, the simulated results obtained for two wind directions are compared against the field measurements. It is shown that the present LES can reproduce the complex turbulent wind flow structures over a complicated terrain such as the Bolund hill. Especially, the present LES results show the best prediction of the turbulent kinetic energy with an average error of 24.1%, which is a 43% smaller than any other model results reported in the past for the Bolund case. Finally, the validated LES methodology is demonstrated to simulate the wind flow over the existing Muukko wind farm located in South-Eastern Finland. The simulation is carried out only for one wind direction and the results on the instantaneous and time-averaged wind speeds are briefly reported. The demonstration case is followed by discussions on the practical aspects of LES for the wind resource assessment over a realistic inland wind farm.

Espoon puhdistamolietteen käsittelymenetelmien selvitys Ämmässuon käsittelykeskukselle

Ravinteiden kierrätys yhä enemmän luonnonvaroja kuluttavassa yhteiskunnassa on noussut tärkeäksi. Hyötykäyttämällä luonnollisesti syntyvät orgaaniset virrat yhteiskunnasta voidaan kierrättää pellolta nostetut ravinteet takaisin. Tässä diplomityössä on tarkasteltu Helsingin seudun ympäristöpalveluiden toteuttamaa mädätetyn puhdistamolietteen käsittelyä ja pohdittu sen käsittelymahdollisuuksia erityisesti Ämmässuon käsittelykeskuksessa. Tarkastelussa on keskitytty kahteen yleisimpään kompostointimenetelmään, eli auma- ja tunnelikompostointiin. Kompostointimenetelmien tarkastelussa on pyritty löytämään optimaalinen tukiaine sekä prosessiaika käsiteltävälle lietteelle Ämmässuolla olemassa olleilla resursseilla. Tutkimuksessa toteutettiin käytännön pilot -kokeita sekä toteutettiin teoreettista prosessin parametrien laskentaa mukaan lukien kustannukset ja resurssit. Tutkimuksen perusteella voidaan todeta, että kyseisen tyyppisellä käsittelyalueella tunneli-kompostointi on näillä parametreilla testattuna paras vaihtoehto. Menetelmällä voidaan saavuttaa nopeampi käsittelyaika sekä vähentää kompostoinnista aiheutuvia hajuhaittoja ympäristöön. Tukiaineseoksessa voidaan korvata yleisestikin käytetty turve vastaanotettavalla lehti-puru viherjäteseoksella sekä risuhakkeella ja prosessiaika on noin 12 -16 viikkoa. Käsittelyn kustannukset verrattuna tämän hetkiseen käsittelyyn Metsäpirtissä eivät olleet merkittävästi suuremmat ja ympäristövaikutusarvioinnin perusteella syntypaikkaa lähellä oleva käsittelyalue sekä prosessin pienemmät vaikutukset tuovat ympäristöedun käsittelylle Ämmässuolla.

Modifications of the sympathetic skin response in workers chronically exposed to lead

The long-term effects of low-level lead intoxication are not known. The sympathetic skin response (SSR) was evaluated in a group of 60 former workers of a primary lead smelter, located in Santo Amaro, BA, Brazil. The individuals participating in the study were submitted to a clinical-epidemiological evaluation including questions related to potential risk factors for intoxication, complaints related to peripheral nervous system (PNS) involvement, neurological clinical examination, and also to electromyography and nerve conduction studies and SSR evaluation. The sample consisted of 57 men and 3 women aged 34 to 69 years (mean ± SD: 46.8 ± 6.9). The neurophysiologic evaluation showed the presence of lumbosacral radiculopathy in one of the individuals (1.7%), axonal sensorimotor polyneuropathy in 2 (3.3%), and carpal tunnel syndrome in 6 (10%). SSR was abnormal or absent in 12 cases, representing 20% of the sample. More than half of the subjects (53.3%) reported a history of acute abdominal pain requiring hospitalization during the period of work at the plant. A history of acute palsy of radial and peroneal nerves was reported by about 16.7 and 8.3% of the individuals, respectively. Mean SSR amplitude did not differ significantly between patients presenting or not the various characteristics in the current neurological situation, except for diaphoresis. The results suggest that chronic lead intoxication induces PNS damage, particularly affecting unmyelinated small fibers. Further systematic study is needed to more precisely define the role of lead in inducing PNS injury.

Vasa hamnvägs MKB : Bedömningsprogram för miljökonsekvenser

Det aktuella projektet ligger vid och i närheten av Södra Stadsfjärden i Vasa stad och Korsholms kommun. Målsättningen med Vasa hamnvägsprojektet är att finna en förbindelse till Vasa hamn som fungerar bättre en den nuvarande och i och med detta utveckla Vasa regionens logistiksystem. Avsikten är att skapa en vägförbindelse som effektivt förenar Vasklots hamn med logistikcentret och riksvägarna 3 och 8 som är en del av det riksomfattande huvudvägnätet. Vasa hamnvägsprojektet omfattar även en sammankoppling av logistikcentret och Lillkyrovägen (landsväg 717) som en del av förbindelsen Helsingby – Vassor av riksväg 8. Utvecklandet av Vasa hamnvägsförbindelsen har varit ett aktuellt och viktigt projekt redan länge. Bakgrunden till planerandet av vägen är de störningsfaktorer som förekommit i Vasa centrum av den genomgående lastbilstrafiken samt de problem som uppkommer för dessa stora transporter till följd av det trånga gatunätverket och broarnas begränsade bärförmåga. Vasa hamnvägs sträckningsalternativ har utretts omfattande i lokaliseringsutredning (2010) och till den fortsatta planeringen har tre vägkorridorsalternativ valts. Dessutom utreds anläggning av tunnel under Vasa centrum, och förbättring av befintlig förbindelse. (Alt 0+) Syftet med denna miljökonsekvensbedömning är att bedöma de miljökonsekvenser som byggandet av Vasa hamnväg orsakar. Bedömningen görs på ett sådant sätt och med en sådan noggrannhet som det förutsätts i MKB-lagen och -förordningen. De konsekvenser som kommer att bedömas: • människors hälsa, levnadsförhållanden och trivsel • jordmån, vatten, luft, klimat, växtlighet, organismer och naturens mångfald • Samhällsstrukturerna, byggnader, landskap, stadsbilden och kulturarvet • trafik och trafiksäkerhet • utnyttjande av naturresurser. Områdets invånare och andra intressenter har möjlighet att delta i planeringen och konsekvensbedömningen. Synpunkter om projektet och dess alternativ insamlas under hela planeringen. Invånare och andra intressenter informeras om projektet genom kungörelser, med hjälp av postningslista samt på internet.

Reduction of carcass weight loss in swines

The carcass fast freezing is one of the aspects of great prominence to the final quality of pork. In order to reduce weight loss, two experiments were performed, in which the carcasses were monitored during 20 hours to evaluate the main variables involved during two different freezing processes (standard and proposed) as follows: microbiological quality, storage temperature, relative humidity (RH) and air velocity. In experiment I, the carcasses were submitted to a system using heat shock (2 hours in static tunnel at - 25 °C) and subsequently sent to the equalization chamber. In experiment II, the carcasses were submitted to the heat shock and stored in a chamber with RH between 80-85%. The chambers used in both experiments showed no change in the variables studied (internal temperature of 5 °C and air velocity of approximately 0.3 m/s). However, the relative humidity in the three chambers was evaluated and significant differences were found; as a consequence, high levels of weight loss were observed in both chambers In experiment II there was an increase of RH, which reduced the weight loss of the carcasses.

Air Ducts

Air duct running along the utilities tunnel.

Early Construction

View of the service tunnel and foundation construction in October 1971.

The finite bandwidth model for spin fluctuations in Pd

The frequency dependence of the electron-spin fluctuation spectrum, P(Q), is calculated in the finite bandwidth model. We find that for Pd, which has a nearly full d-band, the magnitude, the range, and the peak frequency of P(Q) are greatly reduced from those in the standard spin fluctuation theory. The electron self-energy due to spin fluctuations is calculated within the finite bandwidth model. Vertex corrections are examined, and we find that Migdal's theorem is valid for spin fluctuations in the nearly full band. The conductance of a normal metal-insulator-normal metal tunnel junction is examined when spin fluctuations are present in one electrode. We find that for the nearly full band, the momentum independent self-energy due to spin fluctuations enters the expression for the tunneling conductance with approximately the same weight as the self-energy due to phonons. The effect of spin fluctuations on the tunneling conductance is slight within the finite bandwidth model for Pd. The effect of spin fluctuations on the tunneling conductance of a metal with a less full d-band than Pd may be more pronounced. However, in this case the tunneling conductance is not simply proportional to the self-energy.

Water-power at Niagara Falls to be successfully utilized : its objects, facilities and resources, inducements for manufacturers, etc. /

The lithograph, "General view of lands, tunnel and docks of Niagara River Hydraulic Tunnel, Power and Sewer Company," called for p. [4] in the Index, has been removed and encapsulated, and is shelved separately.

Niagara Bruce Trail Club Fonds, 1962-2014, n.d.

The Bruce trail is Canada’s longest and oldest continuous footpath. The trail runs along the Niagara Escarpment from Niagara to Tobermory through private and public land. The main trail is 890 km long and the side trails measure 400 km. In 1961, a “Save the Escarpment” conference was held in Hamilton. Gerry Wolfram, a writer for the St. Catharines Standard proposed that a committee be formed to develop a hiking trail. The Peninsula Field Naturalists Club formed a committee and President Bert Lowe contacted landowners along the proposed route to gain permission to cross their properties. Through Bert Lowe’s effort and dedication, the trail was completed in October 1963. The trail was officially opened on May 24th, 1964 in a ceremony at Queenston. The Niagara group joined the Bruce Trail Association in 1968 at which time the Niagara Bruce Trail Club was formed. The Bruce Trail Association is a charitable, membership-based volunteer organization. Their goal is to preserve public access to the Niagara Escarpment while restoring its natural habitat. The head office of the Bruce Trail Association is located in Hamilton, Ontario. The Niagara Bruce Trail Club’s goal is to secure and preserve a natural corridor along the Niagara Escarpment while providing education, awareness, and access for the public and the future. The club has organized many hikes including special hikes such as the one to commemorate the St. Catharines Centennial. The club has also hosted children’s hikes, cross country skiing hikes, wildflower hikes, jogging hikes, snowshoe hikes and bike outings. They hold annual events such as the End to End hike which is a 3 day walk from Grimsby to Queenston and the 30 km Laura Secord hike to commemorate Laura Secord’s famous walk. Charity hikes have also been held for the Heart and Stroke Foundation and the Lung Association as well as other causes. Major changes have taken place along the trail throughout the years, some of these include: a reroute which eliminated the tunnel passage (1976) and a bridge which eliminated the need to walk to Mountain Road to cross the Queen Elizabeth Way (2008). Other major changes and clean-up projects have been undertaken by the club. The Bruce Trail Conservancy (formerly Association) is made up of 9 clubs including: Niagara Bruce Trail Club (Queenston to Grimsby), Iroquia Bruce Trail Club (Grimsby to Kelso), Toronto Bruce Trail Club (Kelso to Cheltenham), Caledon Hills Bruce Trail Club (Cheltenham to Mono Centre), Dufferin Hi-Land Bruce Trail Club (Mono Centre to Lavender), Blue Mountains Bruce Trail Club (Lavender to Craigleath), Beaver Valley Bruce Trail Club (Craigleath to Blantyre), Sydenham Bruce Trail Club (Blantyre to Wiarton) and Peninsula Bruce Trail Club (Wiarton to Tobermory). Sources: http://www.niagarabrucetrail.org/index.html and http://brucetrail.org/