Adaptive Finite Volume numerical method
This work describes a Finite Volume computational method for the parametric study of phenomena in plasmas, i.e. in situations where numerous runs of three-dimensional simulations are required (in many cases, thousands or tens of thousands). Such problems require simple codes, robust, modular (to add...
Autores principales: | , , , |
---|---|
Formato: | CONF |
Materias: | |
Acceso en línea: | http://hdl.handle.net/20.500.12110/paper_17426588_v591_n1_p_Bilbao |
Aporte de: |
id |
todo:paper_17426588_v591_n1_p_Bilbao |
---|---|
record_format |
dspace |
spelling |
todo:paper_17426588_v591_n1_p_Bilbao2023-10-03T16:31:00Z Adaptive Finite Volume numerical method Bilbao, L. Herrera-Velazquez J.J.E. Vargas-Blanco I. Ad Astra Rocket Company; Instituto Tecnologico de Costa Rica; International Atomic Energy Agency (IAEA); Universidad Nacional de Costa Rica Computational efficiency Fluid dynamics Fusion reactors Hydrodynamics Mesh generation Supersonic aircraft Finite-amplitude perturbations Irregular meshes Kelvin- helmholtz instabilities Physical process Spatial resolution Spatial response Temporal response Three dimensional simulations Numerical methods This work describes a Finite Volume computational method for the parametric study of phenomena in plasmas, i.e. in situations where numerous runs of three-dimensional simulations are required (in many cases, thousands or tens of thousands). Such problems require simple codes, robust, modular (to add or remove physical processes) and do not require high precision. The code is based on a complex multi-component species program with transport and radiation terms. The integration domain is represented with a structured irregular mesh, with fixed connectivity. A new algorithm for the hydrodynamics was implemented in order to improve the computational efficiency plus the improved capability of adapting the mesh to the solution. The improved hydrodynamics method worked well in an ample range of Mach number from subsonic (10-3) to supersonic. After each calculation cycle, mesh vertices are moved arbitrary over the fluid. This is done in order to dynamically adapt the mesh to the solution. The adaptive method consists of shifting mesh vertices over the fluid in order to keep a reasonable mesh structure and increase the spatial resolution where the physical solution demands. As an example, we show the results of the development of the Kelvin-Helmholtz instability in local plane slab models of the magnetopause, showing the development and saturation of the instability in an initially unperturbed structure, i.e., the temporal response approach; and the response of a background equilibrium to the excitation by finite amplitude perturbations generated upstream, i.e., the spatial response of the system. © Published under licence by IOP Publishing Ltd. Fil:Bilbao, L. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales; Argentina. CONF info:eu-repo/semantics/openAccess http://creativecommons.org/licenses/by/2.5/ar http://hdl.handle.net/20.500.12110/paper_17426588_v591_n1_p_Bilbao |
institution |
Universidad de Buenos Aires |
institution_str |
I-28 |
repository_str |
R-134 |
collection |
Biblioteca Digital - Facultad de Ciencias Exactas y Naturales (UBA) |
topic |
Computational efficiency Fluid dynamics Fusion reactors Hydrodynamics Mesh generation Supersonic aircraft Finite-amplitude perturbations Irregular meshes Kelvin- helmholtz instabilities Physical process Spatial resolution Spatial response Temporal response Three dimensional simulations Numerical methods |
spellingShingle |
Computational efficiency Fluid dynamics Fusion reactors Hydrodynamics Mesh generation Supersonic aircraft Finite-amplitude perturbations Irregular meshes Kelvin- helmholtz instabilities Physical process Spatial resolution Spatial response Temporal response Three dimensional simulations Numerical methods Bilbao, L. Herrera-Velazquez J.J.E. Vargas-Blanco I. Ad Astra Rocket Company; Instituto Tecnologico de Costa Rica; International Atomic Energy Agency (IAEA); Universidad Nacional de Costa Rica Adaptive Finite Volume numerical method |
topic_facet |
Computational efficiency Fluid dynamics Fusion reactors Hydrodynamics Mesh generation Supersonic aircraft Finite-amplitude perturbations Irregular meshes Kelvin- helmholtz instabilities Physical process Spatial resolution Spatial response Temporal response Three dimensional simulations Numerical methods |
description |
This work describes a Finite Volume computational method for the parametric study of phenomena in plasmas, i.e. in situations where numerous runs of three-dimensional simulations are required (in many cases, thousands or tens of thousands). Such problems require simple codes, robust, modular (to add or remove physical processes) and do not require high precision. The code is based on a complex multi-component species program with transport and radiation terms. The integration domain is represented with a structured irregular mesh, with fixed connectivity. A new algorithm for the hydrodynamics was implemented in order to improve the computational efficiency plus the improved capability of adapting the mesh to the solution. The improved hydrodynamics method worked well in an ample range of Mach number from subsonic (10-3) to supersonic. After each calculation cycle, mesh vertices are moved arbitrary over the fluid. This is done in order to dynamically adapt the mesh to the solution. The adaptive method consists of shifting mesh vertices over the fluid in order to keep a reasonable mesh structure and increase the spatial resolution where the physical solution demands. As an example, we show the results of the development of the Kelvin-Helmholtz instability in local plane slab models of the magnetopause, showing the development and saturation of the instability in an initially unperturbed structure, i.e., the temporal response approach; and the response of a background equilibrium to the excitation by finite amplitude perturbations generated upstream, i.e., the spatial response of the system. © Published under licence by IOP Publishing Ltd. |
format |
CONF |
author |
Bilbao, L. Herrera-Velazquez J.J.E. Vargas-Blanco I. Ad Astra Rocket Company; Instituto Tecnologico de Costa Rica; International Atomic Energy Agency (IAEA); Universidad Nacional de Costa Rica |
author_facet |
Bilbao, L. Herrera-Velazquez J.J.E. Vargas-Blanco I. Ad Astra Rocket Company; Instituto Tecnologico de Costa Rica; International Atomic Energy Agency (IAEA); Universidad Nacional de Costa Rica |
author_sort |
Bilbao, L. |
title |
Adaptive Finite Volume numerical method |
title_short |
Adaptive Finite Volume numerical method |
title_full |
Adaptive Finite Volume numerical method |
title_fullStr |
Adaptive Finite Volume numerical method |
title_full_unstemmed |
Adaptive Finite Volume numerical method |
title_sort |
adaptive finite volume numerical method |
url |
http://hdl.handle.net/20.500.12110/paper_17426588_v591_n1_p_Bilbao |
work_keys_str_mv |
AT bilbaol adaptivefinitevolumenumericalmethod AT herreravelazquezjje adaptivefinitevolumenumericalmethod AT vargasblancoi adaptivefinitevolumenumericalmethod AT adastrarocketcompanyinstitutotecnologicodecostaricainternationalatomicenergyagencyiaeauniversidadnacionaldecostarica adaptivefinitevolumenumericalmethod |
_version_ |
1807320214762684416 |