DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

COACHES, CLIMATES, “FIELD” GOALS, AND EFFICACY: A “DE-CONSTRUCTION” OF THE MASTERY-APPROACH TO COACHING AND EXAMINATION OF RELATIONSHIPS TO PSYCHOSOCIAL OUTCOMES IN A YOUTH FOOTBALL PLAYER DEVELOPMENT PROGRAM.

In support of the achievement goal theory (AGT), empirical research has demonstrated psychosocial benefits of the mastery-oriented learning climate. In this study, we examined the effects of perceived coaching behaviors on various indicators of psychosocial well-being (competitive anxiety, self-esteem, perceived competence, enjoyment, and future intentions for participation), as mediated by perceptions of the coach-initiated motivational climate, achievement goal orientations and perceptions of sport-specific skills efficacy. Using a pre-post test design, 1,464 boys, ages 10-15 (M = 12.84 years, SD = 1.44), who participated in a series of 12 football skills clinics were surveyed from various locations across the United States. Using structural equation modeling (SEM) path analysis and hierarchical regression analysis, the cumulative direct and indirect effects of the perceived coaching behaviors on the psychosocial variables at post-test were parsed out to determine what types of coaching behaviors are more conducive to the positive psychosocial development of youth athletes. The study demonstrated that how coaching behaviors are perceived impacts the athletes’ perceptions of the motivational climate and achievement goal orientations, as well as self-efficacy beliefs. These effects in turn affect the athletes’ self-esteem, general competence, sport-specific competence, competitive anxiety, enjoyment, and intentions to remain involved in the sport. The findings also clarify how young boys internalize and interpret coaches’ messages through modification of achievement goal orientations and sport-specific efficacy beliefs.

A Social Cognitive Approach to Coping with Acculturative Stress in International Students

The present study employed a cross-sectional design to test a model of coping with acculturative stress in an international student sample. Drawing from Lent’s (2004) social cognitive model of restorative well-being, several direct and mediated paths were hypothesized to predict (negatively) acculturative stress and (positively) life satisfaction. Behavioral acculturation and behavioral enculturation (Kim & Omizo, 2006) were also examined as predictors of coping with acculturative stress among international students. Using a self-report survey, participants’ ratings of acculturative stress, life satisfaction, social support, behavioral acculturation, behavioral enculturation, and coping self-efficacy were assessed. The results revealed that the variables of the model explained 16% of the variance in acculturative stress and 27% of the variance in life satisfaction. A final model, including the use of modification indices, provided good fit to the data. Findings also suggested that coping self-efficacy was a direct predictor of acculturative stress, and that behavioral acculturation and coping self-efficacy were direct predictors of students’ life satisfaction. Limitations, future research, and practical implications are discussed.

A Visual Analytics Approach to Comparing Cohorts of Event Sequences

Sequences of timestamped events are currently being generated across nearly every domain of data analytics, from e-commerce web logging to electronic health records used by doctors and medical researchers. Every day, this data type is reviewed by humans who apply statistical tests, hoping to learn everything they can about how these processes work, why they break, and how they can be improved upon. To further uncover how these processes work the way they do, researchers often compare two groups, or cohorts, of event sequences to find the differences and similarities between outcomes and processes. With temporal event sequence data, this task is complex because of the variety of ways single events and sequences of events can differ between the two cohorts of records: the structure of the event sequences (e.g., event order, co-occurring events, or frequencies of events), the attributes about the events and records (e.g., gender of a patient), or metrics about the timestamps themselves (e.g., duration of an event). Running statistical tests to cover all these cases and determining which results are significant becomes cumbersome. Current visual analytics tools for comparing groups of event sequences emphasize a purely statistical or purely visual approach for comparison. Visual analytics tools leverage humans' ability to easily see patterns and anomalies that they were not expecting, but is limited by uncertainty in findings. Statistical tools emphasize finding significant differences in the data, but often requires researchers have a concrete question and doesn't facilitate more general exploration of the data. Combining visual analytics tools with statistical methods leverages the benefits of both approaches for quicker and easier insight discovery. Integrating statistics into a visualization tool presents many challenges on the frontend (e.g., displaying the results of many different metrics concisely) and in the backend (e.g., scalability challenges with running various metrics on multi-dimensional data at once). I begin by exploring the problem of comparing cohorts of event sequences and understanding the questions that analysts commonly ask in this task. From there, I demonstrate that combining automated statistics with an interactive user interface amplifies the benefits of both types of tools, thereby enabling analysts to conduct quicker and easier data exploration, hypothesis generation, and insight discovery. The direct contributions of this dissertation are: (1) a taxonomy of metrics for comparing cohorts of temporal event sequences, (2) a statistical framework for exploratory data analysis with a method I refer to as high-volume hypothesis testing (HVHT), (3) a family of visualizations and guidelines for interaction techniques that are useful for understanding and parsing the results, and (4) a user study, five long-term case studies, and five short-term case studies which demonstrate the utility and impact of these methods in various domains: four in the medical domain, one in web log analysis, two in education, and one each in social networks, sports analytics, and security. My dissertation contributes an understanding of how cohorts of temporal event sequences are commonly compared and the difficulties associated with applying and parsing the results of these metrics. It also contributes a set of visualizations, algorithms, and design guidelines for balancing automated statistics with user-driven analysis to guide users to significant, distinguishing features between cohorts. This work opens avenues for future research in comparing two or more groups of temporal event sequences, opening traditional machine learning and data mining techniques to user interaction, and extending the principles found in this dissertation to data types beyond temporal event sequences.

Effects of Drive Speed Modulation on Dynamics of Slender Rotating Structures

Slender rotating structures are used in many mechanical systems. These structures can suffer from undesired vibrations that can affect the components and safety of a system. Furthermore, since some these structures can operate in a harsh environment, installation and operation of sensors that are needed for closed-loop and collocated control schemes may not be feasible. Hence, the need for an open-loop non-collocated scheme for control of the dynamics of these structures. In this work, the effects of drive speed modulation on the dynamics of slender rotating structures are studied. Slender rotating structures are a type of mechanical rotating structures, whose length to diameter ratio is large. For these structures, the torsion mode natural frequencies can be low. In particular, for isotropic structures, the first few torsion mode frequencies can be of the same order as the first few bending mode frequencies. These situations can be conducive for energy transfer amongst bending and torsion modes. Scenarios with torsional vibrations experienced by rotating structures with continuous rotor-stator contact occur in many rotating mechanical systems. Drill strings used in the oil and gas industry are an example of rotating structures whose torsional vibrations can be deleterious to the components of the drilling system. As a novel approach to mitigate undesired vibrations, the effects of adding a sinusoidal excitation to the rotation speed of a drill string are studied. A portion of the drill string located within a borewell is considered and this rotating structure has been modeled as an extended Jeffcott rotor and a sinusoidal excitation has been added to the drive speed of the rotor. After constructing a three-degree-of-freedom model to capture lateral and torsional motions, the equations of motions are reduced to a single differential equation governing torsional vibrations during continuous stator contact. An approximate solution has been obtained by making use of the Method of Direct Partition of Motions with the governing torsional equation of motion. The results showed that for a rotor undergoing forward or backward whirling, the addition of sinusoidal excitation to the drive speed can cause an increase in the equivalent torsional stiffness, smooth the discontinuous friction force at contact, and reduce the regions of negative slope in the friction coefficient variation with respect to speed. Experiments with a scaled drill string apparatus have also been conducted and the experimental results show good agreement with the numerical results obtained from the developed models. These findings suggest that the extended Jeffcott rotordynamics model can be useful for studies of rotor dynamics in situations with continuous rotor-stator contact. Furthermore, the results obtained suggest that the drive speed modulation scheme can have value for attenuating drill-string vibrations.