This thesis presents a framework for simulating the fluid dynamical behavior of complex moving boundary problems in a high-performance computing environment. The framework is implemented in the pELAFINT3D software package. Moving boundaries are evolved in a seamless fashion through the use of distributed narrow band level set methods and the effect of moving boundaries is incorporated into the flow solution by a novel Cartesian grid method. The proposed Cartesian grid approach builds on the concept of a ghost fluid method where boundary conditions are applied through least-squares polynomial extrapolations. The method is hybridized such that computational cells adjacent to moving boundaries change discretization schemes smoothly in time to avoid the introduction of strong oscillations in the pressure field. The hybridization is shown to have minimal effect on accuracy while significantly suppressing pressure oscillations. The computational capability of the Cartesian grid approach is enhanced with a massively parallel adaptive meshing algorithm. Local mesh enrichment is effected through the use of octree refinement, and a scalable mesh pruning algorithm is used to reduce the memory footprint of the Cartesian grid for geometries which are not well bounded by a rectangular cuboid. The computational work is kept in a well-balanced state through the use of an adaptive repartitioning strategy. The numerical scheme is validated against many benchmark problems and the composite approach is demonstrated to work well on tens of thousands of computational cores. A simulation of the closure phase of a mechanical heart valve was carried out to demonstrate the ability of the pELAFINT3D package to compute high Reynolds number flows with complex moving boundaries and a wide disparity in length scales. Finally, a novel image-to-computation algorithm was implemented to demonstrate the flexibility the current method allows in designing new applications.

Stability against capsizing is one of the most fundamental requirements to design a ship. In this research, for the first time, CFD is performed to predict main modes of capsizing. CFD first is conducted to predict parametric rolling for a naval ship. Then CFD study of parametric rolling is extended for prediction of broaching both by using CFD as input to NDA model of broaching in replacement of EFD inputs or by using CFD for complete simulation of broaching. The CFD resistance, static heel and drift in calm water and static heel in following wave simulations are conducted to estimate inputs for NDA and 6DOF simulation in following waves are conducted for complete modeling of broaching. CFD parametric rolling simulations show remarkably close agreement with EFD. The CFD stabilized roll angle is very close to those of EFD but CFD predicts larger instability zones. The CFD and EFD results are analyzed with consideration ship theory and compared with NDA. NDA predictions are in qualitative agreement with CFD and EFD. CFD and EFD full Fr curve resistance, static heel and drift in calm water, and static heel in following waves results show fairly close agreement. CFD shows reasonable agreement for static heel and drift linear maneuvering derivatives, whereas large errors are indicated for nonlinear derivatives. The CFD and EFD results are analyzed with consideration ship theory and compared with NDA models. The surge force in following wave is also estimated from Potential Theory and compared with CFD and EFD. It is shown that CFD reproduces the decrease of the surge force near the Fr of 0.2 whereas Potential Theory fails. The CFD broaching simulations are performed for series of heading and Fr and results are compared with the predictions of NDA based on CFD, EFD, and Potential Theory inputs. CFD free model simulations show promising results predicting the instability boundary accurately. CFD calculation of wave and rudders yaw moment explains the processes of surf-riding, broaching, and periodic motion. The NDA simulation using CFD and Potential Flow inputs suggests that CFD/ Potential Flow can be considered as replacement for EFD inputs.

The Boosted Fast Flux Loop (BFFL) will expand the Advanced Test Reactor (ATR) at Idaho National Laboratory. Part of the BFFL is a new corrosion test cap section for testing in the ATR. The corrosion test cap section was designed with parallel channels to reduce the pressure drop and allow coolant contact with specimens. The fluid experiment conducted by Idaho State University found the pressure drop not characteristic of parallel channel flow but greater than without parallel channels. A Computation Fluid Dynamics simulation using STAR-CCM+ was conducted with the objectives of showing sufficient flow through the test cap section for a corrosion test, verifying the fluid experiment's validity, and explaining the abnormal pressure drop. The simulation used a polyhedral volume mesh and the k-e turbulent model with segregated equations. Convergence depended on a low continuity residual and an unchanging pressure drop result. The simulation showed the same pattern as the fluid experiment. The simulation provided evidence of flow through the test cap section needed for a corrosion test. The specimen holding assembly was found to be a small contributor to the pressure drop. The counterintuitive pressure drop was found to be the sum of many factors produced from the geometry of the test cap section. The inlet of the test cap section behaved as a diverging nozzle before a sudden expansion into the test cap section chamber with both creating a pressure drop. The chaotic flow inside the chamber gave rise to pressure loss from mixing. The fluid exited the chamber through a sudden contraction to a converging nozzle behaving exit, again, producing a pressure drop. By varying the flow rate in the simulation, a disturbance in the flow where the gap fluid separated into the parallel channels was found at high flow rates. At low flow rates the pressure drop anomaly was not found. The corrosion test cap section could be used in the ATR but with a higher pressure drop than desirable. The design of the corrosion test cap section created the abnormal pressure drop.

Natural convection is a phenomenon that occurs in a wide range of applications such as cooling towers, air conditioners, and power plants. Natural convection may be used in decay heat removal systems such as spent fuel casks, where the higher reliability inherent of natural convection is more desirable than forced convection. Passive systems, such as natural convection, may provide better safety, and hence have received much attention recently. Cooling of spent fuel rods is conventionally done using water as the coolant. However, it involves contaminating the water with radiation from the fuel rods. Contamination becomes dangerous and difficult for humans to handle. Further, the recent nuclear tragedy in Fukushima, Japan has taught us the dangers of contamination of water with nuclear radiation. Natural convection can perhaps significantly reduce the risk since it is self-sufficient and does not rely on other secondary system such as a blower as in cases of forced convection. The Utah State University Experimental Fluid Dynamics lab has recently designed an experiment that models natural convection using heated rod bundles enclosed in a rectangular cavity. The data available from this experiment provides and opportunity to study and validate computational fluid dynamics(CFD)models. The validated CFD models can be used to study multiple configurations, boundary conditions, and changes in physics(natural and/or forced convection). The results are to be validated using experimental data such as the velocity field from particle image velocimetry (PIV), pressure drops across various sections of the geometry, and temperature distributions along the vertically heated rods. This research work involves modeling natural convection using two-layer turbulence models such as k - ε and RST (Reynolds stress transport) using both shear driven (Wolfstein) and buoyancy driven (Xu) near-wall formulations. The interpolation scheme employed is second-order upwinding using the general purpose code STAR-CCM+. The pressure velocity coupling is done using the SIMPLE method. It is ascertained that turbulence models with two-layer formulations are well suited for modeling natural convection. Further it is established that k - ε and Reynolds stress turbulence models with the buoyancy driven (Xu)formulation are able to accurately predict the flow rate and temperature distribution.

A parallel numerical method has been implemented for solving the Navier Stokes equations on Cartesian and non-orthogonal meshes. To ensure the accuracy of the code first, second and third order differencing schemes, with and without flux-limiters, have been implemented and tested. The most computationally expensive task in the code is the solution of linear equations, and a number of linear solvers have been tested to determine the most efficient. Krylov space, incomplete factorisation, and other iterative and direct solvers from the literature have been implemented, and have been compared with a novel black-box multigrid linear solver that has been developed both as a solver and as a preconditioner for the Krylov space methods. To further reduce execution time the code was parallelised, after a series of experiments comparing the suitability of different parallelisation techniques and computer architectures for the Navier Stokes solver. The code has been applied to the solution of two classes of problem. Two natural convection flows were studied, with an initial study of two dimensional Rayleigh Benard convection being followed by a study of a transient three dimensional flow, in both cases the results being compared with experiment. The second class of problems modelled were free surface flows. A two dimensional free surface driven cavity, and a two dimensional flume flow were modelled, the latter being compared with analytic theory. Finally a three dimensional ship flow was modelled, with the flow about a Wigley hull being simulated for a range of Reynolds and Froude numbers.

There is a growing recognition worldwide of the need for more powerful, smaller petrol engines, capable of delivering the higher picking power of larger engines, yet still being economical and environmentally friendly when used for day-to-day driving. An engineering solution for more efficient engines has been considered by research so far. It has been identified that superchargers can potentially improve the performance of automotive engines; therefore research has focused on developing superchargers and supercharger components with higher efficiency. Of particular interest to the research presented in this thesis has been the twin-screw supercharging compressor with design adapted for automotive use (the twin-screw supercharger). The performance of this supercharger type depends on the volume and total losses of the air flow through the supercharger rotors more than on any other aspects of its behaviour. To accurately predict the efficiency of the twin-screw su percharger for matching a particular engine system, accurate supercharger design is required. The main objective of this research was the investigation of the existing limitations of twin-screw superchargers, in particular leakage and reduced efficiency, leading to the development of optimal asymmetric rotor profiles. This research has been completed in four stages defining an innovative rotor design method. The parametric three-dimensional geometric model of twin-screw supercharger rotors of any aspect ratio was developed. For model validation through visualisation, CAD rotor models with scalable data were generated in commercial CAD software and calibrated experimentally by Laser Doppler Velocimetry (LDV) tests. Calibrated rotor profile data can be transferred into CAD-CFD interface for flow simulation and performance optimisation. Through the application of this new rotor design method, new opportunities are created for the twin-screw supercharger design practice, making it a part of the engineering solution for more efficient engines.

<p>This thesis is performed at Saab Bofors Dynamics in Karlskoga and investigates a technique for ranging with laser pulses. The investigated technique is called Constant Fraction Detection (CFD). Described briefly, the received laser pulse is split into two equal parts, where one part is delayed half the pulse width and inverted. This signal is added to the original pulse. The resulting curve has the shape of a laying S and the detection of the zero level is used to stop the time measurement. The time measurement will be independent of the incoming signal’s amplitude. The CFD technique has the advantage of collecting accurate data for each send pulse, which results in an ability to collect values of measurement with a high frequency. The theses investigates a measurement frequency of 10 kHz that will give an opportunity to implement a scanning function with the possibility to, for example, reproduce a ground structure from a flying object. </p><p>The theses include both digital and analog electronics, which makes it a complex design task. The detector was constructed using analog circuits, from the signal processing of the incoming reflected pulse to the generation of a voltage level as a representation of the distance. The analog part is controlled by digital signals generated by a FPGA, which also performs calculations to convert the voltage level into a distance displayed on a LCD. </p><p>A large part of the work was dedicated to designing a layout and constructing a surface mounted printed circuit board (PCB) and therefor the report treats the whole development process, from technical requirement to construction and verification of a prototype. </p><p>The conclusion states that the CFD technique is a suitable technique for ranging with demands on fast collection of data. The prototype has sufficient accuracy at constant amplitude and was at the time of presentation shown as a prototype for demonstration. The independence of amplitude on the incoming signal was never accomplished and the reason for this is stated in the report. However, further development should solve the problem.</p>

In this study, commercially available Computational Fluid Dynamics (CFD) software, CFX-4.4 has been used for the simulations of aerosol transport through various aerosol-sampling devices. Aerosol transport was modeled as a classical dilute and dispersed two-phase flow problem. Eulerian-Lagrangian framework was adopted wherein the fluid was treated as the continuous phase and aerosol as the dispersed phase, with a one-way coupling between the phases. Initially, performance of the particle transport algorithm implemented in the code was validated against available experimental and numerical data in the literature. Code predictions were found to be in good agreement against experimental data and previous numerical predictions. As a next step, the code was used as a tool to optimize the performance of a virtual impactor prototype. Suggestions on critical geometrical details available in the literature, for a virtual impactor, were numerically investigated on the prototype and the optimum set of parameters was determined. Performance curves were generated for the optimized design at various operating conditions. A computational model of the Linear Slot Virtual Impactor (LSVI) fabricated based on the optimization study, was constructed using the worst-case values of the measured geometrical parameters, with offsets in the horizontal and vertical planes. Simulations were performed on this model for the LSVI operating conditions. Behavior of various sized particles inside the impactor was illustrated with the corresponding particle tracks. Fair agreement was obtained between code predictions and experimental results. Important information on the virtual impactor performance, not known earlier, or, not reported in the literature in the past, obtained from this study, is presented. In the final part of this study, simulations on aerosol deposition in turbulent pipe flow were performed. Code predictions were found to be completely uncorrelated to experimental data. The discrepancy was traced to the performance of the code's turbulent dispersion model. A detailed literature survey revealed the inherent technical deficiencies in the model, even for particle dispersion. Based on the results of this study, it was determined that while the code can be used for simulating aerosol transport under laminar flow conditions, it is not capable of simulating aerosol transport under turbulent flow conditions.

Cette étude porte sur l'utilisation conjointe (appelée ” couplage ”) de modèle de turbulence à grandes échelles LES (Large Eddy Simulation) et du modèle multiphasique VOF (Volume of Fluid). Cette utilisation conjointe est nécessaire dans de nombreuses applications industrielles comme celles de l'automobile où l'on recherche par exemple à évaluer les prestations diphasiques de l'auvent liées au phénomène d'entraînement et d'ingestion d'une nappe d'eau par le HVAC (système d'air conditionné). Cependant, l'utilisation conjointe de ces méthodes nécessite un traitement particulier de la turbulence proche de la surface liquide afin de reproduire convenablement la quantité de mouvement transmise depuis la phase gazeuse, motrice, jusque dans la phase liquide. Basée sur une approche numérique et expérimentale, notre étude est articulée autour de trois axes. Tout d'abord, le développement de techniques de mesures spécifiques pour l'étude expérimentale de notre problématique: le LeDaR pour mesurer les déformées d'une interface et la PIV d'interface afin d'accéder aux champs de vitesse et de turbulence dans chacune des deux phases. Le second axe est la constitution d'une base de données expérimentales sur une configuration de type jet impactant sur une surface liquide représentative des phénomènes rencontrés dans l'auvent. Enfin, le troisième axe de travail est l'évaluation des modèles existants dans le code Ansys Fluent et à partir de cette analyse le développement et la validation de modèles de couplage LES-VOF. L'évaluation des modèles développés a permis de valider une stratégie de calcul adaptée aux simulations de l'entraînement d'une surface d'eau par un écoulement d'air turbulent.

This research presents the results of numerical study of heat transfer and pressure drop in a heat exchanger that is designed with different shape pin fins. The heat exchanger used for this research consists of a rectangular duct fitted with different shape pin fins, and is heated from the lower plate. The pin shape and the compact heat exchanger (CHE) configuration were numerically studied to maximize the heat transfer and minimize the pressure drop across the heat exchanger. A three dimensional finite volume based numerical model using FLUENT© was used to analyze the heat transfer characteristics of various pin fin heat exchangers. The simulation applied to estimate the heat transfer coefficient and pressure drop for a wide range of Reynolds numbers with different pin fins. Circular pin configuration variations included changes in pin spacing, axial pitch and pin height ratio. Rectangular and drop-shaped pin variations also included changes in length and aspect ratio. Correlations for Nusselt number and friction factor were developed. The optimum drop shaped pin array was shown to match the heat transfer rates obtained by the optimum circular pin configuration while incurring less than one third the specific fluid friction power losses. The data and conclusions of this study can be applied to the optimization of different heat exchangers which are used in industry, especially oil cooler in power transformers which are currently working with low cooling efficiency. It can also be used in the design of electronic components, turbine blade cooling or in other high heat flux dissipation applications requiring a low-profile, high area-density based micro-heat exchanger design. This study also shows that numerical models backed with experimental analysis can reduce both the time and money required to create and evaluate engineering concepts, especially those that deal with fluid flow and heat transfer. In the following chapters, first the problems which are encountered by power transformer suppliers are described. Then pin fin technology is studied with more details as a novel solution to the oil cooling problem. Some studies on behavior of power transformer coolers are also conducted to make their problems more clear. Available experimental data in the Iran Transfo company have been used for validation of these studies. They are presented as separated papers at the end of thesis. Finally the results of pin fin studies are presented and horizontal continuous casting (HCC) is explained as a manufacturing method for pin fins production. A separate paper which is based on experimental study on HCC is also included at the end of thesis. / Forskningen presenterad är ett resultat av en numerisk studie av värmeöverföring och tryckfall i en värmeväxlare designad med olika former av Kylflänsar. Värmeväxlaren består av ett rektangulärt kanal utrustat med olika former av Kylflänsar och är uppvärmd underifrån. Kylflänsar forma och den kompakta värmeväxlarens utformning är studerade numeriskt för att maximera värmeöverföringen och minimera tryckfallet över värmeväxlaren.En tredimensionell finit volym baserad på en numerisk modell i FLUENT© användes för att analysera värmeöverföringsegenskaper för olika Kylflänsar konfigurationer. Genom simuleringar uppskattades värmegenomgångstalet och tryckfallet för olika Reynolds tal och Kylflänsar konfigureringar. Cirkulära Kylflänsar konfigurationer inkluderar variation av avstånd mellan Kylflänsar, och förhållandet mellan axiellt avstånd och höjd. Rektangulära och droppformade Kylflänsar inkluderade även variation för längd och aspekt förhållande. Korrelation mellan Nusselts tal och friktionsfaktor utvecklades. Optimal matris för hur droppformade Kylflänsar placerades visades överensstämma med optimal överföring för cirkulära Kylflänsar men bara med en tredjedel av friktionsförlusterna för fluiden. Data och slutsatser från studien kan användas inom för optimering av värmeväxlare använda i industrin, speciellt oljekylda högspänningstransformatorer som har låg effektivitet i kylningen. Resultaten kan även användas inom design av elektronikkomponenter, kylning av turbinblad eller andra komponenter med högt värmeflöde där låg profil, och stor ytdensitet behövs. Studien visar att kombinationen av numeriska modeller som valideras genom experiment kan reducera både tid och kostnad vid utveckling och utvärdering av ingenjörsverktyg, speciellt inom fluidmekanik och värmeöverföring. I följande kapitel beskrivs först problem som identifierats av tillverkare av högspänningstransformatorer. Kylflänsar studeras i detalj som en ny lösning till de identifierade problemen med oljekylning. Några studier har genomförts för att ytterligare belysa problemen kring högspänningstransformatorers kylning. Tillgängliga data från Iran Transfo company har använts för validering av resultat från studierna. Studierna presenteras som separata artiklar i slutet av avhandlingen. Avslutningsvis presenteras resultat från studierna av Kylflänsar och en horisontell kontinuerlig gjutprocess (HCC) för tillverkning av Kylflänsar. HCC-studien presenteras som en separat artikel inkluderad sist i avhandlingen.