Magíster en Ciencias de la Ingeniería, Mención Matemáticas Aplicadas. Ingeniero Civil Matemático / Si se representa un medio poroso por medio de una red de tubos por los cuales puede circular un fluido, el problema de escurrimiento a través de un medio poroso se transforma en un problema de cálculo de caudales en una red de tubos interconectadas. La solución de este último problema hace uso de las leyes de conservación clásicas en mecánica de fluidos, complementadas con leyes de comportamiento adecuadas, como lo son los fluidos newtonianos y fluidos tipo Bingham. En cada tubo se hacen simplificaciones tales como suponer un problema cuasiestático, lo que permite trabajar con problemas unidimensionales (acoplados en las conexiones de los tubos) de resolución numérica simple. A continuación, se programa un algoritmo que modela el flujo por una red de tubos capilares que simula un medio poroso. El programa se caracteriza por su versatilidad para poder simular distintos fluidos, y mallas con distintas geometrías. Se realizan distintos experimentos donde se compara con datos experimentales y con la literatura asociada, con el fin de darle mayor robustez al programa.

Publicación a texto completo no autorizada por el autor / El documento digital no refiere asesor / Establece la dinámica del oleaje en la bahía de Miraflores utilizando formulaciones analíticas y empíricas para estudiar la transformación del oleaje desde aguas profundas hasta agua poco profundas. Calcula todos los parámetros que caracteriza a las rompientes tales como: altura de ola y profundidad rompiente, ángulo de incidencia de la ola rompiente con la línea de costa, y tipos de ola rompiente: spilling (derrumbante), plunging (precipitante) y surging (surgente). Determina la variación de la altura de ola después de romper. Se realiza una comparación entre distintos métodos y criterios para determinar estos parámetros. Describe la transformación del oleaje desde aguas profundas hasta zonas poco profundas: calcula el coeficiente de asomeramiento mediante la teoría lineal y cnoidal, calcula el coeficiente de refracción mediante la teoría lineal, calcula el coeficiente de pérdida de energía por fricción de fondo. Estima la tasa de transporte longitudinal de sedimentos. / Tesis

La Madera es un material natural utilizado en elementos estructurales dentro del ámbito de la Ingeniería Civil. Por tratarse de un biomaterial, tiene características complejas y variables. Particularmente, la madera proveniente de la especie Eucalyptus grandis cultivado en la Mesopotamia Argentina es una de las aceptadas para el uso estructural en el país por el Reglamento Argentino de Estructuras de Madera CIRSOC 601. Las piezas de este material presentan una gran variabilidad en sus propiedades en relación a otros materiales de uso estructural, la cual puede ser explicada tanto por las características de su micro y macro estructura, la incidencia del clima y el suelo, los procesos de producción, y la presencia de defectos que afectan su comportamiento mecánico. Dentro de estos últimos, de singular importancia para el uso estructural de este material, se destacan la presencia de médula y los nudos. Ambos, se encuentran limitados en los métodos de clasificación por resistencia de piezas estructurales. Debido a estos factores, la cuantificación de los efectos de la variabilidad de las propiedades materiales en la respuesta estructural es de interés dentro del ámbito del diseño ingenieril. Para ello, en el presente trabajo, se combinaron tres ejes temáticos principales: el material, el modelado estocástico y el modelado numérico. Estos ejes temáticos fueron empleados en la construcción de modelos numérico-estocásticos aplicables a elementos estructurales desarrollados en el entorno del Método de los Elementos Finitos. Estos modelos han sido calibrados y validados con resultados experimentales existentes mediante simulaciones realizadas a través del Método de Monte Carlo. Luego, fueron utilizados para la realización de estudios de propagación de incertidumbres en la respuesta estructural. Se muestra que un modelo estocástico del material que representa la variabilidad longitudinal de las propiedades del mismo provee una respuesta estructural más próxima a la observada en resultados experimentales. Además, como complemento al desarrollo numérico, se realizaron nuevos ensayos experimentales, dinámicos y estáticos, tanto en tablas como en vigas aserradas de dimensiones estructurales. Con el objetivo de validar las hipótesis asumidas para el desarrollo de los modelos estocásticos. Sumado a esto, dado que no se conocen valores precisos para este material se obtuvieron coeficientes de amortiguamiento, los cuales son necesarios para el estudio de la respuesta estructural bajo acciones dinámicas. Finalmente, y a los efectos de aplicar el modelo desarrollado a una estructura, se realizó un estudio numérico sobre un modelo de puente peatonal de madera con vigas laminadas bajo cargas dinámicas determinísticas y estocásticas producidas por peatones. El mismo, se llevó a cabo para el estudio de las condiciones de servicio. La cuantificación de incertidumbre pretende ampliar la información resultante en la predicción de la respuesta así como la sensibilidad de la misma a la variación de uno o más parámetros del sistema. En esta investigación, esto se logra a través del desarrollo y aplicación de modelos numérico-estocásticos aplicables a elementos estructurales y estructuras de madera aserrada. / Timber is a natural material employed for structural elements in the field of the Civil Engineering. Being a biomaterial, it has complex and random characteristics. Particularly, the Eucalyptus grandis cultivated in the Mesopotamic provinces of Argentine is one of the timber species accepted for the structural use by the Argentinean Standard of Timber Structures CIRSOC 601. The pieces of this material have a high variability in their properties in relation to other materials for structural use, which can be explained by the characteristics of its micro and macro structure, the incidence of the climate and soil, the production processes, and the presence of defects that affect its mechanical behavior. Within the latter, of singular importance for the structural use of this material, the presence of pith and knots are emphasized. Both types of defects are limited in the strength grading classification methods of structural elements. Due to these factors, the quantification of the influence of the timber material properties variability in the structural response is an important issue in the structural design. In order to quantify this influence, three thematic areas are combined: the material, the stochastic modelling and the numerical modelling. These areas were employed in the construction of numerical-stochastic models applicable to structural elements developed in the environment of the Finite Element Method. These models have been calibrated and validated with existing experimental results through simulations performed through the Monte Carlo Method. Then, they were used to carry out uncertainty propagation studies on the structural response. It is shown that a stochastic model of the material which represents the longitudinal variability of the properties provides a structural response closer to the experimental results. In addition, as a complement to the numerical development, new experimental dynamic and static tests were performed in beams and boards of structural dimensions. The aim of these tests is the validation of the assumptions taken into account in the development of the stochastic models. In addition to this, since no precise values are known for this material, damping coefficients were obtained, which are necessary for the study of the structural response under dynamic actions. Finally, and in order to apply the developed model to a structure, a numerical study of the structural behaviour of a timber footbridge under the action of deterministic and stochastic models of walking loads induced by pedestrians was performed. It was carried out in order to evaluate the serviceability conditions. The uncertainty quantification aims to extend the information resulting in the prediction of the response as well as the sensitivity to the variation of one or more parameters of the system. In this work, this is achieved through the development and application of numerical-stochastic models applicable to structural timber elements and structures.

Orientador: Amauri Hassui / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-21T02:13:50Z (GMT). No. of bitstreams: 1 Castro_MarceloHirai_M.pdf: 4008200 bytes, checksum: a0c71450e4b5c3ac934e22f6e3cc9f4b (MD5) Previous issue date: 2012 / Resumo: A inovação é a chave para o diferencial da indústria no mercado competitivo atual. Fazer diferente, melhor e mais barato sempre agrega valor ao produto. Por outro lado, a resistência às mudanças faz parte do comportamento humano, seja na sociedade, ou na indústria. O objetivo principal deste trabalho é estudar a viabilidade de se fazer furos não cilíndricos com uma ferramenta em rotação utilizando máquinas-ferramenta CNC. Estes sempre foram obstáculos na indústria de fabricação pela dificuldade de obtenção. Sendo assim, o processo apresentado nesse trabalho busca usinar um furo quadrado com a mesma configuração de ferramentas e fixação de peças utilizadas para um furo cilíndrico, tão comum à indústria. Isso visa reduzir custos de fabricação e facilitar a obtenção de furos não cilíndricos com os recursos de equipamentos programáveis existentes no chão de fábrica. A partir de um modelo matemático, a máquina é programada para executar movimentos sincronizados entre a rotação e os deslocamentos do eixo da ferramenta. Diversas programações foram testadas em diferentes máquinas ferramenta CNC, com diferentes comandos, com e sem recursos de alta velocidade, com diferentes quantidades de pontos programados, interpolações lineares e circulares. A partir dos experimentos, foram analisadas as configurações das máquinas, a capacidade de execução do controlador numérico instalado, a velocidade desenvolvida pela máquina na operação e a geometria obtida dos furos quadrados. Os experimentos apontam que para executar o furo quadrado, os comandos devem ser capazes de controlar um eixo rotativo como eixo propriamente dito, além de sincronizarem mais três eixos lineares. O crítico para este tipo de operação é a velocidade desempenhada pela máquina, sendo o fator mais significativo, o tempo que o comando precisa para o processamento do bloco de programação, independentemente da ativação de recursos de alta velocidade, interpolações lineares ou circulares. Verificou-se também que quanto maior a quantidade de pontos para compor a trajetória do centro da ferramenta, melhor será a geometria obtida do furo. Sendo assim, as máquinas ferramenta CNC que estão em operação na indústria são capazes de realizar a furação não cilíndrica de forma muito próxima à furação cilíndrica quanto ao tempo de usinagem e ao sistema de fixação da ferramenta e da peça, sendo necessária uma programação sincronizada do eixo spindle presente em toda e qualquer máquina ferramenta / Abstract: The innovation is the key for the differential of the industry on the competitive market nowadays. To make different, better and cheaper always it is worth a lot to the product. By the other side, the resistance to changes makes part of the human behavior, in society, or in the industry. The main goal of this work is to study the feasibility to machining of non-cylindrical holes by a tool in rotation using CNC machine tools. That always has been offering obstacles to manufacturing industry by the difficulty of fabrication. Thus, the process presented on this work will machining a square hole with the same configuration of tools and systems for fixation of the parts, used for a cylindrical drilling, so common to industry. This intends to reduce the costs of fabrication and make easier to get non-cylindrical drilling using resources of the programmable equipment which are working at the shops. From a math model on, the machine is programmed to execute synchronized movements between the rotation and the displacements of the axis of the tool. Different programs were tested in different CNC machine tools, with different controllers, with and no resources for high speed, with different quantity of programmed points, with linear and circular interpolations. From the experiments, it was analyzed the configuration of the machines, the capacity for execution of the controller installed, the speed developed by machine at the operation and the geometry of the square holes. The experiments show that to machine a square hole, the controllers must be able to control a rotating axis itself, beyond to synchronize three axis. The critical to this kind of operation is the speed developed by the machine, being the most significant factor, the time the controller needs to process the block of programming, independently the activation of resources for high speed, linear or circular interpolations. It was verified than as larger the amount of points to compose the trajectory of the center of tool, better will be the geometry of the hole. Thus, the CNC machine tools in operation in the industry are able to machine the non-cylindrical drilling in a way very closely the cylindrical drilling as the time of machining and the fixing system of the tool and parts, requiring a programming that synchronizes the main spindle installed in any machine tool / Mestrado / Materiais e Processos de Fabricação / Mestre em Engenharia Mecânica

Visualizações imersivas são muito valiosas para melhorar a compreensão de uma variedade de fenômenos físicos, que podem ser eventualmente modelados na forma discreta e simulados por computador. Dentre possíveis aplicações podemos utilizar a visualização imersiva como ferramenta pedagógica para percepção aumentada de tópicos complexos, ou como uma poderosa ferramenta de apoio analítico para ajudar os engenheiros a interpretarem os resultados de projetos. Este projeto de pesquisa aborda o uso da visualização imersiva de campos eletromagnéticos, especificamente os campos gerados por redes sem fio, largamente utilizadas no cotidiano como é o caso das redes IEEE 802.11 (Wi-Fi). Para tanto este trabalho propõe métodos novos para visualizar, em tres dimensões, campos eletromagnéticos variantes no tempo e distribuições de parâmetros interessantes relacionados a redes sem fio. Para atingir este objetivo, uma versão aprimorada do método de diferenças finitas no domínio do tempo (FDTD) é desenvolvido: o método FDTD de alta ordem e malha grosseira (Coarse Grid Higher Order FDTD, CGHO-FDTD). Portanto, soluções numéricas muito precisas, mais rápidas, e computacionalmente mais eficientes das equações de Maxwell podem ser obtidas. Os cálculos numéricos podem ser ainda mais rápidos pelo uso de computação paralela em um aglomerado de computadores. As características de domínio de tempo facilitam a criação de instantâneos de campos eletromagnéticos que estão se propagando, e desta maneira é possível criar figuras e animações tridimensionais que podem ser usadas para explicar alguns dos seguintes fenômenos físicos comuns em redes sem fio: difração, reflexão, e atenuação. Para que aumente a percepção física ainda mais, visualizações imersivas são feitas em um ambiente de realidade virtual. Por fim, a ferramenta desenvolvida também pode ser usada para criar distribuições muito detalhadas de parâmetros importantes que afetam o desempenho em uma rede sem fio. É mostrado que simulações de um ambiente fechado para prever a distribuição de potência de uma rede sem fio do tipo IEEE 802.11 (Wi-Fi), estão de acordo com as medidas. / Immersive visualizations are very valuable in order to improve the understanding of a variety of physical phenomena that can be modeled numerically and simulated by computers. Amongst the possible applications, we could utilize immersive visualizations as a pedagogical tool for enhanced perception of complex topics, or as a powerful tool that helps engineers interpret the outcome of simulations. This research project approaches the use of immersive visualizations of electromagnetic fields, especially fields generated by wireless networks widely utilized in the everyday life, as is the case for networks of the type IEEE 802.11 (Wi-Fi). For such a purpose this work proposes new methods to three-dimensionally visualize time-varying electromagnetic fields, and distributions of interesting parameters related to wireless networks. To achieve these objectives, a better version of the finite-difference time-domain (FDTD) method is developed: the Coarse Grid Higher Order FDTD (CGHO-FDTD) method. Thus highly accurate, faster and more computationally efficient numerical solutions of Maxwells equations can be obtained. The numerical calculations are made even faster by the use of parallel computing on a cluster of computers. The characteristics of the time domain facilitate the creation of snapshots of the propagating electromagnetic fields, and in this manner it is possible to create three-dimensional figures and animations that can be used to explain some of the following common physical effects found in wireless networks: diffraction, reflection, and attenuation. To further enhance the perception of the physics, immersive visualizations are carried out in a virtual reality environment. Finally, the developed tool can also be used to create highly detailed distributions of important parameters that affect the performance in wireless networks. It is shown that simulations to predict the power distribution in an indoor wireless network of the type IEEE 802.11 (Wi-Fi), agree very well with measurements.

Este trabalho apresenta um tratamento do Grupo de Renormalização Numérico (GRN) que permite a inclusão de assimetrias na banda de condução, como a provocada por campo magnetostático aplicado a sistemas de spins localizados em metais. Resultados para calor específico e suscetibilidade magnética de spin são obtidos para os modelos de Kondo de dois canais (MK2C), de Falikov-Kimball (MFK) e o modelo de nível ressonante na presença dessas assimetrias. São apresentadas soluções para lidar com as oscilações espúrias das curvas do calor específico, inerentes ao GRN. Tais abordagens, contribuíram de forma determinante para a viabilidade dos resultados apresentados aqui para essas grandezas. / This work presents an approach to the Numerical Renormalization Group which allows for asymmetries in the conduction band, as those brought about by magnetostatic field applied to a system of localized spin in metal. Results for specific heat and magnetic susceptibility are shown for the two-channel Kondo model, as well for the Falikov-Kimball and resonant level models, in the presence of such asymmetry. In addition, novel treatments were described in order to smooth out the curves of specific heat from the oscilating profile that comes along from NRG calculations. This was of outmost importance for the viability of the results presented here.

Atualmente, a competição entre empresas tem sido substituída pela competição entre cadeias de suprimentos, formadas por várias empresas. Um fator chave para garantir a rápida resposta da cadeia de suprimentos é a integração de sistemas. Esta integração deve acontecer não apenas entre os elos da cadeia, mas também dentro de cada empresa entre o chão de fábrica e o gerenciamento da produção. Desta forma, para obter sucesso num ambiente globalizado, torna-se necessário integrar o chão de fábrica a aplicativos MES (Manufacturing Execution System), ERP (Enterprise Resource Planning) e SCM (Supply Chain Management), possibilitando assim o compartilhamento de dados de produção entre empresas de uma mesma cadeia e garantindo a transparência das informações. Neste contexto, a proposta deste trabalho é aplicar o conceito de monitoramento remoto para o monitoramento da execução de pedidos de manufatura em máquinas CNC. Como resultado, a solução proposta proporciona a integração entre empresas, garantindo a transparência dos dados e informações sobre a produção.

Este projeto apresenta a construção de um modelo numérico do Sistema Aqüífero Guarani pelo método de elementos finitos para auxiliar o gerenciamento de recursos hídricos subterrâneos na área do projeto piloto de Ribeirão Preto, localizada na região nordeste do Estado de São Paulo. Dados topográficos, hidrológicos e hidrogeológicos da região de interesse foram obtidos principalmente de mapas (impressos e digitais) e poços cadastrados junto a órgãos governamentais, em cooperação com o grupo coordenador do projeto piloto de Ribeirão Preto. Os dados levantados foram armazenados, processados e analisados utilizando ferramenta SIG, seguida da geração de mapas interpolados com características hidrogeológicas. Esses dados foram transferidos para o software de simulação SPA, onde se deu a construção da malha de elementos finitos e a atribuição dos parâmetros e condições de contorno ao modelo. Em seguida, foi feita uma análise de sensibilidade do modelo, e testadas algumas configurações de parâmetros hidrogeológicos e de taxas de recarga. O modelo apresentou bons resultados e desempenho, e, como parte de uma constante evolução, são propostas visitas de campo para verificação de dados para o seguimento desse projeto. / This work presents a numerical model for the Guarani Aquifer System management in the Ribeirão Preto pilot project area, located in the Northeast of São Paulo State. Topographical, hydrological and hydrogeological data of the region were obtained from printed and digital maps and well database recorded in governmental departments, in cooperation with the Ribeirão Preto coordinating group of the pilot project. The collected data were stored, processed and analyzed making use of GIS tool, followed by the generation of interpolated maps with hydrogeological parameters. These data were transferred to the SPA simulation software, which is based on the finite element method. After mesh construction, and attribution of parameters and boundary conditions, a sensitivity analysis of the model was performed. Different hydrogeological parameter configurations were tested. The model proved to represent the established conceptual model adequately. However base flow data and field visits are necessary, in order to improve the model representativity.

Resultados do Grupo de Renormalização Numérico (GRN) para a condutância linear dependente da temperatura associada a corrente de tunelamento através de uma ponta de prova nas proximidades de uma impureza magnética são apresentados. Nós usamos o Modelo de Anderson de uma impureza para descrever o metal hospedeiro e um Hamiltoniano livre para simular a ponta de prova do MVT (Microscópio de Varredura por Tunelamento). O cálculo da condutância é obtida a partir da fórmula de Kubo com o Hamiltoniano de tunelamento tratado como uma perturbação com dois canais de tunelamento, ponta-impureza e ponta-substrato, com o objetivo de descrever esse sistema que está totalmente fora do equilíbrio. Esse cálculo é guiado pelo GRN de Wilson para determinar a fórmula da condutância em termos de densidades espectrais: a densidade local da impureza e a densidade relativa ao primeiro sítio de condução da rede tight-binding do GRN. Esse resultado para o operador do GRN transforma esse objeto teórico em uma quantidade mensurável. Mostramos sob condições especiais, que o gráfico da condutância em função da temperatura é uma curva universal. Como função da posição ponta-impureza, as correntes de tunelamento mostram oscilações de Friedel, que determinam o tamanho da nuvem Kondo. Finalmente, mostramos como função da energia da impureza, a corrente da impureza para a ponta mostra um platô de Kondo. A interferência entre essa corrente e a que flui da banda de condução para a ponta exibe anti-ressonâncias de Fano como as observadas em medidas espectroscópicas. / Numerical Renormalization Group (NRG) results for the temperature dependent linear conductance associated with the scanning-tunneling current through a probe near a magnetic impurity are reported. We used the Single Impurity Anderson Model to describe the host metal and a free electron Hamiltonian to simulate a STM (Scanning Tunneling Microscope) biased tip. The calculation of the conductance is obtained from the Kubo Formula with the Tunneling Hamiltonian treated as a perturbation with two tunneling channels, STM tip-impurity and STM tip-host metal, with the objective to describe this fully nonequilibrium system. This calculation is guided by Wilson\'s NRG to determine a conductance formula as a funciton of spectral densities: the local impurity density and the density relative to the first conduction site of the NRG tight-binding chain. This result for the NRG operator transforms this theoretical object into a measurable quantity. We show that, under special conditions, plotted as a function of temperature, this zero-bias conductance follows a universal curve. As a function of tip-impurity separation, the tunneling currents display Friedel Oscilations, which determine the size of the Kondo cloud. Finally, plotted as a function of impurity energy, the current from the impurity to the tip displays a Kondo plateau. The inferference between this current and that flowing from the conduction band to the tip displays Fano anti-ressonances analogous to those seen in spectroscopic measurements.

Dados de altura da superfície do mar medidos a partir de satélites altimétricos e saídas de dois modelos de circulação geral dos oceanos (OGCM for the Earth Simulator - OFES e Community Earth System Model - CESM) foram utilizados com o objetivo de verificar se modelos numéricos climáticos reproduzem ondas de Rossby de forma similar à observada em dados altimétricos no Atlântico Sul. Os modelos diferem quanto à forçante, o modelo OFES é forçado com dados do NECP-NCAR enquanto que o modelo CESM é acoplado com modelo atmosférico, de gelo e terrestre. Uma vez que essas ondas dependem da estrutura interna de densidade e da forçante do vento, podemos fazer inferências sobre a adequação desses dois fatores nos modelos à realidade observada pelo altímetro. Uma série de filtros de resposta impulsiva finita 2D (FIR-2D) foi aplicada aos dados de anomalia da altura da superfície do mar com o propósito de detectar as ondas de Rossby e as componentes de onda encontradas nos dados dos modelos foram comparadas com as do altímetro. Os dois modelos são capazes de reproduzir ondas de Rossby e apresentam melhores resultados em baixas latitudes. Porém o modelo OFES apresentou mais dificuldades para reproduzir os parâmetros de onda encontrados no altímetro apresentando diferenças médias de até 68% para a amplitude, 34% para o comprimento de onda e 38% para a velocidade de fase. Em contrapartida as maiores diferenças médias entre os parâmetros de onda calculados a partir dos dados altimétricos e do modelo CESM foram de 32% para a amplitude, 20% para o comprimento de onda e 20% para a velocidade de fase. Além disso, o modelo CESM foi capaz de reproduzir o sinal sazonal com uma correlação média de aproximadamente 0,7 com o sinal sazonal encontrado pelo altímetro em toda a bacia do oceano Atlântico Sul enquanto que o sinal sazonal do modelo OFES apresentou uma correlação média de 0,4 com o sinal encontrado pelo satélite. Estes resultados mostram que os dois modelos reproduzem o fenômeno satisfatoriamente, sendo o CESM melhor que o OFES. A diferença dos resultados deve estar ligada aos aspectos supracitados, especificamente à forçante atmosférica e estrutura de densidade na coluna d\'água nas latitudes ao sul de 20ºS. / Sea surface height data measured from altimetry satellites and outputs of two ocean general circulation models (OGCM for the Earth Simulator - OFES and Community Earth System Model - CESM) were used to determine whether climate numerical models reproduce the Rossby waves in a manner similar to those observed in altimetry record in the South Atlantic. The models differs in the forcing, OFES is forced with NCEP-NCAR data while CESM is coupled with atmospheric model, ice and land model. Because these waves depend on the internal density structure and wind forcing, we can make inferences about the suitability of these two factors in the models in comparison with the altimetry data used as a reference. A series of finite impulse response band-pass filters (FIR-2D) was applied to isolate the westward propagating signals corresponding to Rossby waves in the altimeter. Both models are able to reproduce Rossby waves and show better results in low latitudes. However, the OFES model presented more difficulty to reproduce the wave parameters found in altimeter with differences of up to 68% for amplitude, 34% for the wavelength and 38% for phase velocity. By contrast the greatest differences between the wave parameters computed from the altimeter data and the CESM model were 32% for amplitude, 20% for the wavelength and 20% for phase velocity. Furthermore, the CESM model was capable of reproducing the seasonal signal correlation with an average of approximately 0.7 with altimeter\'s seasonal signal found throughout the South Atlantic basin, while the OFES\'s seasonal signal showed an average correlation 0.4 with the signal found by the satellite. These results indicate that both models can reproduce the phenomenon satisfactorily, the CESM better than OFES. The difference between the results should be related with the aspects cited above, specifically with the atmospheric forcing and the density structure in the water column in latitudes southern then 20ºS.