Partial Differential Equations Commons^™

Enhanced Stability Of A Dewetting Thin Liquid Film In A Single-Frequency Vibration Field, Sergey Shklyaev, Mikhail Khenner, Alexei Alabuzhev

Mathematics Faculty Publications

Dynamics of a thin dewetting liquid film on a vertically oscillating substrate is considered. We assume moderate vibration frequency and large (compared to the mean film thickness) vibration amplitude. Using the lubrication approximation and the averaging method, we formulate the coupled sets of equations governing the pulsatile and the averaged fluid flows in the film, and then derive the nonlinear amplitude equation for the averaged film thickness. We show that there exists a window in the frequency-amplitude domain where the parametric and shear-flow instabilities of the pulsatile flow do not emerge. As a consequence, in this window the averaged description …

Morphologies And Kinetics Of A Dewetting Ultrathin Solid Film, Mikhail Khenner

Mathematics Faculty Publications

The surface evolution model based on geometric partial differential equation is used to numerically study the kinetics of dewetting and dynamic morphologies for the localized pinhole defect in the surface of the ultrathin solid film with the strongly anisotropic surface energy. Depending on parameters such as the initial depth and width of the pinole, the strength of the attractive substrate potential and the strength of the surface energy anisotropy, the pinhole may either extend to the substrate and thus rupture the film, or evolve to the quasiequilibrium shape while the rest of the film surface undergoes phase separation into a …

Numerical Simulation Of Waves And Fronts In Inhomogeneous Solids, A. Berezovski, M. Berezovski, J. Engelbrecht, G. A. Maugin

Publications

Dynamic response of inhomogeneous materials exhibits new effects, which often do not exist in homogeneous media. It is quite natural that most of studies of wave and front propagation in inhomogeneous materials are associated with numerical simulations. To develop a numerical algorithm and to perform the numerical simulations of moving fronts we need to formulate a kinetic law of progress relating the driving force and the velocity of the discontinuity. The velocity of discontinuity is determined by means of the non-equilibrium jump relations at the front. The obtained numerical method generalizes the wave-propagation algorithm to the case of moving discontinuities …

Numerical Simulation Of Nonlinear Elastic Wave Propagation In Piecewise Homogeneous Media, Arkadi Berezovski, Mihhail Berezovski, Juri Engelbrecht

Publications

Systematic experimental work [S. Zhuang, G. Ravichandran, D. Grady, J. Mech. Phys. Solids 51 (2003) 245–265] on laminated composites subjected to high velocity impact loading exhibits the dispersed wave field and the oscillatory behavior of waves with respect to a mean value. Such a behavior is absent in homogeneous solids. An approximate solution to the plate impact in layered heterogeneous solids has been developed in [X. Chen, N. Chandra, A.M. Rajendran, Int. J. Solids Struct. 41 (2004) 4635–4659]. The influence of the particle velocity on many process characteristics was demonstrated. Based on earlier results [A. Berezovski, J. Engelbrecht, G.A. Maugin, …

Lattice Quantum Algorithm For The Schrodinger Wave Equation In 2+1 Dimensions With A Demonstration By Modeling Soliton Instabilities, Jeffrey Yepez, George Vahala, Linda L. Vahala

Electrical & Computer Engineering Faculty Publications

A lattice-based quantum algorithm is presented to model the non-linear Schrödinger-like equations in 2 + 1 dimensions. In this lattice-based model, using only 2 qubits per node, a sequence of unitary collide (qubit-qubit interaction) and stream (qubit translation) operators locally evolve a discrete field of probability amplitudes that in the long-wavelength limit accurately approximates a non-relativistic scalar wave function. The collision operator locally entangles pairs of qubits followed by a streaming operator that spreads the entanglement throughout the two dimensional lattice. The quantum algorithmic scheme employs a non-linear potential that is proportional to the moduli square of the wave function. …

Transient Non-Linear Heat Conduction Solution By A Dual Reciprocity Boundary Element Method With An Effective Posteriori Error Estimator, Eduardo Divo, Alain J. Kassab

Publications

A Dual Reciprocity Boundary Element Method is formulated to solve non-linear heat conduction problems. The approach is based on using the Kirchhoff transform along with lagging of the effective non-linear thermal diffusivity. A posteriori error estimate is used to provide effective estimates of the temporal and spatial error. A numerical example is used to demonstrate the approach.

Dynamics And Control Of A Three Dimensional Gantry Crane With Cable Flexibility, Uchendu H. Eke

Mechanical & Aerospace Engineering Theses & Dissertations

The control of payload swing in industrial gantry cranes is a topic of widespread interest. In this thesis, the control of the swing of a payload idealized by a point mass assumption is investigated utilizing s, compact robotics matrix/vector representation of the dynamical model. In this form, this particular system can be viewed as sn underactuated manipulator with equal numbers of active and passive degrees of freedom. The full three dimensional model is developed and linearized by modelling the cable as a single thread of fiexible wire and the degree of dynamic coupling between the active and passive coordinates is …

Droplet Evaporation And Deformations In An Amplitude Modulated Ultrasonic Field, Nihad E. Daidzic, Rene Stadler, Adrian Melling

Aviation Department Publications

The aim of the report presented is the measurements of droplet oscillations.

Wing-Section Optimization For Supersonic Viscous Flow, Cem Cihan Item

Mechanical & Aerospace Engineering Theses & Dissertations

The recent interest in the High Speed Civil Transport (HSCT) has resulted in renewed research studies of optimized supersonic cruise transport configurations. Incorporation of flow viscosity effects in the design process of such a supersonic wing is currently under investigation. This may lead to more accurate problem formulations and, in tum, greater aerodynamic efficiency than can be obtained by the traditional, inviscid, linear theories. In this context, for a design code to be a candidate for a complex optimization problem, such as three-dimensional viscous supersonic wing design, it should be validated using simpler building-block shapes.

To optimize the shape of …

A Numerical Analysis Of Smoothed Particle Hydrodynamics, David A. Fulk

Theses and Dissertations

This dissertation studies the numerical method of Smoothed Particle Hydrodynamics SPH as a technique for solving systems of conservation equations. The research starts with a detailed consistency analysis of the method. Higher dimensions and non-smooth functions are considered in addition to the smooth one dimensional case. A stability analysis is then performed. Using a linear technique, an instability is found. Solutions are proposed to resolve the instability. Also a total variation stability analysis is performed leading to a monotone form of SPH. The concepts of consistency and stability are then used in a convergence proof. This proof uses lemmas derived …

Implementation Of A Multiblock Sensitivity Analysis Method In Numerical Aerodynamic Shape Optimization, James Matthew Lacasse

Mechanical & Aerospace Engineering Theses & Dissertations

A multiblock sensitivity analysis method is applied in a numerical aerodynamic shape optimization technique. The Sensitivity Analysis Domain Decomposition (SADD) scheme which is implemented in this study was developed to reduce the computer memory requirements resulting from the aerodynamic sensitivity analysis equations, Discrete sensitivity analysis offers the ability to compute quasi-analytical derivatives in a more efficient manner than traditional finite-difference methods, which tend to be computationally expensive and prone to inaccuracies.

The direct optimization procedure couples CFD analysis based on the two-dimensional thin-layer Navier-Stokes equations with a gradient-based numerical optimization technique. The linking mechanism is the sensitivity equation derived from …

A Direct Numerical Simulation Of Vortex Breakdown Induced Tail Buffet, Steven James Massey

Mechanical & Aerospace Engineering Theses & Dissertations

A simulation of tail buffet is presented for a delta wing-vertical tail configuration. Flow conditions are chosen such that the wing primary-vortex cores experience vortex breakdown and the resulting turbulent wake flow impinges on the vertical tail. The dimensions and material properties of the vertical tail are chosen such that the deflections are large enough to insure interaction with the flow, and the natural frequencies are high enough to facilitate a practical computational solution. This multidisciplinary problem is solved sequentially for the fiuid flow, the elastic deformations and the grid displacements. The fluid flow is simulated by time accurately solving …

Electrostatic Positioning Of Droplets In Turbulent Flows (Lstm 375/Te/93), Nihad E. Daidzic, Adrian Melling

Aviation Department Publications

Report LSTM 375/TE/93, Lehrstuhl fuer Stroemungsmechanik Universitaet Erlangen-Nuernberg Cauerstr. 4, 8520 Erlangen Germany.

A Mathematical Model And Simulation Of Natural Circulation Loop Thermalhydraulics, Nihad E. Daidzic

Aviation Department Publications

In this work, a simple, nonlinear, seventh-order mathematical model with lumped parameter analysis of the thermalhydraulic processes in a natural circulation steam boiler loop is presented. Digital simulation has been performed in order to predict open loop transients in severe accidental situations. Parts of the circulation loop such as the steam drum, downcomer, and riser were modelled separately on the basis of the integral balance equations. The two-phase flow in the steam drum and in the riser was modelled with the aid of a homogeneous-equilibrium, two-phase flow model. Special care was taken to ensure the stability of the numerical scheme, …

Viscous Modeling And Computation Of Leading-And Trailing-Edge Vortex Cores Of Delta Wings, Balakrishnan Lakshmanan

Mechanical & Aerospace Engineering Theses & Dissertations

A Finite-Difference method is presented for calculating steady quasi-axisymmetric flow of an incompressible fluid at large Reynolds number. Approximations of the boundary-layer type are employed to reduce the Navier-Stokes equations to a pair of non-linear parabolic equations. Along with the governing equations, initial conditions are specified at some upstream cross section and boundary conditions are specified at the axis of symmetry and on the outer bounding surface.

The governing equations are replaced by a set of quasilinear finite-difference equations. The solution is obtained by a marching technique, which proceeds step-by-step in the axial direction. At each axial station, an iterative …