Source author record

Taha Sochi

Taha Sochi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

52works
18topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

52 published item(s)

preprint2016arXiv

Collision strengths and transition probabilities for Co III forbidden lines

In this paper we compute the collision strengths and their thermally-averaged Maxwellian values for electron transitions between the fifteen lowest levels of doubly-ionised cobalt, Co^{2+}, which give rise to forbidden emission lines in the visible and infrared region of spectrum. The calculations also include transition probabilities and predicted relative line emissivities. The data are particularly useful for analysing the thermodynamic conditions of supernova ejecta.

preprint2016arXiv

Principles of Differential Geometry

The present text is a collection of notes about differential geometry prepared to some extent as part of tutorials about topics and applications related to tensor calculus. They can be regarded as continuation to the previous notes on tensor calculus as they are based on the materials and conventions given in those documents. They can be used as a reference for a first course on the subject or as part of a course on tensor calculus.

preprint2016arXiv

Special Relativity: Scientific or Philosophical Theory?

In this article, we argue that the theory of special relativity, as formulated by Einstein, is a philosophical rather than a scientific theory. What is scientific and experimentally supported is the formalism of the relativistic mechanics embedded in the Lorentz transformations and their direct mathematical, experimental and observational consequences. This is in parallel with the quantum mechanics where the scientific content and experimental support of this branch of physics is embedded in the formalism of quantum mechanics and not in its philosophical interpretations such as the Copenhagen school or the parallel worlds explanations. Einstein theory of special relativity gets unduly credit from the success of the relativistic mechanics of Lorentz transformations. Hence, all the postulates and consequences of Einstein interpretation which have no direct experimental or observational support should be reexamined and the relativistic mechanics of Lorentz transformations should be treated in education, academia and research in a similar fashion to that of quantum mechanics.

preprint2016arXiv

Tensor Calculus

These notes are the second part of the tensor calculus documents which started with the previous set of introductory. In the present text, we continue the discussion of selected topics of the subject at a higher level expanding, when necessary, some topics and developing further concepts and techniques. Unlike the previous notes which are largely based on a Cartesian approach, the present notes are essentially based on assuming an underlying general curvilinear coordinate system.

preprint2015arXiv

Analytical solutions for the flow of Carreau and Cross fluids in circular pipes and thin slits

In this paper, analytical expressions correlating the volumetric flow rate to the pressure drop are derived for the flow of Carreau and Cross fluids through straight rigid circular uniform pipes and long thin slits. The derivation is based on the application of Weissenberg-Rabinowitsch-Mooney-Schofield method to obtain flow solutions for generalized Newtonian fluids through pipes and our adaptation of this method to the flow through slits. The derived expressions are validated by comparing their solutions to the solutions obtained from direct numerical integration. They are also validated by comparison to the solutions obtained from the variational method which we proposed previously. In all the investigated cases, the three methods agree very well. The agreement with the variational method also lends more support to this method and to the variational principle which the method is based upon.

preprint2015arXiv

Collision strengths and transition probabilities for Co II infrared forbidden lines

We calculate collision strengths and their thermally-averaged Maxwellian values for electron excitation and de-excitation between the fifteen lowest levels of singly-ionised cobalt, Co+, which give rise to emission lines in the near- and mid-infrared. Transition probabilities are also calculated and relative line intensities predicted for conditions typical of supernova ejecta. The diagnostic potential of the 10.52, 15.46 and 14.74 micro-metre transition lines is briefly discussed.

preprint2015arXiv

Effective collision strengths for excitation and de-excitation of nebular [O III] optical and infrared lines with kappa distributed electron energies

We present effective collision strengths for electron excitation and de-excitation of the ten forbidden transitions between the five lowest energy levels of the astronomically abundant doubly-ionised oxygen ion, O^{2+}. The raw collision strength data were obtained from an R-matrix intermediate coupling calculation using the Breit-Pauli relativistic approximation published previously by the authors. The effective collision strengths were calculated with kappa-distributed electron energies and are tabulated as a function of the electron temperature and kappa.

preprint2015arXiv

Further validation to the variational method to obtain flow relations for generalized Newtonian fluids

We continue our investigation to the use of the variational method to derive flow relations for generalized Newtonian fluids in confined geometries. While in the previous investigations we used the straight circular tube geometry with eight fluid rheological models to demonstrate and establish the variational method, the focus here is on the plane long thin slit geometry using those eight rheological models, namely: Newtonian, power law, Ree-Eyring, Carreau, Cross, Casson, Bingham and Herschel-Bulkley. We demonstrate how the variational principle based on minimizing the total stress in the flow conduit can be used to derive analytical expressions, which are previously derived by other methods, or used in conjunction with numerical procedures to obtain numerical solutions which are virtually identical to the solutions obtained previously from well established methods of fluid dynamics. In this regard, we use the method of Weissenberg-Rabinowitsch-Mooney-Schofield (WRMS), with our adaptation from the circular pipe geometry to the long thin slit geometry, to derive analytical formulae for the eight types of fluid where these derived formulae are used for comparison and validation of the variational formulae and numerical solutions. Although some examples may be of little value, the optimization principle which the variational method is based upon has a significant theoretical value as it reveals the tendency of the flow system to assume a configuration that minimizes the total stress. Our proposal also offers a new methodology to tackle common problems in fluid dynamics and rheology.

preprint2015arXiv

Reply to "Comment on Sochi's variational method for generalised Newtonian flow" by Pritchard and Corson

In this article we challenge the claim that the previously proposed variational method to obtain flow solutions for generalized Newtonian fluids in circular tubes and plane slits is exact only for power law fluids. We also defend the theoretical foundation and formalism of the method which is based on minimizing the total stress through the application of the Euler-Lagrange principle.

preprint2015arXiv

Using the stress function in the flow of generalized Newtonian fluids through conduits with non-circular or multiply connected cross sections

We investigate the possibility that the spatial dependency of stress in generalized Newtonian flow systems is a function of the applied pressure field and the conduit geometry but not of the fluid rheology. This possibility is well established for the case of a one-dimensional flow through simply connected regions, specifically tubes of circular uniform cross sections and plane thin slits. If it can also be established for the more general case of generalized Newtonian flow through non-circular or multiply connected geometries, such as the two-dimensional flow through conduits of rectangular or elliptical cross sections or the flow through annular circular pipes, then analytical or semi-analytical or highly accurate numerical solutions; regarding stress, rate of strain, velocity profile and volumetric flow rate; for these geometries can be obtained from the stress function, which can be easily obtained from the Newtonian case, in combination with the constitutive rheological relation for the particular non-Newtonian fluid, as done previously for the case of the one-dimensional flow through simply connected regions.

preprint2015arXiv

Using the stress function in the flow of generalized Newtonian fluids through pipes and slits

We use a generic and general numerical method to obtain solutions for the flow of generalized Newtonian fluids through circular pipes and plane slits. The method, which is simple and robust can produce highly accurate solutions which virtually match any analytical solutions. The method is based on employing the stress, as a function of the pipe radius or slit thickness dimension, combined with the rate of strain function as represented by the fluid rheological constitutive relation that correlates the rate of strain to stress. Nine types of generalized Newtonian fluids are tested in this investigation and the solutions obtained from the generic method are compared to the analytical solutions which are obtained from the Weissenberg-Rabinowitsch-Mooney-Schofield method. Very good agreement was obtained in all the investigated cases. All the required quantities of the flow which include local viscosity, rate of strain, flow velocity profile and volumetric flow rate, as well as shear stress, can be obtained from the generic method. This is an advantage as compared to some traditional methods which only produce some of these quantities. The method is also superior to the numerical meshing techniques which may be used for resolving the flow in these systems. The method is particularly useful when analytical solutions are not available or when the available analytical solutions do not yield all the flow parameters.

preprint2015arXiv

Variational approach for resolving the flow of generalized Newtonian fluids in circular pipes and plane slits

In this paper, we use a generic and general variational method to obtain solutions to the flow of generalized Newtonian fluids through circular pipes and plane slits. The new method is not based on the use of the Euler-Lagrange variational principle and hence it is totally independent of our previous approach which is based on this principle. Instead, the method applies a very generic and general optimization approach which can be justified by the Dirichlet principle although this is not the only possible theoretical justification. The results that were obtained from the new method using nine types of fluid are in total agreement, within certain restrictions, with the results obtained from the traditional methods of fluid mechanics as well as the results obtained from the previous variational approach. In addition to being a useful method in its own for resolving the flow field in circular pipes and plane slits, the new variational method lends more support to the old variational method as well as for the use of variational principles in general to resolve the flow of generalized Newtonian fluids and obtain all the quantities of the flow field which include shear stress, local viscosity, rate of strain, speed profile and volumetric flow rate. The theoretical basis of the new variational method, which rests on the use of the Dirichlet principle, also provides theoretical support to the former variational method.

preprint2015arXiv

Variational approach for the flow of Ree-Eyring and Casson fluids in pipes

The flow of Ree-Eyring and Casson non-Newtonian fluids is investigated using a variational principle to optimize the total stress. The variationally-obtained solutions are compared to the analytical solutions derived from the Weissenberg-Rabinowitsch-Mooney equation and the results are found to be identical within acceptable numerical errors and modeling approximations.

preprint2014arXiv

Deterministic and stochastic algorithms for resolving the flow fields in ducts and networks using energy minimization

Several deterministic and stochastic multi-variable global optimization algorithms (Conjugate Gradient, Nelder-Mead, Quasi-Newton, and Global) are investigated in conjunction with energy minimization principle to resolve the pressure and volumetric flow rate fields in single ducts and networks of interconnected ducts. The algorithms are tested with seven types of fluid: Newtonian, power law, Bingham, Herschel-Bulkley, Ellis, Ree-Eyring and Casson. The results obtained from all those algorithms for all these types of fluid agree very well with the analytically derived solutions as obtained from the traditional methods which are based on the conservation principles and fluid constitutive relations. The results confirm and generalize the findings of our previous investigations that the energy minimization principle is at the heart of the flow dynamics systems. The investigation also enriches the methods of Computational Fluid Dynamics for solving the flow fields in tubes and networks for various types of Newtonian and non-Newtonian fluids.

preprint2014arXiv

Emission and recombination coefficients for hydrogen with kappa-distributed electron energies

We provide a data set of emission and recombination coefficients of hydrogen using a kappa-distribution of electron energies rather than the more traditional Maxwell-Boltzmann (MB) distribution. The data are mainly relevant to thin and relatively cold plasma found in planetary nebulae and H II regions. The data set extends the previous data sets provided by Storey and Hummer which were computed using a MB distribution. The data set, which is placed in the public domain, is structured as a function of electron number density, temperature and kappa. Interactive fortran 77 and C++ data servers are also provided as an accessory to probe the data and obtain Lagrange-interpolated values for any choice of all three variables between the explicitly computed values.

preprint2014arXiv

Energy minimization for the flow in ducts and networks

The present paper is an attempt to demonstrate how the energy minimization principle may be considered as a governing rule for the physical equilibrium that determines the flow fields in tubes and networks. We previously investigated this issue using a numerical stochastic method, specifically simulated annealing, where we demonstrated the problem by some illuminating examples and concluded that energy minimization principle can be a valid hypothesis. The investigation in this paper is more general as it is based to a certain extent on an analytical approach.

preprint2014arXiv

Solving the flow fields in conduits and networks using energy minimization principle with simulated annealing

In this paper, we propose and test an intuitive assumption that the pressure field in single conduits and networks of interconnected conduits adjusts itself to minimize the total energy consumption required for transporting a specific quantity of fluid. We test this assumption by using linear flow models of Newtonian fluids transported through rigid tubes and networks in conjunction with a simulated annealing (SA) protocol to minimize the total energy cost. All the results confirm our hypothesis as the SA algorithm produces very close results to those obtained from the traditional deterministic methods of identifying the flow fields by solving a set of simultaneous equations based on the conservation principles. The same results apply to electric ohmic conductors and networks of interconnected ohmic conductors. Computational experiments conducted in this regard confirm this extension. Further studies are required to test the energy minimization hypothesis for the non-linear flow systems.

preprint2014arXiv

The Continuum Emission Spectrum of Hf 2-2 near the Balmer Limit and the ORL versus CEL abundance and temperature Discrepancy

The continuum spectrum of the planetary nebula Hf 2-2 close to the Balmer discontinuity is modeled in the context of the long standing problem of the abundance and temperature discrepancy found when analyzing optical recombination lines and collisionally excited forbidden lines in nebulae. Models are constructed using single and double Maxwell-Boltzmann distributions as well as kappa-distributions for the energies of the free electrons. New results for the necessary continuum and line emission coefficients are presented calculated with kappa-distributed energies. The best fit to the observed continuum spectrum is found to be a model comprising two components with dramatically different temperatures and with a Maxwell-Boltzmann distribution of electron energies. On the basis of a chi-squared analysis, this model is strongly favored over a model with kappa-distributed electron energies.

preprint2014arXiv

The Flow of Newtonian and power law fluids in elastic tubes

We derive analytical expressions for the flow of Newtonian and power law fluids in elastic circularly-symmetric tubes based on a lubrication approximation where the flow velocity profile at each cross section is assumed to have its axially-dependent characteristic shape for the given rheology and cross sectional size. Two pressure-area constitutive elastic relations for the tube elastic response are used in these derivations. We demonstrate the validity of the derived equations by observing qualitatively correct trends in general and quantitatively valid asymptotic convergence to limiting cases. The Newtonian formulae are compared to similar formulae derived previously from a one-dimensional version of the Navier-Stokes equations.

preprint2014arXiv

The flow of power law fluids in elastic networks and porous media

The flow of power law fluids, which include shear thinning and shear thickening as well as Newtonian as a special case, in networks of interconnected elastic tubes is investigated using a residual based pore scale network modeling method with the employment of newly derived formulae. Two relations describing the mechanical interaction between the local pressure and local cross sectional area in distensible tubes of elastic nature are considered in the derivation of these formulae. The model can be used to describe shear dependent flows of mainly viscous nature. The behavior of the proposed model is vindicated by several tests in a number of special and limiting cases where the results can be verified quantitatively or qualitatively. The model, which is the first of its kind, incorporates more than one major non-linearity corresponding to the fluid rheology and conduit mechanical properties, that is non-Newtonian effects and tube distensibility. The formulation, implementation and performance indicate that the model enjoys certain advantages over the existing models such as being exact within the restricting assumptions on which the model is based, easy implementation, low computational costs, reliability and smooth convergence. The proposed model can therefore be used as an alternative to the existing Newtonian distensible models; moreover it stretches the capabilities of the existing modeling approaches to reach non-Newtonian rheologies.

preprint2013arXiv

Collision Strengths for [O III] Optical and Infrared Lines

We present electron collision strengths and their thermally averaged values for the forbidden lines of the astronomically abundant doubly-ionized oxygen ion, O^{2+}, in an intermediate coupling scheme using the Breit-Pauli relativistic terms as implemented in an R-matrix atomic scattering code. We use several atomic targets for the R-matrix scattering calculations including one with 72 atomic terms. We also compare with new results obtained using the intermediate coupling frame transformation method. We find spectroscopically significant differences against a recent Breit-Pauli calculation for the excitation of the [O III] lambda 4363 transition but confirm the results of earlier calculations.

preprint2013arXiv

Comparing Poiseuille with 1D Navier-Stokes Flow in Rigid and Distensible Tubes and Networks

A comparison is made between the Hagen-Poiseuille flow in rigid tubes and networks on one side and the time-independent one-dimensional Navier-Stokes flow in elastic tubes and networks on the other. Analytical relations, a Poiseuille network flow model and two finite element Navier-Stokes one-dimensional flow models have been developed and used in this investigation. The comparison highlights the differences between Poiseuille and one-dimensional Navier-Stokes flow models which may have been unjustifiably treated as equivalent in some studies.

preprint2013arXiv

Flow of Navier-Stokes Fluids in Converging-Diverging Distensible Tubes

We use a method based on the lubrication approximation in conjunction with a residual-based mass-continuity iterative solution scheme to compute the flow rate and pressure field in distensible converging-diverging tubes for Navier-Stokes fluids. We employ an analytical formula derived from a one-dimensional version of the Navier-Stokes equations to describe the underlying flow model that provides the residual function. This formula correlates the flow rate to the boundary pressures in straight cylindrical elastic tubes with constant-radius. We validate our findings by the convergence toward a final solution with fine discretization as well as by comparison to the Poiseuille-type flow in its convergence toward analytic solutions found earlier in rigid converging-diverging tubes. We also tested the method on limiting special cases of cylindrical elastic tubes with constant-radius where the numerical solutions converged to the expected analytical solutions. The distensible model has also been endorsed by its convergence toward the rigid Poiseuille-type model with increasing the tube wall stiffness. Lubrication-based one-dimensional finite element method was also used for verification. In this investigation five converging-diverging geometries are used for demonstration, validation and as prototypes for modeling converging-diverging geometries in general.

preprint2013arXiv

Flow of Navier-Stokes Fluids in Cylindrical Elastic Tubes

Analytical expressions correlating the volumetric flow rate to the inlet and outlet pressures are derived for the time-independent flow of Newtonian fluids in cylindrically-shaped elastic tubes using a one-dimensional Navier-Stokes flow model with two pressure-area constitutive relations. These expressions for elastic tubes are the equivalent of Poiseuille and Poiseuille-type expressions for rigid tubes which were previously derived for the flow of Newtonian and non-Newtonian fluids under various flow conditions. Formulae and procedures for identifying the pressure field and tube geometric profile are also presented. The results are validated by a finite element method implementation. Sensible trends in the analytical and numerical results are observed and documented.

preprint2013arXiv

Flow of non-Newtonian Fluids in Converging-Diverging Rigid Tubes

A residual-based lubrication method is used in this paper to find the flow rate and pressure field in converging-diverging rigid tubes for the flow of time-independent category of non-Newtonian fluids. Five converging-diverging prototype geometries were used in this investigation in conjunction with two fluid models: Ellis and Herschel-Bulkley. The method was validated by convergence behavior sensibility tests, convergence to analytical solutions for the straight tubes as special cases for the converging-diverging tubes, convergence to analytical solutions found earlier for the flow in converging-diverging tubes of Newtonian fluids as special cases for non-Newtonian, and convergence to analytical solutions found earlier for the flow of power-law fluids in converging-diverging tubes. A brief investigation was also conducted on a sample of diverging-converging geometries. The method can in principle be extended to the flow of viscoelastic and thixotropic/rheopectic fluid categories. The method can also be extended to geometries varying in size and shape in the flow direction, other than the perfect cylindrically-symmetric converging-diverging ones, as long as characteristic flow relations correlating the flow rate to the pressure drop on the discretized elements of the lubrication approximation can be found. These relations can be analytical, empirical and even numerical and hence the method has a wide applicability range.

preprint2013arXiv

New program with new approach for spectral data analysis

This article presents a high-throughput computer program, called EasyDD, for batch processing, analyzing and visualizing of spectral data; particularly those related to the new generation of synchrotron detectors and X-ray powder diffraction applications. This computing tool is designed for the treatment of large volumes of data in reasonable time with affordable computational resources. A case study in which this program was used to process and analyze powder diffraction data obtained from the ESRF synchrotron on an alumina-based nickel nanoparticle catalysis system is also presented for demonstration. The development of this computing tool, with the associated protocols, is inspired by a novel approach in spectral data analysis.

preprint2013arXiv

Pore-Scale Modeling of Navier-Stokes Flow in Distensible Networks and Porous Media

In this paper, a pore-scale network modeling method, based on the flow continuity residual in conjunction with a Newton-Raphson non-linear iterative solving technique, is proposed and used to obtain the pressure and flow fields in a network of interconnected distensible ducts representing, for instance, blood vasculature or deformable porous media. A previously derived analytical expression correlating boundary pressures to volumetric flow rate in compliant tubes for a pressure-area constitutive elastic relation has been used to represent the underlying flow model. Comparison to a preceding equivalent method, the one-dimensional Navier-Stokes finite element, was made and the results were analyzed. The advantages of the new method have been highlighted and practical computational issues, related mainly to the rate and speed of convergence, have been discussed.

preprint2013arXiv

The Yield Condition in the Mobilization of Yield-Stress Materials in Distensible Tubes

In this paper we investigate the yield condition in the mobilization of yield-stress materials in distensible tubes. We discuss the two possibilities for modeling the yield-stress materials prior to yield: solid-like materials and highly-viscous fluids and identify the logical consequences of these two approaches on the yield condition. As part of this investigation we derive an analytical expression for the pressure field inside a distensible tube with a Newtonian flow using a one-dimensional Navier-Stokes flow model in conjunction with a pressure-area constitutive relation based on elastic tube wall characteristics.

preprint2013arXiv

Using Euler-Lagrange Variational Principle to Obtain Flow Relations for Generalized Newtonian Fluids

Euler-Lagrange variational principle is used to obtain analytical and numerical flow relations in cylindrical tubes. The method is based on minimizing the total stress in the flow duct using the fluid constitutive relation between stress and rate of strain. Newtonian and non-Newtonian fluid models; which include power law, Bingham, Herschel-Bulkley, Carreau and Cross; are used for demonstration.

preprint2013arXiv

Yield and Solidification of Yield-Stress Materials in Rigid Networks and Porous Structures

In this paper, we address the issue of threshold yield pressure of yield-stress materials in rigid networks of interconnected conduits and porous structures subject to a pressure gradient. We compare the results as obtained dynamically from solving the pressure field to those obtained statically from tracing the path of the minimum sum of threshold yield pressures of the individual conduits by using the threshold path algorithms. We refute criticisms directed recently to our previous findings that the pressure field solution generally produces a higher threshold yield pressure than the one obtained by the threshold path algorithms. Issues related to the solidification of yield stress materials in their transition from fluid phase to solid state have also been investigated and assessed as part of the investigation of the yield point.

preprint2012arXiv

Atomic and Molecular Aspects of Astronomical Spectra

In the first section we present the atomic part where a C2+ atomic target was prepared and used to generate theoretical data to investigate recombination lines arising from electron-ion collisions in thin plasma. R-matrix method was used to describe the C2+ plus electron system. Theoretical data concerning bound and autoionizing states were generated in the intermediate-coupling approximation. The data were used to generate dielectronic recombination data for C+ which include transition lines, oscillator strengths, radiative transition probabilities, emissivities and dielectronic recombination coefficients. The data were cast in a line list containing 6187 optically-allowed transitions which include many C II lines observed in astronomical spectra. This line list was used to analyze the spectra from a number of astronomical objects, mainly planetary nebulae, and identify their electron temperature. The electron temperature investigation was also extended to include free electron energy analysis to investigate the long-standing problem of discrepancy between the results of recombination and forbidden lines analysis and its possible connection to the electron distribution. In the second section we present the results of our molecular investigation; the generation of a comprehensive, calculated line list of frequencies and transition probabilities for H2D+. The line list contains over 22 million rotational-vibrational transitions occurring between more than 33 thousand energy levels and covers frequencies up to 18500 cm-1. About 15% of these levels are fully assigned with approximate rotational and vibrational quantum numbers. A temperature-dependent partition function and cooling function are presented. Temperature-dependent synthetic spectra for the temperatures T=100, 500, 1000 and 2000 K in the frequency range 0-10000 cm-1 were also generated and presented graphically.

preprint2012arXiv

Dielectronic Recombination Lines of C+

The current paper presents atomic data generated to investigate the recombination lines of C II in the spectra of planetary nebulae. These data include energies of bound and autoionizing states, oscillator strengths and radiative transition probabilities, autoionization probabilities, and recombination coefficients. The R-matrix method of electron scattering theory was used to describe the C2+ plus electron system.

preprint2012arXiv

Electron Temperatures and Free-Electron Energy Distributions of Nebulae from C II Dielectronic Recombination Lines

A recently generated theoretical line list of C II dielectronic recombination lines together with observational data gathered from the literature is used to investigate the electron temperature in a range of astronomical objects, mainly planetary nebulae. The electron temperature is obtained by a least-squares optimisation using all the reliable observed lines in each object. In addition, the subset of lines arising directly from autoionising states is used to directly determine the free-electron energy distribution which is then compared with various theoretical possibilities. The method described here can potentially determine whether there are departures from Maxwell-Boltzmann distributions in some nebulae, as has been recently proposed. Using published observations of the three planetary nebulae where the relevant lines are recorded, we find that the data are best matched by Maxwell-Boltzmann distributions but that the uncertainties are sufficiently large at present that kappa-distributions or two-component nebular models are not excluded.

preprint2012arXiv

Newtonian Flow in Converging-Diverging Capillaries

The one-dimensional Navier-Stokes equations are used to derive analytical expressions for the relation between pressure and volumetric flow rate in capillaries of five different converging-diverging axisymmetric geometries for Newtonian fluids. The results are compared to previously-derived expressions for the same geometries using the lubrication approximation. The results of the one-dimensional Navier-Stokes are identical to those obtained from the lubrication approximation within a non-dimensional numerical factor. The derived flow expressions have also been validated by comparison to numerical solutions obtained from discretization with numerical integration. Moreover, they have been certified by testing the convergence of solutions as the converging-diverging geometries approach the limiting straight geometry.

preprint2012arXiv

Recombination Lines of CII in the Spectra of Planetary Nebulae

The current report presents the work to investigate the recombination lines of CII in the spectra of planetary nebulae. Two CIII targets were prepared and used to generate theoretical data required in the investigation of recombination lines that arise from collisions between electrons and ions in thin plasma found in planetary nebulae and other astrophysical objects. One of these targets contains 9 atomic terms while the other contains 26 terms. For each one of these targets, theoretical data concerning bound and autoionizing states were generated in the intermediate coupling approximation by R-matrix and Autostructure codes and compared to experimental data. The comparison revealed very good agreement. These theoretical data were then used to generate emissivity data and compare it to the carbon recombination lines found in the observational line list of Zhang et al [2005] on the planetary nebula NGC 7027. The main tool used in this analysis is the `Emissivity' code which is a program developed by the author [2010] in the course of this investigation to calculate the emissivity of transition lines from resonances and subsequent cascade decay. The results of the analysis from both targets using a least-squares optimization technique indicate a temperature of 10200-10500K of the line-emitting region.

preprint2011arXiv

Slip at Fluid-Solid Interface

The `no-slip' is a fundamental assumption and generally-accepted boundary condition in rheology, tribology and fluid mechanics with strong experimental support. The violations of this condition, however, are widely recognized in many situations, especially in the flow of non-Newtonian fluids. Wall slip could lead to large errors and flow instabilities, such as sharkskin formation and spurt flow, and hence complicates the analysis of fluid systems and introduces serious practical difficulties. In this article, we discuss slip at fluid-solid interface in an attempt to highlight the main issues related to this diverse complex phenomenon and its implications.

preprint2010arXiv

A computed line list for the H2D+ molecular ion

A comprehensive, calculated line list of frequencies and transition probabilities for the singly deuterated isotopologue of H3+, H2D+, is presented. The line list, called ST1, contains over 22 million rotational-vibrational transitions occurring between more than 33 thousand energy levels; it covers frequencies up to 18500 cm-1. All energy levels with rotational quantum number, J, up to 20 are considered, making the line list useful for temperatures up to at least 3000 K. About 15% of these levels are fully assigned with approximate rotational and vibrational quantum numbers. The list is calculated using a previously proposed, high accuracy, ab initio model and consistency checks are carried out to test and validate the results. These checks confirm the accuracy of the list. A temperature-dependent partition function, valid over a more extended temperature range than those previously published, and cooling function are presented. Temperature-dependent synthetic spectra in the frequency range 0 - 10000 cm-1 are also given.

preprint2010arXiv

Computational Techniques for Efficient Conversion of Image Files from Area Detectors

Area detectors are used in many scientific and technological applications such as particle and radiation physics. Thanks to the recent technological developments, the radiation sources are becoming increasingly brighter and the detectors become faster and more efficient. The result is a sharp increase in the size of data collected in a typical experiment. This situation imposes a bottleneck on data processing capabilities, and could pose a real challenge to scientific research in certain areas. This article proposes a number of simple techniques to facilitate rapid and efficient extraction of data obtained from these detectors. These techniques are successfully implemented and tested in a computer program to deal with the extraction of X-ray diffraction patterns from EDF image files obtained from CCD detectors.

preprint2010arXiv

High Throughput Software for Powder Diffraction and its Application to Heterogeneous Catalysis

In this thesis we investigate high throughput computational methods for processing large quantities of data collected from synchrotrons and their application to spectral analysis of powder diffraction data. We also present the main product of this PhD programme, specifically a software called 'EasyDD' developed by the author. This software was created to meet the increasing demand on data processing and analysis capabilities as required by modern detectors which produce huge quantities of data. Modern detectors coupled with the high intensity X-ray sources available at synchrotrons have led to the situation where datasets can be collected in ever shorter time scales and in ever larger numbers. Such large volumes of datasets pose a data processing bottleneck which augments with current and future instrument development. EasyDD has achieved its objectives and made significant contributions to scientific research. It can also be used as a model for more mature attempts in the future. EasyDD is currently in use by a number of researchers in a number of academic and research institutions to process high-energy diffraction data. These include data collected by different techniques such as Energy Dispersive Diffraction, Angle Dispersive Diffraction and Computer Aided Tomography. EasyDD has already been used in a number of published studies, and is currently in use by the High Energy X-Ray Imaging Technology project. The software was also used by the author to process and analyse datasets collected from synchrotron radiation facilities. In this regard, the thesis presents novel scientific research involving the use of EasyDD to handle large diffraction datasets in the study of alumina-supported metal oxide catalyst bodies. These data were collected using Tomographic Energy Dispersive Diffraction Imaging and Computer Aided Tomography techniques.

preprint2010arXiv

Pore-Scale Modeling of Non-Newtonian Flow in Porous Media

The thesis investigates the flow of non-Newtonian fluids in porous media using pore-scale network modeling. Non-Newtonian fluids show very complex time and strain dependent behavior and may have initial yield stress. Their common feature is that they do not obey the simple Newtonian relation of proportionality between stress and rate of deformation. They are generally classified into three main categories: time-independent, time-dependent and viscoelastic. Two three-dimensional networks representing a sand pack and Berea sandstone were used. An iterative numerical technique is used to solve the pressure field and obtain the flow rate and apparent viscosity. The time-independent category is investigated using two fluid models: Ellis and Herschel-Bulkley. The analysis confirmed the reliability of the non-Newtonian network model used in this study. Good results are obtained, especially for the Ellis model, when comparing the network model results to experimental data sets found in the literature. The yield-stress phenomenon is also investigated and several numerical algorithms were developed and implemented to predict threshold yield pressure of the network. An extensive literature survey and investigation were carried out to understand the phenomenon of viscoelasticity with special attention to the flow in porous media. The extensional flow and viscosity and converging-diverging geometry were thoroughly examined as the basis of the peculiar viscoelastic behavior in porous media. The modified Bautista-Manero model was identified as a promising candidate for modeling the flow of viscoelastic materials which also show thixotropic attributes. An algorithm that employs this model was implemented in the non-Newtonian code and the initial results were analyzed. The time-dependent category was examined and several problems in modeling and simulating the flow of these fluids were identified.

preprint2010arXiv

The Flow of Newtonian Fluids in Axisymmetric Corrugated Tubes

This article deals with the flow of Newtonian fluids through axially-symmetric corrugated tubes. An analytical method to derive the relation between volumetric flow rate and pressure drop in laminar flow regimes is presented and applied to a number of simple tube geometries of converging-diverging nature. The method is general in terms of fluid and tube shape within the previous restrictions. Moreover, it can be used as a basis for numerical integration where analytical relations cannot be obtained due to mathematical difficulties.

preprint2010arXiv

The Flow of Power-Law Fluids in Axisymmetric Corrugated Tubes

In this article we present an analytical method for deriving the relationship between the pressure drop and flow rate in laminar flow regimes, and apply it to the flow of power-law fluids through axially-symmetric corrugated tubes. The method, which is general with regards to fluid and tube shape within certain restrictions, can also be used as a foundation for numerical integration where analytical expressions are hard to obtain due to mathematical or practical complexities. Five converging-diverging geometries are used as examples to illustrate the application of this method.

preprint2009arXiv

Emissivity: A Program for Atomic Emissivity Calculations

In this article we report the release of a new program for calculating the emissivity of atomic transitions. The program, which can be obtained with its documentation from our website www.scienceware.net, passed various rigorous tests and was used by the author to generate theoretical data and analyze observational data. It is particularly useful for investigating atomic transition lines in astronomical context as the program is capable of generating a huge amount of theoretical data and comparing it to observational list of lines. A number of atomic transition algorithms and analytical techniques are implemented within the program and can be very useful in various situations. The program can be described as fast and efficient. Moreover, it requires modest computational resources.

preprint2009arXiv

Modeling the Flow of Yield-Stress Fluids in Porous Media

Yield-stress is a problematic and controversial non-Newtonian flow phenomenon. In this article, we investigate the flow of yield-stress substances through porous media within the framework of pore-scale network modeling. We also investigate the validity of the Minimum Threshold Path (MTP) algorithms to predict the pressure yield point of a network depicting random or regular porous media. Percolation theory as a basis for predicting the yield point of a network is briefly presented and assessed. In the course of this study, a yield-stress flow simulation model alongside several numerical algorithms related to yield-stress in porous media were developed, implemented and assessed. The general conclusion is that modeling the flow of yield-stress fluids in porous media is too difficult and problematic. More fundamental modeling strategies are required to tackle this problem in the future.

preprint2009arXiv

Single-Phase Flow of Non-Newtonian Fluids in Porous Media

The study of flow of non-Newtonian fluids in porous media is very important and serves a wide variety of practical applications in processes such as enhanced oil recovery from underground reservoirs, filtration of polymer solutions and soil remediation through the removal of liquid pollutants. These fluids occur in diverse natural and synthetic forms and can be regarded as the rule rather than the exception. They show very complex strain and time dependent behavior and may have initial yield-stress. Their common feature is that they do not obey the simple Newtonian relation of proportionality between stress and rate of deformation. Non-Newtonian fluids are generally classified into three main categories: time-independent whose strain rate solely depends on the instantaneous stress, time-dependent whose strain rate is a function of both magnitude and duration of the applied stress and viscoelastic which shows partial elastic recovery on removal of the deforming stress and usually demonstrates both time and strain dependency. In this article the key aspects of these fluids are reviewed with particular emphasis on single-phase flow through porous media. The four main approaches for describing the flow in porous media are examined and assessed. These are: continuum models, bundle of tubes models, numerical methods and pore-scale network modeling.