No description available.

O trabalho apresenta uma modificação do modelo do número de onda cumulativo (CW – Cumulative Wavenumber) para atender ao balanço de energia na determinação da transferência radiante de gases participantes não uniformes. De um modo geral, o modelo CW fornece resultados que se mostram precisos para o divergente do fluxo de calor radiante (ou taxa volumétrica de geração de energia) quando comparados à solução benchmark (integração linhapor- linha). No entanto, conforme é apresentado nesta pesquisa, uma hipótese importante feita pelo modelo o torna não conservativo e, não garantindo o balanço de energia radiante quando um meio não uniforme é considerado. Como conseqüência, o fluxo de calor radiante pode apresentar desvios consideráveis em relação à solução exata. O modelo do número de onda cumulativo modificado (CWM) foi desenvolvido de forma a manter o mesmo valor para a taxa volumétrica de geração de energia, mas satisfazendo ao balanço de energia radiante. O modelo CWM é aplicado juntamente com o método das ordenadas discretas para resolver a transferência de calor radiante em sistemas unidimensionais formados por placas planas infinitas e paralelas, contendo camadas uniformes e não uniformes de gases típicos da combustão de metano ou óleo combustível. As paredes são negras para a radiação térmica. A aquisição de dados das linhas de absorção espectrais é feita pela utilização dos bancos de dados HITRAN e HITEMP. São apresentados os bancos de dados, bem como a forma em que é efetuada a extração destes e, posteriormente, a elaboração dos espectros de absorção das espécies químicas. São analisados meios uniformes (isotérmicos e homogêneos) primeiramente, evidenciando que o modelo original (CW) não apresenta problemas nesses casos. Posteriormente, meios não uniformes (nãoisotérmicos e homogêneos; ou não-isotérmicos e não homogêneos) são avaliados. Os resultados para ambos os modelos, CW e CWM, são comparados com a solução linha-por-linha (LBL). Casos com diferentes perfis de temperatura e concentração da espécie química são considerados. A solução a partir da nova metodologia se apresentou mais dispendiosa computacionalmente devido ao acréscimo de novos passos iterativos. Os resultados mostram que o CWM leva a resultados mais precisos para o fluxo de calor radiante, satisfazendo ao balanço de energia com um desvio que se deve principalmente à discretização espacial da equação de transporte radiante. / This work presents a modification of the cumulative wavenumber (CW) method to determine the radiative transfer in non-uniform participating gases to enforce the radiative energy balance to be satisfied. In particular, the CW model leads to results for the divergent radiative heat flux (or radiative volumetric heat source) that proved accurate in comparison to the benchmark solution ( LBL integration) for non-isothermal medium. However, as will be shown in this work, one important assumption of the method prevents it of satisfying the radiative energy balance when a non uniform medium is considered. As consequence, the radiative heat flux can present considerable deviation of the correct solution. The modified cumulative wavenumber (CWM) model was developed to keep the same value of the radiative volumetric heat source, but also to satisfy the radiative energy balance. The CWM model is applied together with the discrete ordinates method to solve the radiation heat transfer in a one-dimentional slab containing a uniform/non-uniform layer of typical gases from the methane or fuel oil combustion. The walls are black for the thermal radiation. The HITRAN and HITEMP are used to extract the spectral lines information that is required for modeling. The database are briefly presented as well as the procedure that is used to extract the data and the spectra drawing. Firstly, uniform media (isothermal and homogeneous) are analyzed proving that the original model (CW) presents accurate results in those cases. Next, non-uniform media (non-isothermal and homogeneous; or non-isothermal and non-homogeneous) are taken into account. The results of both the CW and the CWM modeling are compared to the benchmark line-by-line (LBL) integration. Different temperature and concentration profiles are considered. The solution obtained from the new model presents an increase in the computational time due to the insertion of new iterative loops. The results show that the CWM leads to accurate estimation of both the radiative heat flux and volumetric heat source, satisfying the radiative energy balance with an error that is mostly due to the spatial discretization of the radiative transfer equation.

O trabalho apresenta uma modificação do modelo do número de onda cumulativo (CW – Cumulative Wavenumber) para atender ao balanço de energia na determinação da transferência radiante de gases participantes não uniformes. De um modo geral, o modelo CW fornece resultados que se mostram precisos para o divergente do fluxo de calor radiante (ou taxa volumétrica de geração de energia) quando comparados à solução benchmark (integração linhapor- linha). No entanto, conforme é apresentado nesta pesquisa, uma hipótese importante feita pelo modelo o torna não conservativo e, não garantindo o balanço de energia radiante quando um meio não uniforme é considerado. Como conseqüência, o fluxo de calor radiante pode apresentar desvios consideráveis em relação à solução exata. O modelo do número de onda cumulativo modificado (CWM) foi desenvolvido de forma a manter o mesmo valor para a taxa volumétrica de geração de energia, mas satisfazendo ao balanço de energia radiante. O modelo CWM é aplicado juntamente com o método das ordenadas discretas para resolver a transferência de calor radiante em sistemas unidimensionais formados por placas planas infinitas e paralelas, contendo camadas uniformes e não uniformes de gases típicos da combustão de metano ou óleo combustível. As paredes são negras para a radiação térmica. A aquisição de dados das linhas de absorção espectrais é feita pela utilização dos bancos de dados HITRAN e HITEMP. São apresentados os bancos de dados, bem como a forma em que é efetuada a extração destes e, posteriormente, a elaboração dos espectros de absorção das espécies químicas. São analisados meios uniformes (isotérmicos e homogêneos) primeiramente, evidenciando que o modelo original (CW) não apresenta problemas nesses casos. Posteriormente, meios não uniformes (nãoisotérmicos e homogêneos; ou não-isotérmicos e não homogêneos) são avaliados. Os resultados para ambos os modelos, CW e CWM, são comparados com a solução linha-por-linha (LBL). Casos com diferentes perfis de temperatura e concentração da espécie química são considerados. A solução a partir da nova metodologia se apresentou mais dispendiosa computacionalmente devido ao acréscimo de novos passos iterativos. Os resultados mostram que o CWM leva a resultados mais precisos para o fluxo de calor radiante, satisfazendo ao balanço de energia com um desvio que se deve principalmente à discretização espacial da equação de transporte radiante. / This work presents a modification of the cumulative wavenumber (CW) method to determine the radiative transfer in non-uniform participating gases to enforce the radiative energy balance to be satisfied. In particular, the CW model leads to results for the divergent radiative heat flux (or radiative volumetric heat source) that proved accurate in comparison to the benchmark solution ( LBL integration) for non-isothermal medium. However, as will be shown in this work, one important assumption of the method prevents it of satisfying the radiative energy balance when a non uniform medium is considered. As consequence, the radiative heat flux can present considerable deviation of the correct solution. The modified cumulative wavenumber (CWM) model was developed to keep the same value of the radiative volumetric heat source, but also to satisfy the radiative energy balance. The CWM model is applied together with the discrete ordinates method to solve the radiation heat transfer in a one-dimentional slab containing a uniform/non-uniform layer of typical gases from the methane or fuel oil combustion. The walls are black for the thermal radiation. The HITRAN and HITEMP are used to extract the spectral lines information that is required for modeling. The database are briefly presented as well as the procedure that is used to extract the data and the spectra drawing. Firstly, uniform media (isothermal and homogeneous) are analyzed proving that the original model (CW) presents accurate results in those cases. Next, non-uniform media (non-isothermal and homogeneous; or non-isothermal and non-homogeneous) are taken into account. The results of both the CW and the CWM modeling are compared to the benchmark line-by-line (LBL) integration. Different temperature and concentration profiles are considered. The solution obtained from the new model presents an increase in the computational time due to the insertion of new iterative loops. The results show that the CWM leads to accurate estimation of both the radiative heat flux and volumetric heat source, satisfying the radiative energy balance with an error that is mostly due to the spatial discretization of the radiative transfer equation.

O trabalho apresenta uma modificação do modelo do número de onda cumulativo (CW – Cumulative Wavenumber) para atender ao balanço de energia na determinação da transferência radiante de gases participantes não uniformes. De um modo geral, o modelo CW fornece resultados que se mostram precisos para o divergente do fluxo de calor radiante (ou taxa volumétrica de geração de energia) quando comparados à solução benchmark (integração linhapor- linha). No entanto, conforme é apresentado nesta pesquisa, uma hipótese importante feita pelo modelo o torna não conservativo e, não garantindo o balanço de energia radiante quando um meio não uniforme é considerado. Como conseqüência, o fluxo de calor radiante pode apresentar desvios consideráveis em relação à solução exata. O modelo do número de onda cumulativo modificado (CWM) foi desenvolvido de forma a manter o mesmo valor para a taxa volumétrica de geração de energia, mas satisfazendo ao balanço de energia radiante. O modelo CWM é aplicado juntamente com o método das ordenadas discretas para resolver a transferência de calor radiante em sistemas unidimensionais formados por placas planas infinitas e paralelas, contendo camadas uniformes e não uniformes de gases típicos da combustão de metano ou óleo combustível. As paredes são negras para a radiação térmica. A aquisição de dados das linhas de absorção espectrais é feita pela utilização dos bancos de dados HITRAN e HITEMP. São apresentados os bancos de dados, bem como a forma em que é efetuada a extração destes e, posteriormente, a elaboração dos espectros de absorção das espécies químicas. São analisados meios uniformes (isotérmicos e homogêneos) primeiramente, evidenciando que o modelo original (CW) não apresenta problemas nesses casos. Posteriormente, meios não uniformes (nãoisotérmicos e homogêneos; ou não-isotérmicos e não homogêneos) são avaliados. Os resultados para ambos os modelos, CW e CWM, são comparados com a solução linha-por-linha (LBL). Casos com diferentes perfis de temperatura e concentração da espécie química são considerados. A solução a partir da nova metodologia se apresentou mais dispendiosa computacionalmente devido ao acréscimo de novos passos iterativos. Os resultados mostram que o CWM leva a resultados mais precisos para o fluxo de calor radiante, satisfazendo ao balanço de energia com um desvio que se deve principalmente à discretização espacial da equação de transporte radiante. / This work presents a modification of the cumulative wavenumber (CW) method to determine the radiative transfer in non-uniform participating gases to enforce the radiative energy balance to be satisfied. In particular, the CW model leads to results for the divergent radiative heat flux (or radiative volumetric heat source) that proved accurate in comparison to the benchmark solution ( LBL integration) for non-isothermal medium. However, as will be shown in this work, one important assumption of the method prevents it of satisfying the radiative energy balance when a non uniform medium is considered. As consequence, the radiative heat flux can present considerable deviation of the correct solution. The modified cumulative wavenumber (CWM) model was developed to keep the same value of the radiative volumetric heat source, but also to satisfy the radiative energy balance. The CWM model is applied together with the discrete ordinates method to solve the radiation heat transfer in a one-dimentional slab containing a uniform/non-uniform layer of typical gases from the methane or fuel oil combustion. The walls are black for the thermal radiation. The HITRAN and HITEMP are used to extract the spectral lines information that is required for modeling. The database are briefly presented as well as the procedure that is used to extract the data and the spectra drawing. Firstly, uniform media (isothermal and homogeneous) are analyzed proving that the original model (CW) presents accurate results in those cases. Next, non-uniform media (non-isothermal and homogeneous; or non-isothermal and non-homogeneous) are taken into account. The results of both the CW and the CWM modeling are compared to the benchmark line-by-line (LBL) integration. Different temperature and concentration profiles are considered. The solution obtained from the new model presents an increase in the computational time due to the insertion of new iterative loops. The results show that the CWM leads to accurate estimation of both the radiative heat flux and volumetric heat source, satisfying the radiative energy balance with an error that is mostly due to the spatial discretization of the radiative transfer equation.

[pt] A energia eólica é hoje uma das mais promissoras fontes de energia do mundo por ser limpa e abundante. O estudo de fenômenos que estão relacionados com alterações na circulação atmosférica, como o El Niño, são de extrema importância pela possibilidade de afetarem a geração eólica. A fim de explorar o possível efeito de tais fenômenos nos ventos da região Nordeste do Brasil, a qual possui a maior capacidade eólica instalada, é realizada uma análise estatística para a quantificação desse efeito através dos modelos Generalized Autoregressive Score (GAS). Estes permitem a modelagem de séries temporais para diferentes distribuições de probabilidade. Nesse estudo, o GAS é aplicado às séries de velocidade do vento a partir das distribuições Weibull, Gama e Beta. Os resultados mostram que a distribuição Beta produz o melhor ajuste quando se considera o efeito do El Niño através da Temperatura da Superfície do Mar (TSM), na velocidade do vento de uma usina eólica no Nordeste do Brasil. Em particular estes mostram que, o maior impacto do El Niño ocorre nos meses de abril e maio, onde para um aumento de 1 grau Celcius da TSM é observada uma variação de aproximadamente 1.3 por cento na velocidade do vento. / [en] Wind energy is now one of the most promising energy sources of the world being both clean and abundant. The study of phenomena that are related to changes in atmospheric circulation, such as El NiÑo, are extremely important for its ability to affect wind generation. In order to explore the possible effect of such phenomena in the winds of the Northeast region of Brazil which has the largest installed wind capacity, a statistical analysis is performed to quantify this effect through the Generalized Autoregressive Score (GAS) models. These allows the modeling of time series for different probability distributions. In this study, GAS models are applied to the wind speed series using Weibull, Gamma and Beta distributions. The results shows that the Beta distribution produces the best fit when considering the effect of El Niño through Sea Surface Temperature (SST) in the wind speed of an wind farm in the Northeast of Brazil. In particular these show that the greatest impact of the El Niño occurs in the months of April and May, where for an increase of 1 degrees Celcius of SST is observed a variation of approximately 1.3 per cent in wind speed.

[pt] Neste trabalho usamos o arcabouço dos modelos GAS para gerar previsões conjuntas de fator de capacidade eólico, pertencentes a diferentes usinas localizadas em áreas geográficas distintas. Esses cenários são insumos para gerar uma distribuição de fluxo de caixa associada a um portfólio de contratos atrelados aos parques eólicos em questão. Inicialmente modelamos as densidades marginais via um modelo GAS, supondo densidade Beta. De maneira a capturar a estrutura de dependência entre esses fatores de capacidade, usamos uma cópula t-Student com a matriz de correlação também sendo atualizada via mecanismo GAS. Uma das contribuições importantes desse trabalho para o setor elétrico está na geração de cenários conjuntos apenas em um passo, evitando a necessidade de modelar variáveis transformadas e posteriormente transforma-las para retornar às suas respectivas escalas originais. Assim como é feito no caso supondo normalidade para as marginais. Como é sabido, exponenciar valores simulados a partir de uma densidade normal pode gerar resultados equivocados para fatores de capacidade eólico, e por propagação, isso pode afetar severamente as medidas de risco que são obtidas a partir da distribuição simulada de fluxo de caixa associada com o portfolio das usinas eólicas. Nossos resultados mostram que quando a dependência é levada em consideração, os fluxos de caixa tendem a ser maiores do que quando ignora-se a dependência. / [en] In this work we use the framework of GAS models to generate joint forecasts for capacity factors of several wind plants belonging to different geographical areas. Such scenarios are then used as input to raise the distribution of cash flows associated with a portfolio of contracts attached to these wind plants. We first model the marginal density of each capacity factor using a GAS model with Beta density. In order to capture the observed dependence among these capacity factors, we use a copula t- Student with correlation matrix evolving through a GAS mechanism. One of the important contributions of our framework is that generation of scenarios is accomplished in just one step, avoiding the need of transforming back variables to its original scale, as it is the case under a Gaussian assumption for the marginals. As it is known, exponentiation of simulated Gaussian values can result in unrealistic sampling paths for the wind capacity factor, and by propagation, this can badly a ect the risk measures obtained from the simulated distribution of the cash flows associated with a particular portfolio of wind plants. Our results shows that when taking into account dependence the cash flows are higher than when ignoring dependence.

[pt] Esta tese apresente duas contribuições para a modelagem e previsão de sinistros já ocorridos e ainda não avisados (Incurred But Not Reported – IBNR), quando organizados numa estrutura de dados conhecida como triângulo de run-off. Ambas as contribuições são baseadas em arcabouços gerais para a construção de modelos para séries temporais com coeficientes variantes no tempo. Em nossa primeira contribuição desenvolvemos a extensão multivariada do modelo em espaço de estado proposto por Atherino em 2008. A partir dessa extensão é possível modelas simultaneamente um ou mais triângulos de run-off associados às diversas coberturas de uma seguradora, levando-se em consideração a dependência entre os distintos triângulos, capturada pela estrutura da matriz variância – covariância do modelo SUTSE e a dependência entre as células de cada triângulo de run-off, capturada pelas componentes de nível e de periodicidade, de acordo com a proposta de atherino et al. (2010). Em nossa segunda contribuição desenvolvemos um arcabouço geral para a modelagem univariada de triângulos de run-off a partir da estruturas dos modelos GAS (Generalized Autoregressive Score) desenvolvidos por Creal at al. (2013). Esse arcabouço, bastante flexível, permite a escolha de qualquer distribuição para as entradas do triângulo run-off, considerando que os seus parâmetros variem ao longo período de origem ou de desenvolvimento. Em particular consideramos as distribuições gama e log-normal. Nossos resultados foram comparados com os obtidos através do método chain ladder (Mack, 1993), utilizado como benchmark na indústria de seguros. O teste Diebold e Mariano (1995) evidenciou que os modelos propostos geram melhores previsões, comparadas as previsões do método chain ladder. / [en] This thesis presents two contributions to the modeling and prediction of a type of claims in the insurance industry known as IBNR (Incurred But Not Reported) when these are organized in a data structure known as the run-off triangle. Both contributions are based on general frameworks for building models for time series with time varying coefficients. In our first contribution we developed the multivariate extension of the state space model proposed by Atherino in 2008. From this extension it is possible to model simultaneously one or more run-off triangles associated with different coverages from an insurer, taking into account the dependence between different triangles, captured by the structure of the variance-covariance matrix of the SUTSE model, while the dependence between the cells of each triangle run-off is captured by the components of level and periodicity according to the model proposed by Atherino (2008). In our second contribution we developed a general framework for univariate modeling of run-off triangles using the structure of GAS models (Generalized Autoregressive Score) developed by Creal et al. (2013). This framework, very flexible, allows one to choose any distribution to the inputs of the triangle run-off, considering that its parameters can vary over the period of origin or period of development. In particular we have considered both gamma and lognormal distributions. Our results were compared with those obtained by the chain ladder method (Mack, 1993) used as a benchmark in the insurance industry. The Diebold and Mariano test (1995) showed that the proposed models produced better predictions compared to the predictions of the chain ladder method.

Ultra-relativisitic heavy ion collisions produce quark gluon plasma-a hot and dense soup of deconfined quarks and gluons akin to the early universe. We study two models in the context of these collisions namely, Polyakov Quark Meson Model (PQM) and Hadron Resonance Gas Model (HRGM).The PQM Model provides us with a simple and intuitive understanding of the QCD equation of state and thermodynamics at non zero temperature and baryon density while the HRGM is the principle model to analyse the hadron yields measured in these experiments across the entire range of beam energies. We study the effect of including the commonly neglected fermionic vacuum fluctuations to the (2+1) flavor PQM model. The conventional PQM model suffers from a rapid phase transition contrary to what is found through lattice simulations. Addition of the vacuum term tames the rapid transition and significantly improves the model’s agreement to lattice data. We further investigate the role of the vacuum term on the phase diagram. The smoothening effect of the vacuum term persists even at non zero . Depending on the value of the mass of the sigma meson, including the vacuum term results in either pushing the critical end point into higher values of the chemical potential or excluding the possibility of a critical end point altogether. We compute the fluctuations(correlations) of conserved charges up to sixth(fourth) order. Comparison is made with lattice data wherever available and overall good qualitative agreement is found, more so for the case of the normalised susceptibilities. The model predictions for the ratio of susceptibilities approach to that of an ideal gas of hadrons as in HRGM at low temperatures while at high temperature the values are close to that of an ideal gas of massless quarks. We examine the stability of HRGMs by extending them to take care of undiscovered resonances through the Hagedorn formula. We find that the influence of unknown resonances on thermodynamics is large but bounded. We model the decays of resonances and investigate the ratios of particle yields in heavy-ion collisions. We find that extending these models do not have much effect on hydrodynamics but the hadron yield ratios show better agreement with experiment. In principle HRGMs are internally consistent up to a temperature higher than the cross over temperature in QCD; but by examining quark number susceptibilities we find that their region of applicability seems to end even below the QCD cross over.

[pt] Demanda intermitente é definida por períodos de vendas nulas intercaladas com vendas positivas e de quantidade altamente variável. A maior parte das unidades de manutenção de estoque (stock keeping units, em inglês) ao nível loja pode ser caracterizada como contendo demanda desse tipo. Assim, modelos acurados para prever séries com demanda intermitente trazem grandes impactos em relação à gestão de estoque. Nesta dissertação nós propomos o uso do framework GAS com as distribuições adequadas para dados de contagem, além de suas versões com excesso de zeros, e aplicamos os modelos derivados a dados reais obtidos com uma grande rede varejista brasileira. Nós demonstramos que os modelos com excesso de zeros propostos são estimados de forma consistente por máxima verossimilhança e a distribuição dos estimadores é assintóticamente normal. A performance dos modelos propostos é comparada com benchmarks adequados das literaturas de séries temporais para dados de contagem e previsão de demanda intermitente. A avaliação das previsões é feita com base tanto na precisão da distribuição preditiva quanto na precisão das previsões pontuais. Nossos resultados mostram que os modelos propostos, em especial o modelo derivado sob distribuição hurdle Poisson, performam melhor do que os benchmarks analisados. / [en] Intermittent demand is defined by periods of zero sales interleaved with positive sales with highly variable quantities. Most stock keeping units at the store level can be characterized as containing such demand. Thus, accurate models for predicting series with intermittent demand have major impacts in relation to inventory management. In this dissertation we propose the use of the GAS framework with the appropriate distributions for count data, in addition to their versions with excess of zeroes, and apply the derived models to real data obtained from a large Brazilian retail chain. We demonstrate that the proposed models with excess of zeros are consistently estimated via maximum likelihood and the distribution of the estimator is asymptotically normal. The performance of the proposed models is compared to adequate benchmarks from the time series literature for count data and intermittent demand forecast. Forecasting is evaluated based on the accuracy of both the entire predictive distribution and point forecasts. Our results show that the proposed models, specially the one derived from hurdle Poisson distribution, perform better than the analyzed benchmarks.

[pt] O objetivo desta tese é apresentar novas alternativas para modelagem de variáveis aleatórias no setor de seguros, utilizando o arcabouço dos modelos orientados por score com parâmetros variantes no tempo. No primeiro artigo, propomos um modelo dinâmico para a distribuição do sinistro agregado, que corresponde à soma aleatória dos valores de sinistros (severidade) em determinado período de tempo. A obtenção da distribuição do sinistro agregado é um problema clássico na teoria do risco e fundamental para precificação de seguros, cálculo de provisões e de probabilidade de ruína. No entanto, a obtenção da expressão analítica para essa distribuição de probabilidade é uma tarefa difícil. Neste trabalho, especificamos distribuições não-Gaussianas, tanto para o número de sinistros como para severidade, sob o arcabouço GAS (Generalized Autoregressive Score), e, por meio do uso da Transformada Rápida de Fourier obtemos, numericamente, a distribuição do sinistro agregado. O segundo artigo trata da incorporação do efeito de variáveis macroeconômicas na modelagem de variáveis relevantes no setor de seguros, em linha com requisito internacional de avaliação de provisões de forma consistente com mercado, a qual leva em consideração as informações disponíveis nos mercados financeiros e de capital relevantes, utilizando metodologias e parâmetros consistentes com esses mercados. Modelamos uma série bivariada de número de sinistros (duas linhas de negócios) de seguros financeiros com modelos autorregressivos e utilizamos cópulas para modelar a estrutura de dependência das séries temporais condicionado aos modelos ajustados nas marginais. Com esta abordagem, é possível simular números de sinistros futuros de mais de uma carteira, podendo esse resultado ser utilizado em uma avaliação consistente de provisões e da saúde financeira da seguradora. / [en] The objective of this thesis is to present new alternatives for modeling random variables in the insurance industry, using the framework of the score driven models with time-varying parameters. In the first paper, we propose a dynamic model for the aggregate claims distribution, which corresponds to a random sum of claims severity in a certain period of time. Obtaining the aggregate claims distribution is a classic problem in the Risk Theory and fundamental for premium estimation, measurement of obligations and ruin probability valuation. However, obtaining the analytic expression for this probability distribution is a hard task. In this work, we specify nonGaussian distributions for both the number of claims and for the claims severity, under the GAS framework, and, through the use of the fast Fourier transform, we obtain, numerically, the aggregate claims distribution. The second paper deals with the incorporation of macroeconomic variables on the modeling of relevant variables in the insurance sector, in line with the international requirements for market consistent valuation of insurance liabilities, which means that one should take into account the available information in relevant financial and capital markets, using methodologies and parameters consistent with these markets. We model a bivariate time series (two lines of business) of financial insurance with autoregressive models and use copulas models to consider the dependency structure of the time series conditioned to the fitted models for the marginals. Within this approach, it is possible to simulate the numbers of claims from more than one portfolio, and this result can be used in a consistent valuation of liabilities and of the financial health of an insurer.