Garcia Gonzalez, Fernando
Total activity: 10
Department
Department of Applied Physics
School
Barcelona School of Civil Engineering (ETSECCPB)
E-mail
fernando.garcia-gonzalezestudiant.upc.edu
Contact details
UPC directory Open in new window

Graphic summary
  • Show / hide key
  • Information


Scientific and technological production
  •  

1 to 10 of 10 results
  • A comparison of high-order time integrators for highly supercritical thermal convection in rotating spherical shells

     Garcia Gonzalez, Fernando; Net Marce, Marta; Sanchez Umbria, Juan Jose
    Lecture notes in computational science and engineering
    Date of publication: 2014-01-10
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    The efficiency of implicit and semi-implicit time integration codes based on backward differentiation and extrapolation formulas for the solution of the three-dimensional Boussinesq thermal convection equations in rotating spherical shells was studied in Garcia et al. (J Comput Phys 229:7997¿8010, 2010) at weakly supercritical Rayleigh numbers R, moderate (10-3) and low (10-4) Ekman numbers, E, and Prandtl number s = 1. The results presented here extend the previous study and focus on the effect of s and R by analyzing the efficiency of the methods for obtaining solutions at TeX, s = 0. 1 and low and high supercritical R. In the first case (quasiperiodic solutions) the decrease of one order of magnitude does not change the results significantly. In the second case (spatio-temporal chaotic solutions) the differences in the behavior of the semi-implicit codes due to the different treatment of the Coriolis term disappear because the integration is dominated by the nonlinear terms. As in Garcia et al. (J Comput Phys 229:7997¿8010, 2010), high order methods, either with or without time step and order control, increase the efficiency of the time integrators and allow to obtain more accurate solutions.

  • Exponential versus IMEX high-order time integrators for thermal convection in rotating spherical shells

     Garcia Gonzalez, Fernando; Bonaventura, Luca; Net Marce, Marta; Sanchez Umbria, Juan Jose
    Journal of computational physics
    Date of publication: 2014-01-24
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    We assess the accuracy and efficiency of several exponential time integration methods coupled to a spectral discretization of the three-dimensional Boussinesq thermal convection equations in rotating spherical shells. Exponential methods are compared to implicit¿explicit (IMEX) multi-step methods already studied previously in [1]. The results of a wide range of numerical simulations highlight the superior accuracy of exponential methods for a given time step, especially when employed with large time steps and at low Ekman number. However, presently available implementations of exponential methods appear to be in general computationally more expensive than those of IMEX methods and further research is needed to reduce their computational cost per time step. A physically justified extrapolation argument suggests that some exponential methods could be the most efficient option for integrating flows near Earth's outer core conditions.

  • Access to the full text
    Computation of azimuthal waves and their stability in thermal convection in rotating spherical shells with application to the study of a double-Hopf bifurcation  Open access

     Sanchez Umbria, Juan Jose; Garcia Gonzalez, Fernando; Net Marce, Marta
    Physical review E: statistical, nonlinear, and soft matter physics
    Date of publication: 2013-03-22
    Journal article

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    A methodology to compute azimuthal waves, appearing in thermal convection of a pure fluid contained in a rotating spherical shell, and to study their stability is presented. It is based on continuation, Newton-Krylov, and Arnoldi methods. An application to the study of a double-Hopf bifurcation of the basic state is shown for Ekman and Prandtl numbers E=10-4 and =0.1, respectively, radius ratios n[0.32,0.35], Rayleigh numbers R[1.8105,6105], and nonslip and perfectly conducting boundary conditions. The knowledge of the bifurcation diagrams, including the unstable solutions, allows one to understand the coexistence of stable thermal Rossby waves of different azimuthal wave numbers at some parameter regions, and the origin of some new intermittent solutions found, as trajectories close to heteroclinic chains. Moreover, the structure of the eigenfunctions at the secondary bifurcations explains the existence of the amplitude and shape modulated waves.

    A methodology to compute azimuthal waves, appearing in thermal convection of a pure fluid contained in a rotating spherical shell, and to study their stability is presented. It is based on continuation, Newton-Krylov, and Arnoldi methods. An application to the study of a double-Hopf bifurcation of the basic state is shown for Ekman and Prandtl numbers E=10−4 and σ=0.1, respectively, radius ratios η∈[0.32,0.35], Rayleigh numbers R∈[1.8×105,6×105], and nonslip and perfectly conducting boundary conditions. The knowledge of the bifurcation diagrams, including the unstable solutions, allows one to understand the coexistence of stable thermal Rossby waves of different azimuthal wave numbers at some parameter regions, and the origin of some new intermittent solutions found, as trajectories close to heteroclinic chains. Moreover, the structure of the eigenfunctions at the secondary bifurcations explains the existence of the amplitude and shape modulated waves.

  • Access to the full text
    Numerical study of the onset of thermosolutal convection in rotating spherical shells  Open access

     Net Marce, Marta; Garcia Gonzalez, Fernando; Sanchez Umbria, Juan Jose
    Physics of fluids
    Date of publication: 2012-06-01
    Journal article

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    The influence of an externally enforced compositional gradient on the onset of convection of a mixture of two components in a rotating fluid spherical shell is studied for Ekman numbers E = 10−3 and E = 10−6, Prandtl numbers σ = 0.1, 0.001, Lewis numbers τ = 0.01, 0.1, 0.8, and radius ratio η = 0.35. The Boussinesq approximation of the governing equations is derived by taking the denser component of the mixture for the equation of the concentration. Differential and internal heating, an external compositional gradient, and the Soret and Dufour effects are included in the model. By neglecting these two last effects, and by considering only differential heating, it is found that the critical thermal Rayleigh number Rec depends strongly on the direction of the compositional gradient. The results are compared with those obtained previously for pure fluids of the same σ. The influence of the mixture becomes significant when the compositional Rayleigh number Rc is at least of the same order of magnitude as the known Rec computed without mixture. For positive and sufficiently large compositional gradients, Rec decreases and changes sign, indicating that the compositional convection becomes the main source of instability. Then the critical wave number mc decreases, and the drifting waves slow down drastically giving rise to an almost stationary pattern of convection. Negative gradients delay the onset of convection and determine a substantial increase of mc and ωc for Rc sufficiently high. Potential laws are obtained numerically from the dependence of Rec and of the critical frequency ωc on Rc, for the moderate and small Ekman numbers explored.

  • Thermal convection in rotating spherical shells

     Garcia Gonzalez, Fernando
    Defense's date: 2012-11-30
    Department of Applied Physics, Universitat Politècnica de Catalunya
    Theses

     Share Reference managers Reference managers Open in new window

  • Cálculo numérico de variedades invariantes en EDPS disipativas. Aplicaciones a la convección térmica

     Net Marce, Marta; Garcia Gonzalez, Fernando; Sanchez Umbria, Juan Jose
    Participation in a competitive project

     Share

  • Access to the full text
    A comparison of high-order time integrators for thermal convection in rotating spherical shells  Open access

     Garcia Gonzalez, Fernando; Net Marce, Marta; Garcia-Archilla, Bosco; Sanchez Umbria, Juan Jose
    Journal of computational physics
    Date of publication: 2010-10-01
    Journal article

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    A numerical study of several time integration methods for solving the threedimensional Boussinesq thermal convection equations in rotating spherical shells is presented. Implicit and semi-implicit time integration techniques based on backward differentiation and extrapolation formulae are considered. The use of Krylov techniques allows the implicit treatment of the Coriolis term with low storage requirements. The codes are validated with a known benchmark, and their efficiency is studied. The results show that the use of high order methods, especially those with time step and order control, increase the efficiency of the time integration, and allows to obtain more accurate solutions.

    Postprint (author’s final draft)

  • Access to the full text
    Estudio comparativo de la eficiencia de diversos integradores para la simulación numérica de la convección térmica en geometría esférica y rotación  Open access

     Garcia Gonzalez, Fernando; Net Marce, Marta; Garcia-Archilla, Bosco; Sanchez Umbria, Juan Jose
    Congreso de Ecuaciones Diferenciales y Aplicaciones / Congreso de Matemática Aplicada
    Presentation's date: 2009-09-24
    Presentation of work at congresses

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    The efficiency of implicit and semi-implicit time integration codes based on bac- kward differentiation and extrapolation formulae for the solution of the three-dimensio- nal Boussinesq thermal convection equations in rotating spherical shells is studied. The use of Krylov techniques allows the implicit treatment of the Coriolis term with low storage requirements. The results show that high order methods, either with or without time step and order control, increase the efficiency of the time integrators, and allow to obtain more accurate solutions.

  • Access to the full text
    Antisymmetric polar modes of thermal convection in rotating spherical fluid shells at high Taylor numbers  Open access

     Garcia Gonzalez, Fernando; Net Marce, Marta; Sanchez Umbria, Juan Jose
    Physical review letters
    Date of publication: 2008-11-03
    Journal article

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    The onset of thermal convection in a rotating spherical shell of intermediate radius ratio ¼ 0:4 is studied numerically for Taylor numbers Ta 1011 and the Prandtl number of the liquid sodium ( ¼ 0:01). For the first time, it is shown that at very high Taylor numbers the first unstable mode can be antisymmetric with respect to the equator and confined inside a cylinder tangent to the inner sphere at the equator (polar mode). The exponent of the power law determined from the asymptotic dependence of the critical Rayleigh number for very high Ta is 0.57, lower than 2=3, given theoretically for the spiraling columnar modes, and than 0.63, found numerically for the outer equatorially attached modes

  • A note on Weak Stability Boundaries

     Garcia Gonzalez, Fernando; Gómez Muntané, Gerard
    Celestial mechanics and dynamical astronomy
    Date of publication: 2006-12-08
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    This paper is devoted to clarify the algorithmic definition of the weak stability boundary in the framework of the planar Restricted Three Body Problem. The role of the invariant hyperbolic manifolds associated to the central manifolds of the libration points L1 and L2, as boundary of the weak stability region, is shown