Modelling initial survival and growth of radiata pine in New Zealand.

A sensitive framework has been developed for modelling young radiata pine survival, its growth and its size class distribution, from time of planting to age 5 or 6 years. The data and analysis refer to the Central North Island region of New Zealand. The survival function is derived from a Weibull probability density function, to reflect diminishing mortality with the passage of time in young stands. An anamorphic family of trends was used, as very little between-tree competition can be expected in young stands. An exponential height function was found to fit best the lower portion of its sigmoid form. The most appropriate basal area/ha exponential function included an allometric adjustment which resulted in compatible mean height and basal area/ha models. Each of these equations successfully represented the effects of several establishment practices by making coefficients linear functions of site factors, management activities and their interactions. Height and diameter distribution modelling techniques that ensured compatibility with stand values were employed to represent the effects of management practices on crop variation. Model parameters for this research were estimated using data from site preparation experiments in the region and were tested with some independent data sets.

Nanoscale studies of materials using x-ray synchrotron radiation and mesoscopic modelling

X-ray synchrotron radiation was used to study the nanostructure of cellulose in Norway spruce stem wood and powders of cobalt nanoparticles in cellulose support. Furthermore, the growth of metallic clusters was modelled and simulated in the mesoscopic size scale. Norway spruce was characterized with x-ray microanalysis at beamline ID18F of the European Synchrotron Radiation Facility in Grenoble. The average dimensions and the orientation of cellulose crystallites was determined using x-ray microdiffraction. In addition, the nutrient element content was determined using x-ray fluorescence spectroscopy. Diffraction patterns and fluorescence spectra were simultaneously acquired. Cobalt nanoparticles in cellulose support were characterized with x-ray absorption spectroscopy at beamline X1 of the Deutsches Elektronen-Synchrotron in Hamburg, complemented by home lab experiments including x-ray diffraction, electron microscopy and measurement of magnetic properties with a vibrating sample magnetometer. Extended x-ray absorption fine structure spectroscopy (EXAFS) and x-ray diffraction were used to solve the atomic arrangement of the cobalt nanoparticles. Scanning- and transmission electron microscopy were used to image the surfaces of the cellulose fibrils, where the growth of nanoparticles takes place. The EXAFS experiment was complemented by computational coordination number calculations on ideal spherical nanocrystals. The growth process of metallic nanoclusters on cellulose matrix is assumed to be rather complicated, affected not only by the properties of the clusters themselves, but essentially depending on the cluster-fiber interfaces as well as the morphology of the fiber surfaces. The final favored average size for nanoclusters, if such exists, is most probably a consequence of these two competing tendencies towards size selection, one governed by pore sizes, the other by the cluster properties. In this thesis, a mesoscopic model for the growth of metallic nanoclusters on porous cellulose fiber (or inorganic) surfaces is developed. The first step in modelling was to evaluate the special case of how the growth proceeds on flat or wedged surfaces.

Mathematical modelling of groundwater pollution in a small heterogeneous aquifer at Kärkölä, southern Finland

Yhteenveto: Kärkölän likaantuneen pohjavesialueen matemaattinen mallinnus

Modelling vacuum arcs : from plasma initiation to surface interactions

A better understanding of vacuum arcs is desirable in many of today's 'big science' projects including linear colliders, fusion devices, and satellite systems. For the Compact Linear Collider (CLIC) design, radio-frequency (RF) breakdowns occurring in accelerating cavities influence efficiency optimisation and cost reduction issues. Studying vacuum arcs both theoretically as well as experimentally under well-defined and reproducible direct-current (DC) conditions is the first step towards exploring RF breakdowns. In this thesis, we have studied Cu DC vacuum arcs with a combination of experiments, a particle-in-cell (PIC) model of the arc plasma, and molecular dynamics (MD) simulations of the subsequent surface damaging mechanism. We have also developed the 2D Arc-PIC code and the physics model incorporated in it, especially for the purpose of modelling the plasma initiation in vacuum arcs. Assuming the presence of a field emitter at the cathode initially, we have identified the conditions for plasma formation and have studied the transitions from field emission stage to a fully developed arc. The 'footing' of the plasma is the cathode spot that supplies the arc continuously with particles; the high-density core of the plasma is located above this cathode spot. Our results have shown that once an arc plasma is initiated, and as long as energy is available, the arc is self-maintaining due to the plasma sheath that ensures enhanced field emission and sputtering. The plasma model can already give an estimate on how the time-to-breakdown changes with the neutral evaporation rate, which is yet to be determined by atomistic simulations. Due to the non-linearity of the problem, we have also performed a code-to-code comparison. The reproducibility of plasma behaviour and time-to-breakdown with independent codes increased confidence in the results presented here. Our MD simulations identified high-flux, high-energy ion bombardment as a possible mechanism forming the early-stage surface damage in vacuum arcs. In this mechanism, sputtering occurs mostly in clusters, as a consequence of overlapping heat spikes. Different-sized experimental and simulated craters were found to be self-similar with a crater depth-to-width ratio of about 0.23 (sim) - 0.26 (exp). Experiments, which we carried out to investigate the energy dependence of DC breakdown properties, point at an intrinsic connection between DC and RF scaling laws and suggest the possibility of accumulative effects influencing the field enhancement factor.