Non-Orthogonal Transformation of Irregular Geometry for Particle Based Simulation

Author(s):  
Anurag Kumar ◽  
Eiyad Abu-Nada ◽  
Toru Yamada ◽  
Yutako Asako ◽  
Mohammad Faghri

Simulations of irregular geometries using non-orthogonal transformation is widely used in grid based methodology such as computational fluid dynamics. However, this approach is not utilized for particle based models. In this paper we introduce non-orthogonal transformation to simulate fluid flow in irregular geometry using dissipative particle dynamics (DPD). Applying boundary condition is not trivial in DPD methodology and problem becomes more complicated for irregular boundary. In the present work, irregular (physical) domain is transformed into a rectangular domain and boundary particles are frozen along the wall. Transformation for position and velocity is used to relate physical and computational domains. As particle’s position and velocity change with time, transformation matrices are determined for each DPD particle at every time step. In DPD, forces are function of actual distance between the particles and acts within a cutoff radius, which change in transformed domain at every location. To solve this problem, firstly, interacting particles are identified in the physical domain and then forces are calculated in the transformed domain. This approach is described by simulating fluid flow inside a convergent-divergent nozzle, whose geometry is controlled by the contraction ratio (CR) in the middle of the nozzle. The DPD results were validated against in-house computational fluid dynamic (CFD) finite volume code based on the stream function vorticity approach. The range of Reynolds number and CR, under study here, is Re = 10–200 and CR = 0.8 and 0.6, respectively. The results revealed an excellent agreement between the DPD and CFD. The maximum deviation between the DPD and CFD results is within 2%. It is found that using large values of dissipative force parameter velocity fluctuations are less.

2014 ◽  
Vol 493 ◽  
pp. 80-85 ◽  
Author(s):  
C.L Siow ◽  
Jaswar ◽  
Efi Afrizal

Computational Fluid Dynamics (CFD) software is often used to study fluid flow and structures motion in fluids. The CFD normally requires large size of arrays and computer memory and then caused long execution time. However, Innovation of computer hardware such as multi-cores processor provides an alternative solution to improve this programming performance. This paper discussed loop parallelize multi-cores processor for optimization of sequential looping CFD code. This loop parallelize CFD was achieved by applying multi-tasking or multi-threading code into the original CFD code which was developed by one of the authors. The CFD code was developed based on Reynolds Average Navier-Stokes (RANS) method. The new CFD code program was developed using Microsoft Visual Basic (VB) programming language. In the early stage, the whole CFD code was constructed in a sequential flow before it is modified to parallel flow by using VBs multi-threading library. In the comparison, fluid flow around the hull of round-shaped FPSO was selected to compare the performance of both the programming codes. Besides, executed results of this self-developed code such as pressure distribution around the hull were also presented in this paper.


Author(s):  
M. Benaouicha ◽  
S. Guillou ◽  
A. Santa Cruz ◽  
H. Trigui

The study deals with a 3D Fluid-Structure Interaction (FSI) numerical model of a rectangular cantilevered flexible hydrofoil subjected to a turbulent fluid flow regime. The structural response and dynamic deformations are studied by analyzing the oscillations frequencies and amplitudes, under a hydrodynamics loads. The obtained numerical results are confronted with experimental ones, for validation. The numerical model is performed in the same geometric, physical and material conditions as the experimental set-up carried out in a hydrodynamic tunnel. A polyacetal (POM) flexible hydrofoil NACA0015 with an angle of attack of 8° is considered to be immersed in a fluid flow at a Reynold number of 3 × 105. The structure is initially at rest and then moved by the action of the fluid flow. The numerical model is based on a strong coupling procedure for solving the Fluid-Structure Interaction problem. The Arbitrary Lagrangian-Eulerian (ALE) formulation of the Navier-Stokes equations is used and an anisotropic diffusion equation is solved to compute the fluid mesh velocity and position at each time step. The finite volume method is used for the numerical resolution of the fluid dynamics equations. The structure deformations are described by the linear elasticity equation which is solved by the finite elements method. The Fluid-Structure coupled problem is solved by using the partitioned FSI implicit algorithm. A good agreement between numerical and experimental results for the hydrodynamics coefficients and hydrofoil deformations, maximum deflection and frequencies is obtained. The added mass and damping are analyzed and then the FSI effect on the dynamic deformations of the structure is highlighted.


2016 ◽  
Vol 32 (6) ◽  
pp. 777-785
Author(s):  
M. Besanjideh ◽  
M. Hajabdollahi ◽  
S. A. Gandjalikhan Nassab

AbstractThis paper deals with studying fluid flow and heat transfer of nanofluid through a forward facing step channel which is affected by a uniform magnetic field transverse to fluid flow. All the channel walls are assumed to be in constant temperature and the fluid temperature at the channel inlet is less than that of the walls. Also, the nanofluid is considered as a single-phase Newtonian fluid and the proper correlations were utilized to determine the thermophysical properties of nanofluid. Therefore, a code has been developed and two-dimensional continuity, momentum and energy equations were solved, using CFD technique. The computations were conducted for different values of the Reynolds and Hartmann numbers, and contraction ratio and an extensive range of nanoparticles volume fraction. The results indicated that flow separation and reattachment phenomena, in vicinity of the step edge, could be influenced strongly by magnetic field and the average Nusselt number is increased significantly by increasing nanoparticles volume fraction and Hartmann number.


2000 ◽  
Author(s):  
Zhilong He ◽  
Xueyuan Peng ◽  
Pengcheng Shu

Abstract This paper presents a numerical method for simulating the thermal and fluid-dynamic behavior of hermetic compressors in the whole compressor domain. The model of fluid flow is developed by integrating transient one-dimensional conservation equations of continuity, momentum and energy through all of the elements from suction line to discharge line. The model describing thermal behavior is based on heat balance in the components such as muffler, connecting tubes and orifices. The calculation of the thermodynamic and transport properties for different refrigerants at various conditions has been considered, and some numerical results for a hermetic compressor are presented. The present study has demonstrated that the numerical simulation is a fest and reliable tool for compressor design.


2021 ◽  
Author(s):  
Amir A. Mofakham ◽  
Farid Rousta ◽  
Dustin M. Crandall ◽  
Goodarz Ahmadi

Abstract Hydraulic fracturing or fracking is a procedure used extensively by oil and gas companies to extract natural gas or petroleum from unconventional sources. During this process, a pressurized liquid is injected into wellbores to generate fractures in rock formations to create more permeable pathways in low permeability rocks that hold the oil. To keep the rock fractures open after removing the high pressure, proppant, which typically are sands with different shapes and sizes, are injected simultaneously with the fracking fluid to spread them throughout rock fractures. The extraction productivity from shale reservoirs is significantly affected by the performance and quality of the proppant injection process. Since these processes occur under the ground and in the rock fractures, using experimental investigations to examine the process is challenging, if not impossible. Therefore, employing numerical tools for analyzing the process could provide significant insights leading to the fracking process improvement. Accordingly, in this investigation, a 4-way coupled Computational Fluid Dynamic and Discrete Element Method (CFD-DEM) code was used to simulate proppant transport into a numerically generated realistic rock fracture geometry. The simulations were carried out for a sufficiently long period to reach the fractures’ steady coverage by proppant. The proppant fracture coverage is a distinguishing factor that can be used to assess the proppant injection process quality. A series of simulations with different proppant sizes as well as various fracking fluid flow rates, were performed. The corresponding estimated fracture coverages for different cases were compared. The importance of proppant size as well as the fluid flow rate on the efficiency of the proppant injection process, were evaluated and discussed.


2015 ◽  
Author(s):  
Po-Wen Wang ◽  
Chi-Fang Lee ◽  
Yann Quéméner ◽  
Chien-Hua Huang

The objective of this study was to clarify the theoretical basis of sloshing loads and required plate thickness formulations in the harmonized common structural rules. This study used computational fluid dynamic (CFD) to calculate sloshing loads and used finite element analyses (FEA) to evaluate structural response. The sensitivity of the CFD predictions to the time step and grid size was also investigated. Cargo oil tanks were then selected in a handy size oil tanker and a very large crude carrier to evaluate the longitudinal and transverse sloshing loads on the tank boundaries. The results showed that the sloshing pressures computed at four filling levels were mostly consistent with CSR-H. Afterward, the sloshing pressure produced by CFD was applied to the finite element model by using a fluid-structure interaction technique to obtain the dynamic response of the structure. The dynamic responses were investigated to validate the quasistatic approach for sloshing assessment.


Author(s):  
Clemens Bernhard Domnick ◽  
Friedrich-Karl Benra ◽  
Dieter Brillert ◽  
Hans Josef Dohmen ◽  
Christian Musch

The power output of steam turbines is controlled by steam turbine inlet valves. These valves have a large flow capacity and dissipate in throttled operation a huge amount of energy. Due to that, high dynamic forces occur in the valve which can cause undesired valve vibrations. In this paper, the structural dynamics of a valve are analysed. The dynamic steam forces obtained by previous computational fluid dynamic (CFD) calculations at different operating points are impressed on the structural dynamic finite element model (FEM) of the valve. Due to frictional forces at the piston rings and contact effects at the bushings of the valve plug and the valve stem the structural dynamic FEM is highly nonlinear and has to be solved in the time domain. Prior to the actual investigation grid and time step studies are carried out. Also the effect of the temperature distribution within the valve stem is discussed and the influence of the valve actuator on the vibrations is analysed. In the first step, the vibrations generated by the fluid forces are investigated. The effects of the piston rings on the structural dynamics are discussed. It is found, that the piston rings are able to reduce the vibration significantly by frictional damping. In the second step, the effect of the moving valve plug on the dynamic flow in the valve is analysed. The time dependent displacement of the valve is transferred to CFD calculations using deformable meshes. With this one way coupling method the response of the flow to the vibrations is analysed.


2020 ◽  
Vol 11 (1) ◽  
pp. 13
Author(s):  
Vahid Rezania ◽  
Dennis Coombe ◽  
Jack Tuszynski

Tissue engineering, with the goal of repairing or replacing damaged tissue and organs, has continued to make dramatic science-based advances since its origins in the late 1980’s and early 1990’s. Such advances are always multi-disciplinary in nature, from basic biology and chemistry through physics and mathematics to various engineering and computer fields. This review will focus its attention on two topics critical for tissue engineering liver development: (a) fluid flow, zonation, and drug screening, and (b) biomechanics, tissue stiffness, and fibrosis, all within the context of 3D structures. First, a general overview of various bioreactor designs developed to investigate fluid transport and tissue biomechanics is given. This includes a mention of computational fluid dynamic methods used to optimize and validate these designs. Thereafter, the perspective provided by computer simulations of flow, reactive transport, and biomechanics responses at the scale of the liver lobule and liver tissue is outlined, in addition to how bioreactor-measured properties can be utilized in these models. Here, the fundamental issues of tortuosity and upscaling are highlighted, as well as the role of disease and fibrosis in these issues. Some idealized simulations of the effects of fibrosis on lobule drug transport and mechanics responses are provided to further illustrate these concepts. This review concludes with an outline of some practical applications of tissue engineering advances and how efficient computational upscaling techniques, such as dual continuum modeling, might be used to quantify the transition of bioreactor results to the full liver scale.


1967 ◽  
Vol 7 (02) ◽  
pp. 205-220 ◽  
Author(s):  
H.W. Price ◽  
D.A.T. Donohue

Abstract The system of equations describing displacement of a hydrocarbon liquid by a hydrocarbon vapor in a porous medium where mass transfer takes place between the phases is solved numerically for a variety of gas injection processes. Even though the method of solution is quite general, only systems with three hydrocarbon components are considered. Computer simulations of displacement processes wherein mass transfer between phases is both considered and neglected are compared, and it is shown that neglecting mass transfer can give pessimistic displacement efficiencies. Introduction The role of the gas displacement process in the recovery of petroleum has been subjected to a series of detailed analyses; as a result, a number of predictive models have been published in the literature. However, because of major simplifying assumptions, most of these models do not completely represent the physical system. As a result, the effect of making the simplifying assumptions is unknown. Therefore, a complete representation of this process one without major simplifying assumptions should lead to a full understanding of the process, and perhaps to methods of improving it. The general method of developing a model for two-phase fluid flow in a porous medium is to solve simultaneously the continuity equation, the energy equations and the equation-of-state for each phase under the prescribed initial and boundary conditions. For an isothermal system, the energy equations reduce to the momentum equation, Darcy's law. However, since natural gas is the vapor state of the reservoir liquid, interphase mass transfer may take place with concomitant changes in both the intensive and extensive thermodynamic properties of each phase. It is this phenomenon that has often been omitted in previous mathematical models. An additional relation, then, which accounts for mass transfer between the phases, must be included with the other equations to specify a complete model. Completely formulating the equations to be solved is not a difficult task but obtaining their solution has been intractable up to now. Availability of large-memory, high-speed digital computers now makes an attack on this formidable problem possible. This paper presents a preliminary study of the problem. Since this investigation is intended to be exploratory, it is restricted to the linear, horizontal, isothermal, two-phase viscous flow of oil and gas in an oil reservoir. In the early development of predictive models of this process, the reservoir system was considered as a unit and various forms of the material balance equation were proposed. Pressure and saturation gradients were than added in the Buckley-Leverett model. The Buckley-Leverett formulation considered the fluids to be incompressible; thus, the mathematical model reduces to a steady-state system. In the 1950's, studies incorporating numerical techniques were being published. These mathematical models differed in the efficiency of finite difference techniques, the inclusion or exclusion of capillarity or the number of space dimensions considered. To solve these nonlinear, partial differential equations, each phase was considered to be homogeneous with time; therefore, mass transfer between phases was neglected. The effect of mass transfer on the gas displacement process was first reported by Attra. He simulated the one-dimension flow system by a series of cells in each of which the fluids were equilibrated during a time step. In addition, the pressure throughout the system during each time step was predetermined and constant phase velocities were calculated according to the Buckley-Leverett incompressible fluid flow model. Welge et al. developed a model for the displacement of oil by an enriched gas where composition is considered to be a dependent variable. SPEJ P. 205ˆ


1980 ◽  
Vol 20 (06) ◽  
pp. 533-554 ◽  
Author(s):  
Keith H. Coats

Abstract This paper describes a numerical model forsimulating wet or dry, forward or reverse combustionin one, two, or three dimensions. The formulation isconsiderably more general than any reported to date.The model allows any number and identities ofcomponents. Any component may be distributed inany or all of the four phases (water, oil, gas, andsolid or coke.The formulation allows any number of chemicalreactions. Any reaction may have any number ofreactants, products, and stoichiometry, identifiedthrough input data. The energy balance accounts forheat loss and conduction, conversion, and radiationwithin the reservoir.The model uses no assumptions regarding degreeof oxygen consumption. The oxygen concentration iscalculated throughout the reservoir in accordancewith the calculated fluid flow pattern and reactionkinetics. The model, therefore, simulates the effectsof oxygen bypassing caused by kinetic-limitedcombustion or conformance factors.We believe the implicit model formulation resultsin maximum efficiency (lowest computing cost), andrequired computing times are reported in the paper.The paper includes comparisons of model resultswith reported laboratory adiabatic-tube test results.In addition, the paper includes example field-scalecases, with a sensitivity study showing effects on oilrecovery of uncertainties in rock/fluid properties. Introduction Recent papers by Ali, Crookston et al., andYoungren provide a comprehensive review of earlierwork in numerical modeling of the in-situcombustion process.The trend in this modeling has been toward morerigorous treatment of the fluid flow and interphasemass transfer; inclusion of more components, morecomprehensive reaction kinetics, and stoichiometry;and more implicit treatment of the finite differencemodel equations.The purpose of this work was to extend thegenerality of previous models while preserving orreducing the associated computing-time requirement.The most comprehensive or sophisticated combustionmodels described to date appear to be thoseof Crookston et al. and Youngren. Therefore, wecompare our model formulation and results here withthose models.A common objective of different investigators'efforts in modeling in-situ combustion is developmentof more efficient formulations and methods ofsolution. This is especially important in thecombustion case because of the large number ofcomponents and equations involved. For a given numberof components and reactions, computing time pergrid block per time step will increase rapidly as theformulation is rendered more implicit. However, increasing implicitness tends to allow larger timesteps, which in turn reduces overall computingexpense. To pursue the above objective, then, authorsshould present as completely as possible the details oftheir formulations and the associatedcomputing-time requirements.The thermal model described here simulateswet or dry, forward or reverse combustion in one, two, or three dimensions. The formulation allowsany number and identities of components and anynumber of chemical reactions, with reactants, products, and stoichiometry specified through input products, and stoichiometry specified through input data. SPEJ P. 533


Sign in / Sign up

Export Citation Format

Share Document