Simulation of flow and electrokinetic transport in a micro‐analysis‐system

The aim of this work is the investigation of the flow and electrokinetic transport in a micro‐analysis‐system. A mathematical model for electroosmotic and electrophoretic phenomena is introduced in order to perform two‐dimensional, time‐dependent Finite Element simulations for an existing device. The model includes the feature of the flow field, the mass transfer, the external applied electric field and the involved chemistry. The results of the simulation are validated against experimental data and show good agreement. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Systematic derivation of an asymptotic model for the dynamics of curved viscous fibers

This paper presents a slender body theory for the dynamics of a curved inertial viscous Newtonian fiber. Neglecting surface tension and temperature dependence, the fiber flow is modeled as a three‐dimensional free boundary value problem in terms of instationary incompressible Navier–Stokes equations. From regular asymptotic expansions in powers of the slenderness parameter, leading‐order balance laws for mass (cross‐section) and momentum are derived that combine the unrestricted motion of the fiber centerline with the inner viscous transport. The physically reasonable form of the one‐dimensional fiber model results thereby from the introduction of the intrinsic velocity that characterizes the convective terms. For the numerical investigation of the viscous, gravitational and rotational effects on the fiber dynamics, a finite volume approach on a staggered grid with implicit upwind flux discretization is applied. Copyright © 2007 John Wiley & Sons, Ltd.     

Three‐dimensional solidification with two possible crystallization states: Existence of solutions with flow in the melt

In this paper we discuss the existence of solutions of a system of nonlinear and singular partial differential equations constituting a phase field model with convection for non‐isothermal solidification/melting of certain metallic alloys in the case where two different kinds of crystallization are possible. Each one of these crystallization states is described by its own phase field, while the liquid phase is described by another one. The model also allows the occurrence of fluid flow in non‐solid regions, which are a priori unknown, and then we have a free‐boundary value problem. Thus, the model relates the evolutions of these three phase fields, the temperature of the solidification/melting process and the fluid flow in non‐solid regions. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical simulation of generalized newtonian blood flow past a couple of irregular arterial stenoses

This paper looked at the numerical investigations of the generalized Newtonian blood flow through a couple of irregular arterial stenoses. The flow is treated to be axisymmetric, with an outline of the stenoses obtained from a three dimensional casting of a mild stenosed artery, so that the flow effectively becomes two‐dimensional. The Marker and Cell (MAC) method is developed for the governing unsteady generalized Newtonian equations in staggered grid for viscous incompressible flow in the cylindrical polar co‐ordinates system. The derived pressure‐Poisson equation was solved using Successive‐Over‐Relaxation (S.O.R.) method and the pressure‐velocity correction formulae have been derived. Computations are performed for the pressure drop, the wall shear stress distribution and the separation region. The presented computations show that in comparison to the corresponding Newtonian model the generalized Newtonian fluid experiences higher pressure drop, lower peak wall shear stress and smaller separation region. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 960–981, 2011     

Numerical and asymptotic study of non‐axisymmetric magnetohydrodynamic boundary layer stagnation‐point flows

Both numerical and asymptotic analyses are performed to study the similarity solutions of three‐dimensional boundary‐layer viscous stagnation point flow in the presence of a uniform magnetic field. The three‐dimensional boundary‐layer is analyzed in a non‐axisymmetric stagnation point flow, in which the flow is developed because of influence of both applied magnetic field and external mainstream flow. Two approaches for the governing equations are employed: the Keller‐box numerical simulations solving full nonlinear coupled system and a corresponding linearized system that is obtained under a far‐field behavior and in the limit of large shear‐to‐strain‐rate parameter (λ). From these two approaches, the flow phenomena reveals a rich structure of new family of solutions for various values of the magnetic number and λ. The various results for the wall stresses and the displacement thicknesses are presented along with some velocity profiles in both directions. The analysis discovered that the flow separation occurs in the secondary flow direction in the absence of magnetic field, and the flow separation disappears when the applied magnetic field is increased. The flow field is divided into a near‐field (due to viscous forces) and far‐field (due to mainstream flows), and the velocity profiles form because of an interaction between two regions. The magnetic field plays an important role in reducing the thickness of the boundary‐layer. A physical explanation for all observed phenomena is discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Time‐dependent algorithms for viscoelastic flow: Finite element/volume schemes

Hybrid finite volume/element methods are investigated within the context of transient viscoelastic flows. A finite volume algorithm is proposed for the hyperbolic constitutive equation, of Oldroyd‐form, whereas the continuity/momentum balance is accommodated through a Taylor‐Galerkin finite element method. Various finite volume combinations are considered to derive accurate and stable implementations. Consistency of formulation is key, embracing fluctuation distribution and median‐dual‐cell constructs, within a cell‐vertex discretisation on triangles. In addition, we investigate the effect of treating the time‐term in a finite element fashion, using mass‐matrix iteration instead of the standard finite volume mass‐lumping approach. We devise an accurate transient scheme that captures the analytical solution at short and long time, both in core flow and near shear boundaries. In this respect, some difficulties are highlighted. A new method emerges, with the Low Diffusion B (LDB, with or without mass‐matrix iteration) as the optimal choice. We progress to a complex flow application and demonstrate some provocative features due to the influence of true transient boundary conditions on evolutionary flow‐structure in a 4:1 start‐up rounded‐corner contraction problem. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Gliding motion of bacteria on power‐law slime

The present investigation deals with an undulating surface model for the motility of bacteria gliding on a layer of non‐Newtonian slime. The slime being the viscoelastic material is considered as a power‐law fluid. A hydrodynamical model of motility involving an undulating cell surface which transmits stresses through a layer of exuded slime to the substratum is examined. The non‐linear differential equation resulting from the balance of momentum and mass is solved numerically by a finite difference method with an iteration technique. The manner in which the various exponent values of the power‐law flow affect the structure of the boundary layer is delineated. A comparison is made of the power‐law fluid with the Newtonian fluid. For the power‐law fluid with respect to different power‐law exponent values, shear‐thinning and shear‐thickening effects can be observed, respectively. Copyright © 2004 John Wiley & Sons, Ltd.     

Conservation Laws,Hamiltonian Structure,Modulational Instability Properties and Solitary Wave Solutions for a Higher‐Order Model Describing Nonlinear Internal Waves

Recent theoretical advances in connecting the wave‐induced mean flow with the conserved pseudomomentum per unit mass has permitted the first rational derivation of a model that describes the weakly nonlinear propagation of internal gravity plane waves in a continuously stratified fluid. Depending on the particular parameter regime examined the new model corresponds to an extended bright or dark derivative nonlinear Schrödinger equation or an extended complex‐valued modified Korteweg‐de Vries or Sasa–Satsuma equation. Mass, momentum, and energy conservation laws are derived. A noncanonical infinite‐dimensional Hamiltonian formulation of the model is introduced. The modulational stability characteristics associated with the Stokes wave solution of the model are described. The bright and dark solitary wave solutions of the model are obtained.     

Velocity and Concentration Boundary Layers on Spinning Fibers

In this paper a mathematical model for the viscose wet spinning process is presented: We consider a single fibre, which is produced by pressing a basic solution of viscose into a bath containing sulphuric acid. H₂SO₄ diffuses into the viscose solution and reacts with the natrium hydroxide so that a solidifying fibre is formed which is pulled through the bath by drives. Due to the movement of the fibre and of diffusive transport of sulphuric acid into the fibre velocity and concentration boundary layers develop. Starting from the laminar boundary layer equations we investigate the flow and concentration fields in the bath induced by the fibre utilizing the Local Non‐Similarity method. Mass transfer in the fibre is modelled by transport equations for sulphuric acid and natrium hydroxide taking into account the neutralization reaction. The model of the fibre is coupled to the bath phase model by appropriate boundary conditions for the mass flow density and the chemical potential of sulphuric acid. The non‐constant diameter of the fibre is taken into account by a perturbation approach.     

Forces acting on particles in separated wall‐bounded shear flow

Trajectories of solid particles in laminar and turbulent flow over a backward‐facing step (BFS) were numerically computed by integrating the equation of motion for particles. The various forces acting on the particles [5],[6] were calculated for a variety of flow Reynolds numbers and for different particle characteristics such as the Stokes number and the particle‐to‐fluid density ratio. The investigation was conducted for the distinct flow regimes of the BFS flow separately. Generally, the drag and gravitation were found to be the most significant forces. The lift and history force were the next most important, mostly two orders of magnitude smaller, but in some cases closing up to the other two in importance. The pressure and virtual mass effects were very small for the majority of cases. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Macro‐hybrid mixed variational two‐phase flow through fractured porous media

Macro‐hybrid mixed variational models of two‐phase flow, through fractured porous media, are analyzed at the mesoscopic and macroscopic levels. The mesoscopic models are treated in terms of nonoverlapping domain decompositions, in such a manner that the porous rock matrix system and the fracture network interact across rock–rock, rock–fracture, and fracture–fracture interfaces, with flux transmission conditions dualized. Alternatively, the models are scaled to a macroscopic level via an asymptotic process, where the width of the fractures tends to zero, and the fracture network turns out to be an interface system of one less spatial dimension, with variable high permeability. The two‐phase flow is characterized by a fractional flow dual mixed variational model. Augmented two‐field and three‐field variational reformulations are presented for regularization, internal approximations, and macro‐hybrid mixed finite element implementation. Also abstract proximal‐point penalty‐duality algorithms are derived and analyzed for parallel computing. Copyright © 2016 John Wiley & Sons, Ltd.     

An existence result for nonisothermal immiscible incompressible 2‐phase flow in heterogeneous porous media

In the present article, we study the temperature effects on two‐phase immiscible incompressible flow through a porous medium. The mathematical model is given by a coupled system of 2‐phase flow equations and an energy balance equation. The model consists of the usual equations derived from the mass conservation of both fluids along with the Darcy‐Muskat and the capillary pressure laws. The problem is written in terms of the phase formulation; ie, the saturation of one phase, the pressure of the second phase, and the temperature are primary unknowns. The major difficulties related to this model are in the nonlinear degenerate structure of the equations, as well as in the coupling in the system. Under some realistic assumptions on the data, we show the existence of weak solutions with the help of an appropriate regularization and a time discretization. We use suitable test functions to obtain a priori estimates. We prove a new compactness result to pass to the limit in nonlinear terms.     

Investigations of shear free turbulent diffusion in a rotating frame

We reconsider the problem of shear free turbulent diffusion in a rotating frame, rotating about x₁. Shear free turbulence is generated at a vibrating grid in the x₂–x₃ plane and diffuses away from the grid in x₁ direction. An important property of this flow case is that there is no mean flow‐velocity. With the help of Lie‐group methods Reynolds‐stress transport models can be analyzed for this kind of flow in a rotating frame. From the analysis it can be found, that the turbulent diffusion only influences a finite domain. Implicating this solution in the model equations shows that even fully nonlinear Reynolds‐stress transport models (non‐linear in the Reynolds‐stresses for the pressure‐strain model) are insensitive to rotation for this type of flow. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Tangential fluid flow within 3D narrow fissures: Conservative velocity fields on associated triangulations and transport processes

For a fissured medium with uncertainty in the knowledge of fractures' geometry, a conservative tangential flow field is constructed, which is consistent with the physics of stationary fluid flow in porous media and an interpolated geometry of the cracks. The flow field permits computing preferential fluid flow directions of the medium, rates of mechanical energy dissipations, and a stochastic matrix modeling stream lines and fluid mass transportation, for the analysis of solute/contaminant mass advection‐diffusion as well as drainage times.     

Parallel multigrid method for finite element simulations of complex flow problems on locally refined meshes

We present a parallel multigrid solver on locally refined meshes for solving very complex three‐dimensional flow problems. Besides describing the parallel implementation in detail, we prove the smoothing property of the suggested iteration for a simple model problem. For demonstration of the efficiency and feasibility of the solver, we show a chemically reactive flow simulation for a Methane burner using detailed chemical reaction modeling. Further, we give the results of an ocean flow simulation. All described methods are implemented in the finite element toolbox Gascoigne. Copyright © 2010 John Wiley & Sons, Ltd.     

Damping Of Pedestrian‐Induced Bridge Vibrations By Tuned Liquid Column Dampers

The excessive lateral vibrations of Londons Millenium Bridge and the Toda Park Bridge in Japan due to a large number of crossing pedestrians have raised an unexpected problem in footbridge constructions. Secondary tuned structures, like the conventional tuned mass damper (TMD) or the tuned liquid damper (TLD) were installed to the bridge in order to suppress these vibrations. In the present investigation it is proposed to apply the more efficient and more economic tuned liquid column damper (TLCD), which relies on the motion of a liquid mass in a sealed tube to counteract the external motion, while a built‐in orifice plate induces turbulent damping forces that dissipate kinetic energy. For optimal tuning of TLCDs the natural frequency and equivalent linear damping coefficient have to be chosen suitable, likewise to the conventional TMDs, as given in Den Hartog [1]. The advantages of TLCDs are: simple tuning of natural frequency and damping, low cost of design and maintenance and a simple construction. A mathematical model of a three degree‐of‐freedom (DOF) bridge coupled with an optimal tuned TLCD is derived and analyzed numerically. Furthermore, a small scale experimental model set‐up has been constructed in the laboratory of the TU‐Insitute. The experimental results are in good agreement with the theoretical predictions and indicate that TLCDs are effective damping devices for the undesired pedestrian induced footbridge vibrations. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical Investigation of the Interaction of Turbulent Wind Flow and Elastic Structures

The Interaction between wind flow and structures plays an important role in the computation of civil engineering application. In case of gravity prestressed membrane roofs, the wind lifting forces may exceed the dead load leading to high amplitude structural oscillations, which interact with the flow field. To investigate the interaction a consistent discretization method based on stabilized space‐time finite elements is applied. The flow field is modeled with the incompressible Reynolds Averaged Navier‐Stokes (RANS) equations with an anisotropic eddy‐viscosity turbulence model. The structural motion is described with the theory for geometrically nonlinear elastic deformation behavior, a strong coupling algorithm for the time‐dependent fluid‐structure interaction is implemented. Two applications show the capability of the turbulence model in representing the anisotropic turbulence structure, the differences in the flow field over a bluff body between two configurations representing a rigid and an elastic membrane roof, discusses the structural responses of the roof at a high Reynolds number. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Conservative Velocity Interpolation for PDF Methods

In many particle methods the accurate interpolation of a velocity field represented on a computational grid to arbitrary positions is crucial [2, 5]. Here, the importance of mass conservation and order of the interpolation scheme were analyzed. Initially equally distributed particles were tracked in a stationary, incompressible 2d flow field using different interpolation schemes. It could be demonstrated that especially mass conservation is of great importance, in particular in the case of complex flow patterns. The ideas presented in this paper are more general and the methods can be extended for unsteady, compressible 3d flow problems. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A general theoretical model for the magnetohydrodynamic flow of micropolar magnetic fluids. Application to Stokes flow

Many practical applications, which have an inherent interest of physical and mathematical nature, involve the hydrodynamic flow in the presence of a magnetic field. Magnetic fluids comprise a novel class of engineering materials, where the coexistence of liquid and magnetic properties provides us with the opportunity to solve problems with high mathematical and technical complexity. Here, our purpose is to examine the micropolar magnetohydrodynamic flow of magnetic fluids by considering a colloidal suspension of ferromagnetic material (usually non‐conductive) in a carrier magnetic liquid, which is in general electrically conductive. In this case, the ferromagnetic particles behave as rigid magnetic dipoles. Thus, the application of an external magnetic field, apart from the creation of an induced magnetic field of minor significance, will prevent the rotation of each particle, increasing the effective viscosity of the fluid and will cause the appearance of an additional magnetic pressure. Despite the fact that the general consideration consists of rigid particles of arbitrary shape, the assumption of spherical geometry is a very good approximation as a consequence of their small size. Our goal is to develop a general three‐dimensional theoretical model that conforms to physical reality and at the same time permits the analytical investigation of the partial differential equations, which govern the micropolar hydrodynamic flow in such magnetic liquids. Furthermore, in the aim of establishing the consistency of our proposed model with the principles of both ferrohydrodynamics and magnetohydrodynamics, we take into account both magnetization and electrical conductivity of the fluid, respectively. Under this consideration, we perform an analytical treatment of these equations in order to obtain the three‐dimensional effective viscosity and total pressure in terms of the velocity field, the total (applied and induced) magnetic field and the hydrodynamic and magnetic properties of the fluid, independently of the geometry of the flow. Moreover, we demonstrate the usefulness of our analytical approach by assuming a degenerate case of the aforementioned method, which is based on the reduction of the partial differential equations to a simpler shape that is similar to Stokes flow for the creeping motion of magnetic fluids. In view of this aim, we use the potential representation theory to construct a new complete and unique differential representation of magnetic Stokes flow, valid for non‐axisymmetric geometries, which provides the velocity and total pressure fields in terms of easy‐to‐find potentials, via an analytical fashion. Copyright © 2009 John Wiley & Sons, Ltd.