This thesis is concerned with the study of parallel mixing of two dissimilar gases under compressible conditions in the confined environment. A number of numerical studies are reported in the literature for the compressible mixing of two streams of gases where (1) both the streams are of similar gases at the same temperatures, (2) both the streams are at different temperatures with similar gases, and (3) dissimilar gases are with nearly equal temperatures. The combination of dissimilar gases at large temperature difference, mixing under compressible conditions, as in the case of scramjet propulsion, has not been adequately addressed numerically. Also many of the earlier studies have used two dimensional numerical simulation and showed good match with the experimental results on mixing layers that are inherently three dimensional in nature. In the present study, both two-dimensional (2-d) and three dimensional (3-d) studies are reported and in particular the effect of side wall on the three dimensionality of the flow field is analyzed, and the reasons of the good match of two dimensional simulations with experimental results have been discussed. Both two dimensional and three dimensional model free simulations have been conducted for a flow configuration on which experimental results are available. In this flow configuration, the mixing duct has a rectangular cross section with height to width ratio of 0.5. In the upper part of the duct hydrogen gas at a temperature of 103 K is injected through a single manifold of two Ludweig tubes and in the lower part of the duct nitrogen gas at a temperature of 2436 K is supplied through an expansion tube, both the gases are at Mach numbers of 3.1 and 4.0 respectively. Measurements in the experiment are limited to wall pressures and heat flux. The choice of this experimental condition gives an opportunity to study the effect of large temperature difference on the mixing of two dissimilar gases with large molecular weights under compressible conditions. Both two dimensional and three dimensional model free simulations are carried out using higher order numerical scheme (4th order spatial and 2nd order temporal) to understand the structure and evolution of supersonic confined mixing layer of similar and dissimilar gases. Two dimensional simulations are carried out by both SPARK (finite difference method) and OpenFOAM (finite volume method based open source software that was specially picked out and put together), while 3D model free simulations are carried out by OpenFOAM. A fine grid structure with higher grid resolution near the walls and shear layer is chosen. The effect of forcing of fluctuations on the inlet velocity shows no appreciable change in the fully developed turbulent region of the flow. The flow variables are averaged after the attainment of statistical steady state established through monitoring the concentration of inert species introduced in the initial guess. The effect of side wall on the flow structure on the mixing layer is studied by comparing the simulation results with and without side wall. Two dimensional simulations show a good match for the growth rate of shear layer and experimental wall pressures. Three dimensional simulations without side wall shows 14% higher growth rate of shear layer than that of two dimensional simulations. The wall pressures predicted by these three dimensional simulations are also lower than that predicted using two dimensional simulations (6%) and experimental (9%) results in the downstream direction of the mixing duct. Three dimensionality of the flow is thought of as a cause for these differences. Simulations with the presence of side wall show that there is no remarkable difference of three dimensionality of the flow in terms of the variables and turbulence statistics compared to the case without side walls. However, the growth rate of shear layer and wall surface pressures matches well with that predicted using two dimensional simulations. It has been argued that this good match in shear layer growth rate occurs due to formation of oblique disturbances in presence of side walls that are considered responsible for the decrease in growth rate in 3-d mixing layers. The wall pressure match is argued to be good because of hindrance from side wall in the distribution of momentum in third direction results in higher wall pressure. The effect of dissimilar gases at large temperature difference on the growth rate reduction in compressible conditions is studied. Taking experimental conditions as baseline case, simulations are carried out for a range of convective Mach numbers. Simulations are also carried out for the same range of convective Mach numbers considering the mixing of similar gases at the same temperature. The normalized growth rates with incompressible counterpart for both the cases show that the dissimilar gas combination with large temperature difference shows higher growth rate. This result confirms earlier stability analysis that predicts increased growth rate for such cases. The growth rate reduction of a compressible mixing layer is argued to occur due to reduced pressure strain term in the Reynolds stress equation. This reduction also requires the pressure and density fluctuation correlation to be very near to unity. This holds good for a mixing layer formed between two similar gases at same temperature. For dissimilar gases at different temperatures this assumption does not hold well, and pressure-density correlation coefficient shows departure from unity. Further analysis of temperature density correlation factor, and temperature fluctuations shows that the changes in density occur predominantly due to temperature effects, than due to pressure effects. The mechanism of density variations is found to be different for similar and dissimilar gases, while for similar gases the density variations are due to pressure variations. For dissimilar gases density variation is also affected by temperature variations in addition to pressure variations. It has been observed that the traditional k-ε turbulence model within the RANS (Reynolds Averaged Navier Stokes) framework fails to capture the growth rate reduction for compressible shear layers. The performance of k-ε turbulence model is tested for the mixing of dissimilar gases at large temperature difference. For the experimental test case the shear layer growth rate and wall pressures show good match with other model free simulations. Simulations are further carried out for a range of convective Mach numbers keeping the mixing gases and their temperatures same. It has been observed that a drop in the growth rate is well predicted by RANS simulations. Further, the compressibility option has been removed and it has been observed that for the density and temperature difference, even for incompressible case, the drop in growth rate exists. This behaviour shows that the decrease in growth rate is mainly due to the interaction of temperature and species mass fraction on density. Also it can be inferred that RANS with k-ε turbulence model is able to capture the compressible shear layer growth rate for dissimilar gases at high temperature difference. The mixing of heat and species is governed by the values of turbulent Prandtl and Schmidt numbers respectively. These numbers have been observed to vary for different flow conditions, while affecting the flow field considerable in the form of temperature and species distribution. Model free simulations are carried out on an incompressible convective Mach number mixing layer, and the results are compared with that of a compressible mixing layer to study the effect of compressibility on the values of turbulent Prandtl / Schmidt numbers. It has been observed that both turbulent Prandtl and Schmidt numbers show an almost constant value in the mixing layer region for incompressible case. While, for a compressible case, both turbulent Prandtl and Schmidt numbers show a continuous variation within the mixing layer. However, the turbulent Lewis number is observed to be near unity for both incompressible and compressible cases. The thesis is composed of 8 chapters. An introduction of the subject with critical and relevant literature survey is presented in chapter 1. Chapter 2 describes the mathematical formulation and assumptions along with solution methodology needed for the simulations. Chapter 3 deals with the two and three dimensional model free simulations of the non reacting mixing layer. The effect of the presence of side wall is studied in chapter 4. Chapter 5 deals with the effect of compressibility on the mixing of two dissimilar gases at largely different temperatures. The performance of k-ε turbulence model is checked for dissimilar gases in Chapter 6. Chapter 7 is concerned with the effect of compressibility on turbulent Prandtl and Schmidt numbers. Finally concluding remarks are presented in chapter 8. The main aim of this thesis is the exploration of parallel mixing of dissimilar gases under compressible conditions for both two and three dimensional cases. The outcome of the thesis is (a) a finding that the presence of sidewall in a mixing duct does not make flow field two dimensional, instead it causes the formation of oblique disturbances and the shear layer growth rate is reduced, (b) that it has been shown that the growth rates of dissimilar gases are affected far more by large temperature difference than by compressibility as in case of similar gases, (c) that the growth rates of compressible shear layers formed between dissimilar gases are better predicted using k-εturbulence model and (d) that for compressible mixing conditions the turbulent Prandtl and Schmidt numbers vary continuously in the mixing layer region necessitating the use of some kind of model instead of assuming constant values.

The new legal regulations derived from climate change dictate that hydraulic structures must be designed to handle flood events associated with return periods up to 10,000 years. This obviously involves adapting the existing infrastructure to meet such requirements. In order to avoid risks in the restitution of the flow discharged to rivers, such as bank overflows or streambed erosion and scour processes, hydraulic design must be supported by reliable tools capable of reproducing the behavior of hydraulic structures. In the work presented herein, a fully three-dimensional CFD model to reproduce the behavior of different types of air-water flow in hydraulic structures is presented. The flow is assumed to be turbulent, isotropic and incompressible. Several RANS turbulence models are tested and structured rectangular meshes are employed to discretize the analyzed domain. The presence of two fluids is modeled using different VOF approaches and simulations are run using the PIMPLE algorithm. The model is implemented using the open-source platform OpenFOAM and its performance is compared to the commercial code FLOW-3D. The analysis is conducted separately on two different parts of hydraulic structures, namely: the spillway and the stilling basin. Additionally, a case of practical application, where the model reproduces the flow of a real-life case, is also presented in order to prove the suitability of the model to actual design cases. Mesh independence and model validation using experimental data are checked in the results of all the case studies. The sensitivity of the presented model to certain parameters is extensively discussed using different indicator variables. Among these parameters are turbulence closure, discretization scheme, surface tracking approach, CFD code or boundary conditions. Pros and contras of each of them are addressed. The analyzed turbulence models are the Standard k ¿ ¿, the Realizable k ¿ ¿, the RNG k ¿ ¿, and the SST k ¿ ¿. The discretization schemes under study are: a first-order upwind method, the second-order limited Van Leer method, and a second-order limited central difference method. The VOF approaches analyzed are the Partial VOF, as implemented in OpenFOAM, and the TruVOF, as implemented in FLOW-3D. In most cases, the Standard k ¿ ¿ model provides the most accurate estimations of water free surface profiles, although the rest of variables, with few exceptions, are better predicted by the RNG k ¿ ¿. The latter model generally requires slightly longer computation times. The SST k ¿ ¿ reproduces correctly the phenomena under study, although it generally turned out to be less accurate than its k ¿ ¿ counterparts. As regards the comparison among VOF approaches and codes, it is impossible to determine which one performs best. E.g. OpenFOAM, using the Partial VOF, managed to reproduce the in- ternal hydraulic jump structure and all derived variables better than FLOW-3D, using the TruVOF, although the latter seems to capture better the momentum transfer and so all derived variables. In the case of flow in stepped spillways, OpenFOAM captures better the velocity profiles, although FLOW-3D is more accurate when estimating the water free surface profile. It is worth remark- ing that not even their response to certain model parameters is comparable. E.g. FLOW-3D is significantly less sensitive to mesh refinement than OpenFOAM. Given the result accuracy achieved in all cases, the proposed model is fully applicable to more complex design cases, where stilling basins, stepped spillways and hydraulic structures in general must be investigated. / Las nuevas disposiciones legales derivadas del cambio climático dictaminan que las estructuras hidráulicas sean capaces de funcionar correctamente con eventos de inundación asociados a periodos de retorno de hasta 10,000 años. Esto, obviamente, implica adaptar la infraestructura existente para satisfacer dichos requerimientos. A fin de evitar riesgos en la restitución de los caudales vertidos al río, como desbordamientos o procesos erosivos y de socavación, el diseño hidráulico ha de sustentarse en herramientas fiables capaces de reproducir el comportamiento de las estructuras hidráulicas. En este trabajo, se presenta un modelo numérico CFD completamente tridimensional para reproducir el comportamiento de diferentes tipos de flujo aire-agua en estructuras hidráulicas. Se asume que el flujo es turbulento, isotrópico e incompresible. Diversos modelos de turbulencia RANS son contrastados y se emplean mallas estructuradas rectanuglares para discretizar el dominio analizado. La presencia de dos fluidos es modelada utilizando diferentes enfoques VOF y las simulaciones son ejecutadas empleando el algoritmo PIMPLE. El modelo es implementado mediante la plataforma de código abierto OpenFOAM y su respuesta es comparada con la del modelo comercial FLOW-3D. El análisis se lleva a cabo sobre dos partes diferentes de una estructura hidráulica, a saber, el aliviadero y el cuenco amortiguador, de forma separada. Además, un caso de aplicación práctica, donde el modelo reproduce el flujo en una estructura real, es presentado también a fin de probar la adecuación del modelo a casos de diseño aplicado. Se comprueban la independencia de la malla y la validación con datos experimentales de los resultados de todos los casos de estudio. La sensibilidad del modelo presentado a ciertos parámetros es analizada de forma exhaustiva empleando diferentes variables indicadoras. Los pros y contras de cada uno de éstos son planteados. Los modelos de turbulencia analizados son el Standard k-epsilon, el Realizable k-epsilon, el RNG k-epsilon y el SST k-omega. Los esquemas de discretización estudiados son: un método de primer orden upwind, uno de Van Leer de segundo orden y un esquema de segundo orden limitado de diferencias centradas. Los enfoques VOF analizados son el Partial VOF, implementado en OpenFOAM, y el TruVOF, implementado en FLOW-3D. En la mayoría de casos, el modelo k-epsilon aporta las estimaciones más precisas de perfiles de lámina libre de agua, pese a que el resto de variables, con alguna excepción, son mejor predichas por el RNG k-epsilon. Este modelo generalmente requiere mayores tiempos de cálculo. El k-omega reproduce correctamente los fenómenos bajo estudio, pese a que su precisión es generalmente más baja que la de los modelos k-epsilon. En lo que respecta a la comparación entre enfoques VOF y códigos, es imposible determinar cuál es el mejor. Por ejemplo, OpenFOAM, empleando el Partial VOF, logra reproducir la estructura interna del resalto hidráulico y todas las variables derivadas mejor que FLOW-3D, empleando el TruVOF, a pesar de que este último parece capturar mejor la transferencia de cantidad de movimiento y, por tanto, todas las variables derivadas. En el caso del flujo en aliviaderos escalonados, OpenFOAM captura mejor los perfiles de velocidad, pese a que FLOW-3D es más preciso en la estimación de los perfiles de lámina libre de agua. Conviene recalcar que ni tan sólo su respuesta a ciertos parámetros del modelo es comparable. Por ejemplo, FLOW-3D es significativamente menos sensible al refinado de malla que OpenFOAM. A la luz de la precisión de los resultados obtenidos en todos los casos, el modelo propuesto es completamente aplicable a casos de diseño más complejos, donde cuencos amortiguadores, aliviaderos escalonados y estructuras hidráulicas en general han de ser investigadas. / Les noves disposicions legals derivades del canvi climàtic dictaminen que cal que les estructures hidràuliques siguen capaces de funcionar correctament amb esdeveniments d'inundació associats a períodes de retorn de fins a 10,000 anys. Això, òbviament, implica adaptar la infraestrctura existent per satisfer aquests requeriments. A fi d'evitar riscs en la restitució dels cabals vessats al riu, com desbordaments o processos erosius i de socavació, el disseny hidràulic ha de recolzar-se en ferramentes fiables capaces de reproduir el comportament de les estructures hidràuliques. En aquest treball, es prsenta un model numèric CFD completament tridimensional per a reproduir el comportament de diferents tipus de flux aire-aigua en estructures hidràuliques. S'assumeix que el flux és turbulent, isotròpic i incompressible. Diferents models de turbulència RANS són contrastats i s'empren malles estructurades rectangulars per discretitzar el domini analitzat. La presència de dos fluids és modelada utilitzant diferents enfocaments VOF i les simulacions són executades emprant l'algorisme PIMPLE. El model és implementat mitjançant la plataforma de codi obert OpenFOAM i la seua resposta és comparada amb la del codi comercial FLOW-3D. L'anàlisi es du a terme sobre les diferents parts d'una estructura hidràulica, a saber, sobreeixidors esgraonats i vas esmorteïdor, de forma separada. A més, un cas d'aplicació pràctica, on el model reprodueix el flux a una estructura real, és presentat també a fi de provar l'adequació del model a casos de disseny aplicat. Es comproven la independència de la malla i la validació amb dades experimentals dels resultats de tots els casos d'estudi. La sensibilitat del model presentat a certs paràmetres és analitzada de forma exhaustiva emprant diferents variables indicadores. Els pros i contres de cadascun d'aquests són plantejats. Els models de turbulència analitzats són l'Standard k-epsilon, el Realizable k-epsilon, el RNG k-epsilon i l'SST k-omega. Els esquemes de discretització estudiats són: un mètode de primer ordre upwind, un de Van Leer de segon ordre i un esquema de segon ordre limitat de diferències centrades. Els enfocaments VOF analitzats són el Partial VOF, implementat en OpenFOAM, i el TruVOF, implementat en FLOW-3D. En la majoria de casos, el model Standard k-epsilon aporta les estimacions més precises de perfils de làmina lliure d'aigua, tot i que la resta de variables, amb alguna excepció, són millor predites pel RNG k-epsilon. Aquest model generalment requereix majors temps de càlcul. El k-omega reprodueix correctament els fenòmens sota estudi, tot i que la seua precisió és generalment més baixa que la dels models k-epsilon. Pel que fa la comparació entre enfocaments VOF i codis, és impossible determinar quin és el millor. Per exemple, OpenFOAM, emprant el Partial VOF, aconsegueix reproduir l'estructura interna del ressalt hidràulic i totes les variables derivades millor que FLOW-3D, emprant el TruVOF, tot i que aquest últim pareix capturar millor la transferència de quantitat de moviment i, per tant, totes les variables derivades. En el cas del flux en sobreeixidors esgraonats, OpenFOAM captura millor els perfils de velocitat, tot i que FLOW-3D és més precís en estimar els perfils de làmina lliure d'aigua. Cal deixar palès que ni tan sols la seua resposta a certs paràmetres del model és comparable. Per exemple, FLOW-3D és significativament menys sensible al refinament de malla que OpenFOAM. En base a la precisió dels resultats obtinguts en tots els casos, el model proposat és completament aplicable a casos de disseny més complexos, on vassos esmorteïdors, sobreeixidors esgraonats i estructures hidràuliques en general han de ser investigades. / Bayón Barrachina, A. (2017). Numerical analysis of air-water flows in hydraulic structures using computational fluid dynamics (CFD) [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/90440 / TESIS

NOWADAYS, a number of studies keep on demonstrating the existence of a strong relation between high concentrations of particulate matter (PM) and the prevalence of human morbidity and mortality. Large particles can be filtered in the nose or in the throat, while fine particles (about10 micrometer) can settle in the bronchi and lungs, leading to more serious consequences. According to Karagulian et al. the major sources of urban air pollution are traffic (25%), combustion and agriculture (22%), domestic fuel burning (20%), natural dust (18%) and industrial activities (15%).As a consequence, the detailed study of dispersion phenomena within the urban canopy becomes a target of great interest. To this end, Computational Fluid Dynamics (CFD) can be successfully employed to predict turbulence and dispersion patterns, accounting for a detailed characterization of the pollutant sources, complex obstacles and atmospheric stability classes.Despite being intrinsically different phenomena, turbulence and dispersion are closely related. It is universally accepted that, to reach accurate prediction of the concentration field, it is necessary to properly reproduce the turbulence one. For this reason, the present PhD thesis is split into two main Sections: one focused on turbulence modelling and the subsequent, centered on the dispersion modelling.Thanks to its good compromise between accuracy of results and calculation time, Reynolds-averaged Navier-Stokes (RANS) still represents a valid alternative to more resource-demanding methods. However, focusing on the models’ performance in urban studies, Large Eddy Simulation (LES) generally outperforms RANS results, even if the former is at least one order of magnitude more expensive. Stemming from this consideration, the aim of this work is to propose a variety of approaches meant to solve some of the major limitations linked to standard RANS simulation and to further improve its accuracy in disturbed flow fields, without renouncing to its intrinsic feasibility. The proposed models are suitable for the urban context, being capable of automatically switching from a formulation proper for undisturbed flow fields to one suitable for disturbed areas. For neutral homogeneous atmospheric boundary layer (ABL), a comprehensive approach is adopted, solving the issue of the erroneous stream-wise gradients affecting the turbulent profiles and able to correctly represent the various roughness elements. Around obstacles, more performing closures are employed. The transition between the two treatments is achieved through the definition of a Building Influence Area (BIA). The finalgoal is to offer more affordable alternatives to LES simulations without sacrificing a good grade of accuracy.Focusing on the dispersion modelling framework, there exists a number of parameters which have to be properly specified. In particular, the definition of the turbulent Schmidt number Sct, expressing the ratio of turbulent viscosity to turbulent mass diffusivity, is imperative. Despite its relevance, the literature does not report a clear guideline on the definition of this quantity. Nevertheless, the importance of Sct with respect to dispersion is undoubted and further demonstrated in the works of different authors. For atmospheric boundary layer flows, typical constant values range between 0.2 and 1.3. As a matter of fact, the local variability of Sct is supported by experimental evidence and by direct numerical simulations (DNS). These observations further suggest that the turbulent Schmidt number should be prescribed as a dynamic variable. Following these observations a variable turbulent Schmidt number formulation is proposed in this work. The latter stems from the same hypothesis of the variable formulation developed by Gorlé et al. Moreover, the relevant uncertain model parameters are optimized through uncertainty quantification (UQ). This formulation further increased the accuracy of the predictions, and was successfully verified by Di Bernardino et al. However, the turbulent Schmidt number resulting from this formulation is still intrinsically linked to the turbulence model employed, i.e. to the Cμ coefficient. To overcome this constraint, the nature and the dependencies of Sct were further analyzed through correlation studies and employing principal component analysis (PCA) on data obtained through the proposed ABL RANS model. Subsequently, the same data-driven technique was employed based on the high-fidelity outcomes of a delayed Detached Eddy Simulation (dDES) to derive a generalized turbulentSchmidt number formulation. The latter can be employed within a wide range of turbulence models, without limiting its variability. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

No description available.

In this thesis, we have evaluated a stochastic Lagrangian model for computing particle deposition rates with prospects to use for atmospheric deposition applications.  The model is one-dimensional and models the particle dynamics in the boundary layers near walls and obstacles by simulating the coherent turbulent structures and Brownian motion governing the wall-normal transport. The deposition model is used with a hybrid deterministic/stochastic particle dispersion model governing the dynamics in the turbulent bulk flow. We used a steady-state RANS k-ϵ turbulence model to simulate the turbulent fluid flow in a neutral atmospheric boundary layer (ABL) using the with inflow boundary conditions by Richards & Hoxey (1993). The turbulence model is solved with the SIMPLE algorithm using the OpenFOAM software. The mean-field characteristic of the turbulent flow in the computational domain is exported and used for the particle model. The particle model is a Lagrangian Langevin-type model, consisting of a system of stochastic differential equations. The particle model was solved using a weakly first order a-stable scheme. We evaluated the deposition model by computing the deposition rate for a range of particle sizes and compared our results with collected experimental wind tunnel data. The numerical experiment was done in a computational domain based on the ABL model by Hargreaves & Wright (2007), a rectangular domain with a logarithmic wind profile. We used a particle source near the inflow boundary with an instantaneously release at the initial time. Results showed disagreement with the experimental data and was only valid for medium sized particles. However, time restrictions led to the analysis being cut short and only a single simulation was conducted. A definite conclusion on the suitability of the method could not be made based solely on this single results. Some uncertainties were identified and discussed for further potential work on the evaluation of the method. However, one conclusion was drawn on the performance of the method. The computational cost was concluded to be too high with the first order particle scheme used and higher order schemes is required for any practical use of the method for atmospheric deposition applications.

Two computational fluid dynamics models are developed to predict the performance of a solid oxide fuel cell stack, a detailed and a simplified model. In the detailed model, the three dimensional momentum, heat, and species transport equations are coupled with electrochemistry. In the simplified model, the diffusion terms in the transport equations are selectively replaced by rate terms within the core region of the stack. This allows much coarser meshes to be employed at a fraction of the computational cost. Following the mathematical description of the problem, results for single-cell and multi-cell stacks are presented. Comparisons of local current density, temperature, and cell voltage indicate that good agreement is obtained between the detailed and simplified models, verifying the latter as a practical option in stack design. Then, the simplified model is used to determine the effects of utilization on the electrochemical performance and temperature distributions of a 10 cell stack. The results are presented in terms of fluid flow, pressure, species mass fraction, temperature, voltage and current density distributions. The effects of species and flow distributions on electrochemical performance and temperature are then analyzed for a 100 cell stack. The discussion highlights the importance of manifold design on performance and thermal management of large stacks. / Thesis (Master, Mechanical and Materials Engineering) -- Queen's University, 2013-09-30 15:55:18.627

A breed-and-burn molten salt reactor (BBMSR) concept is proposed to address the Generation IV fuel cycle sustainability objective in a once-through cycle with low enrichment and no reprocessing. The BBMSR uses separate fuel and coolant molten salts, with the fuel contained in assemblies of individual tubes that can be shuffled and reclad periodically to enable high burnup. In this dual-salt configuration, the BBMSR may overcome several limitations of previous breed-and-burn (B$\&$B) designs to achieve high uranium utilisation with a simple, passively safe design. A central challenge in design of the BBMSR fuel is balancing the neutronic requirement of large fuel volume fraction for B$\&$B mode with the thermal-hydraulic requirements for safe and economically competitive reactor operation. Natural convection of liquid fuel within the tubes aids heat transfer to the coolant, and a systematic approach is developed to efficiently model this complex effect. Computational fluid dynamics modelling is performed to characterise the unique physics of the system and produce a new heat transfer correlation, which is used alongside established correlations in a numerical model. A design framework is built around this numerical model to iteratively search for the limiting power density of a given fuel and channel geometry, applying several defined temperature and operational constraints. It is found that the trade-offs between power density, core pressure drop, and pumping power are lessened by directing the flow of coolant downwards through the channel. Fuel configurations that satisfy both neutronic and thermal-hydraulic objectives are identified for natural, 5$\%$ enriched, and 20$\%$ enriched uranium feed fuel. B$\&$B operation is achievable in the natural and 5$\%$ enriched versions, with power densities of 73 W/cm$^3$ and 86 W/cm$^3$, and theoretical uranium utilisations of 300 $\mathrm{MWd/kgU_{NAT}}$ and 25.5 $\mathrm{MWd/kgU_{NAT}}$, respectively. Using 20$\%$ enriched feed fuel relaxes neutronic constraints so a wider range of fuel configurations can be considered, but there is a strong inverse correlation between power density and uranium utilisation. The fuel design study demonstrates the flexibility of the BBMSR concept to operate along a spectrum of modes ranging from high fuel utilisation at moderate power density using natural uranium feed fuel, to high power density and moderate utilisation using 20$\%$ uranium enrichment.

This paper covers information research of basic design rules of industrial wet sumps. It describes mathematical models of vortices and method for their identification and visualization. Then the author focuses on CFD modeling of surface vortices with single phase and multiphase approach with Volume of Fluid method. Basic principles of multiphase CFD modelling in OpenFOAM and ANSYS Fluent are given. Description and benchmarking of suitable turbulence models is also present. The single phase and multiphase approach were successfully validated for a simple test case of bathtub surface vortex. Satisfactory agreement with experimental data was achieved. The accuracy and behavior of both solvers were compared between each other. This gives us useful tool for evaluation of inflow condition and danger of surface vortex occurrence in wet sumps. The acquired knowledges were used to design an experimental test case with geometry similar to industrial wet sump. A map of surface vortex occurrence has been created for different operating points. One of the operating point has been used for numerical simulation (both single phase and multiphase). Partial agreement with experimental observation has been achieved.

In the event of a radioactive release in an urban environment the consequent response mustbe swift and precise. As soon as first responders have correct information, they can make anaccurate risk assessment. However, if the position, release rate and time of the radioactiverelease is unknown it is hard to know how the pollutant will spread. This thesis aims to testa model which approximates these three unknowns using weather data (wind and rain) as wellas measurement data collected at sensors placed around an urban environment. An atmospheric dispersion model based on an existing Reynolds Averaged Navier-Stokes modelis set up in two geometries of different complexity to create forward mode synthetic depositiondata and adjoint mode concentration fields resulting from a fixed dry deposition velocity andscavenging effect for wet deposition. Variations of time- and space-dependent rainfall is simu-lated. The resulting data is used in an existing optimization model, where a parameter studyis conducted regarding regularization coefficients. This thesis shows that the optimization model accurately estimates position and its approximaterelease rate of a 2D geometry of radioactive releases using a logarithmic optimization approach,and fail to do so using a linear optimization approach. The logarithmic optimization model alsoapproximately estimates position and release rate in a 3D geometry. Regularization parametersshould be within the range of 0.1 and 1.2 depending on rain. More rain requires smallerparameters and will estimate a lower release rate. Time-dependent rainfall is shown to have amajor negative effect on simulation time.iii

The study and modelling of two-phase flow, even the simplest ones such as the bubbly flow, remains a challenge that requires exploring the physical phenomena from different spatial and temporal resolution levels. CFD (Computational Fluid Dynamics) is a widespread and promising tool for modelling, but nowadays, there is no single approach or method to predict the dynamics of these systems at the different resolution levels providing enough precision of the results. The inherent difficulties of the events occurring in this flow, mainly those related with the interface between phases, makes that low or intermediate resolution level approaches as system codes (RELAP, TRACE, ...) or 3D TFM (Two-Fluid Model) have significant issues to reproduce acceptable results, unless well-known scenarios and global values are considered. Instead, methods based on high resolution level such as Interfacial Tracking Method (ITM) or Volume Of Fluid (VOF) require a high computational effort that makes unfeasible its use in complex systems. In this thesis, an open-source simulation framework has been designed and developed using the OpenFOAM library to analyze the cases from microescale to macroscale levels. The different approaches and the information that is required in each one of them have been studied for bubbly flow. In the first part, the dynamics of single bubbles at a high resolution level have been examined through VOF. This technique has allowed to obtain accurate results related to the bubble formation, terminal velocity, path, wake and instabilities produced by the wake. However, this approach has been impractical for real scenarios with more than dozens of bubbles. Alternatively, this thesis proposes a CFD Discrete Element Method (CFD-DEM) technique, where each bubble is represented discretely. A novel solver for bubbly flow has been developed in this thesis. This includes a large number of improvements necessary to reproduce the bubble-bubble and bubble-wall interactions, turbulence, velocity seen by the bubbles, momentum and mass exchange term over the cells or bubble expansion, among others. But also new implementations as an algorithm to seed the bubbles in the system have been incorporated. As a result, this new solver gives more accurate results as the provided up to date. Following the decrease on resolution level, and therefore the required computational resources, a 3D TFM have been developed with a population balance equation solved with an implementation of the Quadrature Method Of Moments (QMOM). The solver is implemented with the same closure models as the CFD-DEM to analyze the effects involved with the lost of information due to the averaging of the instantaneous Navier-Stokes equation. The analysis of the results with CFD-DEM reveals the discrepancies found by considering averaged values and homogeneous flow in the models of the classical TFM formulation. Finally, for the lowest resolution level approach, the system code RELAP5/MOD3 is used for modelling the bubbly flow regime. The code has been modified to reproduce properly the two-phase flow characteristics in vertical pipes, comparing the performance of the calculation of the drag term based on drift-velocity and drag coefficient approaches. / El estudio y modelado de flujos bifásicos, incluso los más simples como el bubbly flow, sigue siendo un reto que conlleva aproximarse a los fenómenos físicos que lo rigen desde diferentes niveles de resolución espacial y temporal. El uso de códigos CFD (Computational Fluid Dynamics) como herramienta de modelado está muy extendida y resulta prometedora, pero hoy por hoy, no existe una única aproximación o técnica de resolución que permita predecir la dinámica de estos sistemas en los diferentes niveles de resolución, y que ofrezca suficiente precisión en sus resultados. La dificultad intrínseca de los fenómenos que allí ocurren, sobre todo los ligados a la interfase entre ambas fases, hace que los códigos de bajo o medio nivel de resolución, como pueden ser los códigos de sistema (RELAP, TRACE, etc.) o los basados en aproximaciones 3D TFM (Two-Fluid Model) tengan serios problemas para ofrecer resultados aceptables, a no ser que se trate de escenarios muy conocidos y se busquen resultados globales. En cambio, códigos basados en alto nivel de resolución, como los que utilizan VOF (Volume Of Fluid), requirieren de un esfuerzo computacional tan elevado que no pueden ser aplicados a sistemas complejos. En esta tesis, mediante el uso de la librería OpenFOAM se ha creado un marco de simulación de código abierto para analizar los escenarios desde niveles de resolución de microescala a macroescala, analizando las diferentes aproximaciones, así como la información que es necesaria aportar en cada una de ellas, para el estudio del régimen de bubbly flow. En la primera parte se estudia la dinámica de burbujas individuales a un alto nivel de resolución mediante el uso del método VOF (Volume Of Fluid). Esta técnica ha permitido obtener resultados precisos como la formación de la burbuja, velocidad terminal, camino recorrido, estela producida por la burbuja e inestabilidades que produce en su camino. Pero esta aproximación resulta inviable para entornos reales con la participación de más de unas pocas decenas de burbujas. Como alternativa, se propone el uso de técnicas CFD-DEM (Discrete Element Methods) en la que se representa a las burbujas como partículas discretas. En esta tesis se ha desarrollado un nuevo solver para bubbly flow en el que se han añadido un gran número de nuevos modelos, como los necesarios para contemplar los choques entre burbujas o con las paredes, la turbulencia, la velocidad vista por las burbujas, la distribución del intercambio de momento y masas con el fluido en las diferentes celdas por cada una de las burbujas o la expansión de la fase gaseosa entre otros. Pero también se han tenido que incluir nuevos algoritmos como el necesario para inyectar de forma adecuada la fase gaseosa en el sistema. Este nuevo solver ofrece resultados con un nivel de resolución superior a los desarrollados hasta la fecha. Siguiendo con la reducción del nivel de resolución, y por tanto los recursos computacionales necesarios, se efectúa el desarrollo de un solver tridimensional de TFM en el que se ha implementado el método QMOM (Quadrature Method Of Moments) para resolver la ecuación de balance poblacional. El solver se desarrolla con los mismos modelos de cierre que el CFD-DEM para analizar los efectos relacionados con la pérdida de información debido al promediado de las ecuaciones instantáneas de Navier-Stokes. El análisis de resultados de CFD-DEM permite determinar las discrepancias encontradas por considerar los valores promediados y el flujo homogéneo de los modelos clásicos de TFM. Por último, como aproximación de nivel de resolución más bajo, se investiga el uso uso de códigos de sistema, utilizando el código RELAP5/MOD3 para analizar el modelado del flujo en condiciones de bubbly flow. El código es modificado para reproducir correctamente el flujo bifásico en tuberías verticales, comparando el comportamiento de aproximaciones para el cálculo del término d / L'estudi i modelatge de fluxos bifàsics, fins i tot els més simples com bubbly flow, segueix sent un repte que comporta aproximar-se als fenòmens físics que ho regeixen des de diferents nivells de resolució espacial i temporal. L'ús de codis CFD (Computational Fluid Dynamics) com a eina de modelatge està molt estesa i resulta prometedora, però ara per ara, no existeix una única aproximació o tècnica de resolució que permeta predir la dinàmica d'aquests sistemes en els diferents nivells de resolució, i que oferisca suficient precisió en els seus resultats. Les dificultat intrínseques dels fenòmens que allí ocorren, sobre tots els lligats a la interfase entre les dues fases, fa que els codis de baix o mig nivell de resolució, com poden ser els codis de sistema (RELAP,TRACE, etc.) o els basats en aproximacions 3D TFM (Two-Fluid Model) tinguen seriosos problemes per a oferir resultats acceptables , llevat que es tracte d'escenaris molt coneguts i se persegueixen resultats globals. En canvi, codis basats en alt nivell de resolució, com els que utilitzen VOF (Volume Of Fluid), requereixen d'un esforç computacional tan elevat que no poden ser aplicats a sistemes complexos. En aquesta tesi, mitjançant l'ús de la llibreria OpenFOAM s'ha creat un marc de simulació de codi obert per a analitzar els escenaris des de nivells de resolució de microescala a macroescala, analitzant les diferents aproximacions, així com la informació que és necessària aportar en cadascuna d'elles, per a l'estudi del règim de bubbly flow. En la primera part s'estudia la dinàmica de bambolles individuals a un alt nivell de resolució mitjançant l'ús del mètode VOF. Aquesta tècnica ha permès obtenir resultats precisos com la formació de la bambolla, velocitat terminal, camí recorregut, estela produida per la bambolla i inestabilitats que produeix en el seu camí. Però aquesta aproximació resulta inviable per a entorns reals amb la participació de més d'unes poques desenes de bambolles. Com a alternativa en aqueix cas es proposa l'ús de tècniques CFD-DEM (Discrete Element Methods) en la qual es representa a les bambolles com a partícules discretes. En aquesta tesi s'ha desenvolupat un nou solver per a bubbly flow en el qual s'han afegit un gran nombre de nous models, com els necessaris per a contemplar els xocs entre bambolles o amb les parets, la turbulència, la velocitat vista per les bambolles, la distribució de l'intercanvi de moment i masses amb el fluid en les diferents cel·les per cadascuna de les bambolles o els models d'expansió de la fase gasosa entre uns altres. Però també s'ha hagut d'incloure nous algoritmes com el necessari per a injectar de forma adequada la fase gasosa en el sistema. Aquest nou solver ofereix resultats amb un nivell de resolució superior als desenvolupat fins la data. Seguint amb la reducció del nivell de resolució, i per tant els recursos computacionals necessaris, s'efectua el desenvolupament d'un solver tridimensional de TFM en el qual s'ha implementat el mètode QMOM (Quadrature Method Of Moments) per a resoldre l'equació de balanç poblacional. El solver es desenvolupa amb els mateixos models de tancament que el CFD-DEM per a analitzar els efectes relacionats amb la pèrdua d'informació a causa del promitjat de les equacions instantànies de Navier-Stokes. L'anàlisi de resultats de CFD-DEM permet determinar les discrepàncies ocasionades per considerar els valors promitjats i el flux homogeni dels models clàssics de TFM. Finalment, com a aproximació de nivell de resolució més baix, s'analitza l'ús de codis de sistema, utilitzant el codi RELAP5/MOD3 per a analitzar el modelatge del fluxos en règim de bubbly flow. El codi és modificat per a reproduir correctament les característiques del flux bifàsic en canonades verticals, comparant el comportament d'aproximacions per al càlcul del terme de drag basades en velocitat de drift flux model i de les basades en coe / Peña Monferrer, C. (2017). Computational fluid dynamics multiscale modelling of bubbly flow. A critical study and new developments on volume of fluid, discrete element and two-fluid methods [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/90493 / TESIS