Soil carbon (C) is important because even small changes in soil C can affect atmospheric concentrations of CO₂, which in turn can influence global climate. Adequate soil carbon is also required to maintain soil quality, which is important to if agricultural production is to be sustained. The soil carbon balance of New Zealand's pastoral soils is poorly understood, with recent research showing that soils under dairy pasture have lost large amounts of C during the past few decades. The main objective of this research was to determine an annual farm scale C budget for an intensively grazed dairy farm, with a second objective being to determine the amount of CO₂-C lost following cultivation for pasture renewal, and soil pugging by dairy cattle. A third objective was to investigate the environmental controls of CO₂ exchange in a dairy farm pasture system. Net ecosystem exchange (NEE) of CO₂ was measured using an eddy covariance (EC) system from 15 December 2007 to 14 December 2008. Closed chamber techniques were used to measure CO₂ emissions from three cultivated paddocks and three adjacent pasture paddocks between 26 January 2008 and 5 March 2008. CO₂ emissions were also measured using chambers from pugged and control plots between 25 June and 5 August. Coincidentally this research was carried out in a year with a severe summer/autumn drought and a wetter than usual winter. Annual NEE measured with the eddy covariance system was -1,843 kg C ha⁻¹ (a C gain by the land surface). Accounting for C in supplement import, milk export, pasture export and losses in methane, the dairy pasture system was a net sink of -880±500 kg C ha⁻¹. This C sequestration occurred despite severe drought during the study, which was in contrast to other studies of grasslands during drought. Cultivation under dry conditions did not increase cumulative CO₂-C emissions compared to adjacent pasture paddocks. However, when C inputs to pasture paddocks via photosynthesis were included in calculations, net C loss from the cultivated paddocks (during the 39 day study) was estimated to be 622 kg C ha⁻¹ more than the pasture paddocks. CO₂ emissions were lower from pugged plots compared to control plots, probably caused by decreased microbial and root respiration due to wetter soil conditions, and lowered root respiration as a result of lower pasture production. Volumetric soil moisture content (soil moisture) had a dominant effect on CO₂ exchange at a range of temporal scales. Respiration and photosynthesis were both reduced when soil moisture was below 43% (~the lower limit of readily available water) and photosynthesis virtually ceased when soil moisture declined below 24% (~wilting point). Soil moisture also influenced the relationship between temperature and respiration and photosynthetic flux density (PPFD) and NEE. These results suggest that management related soil disturbances of occasional cultivation for pasture renewal and soil pugging, are unlikely to cause large losses of soil C. Further, a severe drought also did not cause CO₂-C losses from the land surface to the atmosphere on an annual scale, in contrast to previous studies.

Graduation date: 1991

Turbulent flows present structures with a wide range of scales. The computation of the complete physics of a turbulent flow (termed DNS) is very expensive and is, for the time being, limited to low and medium Reynolds number flows. As a way to capture high Reynolds number flows, a part of the physics complexity has to be modeled. Large eddy simulation (LES) is a simulation strategy where the large turbulent eddies present on a given mesh are captured and the influence of the non-resolved scales onto the resolved ones is modeled. The present thesis reports on the development and validation of a methodology in order to apply LES for complex wall-bounded flows. Discretization methods and LES models, termed subgrid scale models (SGS), compatible with such a geometrical complexity are discussed. It is proved that discrete a kinetic energy conserving discretization of the convective term is an attractive solution to perform stable simulations without the use of an artificial dissipation, as upwinding. The dissipative effect of the SGS model is thus unaffected by any additional dissipation process. The methodology is first applied to a developed parallel fourth-order incompressible flow solver for cartesian non-uniform meshes. In order to solve the resulting Poisson equation, an efficient multigrid solver is also developed. The code is first validated using DNS (Taylor-Green vortex, channel flow, four-vortex system) and LES (channel flow), and finally applied to the investigation of an aircraft two-vortex system in ground effect. The methodology is then applied to improve a RANS-based industrial unstructured compressible flow solver, developed at CENAERO, to perform well for LES applications. The proposed modifications are tested successfully on the unsteady flow past a sphere at Reynolds of 300 and 10000, corresponding to the subcritical regime.

Road and rail vehicles and aircraft rely mainly or solely on friction brakes. These brakes pose several problems, especially in hybrid vehicles: significant wear, fading, complex and slow actuation, lack of fail-safe features, increased fuel consumption due to power assistance, and requirement for anti-lock controls. To solve these problems, a contactless magnetic brake has been developed. This concept includes a novel flux-shunting structure to control the excitation flux generated by permanent magnets. This brake is wear-free, less-sensitive to temperature than friction brakes, has fast and simple actuation, and has a reduced sensitivity to wheel-lock. The present dissertation includes an introduction to friction braking, a theory of eddy-current braking, analytical and numerical models of the eddy-current brake, its excitation and power generation, record of experimental validation, investigation and simulation of the integration of the brake in conventional and hybrid vehicles.

Turbulent flows present structures with a wide range of scales. The computation of the complete physics of a turbulent flow (termed DNS) is very expensive and is, for the time being, limited to low and medium Reynolds number flows. As a way to capture high Reynolds number flows, a part of the physics complexity has to be modeled. Large eddy simulation (LES) is a simulation strategy where the large turbulent eddies present on a given mesh are captured and the influence of the non-resolved scales onto the resolved ones is modeled. The present thesis reports on the development and validation of a methodology in order to apply LES for complex wall-bounded flows. Discretization methods and LES models, termed subgrid scale models (SGS), compatible with such a geometrical complexity are discussed. It is proved that discrete a kinetic energy conserving discretization of the convective term is an attractive solution to perform stable simulations without the use of an artificial dissipation, as upwinding. The dissipative effect of the SGS model is thus unaffected by any additional dissipation process. The methodology is first applied to a developed parallel fourth-order incompressible flow solver for cartesian non-uniform meshes. In order to solve the resulting Poisson equation, an efficient multigrid solver is also developed. The code is first validated using DNS (Taylor-Green vortex, channel flow, four-vortex system) and LES (channel flow), and finally applied to the investigation of an aircraft two-vortex system in ground effect. The methodology is then applied to improve a RANS-based industrial unstructured compressible flow solver, developed at CENAERO, to perform well for LES applications. The proposed modifications are tested successfully on the unsteady flow past a sphere at Reynolds of 300 and 10000, corresponding to the subcritical regime.

No description available.

The average and fluctuating components in a model ring-type gas turbine combustor are characterized using a Large Eddy Simulation at a Reynolds number of 11,000, based on the bulk velocity and the mean channel height. A spatial filter is applied to the incompressible Navier-Stokes equations, and a high pass filtered Smagorinsky model is used to model the sub-grid scales. Two cases are studied: one with only the swirler inlet active, and one with a single row of dilution jets activated, operating at a momentum flux ratio J of 100. The goal of both of these studies is to validate the capabilities of the solver NEK5000 to resolve important flow features inherent to gas turbine combustors by comparing qualitatively to the work of Jakirlic. Both cases show strong evidence of the Precessing Vortex Core, an essential flow feature in gas turbine combustors. Each case captures other important flow characteristics, such as corner eddies, and in general predicts bulk flow movements well. However, the simulations performed quite poorly in terms of predicting turbulence shear stress quantities. Difficulties in properly emulating the turbulent velocity entering the combustor for the swirl, as well as mesh quality concerns, may have skewed the results. Overall, though small length scale quantities were not accurately captured, the large scale quantities were, and this stress test on the HPF LES model will be built upon in future work that looks at more complex combustors.

In the context of Large Eddy Simulation (LES) of turbulent flows, there is a current need to compare and evaluate different proposed subfilter-scale models. In order to carefully compare subfilter-scale models and compare LES predictions to Direct Numerical Simulation (DNS) results (the latter would be helpful in the comparison and validation of models), there is a real need for a "grid-independent" LES capability and explicit filtering methods offer one means by which this may be achieved. Advantages of explicit filtering are that it provides a means for eliminating aliasing errors, allows for the direct control of commutation errors, and most importantly allows a decoupling between the mesh spacing and the filter width which is the primary reason why there are difficulties in comparing LES solutions obtained on different grids. This thesis considers the design and assessment of discrete explicit filters and their application to isotropic turbulence prediction.

In the context of Large Eddy Simulation (LES) of turbulent flows, there is a current need to compare and evaluate different proposed subfilter-scale models. In order to carefully compare subfilter-scale models and compare LES predictions to Direct Numerical Simulation (DNS) results (the latter would be helpful in the comparison and validation of models), there is a real need for a "grid-independent" LES capability and explicit filtering methods offer one means by which this may be achieved. Advantages of explicit filtering are that it provides a means for eliminating aliasing errors, allows for the direct control of commutation errors, and most importantly allows a decoupling between the mesh spacing and the filter width which is the primary reason why there are difficulties in comparing LES solutions obtained on different grids. This thesis considers the design and assessment of discrete explicit filters and their application to isotropic turbulence prediction.

Large Eddy Simulation is a method of obtaining high accuracy computational results for modelling fluid flow. Unfortunately it is computationally expensive limiting it to users of large parallel machines. However, it may be that the use of LES leads to an over-resolution of the problem because the bulk of the computational domain could be adequately modelled using the Reynolds averaged approach. A study has been undertaken to assess the feasibility, both in accuracy and computational efficiency of using a parallel computer to solve both LES and RANS type turbulence models on the same domain for the problem flow over a circular cylinder at Reynolds number 3 900 To do this the domain has been created and then divided into two sub-domains, one for the LES model and one for the kappa - epsilon turbulence model. The hybrid model has been developed specifically for a parallel computing environment and the user is able to allocate modelling techniques to processors in a way which enables expansion of the model to any number of processors. Computational experimentation has shown that the combination of the Smagorinsky model can be used to capture the vortex shedding from the cylinder and the information successfully passed to the kappa - epsilon model for the dissipation of the vortices further downstream. The results have been compared to high accuracy LES results and with both kappa - epsilon and Smagorinsky LES computations on the same domain. The hybrid models developed compare well with the Smagorinsky model capturing the vortex shedding with the correct periodicity. Suggestions for future work have been made to develop this idea further, and to investigate the possibility of using the technology for the modelling of mixing and fast chemical reactions based on the more accurate prediction of the turbulence levels in the LES sub-domain.