The first two steps in the development of an integrated multidisciplinary design optimization procedure capable of analyzing the nonlinear fluid flow about geometrically complex aeroelastic configurations have been accomplished in the present work. For the first step, a three-dimensional unstructured grid approach to aerodynamic shape sensitivity analysis and design optimization has been developed. The advantage of unstructured grids, when compared with a structured-grid approach, is their inherent ability to discretize irregularly shaped domains with greater efficiency and less effort. Hence, this approach is ideally suited fro geometrically complex configurations of practical interest. In this work the time-dependent, nonlinear Euler equations are solved using an upwind, cell-centered, finite-volume scheme. The discrete, linearized systems which result from this scheme are solved iteratively by a preconditioned conjugate-gradient-like algorithm known as GMRES for the two-dimensional cases and a Gauss-Seidel algorithm for the three-dimensional; at steady-state, similar procedures are used to solve the accompanying linear aerodynamic sensitivity equations in incremental iterative form. As shown, this particular form of the sensitivity equation makes large-scale gradient-based aerodynamic optimization possible by taking advantage of memory efficient methods to construct exact Jacobian matrix-vector products. Various surface parameterization techniques have been employed in the current study to control the shape of the design surface. Once this surface has been deformed, the interior volume of the unstructured grid is adapted by considering the mesh as a system of interconnected tension springs. Grid sensitivities are obtained by differentiating the surface parameterization and the grid adaptation algorithms with ADIFOR, an advanced automatic-differentiation software tool. To demonstrate the ability of this procedure to analyze and design complex configurations of practical interest, the sensitivity analysis and shape optimization has been performed for several two- and three-dimensional cases. In two-dimensions, an initially symmetric NACA-0012 airfoil and a high-lift multi-element airfoil were examined. For the three-dimensional configurations, an initially rectangular wing with uniform NACA-0012 cross-sections was optimized; in additions, a complete Boeing 747-200 aircraft was studied. Furthermore, the current study also examines the effect of inconsistency in the order of spatial accuracy between the nonlinear fluid and linear shape sensitivity equations. The second step was to develop a computationally efficient, high-fidelity, integrated static aeroelastic analysis procedure. To accomplish this, a structural analysis code was coupled with the aforementioned unstructured grid aerodynamic analysis solver. The use of an unstructured grid scheme for the aerodynamic analysis enhances the interactions compatibility with the wing structure. The structural analysis utilizes finite elements to model the wing so that accurate structural deflections may be obtained. In the current work, parameters have been introduced to control the interaction of the computational fluid dynamics and structural analyses; these control parameters permit extremely efficient static aeroelastic computations. To demonstrate and evaluate this procedure, static aeroelastic analysis results for a flexible wing in low subsonic, high subsonic (subcritical), transonic (supercritical), and supersonic flow conditions are presented. / Ph. D.

The primary purpose of this dissertation is to investigate the stability of reticulated domes under multiple static and dynamic loads. Two elastic geometrically nonlinear structural models of a reticulated dome with 21 and 39 degrees of freedom are considered. The nonlinear response of the system to static loads is obtained using nonlinear programming and discrete perturbation techniques. The nonlinear programming technique is used to obtain a starting solution for the discrete perturbation technique and to optimize the choice of the perturbation parameter. Convergence criteria and error estimates to limit errors in a perturbation scheme are developed. A method for selecting a "suitable" perturbation parameter for imperfection sensitive systems is proposed. The investigation of stability of equilibrium of the system subjected to finite disturbances is based on the concept of "degree of stability" and the associated sufficient stability condition. The stability condition is derived from a theorem on extent of asymptotic stability of Liapunov's direct method of the theory of stability of motion. Its application requires the determination of the nonlinear fundamental path and the "nearest" unstable post-buckling path. This is obtained via static analysis. The perturbed motion of the system under a given set of perturbations is obtained by numerically integrating the nonlinear equations of motion. The dynamic stability tests confirm the sufficiency of the dynamic stability condition. However, they also indicate that there is a dynamic disturbance with a specific spatial distribution for which the sufficient condition of stability is also a necessary condition for each equilibrium state tested. Since in practice, the spatial distribution of the disturbances cannot be controlled, the sufficient dynamic stability condition employed is practical for the design of reticulated domes. The stability boundaries corresponding to two independent loads on the models are presented. Limit points lie on a boundary which is convex towards the region of stability. Bifurcation points lie on a continuous but piecewise differentiable boundary. Each piece of the boundary containing bifurcation points appears to be convex towards the region of stability. / Ph. D.

Wetting, which addresses either spontaneous or forced spreading of liquids on a solid surface, is a ubiquitous phenomenon in nature and can be observed by us on a daily basis, e.g., rain drops falling on a windshield and lubricants protecting our corneas. The study of wetting phenomena can be traced back to the observation of water rising in a capillary tube by Hauksbee in 1706 and still remains as a hot topic, since it lays the foundation for a wide spectrum of applications, such as fluid mechanics, surface chemistry, micro/nanofluidic devices, and phase change heat transfer enhancement. Generally, wetting is governed by the dynamic and static behaviors of the three-phase contact line. Therefore, a deep insight into the dynamics and statics of three-phase contact line at nanoscale is necessary for the technological advancement in nanotechnology and nanoscience. This dissertation aims to understand the dynamic wetting under a molecular kinetic framework and resolve the reconfiguration of liquid molecules at the molecular region of contact line. Water spreading on polytetrafluoroethylene surfaces is selected as a classical example to study the dynamic behaviors of three-phase contact line. To accommodate the moving contact line paradox, the excess free energy is considered to be dissipated in the form of molecular dissipation. As-formed contact line friction/dissipation coefficient is calculated for water interacting with PTFE surfaces with varying structures and is found to be on the same order of magnitude with dynamic viscosity. From an ab initio perspective, contact line friction is decomposed into contributions from solid-liquid retarding and viscous damping. A mathematical model is established to generalize the overall friction between a droplet and a solid surface, which is able to clarify the static-to-kinetic transition of solid-liquid friction without introducing contact angle hysteresis. Moreover, drag reduction on lotus-leaf-like surface is accounted for as well. For the first time, the concept of contact line friction is used in the rational design of a superhydrophobic condenser surface for continuous dropwise condensation. We focus on the transport and reconfiguration of liquid molecules confined by a solid wall to shed light on the morphology of the molecular region of a three-phase contact line. A governing equation, which originates from the free energy analysis of a nonuniform monocomponent system, is derived to describe the patterned oscillations of liquid density. By comparing to the Reynolds transport theorem, we find that the oscillatory profiles of interfacial liquids are indeed governed in a combined manner by self-diffusion, surface-induced convection and shifted glass transition. Particularly for interfacial water, the solid confining effects give rise to a bifurcating configuration of hydrogen bonds. Such unique configuration consists of repetitive layer-by-layer water sheets with intra-layer hydrogen bonds and inter-layer defects. Molecular dynamics simulations on the interfacial configuration of water on solid surfaces reveal a quadratic dependence of adhesion on solid-liquid affinity, which bridges the gap between macroscopic interfacial properties and microscopic parameters. / Doctor of Philosophy / The study of wetting phenomena can be traced back to the observation of water rising in a capillary tube by Hauksbee in 1706 and still remains as a hot topic, since it lays the foundation for a wide spectrum of applications, such as fluid mechanics, surface chemistry, micro/nanofluidic devices, and phase change heat transfer enhancement. The conventional hydrodynamic analysis with no slip boundary condition predicts a diverging shear stress at the contact line as well as an unbounded shear force exerted on the solid surface. To accommodate this paradox, different mechanism and models have been proposed to clarify the slip between a moving contact line and a solid surface. Although almost all models yield reasonable agreement with experimental observations or numerical simulations, it is still difficult to pick up a specific model using only macroscopic properties and experiment-accessible quantities, because the energy dissipation mechanism during dynamic wetting is not identified and the contact line deforms over different length scales. In this dissertation, we ascribe the energy dissipation in dynamic wetting to contact line friction/dissipation under the framework of molecular kinetic theory, as it is assumed that the contact line is constantly oscillating around its equilibrium position. By decomposing contact line friction into two contributions: solid-liquid retarding and viscous damping, we are able to derive a universal model for the contact line friction. This model predicts a decaying solid-liquid friction on superhydrophobic surfaces, corresponding to the lotus effect. In the meantime, this model is able to clarify the recently-discovered static-to-kinetic transition of frictional force between a sessile drop and a solid surface. Later, we used the concept of contact line friction in the droplet growth process in dropwise condensation so as to promote the rational design of superhydrophobic condenser surfaces for sustainable dropwise condensation. As the morphology of a contact line is dependent on the length scale of interest, we focus on the molecular region of contact line. We study the transport and structural change of liquid molecules that are several molecular layers away from the solid surface. It is found that liquid molecules in this region experience patterned density oscillations, which cannot be simply attributed to the random deviations from continuum limit. By combining free energy analysis and Reynolds transport theorem, it is demonstrated that the omnipresent density oscillations arise from the collective effects of self-diffusion, surface-induced convection and shifted glass transition. For liquid water, we propose a bifurcating hydrogen bonding network in contrast to its tetrahedron configuration in bulk water.

Software testing is a fundamental part of the software development process. However, testing is still costly and consumes about half of the development cost. The path explosion problem often necessitates one to consider an extremely large number of paths in order to reach a specific target. Symbolic execution can reduce this cost by using symbolic values and heuristic exploration strategies. Although various exploration strategies have been proposed in the past, the number of Satisfiability Modulo Theories (SMT) solver calls for reaching a target is still large, resulting in longer execution times for programs containing many paths. In this paper, we present two partitioning strategies in order to mitigate this problem, consequently reducing unnecessary SMT solver calls as well. In sequential partitioning, code sections are analyzed sequentially to take advantage of infeasible paths discovered in earlier sections. On the other hand, using dynamic partitioning on SSA-applied code, the code sections are analyzed in a non-consecutive order guided by data dependency metrics within the sections. Experimental results show that both strategies can achieve significant speedup in reducing the number of unnecessary solver calls in large programs. More than 1000x speedup can be achieved in large programs over conflict-driven learning. / Master of Science

High transmission data rate, spectral efficiency, and reliability are necessary for future wireless communications systems. In a multipath-rich wireless channel, deploying multiple antennas at both the transmitter and receiver achieves high data rate, without increasing the total transmission power or bandwidth. When perfect knowledge of the wireless channel conditions is available at the receiver, the capacity has been shown to grow linearly with the number of antennas. However, the channel conditions must be estimated since perfect channel knowledge is never known a priori. In practice, the channel estimation procedure can be aided by transmitting pilot symbols that are known at the receiver. System performance depends on the quality of channel estimate, and the number of pilot symbols. It is desirable to limit the number of transmitted pilot symbols because pilot symbols reduce spectral efficiency. This thesis analyzes the system performance of coded multiple-input multiple-output (MIMO) systems for the quasi-static fading channel. The assumption that perfect channel knowledge is available at the receiver must be removed, in order to more accurately examine the system performance. Emphasis is placed on developing channel estimation strategies for an iterative Vertical Bell-Labs Layered Space Time (V-BLAST) architecture. The channel estimate can be sequentially improved between successive iterations of the iterative V-BLAST algorithm. For both the coded and uncoded systems, at high signal to noise ratio only a minimum number of pilot symbols per transmit antenna are required to achieve perfect channel knowledge performance. / Master of Science

Static analysis tools can be used to help programmers identify problems in their code. However, these tools often assume that developers have some programming background knowledge, so they can be hard to use in an educational context. We investigated the most common FindBugs errors from student code submissions and determined those errors that were related to incorrect solutions to problems and potential struggling for students. FindBugs is a static analysis tool that looks for incorrect patterns in Java bytecode analysis to identify potential coding flaws. For the common errors, we rewrote some of the original FindBugs messages to help students more easily understand the problems with their code. We found that students with at least one FindBugs warning in their final submission to an assignment had more submissions, longer work times, and lower correctness scores than students who did not have a FindBugs warning in their final submission. Adding modified FindBugs feedback to the automated grader resulted in students making fewer submissions and decreasing the length of time required to complete assignments. / Master of Science / Professional software developers use automated tools when they code to help them catch potential coding problems. These tools are difficult for novice student programmers because they do not have the same level of background as professionals. In this work, we attempted to change the feedback given by these tools so that students could understand it and use it to fix their code. We found that, across all of the undergraduate courses in this study, FindBugs warnings were associated with students having more trouble with assignments. When students could see FindBugs warnings, their time to complete assignments and the number of attempts they made both went down.

No / To explore the performance of patient-specific prior information, for example, from structural imaging, in improving perimetric procedures. Computer simulation was used to determine the error distribution and presentation count for Structure–Zippy Estimation by Sequential Testing (ZEST), a Bayesian procedure with prior distribution centered on a threshold prediction from structure. Structure-ZEST (SZEST) was trialled for single locations with combinations of true and predicted thresholds between 1 to 35 dB, and compared with a standard procedure with variability similar to Swedish Interactive Thresholding Algorithm (SITA) (Full-Threshold, FT). Clinical tests of glaucomatous visual fields (n = 163, median mean deviation −1.8 dB, 90% range +2.1 to −22.6 dB) were also compared between techniques. For single locations, SZEST typically outperformed FT when structural predictions were within ± 9 dB of true sensitivity, depending on response errors. In damaged locations, mean absolute error was 0.5 to 1.8 dB lower, SD of threshold estimates was 1.2 to 1.5 dB lower, and 2 to 4 (29%–41%) fewer presentations were made for SZEST. Gains were smaller across whole visual fields (SZEST, mean absolute error: 0.5 to 1.2 dB lower, threshold estimate SD: 0.3 to 0.8 dB lower, 1 [17%] fewer presentation). The 90% retest limits of SZEST were median 1 to 3 dB narrower and more consistent (interquartile range 2–8 dB narrower) across the dynamic range than those for FT. Seeding Bayesian perimetric procedures with structural measurements can reduce test variability of perimetry in glaucoma, despite imprecise structural predictions of threshold. Structural data can reduce the variability of current perimetric techniques. A strong structure–function relationship is not necessary, however, structure must predict function within ±9 dB for gains to be realized.

The energy absorption capacity of a series of axially crushed composite tubes fabricated from high tow count graphite fiber is compared with those of similar tubes fabricated from aerospace-grade fiber to determine the viability of considering the use of such fibers in automotive applications. To that end, graphite-epoxy tubular specimens with circular and square cross-sectional geometries; stacking sequences with ±45° fibers and with both ±45° and 0° fibers; and two different fiber types were fabricated and crushed both statically and dynamically to examine the energy absorption characteristics. The fiber types, in the form of preimpregnated tow (towpreg) from Thiokol, were Akzo Fortafil 50k fiber and aerospace-grade T300 12k fiber. Using the towpreg, sixteen tubes were filament wound on aluminum mandrels. Three specimens were cut from each of these tubes for a total of forty-eight specimens. Twenty-four of these specimens were crushed statically in a load frame and twenty-four were crushed dynamically in a drop fixture. In order to characterize the tubes and specimens, a number of measurements were taken. These measurements included length, wall thickness, cross-sectional dimensions, volume, and mass. Two important energy absorption measures were examined: the specific energy absorption (SEA) and the ratio of the peak load to the average load. The geometry had a significant effect on the energy absorption but the stacking sequence did not. It was also found that the 50k material was less effective at absorbing energy than the 12k material, but the 50k still may be acceptable. / Master of Science

The present study is part of the Long Term Bridge Performance Program (LTBP) funded by the Federal Highway Administration. The objectives of this program are to create a comprehensive database of quantitative information of the long-term performance of selected pilot bridges and to develop a methodology to assess bridge performance. Finite element (FE) modeling of the pilot bridges is an intrinsic part of the LTBP program and is intended to not only assist with instrumentation decisions, but also to provide further insight into the behavior of these bridges, which cannot be achieved solely from field testing of the bridges. This thesis provides a comprehensive study of a plethora of issues associated with the development of reliable and accurate FE models of bridges. The first objective of this investigation was to develop reliable finite element models with a variety of levels of refinement and to study the effect of the inclusion of various bridge parameters in the model, such as bridge skew, degree of composite action, thermal gradient and level of support restraint, on the response of bridges. First, the suitability of different modeling techniques and of elements used to model the primary bridge components was assessed using simple models for which analytical solutions are readily available. From these studies, it was concluded that shell elements are adequate to model the bridge deck, and beam and shell elements are both satisfactory to model the bridge girders. From the dynamic analyses of theWildcat Creek River Bridge and the Colquitz River Bridge, flexural modes of vibration were found to be highly sensitive to support restraints and to how the guardrails were modeled and less sensitive to the inclusion of bracing and thermal gradients in the model. The finite element models using extreme boundary conditions were successful in bracketing the field response. The factors identified from these analyses were considered in the analysis of the Virginia pilot bridge. Different support restraints, and the inclusion of skew and level of composite action in the model had noticeable impact on both the static and dynamic responses of the bridge. The results from these analyses were used to assist with instrumentation decisions prior to field-testing. The developed model will also be used to help researchers further understand the bridge's behavior and to help explain a variety of phenomena observed in the field. / Master of Science

Considerable attention has been paid to the design of low-power, high-performance SRAMs (Static Random Access Memories) since they are a critical component in both hand-held devices and high-performance processors. A key in improving the performance of the system is to use an optimum sized SRAM. In this thesis, an SRAM compiler has been developed for the automatic layout of memory elements in the ASIC environment. The compiler generates an SRAM layout based on a given SRAM size, input by the user, with the option of choosing between fast vs. low-power SRAM. Array partitioning is used to partition the SRAM into blocks in order to reduce the total power consumption. Experimental results show that the low-power SRAM is capable of functioning at a minimum operating voltage of 2.1 V and dissipates 17.4 mW of average power at 20 MHz. In this report, we discuss the implementation of the SRAM compiler from the basic component to the top-level SKILL code functions, as well as simulation results and discussion. / Master of Science