Theses and Dissertations from UMD
Permanent URI for this communityhttp://hdl.handle.net/1903/2
New submissions to the thesis/dissertation collections are added automatically as they are received from the Graduate School. Currently, the Graduate School deposits all theses and dissertations from a given semester after the official graduation date. This means that there may be up to a 4 month delay in the appearance of a give thesis/dissertation in DRUM
More information is available at Theses and Dissertations at University of Maryland Libraries.
Browse
3 results
Search Results
Item Numerical Simulations of Granular Physics in the Solar System(2017) Ballouz, Ronald; Richardson, Derek C; Astronomy; Digital Repository at the University of Maryland; University of Maryland (College Park, Md.)Granular physics is a sub-discipline of physics that attempts to combine principles that have been developed for both solid-state physics and engineering (such as soil mechanics) with fluid dynamics in order to formulate a coherent theory for the description of granular materials, which are found in both terrestrial (e.g., earthquakes, landslides, and pharmaceuticals) and extra-terrestrial settings (e.g., asteroids surfaces, asteroid interiors, and planetary ring systems). In the case of our solar system, the growth of this sub-discipline has been key in helping to interpret the formation, structure, and evolution of both asteroids and planetary rings. It is difficult to develop a deterministic theory for granular materials due to the fact that granular systems are composed of a large number of elements that interact through a non-linear combination of various forces (mechanical, gravitational, and electrostatic, for example) leading to a high degree of stochasticity. Hence, we study these environments using an N-body code, pkdgrav, that is able to simulate the gravitational, collisional, and cohesive interactions of grains. Using pkdgrav, I have studied the size segregation on asteroid surfaces due to seismic shaking (the Brazil-nut effect), the interaction of the OSIRIS-REx asteroid sample-return mission sampling head, TAGSAM, with the surface of the asteroid Bennu, the collisional disruptions of rubble-pile asteroids, and the formation of structure in Saturn's rings. In all of these scenarios, I have found that the evolution of a granular system depends sensitively on the intrinsic properties of the individual grains (size, shape, sand surface roughness). For example, through our simulations, we have been able to determine relationships between regolith properties and the amount of surface penetration a spacecraft achieves upon landing. Furthermore, we have demonstrated that this relationship also depends on the strength of the local gravity. By comparing our numerical results to laboratory experiments and observations by spacecraft we can begin to understand which microscopic properties (i.e., grain properties) control the macroscopic properties of the system. For example, we can compare the mechanical response of a spacecraft to landing or Cassini observations of Saturn's ring to understand how the penetration depth of a spacecraft or the complex optical depth structure of a ring system depends on the size and surface properties of the grains in those systems.Item A FLUID STRUCTURE INTERACTION STRATEGY WITH APPLICATION TO LOW REYNOLDS NUMBER FLAPPING FLIGHT(2010) Vanella, Marcos; Balaras, Elias; Mechanical Engineering; Digital Repository at the University of Maryland; University of Maryland (College Park, Md.)In this work a structured adaptive mesh refinement (S-AMR) strategy for fluid-structure interaction (FSI) problems in laminar and turbulent incompressible flows is developed. The Eulerian computational grid consists of nested grid blocks at different refinement levels. The grid topology and data-structure is managed by using the Paramesh© toolkit. The filtered Navier-Stokes equations are evolved in time by means of an explicit second-order projection scheme, where spatial derivatives are approximated with second order central differences on a staggered grid. The level of accuracy of the required variable interpolation operators is studied, and a novel divergence-preserving prolongation scheme for velocities is evolved. A novel direct-forcing embedded-boundary method is developed to enforce boundary conditions on a complex moving body not aligned with the grid lines. In this method, the imposition of no-slip conditions on immersed bodies is done on the Lagrangian markers that represent their wet surfaces, and the resulting force is transferred to the surrounding Eulerian grid points by a moving least squares formulation. Extensive testing and validation of the resulting strategy is done on a numerous set of problems. For transitional and turbulent flow regimes the large-eddy simulation (LES) approach is used. The grid discontinuities introduced in AMR methods lead to numerical errors in LES, especially if non-dissipative, centered schemes are used. A simple strategy is developed to vary the filter size for filtered variables around grid discontinuities. A strategy based on explicit filtering of the advective term is chosen to effectively reduce the numerical errors across refinement jumps. For all the FSI problems reported, the complete set of equations governing the dynamics of the flow and the structure are simultaneously advanced in time by using a predictor-corrector strategy. Dynamic fluid grid adaptation is implemented to reduce the number of grid points and computation costs. Applications to flapping flight comprise the study of flexibility effects on the aerodynamic performance of a hovering airfoil, and simulation of the flow around an insect model under prescribed kinematics and free longitudinal flight. In the airfoil simulations, it is found that peak performance is located in structural flexibility-inertia regions where non-linear resonances are present.Item Numerical Investigations of Gaseous Spherical Diffusion Flames(2009) Lecoustre, Vivien Renaud Francis; Sunderland, Peter B.; Mechanical Engineering; Digital Repository at the University of Maryland; University of Maryland (College Park, Md.)Spherical diffusion flames have several unique characteristics that make them attractive from experimental and theoretical perspectives. They can be modeled with one spatial dimension, which frees computational resources for detailed chemistry, transport, and radiative loss models. This dissertation is a numerical study of two classes of spherical diffusion flames: hydrogen micro-diffusion flames, emphasizing kinetic extinction, and ethylene diffusion flames, emphasizing sooting limits. The flames were modeled using a one-dimensional, time-accurate diffusion flame code with detailed chemistry and transport. Radiative losses from products were modeled using a detailed absorption/emission statistical narrow band model and the discrete ordinates method. During this work the code has been enhanced by the implementation of a soot formation/oxidation model using the method of moments. Hydrogen micro-diffusion flames were studied experimentally and numerically. The experiments involved gas jets of hydrogen. At their quenching limits, these flames had heat release rates of 0.46 and 0.25 W in air and in oxygen, respectively. These are the weakest flames ever observed. The modeling results confirmed the quenching limits and revealed high rates of reactant leakage near the limits. The effects of the burner size and mass flow rate were predicted to have a significant impact on the flame chemistry and species distribution profiles, favoring kinetic extinction. Spherical ethylene diffusion flames at their sooting limits were also examined. Seventeen normal and inverse spherical flames were considered. Initially sooty, these flames were experimentally observed to reach their sooting limits 2 s after ignition. Structure of the flames at 2 s was considered, with an emphasis on the relationships among local temperature, carbon to oxygen atom ratio (C/O), and scalar dissipation rate. A critical C/O ratio was identified, along with two different sooting limit regimes. Diffusion flames with local scalar dissipation rates below 2 1/s were found to have temperatures near 1410 K at the location of the critical C/O ratio, whereas flames with greater local scalar dissipation rate exhibited increased temperatures. The present work sheds light on important combustion phenomenon related to flame extinction and soot formation. Applications to energy efficiency, pollutant reduction, and fire safety are expected.