Automating embedded analysis capabilities and managing software complexity in multiphysics simulation part I: template-based generic programming

An approach for incorporating embedded simulation and analysis capabilities in complex simulation codes through template-based generic programming is presented. This approach relies on templating and operator overloading within the C++ language to transform a given calculation into one that can compute a variety of additional quantities that are necessary for many state-of-the-art simulation and analysis algorithms. An approach for incorporating these ideas into complex simulation codes through general graph-based assembly is also presented. These ideas have been implemented within a set of packages in the Trilinos framework and are demonstrated on a simple problem from chemical engineering

Automating embedded analysis capabilities and managing software complexity in multiphysics simulation part II: application to partial differential equations

A template-based generic programming approach was presented in a previous paper that separates the development effort of programming a physical model from that of computing additional quantities, such as derivatives, needed for embedded analysis algorithms. In this paper, we describe the implementation details for using the template-based generic programming approach for simulation and analysis of partial differential equations (PDEs). We detail several of the hurdles that we have encountered, and some of the software infrastructure developed to overcome them. We end with a demonstration where we present shape optimization and uncertainty quantification results for a 3D PDE application

Telomere maintenance is dependent on activities required for end repair of double-strand breaks

AbstractTelomeres are functionally distinct from ends generated by chromosome breakage, in that telomeres, unlike double-strand breaks, are insulated from recombination with other chromosomal termini [1]. We report that the Ku heterodimer and the Rad50/Mre11/Xrs2 complex, both of which are required for repair of double-strand breaks [2–5], have separate roles in normal telomere maintenance in yeast. Using epistasis analysis, we show that the Ku end-binding complex defined a third telomere-associated activity, required in parallel with telomerase [6] and Cdc13, a protein binding the single-strand portion of telomere DNA [7,8]. Furthermore, loss of Ku function altered the expression of telomere-located genes, indicative of a disruption of telomeric chromatin. These data suggest that the Ku complex and the Cdc13 protein function as terminus-binding factors, contributing distinct roles in chromosome end protection. In contrast, MRE11 and RAD50 were required for the telomerase-mediated pathway, rather than for telomeric end protection; we propose that this complex functions to prepare DNA ends for telomerase to replicate. These results suggest that as a part of normal telomere maintenance, telomeres are identified as double-strand breaks, with additional mechanisms required to prevent telomere recombination. Ku, Cdc13 and telomerase define three epistasis groups required in parallel for telomere maintenance

Albany: Using Component-based Design to Develop a Flexible, Generic Multiphysics Analysis Code

Abstract: Albany is a multiphysics code constructed by assembling a set of reusable, general components. It is an implicit, unstructured grid finite element code that hosts a set of advanced features that are readily combined within a single analysis run. Albany uses template-based generic programming methods to provide extensibility and flexibility; it employs a generic residual evaluation interface to support the easy addition and modification of physics. This interface is coupled to powerful automatic differentiation utilities that are used to implement efficient nonlinear solvers and preconditioners, and also to enable sensitivity analysis and embedded uncertainty quantification capabilities as part of the forward solve. The flexible application programming interfaces in Albany couple to two different adaptive mesh libraries; it internally employs generic integration machinery that supports tetrahedral, hexahedral, and hybrid meshes of user specified order. We present the overall design of Albany, and focus on the specifics of the integration of many of its advanced features. As Albany and the components that form it are openly available on the internet, it is our goal that the reader might find some of the design concepts useful in their own work. Albany results in a code that enables the rapid development of parallel, numerically efficient multiphysics software tools. In discussing the features and details of the integration of many of the components involved, we show the reader the wide variety of solution components that are available and what is possible when they are combined within a simulation capability. Key Words: partial differential equations, finite element analysis, template-based generic programmin

A mouse chromosome 4 balancer ENU-mutagenesis screen isolates eleven lethal lines.

BACKGROUND: ENU-mutagenesis is a powerful technique to identify genes regulating mammalian development. To functionally annotate the distal region of mouse chromosome 4, we performed an ENU-mutagenesis screen using a balancer chromosome targeted to this region of the genome. RESULTS: We isolated 11 lethal lines that map to the region of chromosome 4 between D4Mit117 and D4Mit281. These lines form 10 complementation groups. The majority of lines die during embryonic development between E5.5 and E12.5 and display defects in gastrulation, cardiac development, and craniofacial development. One line displayed postnatal lethality and neurological defects, including ataxia and seizures. CONCLUSION: These eleven mutants allow us to query gene function within the distal region of mouse chromosome 4 and demonstrate that new mouse models of mammalian developmental defects can easily and quickly be generated and mapped with the use of ENU-mutagenesis in combination with balancer chromosomes. The low number of mutations isolated in this screen compared with other balancer chromosome screens indicates that the functions of genes in different regions of the genome vary widely

Standardized Treatment of Active Tuberculosis in Patients with Previous Treatment and/or with Mono-resistance to Isoniazid: A Systematic Review and Meta-analysis

Performing a systematic review of studies evaluating retreatment of tuberculosis or treatment of isoniazid mono-resistant infection, Dick Menzies and colleagues find a paucity of evidence to support the WHO-recommended regimen

The DOE E3SM Coupled Model Version 1: Overview and Evaluation at Standard Resolution

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean and sea ice (60 km in the midlatitudes and 30 km at the equator and poles), and river transport (55 km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model’s strong aerosol-related effective radiative forcing (ERFari+aci = -1.65 W/m2) and high equilibrium climate sensitivity (ECS = 5.3 K).Plain Language SummaryThe U.S. Department of Energy funded the development of a new state-of-the-art Earth system model for research and applications relevant to its mission. The Energy Exascale Earth System Model version 1 (E3SMv1) consists of five interacting components for the global atmosphere, land surface, ocean, sea ice, and rivers. Three of these components (ocean, sea ice, and river) are new and have not been coupled into an Earth system model previously. The atmosphere and land surface components were created by extending existing components part of the Community Earth System Model, Version 1. E3SMv1’s capabilities are demonstrated by performing a set of standardized simulation experiments described by the Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima protocol at standard horizontal spatial resolution of approximately 1° latitude and longitude. The model reproduces global and regional climate features well compared to observations. Simulated warming between 1850 and 2015 matches observations, but the model is too cold by about 0.5 °C between 1960 and 1990 and later warms at a rate greater than observed. A thermodynamic analysis of the model’s response to greenhouse gas and aerosol radiative affects may explain the reasons for the discrepancy.Key PointsThis work documents E3SMv1, the first version of the U.S. DOE Energy Exascale Earth System ModelThe performance of E3SMv1 is documented with a set of standard CMIP6 DECK and historical simulations comprising nearly 3,000 yearsE3SMv1 has a high equilibrium climate sensitivity (5.3 K) and strong aerosol-related effective radiative forcing (-1.65 W/m2)Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/151288/1/jame20860_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/151288/2/jame20860.pd