Ice accretion is a purely unsteady phenomenon that is presently approximated by most icing codes using quasi-steady modeling. The accuracy of ice prediction is thus directly related to the arbitrarily prescribed time span during which the impact of ice growth on both flow and droplets is neglected. The objective of this work is to remove this limitation by implementing a cost-effective unsteady approach. This is done by fully coupling, in time, a diphasic flow (interacting air and droplet particles) with the ice accretion model. The two-phase flow is solved using the Navier-Stokes and Eulerian droplet equations with dual-time stepping in order to improve computational time. The ice shape is either obtained from the conservation of mass and energy within a thin film layer for glaze and mixed icing conditions, or from a mass balance between water droplets impingement and mass flux of ice for rime icing conditions. The iced surface being constantly displaced in time, Arbitrary Lagrangian-Eulerian terms are added to the governing equations to account for mesh movement. Moreover, surface smoothing techniques are developed to prevent degradation of the iced-surface geometric discretization. For rime ice, the numerical results clearly show that the new full unsteady modeling improves the accuracy of ice prediction, compared to the quasi-steady approach, while in addition ensuring time span independence. The applicability of the unsteady icing model for predicting glaze ice accretion is also demonstrated by coupling the diphasic model to the Shallow Water Icing Model. A more rigorous analysis reveals that this model requires the implementation of local surface roughness and that previous quasi-steady validations cannot be carried out using a small number of shots, therefore the need for unsteady simulation.

Even though the power of supercomputers has increased extraordinarily, there is still an insatiable need for more advanced multi-disciplinary CFD simulations in the aircraft analysis and design fields. A particular current interest is in the realistic three-dimensional fully viscous turbulent flow simulation of the highly non-linear aspects of aero-icing. This highly complex simulation is still computationally too demanding in industry, especially when several runs, such as parametric studies, are needed. In order to make such compute-intensive simulations more affordable, this work presents a reduced order modeling approach, based on the "Proper Orthogonal Decomposition", (POD), method to predict a wider swath of flow fields and ice shapes based on a limited number of "snapshots" obtained from complete high-fidelity CFD computations. The procedure of the POD approach is to first decompose the fields into modes, using a limited number of full-calculations snapshots, and then to reconstruct the field and/or ice shapes using those decomposed modes for other conditions, leading to reduced order calculations. The use of the POD technique drastically reduces the computational cost and can provide a more complete map of the performance degradation of an iced aircraft over a wide range of flight and weather conditions.

This work is concerned with the characterization of slowly moving fluids and was carried out on the flow of water through a narrow sedimentation tank. Dispersion in the type of flow structure involved is caused mainly by the presence of large eddies and, due to the fact that shear stresses are small, these eddies persist for a considerable period of time. Two flow models are presented : The first model assumes the X- Y- velocity component pair to form a discrete state Markov process in time and dispersion equations for the mean concentration at a point, the variance as well as concentration cross correlations are generated. In the second model the velocity fluctuation components are assumed to be independent, time-stationary Markov processes with normal probability density functions. The stochastic differential equation describing dispersion of tracer is formulated with and without the effect of molecular diffusion and solutions to both cases are presented. Comparison of the model with experimental data obtained from tracer and anemometer measurements show that the model is capable of describing mean dispersion in a relatively small region of the tank and that the tracer experiments were insensitive to molecular diffusion. / Thesis (Ph.D.)-University of Natal, Durban, 1970.

No description available.

No description available.

No description available.

A three-dimensional numerical model (PHOENICS) was used to investigate the role of stream deflector angle and length on the flow field in a rectangular laboratory flume. Subsequent bed topography surveys were performed to examine the role of obstruction angle on scour hole development over time. The model was capable of predicting laboratory velocity and turbulent kinetic energy measurements, performing better for flow over a flat stable bed than over a deformed sand bed. A new method of incorporating complex bed topography into a structured Cartesian mesh was developed in the process. Flow field properties such as dynamic pressure, velocity amplification, separation zone length and width, and downwelling extent and magnitude were found to be strongly dependent on deflector geometry. Equilibrium scour hole depths and geometry are also angle-dependent. A predictive equation was produced explaining the rate at which scour holes reach equilibrium, and compared well with existing literature. Finally, a method was developed whereby characteristics of the flow field over a flat, stable bed could be used to predict equilibrium scour hole geometry.

No description available.

Thesis (Ph.D)--Mathematics, Georgia Institute of Technology, 2008. / Committee Chair: Goldsztein, Guillermo; Committee Member: Dieci, Luca; Committee Member: McCuan, John; Committee Member: Santamarina, Juan; Committee Member: Zhou, Haomin.

This work is focused on the numerical modelling of fluid-structure interaction in three dimensions. Both internal and external laminar flow around flexible bodies are considered. The fluid flow simulated is based on the incompressible Navier-Stokes equations and the general focus is on laminar Newtonian flow. The streamline upwind/ pressure stabilising Petrov-Galerkin (SUPG/PSPG) method is employed to achieve a stable low order finite element discretisation of the fluid, while the solid is discretised spatially by a standard Galerkin finite element approach. The behavior of the solid is governed by Neo-Hooke elasticity. For temporal discretisation the discrete implicit generalised-alpha method is employed for both the fluid and the solid domains. The motion of the fluid mesh is solved using an arbitrary Lagrangian-Eulerian (ALE) scheme employing a nonlinear pseudo-elastic mesh update method. The fluid-solid interface is modelled using a finite element interpolation method that allows for non-matching meshes and satisfies the required conservation laws. The resulting sets of fully implicit strongly coupled nonlinear equations are then decomposed into a general framework consisting of fluid, interface and solid domains. These equations are then solved using different solution techniques consisting of strongly coupled monolithic Newton and block Gauss-Seidel methods as well as a weakly coupled novel staggered scheme. These solvers are employed to solve a number of three dimensional numerical examples consisting of: External flow: o a soft elastic beam fixed at both ends o a thin cantilever plate.