VTechWorks Repository :: Browsing by Author "Sun, Shu Ming"

Browsing by Author "Sun, Shu Ming"

Now showing 1 - 13 of 13

Almost well-posedness of the full water wave equation on the finite stripe domain
Zhu, Benben (Virginia Tech, 2023-08-18)
The dissertation gives a rigorous study of surface waves on water of finite depth subjected to gravitational force. As for `water', it is an inviscid and incompressible fluid of constant density and the flow is irrotational. The fluid is bounded above by a free surface separating the fluid from the air above (assumed to be a vacuum) and below by a rigid flat bottom. Then, the governing equations for the motion of the fluid flow are called Euler equations. If the initial fluid flow is prescribed at time zero, i.e., mathematically the initial condition for the Euler equations is given, the long-time existence of a unique solution for the Euler equations is still an open problem, even if the initial condition is small (or initial flow is almost motionless). The dissertation tries to make some progress for proving the long-time existence and show that the time interval of the existence is exponentially long, called almost global well-posedness, if the initial condition is small and satisfies some conditions. The main ideas for the study are from the corresponding almost global well-posedness result for surface waves on water of infinite depth.
Dynamics of High-Speed Planetary Gears with a Deformable Ring
Wang, Chenxin (Virginia Tech, 2019-10-17)
This work investigates steady deformations, measured spectra of quasi-static ring deformations, natural frequencies, vibration modes, parametric instabilities, and nonlinear dynamics of high-speed planetary gears with an elastically deformable ring gear and equally-spaced planets. An analytical dynamic model is developed with rigid sun, carrier, and planets coupled to an elastic continuum ring. Coriolis and centripetal acceleration effects resulting from carrier and ring gear rotation are included. Steady deformations and measured spectra of the ring deflections are examined with a quasi-static model reduced from the dynamic one. The steady deformations calculated from the analytical model agree well with those from a finite element/contact mechanics (FE/CM) model. The spectra of ring deflections measured by sensors fixed to the rotating ring, space-fixed ground, and the rotating carrier are much different. Planet mesh phasing significantly affects the measured spectra. Simple rules are derived to explain the spectra for all three sensor locations for in-phase and out-of-phase systems. A floating central member eliminates spectral content near certain mesh frequency harmonics for out-of-phase systems. Natural frequencies and vibration modes are calculated from the analytical dynamic model, and they compare well with those from a FE/CM model. Planetary gears have structured modal properties due to cyclic symmetry, but these modal properties are different for spinning systems with gyroscopic effects and stationary systems without gyroscopic effects. Vibration modes for stationary systems are real-valued standing wave modes, while those for spinning systems are complex-valued traveling wave modes. Stationary planetary gears have exactly four types of modes: rotational, translational, planet, and purely ring modes. Each type has distinctive modal properties. Planet modes may not exist or have one or more subtypes depending on the number of planets. Rotational, translational, and planet modes persist with gyroscopic effects included, but purely ring modes evolve into rotational or one subtype of planet modes. Translational and certain subtypes of planet modes are degenerate with multiplicity two for stationary systems. These modes split into two different subtypes of translational or planet modes when gyroscopic effects are included. Parametric instabilities of planetary gears are examined with the analytical dynamic model subject to time-varying mesh stiffness excitations. With the method of multiple scales, closed-form expressions for the instability boundaries are derived and verified with numerical results from Floquet theory. An instability suppression rule is identified with the modal structure of spinning planetary gears with gyroscopic effects. Each mode is associated with a phase index such that the gear mesh deflections between different planets have unique phase relations. The suppression rule depends on only the modal phase index and planet mesh phasing parameters (gear tooth numbers and the number of planets). Numerical integration of the analytical model with time-varying mesh stiffnesses and tooth separation nonlinearity gives dynamic responses, and they compare well with those from a FE/CM model. Closed-form solutions for primary, subharmonic, superharmonic, and second harmonic resonances are derived with a perturbation analysis. These analytical results agree well with the results from numerical integration. The analytical solutions show suppression of certain resonances as a result of planet mesh phasing. The tooth separation conditions are analytically determined. The influence of the gyroscopic effects on dynamic response is examined numerically and analytically.
Facilitating FPGA Reconfiguration through Low-level Manipulation
Zha, Wenwei (Virginia Tech, 2014-03-24)
The process of FPGA reconfiguration is to recompile a design and then update the FPGA configuration correspondingly. Traditionally, FPGA design compilation follows the way how hardware is compiled for achieving high performance, which requires a long computation time. How to efficiently compile a design becomes the bottleneck for FPGA reconfiguration. It is promising to apply some techniques or concepts from software to facilitate FPGA reconfiguration. This dissertation explores such an idea by utilizing three types of low-level manipulation on FPGA logic and routing resources, i.e. relocating, mapping/placing, and routing. It implements an FMA technique for "fast reconfiguration". The FMA makes use of the software compilation technique of reusing pre-compiled libraries for explicitly reducing FPGA compilation time. Based the software concept of Autonomic Computing, this dissertation proposes to build an Autonomous Adaptive System (AAS) to achieve "self-reconfiguration". An AAS absorbs the computing complexity into itself and compiles the desired change on its own. For routing, an FPGA router is developed. This router is able to route the MCNC benchmark circuits on five Xilinx devices within 0.35 ~ 49.05 seconds. Creating a routing-free sandbox with this router is 1.6 times faster than with OpenPR. The FMA uses relocating to load pre-compiled modules and uses routing to stitch the modules. It is an essential component of TFlow, which achieves 8 ~ 39 times speedup as compared to the traditional ISE flow on various test cases. The core part of an AAS is a lightweight embedded version of utilities for managing the system's hardware functionality. Two major utilities are mapping/placing and routing. This dissertation builds a proof-of-concept AAS with a universal UART transmitter. The system autonomously instantiates the circuit for generating the desired BAUD rate to adapt to the requirement of a remote UART receiver.
Free and Forced Vibration of Linearly Elastic and St. Venant-Kirchhoff Plates using the Third Order Shear and Normal Deformable Theory
Chattopadhyay, Arka Prabha (Virginia Tech, 2019-09-18)
Employing the Finite Element Method (FEM), we numerically study three problems involving free and forced vibrations of linearly and nonlinearly elastic plates with a third order shear and normal deformable theory (TSNDT) and the three dimensional (3D) elasticity theory. We used the commercial software ABAQUS for analyzing 3D deformations, and an in-house developed and verified software for solving the plate theory equations. In the first problem, we consider trapezoidal load-time pulses with linearly increasing and affinely decreasing loads of total durations equal to integer multiples of the time period of the first bending mode of vibration of a plate. For arbitrary spatial distributions of loads applied to monolithic and laminated orthotropic plates, we show that plates' vibrations become miniscule after the load is removed. We call this phenomenon as vibration attenuation. It is independent of the dwell time during which the load is a constant. We hypothesize that plates exhibit this phenomenon because nearly all of plate's strain energy is due to deformations corresponding to the fundamental bending mode of vibration. Thus taking the 1st bending mode shape of the plate vibration as the basis function, we reduce the problem to that of solving a single second-order ordinary differential equation. We show that this reduced-order model gives excellent results for monolithic and composite plates subjected to different loads. Rectangular plates studied in the 2nd problem have points on either one or two normals to their midsurface constrained from translating in all three directions. We find that deformations corresponding to several modes of vibration are annulled in a region of the plate divided by a plane through the constraining points; this phenomenon is termed mode localization. New results include: (i) the localization of both in-plane and out-of-plane modes of vibration, (ii) increase in the mode localization intensity with an increase in the length/width ratio of a rectangular plate, (iii) change in the mode localization characteristics with the fiber orientation angle in unidirectional fiber- reinforced laminae, (iv) mode localization due to points on two normals constrained, and (iv) the exchange of energy during forced harmonic vibrations between two regions separated by the line of nearly stationary points that results in a beating-like phenomenon in a sub-region of the plate. This technique can help design a structure with vibrations limited to its small sub-region, and harvesting energy of vibrations of the sub-region. In the third problem, we study finite transient deformations of rectangular plates using the TSNDT. The mathematical model includes all geometric and material nonlinearities. We compare the results of linear and nonlinear TSNDT FEM with the corresponding 3D FEM results from ABAQUS and note that the TSNDT is capable of predicting reasonably accurate results of displacements and in-plane stresses. However, the errors in computing transverse stresses are larger and the use of a two point stress recovery scheme improves their accuracy. We delineate the effects of nonlinearities by comparing results from the linear and the nonlinear theories. We observe that the linear theory over-predicts the deformations of a plate as compared to those obtained with the inclusion of geometric and material nonlinearities. We hypothesize that this is an effect of stiffening of the material due to the nonlinearity, analogous to the strain hardening phenomenon in plasticity. Based on this observation, we propose that the consideration of nonlinearities is essential in modeling plates undergoing large deformations as linear model over-predicts the deformation resulting in conservative design criteria. We also notice that unlike linear elastic plate bending, the neutral surface of a nonlinearly elastic bending plate, defined as the plane unstretched after the deformation, does not coincide with the mid-surface of the plate. Due to this effect, use of nonlinear models may be of useful in design of sandwich structures where a soft core near the mid-surface will be subjected to large in-plane stresses.
Galerkin Approximations of General Delay Differential Equations with Multiple Discrete or Distributed Delays
Norton, Trevor Michael (Virginia Tech, 2018-06-29)
Delay differential equations (DDEs) are often used to model systems with time-delayed effects, and they have found applications in fields such as climate dynamics, biosciences, engineering, and control theory. In contrast to ordinary differential equations (ODEs), the phase space associated even with a scalar DDE is infinite-dimensional. Oftentimes, it is desirable to have low-dimensional ODE systems that capture qualitative features as well as approximate certain quantitative aspects of the DDE dynamics. In this thesis, we present a Galerkin scheme for a broad class of DDEs and derive convergence results for this scheme. In contrast to other Galerkin schemes devised in the DDE literature, the main new ingredient here is the use of the so called Koornwinder polynomials, which are orthogonal polynomials under an inner product with a point mass. A main advantage of using such polynomials is that they live in the domain of the underlying linear operator, which arguably simplifies the related numerical treatments. The obtained results generalize a previous work to the case of DDEs with multiply delays in the linear terms, either discrete or distributed, or both. We also consider the more challenging case of discrete delays in the nonlinearity and obtain a convergence result by assuming additional assumptions about the Galerkin approximations of the linearized systems.
General boundary value problems of a class of fifth order KdV equations on a bounded interval
Sriskandasingam, Mayuran; Sun, Shu Ming; Zhang, Bing-yu Y. (2024)
Immersed Finite Elements for a Second Order Elliptic Operator and Their Applications
Zhuang, Qiao (Virginia Tech, 2020-06-17)
This dissertation studies immersed finite elements (IFE) for a second order elliptic operator and their applications to interface problems of related partial differential equations. We start with the immersed finite element methods for the second order elliptic operator with a discontinuous coefficient associated with the elliptic interface problems. We introduce an energy norm stronger than the one used in [111]. Then we derive an estimate for the IFE interpolation error with this energy norm using patches of interface elements. We prove both the continuity and coercivity of the bilinear form in a partially penalized IFE (PPIFE) method. These properties allow us to derive an error bound for the PPIFE solution in the energy norm under the standard piecewise $H^2$ regularity assumption instead of the more stringent $H^3$ regularity used in [111]. As an important consequence, this new estimation further enables us to show the optimal convergence in the $L^2$ norm which could not be done by the analysis presented in [111]. Then we consider applications of IFEs developed for the second order elliptic operator to wave propagation and diffusion interface problems. The first application is for the time-harmonic wave interface problem that involves the Helmholtz equation with a discontinuous coefficient. We design PPIFE and DGIFE schemes including the higher degree IFEs for Helmholtz interface problems. We present an error analysis for the symmetric linear/bilinear PPIFE methods. Under the standard piecewise $H^2$ regularity assumption for the exact solution, following Schatz's arguments, we derive optimal error bounds for the PPIFE solutions in both an energy norm and the usual $L^2$ norm provided that the mesh size is sufficiently small. {In the second group of applications, we focus on the error analysis for IFE methods developed for solving typical time-dependent interface problems associated with the second order elliptic operator with a discontinuous coefficient.} For hyperbolic interface problems, which are typical wave propagation interface problems, we reanalyze the fully-discrete PPIFE method in [143]. We derive the optimal error bounds for this PPIFE method for both an energy norm and the $L^2$ norm under the standard piecewise $H^2$ regularity assumption in the space variable of the exact solution. Simulations for standing and travelling waves are presented to corroborate the results of the error analysis. For parabolic interface problems, which are typical diffusion interface problems, we reanalyze the PPIFE methods in [113]. We prove that these PPIFE methods have the optimal convergence not only in an energy norm but also in the usual $L^2$ norm under the standard piecewise $H^2$ regularity.
Mathematical Analysis on the PEC model for Thixotropic Fluids
Wang, Taige (Virginia Tech, 2016-05-03)
A lot of fluids are more complex than water: polymers, paints, gels, ketchup etc., because of big particles and their complicated microstructures, for instance, molecule entanglement. Due to this structure complexity, some material can display that it is still in yielded state when the imposed stress is released. This is referred to as thixotropy. This dissertation establishes mathematical analysis on a thixotropic yield stress fluid using a viscoelastic model under the limit that the ratio of retardation time versus relaxation time approaches zero. The differential equation model (the PEC model) describing the evolution of the conformation tensor is analyzed. We model the flow when simple shearing is imposed by prescribing a total stress. One part of this dissertation focuses on oscillatory shear stresses. In shear flow, different fluid states corresponding to yielded and unyielded phases occur. We use asymptotic analysis to study transition between these phases when slow oscillatory shearing is set up. Simulations will be used to illustrate and supplement the analysis. Another part of the dissertation focuses on planar Poiseuille flow. Since the flow is spatially inhomogeneous in this situation, shear bands are observed. The flow is driven by a homogeneous pressure gradient, leading to a variation of stress in the cross-stream direction. In this setting, the flow would yield in different time scales during the evolution. Formulas linking the yield locations, transition width, and yield time are obtained. When we introduce Korteweg stress in the transition, the yield location is shifted. An equal area rule is identified to fit the shifted locations.
Mathematical models of immune responses following vaccination with application to Brucella infection
Kadelka, Mirjam Sarah (Virginia Tech, 2015-06-17)
For many years bovine brucellosis was a zoonosis endemic in large parts of the world. While it is still endemic in some parts, such as the Middle East or India, several countries such as Australia and Canada have successfully eradicated brucellosis in cattle by applying vaccines, improving the hygienic standards in cattle breeding, and slaughtering or quarantining infected animals. The large economical impact of bovine brucellosis and its virulence for humans, coming in direct contact to fluid discharges from infected animals, makes the eradication of bovine brucellosis important to achieve. To achieve this goal several vaccines have been developed in the past decades. Today the two most commonly used vaccines are Brucella abortus vaccine strain 19 and strain RB51. Both vaccines have been shown to be effective, but the mechanisms of immune responses following vaccination with either of the vaccines are not understood yet. In this thesis we analyze the immunological data obtained through vaccination with the two strains using mathematical modeling. We first design a measure that allows us to separate the subjects into good and bad responders. Then we investigate differences in the immune responses following vaccination with strain 19 or strain RB51 and boosting with strain RB51. We develop a mathematical model of immune responses that accounts for formation of antagonistic pro and anti-inflammatory and memory cells. We show that different characteristics of pro-inflammatory cell development and activity have an impact on the number of memory cells obtained after vaccination.
Modeling, Discontinuous Galerkin Approximation and Simulation of the 1-D Compressible Navier Stokes Equations
Grigorian, Zachary (Virginia Tech, 2019-08-20)
In this thesis we derive time dependent equations that govern the physics of a thermal fluid flowing through a one dimensional pipe. We begin with the conservation laws described by the 3D compressible Navier Stokes equations. This model includes all residual terms resulting from the 1D flow approximations. The final model assumes that all the residual terms are negligible which is a standard assumption in industry. Steady state equations are obtained by assuming the temporal derivatives are zero. We develop a semi-discrete model by applying a linear discontinuous Galerkin method in the spatial dimension. The resulting finite dimensional model is a differential algebraic equation (DAE) which is solved using standard integrators. We investigate two methods for solving the corresponding steady state equations. The first method requires making an initial guess and employs a Newton based solver. The second method is based on a pseudo-transient continuation method. In this method one initializes the dynamic model and integrates forward for a fixed time period to obtain a profile that initializes a Newton solver. We observe that non-uniform meshing can significantly reduce model size while retaining accuracy. For comparison, we employ the same initialization for the pseudo-transient algorithm and the Newton solver. We demonstrate that for the systems considered here, the pseudo-transient initialization algorithm produces initializations that reduce iteration counts and function evaluations when compared to the Newton solver. Several numerical experiments were conducted to illustrate the ideas. Finally, we close with suggestions for future research.
The Reflected Quasipotential: Characterization and Exploration
Farlow, Kasie Geralyn (Virginia Tech, 2013-05-06)
The Reflected Quasipotential V(x) is the solution to a variational problem that arises in the study of reflective Brownian motion. Specifically, the stationary distributions of reflected Brownian motion satisfy a large deviation principle (with respect to a spatial scaling parameter) with V(x) as the rate function. The Skorokhod Problem is an essential device in the construction and analysis of reflected Brownian motion and our value function V(x). Here we characterize V(x) as a solution to a partial differential equation H(DV(x))=0 in the positive n-dimensional orthant with appropriate boundary conditions. H(p) is the Hamiltonian and DV(x) is the gradient of V(x). V(x) is continuous but not differentiable in general. The characterization will need to be in terms of viscosity solutions. Solutions are not unique, thus additional qualifications will be needed for uniqueness. In order to prove our uniqueness result we consider a discounted version of V(x) in a truncated region and pass to the limit. In addition to this characterization of V(x) we explore the possibility of cyclic optimal paths in 3 dimensions.
A Two-Level Galerkin Reduced Order Model for the Steady Navier-Stokes Equations
Park, Dylan (Virginia Tech, 2023-05-15)
In this thesis we propose, analyze, and investigate numerically a novel two-level Galerkin reduced order model (2L-ROM) for the efficient and accurate numerical simulation of the steady Navier-Stokes equations. In the first step of the 2L-ROM, a relatively low-dimensional nonlinear system is solved. In the second step, the Navier-Stokes equations are linearized around the solution found in the first step, and a higher-dimensional system for the linearized problem is solved. We prove an error bound for the new 2L-ROM and compare it to the standard Galerkin ROM, or one-level ROM (1L-ROM), in the numerical simulation of the steady Burgers equation. The 2L-ROM significantly decreases (by a factor of 2 and even 3) the 1L-ROM computational cost, without compromising its numerical accuracy.
A viscoelastic constitutive model for thixotropic yield stress fluids: asymptotic and numerical studies of extension
Grant, Holly Victoria (Virginia Tech, 2017-11-17)
This dissertation establishes a mathematical framework for analyzing a viscoelastic model that displays thixotropic behavior as a model parameter gets very small. The model is the partially extending strand convection model, originally derived for polymeric melts that have long strands that get in the way of fully retracting. A Newtonian solvent is added. The uniaxial and equibiaxial extensional flows are studied using combined asymptotic analysis and numerical simulations. An initial value problem with a prescribed elongational stress is solved in the limit of large relaxation time. This gives rise to multiple time scales. If the initial stress is less than a critical value, the initial elastic elongation is followed by settling to an unyielded state at the slow time scale. If the initial stress is larger than the critical value, then yielding ensues. The extensional flows produce delayed yielding and hysteresis, both associated with thixotropy in complex fluids.

Browsing by Author "Sun, Shu Ming"

Results Per Page

Sort Options