Pub Date : 2019-03-31DOI: 10.2495/CMEM-V7-N2-130-141
R. Ravendran, B. Endelt, J. Christiansen, P. Jensen, Martin Theile, I. Najjar
Understanding the disturbances introduced by cavitation inside spray nozzles is important, when simulating the spray formation of highly viscous liquids. In this paper, a new model for cavitation-induced primary break-up is proposed, which is able to map the influence of cavitating nozzle flow on spray formation. Detailed experimental and numerical investigations of the viscous nozzle flow have been performed in order to develop an improved primary break-up model [1]. The proposed model describes the transition from the flow inside the nozzle, modelled using a homogeneous equilibrium model (HEM) method, to the first primary droplets modelled using a Eulerian–Lagrangian method. Thus, providing the boundary conditions for the calculation of the secondary break-up and spray formation. The nozzle exit is divided into a definite number of patches. Liquid momentum and vapor volume fraction from each patch are used to initialize the primary droplets. The model has been implemented in the open-source CFD software package OpenFOAM and validation has been done using high-speed shadow graphic imaging. The simulated spray tip penetration and spray cone angle at the near-nozzle region show a good agreement with the experiment results.
{"title":"Coupling method for internal nozzle flow and the spray formation for viscous liquids","authors":"R. Ravendran, B. Endelt, J. Christiansen, P. Jensen, Martin Theile, I. Najjar","doi":"10.2495/CMEM-V7-N2-130-141","DOIUrl":"https://doi.org/10.2495/CMEM-V7-N2-130-141","url":null,"abstract":"Understanding the disturbances introduced by cavitation inside spray nozzles is important, when simulating the spray formation of highly viscous liquids. In this paper, a new model for cavitation-induced primary break-up is proposed, which is able to map the influence of cavitating nozzle flow on spray formation. Detailed experimental and numerical investigations of the viscous nozzle flow have been performed in order to develop an improved primary break-up model [1]. The proposed model describes the transition from the flow inside the nozzle, modelled using a homogeneous equilibrium model (HEM) method, to the first primary droplets modelled using a Eulerian–Lagrangian method. Thus, providing the boundary conditions for the calculation of the secondary break-up and spray formation. The nozzle exit is divided into a definite number of patches. Liquid momentum and vapor volume fraction from each patch are used to initialize the primary droplets. The model has been implemented in the open-source CFD software package OpenFOAM and validation has been done using high-speed shadow graphic imaging. The simulated spray tip penetration and spray cone angle at the near-nozzle region show a good agreement with the experiment results.","PeriodicalId":36958,"journal":{"name":"International Journal of Computational Methods and Experimental Measurements","volume":"11 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-03-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81457387","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2019-03-31DOI: 10.2495/CMEM-V7-N2-142-153
Shanon M. Reckinger, Thomas H. Gibson, Fred Hohman, T. Morrison, Scott J. Reckinger, M. Carvalho
Overflows in the ocean occur when dense water flows down a continental slope into less dense ambient water. It is important to study idealized and small-scale models, which allow for confidence and control of parameters. The work presented here is a direct qualitative and quantitative comparison between physical laboratory experiments and lab-scale numerical simulations. Physical parameters are varied, including the Coriolis parameter, the inflow density, and the inflow volumetric flow rate. Laboratory experiments are conducted using a rotating square tank and high-resolution camera mounted on the table in the rotating reference frame. Video results are digitized in order to compare directly to numerical simulations. The MIT General Circulation Model (MITgcm), a three-dimensional ocean model, is used for the direct numerical simulations corresponding to the specific laboratory experiments. It was found that the MITgcm was not a good match to laboratory experiments when physical parameters fell within the high eddy activity regime. However, a more extensive resolution study is needed to understand this fully. The MITgcm simulations did provide a good qualitative and quantitative match to laboratory experiments run in a low eddy activity regime. In all cases, the MITgcm simulations had more eddy activity than the laboratory experiments.
{"title":"The effect of numerical parameters on eddies in oceanic overflows: a laboratory and numerical study","authors":"Shanon M. Reckinger, Thomas H. Gibson, Fred Hohman, T. Morrison, Scott J. Reckinger, M. Carvalho","doi":"10.2495/CMEM-V7-N2-142-153","DOIUrl":"https://doi.org/10.2495/CMEM-V7-N2-142-153","url":null,"abstract":"Overflows in the ocean occur when dense water flows down a continental slope into less dense ambient water. It is important to study idealized and small-scale models, which allow for confidence and control of parameters. The work presented here is a direct qualitative and quantitative comparison between physical laboratory experiments and lab-scale numerical simulations. Physical parameters are varied, including the Coriolis parameter, the inflow density, and the inflow volumetric flow rate. Laboratory experiments are conducted using a rotating square tank and high-resolution camera mounted on the table in the rotating reference frame. Video results are digitized in order to compare directly to numerical simulations. The MIT General Circulation Model (MITgcm), a three-dimensional ocean model, is used for the direct numerical simulations corresponding to the specific laboratory experiments. It was found that the MITgcm was not a good match to laboratory experiments when physical parameters fell within the high eddy activity regime. However, a more extensive resolution study is needed to understand this fully. The MITgcm simulations did provide a good qualitative and quantitative match to laboratory experiments run in a low eddy activity regime. In all cases, the MITgcm simulations had more eddy activity than the laboratory experiments.","PeriodicalId":36958,"journal":{"name":"International Journal of Computational Methods and Experimental Measurements","volume":"8 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-03-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78553502","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2018-07-01DOI: 10.2495/CMEM-V7-N1-57-67
S. D’Alessio
This research investigates the unsteady free convective flow of a viscous incompressible fluid from a differentially heated rotating sphere. The flow is assumed to remain laminar and to possess equatorial and azimuthal symmetry. The governing Navier-Stokes and energy equations are posed in terms of a scaled stream function vorticity formulation and are solved subject to no-slip and specified surface temperature conditions. At t = 0 an impulsive heat flux is applied in the form of a jump in surface temperature. An asymptotic solution valid for large Grashof numbers and small times following the impulsive startup is constructed. Two small parameters have been identified and based on this the flow variables are expanded in a double series in powers of these parameters. The non-zero leading-order terms in the asymptotic expansions have been determined analytically and the corresponding heat transfer coefficient has been found. Future work will involve obtaining numerical solutions.
{"title":"An analytical study of the early stages of unsteady free convective flow from a differentially heated rotating sphere at large Grashof numbers","authors":"S. D’Alessio","doi":"10.2495/CMEM-V7-N1-57-67","DOIUrl":"https://doi.org/10.2495/CMEM-V7-N1-57-67","url":null,"abstract":"This research investigates the unsteady free convective flow of a viscous incompressible fluid from a differentially heated rotating sphere. The flow is assumed to remain laminar and to possess equatorial and azimuthal symmetry. The governing Navier-Stokes and energy equations are posed in terms of a scaled stream function vorticity formulation and are solved subject to no-slip and specified surface temperature conditions. At t = 0 an impulsive heat flux is applied in the form of a jump in surface temperature. An asymptotic solution valid for large Grashof numbers and small times following the impulsive startup is constructed. Two small parameters have been identified and based on this the flow variables are expanded in a double series in powers of these parameters. The non-zero leading-order terms in the asymptotic expansions have been determined analytically and the corresponding heat transfer coefficient has been found. Future work will involve obtaining numerical solutions.","PeriodicalId":36958,"journal":{"name":"International Journal of Computational Methods and Experimental Measurements","volume":"61 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2018-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86821239","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2018-07-01DOI: 10.2495/CMEM-V7-N1-68-78
T. Moulden
Vortical structures have been observed to develop in electrically driven fluid motion at the micro/nano scale, but no coherent theory has been put foreword in the literature to explain such a development. The present paper gives several results in a theory based upon the classical field equations. In particular, it is shown that the origin of vorticity production resides in the applied electric field interacting with any ion concentration gradients present in the fluid as defined by the vorticity equation. This is in addition to any viscous layer vorticity diffusion that may also exist in the flow.
{"title":"Micro/nano flows: vorticity generation","authors":"T. Moulden","doi":"10.2495/CMEM-V7-N1-68-78","DOIUrl":"https://doi.org/10.2495/CMEM-V7-N1-68-78","url":null,"abstract":"Vortical structures have been observed to develop in electrically driven fluid motion at the micro/nano scale, but no coherent theory has been put foreword in the literature to explain such a development. The present paper gives several results in a theory based upon the classical field equations. In particular, it is shown that the origin of vorticity production resides in the applied electric field interacting with any ion concentration gradients present in the fluid as defined by the vorticity equation. This is in addition to any viscous layer vorticity diffusion that may also exist in the flow.","PeriodicalId":36958,"journal":{"name":"International Journal of Computational Methods and Experimental Measurements","volume":"12 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2018-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88013945","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2018-07-01DOI: 10.2495/CMEM-V7-N1-14-21
Salima Djehaichia, R. Lassoued
The vulnerability of reinforced concrete structures, which were built in the 1970s, under the effects of an earthquake is one of the major concerns of researchers, because the designs of these structures have been based on regulations, which are no longer valid. The parameters taken into account in this study to idealize the regulatory shortcomings are: low ratio of reinforcement, type of reinforcement and moderate resistance of concrete. The analysis to test these altered structures with one or both of the above parameters is carried out in the non-linear domain from the perspective of analysing their behaviour in an earthquake. In this paper, the modelling strategy is based on finite elements combined with a discretization of the shear wall most stressed by successive thin layers. The estimation of level of performance is achieved using capacity curves under increasing incremental loads; a non-linear characteristic force-displacement relationship can be determined. The results of the numerical model are compared with those of the Algerian seismic code (RPA). Through this comparison, it was found that there is an improvement in terms of displacement, shearing action and ductility. The introduction of confining as a local model makes it possible to refine the numerical model and improve the total behaviour of the structure. A parametric analysis is carried out to highlight the obvious weakness of structures designed and built in the 1970s.
{"title":"Non-linear behaviour of structural walls","authors":"Salima Djehaichia, R. Lassoued","doi":"10.2495/CMEM-V7-N1-14-21","DOIUrl":"https://doi.org/10.2495/CMEM-V7-N1-14-21","url":null,"abstract":"The vulnerability of reinforced concrete structures, which were built in the 1970s, under the effects of an earthquake is one of the major concerns of researchers, because the designs of these structures have been based on regulations, which are no longer valid. The parameters taken into account in this study to idealize the regulatory shortcomings are: low ratio of reinforcement, type of reinforcement and moderate resistance of concrete. The analysis to test these altered structures with one or both of the above parameters is carried out in the non-linear domain from the perspective of analysing their behaviour in an earthquake. In this paper, the modelling strategy is based on finite elements combined with a discretization of the shear wall most stressed by successive thin layers. The estimation of level of performance is achieved using capacity curves under increasing incremental loads; a non-linear characteristic force-displacement relationship can be determined. The results of the numerical model are compared with those of the Algerian seismic code (RPA). Through this comparison, it was found that there is an improvement in terms of displacement, shearing action and ductility. The introduction of confining as a local model makes it possible to refine the numerical model and improve the total behaviour of the structure. A parametric analysis is carried out to highlight the obvious weakness of structures designed and built in the 1970s.","PeriodicalId":36958,"journal":{"name":"International Journal of Computational Methods and Experimental Measurements","volume":"148 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2018-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"85363314","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2018-07-01DOI: 10.2495/CMEM-V7-N1-45-56
B. Wilkins, T. Hromadka, A. Johnson, Randy Boucher, H. D. McInvale, S. Horton
Solving potential problems, such as those that occur in the analysis of steady-state heat transfer, electrostatics, ideal fluid flow, and groundwater flow, is important in several fields of engineering, science, and applied mathematics. Numerical solution of the relevant governing equations typically involves using techniques such as domain methods (including finite element, finite difference, or finite volume), or boundary element methods (using either real or complex variables). In this paper, the Complex Variable Boundary Element method (“CVBEM”) is examined with respect to the use of different types of basis functions in the CVBEM approximation function. Four basis function families are assessed in their solution success in modeling an important benchmark problem in ideal fluid flow; namely, flow around a 90 degree bend. Identical problem domains are used in the examination, and identical degrees of freedom are used in the CVBEM approximation functions. Further, a new computational modeling error is defined and used to compare the results herein; specifically, M = E / N where M is the proposed computational error measure, E is the maximum difference (in absolute value) between approximation and boundary condition value, and N is the number of degrees of freedom used in the approximation.
解决潜在的问题，例如在分析稳态传热、静电学、理想流体流动和地下水流动时出现的问题，在工程、科学和应用数学的几个领域都是重要的。相关控制方程的数值解通常涉及使用诸如域方法(包括有限元、有限差分或有限体积)或边界元方法(使用实变量或复变量)等技术。本文研究了复杂变量边界元法(CVBEM)在CVBEM近似函数中不同类型基函数的使用。评价了四个基函数族对理想流体流动中一个重要基准问题的求解成功程度;也就是说，绕90度弯道流动。在检验中使用了相同的问题域，并且在CVBEM近似函数中使用了相同的自由度。此外，定义了一个新的计算建模误差，并将其用于比较本文的结果;其中M = E / N，其中M为提出的计算误差度量，E为近似值与边界条件值之间的最大差值(绝对值)，N为近似值所使用的自由度数。
{"title":"Assessment of complex variable basis functions in the approximation of ideal fluid flow problems","authors":"B. Wilkins, T. Hromadka, A. Johnson, Randy Boucher, H. D. McInvale, S. Horton","doi":"10.2495/CMEM-V7-N1-45-56","DOIUrl":"https://doi.org/10.2495/CMEM-V7-N1-45-56","url":null,"abstract":"Solving potential problems, such as those that occur in the analysis of steady-state heat transfer, electrostatics, ideal fluid flow, and groundwater flow, is important in several fields of engineering, science, and applied mathematics. Numerical solution of the relevant governing equations typically involves using techniques such as domain methods (including finite element, finite difference, or finite volume), or boundary element methods (using either real or complex variables). In this paper, the Complex Variable Boundary Element method (“CVBEM”) is examined with respect to the use of different types of basis functions in the CVBEM approximation function. Four basis function families are assessed in their solution success in modeling an important benchmark problem in ideal fluid flow; namely, flow around a 90 degree bend. Identical problem domains are used in the examination, and identical degrees of freedom are used in the CVBEM approximation functions. Further, a new computational modeling error is defined and used to compare the results herein; specifically, M = E / N where M is the proposed computational error measure, E is the maximum difference (in absolute value) between approximation and boundary condition value, and N is the number of degrees of freedom used in the approximation.","PeriodicalId":36958,"journal":{"name":"International Journal of Computational Methods and Experimental Measurements","volume":"9 7 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2018-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74494246","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
N. Demoes, G. Bann, B. Wilkins, T. Hromadka, Randy Boucher
The Complex Variable Boundary Element Method, or CVBEM, was first published in Journal of Numerical Methods in Engineering in year 1984 by authors Hromadka and Guymon [1]. Since that time, several papers and books have been published that present various aspects of the numerical technique as well as advances in the computational method such as extension to three or higher dimensions for arbitrary geometries, nonhomogeneous domains, extension to use of a Hilbert Space setting as well as collocation methods, inclusion of the time derivative via coupling to generalized Fourier series techniques, examination of various families of basis functions including complex monomials, the product of complex polynomials with complex logarithm functions (i.e., the usual CVBEM basis functions), Laurent series expansions, reciprocal of complex monomials, other complex variable analytic functions including exponential and others, as well as linear combinations of these families. Other topics studied and developed include rotation of complex logarithm branch-cuts for extension of the problem computational domain to the exterior of the problem geometry, depiction of computational error in achieving problem boundary conditions by means of the approximate boundary technique, mixed boundary value problems, flow net development and visualization, display of flow field trajectory vectors in two and three dimensions for use in depicting streamlines and flow paths, among other topics. The CVBEM approach has also been extended to solving partial differential equations such as Laplace’s equation, Poisson’s equation, unsteady flow equation, and the wave equation, among other formulations that include sources, sinks and combinations of these equations with mixed boundary conditions. In the current paper, a detailed examination is made of the performance between four families of basis functions in order to assess computational efficiency in problem solving of two dimensional potential problems in a high aspect ratio geometric problem domain. Two selected problems are presented as case studies to demonstrate the different levels of success by each of the four families of examined basis functions. All four families involve basis functions that solve the governing partial differential equation, leaving only the goodness of fit in matching boundary conditions of the boundary value problem as the computational optimization goal. The modeling technique is implemented in computer programs Mathematica and MATLAB. Recommendations are made for future research directions and lessons learned from the current study effort.
复杂变量边界元法(CVBEM)由Hromadka和Guymon于1984年首次发表在《Journal of Numerical Methods in Engineering》上[1]。从那时起，已经出版了一些论文和书籍，介绍了数值技术的各个方面以及计算方法的进展，例如任意几何的三维或更高维度的扩展，非齐次域，希尔伯特空间设置的扩展以及搭配方法，通过与广义傅里叶级数技术的耦合包含时间导数，考察各种基函数族，包括复单项式，复多项式与复对数函数的乘积(即通常的CVBEM基函数)，洛朗级数展开式，复单项式的倒数，其他复变量解析函数，包括指数和其他，以及这些族的线性组合。研究和开发的其他主题包括复对数分支切割的旋转，将问题计算域扩展到问题几何的外部，通过近似边界技术描述实现问题边界条件的计算误差，混合边界值问题，流网的开发和可视化，用于描绘流线和流路的二维和三维流场轨迹向量的显示，在其他话题中。CVBEM方法还被扩展到求解偏微分方程，如拉普拉斯方程、泊松方程、非定常流动方程和波动方程，以及其他公式，包括源、槽和混合边界条件下这些方程的组合。在本文中，为了评估在高纵横比几何问题域中求解二维潜在问题的计算效率，对四种基函数族之间的性能进行了详细的研究。两个选定的问题作为案例研究提出，以证明四个检验基函数族中的每一个的不同程度的成功。这四个族都涉及求解控制偏微分方程的基函数，只留下边值问题的匹配边界条件的拟合优度作为计算优化目标。建模技术在计算机程序Mathematica和MATLAB中实现。最后对今后的研究方向提出了建议，并从目前的研究工作中总结了一些经验教训。
{"title":"35 years of advancements with the complex variable boundary element method","authors":"N. Demoes, G. Bann, B. Wilkins, T. Hromadka, Randy Boucher","doi":"10.2495/CMEM-V7-N1-1-13","DOIUrl":"https://doi.org/10.2495/CMEM-V7-N1-1-13","url":null,"abstract":"The Complex Variable Boundary Element Method, or CVBEM, was first published in Journal of Numerical Methods in Engineering in year 1984 by authors Hromadka and Guymon [1]. Since that time, several papers and books have been published that present various aspects of the numerical technique as well as advances in the computational method such as extension to three or higher dimensions for arbitrary geometries, nonhomogeneous domains, extension to use of a Hilbert Space setting as well as collocation methods, inclusion of the time derivative via coupling to generalized Fourier series techniques, examination of various families of basis functions including complex monomials, the product of complex polynomials with complex logarithm functions (i.e., the usual CVBEM basis functions), Laurent series expansions, reciprocal of complex monomials, other complex variable analytic functions including exponential and others, as well as linear combinations of these families. Other topics studied and developed include rotation of complex logarithm branch-cuts for extension of the problem computational domain to the exterior of the problem geometry, depiction of computational error in achieving problem boundary conditions by means of the approximate boundary technique, mixed boundary value problems, flow net development and visualization, display of flow field trajectory vectors in two and three dimensions for use in depicting streamlines and flow paths, among other topics. The CVBEM approach has also been extended to solving partial differential equations such as Laplace’s equation, Poisson’s equation, unsteady flow equation, and the wave equation, among other formulations that include sources, sinks and combinations of these equations with mixed boundary conditions. In the current paper, a detailed examination is made of the performance between four families of basis functions in order to assess computational efficiency in problem solving of two dimensional potential problems in a high aspect ratio geometric problem domain. Two selected problems are presented as case studies to demonstrate the different levels of success by each of the four families of examined basis functions. All four families involve basis functions that solve the governing partial differential equation, leaving only the goodness of fit in matching boundary conditions of the boundary value problem as the computational optimization goal. The modeling technique is implemented in computer programs Mathematica and MATLAB. Recommendations are made for future research directions and lessons learned from the current study effort.","PeriodicalId":36958,"journal":{"name":"International Journal of Computational Methods and Experimental Measurements","volume":"28 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2018-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"89996367","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2018-07-01DOI: 10.2495/CMEM-V7-N1-33-44
Rasel A. Sultan, Serag Alfarek, M. A. Rahman, S. Zendehboudi
This study analyses three dimensional fluid flow through horizontal pipelines with three-phase gasliquid-solid Newtonian fluids by Computational Fluid Dynamics (CFD) simulation. Validating the simulation with experimental data, the study aims to develop a versatile acceptable simulation model that can be used further for different applied cases. An experimental setup is developed in our laboratory to determine slug flow (air-water) through a horizontal pipeline. Air as gas, water as liquid and silica as solid particle is used in this work. ANSYS Fluent version 16.2 is employed to perform the simulation. The Eulerian multiphase model with the Reynolds Stress Model (RSM) turbulence closure is adopted to analyse multiphase fluid flow. Parameters are selected from experimental works to validate the simulation. After a good agreement with experimental data, sensitivity analysis is conducted to observe the three phase fluid flow characteristics through horizontal flow. Pressure gradient (pressure drop per unit length) and in situ concentration profile are used as primary parameters. This article provides a clear relationship between the different parameters of three-phase fluid flow through a horizontal pipeline.
{"title":"CFD and experimental approach on three phase gas-liquid-solid Newtonian fluid flow in horizontal pipes","authors":"Rasel A. Sultan, Serag Alfarek, M. A. Rahman, S. Zendehboudi","doi":"10.2495/CMEM-V7-N1-33-44","DOIUrl":"https://doi.org/10.2495/CMEM-V7-N1-33-44","url":null,"abstract":"This study analyses three dimensional fluid flow through horizontal pipelines with three-phase gasliquid-solid Newtonian fluids by Computational Fluid Dynamics (CFD) simulation. Validating the simulation with experimental data, the study aims to develop a versatile acceptable simulation model that can be used further for different applied cases. An experimental setup is developed in our laboratory to determine slug flow (air-water) through a horizontal pipeline. Air as gas, water as liquid and silica as solid particle is used in this work. ANSYS Fluent version 16.2 is employed to perform the simulation. The Eulerian multiphase model with the Reynolds Stress Model (RSM) turbulence closure is adopted to analyse multiphase fluid flow. Parameters are selected from experimental works to validate the simulation. After a good agreement with experimental data, sensitivity analysis is conducted to observe the three phase fluid flow characteristics through horizontal flow. Pressure gradient (pressure drop per unit length) and in situ concentration profile are used as primary parameters. This article provides a clear relationship between the different parameters of three-phase fluid flow through a horizontal pipeline.","PeriodicalId":36958,"journal":{"name":"International Journal of Computational Methods and Experimental Measurements","volume":"68 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2018-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80053331","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2018-07-01DOI: 10.2495/CMEM-V7-N1-22-32
A. Malesińska, M. Rogulski, P. Puntorieri, G. Barbaro, B. Kowalska
Pipe lines are useful for transporting water for drinking, irrigation and for fire-ing over long distances, this pipe lines are called “Transmission line” and are used to carry conveying raw or treated water from a well field or remote storage (large lake, reservoir, etc.,) facility to a treatment plant and/or distribution storage tank. In water-carrying piping systems, dangerous phenomena may occur. One such phenomenon is water hammer. The water hammer has always been an area of study, which has captivated the minds of researchers due to its complex and challenging phenomena. Modeling the phenomenon in real conditions is extremely difficult. Due to the dimensions of the piping systems, conducting research at real scales is impossible. However, thanks to the development of numerical methods, the study of water hammer and its effects can be performed using simulation programs. Unfortunately, the simulation results are not always consistent with the actual course of the phenomenon. One of the parameters that describes the nature of the course of a water hammer is the velocity of propagation of the pressure wave, c, which is called celerity. The transient surge pressure, p, may be calculated from the pressure celerity c, and the sudden change in fluid flow velocity, ∆ v. In a piping system, the value of the pressure wave celerity is not equal to the individual celerity, c, for a single pipeline. Therefore for piping systems for ∆p calculations the equivalent celerity shell be used. This article presents value of the equivalent celerity calculated from equations derived using linear analysis of natural vibrations of the system. For implement of the equations, an algorithm in MATLAB has been developed that allows one to easily calculate the equivalent celerity, ce, for N pipelines connected in series with varying diameter, length and material composition.
{"title":"Use of equivalent celerity to estimate maximum pressure increase in serial pipes during water hammer - numerical simulations in MATLAB","authors":"A. Malesińska, M. Rogulski, P. Puntorieri, G. Barbaro, B. Kowalska","doi":"10.2495/CMEM-V7-N1-22-32","DOIUrl":"https://doi.org/10.2495/CMEM-V7-N1-22-32","url":null,"abstract":"Pipe lines are useful for transporting water for drinking, irrigation and for fire-ing over long distances, this pipe lines are called “Transmission line” and are used to carry conveying raw or treated water from a well field or remote storage (large lake, reservoir, etc.,) facility to a treatment plant and/or distribution storage tank. In water-carrying piping systems, dangerous phenomena may occur. One such phenomenon is water hammer. The water hammer has always been an area of study, which has captivated the minds of researchers due to its complex and challenging phenomena. Modeling the phenomenon in real conditions is extremely difficult. Due to the dimensions of the piping systems, conducting research at real scales is impossible. However, thanks to the development of numerical methods, the study of water hammer and its effects can be performed using simulation programs. Unfortunately, the simulation results are not always consistent with the actual course of the phenomenon. One of the parameters that describes the nature of the course of a water hammer is the velocity of propagation of the pressure wave, c, which is called celerity. The transient surge pressure, p, may be calculated from the pressure celerity c, and the sudden change in fluid flow velocity, ∆ v. In a piping system, the value of the pressure wave celerity is not equal to the individual celerity, c, for a single pipeline. Therefore for piping systems for ∆p calculations the equivalent celerity shell be used. This article presents value of the equivalent celerity calculated from equations derived using linear analysis of natural vibrations of the system. For implement of the equations, an algorithm in MATLAB has been developed that allows one to easily calculate the equivalent celerity, ce, for N pipelines connected in series with varying diameter, length and material composition.","PeriodicalId":36958,"journal":{"name":"International Journal of Computational Methods and Experimental Measurements","volume":"27 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2018-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"73798067","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：