The practice of building simulation is split between domains such as energy, multizone airflow, computational fluid dynamics (CFD) airflow, and controls analysis, as well as between the tools which conduct these analyses. Previous work in the integration of these analyses and tools have focused on linking existing tools, written in algorithmic programming languages, together by interfacing them using coupling mechanisms implemented in algorithmic programming languages. This thesis takes a different approach, using the equation-based, object oriented modeling language Modelica to create models in different domains and interfaces between those models within a single framework which has benefits to the modeler/analyst in terms of both representation of physical processes and flexibility in modeling systems composed of many interacting components. Specifically, the simulation of airflows within buildings has historically been compartmentalized into distinct domains such as nodal network (multizone) simulations and CFD. Such airflow simulations are also often treated independently of building energy simulations (via heat transfer) despite their interrelation. Recent work has reported on combining these types of analyses by linking pre-existing simulation software together. Here a prototype CFD package of models is built in Modelica and coupled to models of conductive heat transfer and controls. Comparisons of results of simulations so constituted to analytical solutions and benchmark data available in the literature show good agreement, indicating the technical viability of this approach. Limitations include the absence of turbulence modeling and the lack of modeling features which improve computational efficiency, such as non-uniform grids.

Accurate computational representations of highly deformable surfaces are indispensable in the fields of computer animation, medical simulation, computer vision, digital modeling, and computational physics. The focus of this dissertation is on the animation of physics-based phenomena with highly detailed deformable surfaces represented by triangle meshes. We first present results from an algorithm that generates continuum mechanics animations with intricate surface features. This method combines a finite element method with a tetrahedral mesh generator and a high resolution surface mesh, and it is orders of magnitude more efficient than previous approaches. Next, we present an efficient solution for the challenging problem of computing topological changes in detailed dynamic surface meshes. We then introduce a new physics-inspired surface tracking algorithm that is capable of preserving arbitrarily thin features and reproducing realistic fine-scale topological changes like Rayleigh-Plateau instabilities. This physics-inspired surface tracking technique also opens the door for a unique coupling between surficial finite element methods and volumetric finite difference methods, in order to simulate liquid surface tension phenomena more efficiently than any previous method. Due to its dramatic increase in computational resolution and efficiency, this method yielded the first computer simulations of a fully developed crown splash with droplet pinch off.

Horizontal-axis wind turbines operate in a complex, inherently unsteady aerodynamic environment. The flow over the blades is dominated by 3-D effects, particularly during stall, which is accompanied by massive flow separation and vortex shedding. There is always bluff-body shedding from the turbine nacelle and support structure which interacts with the rotor wake. In addition, the high aspect ratios of wind turbine blades make them very flexible, leading to substantial aeroelastic deformation of the blades, altering the aerodynamics. Finally, when situated in a wind farm, turbines must operate in the unsteady wake of upstream neighbors. Though computational fluid dynamics (CFD) has made significant inroads as a research tool, simple, inexpensive methods, such as blade element momentum theory, are still the workhorses in wind turbine design and aeroelasticity applications. These methods are unable to accurately predict rotor loads near the edges of the operating envelope. In this work, a range of unstructured grid CFD techniques for predicting wind turbine loads and aeroelasticity has been developed and applied to the NREL Unsteady Aerodynamics Experiment Phase VI rotor. First, a kd-tree based nearest neighbor search algorithm was used to improve the computational efficiency of an approximate unsteady actuator blade method. This method was then shown to predict root and tip vortex locations and strengths similar to an overset method, but without the computational expense of modeling the blade surfaces. A hybrid Reynolds-averaged Navier-Stokes / Large Eddy Simulation (HRLES) turbulence model was extended to an unstructured grid framework and demonstrated to improve predictions of unsteady loading and shedding frequency in massively separated cases. For aeroelastic predictions, a methodology for tight coupling between an unstructured CFD solver and a computational structural dynamics tool was developed. Finally, time-accurate overset rotor simulations of a complete turbine---blades, nacelle, and tower---were conducted using both RANS and HRLES turbulence models. The HRLES model was able to accurately predict rotor loads when stalled. In yawed flow, excellent correlations of mean blade loads with experimental data were obtained across the span, and wake asymmetry and unsteadiness were also well-predicted.

A comprehensive understanding of the hydrodynamics in vegetated open-channels and flow-vegetation interaction is of high interest to researchers and practitioners alike for instance in the content of river and coastal restoration schemes. The focus of this study was to investigate the effect of the presence of vegetation on flow resistance, turbulence statistics, and the instantaneous flow in open channels by performing three-dimensional computational-fluid-dynamics (CFD) simulations. Firstly, fully developed turbulent flow in fully-vegetated channel was analyzed by employing the method of high-resolution Large-Eddy Simulation (LES). Flow through a staggered array of rigid, emergent cylinders was simulated and the LES was validated through experiments. After validation, numerical simulations were performed at an extended parameter range of two different cylinder Reynolds numbers (ReD = 500 and 1340) and three different vegetation densities (φ = 0.016, 0.063, and 0.251). Flow structures and statistics were analyzed on the instantaneous flow and the effect of the vegetation density and cylinder Reynolds number was assessed. Moreover, drag forces exerted on the cylinders were calculated explicitly, and the effect of both ReD and φ on the drag coefficient was quantified. Secondly, two new alternative simulation strategies, a RANS based strategy with a vegetative closure model and a low-resolution Large-Eddy Simulation, were devised. They were evaluated by simulating several experimental cases with diverse conditions of the cylinder arrangement (i.e., staggered vs. random distribution), vegetation densities (φ = 0.016, 0.022, 0.063, 0.087, 0.091, 0.150, and 0.251), and cylinder Reynolds number (ReD = 170 - 1700). For the RANS based strategy, the importance of a-priori knowledge was assessed, and for the low-resolution LES, the efficiency and accuracy was demonstrated. Finally, a numerical strategy based on a porosity approach was developed and applied to open-channel flow through a natural plant. The simulated velocities were compared with experimentally acquired ones and results showed reasonable agreement. The results obtained in this research contribute to the understanding of fundamental mechanism of flow-vegetation interaction in vegetated open-channels, resolving turbulent flow-vegetation interaction explicitly. In addition, the new numerical strategies developed as part of this research are expected to allow describing the behavior of turbulent flow through artificial and natural vegetation with high efficiency and accuracy.

The project was aimed at describing the mixing phenomena relevant for both safety analysis, particularly in steam line break and boron dilution scenarios, and mixing phenomena of interest for economical operation and the structural integrity. Measurement data from a set of mixing experiments, gained by using advanced measurement techniques with enhanced resolution in time and space help to improve the basic understanding of turbulent mixing and to provide data for Computational Fluid Dynamics (CFD) code validation. Slug mixing tests simulating the start-up of the first main circulation pump are performed with two 1:5 scaled facilities: The Rossendorf coolant mixing model ROCOM and the VATTENFALL test facility, modelling a German Konvoi type and a Westinghouse type three-loop PWR, respectively. Additional data on slug mixing in a VVER-1000 type reactor gained at a 1:5 scaled metal mock-up at EDO Gidropress are provided. Experimental results on mixing of fluids with density differences obtained at ROCOM and the FORTUM PTS test facility are made available. Concerning mixing phenomena of interest for operational issues and thermal fatigue, flow distribution data available from commissioning tests (Sizewell-B for PWRs, Loviisa and Paks for VVERs) are used together with the data from the ROCOM facility as a basis for the flow distribution studies. The test matrix on flow distribution and steady state mixing performed at ROCOM comprises experiments with various combinations of running pumps and various mass flow rates in the working loops. Computational fluid dynamics calculations are accomplished for selected experiments with two different CFD codes (CFX-5, FLUENT). Best practice guidelines (BPG) are applied in all CFD work when choosing computational grid, time step, turbulence models, modelling of internal geometry, boundary conditions, numerical schemes and convergence criteria. The BPG contain a set of systematic procedures for quantifying and reducing numerical errors. The knowledge of these numerical errors is a prerequisite for the proper judgement of model errors. The strategy of code validation based on the BPG and a matrix of CFD code validation calculations have been elaborated. Besides of the benchmark cases, additional experiments were calculated by new partners and observers, joining the project later. Based on the "best practice solutions", conclusions on the applicability of CFD for turbulent mixing problems in PWR were drawn and recommendations on CFD modelling were given. The high importance of proper grid generation was outlined. In general, second order discretization schemes should be used to minimise numerical diffusion. First order schemes can provide physically wrong results. With optimised "production meshes" reasonable results were obtained, but due to the complex geometry of the flow domains, no fully grid independent solutions were achieved. Therefore, with respect to turbulence models, no final conclusions can be given. However, first order turbulence models like K-e or SST K-w are suitable for momentum driven slug mixing. For buoyancy driven mixing (PTS scenarios), Reynolds stress models provide better results.

<p>Conventional one-point turbulence closures have beenextended with an additional transported scalar for modeling ofmagnetohydrodynamic (MHD) turbulence. The new scalar, α ,captures the length scale anisotropy and tendency towardstwo-dimensionality, which is characteristic feature of MHDturbulence, and allows accurate modeling of the Jouledissipation of turbulence. The concept has been used for both afull Reynolds stress closure, and a three-equation<i>K-ε -α</i>model. An exact transport equation forαwas derived from the governing equations. All terms inthe equation require modeling, however. The proposed modeltransport equation for α includes terms for magneticdissipation, nonlinear energy transfer, and effects of meanshear and strain. Modeling of the magnetic and strain-relatedterms was based on rapid distortion analysis of the linearizedequations, while modeling of nonlinear effects isphenomenological in nature. For homogeneous turbulence, themodel was compared with linear theory, direct numericalsimulations and experiments. For turbulence subjected to astrong magnetic field, the model reproduces the energy andlength scale evolution predicted by linear theory. Whennonlinear effects are of importance, it predicts energy decayand length scale evolution in agreement with experiments. Theeddy viscosity and Reynolds stress versions of the modelcoincide with the respective conventional models in the absenceof a magnetic field. The objective of this project has been todevelop efficient MHD turbulence models for engineeringapplications, especially for modeling of continuous steelcasting. The novel MHD turbulence models appear to benumerically robust, and they have been implemented in acommercial flow solver, together with electromagnetic equationsfor the Lorentz forces in the mean momentum equations.</p><p><b>Keywords:</b>Turbulence model, magnetohydrodynamics, MHD,magnetohydrodynamic turbulence, computational fluid dynamics,continuous casting, dimensionality, Reynolds stresses, eddyviscosity</p>

No description available.

<p>With this contribution to the AVT-113/VFE-2 task group it was possible to prove the feasibility of high Reynolds number CFD computations to resolve and thus better understand the peculiar dual vortex system encountered on the VFE-2 blunt leading edge delta wing. Initial investigations into this phenomenon seemed to undermine the hypothesis, that the formation of the inner vortex system relies on the laminar state of the boundary layer at separation onset. As a result of this research, this initial hypothesis had to be expanded to account also for high Reynolds number cases, where a laminar boundary layer status at separation onset could be excluded. Furthermore, the data published in the same context shows evidence of secondary separation under the inner primary vortex. This further supports the supposition of a different generation mechanism of the inner vortical system other than a pure development out of a possibly laminar separation bubble. The unsteady computations performed on numerical grids with different levels of refinement led furthermore to the establishment of internal guidelines specific to the DES approach.</p>

No description available.

Thesis (M.S.)--Mississippi State University. Computational Engineering Program. / Title from title screen. Includes bibliographical references.