Search results for: non Newtonian
91 Coding Considerations for Standalone Molecular Dynamics Simulations of Atomistic Structures
Authors: R. O. Ocaya, J. J. Terblans
Abstract:
The laws of Newtonian mechanics allow ab-initio molecular dynamics to model and simulate particle trajectories in material science by defining a differentiable potential function. This paper discusses some considerations for the coding of ab-initio programs for simulation on a standalone computer and illustrates the approach by C language codes in the context of embedded metallic atoms in the face-centred cubic structure. The algorithms use velocity-time integration to determine particle parameter evolution for up to several thousands of particles in a thermodynamical ensemble. Such functions are reusable and can be placed in a redistributable header library file. While there are both commercial and free packages available, their heuristic nature prevents dissection. In addition, developing own codes has the obvious advantage of teaching techniques applicable to new problems.Keywords: C language, molecular dynamics, simulation, embedded atom method
Procedia PDF Downloads 27690 Rheological Evaluation of a Mucoadhesive Precursor of Based-Poloxamer 407 or Polyethylenimine Liquid Crystal System for Buccal Administration
Authors: Jéssica Bernegossi, Lívia Nordi Dovigo, Marlus Chorilli
Abstract:
Mucoadhesive liquid crystalline systems are emerging how delivery systems for oral cavity. These systems are interesting since they facilitate the targeting of medicines and change the release enabling a reduction in the number of applications made by the patient. The buccal mucosa is permeable besides present a great blood supply and absence of first pass metabolism, it is a good route of administration. It was developed two systems liquid crystals utilizing as surfactant the ethyl alcohol ethoxylated and propoxylated (30%) as oil phase the oleic acid (60%), and the aqueous phase (10%) dispersion of polymer polyethylenimine (0.5%) or dispersion of polymer poloxamer 407 (16%), with the intention of applying the buccal mucosa. Initially, was performed for characterization of systems the conference by polarized light microscopy and rheological analysis. For the preparation of the systems the components described was added above in glass vials and shaken. Then, 30 and 100% artificial saliva were added to each prepared formulation so as to simulate the environment of the oral cavity. For the verification of the system structure, aliquots of the formulations were observed in glass slide and covered with a coverslip, examined in polarized light microscope (PLM) Axioskop - Zeizz® in 40x magnifier. The formulations were also evaluated for their rheological profile Rheometer TA Instruments®, which were obtained rheograms the selected systems employing fluency mode (flow) in temperature of 37ºC (98.6ºF). In PLM, it was observed that in formulations containing polyethylenimine and poloxamer 407 without the addition of artificial saliva was observed dark-field being indicative of microemulsion, this was also observed with the formulation that was increased with 30% of the artificial saliva. In the formulation that was increased with 100% simulated saliva was shown to be a system structure since it presented anisotropy with the presence of striae being indicative of hexagonal liquid crystalline mesophase system. Upon observation of rheograms, both systems without the addition of artificial saliva showed a Newtonian profile, after addition of 30% artificial saliva have been given a non-Newtonian behavior of the pseudoplastic-thixotropic type and after adding 100% of the saliva artificial proved plastic-thixotropic. Furthermore, it is clearly seen that the formulations containing poloxamer 407 have significantly larger (15-800 Pa) shear stress compared to those containing polyethyleneimine (5-50 Pa), indicating a greater plasticity of these. Thus, it is possible to observe that the addition of saliva was of interest to the system structure, starting from a microemulsion for a liquid crystal system, thereby also changing thereby its rheological behavior. The systems have promising characteristics as controlled release systems to the oral cavity, as it features good fluidity during its possible application and greater structuring of the system when it comes into contact with environmental saliva.Keywords: liquid crystal system, poloxamer 407, polyethylenimine, rheology
Procedia PDF Downloads 42789 Stabilizing Effect of Magnetic Field in a Thermally Modulated Porous Layer
Authors: M. Meenasaranya, S. Saravanan
Abstract:
Nonlinear stability analysis is carried out to determine the effect of surface temperature modulation in an infinite horizontal porous layer heated from below. The layer is saturated by an electrically conducting, viscous, incompressible and Newtonian fluid. The Brinkman model is used for momentum equation, and the Boussinesq approximation is invoked. The system is assumed to be bounded by rigid boundaries. The energy theory is implemented to find the global exponential stability region of the considered system. The results are analysed for arbitrary values of modulation frequency and amplitude. The existence of subcritical instability region is confirmed by comparing the obtained result with the known linear result. The vertical magnetic field is found to stabilize the system.Keywords: Brinkman model, energy method, magnetic field, surface temperature modulation
Procedia PDF Downloads 36888 Yarkovsky Effect on the Orbital Dynamics of the Asteroid (101955) Bennu
Authors: Sanjay Narayan Deo, Badam Singh Kushvah
Abstract:
Bennu(101955) is a half kilometer potentially hazardous near-Earth asteroid. We analyze the influence of Yarkovsky effect and relativistic effect of the Sun on the motion of the asteroid Bennu. The transverse model is used to compute Yarkovsky force on asteroid Bennu. Our dynamical model includes Newtonian perturbations of eight planets, the Moon, the Sun and three massive asteroid (1Ceres, 2Palas and 4Vesta). We showed the variation in orbital elements of nominal orbit of the asteroid. In the presence of Yarkovsky effect, the Semi-major axis of the orbit of the asteroid is decreases by 350 m over one period of orbital motion. The magnitude of Yarkovsky force is computed. We find that maximum magnitude of Yarkovsky force is 0.09 N at the perihelion . We also found that the magnitude of the Sun relativity effect is greater than the Yarkovsky effect on the motion the asteroid Bennu.Keywords: Bennu, orbital elements, relativistic effect, Yarkovsky effect
Procedia PDF Downloads 27087 Temporal and Spatio-Temporal Stability Analyses in Mixed Convection of a Viscoelastic Fluid in a Porous Medium
Authors: P. Naderi, M. N. Ouarzazi, S. C. Hirata, H. Ben Hamed, H. Beji
Abstract:
The stability of mixed convection in a Newtonian fluid medium heated from below and cooled from above, also known as the Poiseuille-Rayleigh-Bénard problem, has been extensively investigated in the past decades. To our knowledge, mixed convection in porous media has received much less attention in the published literature. The present paper extends the mixed convection problem in porous media for the case of a viscoelastic fluid flow owing to its numerous environmental and industrial applications such as the extrusion of polymer fluids, solidification of liquid crystals, suspension solutions and petroleum activities. Without a superimposed through-flow, the natural convection problem of a viscoelastic fluid in a saturated porous medium has already been treated. The effects of the viscoelastic properties of the fluid on the linear and nonlinear dynamics of the thermoconvective instabilities have also been treated in this work. Consequently, the elasticity of the fluid can lead either to a Hopf bifurcation, giving rise to oscillatory structures in the strongly elastic regime, or to a stationary bifurcation in the weakly elastic regime. The objective of this work is to examine the influence of the main horizontal flow on the linear and characteristics of these two types of instabilities. Under the Boussinesq approximation and Darcy's law extended to a viscoelastic fluid, a temporal stability approach shows that the conditions for the appearance of longitudinal rolls are identical to those found in the absence of through-flow. For the general three-dimensional (3D) perturbations, a Squire transformation allows the deduction of the complex frequencies associated with the 3D problem using those obtained by solving the two-dimensional one. The numerical resolution of the eigenvalue problem concludes that the through-flow has a destabilizing effect and selects a convective configuration organized in purely transversal rolls which oscillate in time and propagate in the direction of the main flow. In addition, by using the mathematical formalism of absolute and convective instabilities, we study the nature of unstable three-dimensional disturbances. It is shown that for a non-vanishing through-flow, general three-dimensional instabilities are convectively unstable which means that in the absence of a continuous noise source these instabilities are drifted outside the porous medium, and no long-term pattern is observed. In contrast, purely transversal rolls may exhibit a transition to absolute instability regime and therefore affect the porous medium everywhere including in the absence of a noise source. The absolute instability threshold, the frequency and the wave number associated with purely transversal rolls are determined as a function of the Péclet number and the viscoelastic parameters. Results are discussed and compared to those obtained from laboratory experiments in the case of Newtonian fluids.Keywords: instability, mixed convection, porous media, and viscoelastic fluid
Procedia PDF Downloads 31486 An Implicit Methodology for the Numerical Modeling of Locally Inextensible Membranes
Authors: Aymen Laadhari
Abstract:
We present in this paper a fully implicit finite element method tailored for the numerical modeling of inextensible fluidic membranes in a surrounding Newtonian fluid. We consider a highly simplified version of the Canham-Helfrich model for phospholipid membranes, in which the bending force and spontaneous curvature are disregarded. The coupled problem is formulated in a fully Eulerian framework and the membrane motion is tracked using the level set method. The resulting nonlinear problem is solved by a Newton-Raphson strategy, featuring a quadratic convergence behavior. A monolithic solver is implemented, and we report several numerical experiments aimed at model validation and illustrating the accuracy of the proposed method. We show that stability is maintained for significantly larger time steps with respect to an explicit decoupling method.Keywords: finite element method, level set, Newton, membrane
Procedia PDF Downloads 30585 Pressure Distribution, Load Capacity, and Thermal Effect with Generalized Maxwell Model in Journal Bearing Lubrication
Authors: M. Guemmadi, A. Ouibrahim
Abstract:
This numerical investigation aims to evaluate how a viscoelastic lubricant described by a generalized Maxwell model, affects the pressure distribution, the load capacity and thermal effect in a journal bearing lubrication. We use for the purpose the CFD package software completed by adapted user define functions (UDFs) to solve the coupled equations of momentum, of energy and of the viscoelastic model (generalized Maxwell model). Two parameters, viscosity and relaxation time are involved to show how viscoelasticity substantially affect the pressure distribution, the load capacity and the thermal transfer by comparison to Newtonian lubricant. These results were also compared with the available published results.Keywords: journal bearing, lubrication, Maxwell model, viscoelastic fluids, computational modelling, load capacity
Procedia PDF Downloads 51884 Numerical Investigation of Flow Behaviour Across a Trapezoidal Bluff Body at Low Reynolds Number
Authors: Zaaraoui Abdelkader, Kerfah Rabeh, Noura Belkheir, Matene Elhacene
Abstract:
The trapezoidal bluff body is a typical configuration of vortex shedding bodies. The aim of this work is to study flow behaviour over a trapezoidal cylinder at low Reynolds number. The geometry was constructed from a prototype device for measuring the volumetric flow-rate by counting vortices. Simulations were run for this geometry under steady and unsteady flow conditions using finite volume discretization. Laminar flow was investigated in this model with rigid walls and homogeneous incompressible Newtonian fluid. Calculations were performed for Reynolds number range 5 ≤ Re ≤ 180 and several flow parameters were documented. The present computations are in good agreement with the experimental observations and the numerical calculations by several investigators.Keywords: bluff body, confined flow, numerical calculations, steady and unsteady flow, vortex shedding flow meter
Procedia PDF Downloads 25783 Magneto-Solutal Convection in Newtonian Fluid Layer with Modulated Gravity
Authors: Om Prakash Keshri, Anand Kumar, Vinod K. Gupta
Abstract:
In the present study, the effect of gravity modulation on the onset of convection in viscous fluid layer under the influence of induced magnetic field, salted from above on the boundaries, has been investigated. Linear and nonlinear stability analysis has been performed. A linear stability analysis is performed to show that the gravity modulation can significantly affect the stability limits of the system. A method based on small amplitude of the modulation is used to compute the critical value of Rayleigh number and wave number. The effect of Smith number, salute Rayleigh number and magnetic Prandtl number on the stability of the system is investigated.Keywords: viscous fluid, induced magnetic field, gravity modulation, salute convection
Procedia PDF Downloads 16482 Fiber Orientation Measurements in Reinforced Thermoplastics
Authors: Ihsane Modhaffar
Abstract:
Fiber orientation is essential for the physical properties of composite materials. The theoretical parameters of a given reinforcement are usually known and widely used to predict the behavior of the material. In this work, we propose an image processing approach to estimate true principal directions and fiber orientation during injection molding processes of short fiber reinforced thermoplastics. Generally, a group of fibers are described in terms of probability distribution function or orientation tensor. Numerical techniques for the prediction of fiber orientation are also considered for concentrated situations. The flow was considered to be incompressible, and behave as Newtonian fluid containing suspensions of short-fibers. The governing equations, of this problem are: the continuity, the momentum and the energy. The obtained results were compared to available experimental findings. A good agreement between the numerical results and the experimental data was achieved.Keywords: injection, composites, short-fiber reinforced thermoplastics, fiber orientation, incompressible fluid, numerical simulation
Procedia PDF Downloads 50281 Forced Convection Boundary Layer Flow of a Casson Fluid over a Moving Permeable Flat Plate beneath a Uniform Free Stream
Authors: N. M. Arifin, S. P. M. Isa, R. Nazar, N. Bachok, F. M. Ali, I. Pop
Abstract:
In this paper, the steady forced convection boundary layer flow of a Casson fluid past a moving permeable semi-infinite flat plate beneath a uniform free stream is investigated. The mathematical problem reduces to a pair of noncoupled ordinary differential equations by similarity transformation, which is then solved numerically using the shooting method. Both the cases when the plate moves into or out of the origin are considered. Effects of the non-Newtonian (Casson) parameter, moving parameter, suction or injection parameter and Eckert number on the flow and heat transfer characteristics are thoroughly examined. Dual solutions are found to exist for each value of the governing parameters.Keywords: forced convection, Casson fluids, moving flat plate, boundary layer
Procedia PDF Downloads 43380 Vortex Control by a Downstream Splitter Plate in Psudoplastic Fluid Flow
Authors: Sudipto Sarkar, Anamika Paul
Abstract:
Pseudoplastic (n<1, n is the power index) fluids have great importance in food, pharmaceutical and chemical process industries which require a lot of attention. Unfortunately, due to its complex flow behavior inadequate research works can be found even in laminar flow regime. A practical problem is solved in the present research work by numerical simulation where we tried to control the vortex shedding from a square cylinder using a horizontal splitter plate placed at the downstream flow region. The position of the plate is at the centerline of the cylinder with varying distance from the cylinder to calculate the critical gap-ratio. If the plate is placed inside this critical gap, the vortex shedding from the cylinder suppressed completely. The Reynolds number considered here is in unsteady laminar vortex shedding regime, Re = 100 (Re = U∞a/ν, where U∞ is the free-stream velocity of the flow, a is the side of the cylinder and ν is the maximum value of kinematic viscosity of the fluid). Flow behavior has been studied for three different gap-ratios (G/a = 2, 2.25 and 2.5, where G is the gap between cylinder and plate) and for a fluid with three different flow behavior indices (n =1, 0.8 and 0.5). The flow domain is constructed using Gambit 2.2.30 and this software is also used to generate the mesh and to impose the boundary conditions. For G/a = 2, the domain size is considered as 37.5a × 16a with 316 × 208 grid points in the streamwise and flow-normal directions respectively after a thorough grid independent study. Fine and equal grid spacing is used close to the geometry to capture the vortices shed from the cylinder and the boundary layer developed over the flat plate. Away from the geometry meshes are unequal in size and stretched out. For other gap-ratios, proportionate domain size and total grid points are used with similar kind of mesh distribution. Velocity inlet (u = U∞), pressure outlet (Neumann condition), symmetry (free-slip boundary condition) at upper and lower domain boundary conditions are used for the simulation. Wall boundary condition (u = v = 0) is considered both on the cylinder and the splitter plate surfaces. Discretized forms of fully conservative 2-D unsteady Navier Stokes equations are then solved by Ansys Fluent 14.5. SIMPLE algorithm written in finite volume method is selected for this purpose which is a default solver inculcate in Fluent. The results obtained for Newtonian fluid flow agree well with previous works supporting Fluent’s usefulness in academic research. A thorough analysis of instantaneous and time-averaged flow fields are depicted both for Newtonian and pseudoplastic fluid flow. It has been observed that as the value of n reduces the stretching of shear layers also reduce and these layers try to roll up before the plate. For flow with high pseudoplasticity (n = 0.5) the nature of vortex shedding changes and the value of critical gap-ratio reduces. These are the remarkable findings for laminar periodic vortex shedding regime in pseudoplastic flow environment.Keywords: CFD, pseudoplastic fluid flow, wake-boundary layer interactions, critical gap-ratio
Procedia PDF Downloads 7879 Dam Break Model Using Navier-Stokes Equation
Authors: Alireza Lohrasbi, Alireza Lavaei, Mohammadali M. Shahlaei
Abstract:
The liquid flow and the free surface shape during the initial stage of dam breaking are investigated. A numerical scheme is developed to predict the wave of an unsteady, incompressible viscous flow with free surface. The method involves a two dimensional finite element (2D), in a vertical plan. The Naiver-Stokes equations for conservation of momentum and mass for Newtonian fluids, continuity equation, and full nonlinear kinematic free-surface equation were used as the governing equations. The mapping developed to solve highly deformed free surface problems common in waves formed during wave propagation, transforms the run up model from the physical domain to a computational domain with Arbitrary Lagrangian Eulerian (ALE) finite element modeling technique.Keywords: dam break, Naiver-Stokes equations, free-surface flows, Arbitrary Lagrangian-Eulerian
Procedia PDF Downloads 30078 Chip Less Microfluidic Device for High Throughput Liver Spheroid Generation
Authors: Sourita Ghosh, Falguni Pati, Suhanya Duraiswamy
Abstract:
Spheroid, a simple three-dimensional cellular aggregate, allows us to simulate the in-vivo complexity of cellular signaling and interactions in greater detail than traditional 2D cell culture. It can be used as an in-vitro model for drug toxicity testing, tumor modeling and many other such applications specifically for cancer. Our work is focused on the development of an affordable, user-friendly, robust, reproducible, high throughput microfluidic device for water in oil droplet production, which can, in turn, be used for spheroids manufacturing. Here, we have investigated the droplet breakup between two non-Newtonian fluids, viz. silicone oil and decellularized liver matrix, which acts as our extra cellular matrix (ECM) for spheroids formation. We performed some biochemical assays to characterize the liver ECM, as well as rheological studies on our two fluids and observed a critical dependence of capillary number (Ca) on droplet breakup and homogeneous drop formationKeywords: chip less, droplets, extracellular matrix, liver spheroid
Procedia PDF Downloads 6077 Numerical Investigation of Pressure and Velocity Field Contours of Dynamics of Drop Formation
Authors: Pardeep Bishnoi, Mayank Srivastava, Mrityunjay Kumar Sinha
Abstract:
This article represents the numerical investigation of the pressure and velocity field variation of the dynamics of pendant drop formation through a capillary tube. Numerical simulations are executed using volume of fluid (VOF) method in the computational fluid dynamics (CFD). In this problem, Non Newtonian fluid is considered as dispersed fluid whereas air is considered as a continuous fluid. Pressure contours at various time steps expose that pressure varies nearly hydrostatically at each step of the dynamics of drop formation. A result also shows the pressure variation of the liquid droplet during free fall in the computational domain. The evacuation of the fluid from the necking region is also shown by the contour of the velocity field. The role of surface tension in the Pressure contour of the dynamics of drop formation is also studied.Keywords: pressure contour, surface tension, volume of fluid, velocity field
Procedia PDF Downloads 37676 Effects of G-jitter Combined with Heat and Mass Transfer by Mixed Convection MHD Flow of Maxwell Fluid in a Porous Space
Authors: Faisal Salah, Z. A. Aziz, K. K. Viswanathan
Abstract:
In this article, the effects of g-jitter induced and combined with heat and mass transfer by mixed convection of MHD Maxwell fluid in microgravity situation is investigated for a simple system. This system consists of two heated vertical parallel infinite flat plates held at constant but different temperatures and concentrations. By using modified Darcy’s law, the equations governing the flow are modelled. These equations are solved analytically for the induced velocity, temperature and concentration distributions. Many interesting available results in the relevant literature (i.e. Newtonian fluid) is obtained as the special case of the present general analysis. Finally, the graphical results for the velocity profile of the oscillating flow in the channel are presented and discussed for different values of the material constants.Keywords: g-jitter, heat and mass transfer, mixed convection, Maxwell fluid, porous medium
Procedia PDF Downloads 45575 Modeling Study of Short Fiber Orientation in Simple Injection Molding Processes
Authors: Ihsane Modhaffar, Kamal Gueraoui, Abouelkacem Qais, Abderrahmane Maaouni, Samir Men-La-Yakhaf, Hamid Eltourroug
Abstract:
The main objective of this paper is to develop a Computational Fluid Dynamics (CFD) model to simulate and characterize the fiber suspension in flow in rectangular cavities. The model is intended to describe the velocity profile and to predict the fiber orientation. The flow was considered to be incompressible, and behave as Newtonian fluid containing suspensions of short-fibers. The numerical model for determination of velocity profile and fiber orientation during mold-filling stage of injection molding process was solved using finite volume method. The governing equations of this problem are: the continuity, the momentum and the energy. The obtained results were compared to available experimental findings. A good agreement between the numerical results and the experimental data was achieved.Keywords: injection, composites, short-fiber reinforced thermoplastics, fiber orientation, incompressible fluid, numerical simulation
Procedia PDF Downloads 44274 Implicit Eulerian Fluid-Structure Interaction Method for the Modeling of Highly Deformable Elastic Membranes
Authors: Aymen Laadhari, Gábor Székely
Abstract:
This paper is concerned with the development of a fully implicit and purely Eulerian fluid-structure interaction method tailored for the modeling of the large deformations of elastic membranes in a surrounding Newtonian fluid. We consider a simplified model for the mechanical properties of the membrane, in which the surface strain energy depends on the membrane stretching. The fully Eulerian description is based on the advection of a modified surface tension tensor, and the deformations of the membrane are tracked using a level set strategy. The resulting nonlinear problem is solved by a Newton-Raphson method, featuring a quadratic convergence behavior. A monolithic solver is implemented, and we report several numerical experiments aimed at model validation and illustrating the accuracy of the presented method. We show that stability is maintained for significantly larger time steps.Keywords: finite element method, implicit, level set, membrane, Newton method
Procedia PDF Downloads 28073 Unsteady Forced Convection Flow and Heat Transfer Past a Blunt Headed Semi-Circular Cylinder at Low Reynolds Numbers
Authors: Y. El Khchine, M. Sriti
Abstract:
In the present work, the forced convection heat transfer and fluid flow past an unconfined semi-circular cylinder is investigated. The two-dimensional simulation is employed for Reynolds numbers ranging from 10 ≤ Re ≤ 200, employing air (Pr = 0.71) as an operating fluid with Newtonian constant physics property. Continuity, momentum, and energy equations with appropriate boundary conditions are solved using the Computational Fluid Dynamics (CFD) solver Ansys Fluent. Various parameters flow such as lift, drag, pressure, skin friction coefficients, Nusselt number, Strouhal number, and vortex strength are calculated. The transition from steady to time-periodic flow occurs between Re=60 and 80. The effect of the Reynolds number on heat transfer is discussed. Finally, a developed correlation of Nusselt and Strouhal numbers is presented.Keywords: forced convection, semi-circular cylinder, Nusselt number, Prandtl number
Procedia PDF Downloads 8772 The Effect of Incorporation of Inulin as a Fat Replacer on the Quality of Milk Products Vis-À-Vis Ice Cream
Authors: Harish Kumar Sharma
Abstract:
The influence of different levels of inulin as a fat replacer on the quality of ice cream was investigated. The physicochemical, rheological and textural properties of control ice cream and ice cream prepared with inulin in different proportions were determined and correlated to the different parameters using Pearson correlation and Principle Component Analysis (PCA). Based on the overall acepectability, ice cream with 4% inulin was found best and was selected for preparation of ice cream with inulin:SPI in different proportions. Compared with control ice cream, Inulin:SPI showed different rheological properties, resulting in significantly higher apparent viscosities, consistency coefficient and greater deviations from Newtonian flow. In addition, both hardness and melting resistance significantly increased with increase in the SPI content in ice cream prepared with inulin: SPI. Also hardness value increased for inulin based ice cream compared to control ice cream but it melted significantly faster than the latter. Colour value significantly decreased in both the cases compared to the control sample. The deliberation shall focus on the effect of incorporation of inulin on the quality of ice-cream.Keywords: fat replacer, inulin, ice cream, viscosity, principal component analysis
Procedia PDF Downloads 35571 Numerical Simulation of Magnetohydrodynamic (MHD) Blood Flow in a Stenosed Artery
Authors: Sreeparna Majee, G. C. Shit
Abstract:
Unsteady blood flow has been numerically investigated through stenosed arteries to achieve an idea about the physiological blood flow pattern in diseased arteries. The blood is treated as Newtonian fluid and the arterial wall is considered to be rigid having deposition of plaque in its lumen. For direct numerical simulation, vorticity-stream function formulation has been adopted to solve the problem using implicit finite difference method by developing well known Peaceman-Rachford Alternating Direction Implicit (ADI) scheme. The effects of magnetic parameter and Reynolds number on velocity and wall shear stress are being studied and presented quantitatively over the entire arterial segment. The streamlines have been plotted to understand the flow pattern in the stenosed artery, which has significant alterations in the downstream of the stenosis in the presence of magnetic field. The results show that there are nominal changes in the flow pattern when magnetic field strength is enhanced upto 8T which can have remarkable usage to MRI machines.Keywords: magnetohydrodynamics, blood flow, stenosis, energy dissipation
Procedia PDF Downloads 25570 The Kinks, the Solitons, and the Shocks in Series Connected Discrete Josephson Transmission Lines
Authors: Eugene Kogan
Abstract:
We analytically study the localized running waves in the discrete Josephson transmission lines (JTL), constructed from Josephson junctions (JJ) and capacitors. The quasi-continuum approximation reduces the calculation of the running wave properties to the problem of equilibrium of an elastic rod in the potential field. Making additional approximations, we reduce the problem to the motion of the fictitious Newtonian particle in the potential well. We show that there exist running waves in the form of supersonic kinks and solitons and calculate their velocities and profiles. We show that the nonstationary smooth waves, which are small perturbations on the homogeneous non-zero background, are described by Korteweg-de Vries equation, and those on zero background -by the modified Korteweg-de Vries equation. We also study the effect of dissipation on the running waves in JTL and find that in the presence of the resistors, shunting the JJ and/or in series with the ground capacitors, the only possible stationary running waves are the shock waves, whose profiles are also found.Keywords: Josephson transmission line, shocks, solitary waves, nonlinear waves
Procedia PDF Downloads 9069 Flow of a Second Order Fluid through Constricted Tube with Slip Velocity at Wall Using Integral Method
Authors: Nosheen Zareen Khan, Abdul Majeed Siddiqui, Muhammad Afzal Rana
Abstract:
The steady flow of a second order fluid through constricted tube with slip velocity at wall is modeled and analyzed theoretically. The governing equations are simplified by implying no slip in radial direction. Based on Karman Pohlhausen procedure polynomial solution for axial velocity profile is presented. An expressions for pressure gradient, shear stress, separation and reattachment points and radial velocity are also calculated. The effect of slip and no slip velocity on velocity, shear stress, pressure gradient are discussed and depicted graphically. It is noted that when Reynolds number increases velocity of the fluid decreases in both slip and no slip conditions. It is also found that the wall shear stress, separation and reattachment points are strongly effected by Reynolds number.Keywords: approximate solution, constricted tube, non-Newtonian fluids, Reynolds number
Procedia PDF Downloads 37168 On Energy Condition Violation for Shifting Negative Mass Black Holes
Authors: Manuel Urueña Palomo
Abstract:
In this paper, we introduce the study of a new solution to gravitational singularities by violating the energy conditions of the Penrose Hawking singularity theorems. We consider that a shift to negative energies, and thus, to negative masses, takes place at the event horizon of a black hole, justified by the original, singular and exact Schwarzschild solution. These negative energies are supported by relativistic particle physics considering the negative energy solutions of the Dirac equation, which states that a time transformation shifts to a negative energy particle. In either general relativity or full Newtonian mechanics, these negative masses are predicted to be repulsive. It is demonstrated that the model fits actual observations, and could possibly clarify the size of observed and unexplained supermassive black holes, when considering the inflation that would take place inside the event horizon where massive particles interact antigravitationally. An approximated solution of the model proposed could be simulated in order to compare it with these observations.Keywords: black holes, CPT symmetry, negative mass, time transformation
Procedia PDF Downloads 12167 Electrokinetic Transport of Power Law Fluid through Hydrophobic Micro-Slits
Authors: Ainul Haque, Ameeye Kumar Nayak
Abstract:
Flow enhancement and species transport in a slit hydrophobic microchannel is studied for non-Newtonian fluids with the externally imposed electric field and pressure gradient. The incompressible Poisson-Nernst-Plank equations and the Navier-Stokes equations are approximated by lubrication theory to quantify the flow structure due to hydrophobic and hydrophilic surfaces. The analytical quantification of velocity and pressure of electroosmotic flow (EOF) is made with the numerical results due to the staggered grid based finite volume method for flow governing equations. The resistance force due to fluid friction and shear force along the surface are decreased by the hydrophobicity, enables the faster movement of fluid particles. The resulting flow enhancement factor Ef is increased with the low viscous fluid and provides maximum species transport. Also, the analytical comparison of EOF with pressure driven EOF justifies the flow enhancement due to hydrophobicity and shear impact on flow variation.Keywords: electroosmotic flow, hydrophobic surface, power-law fluid, shear effect
Procedia PDF Downloads 34166 Effect of Slip Condition and Magnetic Field on Unsteady MHD Thin Film Flow of a Third Grade Fluid with Heat Transfer down an Inclined Plane
Authors: Y. M. Aiyesimi, G. T. Okedayo, O. W. Lawal
Abstract:
The analysis has been carried out to study unsteady MHD thin film flow of a third grade fluid down an inclined plane with heat transfer when the slippage between the surface of plane and the lower surface of the fluid is valid. The governing nonlinear partial differential equations involved are reduced to linear partial differential equations using regular perturbation method. The resulting equations were solved analytically using method of separation of variable and eigenfunctions expansion. The solutions obtained were examined and discussed graphically. It is interesting to find that the variation of the velocity and temperature profile with the slip and magnetic field parameter depends on time.Keywords: non-Newtonian fluid, MHD flow, thin film flow, third grade fluid, slip boundary condition, heat transfer, separation of variable, eigenfunction expansion
Procedia PDF Downloads 35465 Entropy Generation Analyze Due to the Steady Natural Convection of Newtonian Fluid in a Square Enclosure
Authors: T. T. Naas, Y. Lasbet, C. Kezrane
Abstract:
The thermal control in many systems is widely accomplished applying mixed convection process due to its low cost, reliability and easy maintenance. Typical applications include the aircraft electronic equipment, rotating-disc heat exchangers, turbo machinery, and nuclear reactors, etc. Natural convection in an inclined square enclosure heated via wall heater has been studied numerically. Finite volume method is used for solving momentum and energy equations in the form of stream function–vorticity. The right and left walls are kept at a constant temperature, while the other parts are adiabatic. The range of the inclination angle covers a whole revolution. The method is validated for a vertical cavity. A general power law dependence of the Nusselt number with respect to the Rayleigh number with the coefficient and exponent as functions of the inclination angle is presented. For a fixed Rayleigh number, the inclination angle increases or decreases is found.Keywords: natural convection in enclosure, inclined enclosure, Nusselt number, entropy generation analyze
Procedia PDF Downloads 22764 Numerical Analysis of the Computational Fluid Dynamics of Co-Digestion in a Large-Scale Continuous Stirred Tank Reactor
Authors: Sylvana A. Vega, Cesar E. Huilinir, Carlos J. Gonzalez
Abstract:
Co-digestion in anaerobic biodigesters is a technology improving hydrolysis by increasing methane generation. In the present study, the dimensional computational fluid dynamics (CFD) is numerically analyzed using Ansys Fluent software for agitation in a full-scale Continuous Stirred Tank Reactor (CSTR) biodigester during the co-digestion process. For this, a rheological study of the substrate is carried out, establishing rotation speeds of the stirrers depending on the microbial activity and energy ranges. The substrate is organic waste from industrial sources of sanitary water, butcher, fishmonger, and dairy. Once the rheological behavior curves have been obtained, it is obtained that it is a non-Newtonian fluid of the pseudoplastic type, with a solids rate of 12%. In the simulation, the rheological results of the fluid are considered, and the full-scale CSTR biodigester is modeled. It was coupling the second-order continuity differential equations, the three-dimensional Navier Stokes, the power-law model for non-Newtonian fluids, and three turbulence models: k-ε RNG, k-ε Realizable, and RMS (Reynolds Stress Model), for a 45° tilt vane impeller. It is simulated for three minutes since it is desired to study an intermittent mixture with a saving benefit of energy consumed. The results show that the absolute errors of the power number associated with the k-ε RNG, k-ε Realizable, and RMS models were 7.62%, 1.85%, and 5.05%, respectively, the numbers of power obtained from the analytical-experimental equation of Nagata. The results of the generalized Reynolds number show that the fluid dynamics have a transition-turbulent flow regime. Concerning the Froude number, the result indicates there is no need to implement baffles in the biodigester design, and the power number provides a steady trend close to 1.5. It is observed that the levels of design speeds within the biodigester are approximately 0.1 m/s, which are speeds suitable for the microbial community, where they can coexist and feed on the substrate in co-digestion. It is concluded that the model that more accurately predicts the behavior of fluid dynamics within the reactor is the k-ε Realizable model. The flow paths obtained are consistent with what is stated in the referenced literature, where the 45° inclination PBT impeller is the right type of agitator to keep particles in suspension and, in turn, increase the dispersion of gas in the liquid phase. If a 24/7 complete mix is considered under stirred agitation, with a plant factor of 80%, 51,840 kWh/year are estimated. On the contrary, if intermittent agitations of 3 min every 15 min are used under the same design conditions, reduce almost 80% of energy costs. It is a feasible solution to predict the energy expenditure of an anaerobic biodigester CSTR. It is recommended to use high mixing intensities, at the beginning and end of the joint phase acetogenesis/methanogenesis. This high intensity of mixing, in the beginning, produces the activation of the bacteria, and once reaching the end of the Hydraulic Retention Time period, it produces another increase in the mixing agitations, favoring the final dispersion of the biogas that may be trapped in the biodigester bottom.Keywords: anaerobic co-digestion, computational fluid dynamics, CFD, net power, organic waste
Procedia PDF Downloads 8463 Numerical Solution of Momentum Equations Using Finite Difference Method for Newtonian Flows in Two-Dimensional Cartesian Coordinate System
Authors: Ali Ateş, Ansar B. Mwimbo, Ali H. Abdulkarim
Abstract:
General transport equation has a wide range of application in Fluid Mechanics and Heat Transfer problems. In this equation, generally when φ variable which represents a flow property is used to represent fluid velocity component, general transport equation turns into momentum equations or with its well known name Navier-Stokes equations. In these non-linear differential equations instead of seeking for analytic solutions, preferring numerical solutions is a more frequently used procedure. Finite difference method is a commonly used numerical solution method. In these equations using velocity and pressure gradients instead of stress tensors decreases the number of unknowns. Also, continuity equation, by integrating the system, number of equations is obtained as number of unknowns. In this situation, velocity and pressure components emerge as two important parameters. In the solution of differential equation system, velocities and pressures must be solved together. However, in the considered grid system, when pressure and velocity values are jointly solved for the same nodal points some problems confront us. To overcome this problem, using staggered grid system is a referred solution method. For the computerized solutions of the staggered grid system various algorithms were developed. From these, two most commonly used are SIMPLE and SIMPLER algorithms. In this study Navier-Stokes equations were numerically solved for Newtonian flow, whose mass or gravitational forces were neglected, for incompressible and laminar fluid, as a hydro dynamically fully developed region and in two dimensional cartesian coordinate system. Finite difference method was chosen as the solution method. This is a parametric study in which varying values of velocity components, pressure and Reynolds numbers were used. Differential equations were discritized using central difference and hybrid scheme. The discritized equation system was solved by Gauss-Siedel iteration method. SIMPLE and SIMPLER were used as solution algorithms. The obtained results, were compared for central difference and hybrid as discritization methods. Also, as solution algorithm, SIMPLE algorithm and SIMPLER algorithm were compared to each other. As a result, it was observed that hybrid discritization method gave better results over a larger area. Furthermore, as computer solution algorithm, besides some disadvantages, it can be said that SIMPLER algorithm is more practical and gave result in short time. For this study, a code was developed in DELPHI programming language. The values obtained in a computer program were converted into graphs and discussed. During sketching, the quality of the graph was increased by adding intermediate values to the obtained result values using Lagrange interpolation formula. For the solution of the system, number of grid and node was found as an estimated. At the same time, to indicate that the obtained results are satisfactory enough, by doing independent analysis from the grid (GCI analysis) for coarse, medium and fine grid system solution domain was obtained. It was observed that when graphs and program outputs were compared with similar studies highly satisfactory results were achieved.Keywords: finite difference method, GCI analysis, numerical solution of the Navier-Stokes equations, SIMPLE and SIMPLER algoritms
Procedia PDF Downloads 36462 Estimation of Damping Force of Double Ended Shear Mode Magnetorheological Damper Using Computational Analysis
Authors: Gurubasavaraju T. M.
Abstract:
The magnetorheological (MR) damper could provide variable damping force with respect to the different input magnetic field. The damping force could be estimated through computational analysis using finite element and computational fluid dynamics analysis. The double-ended damper operates without changing the total volume of fluid. In this paper, damping force of double ended damper under different magnetic field is computed. Initially, the magneto-statics analysis carried out to evaluate the magnetic flux density across the fluid flow gap. The respective change in the rheology of the MR fluid is computed by using the experimentally fitted polynomial equation of shear stress versus magnetic field plot of MR fluid. The obtained values are substituted in the Herschel Buckley model to express the non-Newtonian behavior of MR fluid. Later, using computational fluid dynamic (CFD) analysis damping characteristics in terms of force versus velocity and force versus displacement for the respective magnetic field is estimated. The purpose of the present approach is to characterize the preliminary designed MR damper before fabricating.Keywords: MR fluid, double ended MR damper, CFD, FEA
Procedia PDF Downloads 155