A thermal hydraulic analysis has been conducted to investigate the conditions leading to the channel flow reversal and the subsequent effects that may have on the Primary Heat Transfer System(PHTS) thermohydraulic parameters during the natural circulation under the specific accident scenarios for a generic CANDU 900 MW plant model similar to Darlington NGS. The assumed initiating events are the combination of a small Loss Of Coolant Accident (LOCA) with a loss of Class IV power, as well as the unavailability of Emergency Coolant Injection (ECI) system. No makeup inventory is taken into account in this study, and there is no fuel sheath temperature excursion or fuel centerline melting, i.e., the integrity of fuel is always maintained. A one-dimensional quasi-steady state Homogeneous Equilibrium Model(HEM) has been constructed for the study. A specific node-link structure is adopted to represent the primary heat transfer loop: The whole loop and different components in HTS are represented by a series of nodes that have quasi-static thermal hydraulic characteristics such as pressure and enthalpy, etc. Dynamic characteristics are delivered by the links between nodes, e.g., flow rate and pressure drop. The channel powers (decay heat), the secondary side pressure and the pressure at Reactor Inlet Header(RIH) are chosen as boundary conditions to describe the assumed initiating incidents for the iii model. With ongoing loss of inventory and system depressurization, vapor lock occurs in Steam Generator(SG), and it increases the pressure drop from Reactor Outlet Header(ROH) to RIH across SG and forms an increasingly negative RIH-to-ROH pressure differential. Flow reversal occurs in the channel due to the counter force balance between negative RIH-to-ROH pressure difference and the driving force derived from the density difference between the hot and cold legs. It is found that channels in row A have the highest reversal preference, then followed by the channels in row B, and in that order subsequently. Row A reverses when inventory decreases to about 79.5% of initial value, with following boundary conditions: decay heat is 1.5% of Full Power(FP) and secondary side pressure is 5.070 MPa. In addition, it is found that the decrease in channel power accelerates the depressurization process and brings forward the occurrence of flow reversal in fuel channels. / Thesis / Master of Applied Science (MASc)

This study was conducted to ascertain the impacts of bed leveling, following ship channel dredging operations, and to also investigate the hydrodynamic flow field around box bed levelers. Laboratory experiments were conducted with bed levelers operating in the laboratory using video cameras for flow visualization. Computer software and numerical codes, called FANS, were used to validate the laboratory experiments. The study was split into two major parts: laboratory experiments and hydrodynamic modeling. The laboratory experiment was conducted to model how bed levelers interact with the ship channel bottom after hopper dredge dragheads (blades) made passes and created uneven trenches. These interactions were observed using both underwater and hand-held cameras. The hydrodynamic modeling was accomplished using GRIDGEN and PEGSUS commercial software for generating grid and input data files in the pre-processing phase, Finite-Analytic Navier-Stokes (FANS) software for simulation in the processing phase, and two commercial software (Fieldview and Tecplot) for plotting the images and graphs in the post-processing phase. An interesting phenomenon was observed in the laboratory experimental runs. The flow field showed reversed flow in front of the moving bed leveler and the trench parallel to the direction of the bed leveler. The flow in the parallel trench was observed to be in the same direction as the bed leveler movement, and it was expected that the flow would travel under the bed leveler. The bed leveler was towed at two specified constant speeds: 0.25 m/s (0.82 ft/s) and 0.5 m/s (1.64 ft/s) and at a water depth of 1.22 m (4.00 ft) Similarly, the images and plots of the hydrodynamic modeling obtained from FieldView and Tecplot software showed flow reversal, depicted by the negative velocities, within the vicinity of the trench, as the model bed leveler moved past and interacted with the fluid. The negative velocity had a magnitude close to 0.5 m/s (1.64 ft/s), which was the velocity used in running the laboratory experiments. The hydrodynamic simulation matched closely with the experimental observations, and thus, the laboratory observation was confirmed. The final results obtained from the numerical modeling helped to understand the hydrodynamic effects around the box bed leveler.

A pulsating jet generated by a pulse combustor has been experimentally demonstrated as a technique for impingement heat transfer enhancement relative to a steady jet. The enhancement factor was as high as 2.5. Despite such potential, further studies of this technique have been limited, let alone industrial applications. The ultimate goal of the Pulsed Air Drying project at the Institute of Paper Science and Technology is to develop this technique to commercialization for industrial applications such as paper drying. The main objective of the research in this dissertation is to provide a fundamental basis for the development of the technology. Using CFD simulations, the research studied the characteristics of pulsating single-slot-nozzle jet impingement flows and heat transfer on stationary and moving surfaces. In addition, in order to understand basic flow characteristics of pulse-combustor jets, a simplified model of Helmholtz pulse combustors was developed. The model was used to recommend a strategy to generate a pulsating jet having large amplitude of velocity oscillation. And based on this model, pulsating jets in the simulations were characterized as those at the tailpipe exit of a pulse combustor. The impingement conditions were similar to those in conventional impingement hoods for paper drying. Parameter studies included the effects of jet velocity oscillation amplitude, pulsation frequency, mean jet velocity, tailpipe width, and impingement surface velocity. Simulation results showed that the amplitude of jet velocity oscillation was the most important parameter for heat transfer enhancement, in which two mechanisms were identified: high impinging jet velocity during the positive cycle and strong re-circulating flows in the impingement zone during the negative cycle of jet velocity oscillation. As for the improvement by the pulsating jets relative to steady jets, the maximum heat transfer enhancement and energy saving factors were 1.8 and 3.0, respectively, which were very encouraging for further development of the technology.

O acidente vascular cerebral (AVC) e a primeira causa de obito no Brasil e a principal de incapacidade no mundo. A doenca aterosclerotica carotidea esta associada a cerca de 15% dos AVC. A endarterectomia cirurgica e o tratamento padrao-­ouro para as estenoses carotideas sintomaticas. A angioplastia carotidea com stent via endovascular (ACS) e uma alternativa terapeutica a endarterectomia. O desenvolvimento de dispositivos de protecao embolica levou a reducao das taxas de complicacao durante a ACS. A tecnica de protecao mais utilizada na pratica clinica e a protecao distal com filtro. As estrategias de protecao proximal, por bloqueio ou reversao do fluxo carotideo, tem revelado resultados promissores em relacao a protecao distal com filtro durante a ACS. Assim, o objetivo do estudo foi comparar a reversao de fluxo e o filtro de protecao durante a ACS. Os objetivos primarios foram a incidencia, o numero e o tamanho das lesoes cerebrais isquemicas por ressonancia magnetica pela sequencia diffusion-­weighted-­imaging (DWI) apos a ACS. Os objetivos secundarios foram as complicacoes cardiacas e eventos cerebrovasculares (CCEC), ataque isquemico transitorio (AIT) e as lesoes cerebrais isquemicas definitivas na imagem por ressonancia magnetica pela sequencia fluid-­attenuated inversion recovery (FLAIR) no seguimento de tres meses. As lesoes isquemicas foram obtidas por ressonancia magnetica (RM) de campo 3 Tesla. Os pacientes foram avaliados quanto aos desfechos neurologicos atraves das escalas do National Institutes of Health Stroke Scale (NIHSS) e modified Rankin Scale (mRS). Quarenta pacientes foram incluidos randomicamente no estudo. Comparado a reversao de fluxo (n=21) o filtro de protecao (n=19) resultou em uma reducao significativa na incidencia (15,8% vs. 47,6%, p=0,03), no numero (0,73 vs. 2,6, p=0,05) e no tamanho (0,81 vs. 2,23 mm, p=0,05) das novas lesoes cerebrais isquemicas. Dois pacientes, um de cada grupo, apresentaram AIT no seguimento de tres meses. Nao foram observadas CCEC na internacao ou no seguimento de tres meses. Neste estudo, o filtro de protecao mais eficaz que a reversao de fluxo na protecao cerebral durante a angioplastia carotidea com stent por acesso femoral. / Objectives: Carotid artery stenting (CAS) has become an alternative treatment for patients presenting symptomatic carotid artery stenosis. The improvement in clinical outcomes with CAS has been associated with the development of embolic protection devices. The trial aim is to compare flow-­reversal versus filter protection during CAS through femoral access. Methods: Patients were randomly enrolled in CAS using flow-­reversal or filter protection. The primary endpoints were the incidence, number and size of new ischemic brain lesions after CAS. The secondary endpoints included major adverse cardiac and cerebrovascular events (MACCE), transient ischemic attack (TIA) and definitive ischemic brain lesions on fluid-­attenuated inversion recovery magnetic resonance image (FLAIR-­MRI) at a three-­month follow-­up. Ischemic brain lesions were assessed by a 3T-­MRI. Neurological outcomes were evaluated by means of the National Institutes of Health Stroke Scale (NIHSS) and the modified Rankin Scale (mRS). Results: Forty consecutive patients were randomly assigned. Compared to flow-­reversal (n=21), filter protection (n=19) resulted in a significant reduction in the incidence (15.8% vs. 47.6%, p=0.03), number (0.73 vs. 2.6, p=0.05) and size (0.81 vs. 2.23 mm, p=0.05) of new ischemic brain lesions. Two patients, one from each group, presented TIA at the three-­month follow-­up. There were no MACCE in the hospital or at the three-­month follow-­up. Conclusions: In this small sample size trial, filter protection was more effective than flow-­reversal in reducing ischemic brain lesions during CAS through femoral approach.

O acidente vascular cerebral (AVC) e a primeira causa de obito no Brasil e a principal de incapacidade no mundo. A doenca aterosclerotica carotidea esta associada a cerca de 15% dos AVC. A endarterectomia cirurgica e o tratamento padrao-­ouro para as estenoses carotideas sintomaticas. A angioplastia carotidea com stent via endovascular (ACS) e uma alternativa terapeutica a endarterectomia. O desenvolvimento de dispositivos de protecao embolica levou a reducao das taxas de complicacao durante a ACS. A tecnica de protecao mais utilizada na pratica clinica e a protecao distal com filtro. As estrategias de protecao proximal, por bloqueio ou reversao do fluxo carotideo, tem revelado resultados promissores em relacao a protecao distal com filtro durante a ACS. Assim, o objetivo do estudo foi comparar a reversao de fluxo e o filtro de protecao durante a ACS. Os objetivos primarios foram a incidencia, o numero e o tamanho das lesoes cerebrais isquemicas por ressonancia magnetica pela sequencia diffusion-­weighted-­imaging (DWI) apos a ACS. Os objetivos secundarios foram as complicacoes cardiacas e eventos cerebrovasculares (CCEC), ataque isquemico transitorio (AIT) e as lesoes cerebrais isquemicas definitivas na imagem por ressonancia magnetica pela sequencia fluid-­attenuated inversion recovery (FLAIR) no seguimento de tres meses. As lesoes isquemicas foram obtidas por ressonancia magnetica (RM) de campo 3 Tesla. Os pacientes foram avaliados quanto aos desfechos neurologicos atraves das escalas do National Institutes of Health Stroke Scale (NIHSS) e modified Rankin Scale (mRS). Quarenta pacientes foram incluidos randomicamente no estudo. Comparado a reversao de fluxo (n=21) o filtro de protecao (n=19) resultou em uma reducao significativa na incidencia (15,8% vs. 47,6%, p=0,03), no numero (0,73 vs. 2,6, p=0,05) e no tamanho (0,81 vs. 2,23 mm, p=0,05) das novas lesoes cerebrais isquemicas. Dois pacientes, um de cada grupo, apresentaram AIT no seguimento de tres meses. Nao foram observadas CCEC na internacao ou no seguimento de tres meses. Neste estudo, o filtro de protecao mais eficaz que a reversao de fluxo na protecao cerebral durante a angioplastia carotidea com stent por acesso femoral. / Objectives: Carotid artery stenting (CAS) has become an alternative treatment for patients presenting symptomatic carotid artery stenosis. The improvement in clinical outcomes with CAS has been associated with the development of embolic protection devices. The trial aim is to compare flow-­reversal versus filter protection during CAS through femoral access. Methods: Patients were randomly enrolled in CAS using flow-­reversal or filter protection. The primary endpoints were the incidence, number and size of new ischemic brain lesions after CAS. The secondary endpoints included major adverse cardiac and cerebrovascular events (MACCE), transient ischemic attack (TIA) and definitive ischemic brain lesions on fluid-­attenuated inversion recovery magnetic resonance image (FLAIR-­MRI) at a three-­month follow-­up. Ischemic brain lesions were assessed by a 3T-­MRI. Neurological outcomes were evaluated by means of the National Institutes of Health Stroke Scale (NIHSS) and the modified Rankin Scale (mRS). Results: Forty consecutive patients were randomly assigned. Compared to flow-­reversal (n=21), filter protection (n=19) resulted in a significant reduction in the incidence (15.8% vs. 47.6%, p=0.03), number (0.73 vs. 2.6, p=0.05) and size (0.81 vs. 2.23 mm, p=0.05) of new ischemic brain lesions. Two patients, one from each group, presented TIA at the three-­month follow-­up. There were no MACCE in the hospital or at the three-­month follow-­up. Conclusions: In this small sample size trial, filter protection was more effective than flow-­reversal in reducing ischemic brain lesions during CAS through femoral approach.

District heating systems are common in Nordic countries today and accounts for a great portion of the heat demand. In Sweden, total district heating end use in the last years has been around 50 TWh and district heating accounts for roughly 50 % of the total heat demand. Suppliers of district heating must balance demand and supply, often in large and complex networks. Heat propagation can be in the range of hours and it is not known in detail how the heat will propagate during transient conditions. A dynamic model has been developed in OpenModelica and a method for modeling, handling data, simulating and visualizing the results of a district heating network was developed using Matlab as core. Data from Mälarenergi AB, a district heating producer and grid operator, was used for validation of the model. Validation shows that the model works well in predicting heat propagation and temperature distribution in the network and that the model can be scaled up to a large number of heat exchangers and pipes. The model is robust and can handle bi-directional and reversing flows in complex ring structures. It was concluded that OpenModelica together with Matlab is a good combination for creating models of district heating networks, as a high degree of standardization and automation can be achieved. This, together with visualization of the heat propagation, makes it useful for the understanding of the district heating network during transient conditions. / Smarta Flöden

The flow-reversal elliptical cylindrical end chamber mufflers of short length are used often in the modern day automotive exhaust systems. The conventional 1-D axial plane wave theory is not able to predict their acoustical attenuation performance in view of the fact that the chamber length is not enough for the evanescent 3-D modes generated at the junctions to decay sufficiently for frequencies below the cut-off frequency. Also, due to the large area expansion ratio at the inlet, the first few higher order modes get cut on even in the low frequency regime. This necessitates a 3-D FEM or 3-D BEM analysis, which is cumbersome and time consuming. Therefore, an ingenious 1-D transverse plane wave theory is developed by considering plane wave propagation along the major-axis of the elliptical section, whereby a 2-port axially short elliptical and circular chamber muffler is characterized by means of the transfer matrix [T] or impedance matrix [Z]. Two different approaches are followed: (1) a numerical scheme such as the Matrizant approach, and (2) an analytical approach based upon the Frobenius series solution of the Webster’s equation governing the transverse plane wave propagation. The convective effects of mean flow are neglected; however the dissipative effects at the ports are taken into account. The TL predicted by this 1-D transverse plane wave analysis is compared with that obtained by means of the 3-D analytical approach and numerical (FEM/BEM) methods. An excellent agreement is observed between this simplified 1-D approach and the 3-D approaches at least up to the cut-on frequency of the (1, 1) even mode in the case of elliptical cylindrical chambers, or the (1, 0) mode in the case of circular cylindrical chambers, thereby validating this 1-D transverse plane wave theory. The acoustical attenuation characteristics of such short chamber mufflers for various configurations are discussed, qualitatively as well as quantitatively. Moreover, the Frobenius series solution enables one to obtain non-dimensional frequencies for determining the resonance peak and trough in the TL graph. The use of this theory is, however, limited to configurations in which both the ports are located along the major axis in the case of elliptical chambers and along the same diameter for circular chambers. The method of cascading the [T] matrices of the 2-port elements cannot be used to analyze a network arrangement of 2-port elements owing to the non-unique direction of wave propagation in such a network of acoustic elements. Although, a few papers are found in the literature reporting the analysis of a network of 2-port acoustic elements, no work is seen on the analysis of a network of multi-port elements having more than two external ports. Therefore, a generalized algorithm is proposed for analyzing a general network arrangement of linear multi-port acoustic elements having N inlet ports and M outlet ports. Each of these multi-port elements constituting the network may be interconnected to each other in an arbitrary manner. By appropriate book-keeping of the equations obtained by the [Z] matrix characterizing each of the multi-port and 2-port elements along with the junction laws (which imply the equality of acoustic pressure and conservativeness of mass velocity at a multi-port junction), an overall connectivity matrix is obtained, whereupon a global [Z] matrix is obtained which characterizes the entire network. Generalized expressions are derived for the evaluation of acoustic performance evaluation parameters such as transmission loss (TL) and insertion loss (IL) for a multiple inlet and multiple outlet (MIMO) system. Some of the characteristic properties of a general multi-port element are also studied in this chapter. The 1-D axial and transverse plane wave analysis is used to characterize axially long and short chambers, respectively, in terms of the [Z] matrix. Different network arrangements of multi-port elements are constructed, wherein the TL performance of such MIMO networks obtained on the basis of either the 1-D axial or 1-D transverse plane wave theory are compared with 3-D FEA carried on a commercial software. The versatility of this algorithm is that it can deal with more than two external or terminal ports, i.e., one can have multiple inlets and outlets in a complicated acoustic network. A generalized approach/algorithm is presented to characterize rigid wall reactive multi-port chamber mufflers of different geometries by means of a 3-D analytical formulation based upon the modal expansion and the uniform piston-driven model. The geometries analyzed here are rectangular plenum chambers, circular cylindrical chamber mufflers with and without a pass tube, elliptical cylindrical chamber mufflers, spherical and hemispherical chambers, conical chamber mufflers with and without a co-axial pass tube and sectoral cylindrical chamber mufflers of circular and elliptical cross-section as well as sectoral conical chamber mufflers. Computer codes or subroutines have been developed wherein by choosing appropriate mode functions in the generalized pressure response function, one can characterize a multi-port chamber muffler of any of the aforementioned separable geometrical shapes in terms of the [Z] matrix, subsequent to which the TL performance of these chambers is evaluated in terms of the scattering matrix [S] parameters by making use of the relations between [Z] and [S] matrices derived earlier. Interestingly, the [Z] matrix approach combined with the uniform piston-driven model is indeed ideally suited for the 3-D analytical formulation inasmuch as regardless of the number of ports, one deals with only one area discontinuity at a time, thereby making the analysis convenient for a multi-port muffler configuration with arbitrary location of ports. The TL characteristics of SISO chambers corresponding to each of the aforementioned geometries (especially the elliptical cylindrical chamber) are analyzed in detail with respect to the effect of chamber dimensions (chamber length and transverse dimensions), and relative angular and axial location of ports. Furthermore, the analysis of SIDO (i.e., single inlet and double outlet) chamber mufflers is given special consideration. In particular, we examine (1) the effect of additional outlet port (second outlet port), (2) variation in the relative angular or axial location of the additional or second outlet port (keeping the location of the inlet port and the outlet ports of the original SISO chamber to be constant) and (3) the effect of interchanging the location of the inlet and outlet ports on the TL performance of these mufflers. Thus, design guidelines are developed for the optimal location of the inlet and outlet ports keeping in mind the broadband attenuation characteristics for a single inlet and multiple outlet (SIMO) system. The non-dimensional limits up to which a flow-reversal elliptical (or circular) cylindrical end chamber having an end-inlet and end-outlet configuration is acoustically short (so that the 1-D transverse plane wave theory is applicable) and the limits beyond which it is acoustically long (so that the 1-D axial plane wave theory is applicable) is determined in terms of the ratio or equivalently, in terms of the ratio. Towards this end, two different configurations of the elliptical cylindrical chamber are considered, namely, (1) End-Offset Inlet (located along the major-axis of the ellipse) and End-Centered Outlet (2) End-Offset Inlet and End-Offset Outlet (both the ports located on the major-axis of the ellipse and at equal offset distance from the center). The former configuration is analyzed using 3-D FEA simulations (on SYSNOISE) while the 3-D analytical uniform piston-driven model is used to analyze the latter configuration. The existence of the higher order evanescent modes in the axially long reversal chamber at low frequency (before the cut-on frequency of the (1, 1) even mode or (1, 0) mode) causes a shift in the resonance peak predicted by the 1-D axial plane wave theory and 3-D analytical approach. Thus, the 1-D axial plane wave analysis is corrected by introducing appropriate end correction due to the modified or effective length of the elliptical cylindrical chamber. An empirical formulae has been developed to obtain the average non-dimensional end correction for the aforementioned configurations as functions of the expansion ratio, (i.e., ), minor-axis to major-axis ratio, (i.e., ) and the center-offset distance ratio, (i.e., ). The intermediate limits between which the chamber is neither short nor long (acoustically) has also been obtained. Furthermore, an ingenious method (Quasi 1-D approach) of combining the 1-D transverse plane wave model with the 1-D axial plane wave model using the [Z] matrix is also proposed for the end-offset inlet and end-centered outlet configuration. A 3-D analytical procedure has also been developed which also enables one to determine the end-correction in axially long 2-port flow-reversal end chamber mufflers for different geometries such as rectangular, circular and elliptical cylindrical as well as conical chambers, a priori to the computation of TL. Using this novel analytical technique, we determine the end correction for arbitrary locations on the two end ports on the end face of an axially long flow-reversal end chamber. The applicability of this method is also demonstrated for determination of the end corrections for the 2-port circular cylindrical chamber configuration without and with a pass tube, elliptical cylindrical chambers as well as rectangular and conical chambers.