<p> A stream function finite element model is developed to solve for particle paths within a large ice mass. A steady-state primitive variable finite element model, treating ice as an incompressible non-Newtonian fluid, is used to furnish the necessary input velocities and rotations for the stream function finite element model. Time-integration along the particle paths is used to determine the age of the ice within the ice mass.</p> <p> Two ice masses are studied: the Barnes Ice Cap, Baffin Island, N.W.T., and Mount Logan, Yukon Territory. It is shown that if a realistic approximation of the velocity field of an ice mass can be established, the age of ice determined by time-integration along particle paths corresponds to the age determined by standard methods. Results of simulations using a transient model suggest that the elastic response of large ice masses is negligible.</p> / Thesis / Master of Engineering (MEngr)

AA7075 is a precipitation hardening structural aluminum alloy, which has garnered considerable interest in automotive industry, primarily due its lightweighting capacity compared to many other aluminum alloys from 2xxx and 6xxx series. However, the damage evolution in AA7075 is quite complex due to the presence of different second phase particles in the microstructure and their contribution on damage evolution is largely unknown at large plastic strains. The common second phase particles are η precipitates, θ precipitates and Fe-rich intermetallic particles. The current work presents an extensive multiscale numerical framework, which in conjunction with complementary experiments, is applied to study strain localization, void nucleation, growth, and coalescence in a particle rich matrix. Experimentally, void nucleation is observed to be driven by particle decohesion and particle fracture. Nanoscale molecular dynamics (MD) simulation is carried out to estimate interface properties of the three distinct particle types. The extracted properties are used as input for real particle field 2D and 3D microstructure based finite element (FE) models. The stochastic nature of particle fracture is described using a Weibull distribution, while the effect of grains is incorporated in terms of their Taylor factors. Ductile matrix is described using the well known Gurson Tvergaard Needleman (GTN) void damage model. Complementary experiments included uniaxial tensile tests carried out in-situ in Scanning Electron Microscope (SEM) and X-ray Computed Tomography (XCT), ex-situ high resolution XCT and Electron Back Scattered Diffraction (EBSD) tests. The FE models with three distinct particle stoichiometries and three competing damage mechanisms, show good agreement with experimental observations. Particle fracture marginally dominates particle decohesion. At low plastic strains, void nucleation is initiated by decohesion and fracture of larger Fe-rich particles, which facilitate formation of localized deformation bands. At large plastic strain, elevated stresses within the localized bands facilitate decohesion and fracture of more resistant η and θ precipitates. Due to their inherent larger size and more irregular morphology, θ precipitates contribute to voiding more than η precipitates. Under uniaxial tensile loads, void growth takes place in the middle of the specimen, driven by higher triaxiality stress state in the middle, relative to the surface. Void coalescence occurs along deformation bands driven by higher stresses due accumulated plastic strain within the bands, in a process known as void sheeting. / Thesis / Doctor of Philosophy (PhD)

A finite element solver has been developed for performing analysis and sensitivity analysis with Poisson's equation. An application of Poisson's equation in fluid dynamics is that of poential flow, in which case Posson's equaiton reduces to Laplace's equation. The stiffness matrix and sensitivity of the stiffness matrix are evaluated by direct integrations, as opposed to numerical integration. This allows less computational effort and minimizes the sources of computational errors. The capability of evaluating sensitivity derivatives has been added in orde to perform design sensitivity analysis of non-lifting airfoils. The discrete-direct approach to sensitivity analysis is utilized in the current work. The potential flow equations and the sensitivity equations are computed by using a preconditionaed conjugate gradient method. This method greatly reduces the time required to perfomr analysis, and the subsequent design optimization. Airfoil shape is updated at each design iteratioan by using a Bezier-Berstein surface parameterization. The unstrucured grid is adapted considering the mesh as a system of inteconnected springs. Numerical solutions from the flow solver are compared with analytical results obtained for a Joukowsky airfoil. Sensitivity derivaatives are validated using carefully determined central finite difference values. The developed software is then used to perform inverse design of a NACA 0012 and a multi-element airfoil.

No description available.

No description available.

No description available.

No description available.

No description available.

No description available.

No description available.