Production process of Scots pine.

The accompanying collective research report is the result of the research project in 1986­90 between The Finnish Academy and the former Soviet Academy of Sciences. The project was organized around common field work in Finland and in the former Soviet Union and theoretical analyses of tree growth determining processes. Based on theoretical analyses, dynamic stand growth models were made and their parameters were determined utilizing the field results. Annual cycle affects the tree growth. Our theoretical approach was based on adaptation to local climate conditions from Lapland to South Russia. The initiation of growth was described as a simple low and high temperature accumulation driven model. Linking the theoretical model with long term temperature data allowed us to analyze what type of temperature response produced favorable outcome in different climates. Initiation of growth consumes the carbohydrate reserves in plants. We measured the dynamics of insoluble and soluble sugars in the very northern and Karelian conditions. Clear cyclical pattern was observed but the differences between locations were surprisingly small. Analysis of field measurements of CO2 exchange showed that irradiance is the dominating factor causing variation in photosynthetic rate in natural conditions during summer. The effect of other factors is so small that they can be omitted without any considerable loss of accuracy. A special experiment carried out in Hyytiälä showed that the needle living space, defined as the ratio between the shoot cylindric volume and needle surface area, correlates with the shoot photosynthesis. The penetration of irradiance into Scots pine canopy is a complicated phenomenon because of the movement of the sun on the sky and the complicated structure of branches and needles. A moderately simple but balanced forest radiation regime submodel was constructed. It consists of the tree crown and forest structure, the gap probability calculation and the consideration of spatial and temporal variation of radiation inside the forest. The common field excursions in different geographical regions resulted in a lot of experimental data of regularities of woody structures. The water transport seems to be a good common factor to analyse these properties of tree structure. There are evident regressions between cross-sectional areas measured at different locations along the water pathway from fine roots to needles. The observed regressions have clear geographical trends. For example, the same cross-sectional area can support three times higher needle mass in South Russia than in Lapland. Geographical trends can also be seen in shoot and needle structure. Analysis of data published by several Russian authors show, that one ton of needles transpire 42 ton of water a year. This annual amount of transpiration seems to be independent of geographical location, year and site conditions. The produced theoretical and experimental material is utilised in the development of stand growth model that describes the growth and development of Scots pine stands in Finland and the former Soviet Union. The core of the model is carbon and nutrient balances. This means that carbon obtained in photosynthesis is consumed for growth and maintenance and nutrients are taken according to the metabolic needs. The annual photosynthetic production by trees in the stand is determined as a function of irradiance and shading during the active period. The utilisation of the annual photosynthetic production to the growth of different components of trees is based on structural regularities. Since the fundamental metabolic processes are the same in all locations the same growth model structure can be applied in the large range of Scots pine. The annual photosynthetic production and structural regularities determining the allocation of resources have geographical features. The common field measurements enable the application of the model to the analysis of growth and development of stands growing on the five locations of experiments. The model enables the analysis of geographical differences in the growth of Scots pine. For example, the annual photosynthetic production of a 100-year-old stand at Voronez is 3.5 times higher than in Lapland. The share consumed to needle growth (30 %) and to growth of branches (5 %) seems to be the same in all locations. In contrast, the share of fine roots is decreasing when moving from north to south. It is 20 % in Lapland, 15 % in Hyytiälä Central Finland and Kentjärvi Karelia and 15 % in Voronez South Russia. The stem masses (115­113 ton/ha) are rather similar in Hyytiälä, Kentjärvi and Voronez, but rather low (50 ton/ha) in Lapland. In Voronez the height of the trees reach 29 m being in Hyytiälä and Kentjärvi 22 m and in Lapland only 14 m. The present approach enables utilization of structural and functional knowledge, gained in places of intensive research, in the analysis of growth and development of any stand. This opens new possibilities for growth research and also for applications in forestry practice.

Design, modelling and simulation of vibratory micromachined gyroscopes

Among various MEMS sensors, a rate gyroscope is one of the most complex sensors from the design point of view. The gyro normally consists of a proof mass suspended by an elaborate assembly of beams that allow the system to vibrate in two transverse modes. The structure is normally analysed and designed using commercial FEM packages such as ANSYS or MEMS specific commercial tools such as Coventor or Intellisuite. In either case, the complexity in analysis rises manyfolds when one considers the etch hole topography and the associated fluid flow calculation for damping. In most cases, the FEM analysis becomes prohibitive and one resorts to equivalent electrical circuit simulations using tools like SABER in Coventor. Here, we present a simplified lumped parameter model of the tuning fork gyro and show how easily it can be implemented using a generic tool like SIMULINK. The results obtained are compared with those obtained from more elaborate and intense simulations in Coventor. The comparison shows that lumped parameter SIMULINK model gives equally good results with fractional effort in modelling and computation. Next, the performance of a symmetric and decoupled vibratory gyroscope structure is also evaluated using this approach and a few modifications are made in this design to enhance the sensitivity of the device.

3D Acoustic Analysis Of Spherical Chamber Having Single Inlet And Multiple Outlet: An Impedance Matrix Approach

The transmission loss (TL) performance of spherical chambers having single inlet and multiple outlet is obtained analytically through modal expansion of acoustic field inside the spherical cavity in terms of the spherical Bessel functions and Legendre polynomials. The uniform piston driven model based upon the impedance [Z] matrix is used to characterize the multi-port spherical chamber. It is shown analytically that the [Z] parameters are independent of the azimuthal angle (phi) due to the axisymmetric shape of the sphere; rather, they depend only upon the polar angle (theta) and radius of the chamber R(0). Thus, the effects of relative polar angular location of the ports and number of outlet ports are investigated. The analytical results are shown to be in good agreement with the 3D FEA results, thereby validating the procedure suggested in this work.

3-D acoustic analysis of elliptical chamber mufflers having an end-inlet and a side-outlet: An impedance matrix approach

The acoustical behavior of an elliptical chamber muffler having an end-inlet and side-outlet port is analyzed semi-analytically. A uniform piston source is assumed to model the 3-D acoustic field in the elliptical chamber cavity. Towards this end, we consider the modal expansion of acoustic pressure field in the elliptical cavity in terms of angular and radial Mathieu functions, subjected to rigid wall condition, whereupon under the assumption of a point source, Green's function is obtained. On integrating this function over piston area of the side or end port and dividing it by piston area, one obtains the acoustic field, whence one can find the impedance matrix parameters characterizing the 2-port system. The acoustic performance of these configurations is evaluated in terms of transmission loss (TL). The analytical results thus obtained are compared with 3-D HA carried on a commercial software for certain muffler configurations. These show excellent agreement, thereby validating the 3-D semi-analytical piston driven model. The influence of the chamber length as well as the angular and axial location of the end and side ports on TL performance is also discussed, thus providing useful guidelines to the muffler designer. (c) 2011 Elsevier B.V. All rights reserved.

Acoustical behavior of single inlet and multiple outlet elliptical cylindrical chamber muffler

Transmission loss (TL) of an elliptical cylindrical chamber muffler having a single side/end inlet and multiple side/end outlet is analyzed by means of the 3-D semi-analytical formulation based upon the modal expansion (in terms of the angular and radial Mathieu functions) and the Green's function. The acoustic pressure response obtained in terms of Green's function is integrated over surface area of the side/end ports (modeled as rigid pistons) and upon subsequent division by the port area, yields the acoustic pressure response or impedance Z] matrix parameters due to the uniform piston-driven model. The 3-D semi-analytical results are found to be in excellent agreement with the results obtained by means of 3-D FEA (SYSNOISE) simulations, thereby validating the semi-analytical procedure suggested in this work. Parametric studies such as the effect of chamber length (L), angular and axial locations of the ports, interchanging the locations of inlet and outlet ports as well as the addition of an outlet port for double outlet mufflers on the TL performance are reported, thereby leading to the formulation of design guidelines for obtaining muffler configurations exhibiting a broad-band TL spectrum. One such configuration is an axially long chamber having side-inlet and side-outlet ports such that one of the side ports is located at half the axial length on themajor/minor axis and the other side port is located at three-quarters (or one-quarter) of the axial length on the minor/major axis. (C) 2012 Institute of Noise Control Engineering.

Analysis of distance relay trip regions for EHV and UHV transmission lines

This paper presents comparative evaluation of the distance relay characteristics for UHV and EHV transmission lines. Distance protection relay characteristics for the EHV and UHV systems are developed using Electromagnetic Transients (EMT) program. The variation of ideal trip boundaries for both the systems are presented. Unlike the conventional distance protection relay which uses a lumped parameter model, this paper uses the distributed parameter model. The effect of larger shunt susceptance on the trip boundaries is highlighted. Performance of distance relay with ideal trip boundaries for EHV and UHV lines have been tested for various fault locations and fault resistances. Electromagnetic Transients (EMT) program has been developed considering distributed parameter line model for simulating the test systems. The voltage and current phasors are computed from the signals using an improved full cycle DFT algorithm taking 20 samples per cycle. Two practical transmission systems of Indian power grid, namely 765 kV UHV transmission line and SREB 24-bus 400kV EHV system are used to test the performance of the proposed approach.

Support vector machine based high speed protection relay for EHV/UHV transmission line

This paper presents a fast and accurate relaying technique for a long 765kv UHV transmission line based on support vector machine. For a long EHV/UHV transmission line with large distributed capacitance, a traditional distance relay which uses a lumped parameter model of the transmission line can cause malfunction of the relay. With a frequency of 1kHz, 1/4th cycle of instantaneous values of currents and voltages of all phases at the relying end are fed to Support Vector Machine(SVM). The SVM detects fault type accurately using 3 milliseconds of post-fault data and reduces the fault clearing time which improves the system stability and power transfer capability. The performance of relaying scheme has been checked with a typical 765kV Indian transmission System which is simulated using the Electromagnetic Transients Program(EMTP) developed by authors in which the distributed parameter line model is used. More than 15,000 different short circuit fault cases are simulated by varying fault location, fault impedance, fault incidence angle and fault type to train the SVM for high speed accurate relaying. Simulation studies have shown that the proposed relay provides fast and accurate protection irrespective of fault location, fault impedance, incidence time of fault and fault type. And also the proposed scheme can be used as augmentation for the existing relaying, particularly for Zone-2, Zone-3 protection.

Remoulded shear strength at plastic and semi-solid states

The undrained shear strength of remoulded soils is of great concern in geotechnical engineering applications. This study aims to develop a reliable approach for determining the undrained shear strength of remoulded fine-grained soils, through the use of index test results, at both the plastic and semi-solid states of consistency. Experimental investigation and subsequent analysis involving a number of fine-grained soils of widely varying plasticity and geological origin have led to a two-parameter linear model of the relationship between logarithm of remoulded undrained shear strength and liquidity index. The numerical values of the parameters are found to be dependent to a lesser extent on the soil group and to a greater extent on the soil state. Based on the values of regression coefficient, ranking index and ranking distance, it seems that the relationship represents the experimental results well. It may be pointed out that the possibility of such a relationship in the semi-solid state of a soil has not been explored in the past. It is also shown that the shear strength at the plastic limit is about 32-34 times that at the liquid limit.

Discussion: Remoulded shear strength at plastic and semi-solid states

The undrained shear strength of remoulded soils is of great concern in geotechnical engineering applications. This study aims to develop a reliable approach for determining the undrained shear strength of remoulded fine-grained soils, through the use of index test results, at both the plastic and semi-solid states of consistency. Experimental investigation and subsequent analysis involving a number of fine-grained soils of widely varying plasticity and geological origin have led to a two-parameter linear model of the relationship between logarithm of remoulded undrained shear strength and liquidity index. The numerical values of the parameters are found to be dependent to a lesser extent on the soil group and to a greater extent on the soil state. Based on the values of regression coefficient, ranking index and ranking distance, it seems that the relationship represents the experimental results well. It may be pointed out that the possibility of such a relationship in the semi-solid state of a soil has not been explored in the past. It is also shown that the shear strength at the plastic limit is about 32–34 times that at the liquid limit.

Influence of exhaust gas heating and L/D ratios on the discharge efficiencies for an activated carbon natural gas storage system

A transient 2D axi-symmetric and lumped parameter (LP) model with constant outflow conditions have been developed to study the discharge capacity of an activated carbon bed. The predicted discharge times and variations in bed pressure and temperature are in good agreement with experimental results obtained from a 1.82 l adsorbed natural gas (ANG) storage system. Under ambient air conditions, a maximum temperature drop of 29.5 K and 45.5 K are predicted at the bed center for discharge rates of 1.0 l min(-1) and 5.0 l min(-1) respectively. The corresponding discharge efficiencies are 77% and 71.5% respectively with discharge efficiencies improving with decreasing outflow rates. Increasing the LID ratio from 1.9 to 7.8 had only a marginal increase in the discharge efficiency. Forced convection (exhaust gas) heating had a significant effect on the discharge efficiency, leading to efficiencies as high as 92.8% at a discharge of 1.0 l min(-1) and 88.7% at 5 l min(-1). Our study shows that the LP model can be reliably used to obtain discharge times due to the uniform pressure distributions in the bed. Temperature predictions with the LP model were more accurate at ambient conditions and higher discharge rates, due to greater uniformity in bed temperatures. For the low thermal conductivity carbon porous beds, our study shows that exhaust gas heating can be used as an effective and convenient strategy to improve the discharge characteristics and performance of an ANG system. (C) 2013 Elsevier Ltd. All rights reserved.

Robust Sliding Mode Control for Tokamaks

Nuclear fusion has arisen as an alternative energy to avoid carbon dioxide emissions, being the tokamak a promising nuclear fusion reactor that uses a magnetic field to confine plasma in the shape of a torus. However, different kinds of magnetohydrodynamic instabilities may affect tokamak plasma equilibrium, causing severe reduction of particle confinement and leading to plasma disruptions. In this sense, numerous efforts and resources have been devoted to seeking solutions for the different plasma control problems so as to avoid energy confinement time decrements in these devices. In particular, since the growth rate of the vertical instability increases with the internal inductance, lowering the internal inductance is a fundamental issue to address for the elongated plasmas employed within the advanced tokamaks currently under development. In this sense, this paper introduces a lumped parameter numerical model of the tokamak in order to design a novel robust sliding mode controller for the internal inductance using the transformer primary coil as actuator.

Knowledge modelling in aerospace design

This paper addresses the need for computer support in aerospace design. A review of current design methodologies and computer support tools is presented and the need for further support in aerospace design, particularly in the early formative stages of the design process, is discussed. A parameter-based model of design, founded on the assumption that a design process can be constructed from a predefined set of tasks, is proposed for aerospace design. This is supported by knowledge of possible tasks in which the confidence in key design parameters is used as a basis for identifying, or signposting, the next task. A prototype implementation of the signposting model, for use in the design of helicopter rotor blades, is described and results from trials of the tool are presented. Further areas of research are discussed

Studies in nuclear reactor dynamics : I. The accuracy of point kinetics. II. The effect of delayed neutrons on the spectrum of the group-diffusion operator

This thesis is a theoretical work on the space-time dynamic behavior of a nuclear reactor without feedback. Diffusion theory with G-energy groups is used.

In the first part the accuracy of the point kinetics (lumped-parameter description) model is examined. The fundamental approximation of this model is the splitting of the neutron density into a product of a known function of space and an unknown function of time; then the properties of the system can be averaged in space through the use of appropriate weighting functions; as a result a set of ordinary differential equations is obtained for the description of time behavior. It is clear that changes of the shape of the neutron-density distribution due to space-dependent perturbations are neglected. This results to an error in the eigenvalues and it is to this error that bounds are derived. This is done by using the method of weighted residuals to reduce the original eigenvalue problem to that of a real asymmetric matrix. Then Gershgorin-type theorems .are used to find discs in the complex plane in which the eigenvalues are contained. The radii of the discs depend on the perturbation in a simple manner.

In the second part the effect of delayed neutrons on the eigenvalues of the group-diffusion operator is examined. The delayed neutrons cause a shifting of the prompt-neutron eigenvalue s and the appearance of the delayed eigenvalues. Using a simple perturbation method this shifting is calculated and the delayed eigenvalues are predicted with good accuracy.

Micromechanical damage and fracture in elastomeric polymers

This thesis aims at a simple one-parameter macroscopic model of distributed damage and fracture of polymers that is amenable to a straightforward and efficient numerical implementation. The failure model is motivated by post-mortem fractographic observations of void nucleation, growth and coalescence in polyurea stretched to failure, and accounts for the specific fracture energy per unit area attendant to rupture of the material.

Furthermore, it is shown that the macroscopic model can be rigorously derived, in the sense of optimal scaling, from a micromechanical model of chain elasticity and failure regularized by means of fractional strain-gradient elasticity. Optimal scaling laws that supply a link between the single parameter of the macroscopic model, namely the critical energy-release rate of the material, and micromechanical parameters pertaining to the elasticity and strength of the polymer chains, and to the strain-gradient elasticity regularization, are derived. Based on optimal scaling laws, it is shown how the critical energy-release rate of specific materials can be determined from test data. In addition, the scope and fidelity of the model is demonstrated by means of an example of application, namely Taylor-impact experiments of polyurea rods. Hereby, optimal transportation meshfree approximation schemes using maximum-entropy interpolation functions are employed.

Finally, a different crazing model using full derivatives of the deformation gradient and a core cut-off is presented, along with a numerical non-local regularization model. The numerical model takes into account higher-order deformation gradients in a finite element framework. It is shown how the introduction of non-locality into the model stabilizes the effect of strain localization to small volumes in materials undergoing softening. From an investigation of craze formation in the limit of large deformations, convergence studies verifying scaling properties of both local- and non-local energy contributions are presented.

Thermal modelling of a tubular linear machine for marine renewable generation

A lumped parameter thermal model has been constructed for a tubular linear machine that has been designed for use in a marine environment. It shows good correlation to both steady state and transient experimental tests on the machine. The model has been developed for a stationary machine in a laboratory environment - the modelling techniques used and enhancements to enable the application of the model directly to marine scenarios are discussed.