A Krylov-Schur solution of the eigenvalue problem for the neutron diffusion equation discretized with the Raviart-Thomas method

[EN] Mixed-dual formulations of the finite element method were successfully applied to the neutron diffusion equation, such as the Raviart¿Thomas method in Cartesian geometry and the Raviart¿Thomas¿Schneider in hexagonal geometry. Both methods obtain system matrices which are suitable for solving the eigenvalue problem with the preconditioned power method. This method is very fast and optimized, but only for the calculation of the fundamental mode. However, the determination of non-fundamental modes is important for modal analysis, instabilities, and fluctuations of nuclear reactors. So, effective and fast methods are required for solving eigenvalue problems. The most effective methods are those based on Krylov subspaces projection combined with restart, such as Krylov¿Schur. In this work, a Krylov¿Schur method has been applied to the neutron diffusion equation, discretized with the Raviart¿Thomas and Raviart¿Thomas¿Schneider methods.This work has been partially supported by the Spanish Ministerio de Eduacion Cultura y Deporte [grant number FPU13/01009]; Spanish Ministerio de Ciencia e Innovacion [project number ENE2014-59442-P]; Spanish Ministerio de Economia y Competitividad and the European Fondo Europeo de Desarrollo Regional (FEDER) [project number ENE2015-68353-P (MINECO/FEDER)]; Generalitat Valenciana [project number PROMETEOII/2014/008]; Universitat Politecnica de Valencia [project number UPPTE/2012/118]; Spanish Ministerio de Economia y Competitividad [project number TIN2016-75985-P].Bernal-Garcia, A.; Hébert, A.; Roman, JE.; Miró Herrero, R.; Verdú Martín, GJ. (2017). A Krylov-Schur solution of the eigenvalue problem for the neutron diffusion equation discretized with the Raviart-Thomas method. Journal of Nuclear Science and Technology. 54(10):1085-1094. https://doi.org/10.1080/00223131.2017.1344577S108510945410Hébert, A. (1993). Application of a dual variational formulation to finite element reactor calculations. Annals of Nuclear Energy, 20(12), 823-845. doi:10.1016/0306-4549(93)90076-2Hébert, A. (2008). A Raviart–Thomas–Schneider solution of the diffusion equation in hexagonal geometry. Annals of Nuclear Energy, 35(3), 363-376. doi:10.1016/j.anucene.2007.07.016Hébert, A. (1986). Preconditioning the Power Method for Reactor Calculations. Nuclear Science and Engineering, 94(1), 1-11. doi:10.13182/nse86-a17111Verdú, G., Ginestar, D., Vidal, V., & Muñoz-Cobo, J. L. (1994). 3D λ-modes of the neutron-diffusion equation. Annals of Nuclear Energy, 21(7), 405-421. doi:10.1016/0306-4549(94)90041-8Miró, R., Ginestar, D., Verdú, G., & Hennig, D. (2002). A nodal modal method for the neutron diffusion equation. Application to BWR instabilities analysis. Annals of Nuclear Energy, 29(10), 1171-1194. doi:10.1016/s0306-4549(01)00103-7Hébert A. Applied reactor physics. 2nd ed. Montréal: Presses Internationales Polytechnique; 2016. p. 368–369.Döring, M. G., Kalkkuhl, J. C., & Schröder, W. (1993). Subspace Iteration for Nonsymmetric Eigenvalue Problems Applied to the λ-Eigenvalue Problem. Nuclear Science and Engineering, 115(3), 244-252. doi:10.13182/nse93-a24053Modak, R. S., & Jain, V. K. (1996). Sub-space iteration scheme for the evaluation of λ-modes of finite-differenced multi-group neutron diffusion equations. Annals of Nuclear Energy, 23(3), 229-237. doi:10.1016/0306-4549(95)00015-6Singh, K. P., Modak, R. S., Degweker, S. B., & Singh, K. (2009). Iterative schemes for obtaining dominant alpha-modes of the neutron diffusion equation. Annals of Nuclear Energy, 36(8), 1086-1092. doi:10.1016/j.anucene.2009.05.006Gupta, A., & Modak, R. S. (2011). Evaluation of dominant time-eigenvalues of neutron transport equation by Meyer’s sub-space iterations. Annals of Nuclear Energy, 38(7), 1680-1686. doi:10.1016/j.anucene.2011.02.016Kópházi, J., & Lathouwers, D. (2012). Three-dimensional transport calculation of multiple alpha modes in subcritical systems. Annals of Nuclear Energy, 50, 167-174. doi:10.1016/j.anucene.2012.06.021VERDÚ, G., GINESTAR, D., ROMÁN, J., & VIDAL, V. (2010). 3D Alpha Modes of a Nuclear Power Reactor. Journal of Nuclear Science and Technology, 47(5), 501-514. doi:10.1080/18811248.2010.9711641Lathouwers, D. (2003). Iterative computation of time-eigenvalues of the neutron transport equation. Annals of Nuclear Energy, 30(17), 1793-1806. doi:10.1016/s0306-4549(03)00151-8Warsa, J. S., Wareing, T. A., Morel, J. E., McGhee, J. M., & Lehoucq, R. B. (2004). Krylov Subspace Iterations for Deterministick-Eigenvalue Calculations. Nuclear Science and Engineering, 147(1), 26-42. doi:10.13182/nse04-1Verdu, G., Miro, R., Ginestar, D., & Vidal, V. (1999). The implicit restarted Arnoldi method, an efficient alternative to solve the neutron diffusion equation. Annals of Nuclear Energy, 26(7), 579-593. doi:10.1016/s0306-4549(98)00077-2Lehoucq, R. B., Sorensen, D. C., & Yang, C. (1998). ARPACK Users’ Guide. doi:10.1137/1.9780898719628Boer, B., Lathouwers, D., Kloosterman, J. L., Van Der Hagen, T. H. J. J., & Strydom, G. (2010). Validation of the DALTON-THERMIX Code System with Transient Analyses of the HTR-10 and Application to the PBMR. Nuclear Technology, 170(2), 306-321. doi:10.13182/nt10-a9485Hernandez, V., Roman, J. E., & Vidal, V. (2005). SLEPc. ACM Transactions on Mathematical Software, 31(3), 351-362. doi:10.1145/1089014.1089019Hernández, V., Román, J. E., & Vidal, V. (2003). SLEPc: Scalable Library for Eigenvalue Problem Computations. High Performance Computing for Computational Science — VECPAR 2002, 377-391. doi:10.1007/3-540-36569-9_25Stewart, G. W. (2002). A Krylov--Schur Algorithm for Large Eigenproblems. SIAM Journal on Matrix Analysis and Applications, 23(3), 601-614. doi:10.1137/s0895479800371529Chao, Y. A., & Shatilla, Y. A. (1995). Conformal Mapping and Hexagonal Nodal Methods —II: Implementation in the ANC-H Code. Nuclear Science and Engineering, 121(2), 210-225. doi:10.13182/nse95-a2855

Endangered right whales enhance primary productivity in the bay of fundy

This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Marine mammals have recently been documented as important facilitators of rapid and efficient nutrient recycling in coastal and offshore waters. Whales enhance phytoplankton nutrition by releasing fecal plumes near the surface after feeding and by migrating from highly productive, high-latitude feeding areas to low-latitude nutrient-poor calving areas. In this study, we measured NH4 + and PO4 3- release rates from the feces of North Atlantic right whales (Eubalaena glacialis), a highly endangered baleen whale. Samples for this species were primarily collected by locating aggregations of whales in surface- Active groups (SAGs), which typically consist of a central female surrounded by males competing for sexual activity. When freshly collected feces were incubated in seawater, high initial rates of N release were generally observed, which decreased to near zero within 24 hours of sampling, a pattern that is consistent with the active role of gut microflora on fecal particles. We estimate that at least 10% of particulate N in whale feces becomes available as NH4 + within 24 hours of defecation. Phosphorous was also abundant in fecal samples: Initial release rates of PO4 3- were higher than for NH4 +, yielding low N/P nutrient ratios over the course of our experiments. The rate of PO4 3- release was thus more than sufficient to preclude the possibility that nitrogenous nutrients supplied by whales would lead to phytoplankton production limited by P availability. Phytoplankton growth experiments indicated that NH4 + released from whale feces enhance productivity, as would be expected, with no evidence that fecal metabolites suppress growth. Although North Atlantic right whales are currently rare (approximately 450 individuals), they once numbered about 14,000 and likely played a substantial role in recycling nutrients in areas where they gathered to feed and mate. Even though the NH4 + released from fresh whale fecal material is a small fraction of total whale fecal nitrogen, and recognizing the fact that the additional nitrogen released in whale urine would be difficult to measure in a field study, the results of this study support the idea that the distinctive isotopic signature of the released NH4 + could be used to provide a conservative estimate of the contribution of the whale pump to primary productivity in coastal regions where whales congregate

Inferring introduction routes of invasive species using approximate Bayesian computation on microsatellite data

Determining the routes of introduction provides not only information about the history of an invasion process, but also information about the origin and construction of the genetic composition of the invading population. It remains difficult, however, to infer introduction routes from molecular data because of a lack of appropriate methods. We evaluate here the use of an approximate Bayesian computation (ABC) method for estimating the probabilities of introduction routes of invasive populations based on microsatellite data. We considered the crucial case of a single source population from which two invasive populations originated either serially from a single introduction event or from two independent introduction events. Using simulated datasets, we found that the method gave correct inferences and was robust to many erroneous beliefs. The method was also more efficient than traditional methods based on raw values of statistics such as assignment likelihood or pairwise F(ST). We illustrate some of the features of our ABC method, using real microsatellite datasets obtained for invasive populations of the western corn rootworm, Diabrotica virgifera virgifera. Most computations were performed with the DIYABC program (http://www1.montpellier.inra.fr/CBGP/diyabc/)

Secondary contact and admixture between independently invading populations of the Western corn rootworm, diabrotica virgifera virgifera in Europe

The western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is one of the most destructive pests of corn in North America and is currently invading Europe. The two major invasive outbreaks of rootworm in Europe have occurred, in North-West Italy and in Central and South-Eastern Europe. These two outbreaks originated from independent introductions from North America. Secondary contact probably occurred in North Italy between these two outbreaks, in 2008. We used 13 microsatellite markers to conduct a population genetics study, to demonstrate that this geographic contact resulted in a zone of admixture in the Italian region of Veneto. We show that i) genetic variation is greater in the contact zone than in the parental outbreaks; ii) several signs of admixture were detected in some Venetian samples, in a Bayesian analysis of the population structure and in an approximate Bayesian computation analysis of historical scenarios and, finally, iii) allelic frequency clines were observed at microsatellite loci. The contact between the invasive outbreaks in North-West Italy and Central and South-Eastern Europe resulted in a zone of admixture, with particular characteristics. The evolutionary implications of the existence of a zone of admixture in Northern Italy and their possible impact on the invasion success of the western corn rootworm are discussed

On dynamical probabilities, or: how to learn to shoot straight

© IFIP International Federation for Information Processing 2016.In order to support, for example, a quantitative analysis of various algorithms, protocols etc. probabilistic features have been introduced into a number of programming languages and calculi. It is by now quite standard to define the formal semantics of (various) probabilistic languages, for example, in terms of Discrete Time Markov Chains (DTMCs). In most cases however the probabilities involved are represented by constants, i.e. one deals with static probabilities. In this paper we investigate a semantical framework which allows for changing, i.e. dynamic probabilities which is still based on time-homogenous DTMCs, i.e. the transition matrix representing the semantics of a program does not change over time

Monte Carlo simulation of expected outcomes with the AcrySof® toric intraocular lens

<p>Abstract</p> <p>Background</p> <p>To use a Monte Carlo simulation to predict postoperative results with the AcrySof^®Toric lens, evaluating the likelihood of over- or under-correction using various toric lens selection criteria.</p> <p>Methods</p> <p>Keratometric data were obtained from a large patient population with preoperative corneal astigmatism <= 2.50D (2,000 eyes). The probability distributions for toric marking accuracy, surgically induced astigmatism and lens rotation were estimated using available data. Anticipated residual astigmatism was calculated using a Monte Carlo simulation under two different lens selection scenarios.</p> <p>Results</p> <p>This simulation demonstrated that random errors in alignment, surgically induced astigmatism and lens rotation slightly reduced the overall effect of the toric lens. Residual astigmatism was statistically significantly higher under the simulation of surgery relative to an exact calculation (p < 0.05). The simulation also demonstrated that more aggressive lens selection criteria could produce clinically significant reductions in residual astigmatism in a high percentage of patients.</p> <p>Conclusion</p> <p>Monte Carlo simulation suggests that surgical variability and lens orientation/rotation variability may combine to produce small reductions in the correction achieved with the AcrySof^®Toric^®IOL. Adopting more aggressive lens selection criteria may yield significantly lower residual astigmatism values for many patients, with negligible overcorrections. Surgeons are encouraged to evaluate their AcrySof^®Toric^®outcomes to determine if they should modify their individual lens selection criteria, or their default surgically induced astigmatism value, to benefit their patients.</p

Developmental changes in the role of different metalinguistic awareness skills in Chinese reading acquisition from preschool to third grade

Copyright @ 2014 Wei et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.The present study investigated the relationship between Chinese reading skills and metalinguistic awareness skills such as phonological, morphological, and orthographic awareness for 101 Preschool, 94 Grade-1, 98 Grade-2, and 98 Grade-3 children from two primary schools in Mainland China. The aim of the study was to examine how each of these metalinguistic awareness skills would exert their influence on the success of reading in Chinese with age. The results showed that all three metalinguistic awareness skills significantly predicted reading success. It further revealed that orthographic awareness played a dominant role in the early stages of reading acquisition, and its influence decreased with age, while the opposite was true for the contribution of morphological awareness. The results were in stark contrast with studies in English, where phonological awareness is typically shown as the single most potent metalinguistic awareness factor in literacy acquisition. In order to account for the current data, a three-stage model of reading acquisition in Chinese is discussed.National Natural Science Foundation of China and Knowledge Innovation Program of the Chinese Academy of Sciences

Human kallikrein 6: a new potential serum biomarker for uterine serous papillary cancer.

PURPOSE: The discovery of novel biomarkers might greatly contribute to improve clinical management and outcomes in uterine serous papillary carcinoma (USPC), a highly aggressive variant of endometrial cancer. EXPERIMENTAL DESIGN: Human kallikrein 6 (hK6) gene expression levels were evaluated in 29 snap-frozen endometrial biopsies, including 13 USPC, 13 endometrioid carcinomas, and 3 normal endometrial cells by real-time PCR. Secretion of hK6 protein by 14 tumor cultures, including 3 USPC, 3 endometrioid carcinoma, 5 ovarian serous papillary carcinoma, and 3 cervical cancers, was measured using a sensitive ELISA. Finally, hK6 concentration in 79 serum and plasma samples from 22 healthy women, 20 women with benign diseases, 20 women with endometrioid carcinoma, and 17 USPC patients was studied. RESULTS: hK6 gene expression levels were significantly higher in USPC when compared with endometrioid carcinoma (mean copy number by real-time PCR, 1,927 versus 239, USPC versus endometrioid carcinoma; P < 0.01). In vitro hK6 secretion was detected in all primary USPC cell lines tested (mean, 11.5 microg/L) and the secretion levels were similar to those found in primary ovarian serous papillary carcinoma cultures (mean, 9.6 microg/L). In contrast, no hK6 secretion was detectable in primary endometrioid carcinoma and cervical cancer cultures. hK6 serum and plasma concentrations (mean +/- SE) among normal healthy females (2.7 +/- 0.2 microg/L), patients with benign diseases (2.4 +/- 0.2 microg/L), and patients with endometrioid carcinoma (2.6 +/- 0.2 microg/L) were not significantly different. In contrast, serum and plasma hK6 values in USPC patients (6.1 +/- 1.1) were significantly higher than those in the noncancer group (P = 0.006), benign group (P = 0.003), and endometrioid carcinoma patients (P = 0.005). CONCLUSIONS: hK6 is highly expressed in USPC and is released in the plasma and serum of USPC patients. hK6 may represent a novel biomarker for USPC for monitoring early disease recurrence and response to therapy

Induction of tumor-specific cytotoxicity in tumor infiltrating lymphocytes by HPV16 and HPV18 E7-pulsed autologous dendritic cells in patients with cancer of the uterine cervix.

OBJECTIVE: To evaluate the potential of autologous dendritic cells (DC) pulsed with HPV16 and HPV18 E7 oncoprotein in restoring tumor-specific cytotoxicity in populations of tumor infiltrating lymphocytes (TIL) for adoptive immunotherapy of cervical cancer patients. METHODS: Full-length E7-pulsed DC-stimulated CD8(+) T cells derived from peripheral blood (PBL) and from tumor tissues (TIL) were tested and compared for their ability to induce a HLA class-I-restricted cytotoxic T lymphocyte (CTL) response against autologous tumor cells. In addition, in order to correlate cytotoxic activity by CTL with a particular lymphoid subset, analysis of surface antigens and intracellular CD3 zeta chain and two-color flow cytometric analysis of intracellular cytokine expression (IFN-gamma vs IL-4) at the single cell level were performed. RESULTS: DC stimulation induced powerful cytotoxicity against autologous tumor target cells by TIL-derived CD8(+) T cells from all three cervical cancer patients, while autologous Epstein-Barr virus-transformed lymphoblastoid cell lines were not lysed. Killing of autologous tumor cells was higher by CD8(+) T cells from TIL compared to PBL (P &gt; 0.01) and was more strongly inhibited by anti-HLA class I MAb (P &gt; 0.05). Phenotypically, all CTL populations were CD3(+)/CD8(+), with higher levels of CD56 expression by TIL-derived CTL. Finally, although a marked Type 1 cytokine bias (i.e., IFN-gamma(high)/IL-4(low)) was observable in both PBL- and TIL-derived DC-stimulated CD8(+) T cell populations, TIL-derived CD8(+) T cells showed a higher percentage of IFN-gamma-positive cells compared to PBL. CONCLUSIONS: Full-length E7-pulsed DC can consistently restore strong CD8(+) CTL responses against autologous HPV16- and HPV18-infected cervical cancer cells. DC-stimulated TIL may represent a superior source of tumor-specific CTL compared to PBL for adoptive T cell immunotherapy of patients harboring metastatic or recurrent cervical cancer refractory to standard treatment modalities