Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallelization Strategy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

An mRNA vaccine encoding the Ebola virus glycoprotein induces high neutralizing antibody titers and provides strong protection against lethal infections in mouse models

Ebola virus (EBOV) is the causative agent of Ebola disease (EBOD), a viral hemorrhagic fever with a notably high case fatality rate. Current treatments for EBOD are limited to monoclonal antibodies or two licensed viral vector vaccines, a recombinant vesicular stomatitis virus (rVSV)-vectored vaccine or an adenovirus and modified vaccinia Ankara regimen. However, comparisons of protection, efficacy, and durability with alternative nucleotide platforms remain understudied. Here, we evaluated the immunogenicity of an mRNA vaccine expressing the EBOV glycoprotein (GP) in parallel with rVSV- and DNA-based vaccine platforms. The mRNA EBOV-GP vaccine, formulated in lipid nanoparticles, elicited significantly higher levels of total IgG and neutralizing antibody titers compared to the rVSV-EBOV-GP vaccine. Linear antibody epitope analysis indicated a preference for targeting the mucin-like domain in EBOV-GP1 following rVSV-based vaccination, while the mRNA platform distinctly targeted the internal fusion loop of EBOV-GP2. After characterizing the immunogenicity of the mRNA vaccine, two models of EBOD were used to demonstrate its protective efficacy: a surrogate rVSV-based challenge model of EBOD using type-I interferon deficient C57BL/6 mice and infection of BALB/c mice with authentic mouse-adapted EBOV. In both studies, the EBOV mRNA vaccine fully protected the mouse cohorts against morbidity and mortality. Additionally, the EBOV mRNA vaccine produced greater neutralizing antibody titers compared to the DNA EBOV-GP vaccine. These results suggest that an mRNA vaccine expressing EBOV-GP can induce robust, functional humoral responses that are protective against EBOD, warranting further development as an alternative to, or as part of a vaccine strategy including, viral vectored vaccines.

DNA vaccines↗

Instrumentation and Control Needs for Reliable Operation of Lunar Base Surface Nuclear Power Systems

As one of the near-term goals of the President's Vision for Space Exploration, establishment of a multi-person lunar base will require high-endurance power systems which are independent of the sun, and can operate without replenishment for several years. These requirements may be obtained using nuclear power systems specifically designed for use on the lunar surface. While it is envisioned that such a system will generally be supervised by humans, some of the evolutions required maybe semi or fully autonomous. The entire base complement for near-term missions may be less than 10 individuals, most or all of which may not be qualified nuclear plant operators and may be off-base for extended periods thus, the need for power system autonomous operation. Startup, shutdown, and load following operations will require the application of advanced control and health management strategies with an emphasis on robust, supervisory, coordinated control of, for example, the nuclear heat source, energy conversion plant (e.g., Brayton Energy Conversion units), and power management system. Autonomous operation implies that, in addition to being capable of automatic response to disturbance input or load changes, the system is also capable of assessing the status of the integrated plant, determining the risk associated with the possible actions, and making a decision as to the action that optimizes system performance while minimizing risk to the mission. Adapting the control to deviations from design conditions and degradation due to component failures will be essential to ensure base inhabitant safety and mission success. Intelligent decisions will have to be made to choose the right set of sensors to provide the data needed to do condition monitoring and fault detection and isolation because of liftoff weight and space limitations, it will not be possible to have an extensive set of instruments as used for earth-based systems. Advanced instrumentation and control technologies will be needed to enable this critical functionality of autonomous operation. It will be imperative to consider instrumentation and control requirements in parallel to system configuration development so as to identify control-related, as well as integrated system-related, problem areas early to avoid potentially expensive work-arounds . This paper presents an overview of the enabling technologies necessary for the development of reliable, autonomous lunar base nuclear power systems with an emphasis on system architectures and off-the-shelf algorithms rather than hardware. Autonomy needs are presented in the context of a hypothetical lunar base nuclear power system. The scenarios and applications presented are hypothetical in nature, based on information from open-literature sources, and only intended to provoke thought and provide motivation for the use of autonomous, intelligent control and diagnostics.

Turso, James↗

Site-specific regulation of oral mucosa-recruiting CD8{sup +} T cells in a mouse contact allergy model

Contact allergy is a T cell-mediated, delayed-type hypersensitivity generated by contact exposure of an allergen to the skin and mucosal surface. The clinical manifestations of allergic responses between the skin and oral mucosa vary and the differences in immunopathology have not been clarified. We generated hapten-induced contact hypersensitivity (CH) of the buccal mucosa (BM) in parallel studies with ear skin (ES) CH, and observed several characteristic findings of BM CH. The BM challenge induced more rapid and more severe inflammation than the ES challenge, with abundant granulocyte and CD8{sup +} T cell infiltration. However, these inflammatory responses diminished quickly. Recruiting CD8{sup +} T cells in the BM had higher ratios of CD62L{sup −}CD44{sup low-hi} memory-type cells, and showed impaired IFN-γ, greater PD-1, and comparable Ki-67 expression, suggesting that the recruiting-proliferating CD8{sup +} T cells were unable to differentiate into effector T cells and converted into exhausted T cells at the local site. This finding may explain the rapid recovery of the BM from severe inflammation. Preferentially greater expression of PD-1 ligand (B7-H1), was observed in the BM epithelium under the peak inflammation, and the absence of B7-H1 further accelerated CH responses, suggesting the occurrence of PD-1:B7-H1-mediated immune regulation at the local site. Our results may facilitate the understanding of the unique features of contact allergies in the oral mucosa, and guide the development of new strategies for control of contact allergy. - Highlights: • Contact allergy in the oral mucosa shows severe responses, but quick recovery. • Oral mucosa-infiltrating CD8{sup +} T cells have unique activatory and regulatory status. • PD-1:B7-H1-mediated pathway is involved in the site-specific regulation.

60 APPLIED LIFE SCIENCES↗

Planning for Pre-Exascale Platform Environment (Fiscal Year 2015 Level 2 Milestone 5216)

This Plan for ASC Pre-Exascale Platform Environments document constitutes the deliverable for the fiscal year 2015 (FY15) Advanced Simulation and Computing (ASC) Program Level 2 milestone Planning for Pre-Exascale Platform Environment. It acknowledges and quantifies challenges and recognized gaps for moving the ASC Program towards effective use of exascale platforms and recommends strategies to address these gaps. This document also presents an update to the concerns, strategies, and plans presented in the FY08 predecessor document that dealt with the upcoming (at the time) petascale high performance computing (HPC) platforms. With the looming push towards exascale systems, a review of the earlier document was appropriate in light of the myriad architectural choices currently under consideration. The ASC Program believes the platforms to be fielded in the 2020s will be fundamentally different systems that stress ASC’s ability to modify codes to take full advantage of new or unique features. In addition, the scale of components will increase the difficulty of maintaining an errorfree system, thus driving new approaches to resilience and error detection/correction. The code revamps of the past, from serial- to vector-centric code to distributed memory to threaded implementations, will be revisited as codes adapt to a new message passing interface (MPI) plus “x” or more advanced and dynamic programming models based on architectural specifics. Development efforts are already underway in some cases, and more difficult or uncertain aspects of the new architectures will require research and analysis that may inform future directions for program choices. In addition, the potential diversity of system architectures may require parallel if not duplicative efforts to analyze and modify environments, codes, subsystems, libraries, debugging tools, and performance analysis techniques as well as exploring new monitoring methodologies. It is difficult if not impossible to selectively eliminate some of these activities until more information is available through simulations of potential architectures, analysis of systems designs, and informed study of commodity technologies that will be the constituent parts of future platforms.

97 MATHEMATICS AND COMPUTING↗

Nonrecursive formulations of multibody dynamics and concurrent multiprocessing

Since the late 1980's, research in recursive formulations of multibody dynamics has flourished. Historically, much of this research can be traced to applications of low dimensionality in mechanism and vehicle dynamics. Indeed, there is little doubt that recursive order N methods are the method of choice for this class of systems. This approach has the advantage that a minimal number of coordinates are utilized, parallelism can be induced for certain system topologies, and the method is of order N computational cost for systems of N rigid bodies. Despite the fact that many authors have dismissed redundant coordinate formulations as being of order N(exp 3), and hence less attractive than recursive formulations, we present recent research that demonstrates that at least three distinct classes of redundant, nonrecursive multibody formulations consistently achieve order N computational cost for systems of rigid and/or flexible bodies. These formulations are as follows: (1) the preconditioned range space formulation; (2) penalty methods; and (3) augmented Lagrangian methods for nonlinear multibody dynamics. The first method can be traced to its foundation in equality constrained quadratic optimization, while the last two methods have been studied extensively in the context of coercive variational boundary value problems in computational mechanics. Until recently, however, they have not been investigated in the context of multibody simulation, and present theoretical questions unique to nonlinear dynamics. All of these nonrecursive methods have additional advantages with respect to recursive order N methods: (1) the formalisms retain the highly desirable order N computational cost; (2) the techniques are amenable to concurrent simulation strategies; (3) the approaches do not depend upon system topology to induce concurrency; and (4) the methods can be derived to balance the computational load automatically on concurrent multiprocessors. In addition to the presentation of the fundamental formulations, this paper presents new theoretical results regarding the rate of convergence of order N constraint stabilization schemes associated with the newly introduced class of methods.

Kurdila, Andrew J.↗

A Probabilistic Framework for the Validation and Certification of Computer Simulations

The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.

Ghanem, Roger↗

Efficient data IO for a Parallel Global Cloud Resolving Model

Execution of a Global Cloud Resolving Model (GCRM) at target resolutions of 2-4 km will generate, at a minimum, 10s of Gigabytes of data per variable per snapshot. Writing this data to disk without creating a serious bottleneck in the execution of the GCRM code while also supporting efficient post-execution data analysis is a significant challenge. This paper discusses an Input/Output (IO) application programmer interface (API) for the GCRM that efficiently moves data from the model to disk while maintaining support for community standard formats, avoiding the creation of very large numbers of files, and supporting efficient analysis. Several aspects of the API will be discussed in detail. First, we discuss the output data layout which linearizes the data in a consistent way that is independent of the number of processors used to run the simulation and provides a convenient format for subsequent analyses of the data. Second, we discuss the flexible API interface that enables modelers to easily add variables to the output stream by specifying where in the GCRM code these variables are located and to flexibly configure the choice of outputs and distribution of data across files. The flexibility of the API is designed to allow model developers to add new data fields to the output as the model develops and new physics is added and also provides a mechanism for allowing users of the GCRM code itself to adjust the output frequency and the number of fields written depending on the needs of individual calculations. Third, we describe the mapping to the NetCDF data model with an emphasis on the grid description. Fourth, we describe our messaging algorithms and IO aggregation strategies that are used to achieve high bandwidth while simultaneously writing concurrently from many processors to shared files. We conclude with initial performance results.

99 GENERAL AND MISCELLANEOUS↗

Influence of polarization and a source model for dose calculation in MRT

Purpose: Microbeam Radiation Therapy (MRT), an alternative preclinical treatment strategy using spatially modulated synchrotron radiation on a micrometer scale, has the great potential to cure malignant tumors (e.g., brain tumors) while having low side effects on normal tissue. Dose measurement and calculation in MRT is challenging because of the spatial accuracy required and the arising high dose differences. Dose calculation with Monte Carlo simulations is time consuming and their accuracy is still a matter of debate. In particular, the influence of photon polarization has been discussed in the literature. Moreover, it is controversial whether a complete knowledge of phase space trajectories, i.e., the simulation of the machine from the wiggler to the collimator, is necessary in order to accurately calculate the dose. Methods: With Monte Carlo simulations in the Geant4 toolkit, the authors investigate the influence of polarization on the dose distribution and the therapeutically important peak to valley dose ratios (PVDRs). Furthermore, the authors analyze in detail phase space information provided byMartínez-Rovira et al. [“Development and commissioning of a Monte Carlo photon model for the forthcoming clinical trials in microbeam radiation therapy,” Med. Phys. 39(1), 119–131 (2012)] and examine its influence on peak and valley doses. A simple source model is developed using parallel beams and its applicability is shown in a semiadjoint Monte Carlo simulation. Results are compared to measurements and previously published data. Results: Polarization has a significant influence on the scattered dose outside the microbeam field. In the radiation field, however, dose and PVDRs deduced from calculations without polarization and with polarization differ by less than 3%. The authors show that the key consequences from the phase space information for dose calculations are inhomogeneous primary photon flux, partial absorption due to inclined beam incidence outside the field center, increased beam width and center to center distance due to the beam propagation from the collimator to the phantom surface and imperfect absorption in the absorber material of the Multislit Collimator. These corrections have an effect of approximately 10% on the valley dose and suffice to describe doses in MRT within the measurement uncertainties of currently available dosimetry techniques. Conclusions: The source for the first clinical pet trials in MRT is characterized with respect to its phase space and the photon polarization. The results suggest the use of a presented simplified phase space model in dose calculations and hence pave the way for alternative and fast dose calculation algorithms. They also show that the polarization is of minor importance for the clinical important peak and valley doses inside the microbeam field.

61 RADIATION PROTECTION AND DOSIMETRY↗

Mixed Integer Nonlinear Programming Approaches to Enhance Resiliency and Response Strategies in Critical Infrastructure

Motivation. Critical infrastructures are large, complex engineered systems that must be operated robustly under abnormal conditions resulting from natural hazards or intentional acts. For example, electric power systems must be robust to line faults, water utilities must rapidly mitigate contamination incidents, and computing networks must adapt to adversarial intrusions to protect critical information. Problem. It is difficult for decision-makers within resiliency analysis in critical infrastructure to optimize designs and develop effective response strategies that can account for uncertainties. Facing incomplete information and the sheer scope that a natural hazard or attack vector may incorporate, response can be ineffective without reliable, scalable decision support tools. These problems are intrinsically nonlinear and involve discrete decisions, and unfortunately, existing off-the-shelf mathematical programming methods cannot support optimization-based decision-making of these nonlinear at scale. Method/Approach and Results. This project emphasized development of fundamental optimization strategies that supported real- time mitigation and response for critical infrastructures. In particular, the project developed multi- tree approaches based on piecewise outer-approximations for solution of mixed-integer nonlinear programming (MINLP) problems. These techniques alternate between an MILP or MISOCP relaxation to obtain a lower bound and candidate discrete solutions and an NLP subproblem to obtain upper bounds. Using tailored relaxations based on problem structure, these methods were used to solve several key applications in resilience and response of critical infrastructure. This work resulted in two open-source, copyrighted software packages: CORAMIN (https://github.com/Coramin/Coramin) -- an object-oriented mathematical programming framework that supports tailored multi-tree algorithms for solution of large- scale mixed-integer nonlinear programming; and EGRET (Electrical Grid Research and Engineering Toolkit) (https://github.com/grid- parity-exchange/Egret) -- a declarative mathematical programming framework built upon CORAMIN and Pyomo for formulation and solution of resilience and operations problems in power grid systems. Furthermore, these tools resulted in several important published results, including the following: The first known global optimization approach that could solve the unit-commitment problem with nonlinear power flow constraints on medium-sized test problems; Improved parallel optimization-based bounds tightening and strengthening of relaxations of AC power flow constraints; and, Optimization-based approaches for improved grid resilience and use of demand response to improve grid resilience with reduction in capital requirements. Result Implications. This project developed first-of-a-kind algorithms for decision-making in critical infrastructure resilience operations and planning, as well as a next-generation toolkit for MINLP researchers. These approaches leveraged high-performance computing architectures to solve some of the largest, most challenging nonlinear discrete optimization problems to global optimality, and these successes were captured in open-source software to enable optimization-based decision-making, and efficient solution of MINLP formulations for electric power transmission grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Digital Twin Model for Advanced Manufacture of a Rotating Detonation Engine Injector

A digital twin material model (DTMM) of an additive manufacturing (AM) process was created to advance the state of the art in rotating detonation engine (RDE) injector design. Current RDE injectors are designed with large pressure drops, enabling a stable and repeatable combustion process. However, this comes at the cost of system efficiency. For the technology to transition to commercial fossil-based power generation, it is important to develop injectors with reduced flow losses. Low-loss injectors are difficult to design and manufacture with conventional manufacturing techniques. AM enables new design options, but the AM manufacturing process must be thoroughly understood to result in a robust design. A DTMM provides the necessary insight by defining the cause-effect relationships between process parameters, microstructure features, and properties. Therefore, a DTMM to support the design and manufacturing process was developed and applied to the design of a new additively manufactured low-loss injector. The injector combustion behavior was characterized through hot-fire tests, and mechanical performance was compared to the DTMM predictions. The two project goals were the successful development of the DTMM and the demonstration of an improved RDE injector design. The RDE injector design and DTMM developments occurred on parallel but dependent paths. The injector was designed to reduce pressure drop by increasing the cross-sectional flow area ratio between the injector air passages and the combustor annulus. This resulted in less structural material, raising the concern that thin members would be susceptible to high-cycle fatigue (HCF) under the periodic loading inherent to an RDE. It was most important for the DTMM to predict behavior in these features; therefore, the injector design concept guided the material thicknesses used in fatigue tests. The DTMM development started by manufacturing a series of coupons over the range of possible AM process variations. A design-of-experiment approach was used to select which process variable combinations gave the most efficient coverage relevant to the injector design space. The microstructure in each of these coupons was characterized, and then computational methods were used to create a numerical model of the correlation between process variables and microstructure. Next, a set of HCF samples were tested to calibrate existing models that map microstructure to HCF performance. Together, these two links formed the DTMM that calculates HCF behavior from AM process variables. Two injector prototypes were additively manufactured. The first injector design strategy aggressively pursued low-loss performance by substantially increasing the oxidizer flow area. The combination of manufacturing lead times and the fatigue testing schedule meant that the DTMM was not available when building this first prototype. Therefore, its process parameters were chosen based on a manual review of the available coupon data. This prototype was built successfully and evaluated in 58 combustion tests. Sustained detonation was achieved with remarkably reduced pressure loss, and some tests even displayed pressure loss characteristics similar to conventional gas turbine combustors. This achieved the project goal of improving RDE injector design. The second injector was manufactured according to the optimized parameters predicted by the DTMM. The flow area modifications of this injector were less aggressive than the first injector since demonstrating low pressure loss was not an objective of the second hot-fire test series. Rather, the test objective was to cause high cycle fatigue failure in the part due to periodic loading from the rotating detonation wave. The observed number of cycles to failure was to be compared to the number predicted by the DTMM and thereby assess the utility of the DTMM in component design. However, the required level of vibration was not obtained during combustion. Therefore, high cycle fatigue was not experienced in the hot-fire tests of the second injector. Fatigue data was obtained by further testing the second injector in a conventional HCF test apparatus. The injector demonstrated HCF strength above the DTMM prediction. In fact, it did not fail and testing was only discontinued due to reaching the end of the period of performance. This points to some success in the project’s primary goal of successfully developing and applying the DTMM to a component design. Implementing the DTMM recommendations for optimal processing parameters led to a part with acceptable properties. The DTMM was also shown to be an efficient correlator of data and to provide insight into the relationship between process settings, microstructure, and property performance. However, the failure of the DTMM prediction to match the experimental result of the injector fatigue test also points to the need to include significantly more data in the model development. In this project, coupons made with identical processing parameters exhibited drastically different properties from each other and from the injector part, which clearly influences the accuracy of a model that predicts performance based on parameters. Uncertainties in the build process must be quantified to develop more robust models. A denser and broader matrix of coupon process and geometry variations, several repeated builds of every point, more in-situ build process measurements, and direct observation of tensile and HCF sample microstructure (as opposed to separate microstructure specimens) are recommendations to improve future AM modeling efforts.

20 FOSSIL-FUELED POWER PLANTS↗

CO₂ Capture Membrane Process for Power Plant Flue Gas

Because the fleet of coal-fired power plants is of such importance to the nation's energy production while also being the single largest emitter of CO₂, the development of retrofit, post-combustion CO₂ capture technologies for existing and new, upcoming coal power plants will allow coal to remain a major component of the U.S. energy mix while mitigating global warming. Post-combustion carbon capture technologies are an attractive option for coal-fired power plants as they do not require modification of major power-plant infrastructures, such as fuel processing, boiler, and steam-turbine subsystems. In this project, the overall objective was to develop an advanced, hollow-fiber, polymeric membrane process that could be cost-effectively retrofitted into current pulverized coal-fired power plants to capture at least 90% of the CO₂ from plant flue gas with 95% captured CO₂ purity. The approach for this project tackled the technology development on three different fronts in parallel: membrane materials R&D, hollow-fiber membrane module development, and process development and engineering. The project team consisted of RTI (prime) and two industrial partners, Arkema, Inc. and Generon IGS, Inc. Two CO₂-selective membrane polymer platforms were targeted for development in this project. For the near term, a next-generation, high-flux polycarbonate membrane platform was spun into hollow-fiber membranes that were fabricated into both lab-scale and larger prototype (~2,200 ft²) membrane modules. For the long term, a new fluoropolymer membrane platform based on poly(vinylidene fluoride) [PVDF] chemistry was developed using a copolymer approach as improved capture membrane materials with superior chemical resistance to flue-gas contaminants (moisture, SO₂, NO x , etc.). Specific objectives were: - Development of new, highly chemically resistant, fluorinated polymers as membrane materials with minimum selectivity of 30 for CO₂ over N₂ and CO₂ permeance greater than 300 gas permeation units (GPU) targeted; - Development of next-generation polycarbonate hollow-fiber membranes and membrane modules with higher CO₂ permeance than current commercial polycarbonate membranes; - Development and fabrication of membrane hollow fibers and modules from candidate polymers; - Development of a CO₂ capture membrane process design and integration strategy suitable for end-of-pipe, retrofit installation; and - Techno-economic evaluation of the "best" integrated CO₂ capture membrane process design package In this report, the results of the project research and development efforts are discussed and include the post-combustion capture properties of the two membrane material platforms and the hollow-fiber membrane modules developed from them and the multi-stage process design and analysis developed for 90% CO₂ capture with 95% captured CO₂ purity.

20 FOSSIL-FUELED POWER PLANTS↗

ELSI: A unified software interface for Kohn–Sham electronic structure solvers

Solving the electronic structure from a generalized or standard eigenproblem is often the bottleneck in large scale calculations based on Kohn-Sham density-functional theory. This problem must be addressed by essentially all current electronic structure codes, based on similar matrix expressions, and by high-performance computation. We here present a unified software interface, ELSI, to access different strategies that address the Kohn-Sham eigenvalue problem. Currently supported algorithms include the dense generalized eigensolver library ELPA, the orbital minimization method implemented in libOMM, and the pole expansion and selected inversion (PEXSI) approach with lower computational complexity for semilocal density functionals. The ELSI interface aims to simplify the implementation and optimal use of the different strategies, by offering (a) a unified software framework designed for the electronic structure solvers in Kohn-Sham density-functional theory; (b) reasonable default parameters for a chosen solver; (c) automatic conversion between input and internal working matrix formats, and in the future (d) recommendation of the optimal solver depending on the specific problem. As a result, comparative benchmarks are shown for system sizes up to 11,520 atoms (172,800 basis functions) on distributed memory supercomputing architectures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Combining Sparse Approximate Factorizations with Mixed-precision Iterative Refinement

The standard LU factorization-based solution process for linear systems can be enhanced in speed or accuracy by employing mixed-precision iterative refinement. Most recent work has focused on dense systems. We investigate the potential of mixed-precision iterative refinement to enhance methods for sparse systems based on approximate sparse factorizations. In doing so, we first develop a new error analysis for LU- and GMRES-based iterative refinement under a general model of LU factorization that accounts for the approximation methods typically used by modern sparse solvers, such as low-rank approximations or relaxed pivoting strategies. We then provide a detailed performance analysis of both the execution time and memory consumption of different algorithms, based on a selected set of iterative refinement variants and approximate sparse factorizations. Our performance study uses the multifrontal solver MUMPS, which can exploit block low-rank factorization and static pivoting. We evaluate the performance of the algorithms on large, sparse problems coming from a variety of real-life and industrial applications showing that mixed-precision iterative refinement combined with approximate sparse factorization can lead to considerable reductions of both the time and memory consumption.

97 MATHEMATICS AND COMPUTING↗

Managing Low-Level Waste at Dounreay - 16269

Dounreay Site Restoration Limited (DSRL) is the site licence company responsible for the clean-up and demolition of Britain's former centre of fast reactor research and development. The Dounreay site requires facilities to dispose of up to 175,000 cubic metres of low level waste (LLW) that has been and will continue to be generated during the decommissioning and restoration of the site. In May 2014, construction of two new vaults for the disposal of LLW at Dounreay was completed, including the necessary roads, infrastructures and services, followed by an encapsulation plant in September 2014. The vaults are located immediately to the northeast of the Dounreay licensed site, on NDA-owned land. The facilities are for disposal of the solid LLW generated at Dounreay and the adjacent Vulcan site only. In 2004, following consultation with stakeholders and members of the public, Dounreay issued the 'Best Practicable Environmental Option (BPEO) Report' on the management of Dounreay LLW. This report assessed potential management options for LLW against a range of criteria, including technical, environmental, cost, health and safety issues. The study was conducted in line with best practice and involved consultation with a range of people, including the local community and organisations with an interest in radioactive waste management. The recommendations and conclusions from the study provided the basis for the 'Dounreay Solid LLW Overall Strategy', which was published in March 2005. A fundamental component of this strategy was the development of new below-ground disposal vaults for LLW at Dounreay. Siting the vaults on land at Dounreay avoids any need to transport the LLW away from Dounreay on public roads. Disposal at Dounreay therefore satisfies the proximity principle of managing the waste at source. Disposal at Dounreay is also consistent with UK Government Policy on LLW management, the Scottish Government's Radioactive Waste Management Policy and the Nuclear Decommissioning Authority's (NDA's) strategy on management of LLW from the UK nuclear industry. A planning application was initially lodged with the Highland Council in June 2006. The period for the determination of this application was extended to allow the Scottish Environment Protection Agency (SEPA) sufficient time to assess fully the supporting Environmental Safety Case (ESC). The planning application was re-activated in May 2008, and planning permission was granted by Highland Council in April 2009. The planning permission included 26 conditions on DSRL, covering a wide range of environmental issues such as construction noise/dust, traffic, visual impact and also a condition to develop a community fund aimed at socio-economic development of the area around the Dounreay site. After undertaking the necessary detailed design work, construction of the phase one vaults commenced in November 2011 by principal contractor, Graham Construction. Work was undertaken on both vaults simultaneously, excavating around 243,000 m{sup 3} of rock. When excavation was complete, steel reinforcement was placed and concrete poured to construct the floor slabs. This was followed by the construction of the reinforced walls and then the installation of the steel support structure and materials that form the roof. Mechanical and electrical installation then followed, in parallel with the construction of the vault apron, ramps and access roads. A pumping system and control building was also constructed to remove groundwater from around the vaults to allow the vaults to be kept dry while they are being filled with waste. A layer of excavated material, approximately 3 m thick, was placed across the existing ground surface between the vaults and the coast, to ensure the water table lies well below the new ground surface after closure of the vaults, and that any future trace contaminated groundwater flows out to the sea rather than to the soil zone. This layer is referred to in project documents as the 'enhanced geosphere'. After disposals are completed and the vaults are closed, this layer will be joined to the capping over the vaults so that the ensemble merges visually with the surrounding landscape. The original vegetation and soil mix will be returned as far as is practicable. An encapsulation plant has been constructed to fill the LLW containers with grout prior to disposal in the vault, similar to the process used at the national Low-Level Waste Repository facility (LLWR), near Drigg in Cumbria. The plant has the capacity to grout up to four containers per day. The grout fills voids in the container and helps provide chemical conditions in the vault that limit the release of radionuclides in the far future, after the vault has been penetrated by groundwater. On April 24 2015, active commissioning began at the Encapsulation Plant and LLW Vaults. Half-height ISO containers packed with low level waste are currently being filled with grout and disposed of in the LLW Vault. Also packaged Demolition Low Level Waste (DLLW) is being processed and disposed of in the Demolition Low Level Waste Vault. To date over 2500 packages of DLLW and 70 containers LLW have been disposed of. Completion of construction and the start of operating the Phase One vaults and encapsulation plant are fundamental to the successful and efficient decommissioning of the Dounreay site, upon which many millions of pounds of investment and hundreds of jobs rely.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Dynamic Modeling, Trajectory Optimization, and Linear Control of Cable-Driven Parallel Robots for Automated Panelized Building Retrofits

The construction industry faces a growing need for automation to reduce costs, improve accuracy and productivity, and address labor shortages. One area that stands to benefit significantly from automation is panelized prefabricated building envelope retrofits, which can improve a building’s energy efficiency in heating and cooling interior spaces. In this paper, we propose using cable-driven parallel robots (CDPRs), which can effectively lift and handle large objects, to install these panels. However, implementing CDPRs presents significant challenges because of their nonlinear dynamics, complex trajectory planning, and precise control requirements. To tackle these challenges, this work focuses on a new application of established control and trajectory optimization theories in a CDPR simulation of a building envelope retrofit under real-world conditions. We first model the dynamics of CDPRs, highlighting the critical role of damping in system behavior. Building on this dynamic model, we formulate a trajectory optimization problem to generate feasible and efficient motion plans for the robot under operational and environmental constraints. Given the high precision required in the construction industry, accurately tracking the optimized trajectory is essential. However, challenges such as partial observability and external vibrations complicate this task. To address these issues, a Linear Quadratic Gaussian control framework is applied, enabling the robot to track the optimized trajectories with precision. Simulation results show that the proposed controller enables precise end effector positioning with errors under 4 mm, even in the presence of external wind disturbances. Through comprehensive simulations, our approach allows for an in-depth exploration of the system’s nonlinear dynamics, trajectory optimization, and control strategies under controlled yet highly realistic conditions. The results demonstrate the feasibility of CDPRs for automating panel installation and provide insights into their practical deployment.

CDPR↗

β1 integrin mediates an alternative survival pathway in breast cancer cells resistant to lapatinib

The overexpression of human epidermal growth factor receptor (HER)-2 in 20% of human breast cancers and its association with aggressive growth has led to widespread use of HER2-targeted therapies, such as trastuzumab (T) and lapatinib (L). Despite the success of these drugs, their efficacy is limited in patients whose tumors demonstrate de novo or acquired resistance to treatment. The b1 integrin resides on the membrane of the breast cancer cell, activating several elements of breast tumor progression including proliferation and survival. We developed a panel of HER2-overexpressing cell lines resistant to L, T, and the potent LT combination through long-term exposure and validated these models in 3D culture. Parental and L/T/LT-resistant cells were subject to HER2 and b1 integrin inhibitors in 3D and monitored for 12 days, followed by quantification of colony number. Parallel experiments were conducted where cells were either stained for Ki-67 and Terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) or harvested for protein and analyzed by immunoblot. Results were subjected to statistical testing using analysis of variance and linear contrasts, followed by adjustment with the Sidak method. Using multiple cell lines including BT474 and HCC1954, we reveal that in L and LT resistance, where phosphorylation of EGFR/HER1, HER2, and HER3 are strongly inhibited, kinases downstream of b1 integrin– including focal adhesion kinase (FAK) and Src–are up-regulated. Blockade of b1 by the antibody AIIB2 abrogates this up-regulation and functionally achieves significant growth inhibition of L and LT resistant cells in 3D, without dramatically affecting the parental cells. SiRNA against b1 as well as pharmacologic inhibition of FAK achieve the same growth inhibitory effect. In contrast, trastuzumab-resistant cells, which retain high levels of phosphorylated EGFR/HER1, HER2, and HER3, are only modestly growth-inhibited by AIIB2. Our data suggest that HER2 activity, which is suppressed in resistance involving L but not T alone, dictates whether b1 mediates an alternative pathway driving resistance. Our findings justify clinical studies investigating the inhibition of b1 or its downstream signaling moieties as strategies to overcome acquired L and LT resistance.

59 BASIC BIOLOGICAL SCIENCES↗

Direct health-care cost of head and neck cancers: a population-based study in north-eastern Italy

Improvements in prognosis of head-and-neck squamous cell carcinoma (HNSCC) have paralleled with an increase in health-care costs, so that an economic evaluation is of growing importance. Presently, most of the evidence is from insurance-based studies in the USA. Between 2007 and 2010, 879 HNSCC patients were identified through the population-based cancer registry of the Friuli Venezia Giulia region, including 266 oral, 187 oropharyngeal, 136 hypopharyngeal, and 290 laryngeal cancers. Health-care costs from diagnosis to treatment initiation and in the following 2 years were retrieved through a record linkage with the regional health data warehouse. This database collected comprehensive health information on all resident citizens. Generalized linear models with a gamma distribution and log-link function were applied to model costs. The average health-care cost from diagnosis up to 2 years after treatment initiation was €20,184 (95% confidence interval: €19,634 − 20,733). Heterogeneity emerged according to cancer site, elective treatment, and retreatment for cancer persistence/recurrence (no: €13,896; yes: €24,599; p < 0.001). An advanced stage was associated with increased costs stage (I: €12,969; II: €18,276; III: €26,229; IV: €25,574; p < 0.001) as the result of treatment complexity and elevated frequency of patients retreatment due to recurrence. These findings further support strategies to diagnose patients at an earlier cancer stage and the accurate definition of diagnostic and treatment pathways, to start treating patients when radical unimodal approach is still feasible. Besides the advantage in prognosis due to timely curative treatments, this would reduce the economic burden of cancer treatment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Monofluorination of Naphthyls Promotes the Cofacial π–π Stacking and Increases the Electron Mobility of Non-Planar Zinc(II) Complexes of Di(naphthylethynyl)azadipyrromethene

The homoleptic zinc(II) complex of [2,8-di(1-naphthylethynyl) 3,7-diphenyl 1,9-(4-hexylphenyl)azadipyrromethene (ZnL2) 2 ] is a promising non-planar non-fullerene acceptor for organic photovoltaic applications, but it has a relatively low electron mobility that may limit its performance. Here, we explored the fluorination of peripheral aryl groups to increase intermolecular cofacial π–π stacking interactions, which are desirable for electron transport. Here complexes with fluorine on the distal phenyls [Zn(1F-L2) 2 ], on the naphthyls [Zn(2F-L2) 2 ], and on both [Zn(3F-L2) 2 ] were synthesized and characterized. All three complexes had similar optical and electrochemical properties. The crystal packing structure of Zn(2F-L2) 2 and Zn(3F-L2) 2 revealed cofacial parallel-displaced π–π stacking between the fluorinated 1-naphthylethynyl groups. Such a cofacial orientation was not observed in Zn(L2) 2 crystals, suggesting that fluorination of the naphthyl groups promotes the cofacial π–π stacking orientation. The hole mobility increased from 1.0 × 10 –4 cm 2 V –1 s –1 for Zn(L2) 2 to 0.8–1.0 × 10 –3 cm 2 V –1 s –1 for the fluorinated complexes. Fluorination on the naphthyl groups increased the electron mobility from 4.2 × 10 –5 cm 2 V –1 s –1 for Zn(L2) 2 and Zn(1F-L2) 2 to 2.0 × 10 –4 cm 2 V –1 s –1 for Zn(2F-L2) 2 and Zn(3F-L2) 2 , consistent with cofacial π–π stacking being favorable for electron transport. The three complexes were tested in OPVs using regioregular poly(3-hexylthiophene) (P3HT) as the p-type material, and the best power conversion efficiencies were 5.2, 5.4, and 5.8% for Zn(2F-L2) 2 , Zn(1F-L2) 2 , and Zn(3F-L2) 2 , respectively, compared to 5.5% for Zn(L2) 2 . The fluorination combination found in Zn(3F-L2) 2 resulted in the best device performance. This study points to a viable strategy to increase the electron mobility and performance of non-planar zinc(II) complexes of azadipyrromethene.

36 MATERIALS SCIENCE↗