Modelado detallado de una central nuclear CANDU con CITVAP.

La línea de cálculo de INVAP consiste principalmente de los códigos CONDOR y CITVAP. Este último es la versión mejorada del código CITATION II que resuelve la ecuación de difusión neutrónica multigrupo por el método de diferencias finitas. CITVAP es ampliamente usado para estudiar reactores de investigación y reactores de potencia tales como PWR, BWR, VVER y últimamente se implemento nuevas funciones para estudiar una central PHWR tipo Atucha. Siguiendo con la línea de reactores PHWR, en este trabajo se estudian las capacidades y deficiencias del código de núcleo CITVAP para modelar una central nuclear tipo CANDU. Se plantean mejoras a realizar para un manejo mas eficiente desde el punto de vista del usuario, tanto de la gestión de combustibles, movimientos de barras de control y zonas líquidas como mejoras en el modelo termohidraulico. La metodología consiste en validar la línea de cálculo de INVAP, contrastando los resultados con el benchmark IAEA-tecdoc-887. El proceso de validación consiste en cálculos de celda en dos y tres dimensiones usando los códigos CONDOR y SERPENT respectivamente, obtención de secciones eficaces macroscópicas en función del quemado y cálculos de núcleo para distintas configuraciones de los dispositivos de control usando un núcleo fresco y una distribución de quemado en equilibrio. Se analizan las dificultades que se presentan al modelar el núcleo con las capacidades actuales del código y se plantean posibles soluciones a implementar. Para un estudio completo de un reactor CANDU, se estudian tres de la características distintivas de este tipo de reactor: la termohidraulica, la gestión de combustibles y los dispositivos de control de reactividad, distribución de potencia y apagado.

Modelado detallado de una central nuclear CANDU con CITVAP.

La línea de cálculo de INVAP consiste principalmente de los códigos CONDOR y CITVAP. Este último es la versión mejorada del código CITATION II que resuelve la ecuación de difusión neutrónica multigrupo por el método de diferencias finitas. CITVAP es ampliamente usado para estudiar reactores de investigación y reactores de potencia tales como PWR, BWR, VVER y últimamente se implemento nuevas funciones para estudiar una central PHWR tipo Atucha. Siguiendo con la línea de reactores PHWR, en este trabajo se estudian las capacidades y deficiencias del código de núcleo CITVAP para modelar una central nuclear tipo CANDU. Se plantean mejoras a realizar para un manejo mas eficiente desde el punto de vista del usuario, tanto de la gestión de combustibles, movimientos de barras de control y zonas líquidas como mejoras en el modelo termohidraulico. La metodología consiste en validar la línea de cálculo de INVAP, contrastando los resultados con el benchmark IAEA-tecdoc-887. El proceso de validación consiste en cálculos de celda en dos y tres dimensiones usando los códigos CONDOR y SERPENT respectivamente, obtención de secciones eficaces macroscópicas en función del quemado y cálculos de núcleo para distintas configuraciones de los dispositivos de control usando un núcleo fresco y una distribución de quemado en equilibrio. Se analizan las dificultades que se presentan al modelar el núcleo con las capacidades actuales del código y se plantean posibles soluciones a implementar. Para un estudio completo de un reactor CANDU, se estudian tres de la características distintivas de este tipo de reactor: la termohidraulica, la gestión de combustibles y los dispositivos de control de reactividad, distribución de potencia y apagado.

Numerical comparison between Maxwell stress method and equivalent multipole approach for calculation of the dielectrophoretic force in octupolar cell traps

This work presents detailed numerical calculations of the dielectrophoretic force in octupolar traps designed for single-cell trapping. A trap with eight planar electrodes is studied for spherical and ellipsoidal particles using an indirect implementation of the boundary element method (BEM). Multipolar approximations of orders one to three are compared with the full Maxwell stress tensor (MST) calculation of the electrical force on spherical particles. Ellipsoidal particles are also studied, but in their case only the dipolar approximation is available for comparison with the MST solution. The results show that the full MST calculation is only required in the study of non-spherical particles.

Further development of discrete computational techniques for calculation of restricted diffusion propagators in porous media

Magnetic resonance is a well-established tool for structural characterisation of porous media. Features of pore-space morphology can be inferred from NMR diffusion-diffraction plots or the time-dependence of the apparent diffusion coefficient. Diffusion NMR signal attenuation can be computed from the restricted diffusion propagator, which describes the distribution of diffusing particles for a given starting position and diffusion time. We present two techniques for efficient evaluation of restricted diffusion propagators for use in NMR porous-media characterisation. The first is the Lattice Path Count (LPC). Its physical essence is that the restricted diffusion propagator connecting points A and B in time t is proportional to the number of distinct length-t paths from A to B. By using a discrete lattice, the number of such paths can be counted exactly. The second technique is the Markov transition matrix (MTM). The matrix represents the probabilities of jumps between every pair of lattice nodes within a single timestep. The propagator for an arbitrary diffusion time can be calculated as the appropriate matrix power. For periodic geometries, the transition matrix needs to be defined only for a single unit cell. This makes MTM ideally suited for periodic systems. Both LPC and MTM are closely related to existing computational techniques: LPC, to combinatorial techniques; and MTM, to the Fokker-Planck master equation. The relationship between LPC, MTM and other computational techniques is briefly discussed in the paper. Both LPC and MTM perform favourably compared to Monte Carlo sampling, yielding highly accurate and almost noiseless restricted diffusion propagators. Initial tests indicate that their computational performance is comparable to that of finite element methods. Both LPC and MTM can be applied to complicated pore-space geometries with no analytic solution. We discuss the new methods in the context of diffusion propagator calculation in porous materials and model biological tissues.

Estimation of red blood cell lifespan from alveolar carbon monoxide measurements

Measurement of alveolar carbon monoxide (CO) presents a facile technique to estimate the lifespan, L, of red blood cells (RBCs) in vivo. Several recent studies employ this technique and calculate L (in days) using the expression, L = 13.8 (Hb)/P-CO(end), where (Hb) is the concentration (in g/dL) of hemoglobin in blood, and P-CO(end) is the endogenous production of CO (in ppm). Implicit in this calculation is the assumption that the fraction, f, of endogenous CO production due to RBC turnover is a constant equal to 0.7, which yields the expected RBC lifespan, L approximate to 120 days, in normal controls. In anemic patients, however, enhanced RBC turnover may increase f substantially above 0.7. The above expression then overestimates L. Here, we deriv an alternative tive expression, L = 3390[Hb]/322P(CO (end)-110, that accounts explicitly for the dependence of f on the rate of RBC turnover and thereby provides more accurate estimates of L without requiring additional measurements. Using the latter expression, we recalculate L from recent measurements on hepatitis C virus infected patients undergoing treatment with ribavirin. We find that our estimates of L in these patients (39 +/- 13 days) are significantly lower than current estimates (46 +/- 14 days), indicating that ribavirin affects RBC survival more severely than expected from current studies. Our expression for L is simple to employ in a clinical setting and would render the broadly applicable technique of alveolar CO measurement for the estimation of RBC lifespan more accurate.

Methods for melting temperature calculation

Melting temperature calculation has important applications in the theoretical study of phase diagrams and computational materials screenings. In this thesis, we present two new methods, i.e., the improved Widom's particle insertion method and the small-cell coexistence method, which we developed in order to capture melting temperatures both accurately and quickly.

We propose a scheme that drastically improves the efficiency of Widom's particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.

We propose the small-cell coexistence method based on the statistical analysis of small-size coexistence MD simulations. It eliminates the risk of a metastable superheated solid in the fast-heating method, while also significantly reducing the computer cost relative to the traditional large-scale coexistence method. Using empirical potentials, we validate the method and systematically study the finite-size effect on the calculated melting points. The method converges to the exact result in the limit of a large system size. An accuracy within 100 K in melting temperature is usually achieved when the simulation contains more than 100 atoms. DFT examples of Tantalum, high-pressure Sodium, and ionic material NaCl are shown to demonstrate the accuracy and flexibility of the method in its practical applications. The method serves as a promising approach for large-scale automated material screening in which the melting temperature is a design criterion.

We present in detail two examples of refractory materials. First, we demonstrate how key material properties that provide guidance in the design of refractory materials can be accurately determined via ab initio thermodynamic calculations in conjunction with experimental techniques based on synchrotron X-ray diffraction and thermal analysis under laser-heated aerodynamic levitation. The properties considered include melting point, heat of fusion, heat capacity, thermal expansion coefficients, thermal stability, and sublattice disordering, as illustrated in a motivating example of lanthanum zirconate (La₂Zr₂O₇). The close agreement with experiment in the known but structurally complex compound La₂Zr₂O₇ provides good indication that the computation methods described can be used within a computational screening framework to identify novel refractory materials. Second, we report an extensive investigation into the melting temperatures of the Hf-C and Hf-Ta-C systems using ab initio calculations. With melting points above 4000 K, hafnium carbide (HfC) and tantalum carbide (TaC) are among the most refractory binary compounds known to date. Their mixture, with a general formula Ta_xHf_1-xC_y, is known to have a melting point of 4215 K at the composition Ta₄HfC₅, which has long been considered as the highest melting temperature for any solid. Very few measurements of melting point in tantalum and hafnium carbides have been documented, because of the obvious experimental difficulties at extreme temperatures. The investigation lets us identify three major chemical factors that contribute to the high melting temperatures. Based on these three factors, we propose and explore a new class of materials, which, according to our ab initio calculations, may possess even higher melting temperatures than Ta-Hf-C. This example also demonstrates the feasibility of materials screening and discovery via ab initio calculations for the optimization of "higher-level" properties whose determination requires extensive sampling of atomic configuration space.

Quantum Mechanics Studies of Fuel Cell Catalysts and Proton Conducting Ceramics with Validation by Experiment

We carried out quantum mechanics (QM) studies aimed at improving the performance of hydrogen fuel cells. This led to predictions of improved materials, some of which were subsequently validated with experiments by our collaborators.

In part I, the challenge was to find a replacement for the Pt cathode that would lead to improved performance for the Oxygen Reduction Reaction (ORR) while remaining stable under operational conditions and decreasing cost. Our design strategy was to find an alloy with composition Pt3M that would lead to surface segregation such that the top layer would be pure Pt, with the second and subsequent layers richer in M. Under operating conditions we expect the surface to have significant O and/or OH chemisorbed on the surface, and hence we searched for M that would remain segregated under these conditions. Using QM we examined surface segregation for 28 Pt₃M alloys, where M is a transition metal. We found that only Pt₃Os and Pt₃Ir showed significant surface segregation when O and OH are chemisorbed on the catalyst surfaces. This result indicates that Pt₃Os and Pt3Ir favor formation of a Pt-skin surface layer structure that would resist the acidic electrolyte corrosion during fuel cell operation environments. We chose to focus on Os because the phase diagram for Pt-Ir indicated that Pt-Ir could not form a homogeneous alloy at lower temperature. To determine the performance for ORR, we used QM to examine all intermediates, reaction pathways, and reaction barriers involved in the processes for which protons from the anode reactions react with O₂ to form H₂O. These QM calculations used our Poisson-Boltzmann implicit solvation model include the effects of the solvent (water with dielectric constant 78 with pH 7 at 298K). We found that the rate determination step (RDS) was the O_ad hydration reaction (O_ad + H₂O_ad -> OH_ad + OH_ad) in both cases, but that the barrier for pure Pt of 0.50 eV is reduced to 0.48 eV for Pt₃Os, which at 80 degrees C would increase the rate by 218%. We collaborated with the Pu-Wei Wu’s group to carry out experiments, where we found that the dealloying process-treated Pt2Os catalyst showed two-fold higher activity at 25 degrees C than pure Pt and that the alloy had 272% improved stability, validating our theoretical predictions.

We also carried out similar QM studies followed by experimental validation for the Os/Pt core-shell catalyst fabricated by the underpotential deposition (UPD) method. The QM results indicated that the RDS for ORR is a compromise between the OOH formation step (0.37 eV for Pt, 0.23 eV for Pt_2ML/Os core-shell) and H₂O formation steps (0.32 eV for Pt, 0.22 eV for Pt_2ML/Os core-shell). We found that Pt_2ML/Os has the highest activity (compared to pure Pt and to the Pt₃Os alloy) because the 0.37 eV barrier decreases to 0.23 eV. To understand what aspects of the core shell structure lead to this improved performance, we considered the effect on ORR of compressing the alloy slab to the dimensions of pure Pt. However this had little effect, with the same RDS barrier 0.37 eV. This shows that the ligand effect (the electronic structure modification resulting from the Os substrate) plays a more important role than the strain effect, and is responsible for the improved activity of the core- shell catalyst. Experimental materials characterization proves the core-shell feature of our catalyst. The electrochemical experiment for Pt_2ML/Os/C showed 3.5 to 5 times better ORR activity at 0.9V (vs. NHE) in 0.1M HClO₄ solution at 25 degrees C as compared to those of commercially available Pt/C. The excellent correlation between experimental half potential and the OH binding energies and RDS barriers validate the feasibility of predicting catalyst activity using QM calculation and a simple Langmuir–Hinshelwood model.

In part II, we used QM calculations to study methane stream reforming on a Ni-alloy catalyst surfaces for solid oxide fuel cell (SOFC) application. SOFC has wide fuel adaptability but the coking and sulfur poisoning will reduce its stability. Experimental results suggested that the Ni4Fe alloy improves both its activity and stability compared to pure Ni. To understand the atomistic origin of this, we carried out QM calculations on surface segregation and found that the most stable configuration for Ni₄Fe has a Fe atom distribution of (0%, 50%, 25%, 25%, 0%) starting at the bottom layer. We calculated that the binding of C atoms on the Ni4Fe surface is 142.9 Kcal/mol, which is about 10 Kcal/mol weaker compared to the pure Ni surface. This weaker C binding energy is expected to make coke formation less favorable, explaining why Ni₄Fe has better coking resistance. This result confirms the experimental observation. The reaction energy barriers for CHx decomposition and C binding on various alloy surface, Ni₄X (X=Fe, Co, Mn, and Mo), showed Ni₄Fe, Ni₄Co, and Fe₄Mn all have better coking resistance than pure Ni, but that only Ni₄Fe and Fe₄Mn have (slightly) improved activity compared to pure Ni.

In part III, we used QM to examine the proton transport in doped perovskite-ceramics. Here we used a 2x2x2 supercell of perovskite with composition Ba₈X₇M₁(OH)₁O₂₃ where X=Ce or Zr and M=Y, Gd, or Dy. Thus in each case a 4⁺ X is replace by a 3⁺ M plus a proton on one O. Here we predicted the barriers for proton diffusion allowing both includes intra-octahedron and inter-octahedra proton transfer. Without any restriction, we only observed the inter-octahedra proton transfer with similar energy barrier as previous computational work but 0.2 eV higher than experimental result for Y doped zirconate. For one restriction in our calculations is that the O_donor-O_acceptor atoms were kept at fixed distances, we found that the barrier difference between cerates/zirconates with various dopants are only 0.02~0.03 eV. To fully address performance one would need to examine proton transfer at grain boundaries, which will require larger scale ReaxFF reactive dynamics for systems with millions of atoms. The QM calculations used here will be used to train the ReaxFF force field.

Analysis of leakage current in GaAs micro-solar cell arrays

The output characteristics of micro-solar cell arrays are analyzed on the basis of a modified model in which the shunt resistance between cell lines results in current leakage. The modification mainly consists of adding a shunt resistor network to the traditional model. The obtained results agree well with the reported experimental results. The calculation results demonstrate that leakage current in substrate affects seriously the performance of GaAs micro- solar cell arrays. The performance of arrays can be improved by reducing the number of cells per line. In addition, at a certain level of integration, an appropriate space occupancy rate of the single cell is recommended for ensuring high open circuit voltages, and it is more appropriate to set the rates at 80%-90% through the calculation.

Evaluation of thermal radiation dependent performance of GaSb thermophotovoltaic cell based on an analytical absorption coefficient model

Applying the model dielectric function method, we have expressed the absorption coefficient of GaSb analytically at room temperature relating to the contribution of various critical points of its electronic band structure. The calculated absorption spectrum shows good agreement with the reported experimental data obtained by spectral ellipsometry on nominally undoped sample. Based on this analytical absorption spectrum, we have qualitatively evaluated the response of active absorbing layer structure and its photoelectric conversion properties of GaSb thermophotovoltaic device on the perturbation of external thermal radiation induced by the varying radiator temperature or emissivity. Our calculation has demonstrated that desirable thickness to achieve the maximum conversion efficiency should be decreased with the increment of radiator temperature and the performance degradation brought by any structure deviation from its optimal one would be stronger meanwhile. For the popular radiator temperature, no more than 1500 K in a real solar thermophotovoltaic system, and typical doping profile in GaSb cell, a reasonable absorbing layer structure parameter should be controlled within 100-300 nm for the emitter while 3000-5000 nm for the base.

Effect of partial site occupation on the magnetic properties of Nd2Fe17-xSix by supercell calculation

The magnetic properties of the Nd2Fe17-xSix intermetallic compounds are studied by means of spin-polarized supercell calculations in which the selected sites of substitution are close to the situations in real samples. It is shown that the average Fe moment increases with x and saturates near x = 3. This correlates quite well with the experimental dependence of Te on x. The difference between supercell and unit cell calculations are pointed out and the influence of Si atoms on the density of states of the nearby Fe atoms is emphasized. (C) 1997 American Institute of Physics.

Analysis of leakage current in GaAs micro-solar cell arrays

The output characteristics of micro-solar cell arrays are analyzed on the basis of a modified model in which the shunt resistance between cell lines results in current leakage. The modification mainly consists of adding a shunt resistor network to the traditional model. The obtained results agree well with the reported experimental results. The calculation results demonstrate that leakage current in substrate affects seriously the performance of GaAs micro- solar cell arrays. The performance of arrays can be improved by reducing the number of cells per line. In addition, at a certain level of integration, an appropriate space occupancy rate of the single cell is recommended for ensuring high open circuit voltages, and it is more appropriate to set the rates at 80%-90% through the calculation.

INVESTIGATION OF CETYLPYRIDINIUM BROMIDE ADSORPTION AT A GLASSY-CARBON ELECTRODE SURFACE BY SPECTROELECTROCHEMISTRY WITH A LONG OPTICAL-PATH LENGTH THIN-LAYER CELL

The adsorption of cationic surfactant cetylpyridinium bromide (CPB) on a glassy carbon (GC) electrode surface has been studied by spectroelectrochemistry with a long optical path length thin-layer cell (LOPTLC) for the first time. A fine adsorption isotherm of CPB molecules from an aqueous solution containing 0.10 M KBr has been obtained over the range of (1.00-8.00) x 10(-5) M. From theoretical calculation and experimental data, adsorption of CPB on the GC electrode surface shows four distinct orientations and three large orientation transitions. Compared with the ordinary isotherm, the differential isotherm is more characteristic and would be suitable for the study of orientation transitions of organic compounds. With a theoretical treatment of the adsorption isotherm, four orientations of adsorbed CPB on a GC electrode surface coincide with the Frumkin-Langmuir type. From adsorption parameters the Frumkin-Langmuir equations, the adsorption free energy and, therefore, the equilibrium constants of orientation transitions of the CPB molecule can be obtained.

Manufacturing of uniformly-sized silicon particles for solar cell from molten metal jet by electromagnetic pinch force

Spherical silicon solar cells are expected to serve as a technology to reduce silicon usage of photovoltaic (PV) power systems[1, 2, 3]. In order to establish the spherical silicon solar cell, a manufacturing method of uniformly sized silicon particles of 1mm in diameter is required. However, it is difficult to mass-produce the mono-sized silicon particles at low cost by existent processes now. We proposed a new method to generate liquid metal droplets uniformly by applying electromagnetic pinch force to a liquid metal jet[4]. The electromagnetic force was intermittently applied to the liquid metal jet issued from a nozzle in order to fluctuate the surface of the jet. As the fluctuation grew, the liquid jet was broken up into small droplets according to a frequency of the intermittent electromagnetic force. Firstly, a preliminary experiment was carried out. A single pulse current was applied instantaneously to a single turn coil around a molten gallium jet. It was confirmed that the jet could be split up by pinch force generated by the current. And then, electromagnetic pinch force was applied intermittently to the jet. It was found that the jet was broken up into mono-sized droplets in the case of a force frequency was equal to a critical frequency[5], which corresponds to a natural disturbance wave length of the jet. Numerical simulations of the droplet generation from the liquid jet were then carried out, which consisted of an electromagnetic analysis and a fluid flow calculation with a free surface of the jet. The simulation results were compared with the experiments and the agreement between the two was quite good.

Optical and mechanical analysis on a biological cell in optical tweezers

La réponse mécanique d’une cellule à une force externe permet d’inférer sa structure et fonction. Les pinces optiques s’avèrent une approche particulièrement attrayante pour la manipulation et caractérisation biophysique sophistiquée des cellules de façon non invasive. Cette thèse explore l’utilisation de trois types de pinces optiques couramment utilisées : 1) statiques (static), 2) à exposition partagée (time-sharing) et 3) oscillantes (oscillating). L’utilisation d’un code basé sur la méthode des éléments finis en trois dimensions (3DFEM) nous permet de modéliser ces trois types de piégeage optique afin d’extraire les propriétés mécaniques cellulaires à partir des expériences. La combinaison des pinces optiques avec la mécanique des cellules requiert des compétences interdisciplinaires. Une revue des approches expérimentales sur le piégeage optique et les tests unicellulaires est présentée. Les bases théoriques liant l’interaction entre la force radiative optique et la réponse mécanique de la cellule aussi. Pour la première fois, une simulation adaptée (3DFEM) incluant la diffusion lumineuse et la distribution du stress radiatif permet de prédire la déformation d’une cellule biconcave –analogue aux globules rouges—dans un piège statique double (static dual-trap). À l’équilibre, on observe que la déformation finale est donnée par l’espacement entre les deux faisceaux lasers: la cellule peut être étirée ou même comprimée. L’exposition partagée (time-sharing) est la technique qui permet de maintenir plusieurs sites de piégeage simultanément à partir du même faisceau laser. Notre analyse quantitative montre que, même oscillantes, la force optique et la déformation sont omniprésentes dans la cellule : la déformation viscoélastique et la dissipation de l’énergie sont analysées. Une autre cellule-type, la tige cubique, est étudiée : cela nous permet d’élucider de nouvelles propriétés sur la symétrie de la réponse mécanique. Enfin, l’analyse de la déformation résolue en temps dans un piége statique ou à exposition partagée montre que la déformation dépend simultanément de la viscoélasticité, la force externe et sa forme tridimensionnelle. La technique à force oscillante (oscillating tweezers) montre toutefois un décalage temporel, entre la force et la déformation, indépendant de la forme 3D; cette approche donnerait directement accès au tenseur viscoélastique complexe de la cellule.

Calculation of the magnetic field penetration depth for high-Tc cuprate superconductors based on the Interlayer Pair Tunneling model /

In this work, the magnetic field penetration depth for high-Tc cuprate superconductors is calculated using a recent Interlayer Pair Tunneling (ILPT) model proposed by Chakravarty, Sudb0, Anderson, and Strong [1] to explain high temperature superconductivity. This model involves a "hopping" of Cooper pairs between layers of the unit cell which acts to amplify the pairing mechanism within the planes themselves. Recent work has shown that this model can account reasonably well for the isotope effect and the dependence of Tc on nonmagnetic in-plane impurities [2] , as well as the Knight shift curves [3] and the presence of a magnetic peak in the neutron scattering intensity [4]. In the latter case, Yin et al. emphasize that the pair tunneling must be the dominant pairing mechanism in the high-Tc cuprates in order to capture the features found in experiments. The goal of this work is to determine whether or not the ILPT model can account for the experimental observations of the magnetic field penetration depth in YBa2Cu307_a7. Calculations are performed in the weak and strong coupling limits, and the efi"ects of both small and large strengths of interlayer pair tunneling are investigated. Furthermore, as a follow up to the penetration depth calculations, both the neutron scattering intensity and the Knight shift are calculated within the ILPT formalism. The aim is to determine if the ILPT model can yield results consistent with experiments performed for these properties. The results for all three thermodynamic properties considered are not consistent with the notion that the interlayer pair tunneling must be the dominate pairing mechanism in these high-Tc cuprate superconductors. Instead, it is found that reasonable agreement with experiments is obtained for small strengths of pair tunneling, and that large pair tunneling yields results which do not resemble those of the experiments.