Search NASA⌕ Search

Engineering topics

Martin, William R.

Publications and source records attributed to Martin, William R..

Code Verification and Solution Verification framework in pin-resolved neutron transport code MPACT

Program verification in scientific computing encompasses the application of formal and mathematical techniques to a scientific computing code for its credibility, accuracy, and validity. Code Verification identifies bugs and performance issues in the software development stage. Solution Verification assesses the applicability of the code and the accuracy of the solution to problems of interest. Both activities utilize application cases and quantify the error against prescribed acceptance criteria. However, simply executing more application cases does not guarantee stronger or more comprehensive credibility. Here, we establish a verification framework that involves Code Verification and Solution Verification, both of which work together such that the overarching goal of “converge to the correct answer for the intended application” can be reasonably inferred. The application of such a verification framework is demonstrated using the pin-resolved neutron transport code MPACT, where standard unit tests and regression tests are covered, and where the Method of Exact Solutions and the Method of Manufactured Solutions are successfully used. Additionally, the applicability of Method of Manufactured Solutions is extended to the OECD/NEA C5G7 benchmark problems of practical material and geometric configurations. Solution Verification activities are demonstrated on a practical hierarchy of application models of increasing complexity ranging from 2D pin cell problems to 3D assembly problems. The convergence behavior and rate of convergence with respect to each individual variable are studied and provided. This framework can be adapted broadly to other fields involving scientific computing codes.

97 MATHEMATICS AND COMPUTING↗

Improved MPACT energy deposition and explicit heat generation coupling with CTF

The CASL neutronics code MPACT assumes local energy deposition with equilibrium delayed energy for steady-state and transient calculations. These approximations limit the accurate representation of the heat generation in space and its variations over time, both of which are essential for power distribution and thermal–hydraulic coupling analyses. This paper presents an explicit energy deposition model considering the neutron and gamma heating of all regions and the time-dependent delayed energy. The MPACT/CTF interfaces are updated to transfer the explicit heat sources. The new model agrees well with MCNP for problems without thermal–hydraulic feedback. The MPACT/CTF coupled calculation for a hot-full-power quarter-core case using the new model exhibits a reduction of peak pin power by 2.3% and fuel centerline temperature by 17 K. The new model also shows a keff difference up to -100 pcm in depletion and a 30 K peak fuel temperature reduction in a reactivity-initiated-accident problem.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Comprehensive characterization of protein–protein interactions perturbed by disease mutations

Technological and computational advances in genomics and interactomics have made it possible to identify how disease mutations perturb protein–protein interaction (PPI) networks within human cells. Here, we show that disease-associated germline variants are significantly enriched in sequences encoding PPI interfaces compared to variants identified in healthy participants from the projects 1000 Genomes and ExAC. Somatic missense mutations are also significantly enriched in PPI interfaces compared to noninterfaces in 10,861 tumor exomes. We computationally identified 470 putative oncoPPIs in a pan-cancer analysis and demonstrate that oncoPPIs are highly correlated with patient survival and drug resistance/sensitivity. Further, we experimentally validate the network effects of 13 oncoPPIs using a systematic binary interaction assay, and also demonstrate the functional consequences of two of these on tumor cell growth. In summary, this human interactome network framework provides a powerful tool for prioritization of alleles with PPI-perturbing mutations to inform pathobiological mechanism- and genotype-based therapeutic discovery.

59 BASIC BIOLOGICAL SCIENCES↗

Multigroup Constant Calculation with Static α-Eigenvalue Monte Carlo for Time-Dependent Neutron Transport Simulations

Herein, multigroup constants for deterministic methods that preserve the time-dependent physics of the neutron transport equations are derived. Alternative multigroup constant weighting spectra are discussed: (1) the fundamental k-eigenfunction, (2) the fundamental α-eigenfunction, and (3) a composite of several α-modes. To generate the fundamental α-eigenfunction for calculating the multigroup constants, a static fundamental α-eigenvalue method is implemented into the open source Monte Carlo code OpenMC. Several static and kinetic problems are devised to verify the implementations and to investigate the relative performance of the alternative multigroup constant weighting spectra. Results emphasize that as a multigroup constant weighting spectrum, the fundamental α-eigenfunction offers physical characteristics that make it advantageous (in producing accurate solutions) over the typically used fundamental k-eigenfunction.

42 ENGINEERING↗

An Improved Energy Deposition Model in MPACT and Explicit Heat Generation Coupling with CTF

The default energy deposition model in the CASL neutronics code MPACT assumes all fission energy is deposited locally in fuel rods. Furthermore, equilibrium delayed energy release is assumed for both steady-state and transient calculations. These approximations limit the accurate representation of the heat generation distribution in space and its variations over time, which are essential for power distribution and thermal-hydraulic coupling calculations. In this paper, an improved energy deposition model is presented in both the spatial and time domains. Spatially, the energy deposition through fission, neutron capture, and slowing-down reactions are explicitly modeled to account for the heat generation from all regions of a reactor core, and a gamma smearing scheme is developed that utilizes the gamma sources from neutron fission and capture. In the time domain, the delayed energy release is modeled by solving an additional equation of delayed heat emitters, similar to the equation of delayed neutron precursors.To allow the explicit heat generation coupling, the interfaces between MPACT and CTF were updated to transfer separate heat sources for different material regions (fuel, clad, moderator and guide tube). The results show that the distributions of the energy deposition between MPACT and MCNP agree very well for various 2-D assembly and quarter-core problems without TH feedback. The MPACT/CTF coupled calculation for the hot full power quarter-core case exhibited a reduced peak pin power by 2.3% and a reduced peak fuel centerline temperature by 17 K when using the explicit energy deposition and heat transfer. The new model also shows a maximum 100 pcm keff effect on assembly depletion problems and an increased overall energy release by 7% in a PWR reactivity-initiated accident (RIA) problem.

Liu, Yuxuan↗

In Silico Insights into Protein–Protein Interaction Disruptive Mutations in the PCSK9-LDLR Complex

Gain-of-function mutations in PCSK9 (proprotein convertase subtilisin/kexin type 9) lead to reduced uptake of LDL (low density lipoprotein) cholesterol and, therefore, increased plasma LDL levels. However, the mechanism by which these mutants reduce LDL reuptake is not fully understood. Here, we have used molecular dynamics simulations, MM/PBSA (Molecular Mechanics/Poisson–Boltzmann Surface Area) binding affinity calculations, and residue interaction networks, to investigate the protein–protein interaction (PPI) disruptive effects of two of PCSK9's gain-of-function mutations, Ser127Arg and Asp374Tyr on the PCSK9 and LDL receptor complex. In addition to these PPI disruptive mutants, a third, non-interface mutation (Arg496Trp) is included as a positive control. Our results indicate that Ser127Arg and Asp374Tyr confer significantly improved binding affinity, as well as different binding modes, when compared to the wild-type. These PPI disruptive mutations lie between the EGF(A) (epidermal growth factor precursor homology domain A) of the LDL receptor and the catalytic domain of PCSK9 (Asp374Tyr) and between the prodomain of PCSK9 and the β-propeller of the LDL receptor (Ser127Arg). The interactions involved in these two interfaces result in an LDL receptor that is sterically inhibited from entering its closed conformation. This could potentially implicate the prodomain as a target for small molecule inhibitors.

59 BASIC BIOLOGICAL SCIENCES↗

α-weighted Transition Rate Matrix Method

Relative effectiveness of an α-weighted Transition Rate Matrix Method (TRMM) as compared to the normally used κ- weighted, for solving reactor transient problem is investigated. An improved Monte Carlo κ-α iteration with delayed neutron is implemented to the OpenMC to generate the fundamental α-eigenfunction to be used as the weighting function in calculating the TRMM forward and adjoint operators. A prompt supercritical transient problem and a deep subcritical transient problem of a two-region slab system are devised to compare the relative performance of the α- and κ-weighted TRMM.

Variansyah, Ilham↗

A Metaheuristic Optimization Tool for High Flux Isotope Reactor Low-Enriched Uranium Core Design

A metaheuristic optimization tool was developed for design optimization studies of a High Flux Isotope Reactor low-enriched uranium core. Low-enriched uranium core design presents a multidimensional, multi-objective optimization problem that requires a systematic approach to avoid arbitrary design decisions. This tool searches this multidimensional design space to generate an optimum design, considering a weighted design performance metric that includes several parameters. The search methods are discussed herein, and the tool is applied to uranium silicide fuel designs.

Variansyah, I.↗

α-Weighted Transition Rate Matrix Method

A more direct approach in calculating the Transition Rate Matrix Method operators and the significance of this method are presented. Consequences of using α-spectrum as the Transition Rate Matrix Method weighting function instead of the κ-spectrum are investigated. A prompt supercritical infinite medium multigroup system is considered. Results on continuous problem using OpenMC will be provided in the full paper.

Variansyah, I.↗

Implementation of Gamma Transport Capability in MPACT

This paper presents the latest work on the development of the gamma transport capability for the CASL neutronics code MPACT. MPACT is capable of performing pin-resolved neutronics calculations for full-core problems and the purpose of this project is to further expand its capability for gamma transport calculations at the same level of accuracy for potential applications of interests, such as gamma heating and gamma fluence calculations. The development of the gamma transport capability reuses many existing geometrical capabilities and neutronics solvers in MPACT, because essentially the same form of Boltzmann transport equation is solved. First, 2-D gamma transport is implemented in this work by leveraging the method of characteristics (MOC) sweeper. The HELIOS library is used to provide the neutron-induced gamma production and gamma interaction cross sections. The new gamma transport capability is verified with MCNP. Both the gamma spectrum of a pin cell problem and the gamma flux distribution of 2D assembly problems are consistent with reference calculations. A 2D quarter core case is qualitatively verified by comparing MCNP and MPACT predictions of the gamma flux and thermal neutron flux to confirm the gamma smearing effects.

Wang, Xinyan↗

Ultra-precision processes for optics manufacturing

The Optics MODIL (Manufacturing Operations Development and Integration Laboratory) is developing advanced manufacturing technologies for fabrication of ultra precision optical components, aiming for a ten-fold improvement in precision and a shortening of the scheduled lead time. Current work focuses on diamond single point turning, ductile grinding, ion milling, and in/on process metrology.

Martin, William R.↗

Probabilistic structural mechanics research for parallel processing computers

Aerospace structures and spacecraft are a complex assemblage of structural components that are subjected to a variety of complex, cyclic, and transient loading conditions. Significant modeling uncertainties are present in these structures, in addition to the inherent randomness of material properties and loads. To properly account for these uncertainties in evaluating and assessing the reliability of these components and structures, probabilistic structural mechanics (PSM) procedures must be used. Much research has focused on basic theory development and the development of approximate analytic solution methods in random vibrations and structural reliability. Practical application of PSM methods was hampered by their computationally intense nature. Solution of PSM problems requires repeated analyses of structures that are often large, and exhibit nonlinear and/or dynamic response behavior. These methods are all inherently parallel and ideally suited to implementation on parallel processing computers. New hardware architectures and innovative control software and solution methodologies are needed to make solution of large scale PSM problems practical.

Sues, Robert H.↗