Settlement-Induced Bending Moments (SIBM) are an important design condition that must be considered whenever battered piles are placed in settling soils. The objective of this research was to investigate various parameters which can affect SIBM in battered piles within a levee embankment. The results from the current study were compared and verified against those obtained from centrifuge testing and alternative numerical simulations. A series of centrifuge testing as well as finite difference numerical simulations in Fast Langrangian Analysis of Continua (FLAC) were conducted. Different parameters which may affect the bending moments were investigated including pile connection fixity, batter, and stiffness of the pile as well as the magnitude of settlement. The simulations show that these parameters can have large impacts on the magnitude and location of the bending moments. Findings of this research can be used to validate or identify the need for adjustment of the current modeling/design approach.

A numerical simulation of electrowetting on a dielectric was performed in COMSOL to grant insight on various parameters that play a critical role in system performance. The specific system being simulated was the Open Drop experiment and the parameters being investigated were the applied voltage, contact angle at the advancing triple point, and droplet overlap onto neighboring actuated electrodes. These parameters were investigated with respect to their effect on droplet locomotion performance. This performance was quantified by the droplets velocity and the dielectrophortic (DEP) force’s magnitude; the DEP force was calculated from integration of the Maxwell Stress Tensor, however, the force was not integrated into the simulation to assist with droplet movement. It was found that as the droplet overlap onto the neighboring electrode, or droplet radius to electrode size ratio, decreased, the droplet velocity increased. As the applied potential increased, and induced contact angle at the advancing triple point decreased, droplet velocity also increased. Both the decreasing overlap and increasing voltage had a linear effect on droplet velocity. As the droplet overlap increased, the rate of change of droplet velocity decreased as increasing voltages were considered. A 2D DEP calculation illustrated that an increase in voltage induced a tenfold increase in the corresponding DEP force; a linear relationship was found between droplet overlap and DEP force for the Open Drop size regime.

No description available.

Castings of aluminum alloys are widely used in the automotive and aerospace industries since they play a significant role in improving the performance and fuel efficiency. In aluminum industries, sampling is the most common method to evaluate the inclusion levels which is a key indicator for the quality of the aluminum alloys. Since how the filling process and solidification process will influence the inclusion characteristics during the sampling procedure is of great importance, the objectives of this work is to create a the two-phase flow model to simulate the filling process and solidification process, as well as calculate the particles movement in the whole sampling procedure. Computational Fluid Dynamics (CFD) modeling was used and this work was performed in the software ANSYS FLUENT. A numerical two dimensional (2D) axisymmetric model was built to simulate the sampling procedure with the assumption that the filling could be done along the main axis automatically. First, the initial solidification during the filling was taken into account without particle injection. The realizable k − ε turbulence model was used to model the effects of the turbulence. Several simulations with different inlet filling rate, different initial filling temperature and different inlet diameter was calculated to see the influence on the solidification behavior. Then, the whole sampling system was modeled with particle injection. The Discrete Phase Model (DPM) was used to simulate the particle motion in the melt and the focus was on the influence of the initial solidification on the inclusion distributions. Finally, the optimal sampling position inside the aluminum sampler mold was calculated. / Gjutningar av aluminiumlegeringar används ofta inom bil-, och flygindustrin eftersom de spelar en viktig roll för att förbättra prestanda och bränsleeffektivitet. Inom aluminiumindustrin är provtagning den vanligaste metoden att utvärdera mängden inneslutningar i smältan, vilket är en nyckelindikator för kvaliteten på aluminiumlegeringarna. Eftersom både fyllnads- och stelningsprocessen kommer att påverka inneslutningskarakteristiken är provtagningsproceduren av stor betydelse. Syftet med detta arbete är att skapa en två-fasflödesmodell för att simulera fyllnings-, och stelningsprocessen, samt att beräkna partikelrörelserna under provtagningsförfarandet. Computational Fluid Dynamics (CFD) modellering användes och arbetet har utfördes med mjukvaran ANSYS FLUENT. En numerisk tvådimensionell (2D) axisymmetrisk modell byggdes för att simulera provtagningsproceduren med antagandet att påfyllningen kan göras automatiskt längs huvudaxeln. Till att börja med betraktades det första stelnandet under fyllningen utan partikelinjektion. En antagen k - ε turbulensmodell användes för att modellera effekten av turbulens. Flera simuleringar med olika inloppshastighet, påfyllningstemperatur och inloppsdiametrar beräknades för att se påverkan på stelningsbeteendet. Därefter modellerades hela provtagningsmodellen med partikelinjektion. En Diskret Fasmodell (DPM) användes för att simulera partikelrörelsen i smältan och fokus var inverkan av det initiala stelnandet på inneslutningsfördelningen. Slutligen beräknades den optimala provtagningspositionen inuti aluminiumprovformen.

Photonic crystals and nanophotonics have received a great deal of attention over the last decade, largely due to improved numerical modeling and advances in fabrication technologies. To this day, fabrication and optical behavior remain decoupled during the design phase and numerous assumptions are made about "perfect" geometry. As research moves from theory to real devices, predicting device behavior based on realistic geometry becomes critical. In this dissertation, a set of numerical tools was developed to model micro and nano fabrication processes. They were combined with equally capable tools to model optical performance of the simulated structures. Using these tools, it was predicted and demonstrated that 3D nanostructures may be formed on a standard mask aligner. A space-variant photonic crystal filter was designed and optimized based on a simple fabrication method of etching holes through hetero-structured substrates. It was found that hole taper limited their optical performance and a method was developed to compensate. A method was developed to tune the spectral response of guided-mode resonance filters at the time of fabrication using models of etching and deposition. Autocloning was modeled and shown that it could be used to form extremely high aspect ratio structures to improve performance of form-birefringent devices. Finally, the numerical tools were applied to metallic photonic crystal devices.

This thesis focuses on the simulation of astronomic tides of the Pascagoula River. The work is comprised of five steps: 1) Production of a digital elevation model describing the entire Pascagoula River system; 2) Development of an inlet-based, unstructured mesh for inbank flow to better understand the basis of the hydrodynamics within the Pascagoula riverine system. In order to assist in the mesh development, a toolbox was constructed to implement one-dimensional river cross sections into the two-dimensional model; 3) Implementation of a sensitivity analysis of the Pascagoula River two inlet system to examine the inlet effects on tidal propagation; 4) Improvement of the inlet-based model by performing a preliminary assessment of a spatially varied bottom friction; 5) Implementation of an advection analysis to reveal its influence on the flow velocity and water elevation within the domain. The hydrodynamic model employed for calculating tides is ADCIRC-2DDI (ADvanced CIRCulation Model for Shelves, Coasts and Estuaries, Two-Dimensional Depth Integrated). This finite element based model solves the shallow water equations in their full nonlinear form. Boundary conditions including water surface elevation at the off-shore boundary and tidal potential terms allow the full simulation of astronomic tides. The improved astronomic tide model showed strong agreement with the historical data at seven water level monitoring gauge stations. The main conclusions of this research are: 1) The western inlet of the Pascagoula River is more dominant than the eastern inlet; however, it is necessary to include both inlets in the model. 2) Although advection plays a significant role in velocity simulation, water elevations are insensitive to advection. 3) The astronomic model is sensitive to bottom friction (both global and spatial variations); therefore, a spatially varied bottom friction coefficient is suggested. As a result of this successful effort to produce an astronomic tide model of the Pascagoula River, a comprehensive storm surge model can be developed. With the addition of inundation areas the surge model can be expected to accurately predict storm tides generated by hurricanes along the Gulf Coast.

Optical filters encompass a vast array of devices and structures for a wide variety of applications. Generally speaking, an optical filter is some structure that applies a designed amplitude and phase transform to an incident signal. Different classes of filters have vastly divergent characteristics, and one of the challenges in the optical design process is identifying the ideal filter for a given application and optimizing it to obtain a specific response. In particular, it is highly advantageous to obtain a filter that can be seamlessly integrated into an overall device package without requiring exotic fabrication steps, extremely sensitive alignments, or complicated conversions between optical and electrical signals. This dissertation explores three classes of nano-scale optical filters in an effort to obtain different types of dispersive response functions. First, dispersive waveguides are designed using a sub-wavelength periodic structure to transmit a single TE propagating mode with very high second order dispersion. Next, an innovative approach for decoupling waveguide trajectories from Bragg gratings is outlined and used to obtain a uniform second-order dispersion response while minimizing fabrication limitations. Finally, high Q-factor microcavities are coupled into axisymmetric pillar structures that offer extremely high group delay over very narrow transmission bandwidths. While these three novel filters are quite diverse in their operation and target applications, they offer extremely compact structures given the magnitude of the dispersion or group delay they introduce to an incident signal. They are also designed and structured as to be formed on an optical wafer scale using standard integrated circuit fabrication techniques. A number of frequency-domain numerical simulation methods are developed to fully characterize and model each of the different filters. The complete filter response, which includes the dispersion and delay characteristics and optical coupling, is used to evaluate each filter design concept. However, due to the complex nature of the structure geometries and electromagnetic interactions, an iterative optimization approach is required to improve the structure designs and obtain a suitable response. To this end, a Particle Swarm Optimization algorithm is developed and applied to the simulated filter responses to generate optimal filter designs.

Detectors sensitive to thermal and reflected infrared radiation are widely used for night-vision, communications, thermography, and object tracking among other military, industrial, and commercial applications. System requirements for the next generation of ultra-high-performance infrared detectors call for increased functionality such as large formats (> 4K HD) with wide field-of-view, multispectral sensitivity, and on-chip processing. Due to the low yield of infrared material processing, the development of these next-generation technologies has become prohibitively costly and time consuming. In this work, it will be shown that physics-based numerical models can be applied to predictively simulate infrared detector arrays of current technological interest. The models can be used to a priori estimate detector characteristics, intelligently design detector architectures, and assist in the analysis and interpretation of existing systems. This dissertation develops a multi-scale simulation model which evaluates the physics of infrared systems from the atomic (material properties and electronic structure) to systems level (modulation transfer function, dense array effects). The framework is used to determine the electronic structure of several infrared materials, optimize the design of a two-color back-to-back HgCdTe photodiode, investigate a predicted failure mechanism for next-generation arrays, and predict the systems-level measurables of a number of detector architectures.

Modeling of friction stir welding (FSW) is challenging, as there are large gradients in both strain rate and temperature that must be accounted for in the constitutive law of the material being joined. Constitutive laws are most often calibrated using flow stresses from hot compression or hot torsion testing, where strain rates are much lower than those of the FSW process. As such, the current work employed a recently developed method to measure flow stresses in AA 2219-T67 at the high strain rates typical of FSW. These data were used in the development of a finite element simulation of FSW to study the effect of the new flow stress data on temperature, torque, and load predictions, compared to standard material models calibrated with hot compression or hot torsion data. It was found that load predictions were significantly better with the new material law, reducing the error with respect to experimental measurements by approximately 79%. Because heat generation during FSW is primarily a function of friction between the rapidly spinning tool and the workpiece, the choice of friction law, and associated parameters, were also studied with respect to FE model predictions. It was found that the Norton (viscoplastic) friction law was the most appropriate for modeling FSW, because its predictions were more accurate for both the transient and steady-state phases of the FSW plunge experiment. The postulated reason for the superior performance of the Norton law was its ability to account for temperature and rate sensitivity of the workpiece material sheared by the tool, while the Tresca limited Coulomb law favored contact pressure, with essentially no temperature or rate dependence of local material properties. The combination of the new flow stress data and the optimized Norton friction law resulted in a 63% overall reduction in model error, compared to the use of a standard material law and boilerplate friction parameters. The overall error was calculated as an equally weighted comparison of temperatures, torques, and forces with experimentally measured values.

We propose an efficient wideband theory for adjoint variable sensitivity analysis of problems with general anisotropic materials. The method is formulated based on the transmission line numerical modeling technique. The anisotropic material properties of potential interest are the full tensors of permittivity, permeability, electrical conductivity, magnetic resistivity, magnetoelectric coupling, and electromagnetic coupling. The tensors may contain non-diagonal elements. Our method estimates the gradients of the desired response with respect to all designable parameters using at most one extra simulation, regardless of their number. In contrast, in the conventional sensitivity analysis method using central finite differences, the number of the required simulations scales linearly with the number of designable parameters. The theory has been implemented for sensitivity analysis of the two and three-dimensional structures. The available adjoint variable method (AVM) sensitivities enable the optimization-based design of anisotropic and dispersive anisotropic structures. We apply our AVM technique to optimization-based wideband invisibility cloak design of arbitrary-shape objects. Our method optimizes the voxel-by-voxel constitutive parameters of an anisotropic cloak. This results in a large number of optimizable parameters. The associated sensitivities of a wideband cloaking objective function are efficiently estimated using our anisotropic adjoint variable method technique. A gradient-based optimization algorithm utilizes the available sensitivity information to iteratively minimize the visibility objective function and to determine the constitutive parameters of the optimal cloak. / Thesis / Doctor of Philosophy (PhD)