This project solves the problem of cooling the photovoltaic cell. Solar cell was modeled according to a real model in SolidWorks, and subsequently created the model was simulated in SolidWorks Flow Simulation and Ansys Fluent. The use of both systems allow a comparison of their possibilities in the field of heat transfer and their suitability for the case. The conclusion summarizes the first results and outline further developments cooling design to optimize the performance of the solar cell.

In our study, the Finite Volume Method (FVM) is successfully implemented to simulate thermal process of polymerization. This application is verified based on the obtained plots compared with those from other two methods as well as experimental data. After the verification, a method is developed to optimize heat history in order to reduce processing time and in the meantime to maintain the uniformity of cure state. Also sensitivities of cure state to different parameters are examined. Besides, a correlation between temperature and the degree of polymerization profile on sample surface is found using on-line monitoring method.

Consider a bottom-hinged Oscillating Wave Surge Converter (OWSC): This device oscillates due to the hydrodynamic forces applied on it by the action of ocean waves. The focus of this thesis is to build upon the in-house multi-block generalized coordinate finite volume solver GenIDLEST using a collocated grid arrangement within the framework of the fractional-step method to make it compatible to simulate such systems. The first step in this process is to deploy a convection scheme which differentiates between air and water. This process is further complicated by the 1:1000 density and 1:100 viscosity ratio between the two fluids. For this purpose, a phase field method is chosen for its ease of implementation and proven boundedness and conservativeness properties. Extensive validation and verification using standard test cases, such as droplet in shear flow, Rayleigh Taylor instability, and the Dam Break Problem is carried out. This development is then coupled with the present Immersed Boundary Module which is used to simulate the presence of moving bodies and again verified against test cases, such as the Dam Break problem with a vertical obstacle and heave decay of a partially submerged buoyant cylinder. Finally, a relaxation zone technique is used to generate waves and a numerical beach technique is used to absorb them. These are then used to simulate the Oscillating Surge Wave Converter. / Master of Science / An Oscillating Wave Surge Converter can be best described as a rectangular flap, hinged at the bottom, rotating under the influence of ocean waves from which energy is harvested. The singular aim of this thesis is to model this device using Computational Fluid Dynamics (CFD). More specifically, the aim is to model this dynamic device with the full Navier Stokes Equations, which include inertial forces, arising due to the motion of the fluid, viscous forces which dissipate energy, and body forces such as gravity. This involves three key steps: 1. Modelling the air-water interface using a convection scheme. A phase field method is used to differentiate between the two fluids. This task is made more challenging because of the very large density and viscosity differences between air and water. 2. Model dynamic moving geometries in a time-dependent framework. For this, we rely on the Immersed Boundary Method. 3. Develop a numerical apparatus to generate and absorb ocean waves. For this, we rely on the Relaxation Zone and Numerical Beach Method. These developments are validated in different canonical problems and finally applied to a two-dimensional oscillating surge wave energy converter.

No description available.

No description available.

Mathematical modeling of ocean waves is based on the formulation and solution of the appropriate equations of continuity, momentum and the choice of proper initial and boundary conditions. Under the influence of gravity, many free surface water waves can be modeled by the shallow water equations (SWE) with the assumption that the horizontal length scale of the wave is much greater than the depth scale and the wave height is much less than the fluid's mean depth. Furthermore, to describe three dimensional flows in the hydrostatic and Boussinesq limits, the multilayer SWE model is used, where the fluid is discretized horizontally into a set of vertical layers, each having its own height, density, horizontal velocity and geopotential. In this study, we used an explicit staggered finite volume method to solve single and multilayer SWE, with and without density stratification and bathymetry, to understand the dynamic of surface waves and internal waves. We implemented a two-dimensional version of the incompressible DYNAMICO method and compare it with a one-dimensional SWE. For multilayer SWE, we considered both two layer and a linear stratification of density, with very small density gradient, consistent with Boussinesq approximation. We used Lagrangian vertical coordinate which doesn't allow mass to flow across vertical layers. Numerical examples are presented to verify multilayer SWE model against single layer SWE, in terms of the phase speed and the steepness criteria of wave profile. In addition, the phase speed of the barotropic and baroclinic mode of two-layer SWE also verified our multilayer SWE model. We found that, for multilayer SWE, waves move slower than single layer SWE and get steeper than normal when they flow across bathymetry. A series of numerical experiment were carried out to compare 1-D shallow water solutions to 2-D solutions with and without density as well as to explain the dynamics of surface wave and internal wave. We found that, a positive fluctuations on free surface causes water to rise above surface level, gravity pulls it back and the forces that acquired during the falling movement causes the water to penetrate beneath it's equilibrium level, influences the generation of internal waves. Internal waves travel considerably more slowly than surface waves. On the other hand, a bumpy or a slicky formation of surface waves is associated with the propagation of internal waves. The interaction between these two waves is therefore demonstrated and discussed. / Thesis / Master of Science (MSc) / In the modelling of ocean wave, the formulation and solution of appropriate equations and proper initial and boundary conditions are required. The shallow water equations (SWE) are derived from the conservation of mass and momentum equations, in the case where the horizontal length scale of the wave is much greater than the depth scale and the wave height is much less than the fluid's mean depth. In multilayer SWE, the fluid is discretized horizontally into a set of vertical layers, each having its own height, density, horizontal velocity and geopotential. In this study, we used an explicit staggered finite volume method to solve single and multilayer SWE, with and without density stratification and bathymetry, to understand the dynamic of surface waves and internal waves. A series of numerical experiments were carried out to validate our multilayer model. It is found that, in the presence of density differences, surface waves for the multilayer SWE move slowly and get more steep than normal when they flow across bathymetry. Also, a positive fluctuations on free surface generates internal waves at the interior of ocean which propagate along the line of density gradient.

In this dissertation we discuss bilinear immersed finite elements (IFE) for solving interface problems. The related research works can be categorized into three aspects: (1) the construction of the bilinear immersed finite element spaces; (2) numerical methods based on these IFE spaces for solving interface problems; and (3) the corresponding error analysis. All of these together form a solid foundation for the bilinear IFEs. The research on immersed finite elements is motivated by many real world applications, in which a simulation domain is often formed by several materials separated from each other by curves or surfaces while a mesh independent of interface instead of a body-fitting mesh is preferred. The bilinear IFE spaces are nonconforming finite element spaces and the mesh can be independent of interface. The error estimates for the interpolation of a Sobolev function in a bilinear IFE space indicate that this space has the usual approximation capability expected from bilinear polynomials, which is <i>O</i>(<i>h</i>²) in <i>L</i>² norm and <i>O</i>(<i>h</i>) in <i>H</i>¹ norm. Then the immersed spaces are applied in Galerkin, finite volume element (FVE) and discontinuous Galerkin (DG) methods for solving interface problems. Numerical examples show that these methods based on the bilinear IFE spaces have the same optimal convergence rates as those based on the standard bilinear finite element for solutions with certain smoothness. For the symmetric selective immersed discontinuous Galerkin method based on bilinear IFE, we have established its optimal convergence rate. For the Galerkin method based on bilinear IFE, we have also established its convergence. One of the important advantages of the discontinuous Galerkin method is its flexibility for both <i>p</i> and <i>h</i> mesh refinement. Because IFEs can use a mesh independent of interface, such as a structured mesh, the combination of a DG method and IFEs allows a flexible adaptive mesh independent of interface to be used for solving interface problems. That is, a mesh independent of interface can be refined wherever needed, such as around the interface and the singular source. We also develop an efficient selective immersed discontinuous Galerkin method. It uses the sophisticated discontinuous Galerkin formulation only around the locations needed, but uses the simpler Galerkin formulation everywhere else. This selective formulation leads to an algebraic system with far less unknowns than the immersed DG method without scarifying the accuracy; hence it is far more efficient than the conventional discontinuous Galerkin formulations. / Ph. D.

Numerical simulations of fluid-structure interaction (FSI) problems were performed in an effort to verify and validate a commercially available FSI tool. This tool uses an iterative partitioned coupling scheme between CD-adapco's STAR-CCM+ finite volume fluid solver and Simulia's Abaqus finite element structural solver to simulate the FSI response of a system. Preliminary verification and validation work (VandV) was carried out to understand the numerical behavior of the codes individually and together as a FSI tool. Verification and Validation work that was completed included code order verification of the respective fluid and structural solvers with Couette-Pouiselle flow and Euler-Bernoulli beam theory. These results confirmed the 2nd order accuracy of the spatial discretizations used. Following that, a mixture of solution verifications and model calibrations was performed with the inclusion of the physics models implemented in the solution of the FSI problems. Solution verifications were completed for fluid and structural stand-alone models as well as for the coupled FSI solutions. These results re-confirmed the spatial order of accuracy but for more complex flows and physics models as well as the order of accuracy of the temporal discretizations. In lieu of a good material definition, model calibration is performed to reproduce the experimental results. This work used model calibration for both instances of hyperelastic materials which were presented in the literature as validation cases because these materials were defined as linear elastic. Calibrated, three dimensional models of the bow seal on the University of Michigan bow seal test platform showed the ability to reproduce the experimental results qualitatively through averaging of the forces and seal displacements. These simulations represent the only current 3D results for this case. One significant result of this study is the ability to visualize the flow around the seal and to directly measure the seal resistances at varying cushion pressures, seal immersions, forward speeds, and different seal materials. SES design analysis could greatly benefit from the inclusion of flexible seals in simulations, and this work is a positive step in that direction. In future work, the inclusion of more complex seal geometries and contact will further enhance the capability of this tool. / Ph. D.

The solution to partial differential equations generally requires approximations that result in numerical error in the final solution. Of the different types of numerical error in a solution, discretization error is the largest and most difficult error to estimate. In addition, the accuracy of the discretization error estimates relies on the solution (or multiple solutions used in the estimate) being in the asymptotic range. The asymptotic range is used to describe the convergence of a solution, where an asymptotic solution approaches the exact solution at a rate proportional to the change in mesh spacing to an exponent equal to the formal order of accuracy. A non-asymptotic solution can result in unpredictable convergence rates introducing uncertainty in discretization error estimates. To account for the additional uncertainty, various discretization uncertainty estimators have been developed. The goal of this work is to evaluation discretization error and discretization uncertainty estimators based on Richardson extrapolation for computational fluid dynamics problems. In order to evaluate the estimators, the exact solution should be known. A select set of solutions to the 2D Euler equations with known exact solutions are used to evaluate the estimators. Since exact solutions are only available for trivial cases, two applications are also used to evaluate the estimators which are solutions to the Navier-Stokes equations: a laminar flat plate and a turbulent flat plate using the k-Ï SST turbulence model. Since the exact solutions to the Navier-Stokes equations for these cases are unknown, numerical benchmarks are created which are solutions on significantly finer meshes than the solutions used to estimate the discretization error and uncertainty. Metrics are developed to evaluate the accuracy of the error and uncertainty estimates and to study the behavior of each estimator when the solutions are in, near, and far from the asymptotic range. Based on the results, general recommendations are made for the implementation of the error and uncertainty estimators. In addition, a new uncertainty estimator is proposed with the goal of combining the favorable attributes of the discretization error and uncertainty estimators evaluated. The new estimator is evaluated using numerical solutions which were not used for development and shows improved accuracy over the evaluated estimators. / Master of Science

In Diesel engines, the internal flow characteristics in the fuel injection nozzles, such as the turbulence level and distribution, the cavitation pattern and the velocity profile affect significantly the air-fuel mixture in the spray and subsequently the combustion process. Since the possibility to observe experimentally and measure the flow inside real size Diesel injectors is very limited, Computational Fluid Dynamics (CFD) calculations are generally used to obtain the relevant information. The work presented within this thesis is focused on the study of cavitation in real size automotive injectors by using a commercial CFD code. It is divided in three major phases, each corresponding to a different complementary objective. The first objective of the current work is to assess the ability of the cavitation model included in the CFD code to predict cavitating flow conditions. For this, the model is validated for an injector-like study case defined in the literature, and for which experimental data is available in different operating conditions, before and after the start of cavitation. Preliminary studies are performed to analyze the effects on the solution obtained of various numerical parameters of the cavitation model itself and of the solver, and to determine the adequate setup of the model. It may be concluded that overall the cavitation model is able to predict the onset and development of cavitation accurately. Indeed, there is satisfactory agreement between the experimental data of injection rate and choked flow conditions and the corresponding numerical solution.This study serves as the basis for the physical and numerical understanding of the problem. Next, using the model configuration obtained from the previous study, unsteady flow calculations are performed for real-size single and multi-hole sac type Diesel injectors, each one with two types of nozzles, tapered and cylindrical. The objective is to validate the model with real automotive cases and to ununderstand in what way some physical factors, such as geometry, operating conditions and needle position affect the inception of cavitation and its development in the nozzle holes. These calculations are made at full needle lift and for various values of injection pressure and back-pressure. The results obtained for injection rate, momentum flux and effective injection velocity at the exit of the nozzles are compared with available CMT-Motores Térmicos in-house experimental data. Also, the cavitation pattern inside the nozzle and its effect on the internal nozzle flow is analyzed. The model predicts with reasonable accuracy the effects of geometry and operating conditions. / Patouna, S. (2012). A CFD STUDY OF CAVITATION IN REAL SIZE DIESEL INJECTORS [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/14723