Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational methods”

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 595 records · Page 33

Simplifying computational workflows with the Multiscale Atomic Zeolite Simulation Environment (MAZE)

Zeolites, an important class of 3-dimensional nanoporous materials, have been widely explored for a variety of applications including gas storage, separations, and catalysis. As the properties of these aluminosilicate materials depend on a number of factors (e.g., framework topology, Si/Al ratio, extra-framework cations etc.), detailed experiments (e.g., catalytic properties, adsorption capacities etc.) are often limited to only a handful of materials. Computational methods have played an important role in (1) providing molecular level insights to rationalize experimental observations, and (2) screening large libraries of zeolites to identify promising candidates for experimental synthesis and validation. Different levels of theory and computational chemistry codes are necessary to describe the range of relevant phenomena such as adsorption (e.g., grand canonical Monte Carlo), diffusion (e.g., molecular dynamics), and chemical reactions (e.g., density functional theory). Manipulation of atomic structures, handling of input files, and developing robust workflows becomes quite cumbersome. To mitigate these challenges, we describe the development of the Multiscale Atomic Zeolite Simulation Environment (MAZE) – a Python package that simplifies zeolite-specific calculation workflows by providing a user-friendly interface for systematically manipulating zeolite structures

97 MATHEMATICS AND COMPUTING↗

Fibre Channel Development Final Report CRADA No. TC-0340-92-D

The Fibre Channel (FC) Development project was designed to enhance U.S. industrial competitiveness by accelerating the standardization and subsequent industrial acceptance of the American National Standards Institute’s (AIW) Fibre Channel Standard (now X3T1 1). It was also intended to give the principal U.S. participants a one- to two-year lead in developing heterogeneous high-speed computer networks, dktibuted computing through clustering, scientific visualization techniques, and high speed channels. The project was proposed to benefit U.S. industry as these high performance computing methods became generally available. In 1992-93, during the initial project phases, Fibre Channel PC) was emerging as a new ANSI standard. It was envisioned as becoming the mechanism for data transfer between the next generation of computers, and between these computers and their peripherals. Using this channel standard, U.S. computer vendors were expecting to be able to provide high-speed interoperability between heterogeneous computers and their storage devices, with data transfer rates greater than a GBit/sec. Our project was composed of key principal players in the FC standardization effort from both U.S. industry and the Department of Energy’s (DOES) Lawrence Liverrnore National Laboratory (LLNL). The industial participants were developing individual FC interfaces to operate with their platform products. LLNL’s role was to be a neutral testing site and a facilitator for helping to resolve conflicting interpretations of the embryonic Fibre Channel specification. LLNL anticipated that conflicts in the interpretation of such a standard would arise and require neutral resolution and that such resolutions would best be achieved in cooperation with the ANSI FC Working Group. Project plans included transferring to U.S. industry any resilient FC testing capabilities at the project conclusion.

42 ENGINEERING↗

Fibre Channel Development Final Report CRADA No. TC-0340-92-C

The Fibre Channel (FC) Development project was designed to enhance U.S. industrial competitiveness by accelerating the standardization and subsequent industrial acceptance of the American National Standards Institute’s (ANSI) Fibre Channel Standard (now X3T1 1). It was also intended to give the principal U.S. participants a one-to two-year lead in developing heterogeneous high-speed computer networks, distributed computing through clustering, scientific visualization techniques, and high speed channels. The project was proposed to benefit U.S. industry as these high performance computing methods became generally available. In 1992-93, during the initial project phases, Fibre Channel (FC) was emerging as a new ANSI standard. It was envisioned as becoming the mechanism for data transfer between the next generation of computers, and between these computers and their peripherals. Using this channel standard, U.S. computer vendors were expecting to be able to provide high-speed interoperatillity between heterogeneous computers and their storage devices, with data transfer rates greater than a GBit/sec. Our project was composed of key principal players in the FC standardization effort from both U.S. industry and the Department of Energy’s@OEs) Lawrence Liverrnore National Laboratory (LLNL). The industrial participants were developing individual FC interfaces to operate with their platform products. LLNL’s.role was to be a neutral testing site and a facilitator for helping to resolve conflicting interpretations of the embryonic Fibre Channel specification. LLNL anticipated that conflicts in the interpretation of such a standard would arise and require neutral resolution and that such resolutions would best be achieved in cooperation with the ANSI FC Working Group. Project plans included transferring to U.S. industry any resilient FC testing capabilities at the project conclusion.

42 ENGINEERING↗

EBS Task Force: Task 9/FEBEX Modeling Final Report: Thermo-Hydrological Modeling with PFLOTRAN

This report outlines Sandia National Laboratories modeling studies applied to Stage 1 and Stage 2 of the Full-scale Engineered Barriers Experiment in Crystalline Host Rock (FEBEX) in situ test for the SKB EBS Task Force Task 9. The FEBEX test was a full-scale test conducted over ~18 years at the Grimsel, Switzerland Underground Research Laboratory (URL) managed by NAGRA. It involved emplacing simulated waste packages, in the form of welded cylindrical heaters, inside a tunnel in crystalline granitic rock and surrounded by a bentonite barrier and cement plug. Sensors emplaced within the bentonite monitored the wetting-up, heating, and drying out of the bentonite barrier, and the large resulting data set provides an excellent opportunity for validation of multiphysics Thermal-Hydrological (TH), Thermal-Hydrologic-Chemical (THC), and Thermal-Hydrological-Mechanical (THM) modeling approaches for underground nuclear waste storage and the performance of engineered bentonite barriers. The present status of the EBS Task Force is finalizing Task 9, which follows years of modeling studies of the FEBEX test, by many notable modeling teams (Gens et al., 2009; Sanchez et al. 2010; 2012; Samper et al., 2018). These modeling studies generally use two-dimensional axisymmetric meshes, ignoring threedimensional effects, gravity and asymmetric wetting and dry out of the bentonite engineered barrier. This study investigates these effects with use of the PFLOTRAN THC code with massively parallel computational methods in modeling FEBEX Stage 1 and Stage 2 results. The PFLOTRAN numerical code is an open source, state-of-the-art, massively parallel subsurface flow and reactive transport code operating in a high-performance computing environment (Hammond et al., 2014). Section 2 describes the applied partial differential equations describing mass, momentum and energy balance used in this study, considerations derived by assuming phase equilibrium between gas and liquid phases, constitutive equations for granite, cement plug, and bentonite domains, and specific approaches for use inthe PFLOTRAN code. Section 3 describes the geometry, meshing, and model set-up. Section 4 describes modeling results, Section 5 compares modeling results to field testing data, and Section 6 gives conclusions. The Appendix provides detailed information required by the EBSTask Force for final reporting.

42 ENGINEERING↗

Grid Optimization Competition Challenge 3 Problem Formulation

This report contains the problem formulation for the Grid Optimization (GO) Competition Challenge 3. The Grid Optimization Competition is run by a team of researchers from a number of organizations, including the sponsor Advanced Research Projects Agency - Energy (ARPA-E), lead organization Pacific Northwest National Laboratory (PNNL), and technical contributors from Los Alamos National Laboratory (LANL), National Renewable Energy Laboratory (NREL), Texas A&M University (TAMU), Georgia Institute of Technology (GT), University of Wisconsin (UW), and others. The GO Competition poses challenge problems in the field of power grid management, invites entrants to develop solvers for these problems, invokes the solvers on a set of problem instances using common hardware, ranks the solvers according to their performance, and awards prizes according to the rankings. The overall goal of the GO Competition is to spur innovative research on high impact and computationally challenging problems in power grid management from initial development through commercial deployment. Complete information about the GO Competition can be found online at [2]. The webpage covers previous Challenges, rules, timeline, registration information, data formats, scoring methods, computational platform information, information on supported solvers and languages, sponsor information, frequently asked questions, administrator contact information, publicly available problem instances, computer code for reading and evaluating problem and solution data, a sandbox for testing solvers, and a solver submission interface, results, and publications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Natural equilibria in steady-state neutron diffusion with temperature feedback

The critical diffusion equation with feedback is investigated within the context of steady-state multiphysics. It is proposed that for critical configurations there is no need to include the multiplication factor k in the formulation of the diffusion equation. This is notable because exclusion of k from the coupled system of equations precludes the mathematically tenuous notion of a nonlinear eigenvalue problem. On the other hand, it is shown that if the factor k is retained in the diffusion equation, as is currently common practice, then the resulting problem is equivalent to the constrained minimization of a functional representing the critical equilibrium of neutron and temperature distributions. The unconstrained solution corresponding to k = 1 represents the natural equilibrium of a critical system at steady-state. Computational methods for solving the constrained problem (with k) are briefly reviewed from the literature and a method for the unconstrained problem (without k) is outlined. A numerical example is studied to examine the effects of the constraint in the nonlinear system. (authors)

97 MATHEMATICS AND COMPUTING↗

Likelihood-informed dimension reduction for nonlinear inverse problems

The intrinsic dimensionality of an inverse problem is affected by prior information, the accuracy and number of observations, and the smoothing properties of the forward operator. From a Bayesian perspective, changes from the prior to the posterior may, in many problems, be confined to a relatively low-dimensional subspace of the parameter space. We present a dimension reduction approach that defines and identifies such a subspace, called the 'likelihood-informed subspace' (LIS), by characterizing the relative influences of the prior and the likelihood over the support of the posterior distribution. This identification enables new and more efficient computational methods for Bayesian inference with nonlinear forward models and Gaussian priors. In particular, we approximate the posterior distribution as the product of a lower-dimensional posterior defined on the LIS and the prior distribution marginalized onto the complementary subspace. Markov chain Monte Carlo sampling can then proceed in lower dimensions, with significant gains in computational efficiency. We also introduce a Rao–Blackwellization strategy that de-randomizes Monte Carlo estimates of posterior expectations for additional variance reduction. We demonstrate the efficiency of our methods using two numerical examples: inference of permeability in a groundwater system governed by an elliptic PDE, and an atmospheric remote sensing problem based on Global Ozone Monitoring System (GOMOS) observations.

42 ENGINEERING↗

Overcoming the Memory Bottleneck in Auxiliary Field Quantum Monte Carlo Simulations with Interpolative Separable Density Fitting

We investigate the use of interpolative separable density fitting (ISDF) as a means to reduce the memory bottleneck in auxiliary field quantum Monte Carlo (AFQMC) simulations of real materials in Gaussian basis sets. Here, we find that ISDF can reduce the memory scaling of AFQMC simulations from O(M 4 ) to O(M 2 ) . We test these developments by computing the structural properties of carbon in the diamond phase, comparing to results from existing computational methods and experiment. (LLNL-JRNL-759024)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Method and computer-readable model for additively manufacturing ducting arrangement for a gas turbine engine

Method and computer-readable model for additively manufacturing a ducting arrangement (10) for a gas turbine engine are provided. Ducting arrangement (10) may include a duct (18) to be fluidly coupled to receive a cross-flow of combustion gases from a main combustion stage. Duct (18) includes a duct segment (23) with an expanding cross-sectional area (24) where one or more injector assemblies (26) are disposed. Injector assembly (26) includes one or more reactant-guiding structures (27) arranged to deliver a flow of reactants to be mixed with the cross-flow of combustion gases. The ducting arrangement is effective to reduce total pressure loss while providing an effective level of mixing of the injected reactants with the passing cross-flow. Respective duct components or the entire ducting arrangement may be formed as a unitized structure, such as a single piece using a rapid manufacturing technology, such as 3D Printing/Additive Manufacturing (AM) technologies.

Chen, Yuntao↗

Silver diamine fluoride differentially affects dentin and hypomineralized enamel permeabilities

OBJECTIVES: To investigate the physicochemical effect of silver diamine fluoride (SDF) by correlating permeability with mineral density and elemental composition of hypomineralized enamel and carious dentin. METHODS: Enamel and dentin from human carious primary teeth with and without SDF treatment in-vivo, and hypomineralized enamel from permanent molars with and without SDF treatment in-vitro were scanned using micro X-ray computed tomography. Spatial maps of biometals (calcium, zinc), phosphorus, and silver were generated using X-ray fluorescence microprobe. Permeabilities were computed using Porous Microstructure Analysis software. RESULTS: The intrinsic permeability of SDF-treated carious dentin was 14.3 % lower than untreated sound dentin (6.39e-15 ± 3.01e-15 m² vs 7.46e-15 ± 1.82e-15 m²; P < 0.0001), while untreated carious dentin was 98.4 % higher (1.48e-14 ± 7.11e-15 m²; P < 0.0001). SDF-treated and untreated transparent dentin showed similar reduced permeabilities (75.6 % and 78.4 % lower than untreated sound dentin, respectively; P = 0.93). Severely hypomineralized enamel showed permeability reaching 108.1 % of adjacent sound dentin (5.71e-15 ± 2.04e-15 m² vs 5.28e-15 ± 1.30e-15 m²; P = 0.1409) and was significantly higher than mildly hypomineralized enamel (1.39e-15 ± 1.04e-15 m²; P < 0.0001). SDF treatment did not significantly impact the permeability of severely hypomineralized enamel (12.4 % reduction; P = 0.07). Principal component regression identified Zn level as a significant effector of tissue permeabilities in carious primary teeth (P < 0.0001). SIGNIFICANCE: This study introduces a computational method to measure dental tissue permeability, and demonstrates that SDF significantly reduces permeability in carious dentin but not intact hypomineralized enamel. The study reveals biometal Zn localization can alter dentin and enamel permeabilities, providing new insights into pathobiological mechanisms underlying caries and hypomineralization.

Chou, Conrad↗

Hydrodynamic Analysis of a Suspended Cylinder Under Regular Wave Loading Based on Computational Fluid Dynamics

An investigation into the computation of hydrodynamic loads on a suspended cylinder in regular waves is presented. The primary goal was to perform a three-way validation of the loads between experimental measurements and simulations from two computational methods. Experimental measurements of the longitudinal in-line force on a cylinder suspended at a fixed position were available from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project, Phase Ia. These measurements were compared to computational fluid dynamics (CFD) simulations based on the solution of Reynolds-averaged Navier-Stokes (RANS) equations, as implemented in STAR-CCM+. The study encompassed a sensitivity analysis of the loads computed in STAR-CCM+ based on wave modeling, boundary conditions, turbulence modeling, and spatial and temporal discretization. The analysis was supplemented by results generated with the offshore wind turbine engineering software OpenFAST, based on a hybrid combination of second-order potential flow and viscous drag from Morison's equation. The focus of the investigation was on the assessment of the accuracy of the computation of first- and higher-order hydrodynamic loads. Substantial differences were observed in the numerical prediction of the second and third harmonic force contribution. Local flow field analysis with CFD was applied to study the physics of wave run-up and diffraction dynamics to identify the causes.

16 TIDAL AND WAVE POWER↗

Hydrodynamic Analysis of a Suspended Cylinder Under Regular Wave Loading Based on Computational Fluid Dynamics: Preprint

An investigation into the computation of hydrodynamic loads on a suspended cylinder in regular waves is presented. The primary goal was to perform a three-way-validation of the loads between experimental measurements and simulations from two computational methods. Experimental measurements of the longitudinal inline force on a cylinder suspended at a fixed po- sition were available from the Offshore Code Comparison Col- laboration, Continued, with Correlation (OC5) project, Phase Ia. These measurements were compared to Computational Fluid Dynamics (CFD) simulations based on the solution of Reynolds- averaged Navier-Stokes (RANS) equations, as implemented in STAR-CCM plus. The study encompassed a sensitivity analysis of the loads computed in STAR-CCM plus based on wave modeling, boundary conditions, turbulence modeling, and spatial and tem- poral discretization. The analysis was supplemented by results generated with the offshore wind turbine engineering software OpenFAST, based on a hybrid combination of second-order po- tential flow and viscous drag fromMorison's equation. The focus of the investigation was on the assessment of the accuracy of the computation of first- and higher-order hydrodynamic loads. Sub- stantial differences were observed in the numerical prediction of the second and third harmonic force contribution. Local flow field analysis with CFD was applied to study the physics of wave run-up and diffraction dynamics to identify the causes.

16 TIDAL AND WAVE POWER↗

Understanding the influence of capillary waves on solvation at the liquid-vapor interface

This work investigates the question if surface capillary waves (CWs) affect interfacial solvation thermodynamic properties that determine the propensity of small molecules toward the liquid-vapor interface. We focus on (1) the evaluation of these properties from molecular simulations in a practical manner and (2) understanding them from the perspective of theories in solvation thermodynamics, especially solvent reorganization effects. Concerning the former objective, we propose a computational method that exploits the relationship between an external field acting on the liquid-vapor interface and the magnitude of CWs. The system considered contains the solvent, an externally applied field (f) and the solute molecule fixed at a particular location. The magnitude of f is selected to induce changes in CWs. The difference between the solvation free energies computed in the presence and in the absence of f is then shown to quantify the contribution of CWs to interfacial solvation. We describe the implementation of this method in the canonical ensemble by using a Lennard-Jones solvent and a non-ionic solute. Results are shown for three types of solutes that differ in the nature of short-ranged repulsive (hard-core) interactions. Overall, we observe that CWs have a negligible or very small effect on the interfacial solvation free energy of a solute molecule fixed near the liquid-vapor interface for the above systems. We also explain how the effects of pinning or dampening of CWs caused by a fixed solute are effectively compensated and do not contribute to the solvation free energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Long non-coding subgenomic flavivirus RNAs have extended 3D structures and are flexible in solution

Most mosquito-borne flaviviruses, including Zika virus (ZIKV), Dengue virus (DENV), and West Nile virus (WNV), produce long non-coding subgenomic RNAs (sfRNAs) in infected cells that link to pathogenicity and immune evasion. Until now, the structural characterization of these lncRNAs remains limited. Here, we studied the 3D structures of individual and combined subdomains of sfRNAs, and visualized the accessible 3D conformational spaces of complete sfRNAs from DENV2, ZIKV, and WNV by small angle X-ray scattering (SAXS) and computational modeling. The individual xrRNA1s and xrRNA2s adopt similar structures in solution as the crystal structure of ZIKV xrRNA1, and all xrRNA1-2s form compact structures with reduced flexibility. While the DB12 of DENV2 is extended, the DB12s of ZIKV and WNV are compact due to the formation of intertwined double pseudoknots. All 3' stem-loops (3'SLs) share similar rod-like structures. Complete sfRNAs are extended and sample a large conformational space in solution. Our work not only provides structural insight into the function of flavivirus sfRNAs, but also highlights strategies of visualizing other lncRNAs in solution by SAXS and computational methods.

3D structure↗

On Xenon Fission Product Poisoning

We all appreciate the special role Xenon-135 plays in reactor transients. Unless there is sufficient reactivity override, on shutdown a conventional reactor cannot immediately restart because of the buildup of strongly absorbing Xenon from Iodine decay and the absence of neutron flux to continue its elimination. We are also aware of how the flux imbalance in a large reactor core gives rise to Xenon induced power precession around the core in unpredictable ways. In addition, with regard to computational methods, we lack readily accessible numerical algorithms enabling reliable characterization of Xenon oscillations with simple models to explain this fundamentally important phenomenon. Nowhere is this more apparent than in the classroom. In this presentation, we first consider a new 12-element comprehensive decay to model Xenon and its corresponding numerical implementation and verification to establish a predictive model of Xenon buildup in a constant flux. We then follow with a reduced self-consistent flux model in an infinite reactor- a model not found in textbooks. This presentation should appeal to those interested in developing high order numerical methods for reactor transients.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetic structure of UO 2 and NpO 2 by first-principle methods

Here, the magnetic structure of the actinide dioxides (AnO 2 ) remains a field of intense research. A low-temperature experimental investigation of the magnetic ground-state is complicated by thermal energy released from the radioactive decay of the actinide nuclei. To establish the magnetic ground-state, we have employed high-accuracy computational methods to systematically probe different magnetic structures. A transverse 1k antiferromagnetic ground-state with Fmmm (No. 69) crystal symmetry has been established for UO 2 , whereas a ferromagnetic (111) ground-state with R¯3 with combining macron]m (No. 166) has been established for NpO 2 . Band structure calculations have been performed to analyse these results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗

NPIP: A skew line needle configuration optimization system for HDR brachytherapy

Purpose: In this study, the authors introduce skew line needle configurations for high dose rate (HDR) brachytherapy and needle planning by integer program (NPIP), a computational method for generating these configurations. NPIP generates needle configurations that are specific to the anatomy of the patient, avoid critical structures near the penile bulb and other healthy structures, and avoid needle collisions inside the body. Methods: NPIP consisted of three major components: a method for generating a set of candidate needles, a needle selection component that chose a candidate needle subset to be inserted, and a dose planner for verifying that the final needle configuration could meet dose objectives. NPIP was used to compute needle configurations for prostate cancer data sets from patients previously treated at our clinic. NPIP took two user-parameters: a number of candidate needles, and needle coverage radius, {delta}. The candidate needle set consisted of 5000 needles, and a range of {delta} values was used to compute different needle configurations for each patient. Dose plans were computed for each needle configuration. The number of needles generated and dosimetry were analyzed and compared to the physician implant. Results: NPIP computed at least one needle configuration for every patient that met dose objectives, avoided healthy structures and needle collisions, and used as many or fewer needles than standard practice. These needle configurations corresponded to a narrow range of {delta} values, which could be used as default values if this system is used in practice. The average end-to-end runtime for this implementation of NPIP was 286 s, but there was a wide variation from case to case. Conclusions: The authors have shown that NPIP can automatically generate skew line needle configurations with the aforementioned properties, and that given the correct input parameters, NPIP can generate needle configurations which meet dose objectives and use as many or fewer needles than the current HDR brachytherapy workflow. Combined with robot assisted brachytherapy, this system has the potential to reduce side effects associated with treatment. A physical trial should be done to test the implant feasibility of NPIP needle configurations.

60 APPLIED LIFE SCIENCES↗