The present study aims to explore flow behavior in the mold during centrifugal casting process by numerical simulations. The theoretical model comprises two groups of steady conservation equations of mass and momentum and the governing equations are solved numerically with k-£` turbulence model and iterative SIMPLE(Semi-Implicit for Pressure-Linked Equations) algorithm to determine the flow property. The numerical results indicate that the melt liquids of titanium compound metal flow near the walls in the high rotation rate. With the high rotation rate, the outflow velocity is rapid. The products depend on the flow of the melt liquids. The flow is rapid, the filled process is quick, and the temperature is uniform. But if the flow is slow, the temperature is not uniform in the filled process. Therefore, the products may result in faults. When the rotation rate is up to 50 rpm, the melt liquids flow near the walls with the affect of centrifugal force. The velocity is larger than the velocity with zero rotation rate . Thus the flow with rotation rate can help to fill quickly, and reduce the temperature loss.

A FORTRAN code was developed to numerically simulate the mixed convective flow over a three-dimensional horizontal backward-facing step. The momentum and energy equations under the assumption of the Boussinesq approximation were discretized by means of a finite volume technique. The SIMPLE algorithm scheme was applied to link the pressure and velocity fields inside the domain while an OpenMP parallel implementation was proposed to improve the numerical performance and to accelerate the numerical solution. The heating process corresponds to a channel heated from below at constant temperature keeping insulated all the other channel walls. In addition, the back-step was considered as a thermally conducting block and its influence in the heating process was explored by holding different solid to fluid thermal conductivity ratios. The effects over the velocity and temperature distribution of buoyancy forces, acting perpendicular to the mainstream flow, are studied for three different Richardson numbers Ri=3, 2, and 1 and the results are compared against those of pure forced convection Ri=0. In these simulations the Reynolds number is fixed at 200 while the bottom wall temperature is adjusted to fulfill the conditions for the different Ri. Under this assumption, as Ri increases the buoyancy effects are the dominant effects in the mixed convective process. The numerical results indicate that the velocity field and the temperature distribution for pure forced convection are highly distorted if compared with the mixed convective flow. If the Ri parameter is increased, then the primary re-circulation zone is reduced. Similarly, as the buoyancy forces become predominant in the flow, the convective rolls, in the form of spiral-flow structures, become curlier and then higher velocity components are found inside the domain. The temperature field distribution showed that as the Ri is increased a thicker layer of high temperature flow is located at the channel??s top wall as a result of the higher rates of low-density flow moving to the top wall. The flow is ascending by the channel sidewalls, while descending by the channel span-wise central plane. The parallel numerical strategy is presented and some results for the performance of the OpenMP implementation are included. In this sense, linear speedup was obtained when using 16 possessors in parallel.

In this work, numerical simulations of flow separation control are performed. The sep-aration control technique studied is called 'synthetic jet actuation'. The developed code employs a cell centered finite volume scheme which handles viscous, steady and unsteady compressible turbulent flows. The pulsating zero mass jet flow is simulated by imposing a harmonically varying transpiration boundary condition on the airfoil's surface. Turbulence is modeled with the algebraic model of Baldwin and Lomax. The application of synthetic jet actuators is based in their ability to energize the boundary layer, thereby providing signifcant increase in the lift coefficient. This has been corroborated experimentally and it is corroborated numerically in this research. The performed numerical simulation investigates the flow over a NACA0015 air-foil. For this flow Re = 9??105 and the reduced frequency and momentum coefficient are F+ = 1:1 and C?? = 0:04 respectively. The oscillatory injection takes place at 12.27% chord from the leading edge. A maximum increase in the mean lift coefficient of 93% is predicted by the code. A discrepancy of approximately 10% is observed with corresponding experimental data from the literature. The general trend is, how-ever, well captured. The discrepancy is attributed to the modeling of the injection boundary condition and to the turbulence model.A sensitivity analysis of the lift coefficient to different values of the oscillation parameters is performed. It is concluded that tangential injection, F + ?? O(1) and the utilized grid resolution around the site of injection are optimal. Streamline fields ob-tained for different angles of injection are analyzed. Flow separation and attachment as functions of the injection angle and of the velocity of injection can be observed. It is finally concluded that a reliable numerical tool has been developed which can be utilized as a support tool in the optimization of the synthetic jet operation and in the modeling of its operation.

Proton exchange membrane (PEM) fuel cell possesses the characteristics of microminiaturization and low temperature operation. For this reason, the proton exchange membrane fuel cell is very suitable to serve as power source of portable electronic products. In this paper, a three-dimensional numerical model to evaluate the voltage and the total current density of a PEM fuel cell stack was developed. The polarization curves of the PEM fuel cell stack under three different operating temperatures were investigated. In this study, the micro PEM fuel cell stack contains two single cells. Pure H2 gas stream was supplied as the anode inlet flow and air as the cathode inlet flow under constant pressure at 97 kPa and constant cell temperate (298K¡B308K¡B323K) conditions. Because the cell temperature may affect the chemical reaction rate on the cathode side, we discussed the influences of different temperatures on the cell performance. Solutions were compared with the experimental data. Both the value of power density and the tendency of polarization curve are in good agreement with the experimental data.

Enceladus, a tiny moon of Saturn, is found to be geologically active. In 2005, Cassini detected an anomalously warm region and a plume, consisting of mostly water vapor and ice grains, at its south pole. The plume has far-reaching effects on the Saturnian system and offers clues into the moon’s interior, particularly as to whether liquid water exists underground. Consequently, understanding the physics and source conditions of the plume is crucial, which is the focus of this work. The plume is not only two-phase but also multi-regime in nature and can be divided into three distinct regions: a subsurface region, a collisional near-field and a free-molecular far-field. To study it, a hybrid model of the plume, which treats each region separately, is constructed. Two subsurface models are considered. Using the resulting vent conditions from these models, the plume is propagated from the surface vents out to several Enceladus radii using the direct simulation Monte Carlo (DSMC) method in the near-field and a free-molecular model in the far-field. The model is used to examine the plume flow, with and without grains. Collisions are found to be important in various processes, including grain condensation and flow acceleration. Since collision rates drop away from the vent, they must be high enough at the vent to enable significant condensation to occur and the gas to accelerate to the maximum speed possible by allowing energy stored in internal modes to be converted into translational energy as the gas expands. When slower grains are present, however, they may decelerate and push the gas out more laterally. Moreover, grains may form a thick column and restrict the free expansion of the gas. Smaller grains have greater and more extensive effects on the gas, but are also more strongly affected by the gas. Their motions decouple from the gas motions higher above the vent. They are also more likely to spread with the gas and be accelerated to the gas speeds. By constraining the plume far-field using Cassini data, the H2O and grain production rates from the plume are estimated to be ~100–1000 kg/s and < 10 kg/s respectively, which agree with other estimates. Based on fit results, the gas jets appear to be narrow, suggesting high Mach numbers at the vents. / text

Controlled cooling is an important part of steel production industry that affects the properties of the outcome steel. Many of the researches done in controlled cooling are experimental. Due to progress in the numerical techniques and high cost of experimental works in this field the numerical work seems more feasible. Heat transfer analysis is the necessary element of successful controlled cooling and ultimately achievement of novel properties in steel. Heat transfer on the surface of the plate normally contains different regimes such as film boiling, nucleate boiling, transition boiling and radiation heat transfer. This makes the analysis more complicated. In order to perform the heat transfer analysis often empirical correlations are being used. In these correlations the velocity and pressure within the fluid domain is involved. Therefore in order to obtain a better understanding of heat transfer process, study of hydrodynamics of the fluid becomes necessary. Circular jet due to its high efficiency has been used vastly in the industry. Although some experimental studies of round jet arrays have been done, yet the characteristics of a single jet with industrial geometric and flow parameters on the surface of a flat plate is not fully understood. Study of hydrodynamics of the jet impingement is the first step to achieve better understanding of heat transfer process. Finite element method as a popular numerical method has been used vastly to simulate different domains. Traditional approaches of finite element method, Lagrangian and Eulerian, each has its own benefits and drawbacks. Lagrangian approach has been used widely in solid domains and Eulerian approach has been widely used in fluid fields. Jet impingement problem, due to its unknown free surface and the change in the boundary, falls in the category of special problems and none of the traditional approaches is suitable for this application. The Arbitrary Lagrangian Eulerian (ALE) formulation has emerged as a technique that can alleviate many of the shortcomings of the traditional Lagrangian and Eulerian formulations in handling these types of problems. Using the ALE formulation the computational grid need not adhere to the material (Lagrangian) nor be fixed in space (Eulerian) but can be moved arbitrarily. Two distinct techniques are being used to implement the ALE formulation, namely the operator split approach and the fully coupled approach. This thesis presents a fully coupled ALE formulation for the simulation of flow field. ALE form of Navier-Stokes equations are derived from the basic principles of continuum mechanics and conservation laws in the fluid. These formulations are then converted in to ALE finite element equations for the fluid flow. The axi-symmetric form of these equations are then derived in order to be used for jet impingement application. In the ALE Formulation as the mesh or the computational grid can move independent of the material and space, an additional set of unknowns representing mesh movement appears in the equations. Prescribing a mesh motion scheme in order to define these unknowns is problem-dependent and has not been yet generalized for all applications. After investigating different methods, the Winslow method is chosen for jet impingement application. This method is based on adding a specific set of partial differential Equations(Laplace equations) to the existing equations in order to obtain enough equations for the unknowns. Then these set of PDEs are converted to finite element equations and derived in axi-symmetric form to be used in jet impingement application. These equations together with the field equations are then applied to jet impingement problem. Due to the number of equations and nonlinearity of the field equations the solution of the problem faces some challenges in terms of convergence characteristics and modeling strategies. Some suggestions are made to deal with these challenges and convergence problems. Finally the numerical treatment and results of analyzing hydrodynamics of the Jet Impingement is presented. The work in this thesis is confined to the numerical simulation of the jet impingement and the specifications of an industrial test setup only have been used in order to obtain the parameters of the numerical model.

In this work, the lattice Boltzmann method (LBM) was validated for direct numerical simulation (DNS) of wall-bounded turbulent flows. The LBM is a discrete-particle-based method that numerically solves the Boltzmann equation as opposed to conventional DNS methods that are based on the Navier-Stokes (NS) equations. The advantages of the LBM are its simple implementation, its ability to handle complex geometries, and its scalability on modern high-performance computers. An LBM code was developed and used to simulate fully-developed turbulent channel flow. In order to validate the results, the turbulence statistics were compared to those calculated from a conventional NS-based finite difference (FD) simulation. In the present study, special care was taken to make sure the computational domains for LBM and FD simulations were the same. Similar validation studies in the literature have used LBM simulations with smaller computational domains in order to reduce the computational cost. However, reducing the size of the computational domain affects the turbulence statistics and confounds the results of the validation. The turbulence statistics calculated from the LBM and FD simulations were found to agree qualitatively; however, there were several significant deviations, particularly in the variance profiles. The largest discrepancy was in the variance of the pressure fluctuations, which differed by approximately 7%. Given that both the LBM and FD simulations resolved the full range of turbulent scales and no models were used, this error was deemed to be significant. The cause of the discrepancy in the pressure variance was found to be the compressibility of the LBM. The LBM allows the density to vary, while the FD method does not since it solves the incompressible form of the NS equations. The effect of the compressibility could be reduced by lowering the Mach number, but this would come at the cost of significantly increasing the computational cost. Therefore, the conclusion of this work is that, while the LBM is capable of producing accurate solutions for incompressible turbulent flows, it is significantly more expensive than conventional methods for simple wall-bounded turbulent flows. / Thesis (Ph.D, Mechanical and Materials Engineering) -- Queen's University, 2011-09-15 23:24:09.968

No description available.

This thesis is devoted to mathematical modelling, including numerical simulation, in nonlinear laser physics. We study the concrete problem of finding optimal parameters for soliton lasers which is important for many applications. The basic mathematical tool used in our studies is the theory of nonlinear differential equations and numerical methods of their investigations.

abstract: Numerical simulations are very helpful in understanding the physics of the formation of structure and galaxies. However, it is sometimes difficult to interpret model data with respect to observations, partly due to the difficulties and background noise inherent to observation. The goal, here, is to attempt to bridge this gap between simulation and observation by rendering the model output in image format which is then processed by tools commonly used in observational astronomy. Images are synthesized in various filters by folding the output of cosmological simulations of gasdynamics with star-formation and dark matter with the Bruzual- Charlot stellar population synthesis models. A variation of the Virgo-Gadget numerical simulation code is used with the hybrid gas and stellar formation models of Springel and Hernquist (2003). Outputs taken at various redshifts are stacked to create a synthetic view of the simulated star clusters. Source Extractor (SExtractor) is used to find groupings of stellar populations which are considered as galaxies or galaxy building blocks and photometry used to estimate the rest frame luminosities and distribution functions. With further refinements, this is expected to provide support for missions such as JWST, as well as to probe what additional physics are needed to model the data. The results show good agreement in many respects with observed properties of the galaxy luminosity function (LF) over a wide range of high redshifts. In particular, the slope (alpha) when fitted to the standard Schechter function shows excellent agreement both in value and evolution with redshift, when compared with observation. Discrepancies of other properties with observation are seen to be a result of limitations of the simulation and additional feedback mechanisms which are needed. / Dissertation/Thesis / Ph.D. Physics 2012