This thesis aims to provide the thermal condition of melt pool convection by CFD simulation, which is important to the assessment of the invessel melt retention (IVR) strategy widely adopted in Generation III pressurized water reactors (PWRs). As a severe accident mitigation measure, the IVR strategy is realized through external cooling of the lower head of a reactor pressure vessel (RPV). To achieve the coolability and retention of the corium pool in the RPV lower head, the heat flux at the outer surface of the vessel should be less than the critical heat flux (CHF) of boiling around the lower head. Under such condition, the integrity of the RPV is guaranteed by the adequate thickness of the unmelted vessel wall. The thesis work starts from the selection and validation of a turbulence model in the CFD computational tool chosen (Fluent). Afterwards a numerical model is set up for estimation of melt pool heat transfer of a reference PWR with the power capacity of 1000 MWe, including a mesh sensitivity study. Based on the numerical model of a twolayer melt pool, four tasks are carried out to investigate the effects of Zr oxidation ratio, Fe content, and radiation emissivity on heat flux profiles, as well as the focus effect under extreme conditions. Selection and validation of the turbulence model are conducted by comparing the simulation results of different turbulence models with the DNS data on the convection of volumetrically heated fluid layer bounded by rigid isothermal horizontal walls at equal temperature. The internal Rayleigh numbers of the flow reach up to 10e6. The comparison shows a good agreement of the SST k-ω turbulence model results with the DNS data. The simulations with the Zr oxidation ratio of 0, 0.2 and 0.5, correspondingly, the oxide layer of 1.389m, 1.467m and 1.580m, and the metal layer of 0.705m, 0.646m and 0.561m in height, show that, the temperature of the oxide layer will increase with Zr oxidation ratio, while the temperature of the metal layer will decrease resulting in more heat transfer through the oxide layer sidewall and less top radiation. Nevertheless, the effect of the Zr oxidation ratio is not pronounced in the range of 00.5. The simulations with the Fe mass of 22t, 33t and 45t, and respective height of the metal layer of 0.462m, 0.568m and 0.646m, show that, the inner metal layer will significantly increase the temperatures of both the metal layer and the oxide layer. The percentage of heat transfer at the oxide layer sidewall will increase to supplement the reduction of that at the metal layer. The simulations with the radiation emissivity of 0.2, 0.35, 0.45 and 0.7 show that, the emissivity below 0.45 has an impact on heat transfer, and the temperatures and sidewall heat flux of both the oxide layer and the metal layer will increase with decreasing emissivity. The impact is negligible when the emissivity is above 0.45. The simulations under the hypothetically extreme conditions with either an adiabatic top boundary or a very thin metal layer show the focusing effect may occur, i.e., the heat flux through the metal sidewall is larger than that in the oxide layer. But the local high heat flux is flattened by the vessel wall with good heat conductivity. In summary, the simulations demonstrate that, except for the cases under extreme conditions, the heat fluxes of the melt pools in all other cases are significantly lower than the CHF of external cooling of the lower head. Therefore, the safety margin of the IVR strategy of the PWR chosen is seems sufficient. However, due to some limitations (e.g., simplification and assumptions) in the simulation cases and coupling of different influential factors, as indicated by the present study, the precise predictions of heat flux under all scenarios are still difficult. Therefore, the conclusions could not be generalized to the other conditions or other configurations of the molten pools. By discussing the model and simplifications/assumptions adopted in this work, the improvement directions of the numerical model and other perspectives are proposed at the end of the thesis. / Denna avhandling syftar till att tillhandahålla det termiska tillståndet för smältbassängskonvektion genom CFD-simulering, vilket är viktigt för bedömningen av IVR-strategin som allmänt antagits i tryckvattenreaktorer (PWR) i Generation III. Som en åtgärd för att mildra allvarliga olyckor realiseras IVR-strategin genom extern kylning av det nedre huvudet av ett reaktortryckkärl (RPV). För att uppnå kylbarhet och kvarhållning av koriumbassängen i det nedre RPV-huvudet bör värmeflöde vid den yttre ytan av kärlet vara mindre än det kritiska värmeflödet (CHF) som kokar runt det nedre huvudet. Under sådant tillstånd garanteras RPV: s integritet av den osmälta kärlväggens tillräckliga tjocklek. Examensarbetet startar från valet och valideringen av en turbulensmodell i det valda CFD-beräkningsverktyget (Fluent). Därefter sätts en numerisk modell upp för uppskattning av smältbassängens värmeöverföring av en referens PWR med en effektkapacitet på 1000 MWe, inklusive en nätkänslighetsstudie. Baserat på den numeriska modellen för en tvålagers smältbassäng utförs fyra uppgifter för att undersöka effekterna av Zr-oxidationsförhållande, Fe-innehåll och strålningsemissivitet på värmeflödesprofiler, liksom fokuseffekten under extrema förhållanden. Val och validering av turbulensmodellen utförs genom att jämföra simuleringsresultaten för olika turbulensmodeller med DNS-data för konvektionen av volymetriskt uppvärmt fluidskikt avgränsat av styva isoterma horisontella väggar vid lika temperatur. De interna Rayleigh-siffrorna i flödet når upp till 10e6. Jämförelsen visar att SST k-ω turbulensmodellresultaten överensstämmer med DNS-data. Simuleringarna med Zr-oxidationsförhållandet 0, 0,2 och 0,5, motsvarande oxidskiktet på 1,389 m, 1,467 m och 1,580 m, och metallskiktet på 0,705 m, 0,664 m och 0,561 m i höjd, visar att temperaturen av oxidskiktet kommer att öka med Zr-oxidationsförhållandet, medan metallskiktets temperatur kommer att minska vilket resulterar i mer värmeöverföring genom oxidskiktets sidovägg och mindre toppstrålning. Ändå är effekten av Zr-oxidationsförhållandet inte uttalad i intervallet 00,5. Simuleringarna med Fe-massan på 22t, 33t och 45t och respektive höjd av metallskiktet på 0,462m, 0,568m och 0,664m visar att det inre metallskiktet avsevärt kommer att öka temperaturerna för både metallskiktet och oxiden lager. Andelen värmeöverföring vid oxidskiktets sidovägg ökar för att komplettera minskningen av den vid metallskiktet. Simuleringarna med strålningsemissiviteten 0,2, 0,35, 0,45 och 0,7 visar att emissiviteten under 0,45 påverkar värmeöverföringen, och temperaturerna och sidoväggens värmeflöde för både oxidskiktet och metallskiktet kommer att öka med minskande emissivitet. Effekten är försumbar när strålningen är över 0,45. Simuleringarna under de hypotetiskt extrema förhållandena med antingen en adiabatisk övre gräns eller ett mycket tunt metallskikt visar att fokuseringseffekten kan uppstå, dvs. värmeflödet genom metallsidan är större än det i oxidskiktet. Men det lokala höga värmeflödet plattas ut av kärlväggen med god värmeledningsförmåga. Sammanfattningsvis visar simuleringarna att, förutom fall under extrema förhållanden, är värmeflödet från smältpoolerna i alla andra fall betydligt lägre än CHF för extern kylning av nedre huvudet. Därför verkar säkerhetsmarginalen för IVR-strategin för den valda PWR tillräcklig. På grund av vissa begränsningar (t.ex. förenkling och antaganden) i simuleringsfall och koppling av olika inflytelserika faktorer, vilket indikeras av den aktuella studien, är de exakta förutsägelserna av värmeflöde under alla scenarier fortfarande svåra. Därför kunde slutsatserna inte generaliseras till de andra förhållandena eller andra konfigurationer av de smälta poolerna. Genom att diskutera modellen och förenklingar / antaganden som antagits i detta arbete föreslås förbättringsriktningarna för den numeriska modellen och andra perspektiv i slutet av avhandlingen.

The turbine inlet guide vanes of a gas-turbine are subjected to extreme hot gas temperatures which increases the risk of mechanical failure and overall reduces the component lifespan. Hence, it is of great interest for gas-turbine manufacturers to establish methods for accurately estimating the temperature distribution along the vane surface. Due to the three-dimensional nature of turbine flow, it is of interest to establish Computational Fluid Dynamics (CFD) methodology which capture these three-dimensional effects. This thesis is one in a collection of theses conducted at Siemens Energy AB on the subject. Previous studies have investigated and validated the implementation of RANS simulations on non-cooled turbine vanes and endwalls. In this study, the focus is on studying a film cooled vane and establishing one RANS as well as one hybrid modelling strategy for heat transfer coefficient (HTC) predictions. The HTC prediction capabilities are compared and validated against experimental data presented in the doctoral thesis by Fabrizio Fontaneto on the LS94 vane at Von Karman Institute. The chosen RANS modelling method was the Shear Stress Transport (SST) k-ω turbulence model, with γ-Reθ transition modelling, based on the findings by Enico (2021) and Daugulis (2022). The model proved capable in estimating the HTC well on mainly the suction side of the vane. The pressure side HTC was largely under-predicted, a common issue with the SST model also seen in the previous theses as well as the hybrid simulations. The strength of the SST k-ω turbulence model, with γ-Reθ transition modelling, is in accurately capturing the HTC magnitude, most likely due to the well-predicted turbulence intensity decay at the inlet. However, it largely under-predicts the HTC along the suction side film-coolant layer, implying that it may be over-estimating the film-cooling capabilities. The hybrid model chosen was the Scale Resolving Hybrid (SRH) model, with underlying RANS SST k-ω. Compared to RANS, hybrid results were under-estimated, seemingly offset from the experimental data by a constant 200 units along the entire vane midspan. This is likely due to the inaccurate turbulence intensity presented in the SRH simulations, which decays quickly along the inlet compared to RANS and experimental data. Yet still, the hybrid model showed potential in capturing certain results not seen with RANS, such as the secondary flow effects by the vane endwalls, as well as arguably capturing the general HTC trend at midspan seen in the experimental data. Additionally, the section of severely under-predicted HTC by the suction side film-coolant seen with RANS is not present in the hybrid results. Although the hybrid model has proven promising in many aspects, in its current state it is not a viable method for HTC predictions due to its general under-prediction of HTC. Largely, the authors suspect this is due to the undesirably coarse mesh around the cooling holes, which leads to RANS computation in regions where SRH is desired. Thus, improvements would need to be made to the model, where, for example, implementing a zonal hybrid RANS-LES model would be an option. Considering the hybrid model in its current state, RANS is the preferred method, especially when considering the greater computational cost and the labor associated with hybrid simulations which were experienced during this study. In conclusion, it is evident that the correct capture of inlet turbulence intensity decay as well as suitable mesh refinement by the cooling holes are crucial for obtaining the correct magnitudes of HTC, and thus, the capture of it should be of utmost priority in future work within the field.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第21723号 / 工博第4540号 / 新制||工||1708(附属図書館) / 京都大学大学院工学研究科社会基盤工学専攻 / (主査)教授 中川 一, 准教授 竹林 洋史, 准教授 川池 健司 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

<p>More and more attention has been paid to decreasing the number and size of non-metallic inclusions existing in the final products recently in steel industries. Therefore, more efforts have been made to monitor the inclusions' size distributions during the metallurgy process, especially at the secondary steelmaking period. A liquid sampling procedure is one of the commonly applied methods that monitoring the inclusion size distribution in ladles, for example, during the secondary steelmaking. Here, a crucial point is that the steel sampler should be filled and solidified without changing the inclusion characteristics that exist at steel making temperatures. In order to preserve the original size and distributions in the extracted samples, it is important to avoid their collisions and coagulations inside samplers during filling. Therefore, one of the first steps to investigate is the flow pattern inside samplers during filling in order to obtain a more in-depth knowledge of the sampling process to make sure that the influence is minimized.</p><p>The main objective of this work is to fundamentally study the above mentioned sampler filling process. A production sampler employed in the industries has been scaled-up according to the similarity of Froude Number in the experimental study. A Particle Image Velocimetry (PIV) was used to capture the flow field and calculate the velocity vectors during the entire experiment. Also, a mathematical model has been developed to have an in-depth investigate of the flow pattern in side the sampler during its filling. Two different turbulence models were applied in the numerical study, the realizable k-ε model and Wilcox k-ω model. The predictions were compared to experimental results obtained by the PIV measurements. Furthermore, it was illustrated that there is a fairly good agreement between the measurements obtained by PIV and calculations predicted by the Wilcox k-ω model. Thus, it is concluded that the Wilcox k-ω model can be used in the future to predict the filling of steel samplers.</p>

This document is a report summering the master thesis work dealing with the Computational Fluid Dynamic (CFD) study of the flow around a high-speed train. The model is a scaled 1:50 generic train with two cars, one inter-car gap and simplified bogies. A platform is set on the side of the train since one of the aim of the study is to look at the consequences of the phenomena in the wake on people or objects standing on the platform. The slipstream is one of this phenomena, it is due to the fact that the viscous air is dragged when the train is passing. If too strong, it can move or destabilize people or objects on the platform. In addition of the slipstream study, a velocity profile study, a drag and lift coefficients analyze as well as a Q-factor study and a frequency study have been realized. Some results of these different studies are compared with the ones obtained on the same model with a Delayed Detached Eddy Simulation (DDES). Since the flow is turbulent, for those different studies, the flow has been simulated with a Reynolds Averaged Navier-Stokes equation model (RANS) which is the k-ω SST model for the turbulence. The study of the slipstream allowed to calculate the Technical Specification for Interoperability (TSI) which must not be higher that the European Union requirement set at 15.5 m/s, the result obtained is 8.1 m/s which is then lower than the limit. The velocity profile shows similarities with the DDES results even though it is less detailed. The same conclusion is done for the Q-plot where is clearly visible the two counter-rotating vortices in the wake. Finally, a Fast Fourier Transform algorithm has been applied to instantaneous velocity results in the wake of the train in order to get the frequency of the aerodynamic phenomena in that wake. The main frequency is 25 Hz and corresponds to a Strouhal number of 0.1, quite closed to the results obtained with DDES which is 0.085. The results of the RANS and DDES are reasonably similar and by regarding at the large difference between the cell numbers (respectively 8 500 000 and 20 000 000) it can be conclude that in some ways the RANS model can be preferred at the DDES to save time for the computation but it does not contain the small scales resolved by the DDES.

More and more attention has been paid to decreasing the number and size of non-metallic inclusions existing in the final products recently in steel industries. Therefore, more efforts have been made to monitor the inclusions' size distributions during the metallurgy process, especially at the secondary steelmaking period. A liquid sampling procedure is one of the commonly applied methods that monitoring the inclusion size distribution in ladles, for example, during the secondary steelmaking. Here, a crucial point is that the steel sampler should be filled and solidified without changing the inclusion characteristics that exist at steel making temperatures. In order to preserve the original size and distributions in the extracted samples, it is important to avoid their collisions and coagulations inside samplers during filling. Therefore, one of the first steps to investigate is the flow pattern inside samplers during filling in order to obtain a more in-depth knowledge of the sampling process to make sure that the influence is minimized. The main objective of this work is to fundamentally study the above mentioned sampler filling process. A production sampler employed in the industries has been scaled-up according to the similarity of Froude Number in the experimental study. A Particle Image Velocimetry (PIV) was used to capture the flow field and calculate the velocity vectors during the entire experiment. Also, a mathematical model has been developed to have an in-depth investigate of the flow pattern in side the sampler during its filling. Two different turbulence models were applied in the numerical study, the realizable k-ε model and Wilcox k-ω model. The predictions were compared to experimental results obtained by the PIV measurements. Furthermore, it was illustrated that there is a fairly good agreement between the measurements obtained by PIV and calculations predicted by the Wilcox k-ω model. Thus, it is concluded that the Wilcox k-ω model can be used in the future to predict the filling of steel samplers.

El cultivo protegido bajo invernadero tiene altos costes energéticos derivados de la ventilación forzada. En cambio, la ventilación natural puede ser una solución barata que reduzca el consumo energético. No obstante, un diseño que permita la ventilación natural del invernadero supone un reto debido a la complejidad del fenómeno físico. Frente a otros métodos de diseño, como los métodos experimentales o analíticos, en los últimos tiempos la dinámica de fluidos computacional (CFD por sus siglas en inglés) se ha convertido en la herramienta más utilizada para estudiar este tipo de fenómenos, gracias a su relativo bajo coste y a la rapidez en la obtención de resultados. No obstante, los modelos CFD deben refrendarse mediante validaciones realizadas a través de datos experimentales. Un análisis bibliográfico detallado del uso de la CFD aplicada a invernaderos muestra, en general, que los trabajos carecen de suficientes datos experimentales, seguramente debido al alto coste de los sensores para adquirirlos y la dificultad que conlleva el trabajo en condiciones de campo. Además, se observa que las simulaciones CFD en invernaderos no se realizan con un procedimiento sistematizado. Por ello, la presente tesis, por un lado describe un dispositivo y unos métodos sencillos y baratos para obtener datos atmosféricos, y por otro, propone una visión crítica sobre la investigación realizada hasta el momento, con el fin de sistematizar la manera de generar modelos CFD aplicados a la ventilación natural de invernaderos. Finalmente, la tesis se complementa con un ejemplo sobre un caso práctico. Para ello, en primer lugar, se realizó una revisión bibliográfica de las diferentes guías de buenas prácticas en diferentes campos de la tecnología principalmente en edificación, para sistematizar y adaptar las recomendaciones para generar modelos CFD en invernaderos. En segundo lugar, se desarrolló un sistema de adquisición de datos sencillo, consistente en una red de sensores, que permite medir simultáneamente la velocidad y dirección de viento en 20 puntos. Este sistema de sensores fue calibrado y probado en campo satisfactoriamente obteniendo una precisión similar a los anemómetros comerciales con un precio 30 veces superior. En tercer lugar, se generaron 24 modelos del flujo de aire alrededor del invernadero, resultantes de la combinación de cuatro modelos de turbulencia (k- ¿; RNG k-¿; k-¿ y RSM); dos esquemas de discretización (primer orden y segundo orden) y tres velocidades de viento exterior (3; 3,5 y 4 m/s) con el fin de analizar sus diferencias y demostrar sus ventajas e inconvenientes. Por este motivo, en cuarto lugar, se comprobó su capacidad de ajustarse a los datos de campo, validando los modelos con un análisis de regresión lineal sobre los datos experimentales. Con este estudió se reveló que los modelos SST k-¿ y el RSM (segundo orden) son los que mejor representan el flujo de ventilación y se demostró que el modelo k-¿ estándar (primer orden), el más utilizado en la bibliografía, no sólo ofrece resultados diferentes al resto de modelos, sino que su rendimiento es pobre para predecir el flujo de ventilación. / Granell Ruiz, R. (2014). Análisis del flujo ambiental y propuesta metodológica para simulaciones CFD aplicadas a la ventilación natural de invernaderos [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/37194

In recent years, application of confluent jets for design of ventilation supply devices has been studied. Similarly, numerus studies have been made on the potential and application of variable air volume (VAV) in order to reduce the energy demand of ventilation systems. This study investigates the combination of supply devices based on confluent jets and VAV, both in terms of the nearfield flow behavior of the device and the impact on thermal comfort, indoor air quality and energy efficiency on a classroom-level space when the airflow rate is varied. The method used in this study is an experimental field study where the confluent jets-based supply devices were compared to the previously installed displacement ventilation. The field study evaluated the energy efficiency, thermal comfort and indoor air quality of the two systems. In the case of the confluent jets supply devices, airflow rate was varied in order to see what impact the variation had on the performance of the system for each airflow rate. Furthermore, the confluent jets supply devices were investigated both experimentally and numerically in a well insulated test room to get high resolution data on the particular flow characteristics for this type of supply device when the airflow rate is varied. The results from the field study show nearly uniform distribution of the local mean age of air in the occupied zone, even in the cases of relatively low airflow rates. The airflow rates have no significant effect on the degree of mixing. The thermal comfort in the classroom was increased when the airflow rate was adapted to the heat load compared to the displacement system. The results lead to the conclusion that the combination of supply devices based on confluent jets can reduce energy usage in the school while maintaining indoor air quality and increasing the thermal comfort in the occupied zone. The results from the experimental and numerical study show that the flow pattern and velocity in each nozzle is directly dependent on the total airflow rate. However, the flow pattern does not vary between the three different airflow rates. The numerical investigation shows that velocity profiles for each nozzle have the same pattern regardless of the airflow rate, but the magnitude of the velocity profile increases as the airflow increases. Thus, a supply device of this kind could be used for variable air volume and produce confluent jets for different airflow rates. The results from both studies show that the airflow rate does not affect the distribution of the airflow on both near-field and room level. The distribution of air is nearly uniform in the case of the near-field results and the room-level measurement shows a completely uniform degree of mixing and air quality in the occupied zone for each airflow rate. This means that there is potential for combining these two schemes for designing air distribution systems with high energy efficiency and high thermal comfort and indoor air quality. / Under senare tid har applikation av Confluent jets för design av tilluftsdon studerats. Många studier har även utförts över potentialen av att applicera variabelt luftflöde (VAV) för att minska energianvändningen i ventilationssystem. Denna studie undersöker möjligheten att kombinera Confluent jets-don med VAV, både med avseende på den lokala flödesbilden och dess påverkan på termisk komfort, luftkvalitet och energieffektivitet i en klassrumsmiljö där luftflödes varieras. Denna studie baseras dels på en experimentell fältstudie där tilluftsdon baserade på Confluents jets jämfördes med befintliga deplacerande tilluftsdon. Fältstudien utvärderade energieffektiviteten, den termiska komforten och luftkvaliteten för båda typerna av tillluftsdon. Confluent jets-donen testades under varierat luftflöde för att se påverkan av flödesvariationen på ventilationens prestation under de olika flödena. Utöver fältstudien testades Confluent jets-donen experimentellt och numeriskt i ett välisolerat test-rum för få den detaljerade flödeskarakteristiken för den här typen tilluftsdon vid varierat luftflöde. Resultaten från fältstudien visar på en jämn fördelning av den lokala luftsmedelåldern i vistelsezonen, även för fallen med relativt låga luftflöden. Luftflöden har ingen signifikant effekt på omblandningen. Den termiska komforten i klassrummet ökade när luftflödet anpassades efter värmelasten jämfört med de deplacerande donen. Slutsatsen från fältstudien är att kombinationen av VAV och Confluent jets-don kan användas för att minska energianvändningen på skolan och bevara luftkvaliteten och den termiska komforten i vistelsezonen. Resultaten från den experimental och numeriska studien visar luftflödet och lufthastigheten i varje enskild dysa är direkt beroende på det totala luftflödet genom donet. Dock är flödesfördelningen mellan dysorna oberoende av de tre olika luftflödena. Den numeriska undersökningen visar att flödesprofilen för varje dysa är konstant trots att flödet varieras, men amplituden för varje profil ökar med en höjning av luftflödet. Det betyder att tilluftsdon av den här typen kan användas med VAV för att producera Confluent jets för olika luftflöden. Resultaten från båda studierna visar att luftflöde inte påverkar fördelningen av luften vare sig längs luftdonen eller på rumsnivå. Fördelningen av luften är nästan helt jämn längs donen och på rumsnivå är omblandningen och luftkvalitet den samma för varje luftflöde. Det betyder att det finns potential för att kombinera det här två teknikerna för att designa luftdistribueringssystem med hög energieffektivitet och hög termisk komfort med god luftkvalitet.

With steadily increasing computational power, computational fluid dynamics (CFD) can be applied to unsteady problems such as seakeeping simulations. Therefore, a good balance between accuracy and computational speed is required. This thesis investigates the application of CFD to seakeeping performance predictions and aims to propose a best-practice procedure for efficient seakeeping simulations. The widely used KVLCC2 research vessel serves as a test case for this thesis and FINEŠ/Marine software package is used for CFD computations. In order to validate the simulations, results are compared to recent experimental data from SSPA as well as predictions with potential ˛ow code SHIPFLOW® Motions. As for the calm water simulations, both inviscid and viscous ˛ow computations are performed in combination with three mesh refinement levels. Seakeeping simulations with regular head waves of different wavelengths are set-up correspondingly. Furthermore, different strategies for time discretization are investigated. With the given computational resources, it is not feasible to complete seakeeping simulations with a ˝ne mesh. However, already the coarse meshes give good agreement to experiments and SHIPFLOW® Motions' predictions. Viscous ˛ow simulations turn out to be more robust than Euler ˛ow computations and thus should be preferred. Regarding the time discretization, a fixed time discretization of 150 steps per wave period has shown the best balance between accuracy and speed. Based on these findings, a best-practice procedure for seakeeping performance predictions in FINEŠ/Marine is established. Taking the most efficient settings obtained from head wave simulations, the vessel is subjected to oblique waves with 160° encounter angle. Under similar wave conditions, CFD predictions of a similar thesis show close agreement in terms of added wave resistance. Compared to the previous head wave conditions of this study, added resistance in 160° oblique waves is found to be significantly higher. This underlines that oblique bow quartering waves represent a relevant case for determining the maximum required power of a ship. CFD and potential ˛ow show similar accuracy with respect to ship motions and added wave resistance, albeit potential ˛ow outperforms CFD in terms of computational speed. Hence, CFD should be applied in cases where viscous effects are known to have large influence on a vessel's seakeeping behavior. This can be the case if motion control and damping devices are to be evaluated, for instance. / Tack vare den stadigt ökande beräkningskraften kan beräkningsuiddynamik (CFD) idag användas på beräkningsintensiva problem som sjöegenskapssimulationer. Den här rapporten undersöker användning av CFD på sjöegenskapsprestanda och syftar till att foreslå ett best-practice förfaringssätt för effektiv sjöegenskapssimulationer. Forskningsskrovet KVLCC2 fungerar som ett testfall för denna rapport och FINE—/Marine-mjukvarupaketet används för CFD-beräkningar. Viktiga parametrar, såsom ödestyp, beräkningsnät och tidssteg varierars systematiskt. Resultaten jämförs med experiment gjorda vid SSPA. Baserat på resultaten förelås en best-practice. Den föreslagna best-practice användas vidare för berökningar av sjöegenskaper i sneda vågor. Jämförelse av resultaten med liknande studier visar god överensstämmelse. Genom att använda det föreslagna förfarandet för best-practice kan CFD-sjöegenskapssimulationer användas på fall där viskösa krafter måste beaktas, till exempel rörelseregleringsanordningar.

Cette 'Habilitation à Diriger des Recherches' a pour objectif de retracer le parcours et les activités du candidat sur les dix dernières années. Elle donne un aperçu de l'expérience acquise dans les domaines de l'enseignement, de la recherche et de l'implication dans la vie de l'établissement. La première partie récapitule le parcours académique et les activités du candidat ainsi que sa production scientifique. La deuxième partie présente une description synthétique des activités de recherche portant sur le thème principal qui traite de la modélisation analytique des couches limites turbulentes et la dispersion de particules solides en suspension. L'idée principale de cette approche consiste à proposer à l'ingénieur des outils à la fois pratiques et précis. Afin d'améliorer la description des écoulements turbulents proche parois, nous avons commencé par travailler à la résolution analytique de l'équation de l'énergie cinétique turbulente k. Nous avons proposé une solution générale qui a permis d'obtenir une fonction au voisinage immédiat de la paroi. Cette fonction a été d'abord validée par des données de simulation numérique directe (DNS) puis évaluée au niveau du code Fluent. Ces résultats ont permis de développer des modèles analytiques de viscosité turbulente qui ont permis de reconstituer avec précision le champ des vitesses en canal turbulent. A partir d'une analyse des équations de Navier-Stokes moyennées au sens de Reynolds (RANS) appliquées aux écoulements à surface libre, nous avons obtenu une équation différentielle ordinaire (EDO) pour les vitesses ainsi que sa solution semi-analytique. L'application aux canaux étroits (rapport d'aspect < 5) montre des profils de vitesses avec un maximum au-dessous de la surface libre. Nous avons développé des formulations analytiques du nombre de Schmidt turbulent nécessaire à l'évaluation du transfert de masse et des concentrations au sein des milieux fluides en mouvement avec particules. Nous avons considéré le cas de particules de sable en suspension sous l'effet d'un écoulement oscillant sur un fond de rides. Le modèle analytique de viscosité turbulente a été généralisé à ce type d'écoulements et validé par le modèle à deux équations k-ω. Cet exemple a permis en particulier de distinguer l'effet de sédimentation du phénomène de détachement tourbillonnaire. En conclusion, cette présentation sera complétée par des perspectives à court et à moyen terme.