Over the years numerous computational fluid dynamics models have been developed in order to study the fluid flow in gas and induction stirred ladles. These models are used to gain insight in the industrial processes used in ladle treatment of steel. A unified model of an induction stirred Ladle in two and three dimensions is presented. Induction stirring of molten steel is a coupled multi-physics phenomena involving electromagnetic and fluid flow. Models presented in this thesis gives a more accurate description of the real stirring conditions and flow pattern, by taking into account the multi-physics behavior of the induction stirring process in an induction stirred ladle. This thesis presents a formulation of coupled electromagnetic and fluid flow equations. The coupled electromagnetic and fluid flow equations are solved using the finite element method in two and three-dimensions. The simulation model is used to predict values of steel velocities and magnetic flux density. The simulation model is also used to predict the effect of increased current density on flow velocity. Magnetic flux density values obtained from the model are verified against experimental values. / QC 20120615

The ability of Reynolds Averaged Navier Stokes (RANS) models to predict the characteristics of a non-Body of Revolution (non-BOR) Ellipsoidal model is studied to establish the feasibility of utilizing RANS as a non-BOR concept design tool. Data unable to be obtained experimentally, such as streamwise and spanwise pressure gradients and yaw turn boundary layer characteristics, are also established. A range of conditions are studied including ahead, pitched up, steady 10 and 15 degree yaw turns, and unsteady 10 and 15 degree yaw turns. Simulation results show good agreement for ahead and pitched forces and moments. Straight ahead skin friction values also showed good agreement, providing even improved agreement over an LES model which utilized wall functions. Yaw turn conditions also showed good agreement for roll angles up to 10 degrees. Steady maneuvering forces and moments showed good agreement up to 10 degrees roll and separation calculations also showed good agreement up to 10 degrees roll. Unsteady maneuvering characteristics showed mixed results, with the normal force and pitching moment trends generally agreeing with experimental data, whereas the unsteady rolling moment did not tend to follow experimental trends. Two primary conditions, the change in curvature between the mid-body and elliptical ends and the accuracy of modeling of 3D flows with RANS, are discussed as sources of discrepancies between the experimental data and steady simulations greater than 10 degrees roll and unsteady rolling simulations. / Master of Science

<P> The work presented m this thesis focused on using the hybrid Unsteady Reynolds-Averaged Navier-Stokes (URANS)/Large Eddy Simulation (LES) methodology to investigate the flow pulsation phenomenon in compound rectangular channels for isothermal flows. The specific form of the hybrid URANS/LES approach that was used is the Strelets (2001) version of the Detached Eddy Simulation (DES). It is of fundamental interest to study the problem of flow pulsations, as it is one of the most important mechanisms that directly affect the heat transfer occurring in sub-channel geometries such as those in nuclear fuel bundles. The predictions associated with the heat transfer and fluid flow in sub-channel geometry can be used to develop simplified physical models for sub-channel mixing for use in broader safety analysis codes. The primary goal of the current research work was to determine the applicability of the DES approach to predict the flow pulsations in sub-channel geometries. It was of interest to see how accurately the dynamics associated with the flow pulsations can be resolved from a spatial-temporal perspective using the specific DES model. The research work carried out for this thesis was divided into two stages. </p> <p> In the first stage of the research work, effort was concentrated to primarily understand the field of sub-channel flow pulsations and its implications from both an experimental and numerical point of view. It was noted that unsteady turbulence modeling approaches have great potential in providing insights into the fundamentals of sub-channel flow pulsations. It was proposed that for this thesis work, the Shear Stress Transport (SST) based DES model be used to understand the dynamics associated with sub-channel flow pulsations. To the author's knowledge the DES-SST based turbulence model has never been used for resolving the effects of sub-channel flow pulsations. Next, the hybrid URANS/LES turbulence modeling technique was reviewed in great detail to understand the philosophy of the hybrid URANS/LES technique and its ability to resolve fundamental flows of interest. Effort was directed to understand the switching mechanism (which blends the URANS region with the LES region) in the DES-SST model for fully wall bounded turbulent flows without boundary layer separation. To the author's knowledge, the DES-SST model has never been used on a fully wall bounded turbulent flow problem without boundary layer separation. Thus, the DES-SST model was first completely validated for a fully developed turbulent channel flow problem without boundary layer separation. </p> <p> In the second stage of the research work, the DES-SST model was used to study the flow pulsation phenomena on two rectangular sub-channels connected by a gap, on which extensive experiments were conducted by Meyer and Rehme (1994). It was found that the DES-SST model was successful in resolving significant portion of the flow field in the vicinity of the gap region. The span-wise velocity contours, velocity vector plots, and time traces of the velocity components showed the expected cross flow mixing between the sub-channels through the gap. The predicted turbulent kinetic energy showed two clear peaks at the edges of the gap. The dynamics of the flow pulsations were quantitatively described through temporal auto-correlations, spatial cross-correlations and power spectral functions. The numerical predictions were in general agreement with the experiments in terms of the quantitative aspects. From an instantaneous time scale point of view, the DES-SST model was able to identify different flow mixing patterns. The pulsating flow is basically an effect of the variation of the pressure field which is a response to the instability causing the fluid flow pulsations. Coherent structures were identified in the flow field to be comprised of eddies, shear zones and streams. Eddy structures with high vorticity and low pressure cores were found to exist near the vicinity of the gap edge region. A three dimensional vorticity field was identified and found to exist near the gap edge region. The instability mechanism and the probable cause behind the quasi-periodic fluid flow pulsations was identified and related to the inflectional stream-wise velocity profile. Simulations were also performed with two different channel lengths in comparison to the reference channel length. Different channel length studies showed similar statistical description of the flow field. However, frequency independent results were not obtained. In general, simulations performed using the DES-SST model were successful in capturing the effects of the fluid flow pulsations. This modeling technique has great potential to be used for actual rod bundle configurations. </p> / Thesis / Doctor of Philosophy (PhD)

We consider three different phenomena governing the fluid flow in the presence of charged particles: electroconvection in fluids, electroconvection in porous media, and electrodiffusion. Electroconvecton in fluids is mathematically represented by a nonlinear drift-diffusion partial differential equation describing the time evolution of a surface charge density in a two-dimensional incompressible fluid. The velocity of the fluid evolves according to Navier-Stokes equations forced nonlinearly by the electrical forces due to the presence of the charge density. The resulting model is reminiscent of the quasi-geostrophic equation, where the main difference resides in the dependence of the drift velocity on the charge density. When the fluid flows through a porous medium, the velocity and the electrical forces are related according to Darcy’s law, which yields a challenging doubly nonlinear and doubly nonlocal model describing electroconvection in porous media. A different type of particle-fluid interaction, called electrodiffusion, is also considered. This latter phenomenon is described by nonlinearly advected and nonlinearly forced continuity equations tracking the time evolution of the concentrations of many ionic species having different valences and diffusivities and interacting with an incompressible fluid. This work is based on [1, 2, 3] and addresses the global well-posedness, long-time dynamics, and other features associated with the aforementioned three models. REFERENCES:[1] E. Abdo, M. Ignatova, Long time dynamics of a model of electroconvection, Trans. Amer. Math. Soc. 374 (2021), 5849–5875. [2] E. Abdo, M. Ignatova, Long Time Finite Dimensionality in Charged Fluids, Nonlinearity 34 (9) (2021), 6173–6209. [3] E. Abdo, M. Ignatova, On Electroconvection in Porous Media, to appear in Indiana University Mathematics Journal (2023). / Mathematics

No description available.

No description available.

We analyze the mathematical models of two classes of physical phenomena. The first class of phenomena we consider is the interaction between one or more insulating rigid bodies and an electrically conducting fluid, inside of which the bodies are contained, as well as the electromagnetic fields trespassing both of the materials. We take into account both the cases of incompressible and compressible fluids. In both cases our main result yields the existence of weak solutions to the associated system of partial differential equations, respectively. The proofs of these results are built upon hybrid discrete-continuous approximation schemes: Parts of the systems are discretized with respect to time in order to deal with the solution-dependent test functions in the induction equation. The remaining parts are treated as continuous equations on the small intervals between consecutive discrete time points, allowing us to employ techniques which do not transfer to the discretized setting. Moreover, the solution-dependent test functions in the momentum equation are handled via the use of classical penalization methods. The second class of phenomena we consider is the evolution of a magnetoelastic material. Here too, our main result proves the existence of weak solutions to the corresponding system of partial differential equations. Its proof is based on De Giorgi's minimizing movements method, in which the system is discretized in time and, at each discrete time point, a minimization problem is solved, the associated Euler-Lagrange equations of which constitute a suitable approximation of the original equation of motion and magnetic force balance. The construction of such a minimization problem is made possible by the realization that, already on the continuous level, both of these equations can be written in terms of the same energy and dissipation potentials. The functional for the discrete minimization problem can then be constructed on the basis of these potentials. / Wir analysieren die mathematischen Modelle von zwei Arten physikalischer Phänomene. Die erste Art von Phänomenen, die wir betrachten, ist die Wechselwirkung zwischen einem oder mehreren isolierenden starren Körpern und einem elektrisch leitenden Fluid, das die Körper umgibt, sowie den elektromagnetischen Feldern in beiden Materialien. Wir untersuchen sowohl den Fall inkompressibler als auch kompressibler Fluide. In beiden Fällen liefert unser Hauptresultat die Existenz von schwachen Lösungen für das zugehörige System partieller Differentialgleichungen. Die Beweise dieser Resultate beruhen auf hybriden diskret-kontinuierlichen Approximationsmethoden: Teile der Systeme werden in der Zeit diskretisiert, um das Problem der lösungsabhängigen Testfunktionen in der Induktionsgleichung zu bewältigen. Die verbleibenden Gleichungen werden als kontinuierliche Gleichungen auf den kleinen Intervallen zwischen aufeinanderfolgenden diskreten Zeitpunkten behandelt, sodass wir Techniken anwenden können, die sich nicht auf das diskretisierte System übertragen lassen. Darüber hinaus wird das Problem der lösungsabhängigen Testfunktionen in der Impulsgleichung durch die Verwendung klassischer Penalisierungsmethoden gelöst. Die zweite Art von Phänomenen, die wir betrachten, ist die Entwicklung eines magnetoelastischen Materials. Auch hier beweist unser Hauptresultat die Existenz schwacher Lösungen für das zugehörige System partieller Differentialgleichungen. Der Beweis basiert auf der Methode von De Giorgi, bei der das System in der Zeit diskretisiert und in jedem diskreten Zeitpunkt ein Minimierungsproblem gelöst wird, dessen zugehörige Euler-Lagrange-Gleichungen eine geeignete Approximation an die ursprüngliche Bewegungsgleichung und mikromagnetische Gleichung darstellen. Die Konstruktion eines solchen Minimierungsproblems wird durch die Erkenntnis ermöglicht, dass diese beiden Gleichungen bereits im kontinuierlichen System mithilfe derselben Energie- und Dissipationspotenziale ausgedrückt werden können. Das Funktional für das diskrete Minimierungsproblem kann dann auf Grundlage dieser Potenziale konstruiert werden.

This work is concerned with the development of computational methods for approximating sensitivities of solutions to boundary value problems. We focus on the continuous sensitivity equation method and investigate the application of adaptive meshing and smoothing projection techniques to enhance the basic scheme. The fundamental ideas are first developed for a one dimensional problem and then extended to 2-D flow problems governed by the incompressible Navier-Stokes equations. Numerical experiments are conducted to test the algorithms and to investigate the benefits of adaptivity and smoothing. / Ph. D.

We provide an algorithm and analysis of a high order projection scheme for time integration of the incompressible Navier-Stokes equations (NSE). The method is based on a projection onto the subspace of divergence-free (incompressible) functions interleaved with a Krylov-based exponential time integration (KBEI). These time integration methods provide a high order accurate, stable approach with many of the advantages of explicit methods, and can reduce the computational resources over conventional methods. The method is scalable in the sense that the computational costs grow linearly with problem size. Exponential integrators, used typically to solve systems of ODEs, utilize matrix vector products of the exponential of the Jacobian on a vector. For large systems, this product can be approximated efficiently by Krylov subspace methods. However, in contrast to explicit methods, KBEIs are not restricted by the time step. While implicit methods require a solution of a linear system with the Jacobian, KBEIs only require matrix vector products of the Jacobian. Furthermore, these methods are based on linearization, so there is no non-linear system solve at each time step. Differential-algebraic equations (DAEs) are ordinary differential equations (ODEs) subject to algebraic constraints. The discretized NSE constitute a system of DAEs, where the incompressibility condition is the algebraic constraint. Exponential integrators can be extended to DAEs with linear constraints imposed via a projection onto the constraint manifold. This results in a projected ODE that is integrated by a KBEI. In this approach, the Krylov subspace satisfies the constraint, hence the solution at the advanced time step automatically satisfies the constraint as well. For the NSE, the projection onto the constraint is typically achieved by a projection induced by the L2 inner product. We examine this L2 projection and an H1 projection induced by the H1 semi-inner product. The H1 projection has an advantage over the L2 projection in that it retains tangential Dirichlet boundary conditions for the flow. Both the H1 and L2 projections are solutions to saddle point problems that are efficiently solved by a preconditioned Uzawa algorithm. / Ph. D.

This thesis describes a simple efficient and robust numerical technique for solving two-dimensional incompressible laminar steady flows at moderate-to-high Reynolds numbers. The method uses an incremental multigrid method and an extrapolation procedure based on minimum residual concepts to accelerate the convergence rate of a robust block-line-Gauss-Seidel solver for the vorticity-stream function equations. Results are presented for the driven cavity flow problem using uniform and nonuniform grids and for the flow past a backward facing step in a channel. / M.S.