Closed-loop ground source heat pumps (GSHPs) are used to transfer thermal energy between the subsurface and conditioned spaces for heating and cooling applications. A basic GSHP is composed of a ground heat exchanger (GHX), which is a closed loop of pipe buried in the shallow subsurface circulating a heat exchange fluid, connected to a heat pump. These systems offer an energy efficient alternative to conventional heating and cooling systems; however, installation costs are higher due to the additional cost associated with the GHX. By further developing our understanding of how these ground loops interact with the subsurface, it may possible to design them more intelligently, efficiently, and economically. To gain insight into the physical processes occurring between the GHX and the subsurface and to identify efficiencies and inefficiencies in GSHP design and operation, two main research goals were defined: comprehensive monitoring of a fully functioning GSHP and intensive simulation of these systems using computer models. A 6-ton GSHP was installed at a residence in Elora, ON. An array of 64 temperature sensors was installed on and surrounding the GHX and power consumption and temperature sensors were installed on the system inside the residence. The data collected were used to help characterize and understand the function of the system, provide motivation for further investigations, and assess the impact of the time of use billing scheme on GSHP operation costs. To simulate GSHPs, two computer models were utilized. A 3D finite element model was employed to analyse the effects of pipe configuration and pipe spacing on system performance. A unique, transient 1D finite difference heat conduction model was developed to simulate a single pipe in a U-tube shape with inter-pipe interactions and was benchmarked against a tested analytical solution. The model was used to compare quasi-steady state and transient simulation of GSHPs, identify system performance efficiencies through pump schedule optimization, and investigate the effect of pipe length on system performance. A comprehensive comparison of steady state and pulsed simulation concludes that it is possible to simulate transient operation using a steady state assumption for some cases. Optimal pipe configurations are identified for a range of soil thermal properties. Optimized pump schedules are identified and analysed for a specific heat pump and fluid circulation pump. Finally, the effect of pipe spacing and length on system performance is characterized. It was found that there are few design inefficiencies that could be easily addressed to improve general design practice.

This thesis presents an axisymmetric numerical model to simulate interfacial flows near a sharp corner, where contact line pinning occurs. The method has been used to analyze drop imbibition into a capillary. To evaluate the performance of the numerical method, for a liquid drop initially placed partially within a capillary, a thermodynamic model has also been developed, to predict equilibrium states. The first part of this thesis presents an axisymmetric VoF algorithm to simulate interfacial flows near a sharp corner. (1) A new method to exactly calculate the normals and curvatures of any circle with a radius as small as the grid size is presented. This method is a hybrid least squares height function technique which fits a discretized osculating circle to a curve, from which interface normals and curvature can be evaluated. (2) A novel technique for applying the contact angle boundary condition has been devised, based on the definition of an osculating circle near a solid phase. (3) A new flux volume construction technique is presented, which can be applied to any split advection scheme. Unlike the traditional approach where the flux volumes are assumed rectangular, the new flux volumes can be either trapezoidal or triangular. The new technique improves the accuracy and consistency of the advection scheme. (4) Explicit PLIC reconstruction expressions for axisymmetric coordinates have been derived. (5) Finally, a numerical treatment of VoF for contact line motion near a sharp corner is presented, base on the idea of contact line pinning and an edge contact angle. The second part of the thesis is on the imbibition of a drop into an open capillary. A thermodynamic analysis based on minimization of an interfacial surface energy function is presented to predict equilibrium configurations of drops. Based on the drop size compared to the hole size, the equilibrium contact angle, and the geometry of the capillary, the drop can be totally imbibed by the capillary, or may not wet the capillary at all. The thesis concludes with application of the numerical scheme to the same problem, to examine the dynamics of wetting or dewetting of a capillary. All of the simulations yield results that correspond to the equilibrium states predicted by the thermodynamic analysis, but offer additional insight on contact line motion and interface deformation near the capillary edge.

In this thesis, a numerical model for coupled heat and moisture transfer in a run around membrane energy exchanger (RAMEE) with a liquid desiccant as a coupling fluid is developed. The numerical model is two dimensional, transient and is formulated using the finite difference method with an implicit time discretization. The model for the case of only heat transfer for a single heat exchanger is compared to an available analytical solution and good agreement is obtained. It is shown that the discrepancy between the numerical and theoretical dimensionless bulk outlet temperature of the fluids is less than 4% during the transient period. The model is also validated for the case of simultaneous heat and moisture transfer using experimental data measured during the laboratory testing of a RAMEE system. The results for both sensible and latent effectiveness showed satisfactory agreement at different operating conditions. However, there are some discrepancies between the simulation and the experimental data during the transient times. It is proposed that these discrepancies may be due to experimental flow distribution problems within the exchanger. The maximum average absolute differences between the measured and simulated transient effectivenesses were 7.5% and 10.3% for summer and winter operating conditions, respectively.<p> The transient response of the RAMEE system for step changes in the inlet supply air temperature and humidity ratio is presented using the numerical model. In addition, the system quasi steady state operating conditions are predicted as the system approaches its steady state operating condition. The effect of various dimensionless parameters on the transient response is predicted separately. These included: the number of heat transfer units, thermal capacity ratio, heat loss/gain ratio, storage volume ratio and the normalized initial salt solution concentration. It is shown that the initial salt solution concentration and the storage volume of the salt solution have significant impacts on the transient response of the system and the heat loss/gain rates from/to the circulated fluid flow can change the system quasi steady effectiveness substantially. The detailed study of the transient performance of the RAMEE is useful to determine the transient response time of the system under different practical situations.

In this study, a numerical model is developed for determining coupled heat and moisture transfer in a run-around membrane energy exchanger (RAMEE) using two counter/cross flow exchangers and with a salt solution of MgCl2 as the coupling fluid. The counter/cross flow exchanger is a counter-flow exchanger with cross-flow inlet and outlet headers. The model is two-dimensional, steady-state and based on the physical principles of conservation of momentum, energy, and mass. The finite difference method is used in this model to discretize the governing equations.<p> The heat transfer model is validated with effectiveness correlations in the literature. It is shown that the difference between the numerical model and correlations is less than ¡À2% and ¡À2.5% for heat exchangers and run around heat exchangers (RAHE), respectively. The simultaneous heat and moisture transfer model is validated with data from another model and experiments. The inter-model comparison shows a difference of less than 1%. The experimental validation shows an average discrepancy of 1% to 17% between the experimental and numerical data for overall total effectiveness. At lower NTUs the numerical and experimental results show better agreement (e.g. within 1-4% at NTU=4).<p> The model for RAHE is used to develop new effectiveness correlations for the geometrically more complex counter/cross flow heat exchangers and RAHE systems. The correlations are developed to predict the response of the exchangers and overall system to the change of different design characteristics as it is determined by the model. Discrepancies between the simulated and correlated results are within ¡À2% for both the heat exchangers and the RAHE systems.<p> It is revealed by the model that the overall effectiveness of the counter/cross flow RAMEE depends on the entrance ratio (the ratio of the length of the inlet and outlet headers to the length of the exchanger, xi/x0), aspect ratio (the ratio of the height to the length of the exchanger, y0/x0), number of heat transfer units (NTU), heat capacity rate ratio (Cr*), number of mass transfer units (NTUm), and the mass flow rate ratio of pure salt in desiccant solution to dry air (m*). Beside these dimensionless parameters, the performance of the RAMEE system is affected by the liquid-air flow configuration and the operating inlet temperature and humidity.<p> This study concludes that the maximum effectiveness of the RAMEE system with two counter/cross flow exchangers occurs when NTU and NTUm are large (e.g. greater than 10). At any NTU, the overall effectiveness of the RAMEE system increases with Cr* until it reaches a maximum value when Cr*= . Increasing Cr* above causes the overall effectiveness to decrease slightly. Therefore, to achieve the maximum overall effectiveness of the system, Cr* must be close to . is a function of NTU and operating conditions e.g., with NTU=10, and under AHRI summer and winter operating conditions, respectively. The exchangers in the RAMEE system are needed to have a small aspect ratio (e.g. y0/x0<0.2) and small entrance ratio (e.g. xi/x0<0.1) to get the maximum overall effectiveness of a RAMEE system using two counter/cross flow exchangers. Such a RAMEE system has a total effectiveness 6% higher and 1.5% lower compared to the same cross-flow and counter-flow RAMEE, respectively (at NTU=10, Cr*¡Ö3, y0/x0=0.2 and xi/x0=0.1).

It is well established that the temperature gradients in the interiors of internally-heated mantle convection models are subadiabatic (e.g. Parmentier et al., 1994; Bunge et al., 1997, 2001). The subadiabatic gradients have been explained to arise due to a balance between vertical advection and internal heating, however, a detailed analysis of the energy balance in the subadiabatic regions has not been undertaken. In this research, I examine in detail the energy balance in a suite of two-dimensional convection calculations with mixed internal and basal heating, depth-dependent viscosity and continents. I find that there are three causes of subadiabatic gradients. One is the above-mentioned balance, which becomes significant when the ratio of internal heating to surface heat flux is large. The second mechanism involves the growth of the overshoot (maximum and minimum Temperatures along a geotherm) of the geotherm near the lower boundary where the dominant balance is between vertical and horizontal advection. The latter mechanism is significant even in relatively weakly internally heated calculations. For time-dependent calculations, I find that local secular cooling can be a dominant term in the energy equation and can lead to subadiabaticity. However, it does not show its signature on the shape of the time-averaged geotherm. I also compare the basal heat flux with parameterized calculations based on the temperature drop at the core-mantle boundary, calculated both with and without taking the subadiabatic gradient into account and I find a significantly improved fit with its inclusion.<p> I also explore a wide range of parameter space to investigate the dynamical interaction between effects due to surface boundary conditions representing continental and oceanic lithosphere and the endothermic phase boundary at 660 km-depth in two-dimensional Cartesian coordinate convection calculations. I find that phase boundary induced mantle layering is strongly affected by the wavelength of convective flows and mixed surface boundary conditions strongly increase the horizontal wavelength of convection. My study shows that for mixed cases the effects of the surface boundary conditions dominate the effects of the phase boundary. I show that the calculations with complete continental coverage have the most significantly decoupled lower and upper mantle flows and substantial thermal and mechanical layering. Unlike the free-slip case where the surface heat flux decreases substantially with increasing magnitude of the Clapeyron slope, surface heat flux is shown to be almost independent of the Clapeyron slope for mixed boundary condition cases. Although very different when not layered, models with free and mixed surfaces have very similar planforms with very large aspect ratio flows when run with large magnitudes of the Clapeyron slope. I also calculate the critical boundary layer Rayleigh number as a measure of the thermal resistance of the surface boundary layer. My results show that the thermal resistance in the oceanic and the continental regions of the mixed cases are similar to fully free and no-slip cases, respectively. I find that, even for purely basally heated models, the mantle becomes significantly subadiabatic in the presence of partial continental coverage. This is due to the significant horizontal advection of heat that occurs with very large aspect ratio convection cells.

While shoaling from deepwater in a stratified ocean, an interfacial solitary wave (ISW) may experience waveform inversion on a continental margin. Although many oceanographers have believed that the inversion from depression to elevation may commence at the turning point where the upper and bottom layers are equal in depth, this phenomenon has not been fully verified in field observations nor in a laboratory. In this study, a series of laboratory experiments and numerical modeling were conducted on the evolution of an ISW of depression across uniform slope joining a horizontal plateau which resembles pseudo slope-shelf topography, in order to clarify this fascinating phenomenon and the variations of wave properties associated with the process. In the laboratory experiments, a depression ISW was produced by a collapse mechanism in a stratified two-layer fluid system within a steel-framed wave flume (12 m long, 0.7 m high by 0.5 m wide) at the National Sun Yat-sen University in Taiwan. The fluid density in the upper (fresh) and bottom (brine) layers was 996 and 1030 kg/m3, respectively. A series of experiments were conducted upon varying the magnitude of the most important physical factors (i.e., nominal thickness of pycnocline, depth ratio between upper and bottom layer, front gradient and shape of pseudo slope-shelf), from which the results are now discussed in four separate chapters in this thesis. Present laboratory results indicate that the process of waveform inversion took place after an ISW had experienced internal run-down, hydraulic jump, vortex motion and surge-up on the front slope, prior to its propagation onto the plateau. Moreover, the fundamental wave period of leading wave on the plateau was significantly smaller than that in the preceding sections on the front slope and the incident stage earlier, thus representing frequency downshift. Amongst the factors involved, the depth ratio between the upper and bottom layer was the most significant one for waveform inversion. Only when the upper layer was thicker than the bottom layer on the plateau of pseudo slope-shelf, waveform inversion could occur, besides the length of the plateau. On the other hand, the front gradient and shape of pseudo slope-shelf also affected the magnitude of the transmitted wave over the plateau as the wave across this specific topography. In the case of a steeper front gradient, waveform inversion became insignificant due to stronger wave reflection and intense energy dissipation caused by turbulent mixing while a depression ISW propagated over a slope-shelf; particularly against a submerged vertical cliff. As a depression ISW across pseudo slope-shelf with short plateau, intense wave breaking might occur again with vortex motion at its rear end as the newly inversed waveform reentering deep water. In this region, the upper layer was smaller than the bottom layer, hence it could not support the continuous existence of an ISW in elevation. Again, energy dissipation occurred due to turbulent mixing beyond the rear end of a short plateau. Finally, a different mode of ISW appeared within pycnocline, while its nominal thickness was larger than the amplitude of the incident wave. In addition to the laboratory investigations, numerical model was also adopted to study the variations in the flow field as an ISW propagated over a pseudo slope-shelf, in order to complement the experimental results. The results of numerical modeling revealed that the horizontal velocity in the bottom layer increased when the wave encountered the front slope, even if the depth of upper layer was thinner than that of the bottom layer on the plateau. Consequently, the velocity in the upper layer became less than that in the bottom layer when the former was thicker than that of the latter on the plateau. On the other hand, the vertical velocity within the self-generated vortex switched direction as waveform inversion commenced after the wave across the shoulder of pseudo slope-shelf where the local depth of the upper layer was larger than that of bottom part. Overall, the significance of the four pertinent factors (i.e., nominal thickness of pycnocline, water depth ratio, front slope, and plateau length) that affected a depression ISW across pseudo slope-shelf is discussed in detail in this thesis, as well as the variation of flow field calculated by the numerical mode presented.

The sediment resuspension and turbidity created during dredging operations is both an economical and environmental issue. The movement of sediment plumes created from dredging operations has been predicted with numerical modeling, however, these far-field models need a “source term” or near-field model as input. Although data from field tests have been used to create near-field models that predict the amount of material suspended in the water column, these results are skewed due to limitations such as non-uniform sediment distributions, water currents, and water quality issues. Laboratory investigations have obtained data for turbidity during dredging operations, but these results do not take advantage of the most contemporary testing methods. The purpose of this dissertation is to provide an estimation of turbidity created during a cutter suction dredging operation. This estimation was facilited by the development of resuspension measurement and data acquisition techniques in a laboratory setting. Near-field turbidity measurements around the cutter head were measured in the Haynes Coastal Engineering Laboratory at Texas A&M University. The laboratory contains a dredge/tow tank that is ideal for conducting dredging research. A dredge carriage is located in the dredge/tow tank and is composed of a carriage, cradle, and ladder. Acoustic Doppler Velocimetry (ADV) and Optical Backscatter Sensor (OBS) measurements were taken at specific points around the cutter head. The variables of suction flow rate, cutter speed, and the thickness of cut were investigated to understand their specific effect on turbidity generation and turbulence production around the cutter head. A near-field advection diffusion model was created to predict resuspension of sediment from a cutter suction dredge. The model incorporates the laboratory data to determine the velocity field as well as the turbulent diffusion. The model is validated with laboratory testing as well as field data. Conclusions from this research demonstrate undercutting consistently produced larger point specific turbidity maximum than overcutting in the laboratory testing. An increase in suction flow rate was shown to increase production and decrease turbidity around the cutter head. In general, an increase in cutter speed led to an increase in turbidity. The thickness of cut produced less resuspension for a full cut versus a partial cut. Data for a “shallow cut” also produced less turbidity generation than partial cuts. The numerical model was compared to all laboratory testing cases as well as the Calumet Harbor and New Bedford cutter resuspension data and produced suitable MRA values for all tests. The numerical model produced higher point specific regions of turbidity for undercutting but produced larger mean values of turbidity for overcutting.

ABSTRACT Three-dimensional modeling on the aerodynamics of airflow and diffusion of air pollutants in a longitudinal-ventilated traffic tunnel was carried out. The model takes ventilation fans, traffic flow rate, speed, emission factor and piston effect of moving vehicles into consideration. Turbulent flow and dispersion of gaseous pollutants in road tunnels were solved numerically using the finite volume method. Traffic emissions were accordingly modeled as banded line sources along the tunnel floor. The effects of fan ventilation, roughness and piston effect of moving vehicles on the air flow and pollutant dilution are examined. Concentrations of gaseous pollutants CO, NOX, SO2 and THC (total hydrocarbons) at three axial locations in the tunnel, together with traffic flow rate, traffic speed and types of vehicle were measured. Case study was conducted on the Cross-Harbor Tunnel and the Chungcheng Tunnel in which on-site measurements of traffic flow were also conducted concurrently to provide traffic emission data to the tunnel environment for numerical simulation and comparisons. The aim of this study was to understand the spatial variation of air pollutants generated by traffic emissions and evaluation of ventilation performance and piston effect of moving vehicles on dilution of air pollutants in these tunnels. The results show that the major emission sources of CO are passenger cars and motorcycles, while major emission sources of NOx are trucks. Pollutants convect downstream with the wind generated either by longitudinal ventilation fans and/or moving vehicles, thus causing increasing pollutants concentrations with increasing downstream distance. The piston effect of moving vehicle alone can provide 64% ~ 85% increase of wind speed in Chungchen Tunnel and 13% ~ 20% in Cross-Harbor Tunnel. When all fans are on, showing 185% ~ 328% and 120% ~ 182% increases in Chungchen Tunnel and Cross-Harbor Tunnel, respectively. The piston effect of moving vehicle alone can provide 14% ~ 32% dilution of air pollutants in the Chungcheng Tunnel. The piston effect of moving vehicles is compounded with ventilation fans, showing a 47% ~ 66% dilution effect when all fans are on. For the Cross-Harbor Tunnel, the piston effect of moving vehicle alone can provide 9% ~ 23% dilution of air pollutants and 36% ~ 74% dilution effect when all fans are on. The results reveal that cross-sectional concentrations are non-uniformly distributed and that concentrations rise with downstream distance. When all fans were turned off, wind speed in tunnels would be considered as constant, and gaseous pollutants concentration agree with linearly alone the tunnel.

Rod seals are one of the most critical components of hydraulic systems. However, the fundamental physics of seal behavior is still poorly understood and the seal designers have virtually no analytical tools with which to predict the behavior of potential seal designs. In pursuit of a comprehensive physics based seal analysis/ design tool, in this work, a multi-scale multi-physics (MSMP) seal model is developed. The model solves the transient problem involving macro-scale viscoelastic deformation mechanics, macro-scale contact, micro-scale two phase fluid mechanics in the sealing zone, micro-scale asperity contact mechanics and micro-scale deformation mechanics of the sealing edge in a strongly coupled manner. The model takes into account surface roughness, mixed lubrication, cavitation and two phase flow, transient squeeze film effects and the dynamic operation as well as the effect of macro/micro/nano scale viscoelasticity. A hybrid finite element-finite volume-statistical computational framework is developed to solve the highly coupled multi-physics interactions of the MSMP model simultaneously. Surface characterization experiments are performed to extract the parameters like RMS roughness, asperity density, autocorrelation length and asperity radius needed by MSMP. To remove the high frequency noise without removing the high frequency real surface features, a wavelet transform based adaptive surface extraction method is implemented. Dynamic mechanical analysis (DMA) is performed to extract the macro-scale viscoelastic parameters of the seal. Through atomic force microscopy (AFM) experiments, the local micro/nano scale elastic moduli were found to be varying within two orders of magnitude higher than the bulk of the polymer. Significant differences in local stiffness, adhesion and the relaxation time scales of individual surface asperities were also observed. With the MSMP model, dynamic seal performance was analyzed. The results confirmed the mixed lubrication and the effect of surface roughness. Thicker fluid films during instroke and cavitation during the outstroke were found to be important for non-leakage. Seal behavior was a function of the complex dual dependence on the time varying sealed pressure and hydrodynamic effects. Viscoelasticity is seen to critically affect the leakage and friction characteristics. It produces thicker fluid films and produces a significant increase in Poiseuille component of flow during instroke. Ignoring viscoelasticity leads to under-prediction of the time required to reach the zero leakage state. Several high pressure - high frequency sealing applications were analyzed. In such applications, a new phenomenon of "secondary contact" was observed. Viscoelastic creep was seen to critically affect the contact pressure and hence the friction characteristics. In high frequency applications, viscoelasticity induced significant differences in Poiseuille flow and friction force from cycle to cycle. Cycle frequency was seen to play an important role in governing visco-elastohydrodynamics and the leakage of such seals. The seals need to be designed by considering the relationship between relaxation time scales of the polymer and the cycle frequencies. Study also revealed the presence of characteristics like "critical temperature" and "critical frequency". Using the multi-physics modeling capability of MSMP framework, several novel seal designs using smart materials like piezo-ceramic embedded polymers are proposed and analyzed. The MSMP computational framework developed here has a great potential to be used as a stand-alone seal design and analysis software in academic and industrial research.

Impact hammer is the current modal testing way in Dynapac testing department. Due to highly damped characteristic of big construction machines, there are a few weaknesses for modal testing when using hammer, such as short response time, limited frequency resolution, poor quality of frequency response functions. Therefore, a more advanced excitation equipment is needed to improve the measurement quality. The object for this study is to compare two different measuring methods. The thesis will show a comparison between the hammer testing and the shaker MIMO testing compared with analytical model in a highly damped system. It will also give a reference for further highly damped modal analysis and budgetary assessment to decide the budget expenditure. Result from shaker testing shows a little better correlation than hammer testing compared with FEM model. While the correlation between FEM model and measurement is bad due to many reasons, such as many local modes that can not excited, lack of excitation points, unexpected noise and error from the measurement. While considering the compared results obtained from this machine for now, a simpler structure experiment is suggested to be carried on in the future. Shorter length of stinger can be used to enable higher amplitude of force to excite the property on this machine.