El método de las ecuaciones integrales singulares (B.I.E.M.)

El Método de Ecuaciones Integrales Singulares es una alternativa a los métodos dominantes, como el conocido Método de los Elementos Finitos. La idea esencial es la combinación de las clásicas relaciones recíprocas con la filosofía de la discretización del Método de los Elementos Finitos.

Los métodos proyectivos en la mecánica de los medios contínuos

El Método de las Ecuaciones Integrales es una potente alternativa a los Métodos de Dominio tales como el Método de los Elementos Finitos. La idea ensencial es la combinación de la clásica relación de la reciprocidad con la filosofía de la discretización del F.E.M. La aplicación a algunos problemas reales ha demostrado que en ciertos casos el B.I.E.M. es preferiole al F.E.M. y ello es especialmente así cuando los problemas a tratar son tridimensionales y con geometría complicada. En esta ocasión se analizan comparativamente algunos aspectos matemáticos del procedimiento = Boundary integral equation method (B.I.E.M.)is a powerful alternative to the domain methods, as the well know Finite Element Method (F .E.M.) The esential idea, are the combination of the classical reciprocity re!ations with the discretization phylosophy of F.E.M. The reduction in dimension of the domain to be discretized, the easy treatment of infinite domains and the high accuracy of the results are the main adventages of B.I.E.M. Between the drawacks the nonsymetry and non sparseness of the matrices to be treated are worth remembering. Application to several real problems has shown that in certain cases B.I.E.M. is better than F.E.M. and this is specially true when tridimensional problems of complicated geometries have to be treated. Active research is in progress of its extensión to non linear and time dependent problems.

B.I.E.M. y estudio de singularidades en problemas de potencial

El objeto del presente artículo es el estudio de singularidades en problemas de Potencial mediante el uso del Método de las Ecuaciones Integrales sobre el contorno del dominio en estudio. Frente a soluciones basadas en la mejora de la discretización, análisis asintótico o introducción de funciones de forma que representen mejor la evolución de la función, una nueva hipótesis es presentada: el término responsable de la singularidad es incluido en la integral sobre el contorno de la función auxiliar. Los resultados obtenidos mejoran los de soluciones anteriores simplificando también el tiempo de cálculo = The subject of this paper is the modelling of singularities in potential problems, using the Boundary Integral Equation Method. As a logical alternative to classical methods (discretization refinement, asymptotic analysis, high order interpolatory functions) a new hypothesis is presented: the singularity responsible term is included in the interpolatory shape function. As shown by several exemples results are splendid and computer time radically shortened.

Contribución al análisis eficiente y a la mejora de prestaciones de antenas reflectarray

El diseño de una antena reflectarray bajo la aproximación de periodicidad local requiere la determinación de la matriz de scattering de estructuras multicapa con metalizaciones periódicas para un gran número de geometrías diferentes. Por lo tanto, a la hora de diseñar antenas reflectarray en tiempos de CPU razonables, se necesitan herramientas númericas rápidas y precisas para el análisis de las estructuras periódicas multicapa. En esta tesis se aplica la versión Galerkin del Método de los Momentos (MDM) en el dominio espectral al análisis de las estructuras periódicas multicapa necesarias para el diseño de antenas reflectarray basadas en parches apilados o en dipolos paralelos coplanares. Desgraciadamente, la aplicación de este método numérico involucra el cálculo de series dobles infinitas, y mientras que algunas series convergen muy rápidamente, otras lo hacen muy lentamente. Para aliviar este problema, en esta tesis se propone un novedoso MDM espectral-espacial para el análisis de las estructuras periódicas multicapa, en el cual las series rápidamente convergente se calculan en el dominio espectral, y las series lentamente convergentes se calculan en el dominio espacial mediante una versión mejorada de la formulación de ecuaciones integrales de potenciales mixtos (EIPM) del MDM. Esta versión mejorada se basa en la interpolación eficiente de las funciones de Green multicapa periódicas, y en el cálculo eficiente de las integrales singulares que conducen a los elementos de la matriz del MDM. El novedoso método híbrido espectral-espacial y el tradicional MDM en el dominio espectral se han comparado en el caso de los elementos reflectarray basado en parches apilados. Las simulaciones numéricas han demostrado que el tiempo de CPU requerido por el MDM híbrido es alrededor de unas 60 veces más rápido que el requerido por el tradicional MDM en el dominio espectral para una precisión de dos cifras significativas. El uso combinado de elementos reflectarray con parches apilados y técnicas de optimización de banda ancha ha hecho posible diseñar antenas reflectarray de transmisiónrecepción (Tx-Rx) y polarización dual para aplicaciones de espacio con requisitos muy restrictivos. Desgraciadamente, el nivel de aislamiento entre las polarizaciones ortogonales en antenas DBS (típicamente 30 dB) es demasiado exigente para ser conseguido con las antenas basadas en parches apilados. Además, el uso de elementos reflectarray con parches apilados conlleva procesos de fabricación complejos y costosos. En esta tesis se investigan varias configuraciones de elementos reflectarray basadas en conjuntos de dipolos paralelos con el fin de superar los inconvenientes que presenta el elemento basado en parches apilados. Primeramente, se propone un elemento consistente en dos conjuntos apilados ortogonales de tres dipolos paralelos para aplicaciones de polarización dual. Se ha diseñado, fabricado y medido una antena basada en este elemento, y los resultados obtenidos para la antena indican que tiene unas altas prestaciones en términos de ancho de banda, pérdidas, eficiencia y discriminación contrapolar, además de requerir un proceso de fabricación mucho más sencillo que el de las antenas basadas en tres parches apilados. Desgraciadamente, el elemento basado en dos conjuntos ortogonales de tres dipolos paralelos no proporciona suficientes grados de libertad para diseñar antenas reflectarray de transmisión-recepción (Tx-Rx) de polarización dual para aplicaciones de espacio por medio de técnicas de optimización de banda ancha. Por este motivo, en la tesis se propone un nuevo elemento reflectarray que proporciona los grados de libertad suficientes para cada polarización. El nuevo elemento consiste en dos conjuntos ortogonales de cuatro dipolos paralelos. Cada conjunto contiene tres dipolos coplanares y un dipolo apilado. Para poder acomodar los dos conjuntos de dipolos en una sola celda de la antena reflectarray, el conjunto de dipolos de una polarización está desplazado medio período con respecto al conjunto de dipolos de la otra polarización. Este hecho permite usar solamente dos niveles de metalización para cada elemento de la antena, lo cual simplifica el proceso de fabricación como en el caso del elemento basados en dos conjuntos de tres dipolos paralelos coplanares. Una antena de doble polarización y doble banda (Tx-Rx) basada en el nuevo elemento ha sido diseñada, fabricada y medida. La antena muestra muy buenas presentaciones en las dos bandas de frecuencia con muy bajos niveles de polarización cruzada. Simulaciones numéricas presentadas en la tesis muestran que estos bajos de niveles de polarización cruzada se pueden reducir todavía más si se llevan a cabo pequeñas rotaciones de los dos conjuntos de dipolos asociados a cada polarización. ABSTRACT The design of a reflectarray antenna under the local periodicity assumption requires the determination of the scattering matrix of a multilayered structure with periodic metallizations for quite a large number of different geometries. Therefore, in order to design reflectarray antennas within reasonable CPU times, fast and accurate numerical tools for the analysis of the periodic multilayered structures are required. In this thesis the Galerkin’s version of the Method of Moments (MoM) in the spectral domain is applied to the analysis of the periodic multilayered structures involved in the design of reflectarray antennas made of either stacked patches or coplanar parallel dipoles. Unfortunately, this numerical approach involves the computation of double infinite summations, and whereas some of these summations converge very fast, some others converge very slowly. In order to alleviate this problem, in the thesis a novel hybrid MoM spectral-spatial domain approach is proposed for the analysis of the periodic multilayered structures. In the novel approach, whereas the fast convergent summations are computed in the spectral domain, the slowly convergent summations are computed by means of an enhanced Mixed Potential Integral Equation (MPIE) formulation of the MoM in the spatial domain. This enhanced formulation is based on the efficient interpolation of the multilayered periodic Green’s functions, and on the efficient computation of the singular integrals leading to the MoM matrix entries. The novel hybrid spectral-spatial MoM code and the standard spectral domain MoM code have both been compared in the case of reflectarray elements based on multilayered stacked patches. Numerical simulations have shown that the CPU time required by the hybrid MoM is around 60 times smaller than that required by the standard spectral MoM for an accuracy of two significant figures. The combined use of reflectarray elements based on stacked patches and wideband optimization techniques has made it possible to design dual polarization transmit-receive (Tx-Rx) reflectarrays for space applications with stringent requirements. Unfortunately, the required level of isolation between orthogonal polarizations in DBS antennas (typically 30 dB) is hard to achieve with the configuration of stacked patches. Moreover, the use of reflectarrays based on stacked patches leads to a complex and expensive manufacturing process. In this thesis, we investigate several configurations of reflectarray elements based on sets of parallel dipoles that try to overcome the drawbacks introduced by the element based on stacked patches. First, an element based on two stacked orthogonal sets of three coplanar parallel dipoles is proposed for dual polarization applications. An antenna made of this element has been designed, manufactured and measured, and the results obtained show that the antenna presents a high performance in terms of bandwidth, losses, efficiency and cross-polarization discrimination, while the manufacturing process is cheaper and simpler than that of the antennas made of stacked patches. Unfortunately, the element based on two sets of three coplanar parallel dipoles does not provide enough degrees of freedom to design dual-polarization transmit-receive (Tx-Rx) reflectarray antennas for space applications by means of wideband optimization techniques. For this reason, in the thesis a new reflectarray element is proposed which does provide enough degrees of freedom for each polarization. This new element consists of two orthogonal sets of four parallel dipoles, each set containing three coplanar dipoles and one stacked dipole. In order to accommodate the two sets of dipoles in each reflectarray cell, the set of dipoles for one polarization is shifted half a period from the set of dipoles for the other polarization. This also makes it possible to use only two levels of metallization for the reflectarray element, which simplifies the manufacturing process as in the case of the reflectarray element based on two sets of three parallel dipoles. A dual polarization dual-band (Tx-Rx) reflectarray antenna based on the new element has been designed, manufactured and measured. The antenna shows a very good performance in both Tx and Rx frequency bands with very low levels of cross-polarization. Numerical simulations carried out in the thesis have shown that the low levels of cross-polarization can be even made smaller by means of small rotations of the two sets of dipoles associated to each polarization.

Productividad y Factores Humanos. Un Modelo con Ecuaciones Estructurales

Las personas se asocian, colaboran e interactúan para el logro de sus metas y los de la organización. Son las personas las que le otorgan ‘vida’ a las organizaciones. Asimismo, este hecho asociativo genera procesos psicológicos y psicosociales que afectan a los individuos, a los grupos y a los resultados de la organización. El conocimiento de estos procesos es de suma importancia para quienes dirigen las organizaciones, ya que su adecuada gestión permitirá obtener los resultados deseados. A fin de establecer los factores humanos que inciden en la productividad, se realizó un análisis factorial donde se obtuvieron cuatro componentes, lo que permitió clasificarlos en cuatro factores: individual, grupal, organizacional y de resultados. De acuerdo con la fundamentación teórica se encontró que las variables determinaron cuatro constructos latentes subyacentes, por lo cual se planteó un modelo conformado por tres variables latentes exógenas denominadas ‘ factores individuales’, ‘ factores grupales’ y ‘ factores organizacionales’, que se relacionan con una variable latente endógena denominada ‘productividad’. Mediante la técnica de modelos de ecuaciones estructurales (SEM) se constataron las relaciones entre estas variables latentes y sus indicadores asociados. Las relaciones inicialmente propuestas entre las variables latentes fueron corroboradas por los ajustes del modelo SEM

Ecuaciones de cubicación para el pino silvestre (Pinus sylvestris L.) en el monte Aguas Vertientes N.138 del C.U.P. en el municipio de El Espinar (Segovia).

La finalidad de este trabajo es proponer unas nuevas ecuaciones de cubicación, ya que los métodos utilizados por la Administración podrían estar anticuados. Con estas ecuaciones actualizadas se intentan paliar las diferencias entre Administración y aserradero. Desde hace unos 25 años la Administración de Segovia viene utilizando unos valores modulares para estimar el volumen en pié de los árboles del monte "Aguas Vertientes" Estos valores modulares son tablas en las que se da un volumen teórico para cada pié en cada uno de los rodales y según clases diamétricas de 5 cm. Estas tablas fueron calculadas hace tiempo, y, por lo tanto, sería recomendable realizar unas nuevas tarifas de cubicación. Los objetivos del presente trabajo, son, por tanto, construir una nueva tarifa de cubicación de una entrada V=f(dn) y una tabla de cubicación con dos entradas V=f(dn, H)

Evaluación de integrales singulares y cuasisingulares en el método de los elementos de contorno bidimensional mediante transformaciones no lineales

Se presentan en este trabajo distintos niétodos, basados en transformaciones de coordenadas no-lineales, para la evaluación de las integrales singulares y cuasisingulares que aparecen en el Método Directo de los Elementos de Contorno. Se detecta un error inherente a algunas de las transformaciones propuestas y finalmente se sugieren dos nuevas transformaciones que mejoran las actualmente disponibles.

Problemas de matemáticas

Este libro, Problemas de Matemáticas, junto con otros dos, Problemas de Geometría y Problemas de Geometría Analítica y Diferencial, están dedicados a la presentación y resolución de problemas que se planteaban hace unas décadas, en la preparación para ingreso en las carreras de ingeniería técnica superior. Incluye 1578 problemas, de los que 848 se refieren al Álgebra (operaciones algebraicas, divisibilidad, combinatoria, determinantes, ecuaciones e inecuaciones, fracciones continuas, números complejos, límites, sucesiones y series, y algunos sobre vectores y mecánica), 175 a la Trigonometría (plana y esférica), 282 al Cálculo diferencial (funciones de una variable, y de dos o más variables), 246 al Cálculo integral (integrales, integrales definidas, integrales en el campo de dos o más variables y ecuaciones diferenciales) y 27 a la Estadística. Esta tercera edición de Problemas de Matemáticas tiene por objeto su puesta a disposición de la Escuela de Ingenieros de Minas de la Universidad Politécnica de Madrid.

Nueva formulación para la obtención de la energía cinética de la lluvia. Análisis comparativo con ecuaciones empíricas

El presente trabajo describe el estudio de la caracterización cuantitativa de una propiedad de la lluvia muy importante, la energía cinética, y es importante porque es la fuente de entrada de energía en los procesos de erosión hídrica y por tanto, de la pérdida de suelo que sufren nuestros campos.Este trabajo presenta un estudio de la relación existente entre la energía cinética de la lluvia y la intensidad de la lluvia con datos de lluvia registrados de forma automática con el disdrómetro de JOSS y WALDVOGEL (1967). Se ha obtenido una ecuación de regresión entre la energía cinética y la intensi- dad pero con algunas peculiaridades, ya que es una ecuación exponencial que introduce un valor máximo de energía cinética. Esta ecuación exponencial ha proporcionado resultados muy satisfacto- nos, presenta errores relativos bajos y además tiene un significado práctico y físico.Además se ha realizado un análisis comparativo entre los resultados obtenidos con esta ecuación exponencial y ecuaciones empíricas tradicionales en el cálculo de la energía cinética. El criterio de selección de las ecuaciones a comparar fue en primer lugar escoger aquellas ecuaciones que fuesen habituales en los estudios de erosión y también aquellas ecuaciones que por la localización donde fueron obtenidas pudiesen ser interesante para su comparación.Los resultados obtenidos de la comparación conducen a pensar que la relación entre la energía cinética y la intensidad no estaría influenciado por condiciones de localización. Se concluye tambien que el uso de la ecuación exponencial obtenida en este estudio resolvería algunos de los problemas que conlleva el uso de algunas de la ecuaciones empíricas.

Un nuevo procedimiento numérico para la evaluación de integrales singulares de Valor Principal de Cauchy en métodos de contorno

Es indudable que las estrategias para evaluar numéricamente en el computador integrales singulares en el sentido de Valor Principal de Cauchy (VPC) constituyen uno de los aspectos clave en la precisión y confiabilidad del Método de los Elementos de Contorno. Así pues, en este trabajo se presenta un nuevo procedimiento numérico orientado a la evaluación de integrales singulares VPC y de aquellas que contengan otro tipo de singularidades, basado en una transformación bi-cúbica de coordenadas. Se desarrollan las ideas y detalles de la formulación propuesta, obteniéndose expresiones sencillas y atractivas para la cuadratura numérica y rápidamente incorporables a los códigos MEC ya existentes. Se presentan algunos ejemplos numéricos que avalan la estabilidad y precisión de los nuevos algoritmos.

La dimensión de la subjetividad en la formación del valor inmobiliario : aplicación del método de análisis de ecuaciones estructurales al mercado residencial de Lisboa

The discretionality and the appraisers’ subjectivity that characterize traditional real estate valuation are still allowed to take part in the formation of the asset price even when respecting international standards (EVS, IVS) or Appraisal Institution´s regulations (TEGOVA, RICS, etc.). The application of econometric and statistical methods to real estate valuation aims at the elimination of subjectivity on the appraisal process. But the unanswered question underneath this subject is the following: How important is the subjective component on real estate appraisal value formation? On this study Structural Equation Models (SEM) are used to determine the importance of the objective and subjective components on real estate valuation value formation as well as the weight of economic factors and the current economic context on real estate appraisal for mortgage purposes price formation. There were used two latent variables, Objective Component and Subjective Component, witch aggregate objective observed variables and subjective observed and unobserved variables, respectively. Factorial Exploratory Analysis is the statistical technique used in order to link the observed variables extracted from the valuation appraisal reports to the latent constructs derived from the theoretical model. SEM models were used to refine the model, eliminate non‐significant variables and to determine the weight of Objective and Subjective latent variables. These techniques were applied to a sample of over 11.000 real estate assets appraisal reports throughout the time period between November of 2006 and April of 2012. The real assets used on this study are located on Lisbon’s Metropolitan Area – “Grande Lisboa” –, Portugal. From this study, we conclude that Subjective Component has a considerable weight on real estate appraisal value formation and that the external factor Economic Situation has a very small impact on real estate appraisal value formation.

Sobre la formulación de las ecuaciones del campo eléctrico y magnético

Con esta breve nota no intentaré más que dar una pequeña idea de la actual controversia surgida con motivo de la aparición del libro de Chu y la consecuente discusión de los modelos más característicos de formulación de las ecuaciones de Maxwell: el modelo de Minkowski, el de Boffi o Amperiano y el de Chu.

Precisión de las ecuaciones para estimar la frecuencia cardíaca máxima en cicloergómetro

Precisión de las ecuaciones para estimar la frecuencia cardíaca máxima en cicloergómetro

Nuevas ecuaciones de cubicación para pino silvestre en "aguas vertientes" (El Espinar, Segovia).

El monte "Aguas Vertientes" (El Espinar) es el nº 138 del C.U.P. de Segovia. Se trata de una masa natural de pino silvestre, con abundante melojo en sus cotas inferiores. Su primer Proyecto de Ordenación fue redactado en 1904, estando actualmente redactada y en fase de aprobación la novena revisión. La madera de pino obtenida es elaborada en un aserradero cercano, perteneciente a la misma entidad propietaria. Durante los últimos años los responsables de dicho aserradero vienen manifestando que la cubicación en pie realizada por la administración forestal presenta errores de magnitud suficiente como para ser necesaria su corrección. Este trabajo fin de carrera ha consistido en comprobar la fiabilidad de los valores modulares empleados hasta ahora, y elaborar nuevas ecuaciones de cubicación. Para ello se ha manejado una muestra de 234 pinos repartidos en malla cuadrada y cubicados en pie mediante el método de Pressler-Bitterlich. Como resultado, se propone una ecuación de doble entrada (función de diámetro normal y altura) para todo el monte; y una ecuación de una entrada (función de diámetro normal) para cotas superiores e inferiores a 1.500 m. Los nuevos modelos son validados mediante una muestra independiente.

Monte Carlo Neutronics and Thermal-hydraulics coupling applied to Fast Reactors

Un escenario habitualmente considerado para el uso sostenible y prolongado de la energía nuclear contempla un parque de reactores rápidos refrigerados por metales líquidos (LMFR) dedicados al reciclado de Pu y la transmutación de actínidos minoritarios (MA). Otra opción es combinar dichos reactores con algunos sistemas subcríticos asistidos por acelerador (ADS), exclusivamente destinados a la eliminación de MA. El diseño y licenciamiento de estos reactores innovadores requiere herramientas computacionales prácticas y precisas, que incorporen el conocimiento obtenido en la investigación experimental de nuevas configuraciones de reactores, materiales y sistemas. A pesar de que se han construido y operado un cierto número de reactores rápidos a nivel mundial, la experiencia operacional es todavía reducida y no todos los transitorios se han podido entender completamente. Por tanto, los análisis de seguridad de nuevos LMFR están basados fundamentalmente en métodos deterministas, al contrario que las aproximaciones modernas para reactores de agua ligera (LWR), que se benefician también de los métodos probabilistas. La aproximación más usada en los estudios de seguridad de LMFR es utilizar una variedad de códigos, desarrollados a base de distintas teorías, en busca de soluciones integrales para los transitorios e incluyendo incertidumbres. En este marco, los nuevos códigos para cálculos de mejor estimación ("best estimate") que no incluyen aproximaciones conservadoras, son de una importancia primordial para analizar estacionarios y transitorios en reactores rápidos. Esta tesis se centra en el desarrollo de un código acoplado para realizar análisis realistas en reactores rápidos críticos aplicando el método de Monte Carlo. Hoy en día, dado el mayor potencial de recursos computacionales, los códigos de transporte neutrónico por Monte Carlo se pueden usar de manera práctica para realizar cálculos detallados de núcleos completos, incluso de elevada heterogeneidad material. Además, los códigos de Monte Carlo se toman normalmente como referencia para los códigos deterministas de difusión en multigrupos en aplicaciones con reactores rápidos, porque usan secciones eficaces punto a punto, un modelo geométrico exacto y tienen en cuenta intrínsecamente la dependencia angular de flujo. En esta tesis se presenta una metodología de acoplamiento entre el conocido código MCNP, que calcula la generación de potencia en el reactor, y el código de termohidráulica de subcanal COBRA-IV, que obtiene las distribuciones de temperatura y densidad en el sistema. COBRA-IV es un código apropiado para aplicaciones en reactores rápidos ya que ha sido validado con resultados experimentales en haces de barras con sodio, incluyendo las correlaciones más apropiadas para metales líquidos. En una primera fase de la tesis, ambos códigos se han acoplado en estado estacionario utilizando un método iterativo con intercambio de archivos externos. El principal problema en el acoplamiento neutrónico y termohidráulico en estacionario con códigos de Monte Carlo es la manipulación de las secciones eficaces para tener en cuenta el ensanchamiento Doppler cuando la temperatura del combustible aumenta. Entre todas las opciones disponibles, en esta tesis se ha escogido la aproximación de pseudo materiales, y se ha comprobado que proporciona resultados aceptables en su aplicación con reactores rápidos. Por otro lado, los cambios geométricos originados por grandes gradientes de temperatura en el núcleo de reactores rápidos resultan importantes para la neutrónica como consecuencia del elevado recorrido libre medio del neutrón en estos sistemas. Por tanto, se ha desarrollado un módulo adicional que simula la geometría del reactor en caliente y permite estimar la reactividad debido a la expansión del núcleo en un transitorio. éste módulo calcula automáticamente la longitud del combustible, el radio de la vaina, la separación de los elementos de combustible y el radio de la placa soporte en función de la temperatura. éste efecto es muy relevante en transitorios sin inserción de bancos de parada. También relacionado con los cambios geométricos, se ha implementado una herramienta que, automatiza el movimiento de las barras de control en busca d la criticidad del reactor, o bien calcula el valor de inserción axial las barras de control. Una segunda fase en la plataforma de cálculo que se ha desarrollado es la simulació dinámica. Puesto que MCNP sólo realiza cálculos estacionarios para sistemas críticos o supercríticos, la solución más directa que se propone sin modificar el código fuente de MCNP es usar la aproximación de factorización de flujo, que resuelve por separado la forma del flujo y la amplitud. En este caso se han estudiado en profundidad dos aproximaciones: adiabática y quasiestática. El método adiabático usa un esquema de acoplamiento que alterna en el tiempo los cálculos neutrónicos y termohidráulicos. MCNP calcula el modo fundamental de la distribución de neutrones y la reactividad al final de cada paso de tiempo, y COBRA-IV calcula las propiedades térmicas en el punto intermedio de los pasos de tiempo. La evolución de la amplitud de flujo se calcula resolviendo las ecuaciones de cinética puntual. Este método calcula la reactividad estática en cada paso de tiempo que, en general, difiere de la reactividad dinámica que se obtendría con la distribución de flujo exacta y dependiente de tiempo. No obstante, para entornos no excesivamente alejados de la criticidad ambas reactividades son similares y el método conduce a resultados prácticos aceptables. Siguiendo esta línea, se ha desarrollado después un método mejorado para intentar tener en cuenta el efecto de la fuente de neutrones retardados en la evolución de la forma del flujo durante el transitorio. El esquema consiste en realizar un cálculo cuasiestacionario por cada paso de tiempo con MCNP. La simulación cuasiestacionaria se basa EN la aproximación de fuente constante de neutrones retardados, y consiste en dar un determinado peso o importancia a cada ciclo computacial del cálculo de criticidad con MCNP para la estimación del flujo final. Ambos métodos se han verificado tomando como referencia los resultados del código de difusión COBAYA3 frente a un ejercicio común y suficientemente significativo. Finalmente, con objeto de demostrar la posibilidad de uso práctico del código, se ha simulado un transitorio en el concepto de reactor crítico en fase de diseño MYRRHA/FASTEF, de 100 MW de potencia térmica y refrigerado por plomo-bismuto. ABSTRACT Long term sustainable nuclear energy scenarios envisage a fleet of Liquid Metal Fast Reactors (LMFR) for the Pu recycling and minor actinides (MAs) transmutation or combined with some accelerator driven systems (ADS) just for MAs elimination. Design and licensing of these innovative reactor concepts require accurate computational tools, implementing the knowledge obtained in experimental research for new reactor configurations, materials and associated systems. Although a number of fast reactor systems have already been built, the operational experience is still reduced, especially for lead reactors, and not all the transients are fully understood. The safety analysis approach for LMFR is therefore based only on deterministic methods, different from modern approach for Light Water Reactors (LWR) which also benefit from probabilistic methods. Usually, the approach adopted in LMFR safety assessments is to employ a variety of codes, somewhat different for the each other, to analyze transients looking for a comprehensive solution and including uncertainties. In this frame, new best estimate simulation codes are of prime importance in order to analyze fast reactors steady state and transients. This thesis is focused on the development of a coupled code system for best estimate analysis in fast critical reactor. Currently due to the increase in the computational resources, Monte Carlo methods for neutrons transport can be used for detailed full core calculations. Furthermore, Monte Carlo codes are usually taken as reference for deterministic diffusion multigroups codes in fast reactors applications because they employ point-wise cross sections in an exact geometry model and intrinsically account for directional dependence of the ux. The coupling methodology presented here uses MCNP to calculate the power deposition within the reactor. The subchannel code COBRA-IV calculates the temperature and density distribution within the reactor. COBRA-IV is suitable for fast reactors applications because it has been validated against experimental results in sodium rod bundles. The proper correlations for liquid metal applications have been added to the thermal-hydraulics program. Both codes are coupled at steady state using an iterative method and external files exchange. The main issue in the Monte Carlo/thermal-hydraulics steady state coupling is the cross section handling to take into account Doppler broadening when temperature rises. Among every available options, the pseudo materials approach has been chosen in this thesis. This approach obtains reasonable results in fast reactor applications. Furthermore, geometrical changes caused by large temperature gradients in the core, are of major importance in fast reactor due to the large neutron mean free path. An additional module has therefore been included in order to simulate the reactor geometry in hot state or to estimate the reactivity due to core expansion in a transient. The module automatically calculates the fuel length, cladding radius, fuel assembly pitch and diagrid radius with the temperature. This effect will be crucial in some unprotected transients. Also related to geometrical changes, an automatic control rod movement feature has been implemented in order to achieve a just critical reactor or to calculate control rod worth. A step forward in the coupling platform is the dynamic simulation. Since MCNP performs only steady state calculations for critical systems, the more straight forward option without modifying MCNP source code, is to use the flux factorization approach solving separately the flux shape and amplitude. In this thesis two options have been studied to tackle time dependent neutronic simulations using a Monte Carlo code: adiabatic and quasistatic methods. The adiabatic methods uses a staggered time coupling scheme for the time advance of neutronics and the thermal-hydraulics calculations. MCNP computes the fundamental mode of the neutron flux distribution and the reactivity at the end of each time step and COBRA-IV the thermal properties at half of the the time steps. To calculate the flux amplitude evolution a solver of the point kinetics equations is used. This method calculates the static reactivity in each time step that in general is different from the dynamic reactivity calculated with the exact flux distribution. Nevertheless, for close to critical situations, both reactivities are similar and the method leads to acceptable practical results. In this line, an improved method as an attempt to take into account the effect of delayed neutron source in the transient flux shape evolutions is developed. The scheme performs a quasistationary calculation per time step with MCNP. This quasistationary simulations is based con the constant delayed source approach, taking into account the importance of each criticality cycle in the final flux estimation. Both adiabatic and quasistatic methods have been verified against the diffusion code COBAYA3, using a theoretical kinetic exercise. Finally, a transient in a critical 100 MWth lead-bismuth-eutectic reactor concept is analyzed using the adiabatic method as an application example in a real system.