DOCSLIB.ORG

Couple Openfoam with Any Other Solver Using Precice

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Couple Openfoam with Any Other Solver Using Precice Couple OpenFOAM with any other solver using preCICE Gerasimos Chourdakis et al. Technical University of Munich Department of Informatics Chair of Scientific Computing in Computer Science 2nd German OpenFoam User meetiNg TU Braunschweig February 21, 2018 Part II: + Agenda Part I: + Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN22 + Agenda Part I: + Part II: Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN22 How to simulate this heat exchanger? Geometry of a shell-and-tube heat exchanger (Image by L. Cheung Yau, 2016) Surface plot and streamlines of the two fluids colored by temperature. Solid not shown. Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN23 A shell-and-tube heat exchanger with preCICE Coupling interfaces (Image by L. Cheung Yau, 2016) Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN24 preCICE precise Code Interaction Coupling Environment libprecice Communication adapter Data Mapping Coupling Schemes Solver A Solver B • Free (GNU LGPL), developed at TU Munich & Univ. of Stuttgart. • Version 1.0 in November 2017 (10+ years, 3 PhD generations). • Official adapters for CalculiX, Code_Aster, COMSOL, Fluent, OpenFOAM, SU2 • Third-party adapters for Ateles, Alya, Carat++, FASTEST, FEAP, foam-extend, ... • API in C, C++, Fortran, Python Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN25 But why preCICE? libprecice Ope nFOAM CHT adapter CalculiX in-house FSI solver • Pure library approach ! flexibility • Fully parallel, peer-to-peer concept ! scalable and efficient communication • Sophisticated and robust quasi-Newton coupling algorithms • Multi-coupling Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN26 The roles of an adapter Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN27 + Part IIa: previous approach Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN28 Duplicated development effort Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN29 Duplicated development effort Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN29 Duplicated development effort All these adapters are bound to specific solvers! Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN29 Duplicated development effort All these adapters are bound to specific solvers! ! We need an official, general adapter! Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN29 Example of an adapted solver (previous) 1 /* Adapter: Initialize coupling */ 22 /* solve the equations */ 2 adapter.initialize(); 23 #include "rhoEqn.H" 3 24 while (pimple.loop()) 4 Info<<"\nStarting time loop\n" << endl; 25 { 5 26 ... 6 while (adapter.isCouplingOngoing()) { 27 } 7 #include "readTimeControls.H" 28 8 #include "compressibleCourantNo.H" 29 /* Adapter: Write in buffers */ 9 #include "setDeltaT.H" 30 adapter.writeCouplingData(); 10 31 11 /* Adapter: Adjust solver time */ 32 /* Adapter: advance the coupling */ 12 adapter.adjustSolverTimeStep(); 33 adapter.advance(); 13 34 14 /* Adapter: Write checkpoint */ 35 /* Adapter: Read checkpoint */ 15 if(adapter.isWriteCheckptRequired()) 36 if(adapter.isReadCheckptRequired()) 16 adapter.writeCheckpoint(); 37 adapter.readCheckpoint(); 17 38 18 runTime++; 39 if(adapter.isCouplTimeStepComplete()) 19 40 runTime.write(); 20 /* Adapter: Receive coupling data */ 41 21 adapter.readCouplingData(); 42 } Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 10 Before: Working and validated prototypes Image from desertcart.ae. Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 11 Before: Working and validated prototypes Image from desertcart.ae. Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 11 Now: A user-friendly, plug-and-play adapter The human-like figure is a property of ikea.com. Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 12 + Part IIb: a new, official adapter Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 13 Several advantages: • No source code changes • Load at runtime • (mostly) Solver agnostic However: • Still ready-to-run only for CHT • but... Making this a function object Several challenges: • No changes in the source allowed − Cannot use variables directly − Ask the objects’ registry • One adapter for all the solvers and problem types − Some parameters are not available • Only one call to execute() at the end − We may need to reload a checkpoint at the last timestep... − Set the endTime to GREAT and exit when ready. • Collaboration with other function objects − At the end, call any other end() methods explicitly. • Error handling − read() degrades errors to warnings − Catch them and throw them in execute • One adapter for all the OpenFOAM flavors and versions? − E.g. boundaryField() and boundaryFieldRef() − E.g. missing adjustTimeStep() − How to distribute? Branches/Tags? Preprocessor ifdef? • ... Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 14 However: • Still ready-to-run only for CHT • but... Making this a function object Several challenges: Several advantages: • No changes in the source allowed • No source code changes − Cannot use variables directly • Load at runtime − Ask the objects’ registry • (mostly) Solver agnostic • One adapter for all the solvers and problem types − Some parameters are not available • Only one call to execute() at the end − We may need to reload a checkpoint at the last timestep... − Set the endTime to GREAT and exit when ready. • Collaboration with other function objects − At the end, call any other end() methods explicitly. • Error handling − read() degrades errors to warnings − Catch them and throw them in execute • One adapter for all the OpenFOAM flavors and versions? − E.g. boundaryField() and boundaryFieldRef() − E.g. missing adjustTimeStep() − How to distribute? Branches/Tags? Preprocessor ifdef? • ... Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 14 Making this a function object Several challenges: Several advantages: • No changes in the source allowed • No source code changes − Cannot use variables directly • Load at runtime − Ask the objects’ registry • (mostly) Solver agnostic • One adapter for all the solvers and problem types However: − Some parameters are not available • Still ready-to-run only for CHT • Only one call to execute() at the end • but... − We may need to reload a checkpoint at the last timestep... − Set the endTime to GREAT and exit when ready. • Collaboration with other function objects − At the end, call any other end() methods explicitly. • Error handling − read() degrades errors to warnings − Catch them and throw them in execute • One adapter for all the OpenFOAM flavors and versions? − E.g. boundaryField() and boundaryFieldRef() − E.g. missing adjustTimeStep() − How to distribute? Branches/Tags? Preprocessor ifdef? • ... Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 14 An extensible adapter Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 15 OK! I want to use it! Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 16 Properties for incompressible solvers: 1 // constant/transportProperties 2 rho rho [1-300000]1; 3 Cp Cp [02-2-1000] 5000; Properties for basic solvers: 1 // constant/transportProperties 2 k k [11-3-1000] 100; OpenFOAM configuration 1 // system/controlDict 2 functions 3 { 4 preCICE_Adapter 5 { 6 type preciceAdapterFunctionObject; 7 libs ("libpreciceAdapterFunctionObject.so"); 8 } 9 } Set the appropriate boundary condition types: 1 // 0/T 2 interface 3 { 4 type fixedValue; 5 value uniform 300; 6 } 7 8 // other types: fixedGradient, mixed Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 17 OpenFOAM configuration 1 // system/controlDict 2 functions 3 { 4 preCICE_Adapter 5 { 6 type preciceAdapterFunctionObject; 7 libs ("libpreciceAdapterFunctionObject.so"); 8 } 9 } Set the appropriate boundary condition types: Properties for incompressible solvers: 1 // 0/T 1 // constant/transportProperties 2 interface 2 rho rho [1-300000]1; 3 { 3 Cp Cp [02-2-1000] 5000; 4 type fixedValue; Properties for basic solvers: 5 value uniform 300; 6 } 1 // constant/transportProperties 7 2 k k [11-3-1000] 100; 8 // other types: fixedGradient, mixed Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 17 preCICE & adapter configuration libprecice adapter Solver A Solver B YAML YAML XML Adapter A Config Adapter B Config preCICE Config To run the simulation, just execute the solvers as usual. Gerasimos Chourdakis (TUM) | Couple OpenFOAM with any other solver using preCICE | Feb 21, 2018 | GOFUN2 18 Tutorials Gerasimos
Load more

Recommended publications
  • Advanced Molecular Dynamics in Openfoam✩ S.M
    Computer Physics Communications 224 (2018) 1–21 Contents lists available at ScienceDirect Computer Physics Communications journal homepage: www.elsevier.com/locate/cpc Feature article mdFoam+: Advanced molecular dynamics in OpenFOAMI S.M. Longshaw a,*, M.K. Borg b, S.B. Ramisetti b, J. Zhang b, D.A. Lockerby c, D.R. Emerson a, J.M. Reese b a Scientific Computing Department, The Science & Technology Facilities Council, Daresbury Laboratory, Warrington, Cheshire WA4 4AD, UK b School of Engineering, University of Edinburgh, Edinburgh EH9 3FB, UK c School of Engineering, University of Warwick, Coventry, CV4 7AL, UK article info a b s t r a c t Article history: This paper introduces mdFoam+, which is an MPI parallelised molecular dynamics (MD) solver imple- Received 6 March 2017 mented entirely within the OpenFOAM software framework. It is open-source and released under the Received in revised form 17 August 2017 same GNU General Public License (GPL) as OpenFOAM. The source code is released as a publicly open Accepted 25 September 2017 software repository that includes detailed documentation and tutorial cases. Since mdFoam+ is designed Available online 23 October 2017 entirely within the OpenFOAM C++ object-oriented framework, it inherits a number of key features. The Keywords: code is designed for extensibility and flexibility, so it is aimed first and foremost as an MD research tool, OpenFOAM in which new models and test cases can be developed and tested rapidly. Implementing mdFoam+ in Molecular dynamics OpenFOAM also enables easier development of hybrid methods that couple MD with continuum-based MD solvers. Setting up MD cases follows the standard OpenFOAM format, as mdFoam+ also relies upon Computational fluid dynamics the OpenFOAM dictionary-based directory structure.
    [Show full text]
  • Introduction to Nuclear Reactor Modelling Using Openfoam Carlo Fiorina 2
    Introduction to nuclear reactor modelling using OpenFOAM Carlo Fiorina 2 Disclaimer C. Fiorina ▪ Focus on the use of OpenFOAM for multiphysics ○ Use of OpenFOAM as CFD tool widely covered by documentation, forums, courses, etc. ▪ Focus on already existing tool (GeN-Foam) as an example ○ Programming from scratch is not that difficult, but unsuited for a 75 minutes lecture ▪ In the slides, more material than can actually be covered in this lecture ○ Can help better understanding the slides after the lecture 3 Content C. Fiorina ▪ General Introduction ▪ Introduction to the use of OpenFOAM ▪ Basics of GeN-Foam ▪ Short introduction on the use of GeN-Foam 4 Objective C. Fiorina What is it about? ▪ Provide you with general information, references, suggestions, terminology and lessons learnt that can facilitate your approach to the OpenFOAM world ▪ Provide with slides that can help you out orienting yourself if you decide to embrace the use of OpenFOAM What is not about? ▪ Detailed course on the use of OpenFOAM ▪ Hands-on training 5 The IAEA-facilitated ONCORE initiative C. Fiorina ONCORE to support the open-source nuclear community and help addressing typical shortcomings of open-source development (scattered community, documentation, QA, loss of knowledge) ▪ Promote collaboration and facilitate communication (connect the community) ▪ Provide guidelines for code contribution (documentation, QA) ▪ Provide development best practices (QA) ▪ Preserve knowledge ○ Incl. compiling a list of open-source codes https://www.iaea.org/topics/nuclear-power-reactors/open-source-nuclear-cod e-for-reactor-analysis-oncore 6 A first important outcome: list of C. Fiorina available codes ▪ https://nucleus.iaea.org/sites/oncore/SitePages/List%20of%20Codes.aspx ▪ A vibrant community with an impressive R&D output ▪ ~35 codes already identified so far: ○ OpenMC ○ Raven ○ Dragon ○ MOOSE ○ Salome platform (Code_Saturne, Code_Aster) ○ TrioCFD ○ … ○ Several OpenFOAM-based tools 7 A central tool for open-source simulation C.
    [Show full text]
  • 15Th Openfoam Workshop - June 22, 2020 - Technical Session I-D a Flexible and Precice Solver Coupling Ecosystem
    15th OpenFOAM Workshop - June 22, 2020 - Technical Session I-D A flexible and preCICE solver coupling ecosystem Gerasimos Chourdakis, Technical University of Munich [email protected] Benjamin Uekermann, Eindhoven University of Technology [email protected] Find these slides on GitHub https://github.com/MakisH/ofw15-slides The Big Picture e r ic r te c e p e lv a r o d p s a lib structure fluid solver solver The Big Picture e r ic r te c e p e lv a r o d p s a lib structure fluid solver solver in-house commercial solver solver The Big Picture e r ic r te c e p e lv a r o d p s a lib structure fluid solver solver OpenFOAM CalculiX SU2 Code_Aster foam-extend FEniCS deal-ii Nutils MBDyn in-house commercial solver solver API in: C++ Python Matlab ANSYS Fluent C Fortran COMSOL The Big Picture e r ic r te c e p e lv a r o d p s a lib structure fluid solver A Coupling Library for Partitioned Multi-Physics Simulations solver OpenFOAM CalculiX SU2 Code_Aster . foam-extend FEniCS deal-ii . Nutils MBDyn communication data mapping in-house commercial solver solver coupling schemes time interpolation API in: C++ Python Matlab ANSYS Fluent C Fortran COMSOL News preCICE v2.0 Simplified config & API Better building & testing XML reference & visualizer xSDK member Faster initialization Better Python bindings Spack / Debian / AUR New Matlab bindings packages Upgrade guide in the wiki Other news deal.ii adapter new non-linear example for FSI FEniCS adapter new example for FSI code_aster adapter revived for code_aster 14 and preCICE v2 The OpenFOAM adapter
    [Show full text]
  • Model Order Reduction for Aerodynamic Lifting Surfaces Aerospace Engineering
    Model Order Reduction for Aerodynamic Lifting Surfaces Gonçalo da Cunha Laboreiro Mendonça Thesis to obtain the Master of Science Degree in Aerospace Engineering Supervisors: Prof. Fernando José Parracho Lau Dr. Frederico José Prata Rente Reis Afonso Examination Committee Chairperson: Prof. Filipe Szolnoky Ramos Pinto Cunha Supervisor: Prof. Fernando José Parracho Lau Member of the Committee: Prof. Afzal Suleman November 2017 ii Dedicated to my family and friends, who were always there for me. iii iv Acknowledgments I would like to thank dearly my supervisors Prof. Lau and Dr. Afonso who gave me constant support throughout this thesis up until the very end. Their insights on aerodynamics and CFD models guided me in my research of model order reduction methods and allowed me to better understand the models with which I had to work. Their demands for rigor and quality also pushed me to better my work, and in the end write a better thesis. v vi Resumo Nesta tese o tema de reduc¸ao˜ de modelos e a sua aplicac¸ao˜ a Mecanicaˆ de Fluidos Computacional sao˜ abordados. E´ mostrada a necessidade da industria´ aeroespacial, seja nacional ou Europeia, de modelos mais rapidos´ mas fieis´ a` realidade. Isto e´ devido ao elevado tempo de calculo´ associado aos modelos de alta-fidelidade. Estes mostram-se pouco viaveis´ para aplicac¸oes˜ do tipo Optimizac¸ao˜ Multi- disciplinar, como a plataforma de optimizac¸ao˜ NOVEMOR. Tendo por objectivo testar e aplicar reduc¸ao˜ de modelos a modelos CFD de superf´ıcies sustentadoras, uma pesquisa bibliografica´ abrangendo a reduc¸ao˜ de modelos nao-lineares,˜ dinamicosˆ e ou estaticos´ foi feita.
    [Show full text]
  • The One-Way FSI Method Based on RANS-FEM for the Open Water Test of a Marine Propeller at the Different Loading Conditions
    Journal of Marine Science and Engineering Article The One-Way FSI Method Based on RANS-FEM for the Open Water Test of a Marine Propeller at the Different Loading Conditions Mobin Masoomi 1 and Amir Mosavi 2,3,* 1 Department of Mechanical Engineering, Babol Noshirvani University of Technology, Babol, Iran; [email protected] 2 Faculty of Civil Engineering, Technische Universität Dresden, 01069 Dresden, Germany 3 John von Neumann Faculty of Informatics, Obuda University, 1034 Budapest, Hungary * Correspondence: [email protected] Abstract: This paper aims to assess a new fluid–structure interaction (FSI) coupling approach for the vp1304 propeller to predict pressure and stress distributions with a low-cost and high-precision approach with the ability of repeatability for the number of different structural sets involved, other materials, or layup methods. An outline of the present coupling approach is based on an open- access software (OpenFOAM) as a fluid solver, and Abaqus used to evaluate and predict the blade’s deformation and strength in dry condition mode, which means the added mass effects due to propeller blades vibration is neglected. Wherein the imposed pressures on the blade surfaces are extracted for all time-steps. Then, these pressures are transferred to the structural solver as a load condition. Although this coupling approach was verified formerly (wedge impact), for the case in-hand, a further verification case, open water test, was performed to evaluate the hydrodynamic part of the solution with an e = 7.5% average error. A key factor for the current coupling approach is Citation: Masoomi, M.; Mosavi, A.
    [Show full text]
  • Wind Field Simulation in a Wind Farm Using Openfoam and Actuator Line Model
    ParCFD'2019 31st International Conference on Parallel Computational Fluid Dynamics May-14-17 2019, Antalya TURKEY WIND FIELD SIMULATION IN A WIND FARM USING OPENFOAM AND ACTUATOR LINE MODEL Huseyin Can Onel∗ & Dr. Ismail H. Tuncery ∗ Middle East Technical University (METU) Department of Aerospace Engineering 06800 Ankara, TURKEY e-mail: [email protected] yMiddle East Technical University (METU) Department of Aerospace Engineering 06800 Ankara, TURKEY e-mail: [email protected] - Web page: http://www.ae.metu.edu.tr/tuncer/ Key words: Aerospace applications, Wind turbine, HAWT, Actuator Line Model, Wake calculation Abstract. In this study, a horizontal axis wind turbine (HAWT) is modeled using so called Actuator Line Model (ALM), where full resolution of boundary layer over turbine blades is not needed and hence computation is cheaper. Results are validated against other numerical and experimental studies as well as panel method (XFOIL) and Blade Element Momentum Theory (BEMT) results which are still widely employed in today's wind energy industry. Important simulation and operation parameters and their effects on accuracy are discussed. It is concluded that within a certain range of tip speed ratios, ALM gives acceptable results and is a promising model for full-scale wind farm simulations to estimate energy production. 1 INTRODUCTION Market share of renewable energy grows at ever highest rates and wind turbine and wind farm design processes becomes more sophisticated with the advancements in computation technologies. There are two main design problems in wind energy: • Design of an individual wind turbine at its ideal operation conditions, where classical methods like 2D airfoil theory, potential flow theory and Blade Element Momentum Theory (BEMT) are still widely used, • Design of a complete wind farm, in which statistical meteorological data is used for macro-siting and simple analytical or empirical methods are used for micro-siting.
    [Show full text]
  • EOF-Library: Open-Source Elmer FEM and Openfoam Coupler For
    EOF-Library: Open-source Elmer FEM and OpenFOAM coupler for electromagnetics and fluid dynamics (CC BY-SA 4.0) Motivation for MHD - controlling liquid metals Liquid metal pump rotating permanent pipe magnets F1 F2 F FL FL Mechanically FL FL Electromagnetically Required physics Modelling coupled ● Electromagnetics ● Fluid dynamics Extra ● Free surface ● Heat transfer Liquid metal pump - modelling requires two-way ● Solidification / melting coupled electromagnetics and fluid dynamics 1 Required numerical models Fluid dynamics ● 3D, transient Electromagnetics ● High Reynolds ● Complex geometries, turbulence models multi-regions ● Complex A−V ● Volume of Fluid magnetic vector ● Advanced pre- and potential (VOF) for free post-processing tools surface modelling formulation ● Parallelization with ● Conductivity ● Viscosity good scaling dependence on dependence on temperature ● User community, temperature documentation 2 Existing Open Source codes Best for Best for Electromagnetics (EM) Fluid dynamics (FD) FEM FVM getDP OpenFOAM Elmer SU2 ... ... 3 Choosing Electromagnetics solver (1) GetDP + Designed for EM in mind, great for eddy-current problems + Uses PETSc iterative solvers, PETSc is MPI parallelized - Matrix assembly is not parallel (neither MPI nor OpenMP) - Weak integration with external post-processing tools MaxFEM, FreeFem++, deal.II, … Pictures from getdp.info ● I have no opinion 4 Choosing Electromagnetics solver (2) Elmer by CSC (Finland) + Multiphysics out-of-the-box (>40 solvers) + Build-in linear & nonlinear iterative solvers,
    [Show full text]
  • Openfoam and Third Party Structural Solver for Fluid–Structure Interaction Simulations Robert L
    OpenFOAM and Third Party Structural Solver for Fluid–Structure Interaction Simulations Robert L. Campbell Brent A. Craven [email protected] [email protected] Fluids and Structural Mechanics Office Applied Research Laboratory The Pennsylvania State University 6th OpenFOAM Workshop Penn State University University Park, Pa 13-17 June 2011 Objective ● Demonstrate a method to include a third-party, structural solver of choice for FSI simulations (quasi-steady) ● Why important? Structural analysts seem to have strong ties to their structural software and want to carry them forward for FSI simulations 2 Outline ● Brief Introduction to Fluid–Structure Interaction – Monolithic – Partitioned • Loose coupling • Tight coupling ● Solver Implementation – Flow solver – Structure solver – Solver interface ● Example Problem – Setup – Solution 3 FSI Introduction ● Fluid–Structure Interaction (FSI) modeling is a type of multi- physics modeling that combines physical models of fluids and structures ● FSI, in this context, implies two-way coupling ● Two approaches to FSI modeling – Monolithic • Governing equations for both the fluid and solid cast in terms of the same primitive variables • Single discretization scheme applied to entire domain – Partitioned • Fluid and solid domains modeled separately • Separate discretizations of each domain • Stress and displacement communication across the domain interface • Two coupling schemes: Loose and Tight 4 FSI Loose Coupling ● Loose coupling approach does START not check for convergence at each time step: Estimate
    [Show full text]
  • Simulation in CFD of a Pebble Bed: Advanced High Temperature Reactor
    2017 International Nuclear Atlantic Conference - INAC 2017 Belo Horizonte, MG, Brazil, October 22-27, 2017 ASSOCIAC¸ AO˜ BRASILEIRA DE ENERGIA NUCLEAR - ABEN SIMULATION IN CFD OF A PEBBLE BED - ADVANCED HIGH TEMPERATURE REACTOR CORE USING OPENFOAM Pamela M. Dahl1 and Jian Su1 1Programa de Engenharia Nuclear, COPPE Universidade Federal do Rio de Janeiro 21941-972 Rio de Janeiro, RJ [email protected] ABSTRACT Numerical simulations of a Pebble Bed nuclear reactor core are presented using the multi-physics tool-kit OpenFOAM. The HTR-PM is modeled using the porous media approach, accounting both for viscous and inertial effects through the Darcy and Forchheimer model. Initially, cylindrical 2D and 3D simulations are compared, in order to evaluate their differences and decide if the 2D simulations carry enough of the sought information, considering the savings in computational costs. The porous medium is considered to be isotropic, with the whole length of the packed bed occupied homogeneously with the spherical fuel elements. Steady-state simulations for normal equilibrium operation are performed, using a semi- sine function of the power density along the vertical axis as the source term for the energy balance equation.Total pressure drop is calculated and compared with that obtained from literature for a similar case. At a second stage, transient simulations are performed, where relevant parameters are calculated and compared to those of the literature. 1. INTRODUCTION The modular high temperature gas-cooled reactor (HTR) is quickly becoming a strong candidate for the IVth Generation of nuclear energy system technology, because of its inherent safety features, modularity, relatively low cost and a shorter construction time, which are obvious attractive features for a commercial reactor [1].
    [Show full text]
  • Open Source, Component Based Simulation Software Development Using Orcan
    Open Source, component based simulation software development using Orcan Horst Hadler University Erlangen-Nuernberg, Department of Computer Science 9 Thomas Jung Fraunhofer Institute of Integrated Systems and Device Technology, Erlangen (IISB) Michael Kellner University Erlangen-Nuernberg, Department of Materials Science 6 Department of Computer Science 10 Crystal Growth Lab (Prof. G. Mueller): www.cgl-erlangen.com Institute for Material Sciences VI, University Erlangen-Nuremberg Fraunhofer-Institute IISB, Erlangen Development and analysis of crystal growth processes (Si, InP, GaAs, CaF2, oxide crystals, ...) Experiments + Simulation CrysVUn: TMT for the user support program of the MSL Background of ORCAN development: need to be able to simulate radiative heat transfer in complex 3D geometries ± including volume effects (absorption, scattering), and different surface effects (diffuse/ specular reflection, ...), e.g. for optical crystals, szintillators, ... Options: Use a commercial CFD-code ? +? most is already done... (however, radiation models are weak) - you don©t know what is going on behind the scenes - difficult to extend - expensive, take a lot of learning time: have to select one, and stick to it Develop your own proprietary code ? - beyond our possibilities: we cannot do everything needed Use the wealth of existing open source tools: mesh generators, solvers, geometry handling, visualizations, ... + its free .. + source code available -> possibility to check what is really done, and to extend the code - no common interfaces:
    [Show full text]
  • Agenda Item Submittal
    AGENDA ADMINISTRATION/FINANCE ISSUES COMMITTEE MEETING WITH BOARD OF DIRECTORS* ORANGE COUNTY WATER DISTRICT 18700 Ward Street, Fountain Valley, CA 92708 Thursday, March 15, 2018, 8:00 a.m. - Conference Room C-2 *The OCWD Administration and Finance Issues Committee meeting is noticed as a joint meeting with the Board of Directors for the purpose of strict compliance with the Brown Act and it provides an opportunity for all Directors to hear presentations and participate in discussions. Directors receive no additional compensation or stipend as a result of simultaneously convening this meeting. Items recommended for approval at this meeting will be placed on the March 21, 2018 Board meeting Agenda for approval. ROLL CALL ITEMS RECEIVED TOO LATE TO BE AGENDIZED RECOMMENDATION: Adopt resolution determining need to take immediate action on item(s) and that the need for action came to the attention of the District subsequent to the posting of the Agenda (requires two-thirds vote of the Board members present, or, if less than two-thirds of the members are present, a unanimous vote of those members present). VISITOR PARTICIPATION Time has been reserved at this point in the agenda for persons wishing to comment for up to three minutes to the Board of Directors on any item that is not listed on the agenda, but within the subject matter jurisdiction of the District. By law, the Board of Directors is prohibited from taking action on such public comments. As appropriate, matters raised in these public comments will be referred to District staff or placed on the agenda of an upcoming Board meeting.
    [Show full text]
  • Vertical-Axis Wind-Turbine Computations Using a 2D Hybrid Wake Actuator-Cylinder Model
    Wind Energ. Sci., 6, 1061–1077, 2021 https://doi.org/10.5194/wes-6-1061-2021 © Author(s) 2021. This work is distributed under the Creative Commons Attribution 4.0 License. Vertical-axis wind-turbine computations using a 2D hybrid wake actuator-cylinder model Edgar Martinez-Ojeda1, Francisco Javier Solorio Ordaz1, and Mihir Sen2 1Facultad de Ingeniería, Universidad Nacional Autónoma de México, Mexico City, 04510, Mexico 2Department of Aerospace and Mechanical Engineering, University of Notre Dame, Notre Dame, IN 46556, USA Correspondence: Edgar Martinez-Ojeda ([email protected]) Received: 7 May 2021 – Discussion started: 10 May 2021 Revised: 8 July 2021 – Accepted: 19 July 2021 – Published: 16 August 2021 Abstract. The actuator-cylinder model was implemented in OpenFOAM by virtue of source terms in the Navier–Stokes equations. Since the stand-alone actuator cylinder is not able to properly model the wake of a vertical-axis wind turbine, the steady incompressible flow solver simpleFoam provided by OpenFOAM was used to resolve the entire flow and wakes of the turbines. The source terms are only applied inside a certain region of the computational domain, namely a finite-thickness cylinder which represents the flight path of the blades. One of the major advantages of this approach is its implicitness – that is, the velocities inside the hollow cylinder region feed the stand-alone actuator-cylinder model (AC); this in turn computes the volumetric forces and passes them to the OpenFOAM solver in order to be applied inside the hollow cylinder region. The process is repeated in each iteration of the solver until convergence is achieved.
    [Show full text]