Articles de revistahttp://hdl.handle.net/2117/805232021-07-25T23:25:43Z2021-07-25T23:25:43ZOverview of physics studies on ASDEX UpgradeMeyer, H.Aguiam, D.Angioni, C.Albert, C.G.Arden, N.Mantsinen, MerviSchmitz, O.http://hdl.handle.net/2117/3497132021-07-25T20:40:26Z2021-07-20T08:20:24ZOverview of physics studies on ASDEX Upgrade
Meyer, H.; Aguiam, D.; Angioni, C.; Albert, C.G.; Arden, N.; Mantsinen, Mervi; Schmitz, O.
The ASDEX Upgrade (AUG) programme, jointly run with the EUROfusion MST1 task force, continues to significantly enhance the physics base of ITER and DEMO. Here, the full tungsten wall is a key asset for extrapolating to future devices. The high overall heating power, flexible heating mix and comprehensive diagnostic set allows studies ranging from mimicking the scrape-off-layer and divertor conditions of ITER and DEMO at high density to fully non-inductive operation (q95 = 5.5, ) at low density. Higher installed electron cyclotron resonance heating power 6 MW, new diagnostics and improved analysis techniques have further enhanced the capabilities of AUG.
Stable high-density H-modes with MW m−1 with fully detached strike-points have been demonstrated. The ballooning instability close to the separatrix has been identified as a potential cause leading to the H-mode density limit and is also found to play an important role for the access to small edge-localized modes (ELMs). Density limit disruptions have been successfully avoided using a path-oriented approach to disruption handling and progress has been made in understanding the dissipation and avoidance of runaway electron beams. ELM suppression with resonant magnetic perturbations is now routinely achieved reaching transiently . This gives new insight into the field penetration physics, in particular with respect to plasma flows. Modelling agrees well with plasma response measurements and a helically localised ballooning structure observed prior to the ELM is evidence for the changed edge stability due to the magnetic perturbations. The impact of 3D perturbations on heat load patterns and fast-ion losses have been further elaborated.
Progress has also been made in understanding the ELM cycle itself. Here, new fast measurements of and Er allow for inter ELM transport analysis confirming that Er is dominated by the diamagnetic term even for fast timescales. New analysis techniques allow detailed comparison of the ELM crash and are in good agreement with nonlinear MHD modelling. The observation of accelerated ions during the ELM crash can be seen as evidence for the reconnection during the ELM. As type-I ELMs (even mitigated) are likely not a viable operational regime in DEMO studies of 'natural' no ELM regimes have been extended. Stable I-modes up to have been characterised using -feedback.
Core physics has been advanced by more detailed characterisation of the turbulence with new measurements such as the eddy tilt angle—measured for the first time—or the cross-phase angle of and fluctuations. These new data put strong constraints on gyro-kinetic turbulence modelling. In addition, carefully executed studies in different main species (H, D and He) and with different heating mixes highlight the importance of the collisional energy exchange for interpreting energy confinement. A new regime with a hollow profile now gives access to regimes mimicking aspects of burning plasma conditions and lead to nonlinear interactions of energetic particle modes despite the sub-Alfvénic beam energy. This will help to validate the fast-ion codes for predicting ITER and DEMO.
2021-07-20T08:20:24ZMeyer, H.Aguiam, D.Angioni, C.Albert, C.G.Arden, N.Mantsinen, MerviSchmitz, O.The ASDEX Upgrade (AUG) programme, jointly run with the EUROfusion MST1 task force, continues to significantly enhance the physics base of ITER and DEMO. Here, the full tungsten wall is a key asset for extrapolating to future devices. The high overall heating power, flexible heating mix and comprehensive diagnostic set allows studies ranging from mimicking the scrape-off-layer and divertor conditions of ITER and DEMO at high density to fully non-inductive operation (q95 = 5.5, ) at low density. Higher installed electron cyclotron resonance heating power 6 MW, new diagnostics and improved analysis techniques have further enhanced the capabilities of AUG.
Stable high-density H-modes with MW m−1 with fully detached strike-points have been demonstrated. The ballooning instability close to the separatrix has been identified as a potential cause leading to the H-mode density limit and is also found to play an important role for the access to small edge-localized modes (ELMs). Density limit disruptions have been successfully avoided using a path-oriented approach to disruption handling and progress has been made in understanding the dissipation and avoidance of runaway electron beams. ELM suppression with resonant magnetic perturbations is now routinely achieved reaching transiently . This gives new insight into the field penetration physics, in particular with respect to plasma flows. Modelling agrees well with plasma response measurements and a helically localised ballooning structure observed prior to the ELM is evidence for the changed edge stability due to the magnetic perturbations. The impact of 3D perturbations on heat load patterns and fast-ion losses have been further elaborated.
Progress has also been made in understanding the ELM cycle itself. Here, new fast measurements of and Er allow for inter ELM transport analysis confirming that Er is dominated by the diamagnetic term even for fast timescales. New analysis techniques allow detailed comparison of the ELM crash and are in good agreement with nonlinear MHD modelling. The observation of accelerated ions during the ELM crash can be seen as evidence for the reconnection during the ELM. As type-I ELMs (even mitigated) are likely not a viable operational regime in DEMO studies of 'natural' no ELM regimes have been extended. Stable I-modes up to have been characterised using -feedback.
Core physics has been advanced by more detailed characterisation of the turbulence with new measurements such as the eddy tilt angle—measured for the first time—or the cross-phase angle of and fluctuations. These new data put strong constraints on gyro-kinetic turbulence modelling. In addition, carefully executed studies in different main species (H, D and He) and with different heating mixes highlight the importance of the collisional energy exchange for interpreting energy confinement. A new regime with a hollow profile now gives access to regimes mimicking aspects of burning plasma conditions and lead to nonlinear interactions of energetic particle modes despite the sub-Alfvénic beam energy. This will help to validate the fast-ion codes for predicting ITER and DEMO.Wildfires vegetation recovery through satellite remote sensing and functional data analysisSerra Burriel, FeliuDelicado Useros, Pedro FranciscoCucchietti, Fernandohttp://hdl.handle.net/2117/3493672021-07-18T21:01:15Z2021-07-15T07:16:52ZWildfires vegetation recovery through satellite remote sensing and functional data analysis
Serra Burriel, Feliu; Delicado Useros, Pedro Francisco; Cucchietti, Fernando
In recent years, wildfires have caused havoc across the world, which are especially aggravated in certain regions due to climate change. Remote sensing has become a powerful tool for monitoring fires, as well as for measuring their effects on vegetation over the following years. We aim to explain the dynamics of wildfires’ effects on a vegetation index (previously estimated by causal inference through synthetic controls) from pre-wildfire available information (mainly proceeding from satellites). For this purpose, we use regression models from Functional Data Analysis, where wildfire effects are considered functional responses, depending on elapsed time after each wildfire, while pre-wildfire information acts as scalar covariates. Our main findings show that vegetation recovery after wildfires is a slow process, affected by many pre-wildfire conditions, among which the richness and diversity of vegetation is one of the best predictors for the recovery.
2021-07-15T07:16:52ZSerra Burriel, FeliuDelicado Useros, Pedro FranciscoCucchietti, FernandoIn recent years, wildfires have caused havoc across the world, which are especially aggravated in certain regions due to climate change. Remote sensing has become a powerful tool for monitoring fires, as well as for measuring their effects on vegetation over the following years. We aim to explain the dynamics of wildfires’ effects on a vegetation index (previously estimated by causal inference through synthetic controls) from pre-wildfire available information (mainly proceeding from satellites). For this purpose, we use regression models from Functional Data Analysis, where wildfire effects are considered functional responses, depending on elapsed time after each wildfire, while pre-wildfire information acts as scalar covariates. Our main findings show that vegetation recovery after wildfires is a slow process, affected by many pre-wildfire conditions, among which the richness and diversity of vegetation is one of the best predictors for the recovery.Measuring and improving the geometric accuracy of piece-wise polynomial boundary meshesRuiz Gironès, EloiSarrate Ramos, JosepRoca Navarro, Francisco Javierhttp://hdl.handle.net/2117/3490242021-07-18T08:07:05Z2021-07-12T13:16:39ZMeasuring and improving the geometric accuracy of piece-wise polynomial boundary meshes
Ruiz Gironès, Eloi; Sarrate Ramos, Josep; Roca Navarro, Francisco Javier
We present a new disparity functional to measure and improve the geometric accuracy of a curved high-order mesh that approximates a target geometry model. We have devised the disparity to account for compound models, be independent of the entity parameterization, and allow trimmed entities. The disparity depends on the physical mesh and the auxiliary parametric meshes. Since it is two times differentiable on all these variables, we can minimize it with a second-order method. Its minimization with the parametric meshes as design variables measures the geometric accuracy of a given mesh. Furthermore, the minimization with both the physical and parametric meshes as design variables improves the geometric accuracy of an initial mesh. We have numerical evidence that the obtained meshes converge to the target geometry (unitary normal) algebraically, in terms of the element size, with order 2p (2p-1, respectively), where p is the polynomial degree of the mesh. Although we obtain meshes with non-interpolative boundary nodes, we propose a post-process to enforce, if required by the application, meshes with interpolative boundary nodes and featuring the same order of geometric accuracy. In conclusion, we can obtain super-convergent orders, at least for sufficiently smooth parametric curve (surface) entities, for meshes of polynomial degrees up to 4 (3, respectively). In perspective, this super-convergence might enable using a lower polynomial degree to approximate the geometry than to approximate the solution without hampering the required geometric accuracy for high-order analysis.
2021-07-12T13:16:39ZRuiz Gironès, EloiSarrate Ramos, JosepRoca Navarro, Francisco JavierWe present a new disparity functional to measure and improve the geometric accuracy of a curved high-order mesh that approximates a target geometry model. We have devised the disparity to account for compound models, be independent of the entity parameterization, and allow trimmed entities. The disparity depends on the physical mesh and the auxiliary parametric meshes. Since it is two times differentiable on all these variables, we can minimize it with a second-order method. Its minimization with the parametric meshes as design variables measures the geometric accuracy of a given mesh. Furthermore, the minimization with both the physical and parametric meshes as design variables improves the geometric accuracy of an initial mesh. We have numerical evidence that the obtained meshes converge to the target geometry (unitary normal) algebraically, in terms of the element size, with order 2p (2p-1, respectively), where p is the polynomial degree of the mesh. Although we obtain meshes with non-interpolative boundary nodes, we propose a post-process to enforce, if required by the application, meshes with interpolative boundary nodes and featuring the same order of geometric accuracy. In conclusion, we can obtain super-convergent orders, at least for sufficiently smooth parametric curve (surface) entities, for meshes of polynomial degrees up to 4 (3, respectively). In perspective, this super-convergence might enable using a lower polynomial degree to approximate the geometry than to approximate the solution without hampering the required geometric accuracy for high-order analysis.Vortex induced vibrations of a pivoted finite height cylinder at low Reynolds numberCajas García, Juan CarlosPastrana Maldonado, DanielRodríguez Pérez, Ivette MaríaLehmkuhl Barba, OriolHouzeaux, GuillaumeVázquez, MarianoTreviño Treviño, Cesarhttp://hdl.handle.net/2117/3466652021-06-13T19:20:26Z2021-06-07T07:50:16ZVortex induced vibrations of a pivoted finite height cylinder at low Reynolds number
Cajas García, Juan Carlos; Pastrana Maldonado, Daniel; Rodríguez Pérez, Ivette María; Lehmkuhl Barba, Oriol; Houzeaux, Guillaume; Vázquez, Mariano; Treviño Treviño, Cesar
The vortex induced vibrations (VIVs) of a pivoted cylinder with finite height have been numerically investigated. A mathematical model is introduced and described, and the resulting equations are numerically solved for low Reynolds number Re = 100, 200 and several combinations of the governing parameters. Results on the solid body trajectories, the maximum amplitude of the oscillations, the hydrodynamic force coefficients, the wake structure, and details on the vortex shedding near the cylinder are presented and discussed. The numerical results compare reasonably well with the canonical system of VIV of two-degrees of freedom circular cylinder in the laminar regime. Also, qualitative similarities with closely related VIV systems at larger Re suggest interesting lines of future research. Analytical approximations for limiting cases are done and an excellent agreement with the numerical results is obtained.
2021-06-07T07:50:16ZCajas García, Juan CarlosPastrana Maldonado, DanielRodríguez Pérez, Ivette MaríaLehmkuhl Barba, OriolHouzeaux, GuillaumeVázquez, MarianoTreviño Treviño, CesarThe vortex induced vibrations (VIVs) of a pivoted cylinder with finite height have been numerically investigated. A mathematical model is introduced and described, and the resulting equations are numerically solved for low Reynolds number Re = 100, 200 and several combinations of the governing parameters. Results on the solid body trajectories, the maximum amplitude of the oscillations, the hydrodynamic force coefficients, the wake structure, and details on the vortex shedding near the cylinder are presented and discussed. The numerical results compare reasonably well with the canonical system of VIV of two-degrees of freedom circular cylinder in the laminar regime. Also, qualitative similarities with closely related VIV systems at larger Re suggest interesting lines of future research. Analytical approximations for limiting cases are done and an excellent agreement with the numerical results is obtained.A review of recent progress in thermoelectric materials through computational methodsGutiérrez Moreno, J. JulioCao, JiangFronzi, MarcoAssadi, M. Hussein N.http://hdl.handle.net/2117/3465722021-06-06T21:11:11Z2021-06-03T10:43:26ZA review of recent progress in thermoelectric materials through computational methods
Gutiérrez Moreno, J. Julio; Cao, Jiang; Fronzi, Marco; Assadi, M. Hussein N.
Reducing our overwhelming dependence on fossil fuels requires groundbreaking innovations in increasing our efficiency in energy consumption for current technologies and moving towards renewable energy sources. Thermoelectric materials can help in achieving both goals. Moreover, because of recent advances in high-performance computing, researchers more increasingly rely on computational methods in discovering new thermoelectric materials with economically feasible performance. In this article, significant thermoelectric materials discovered through these computational methods are systematically reviewed. Furthermore, the primary computational tools that aid the design of the next-generation thermoelectric materials are introduced and discussed. These techniques include various levels of density functional theory, electronic transport simulations, and phonon calculations.
2021-06-03T10:43:26ZGutiérrez Moreno, J. JulioCao, JiangFronzi, MarcoAssadi, M. Hussein N.Reducing our overwhelming dependence on fossil fuels requires groundbreaking innovations in increasing our efficiency in energy consumption for current technologies and moving towards renewable energy sources. Thermoelectric materials can help in achieving both goals. Moreover, because of recent advances in high-performance computing, researchers more increasingly rely on computational methods in discovering new thermoelectric materials with economically feasible performance. In this article, significant thermoelectric materials discovered through these computational methods are systematically reviewed. Furthermore, the primary computational tools that aid the design of the next-generation thermoelectric materials are introduced and discussed. These techniques include various levels of density functional theory, electronic transport simulations, and phonon calculations.A FE2 multi-scale implementation for modeling composite materials on distributed architecturesGiuntoli, GuidoAguilar Mena, JimmyVázquez, MarianoOller Martínez, Sergio HoracioHouzeaux, Guillaumehttp://hdl.handle.net/2117/3461812021-06-06T06:55:57Z2021-05-27T08:00:59ZA FE2 multi-scale implementation for modeling composite materials on distributed architectures
Giuntoli, Guido; Aguilar Mena, Jimmy; Vázquez, Mariano; Oller Martínez, Sergio Horacio; Houzeaux, Guillaume
This work investigates the accuracy and performance of a FE2 multi-scale implementation used to predict the behavior of composite materials. The equations are formulated assuming the small deformations solid mechanics approach in non-linear material models with hardening plasticity. The uniform strain boundary conditions are applied for the macro-to-micro transitions. A parallel algorithm was implemented in order to solve large engineering problems. The scheme proposed takes advantage of the domain decomposition method at the macro-scale and the coupling between each subdomain with a micro-scale model. The precision of the method is validated with a composite material problem and scalability tests are performed for showing the efficiency.
2021-05-27T08:00:59ZGiuntoli, GuidoAguilar Mena, JimmyVázquez, MarianoOller Martínez, Sergio HoracioHouzeaux, GuillaumeThis work investigates the accuracy and performance of a FE2 multi-scale implementation used to predict the behavior of composite materials. The equations are formulated assuming the small deformations solid mechanics approach in non-linear material models with hardening plasticity. The uniform strain boundary conditions are applied for the macro-to-micro transitions. A parallel algorithm was implemented in order to solve large engineering problems. The scheme proposed takes advantage of the domain decomposition method at the macro-scale and the coupling between each subdomain with a micro-scale model. The precision of the method is validated with a composite material problem and scalability tests are performed for showing the efficiency.A city of cities: Measuring how 15-minutes urban accessibility shapes human mobility in BarcelonaGraells Garrido, EduardoSerra Burriel, FeliuRowe, FranciscoCucchietti, FernandoReyes Valenzuela, Patricio Alejandrohttp://hdl.handle.net/2117/3460972021-05-30T22:24:20Z2021-05-26T11:14:31ZA city of cities: Measuring how 15-minutes urban accessibility shapes human mobility in Barcelona
Graells Garrido, Eduardo; Serra Burriel, Feliu; Rowe, Francisco; Cucchietti, Fernando; Reyes Valenzuela, Patricio Alejandro
As cities expand, human mobility has become a central focus of urban planning and policy making to make cities more inclusive and sustainable. Initiatives such as the “15-minutes city” have been put in place to shift the attention from monocentric city configurations to polycentric structures, increasing the availability and diversity of local urban amenities. Ultimately they expect to increase local walkability and increase mobility within residential areas. While we know how urban amenities influence human mobility at the city level, little is known about spatial variations in this relationship. Here, we use mobile phone, census, and volunteered geographical data to measure geographic variations in the relationship between origin-destination flows and local urban accessibility in Barcelona. Using a Negative Binomial Geographically Weighted Regression model, we show that, globally, people tend to visit neighborhoods with better access to education and retail. Locally, these and other features change in sign and magnitude through the different neighborhoods of the city in ways that are not explained by administrative boundaries, and that provide deeper insights regarding urban characteristics such as rental prices. In conclusion, our work suggests that the qualities of a 15-minutes city can be measured at scale, delivering actionable insights on the polycentric structure of cities, and how people use and access this structure.
2021-05-26T11:14:31ZGraells Garrido, EduardoSerra Burriel, FeliuRowe, FranciscoCucchietti, FernandoReyes Valenzuela, Patricio AlejandroAs cities expand, human mobility has become a central focus of urban planning and policy making to make cities more inclusive and sustainable. Initiatives such as the “15-minutes city” have been put in place to shift the attention from monocentric city configurations to polycentric structures, increasing the availability and diversity of local urban amenities. Ultimately they expect to increase local walkability and increase mobility within residential areas. While we know how urban amenities influence human mobility at the city level, little is known about spatial variations in this relationship. Here, we use mobile phone, census, and volunteered geographical data to measure geographic variations in the relationship between origin-destination flows and local urban accessibility in Barcelona. Using a Negative Binomial Geographically Weighted Regression model, we show that, globally, people tend to visit neighborhoods with better access to education and retail. Locally, these and other features change in sign and magnitude through the different neighborhoods of the city in ways that are not explained by administrative boundaries, and that provide deeper insights regarding urban characteristics such as rental prices. In conclusion, our work suggests that the qualities of a 15-minutes city can be measured at scale, delivering actionable insights on the polycentric structure of cities, and how people use and access this structure.The making of the New European Wind Atlas – Part 2: Production and evaluationDörenkämper, MartinOlsen, Bjarke T.Witha, BjörnHahmann, Andrea N.Barcons, Jordihttp://hdl.handle.net/2117/3460442021-05-30T22:32:28Z2021-05-25T09:29:31ZThe making of the New European Wind Atlas – Part 2: Production and evaluation
Dörenkämper, Martin; Olsen, Bjarke T.; Witha, Björn; Hahmann, Andrea N.; Barcons, Jordi
This is the second of two papers that document the creation of the New European Wind Atlas (NEWA). In Part 1, we described the sensitivity experiments and accompanying evaluation done to arrive at the final mesoscale model setup used to produce the mesoscale wind atlas. In this paper, Part 2, we document how we made the final wind atlas product, covering both the production of the mesoscale climatology generated with the Weather Research and Forecasting (WRF) model and the microscale climatology generated with the Wind Atlas Analysis and Applications Program (WAsP). The paper includes a detailed description of the technical and practical aspects that went into running the mesoscale simulations and the downscaling using WAsP. We show the main results from the final wind atlas and present a comprehensive evaluation of each component of the NEWA model chain using observations from a large set of tall masts located all over Europe. The added value of the WRF and WAsP downscaling of wind climatologies is evaluated relative to the performance of the driving ERA5 reanalysis and shows that the WRF downscaling reduces the mean wind speed bias and spread relative to that of ERA5 from −1.50±1.30 to 0.02±0.78 m s−1. The WAsP downscaling has an added positive impact relative to that of the WRF model in simple terrain. In complex terrain, where the assumptions of the linearized flow model break down, both the mean bias and spread in wind speed are worse than those from the raw mesoscale results.
2021-05-25T09:29:31ZDörenkämper, MartinOlsen, Bjarke T.Witha, BjörnHahmann, Andrea N.Barcons, JordiThis is the second of two papers that document the creation of the New European Wind Atlas (NEWA). In Part 1, we described the sensitivity experiments and accompanying evaluation done to arrive at the final mesoscale model setup used to produce the mesoscale wind atlas. In this paper, Part 2, we document how we made the final wind atlas product, covering both the production of the mesoscale climatology generated with the Weather Research and Forecasting (WRF) model and the microscale climatology generated with the Wind Atlas Analysis and Applications Program (WAsP). The paper includes a detailed description of the technical and practical aspects that went into running the mesoscale simulations and the downscaling using WAsP. We show the main results from the final wind atlas and present a comprehensive evaluation of each component of the NEWA model chain using observations from a large set of tall masts located all over Europe. The added value of the WRF and WAsP downscaling of wind climatologies is evaluated relative to the performance of the driving ERA5 reanalysis and shows that the WRF downscaling reduces the mean wind speed bias and spread relative to that of ERA5 from −1.50±1.30 to 0.02±0.78 m s−1. The WAsP downscaling has an added positive impact relative to that of the WRF model in simple terrain. In complex terrain, where the assumptions of the linearized flow model break down, both the mean bias and spread in wind speed are worse than those from the raw mesoscale results.Specific heat effects in two-dimensional shock refractionsMartínez Ruiz, DanielHuete, CesarMartínez Ferrer, Pedro J.Mira Martinez, Danielhttp://hdl.handle.net/2117/3459652021-05-23T20:11:27Z2021-05-20T14:27:46ZSpecific heat effects in two-dimensional shock refractions
Martínez Ruiz, Daniel; Huete, Cesar; Martínez Ferrer, Pedro J.; Mira Martinez, Daniel
Compressible mixtures in supersonic flows are subject to significant temperature changes via shock waves and expansions, which affect several properties of the flow. Besides the widely studied variable transport effects such as temperature-dependent viscosity and conductivity, vibrational and rotational molecular energy storage is also modified through the variation of the heat capacity cp and heat capacity ratio γ, especially in hypersonic flows. Changes in the composition of the mixture may also modify its value through the species mass fraction Yα, thereby affecting the compression capacity of the flow. Canonical configurations are studied here to explore their sharply conditioned mechanical equilibrium under variations of these thermal models. In particular, effects of cp(T,Yα) and γ(T,Yα) on the stability of shock-impinged supersonic shear and mixing layers are addressed, on condition that a shock wave is refracted. It is found that the limits defining regular structures are affected (usually broadened out) by the dependence of heat capacities with temperature. Theoretical and high-fidelity numerical simulations exhibit a good agreement in the prediction of regular shock reflections and their post-shock aerothermal properties.
2021-05-20T14:27:46ZMartínez Ruiz, DanielHuete, CesarMartínez Ferrer, Pedro J.Mira Martinez, DanielCompressible mixtures in supersonic flows are subject to significant temperature changes via shock waves and expansions, which affect several properties of the flow. Besides the widely studied variable transport effects such as temperature-dependent viscosity and conductivity, vibrational and rotational molecular energy storage is also modified through the variation of the heat capacity cp and heat capacity ratio γ, especially in hypersonic flows. Changes in the composition of the mixture may also modify its value through the species mass fraction Yα, thereby affecting the compression capacity of the flow. Canonical configurations are studied here to explore their sharply conditioned mechanical equilibrium under variations of these thermal models. In particular, effects of cp(T,Yα) and γ(T,Yα) on the stability of shock-impinged supersonic shear and mixing layers are addressed, on condition that a shock wave is refracted. It is found that the limits defining regular structures are affected (usually broadened out) by the dependence of heat capacities with temperature. Theoretical and high-fidelity numerical simulations exhibit a good agreement in the prediction of regular shock reflections and their post-shock aerothermal properties.High order methods for acoustic scattering: Coupling farfield expansions ABC with deferred-correction methodsVillamizar, VianeyGrundvig, DaneRojas, OtilioAcosta, Sebastianhttp://hdl.handle.net/2117/3459282021-05-23T20:11:22Z2021-05-19T16:23:24ZHigh order methods for acoustic scattering: Coupling farfield expansions ABC with deferred-correction methods
Villamizar, Vianey; Grundvig, Dane; Rojas, Otilio; Acosta, Sebastian
Arbitrary high order numerical methods for time-harmonic acoustic scattering problems originally defined on unbounded domains are constructed. This is done by coupling recently developed high order local absorbing boundary conditions (ABCs) with finite difference methods for the Helmholtz equation. These ABCs are based on exact representations of the outgoing waves by means of farfield expansions. The finite difference methods, which are constructed from a deferred-correction (DC) technique, approximate the Helmholtz equation and the ABCs, with the appropriate number of terms, to any desired order. As a result, high order numerical methods with an overall order of convergence equal to the order of the DC schemes are obtained. A detailed construction of these DC finite difference schemes is presented. Additionally, a rigorous proof of the consistency of the DC schemes with the Helmholtz equation and the ABCs in polar coordinates is also given. The results of several numerical experiments corroborate the high order convergence of the novel method.
2021-05-19T16:23:24ZVillamizar, VianeyGrundvig, DaneRojas, OtilioAcosta, SebastianArbitrary high order numerical methods for time-harmonic acoustic scattering problems originally defined on unbounded domains are constructed. This is done by coupling recently developed high order local absorbing boundary conditions (ABCs) with finite difference methods for the Helmholtz equation. These ABCs are based on exact representations of the outgoing waves by means of farfield expansions. The finite difference methods, which are constructed from a deferred-correction (DC) technique, approximate the Helmholtz equation and the ABCs, with the appropriate number of terms, to any desired order. As a result, high order numerical methods with an overall order of convergence equal to the order of the DC schemes are obtained. A detailed construction of these DC finite difference schemes is presented. Additionally, a rigorous proof of the consistency of the DC schemes with the Helmholtz equation and the ABCs in polar coordinates is also given. The results of several numerical experiments corroborate the high order convergence of the novel method.