Search NASA⌕ Search

SEARCH · Search NASA

Results for “discretization 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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Nonlinear discrete wavefront shaping for spatiotemporal pulse compression with multicore fibers

Here in this work, we apply an effective discrete phase front focusing method to control the coherent combining and temporal compression of laser pulses by means of nonlinear multicore fibers. We numerically demonstrate the possibility of combining almost entirely the input optical pulses injected in all cores, into the central core of a hexagonal lattice multicore fiber. We investigate the optimal operational conditions for the most effective pulse combining and analyze the influence of a positive input pulse temporal chirp and spatial phase modulation. We demonstrate that a pulse, which has been focused into an arbitrary core of the fiber, also undergoes temporal compression.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Transfer learning for metamaterial design and simulation

Abstract We demonstrate transfer learning as a tool to improve the efficacy of training deep learning models based on residual neural networks (ResNets). Specifically, we examine its use for study of multi-scale electrically large metasurface arrays under open boundary conditions in electromagnetic metamaterials. Our aim is to assess the efficiency of transfer learning across a range of problem domains that vary in their resemblance to the original base problem for which the ResNet model was initially trained. We use a quasi-analytical discrete dipole approximation (DDA) method to simulate electrically large metasurface arrays to obtain ground truth data for training and testing of our deep neural network. Our approach can save significant time for examining novel metasurface designs by harnessing the power of transfer learning, as it effectively mitigates the pervasive data bottleneck issue commonly encountered in deep learning. We demonstrate that for the best case when the transfer task is sufficiently similar to the target task, a new task can be effectively trained using only a few data points yet still achieve a test mean absolute relative error of 3 % with a pre-trained neural network, realizing data reduction by a factor of 1000.

Peng, Rixi↗
DOE OSTI2015

Proteus-SN user manual

This user manual describes how to set up a neutron transport simulation with the PROTEUS-SN code. A companion methodology manual describes the theory and algorithms within PROTEUS-SN.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

PROTEUS-SN User Manual

PROTEUS-SN is a three-dimensional, highly scalable, high-fidelity neutron transport code developed at Argonne National Laboratory. The code is applicable to all spectrum reactor transport calculations, particularly those in which a high degree of fidelity is needed either to represent spatial detail or to resolve solution gradients. PROTEUS-SN solves the second order formulation of the transport equation using the continuous Galerkin finite element method in space, the discrete ordinates approximation in angle, and the multigroup approximation in energy. PROTEUS-SN’s parallel methodology permits the efficient decomposition of the problem by both space and angle, permitting large problems to run efficiently on hundreds of thousands of cores. PROTEUS-SN can also be used in serial or on smaller compute clusters (10’s to 100’s of cores) for smaller homogenized problems, although it is generally more computationally expensive than traditional homogenized methodology codes. PROTEUS-SN has been used to model partially homogenized systems, where regions of interest are represented explicitly and other regions are homogenized to reduce the problem size and required computational resources. PROTEUS-SN solves forward and adjoint eigenvalue problems and permits both neutron upscattering and downscattering. An adiabatic kinetics option has recently been included for performing simple time-dependent calculations in addition to standard steady state calculations. PROTEUS-SN handles void and reflective boundary conditions. Multigroup cross sections can be generated externally using the MC2-3 fast reactor multigroup cross section generation code or internally using the cross section application programming interface (API) which can treat the subgroup or resonance table libraries. PROTEUS-SN is written in Fortran 90 and also includes C preprocessor definitions. The code links against the PETSc, METIS, HDF5, and MPICH libraries. It optionally links against the MOAB library and is a part of the SHARP multi-physics suite for coupled multi-physics analysis of nuclear reactors. This user manual describes how to set up a neutron transport simulation with the PROTEUS-SN code. A companion methodology manual describes the theory and algorithms within PROTEUS-SN.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

SCEPTRE 1.4 Quick Start Guide

This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general purpose C++ code for solving the Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates. Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Party Libraries (TPL) to be available, and example scripts for building these TPL’s are provided. The TPL’s needed by SCEPTRE are Trilinos, boost, and netcdf. SCEPTRE uses an autoconf build system, and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Tycho 2: A Proxy Application for Kinetic Transport Sweeps

Tycho 2 is a proxy application that implements discrete ordinates (SN) kinetic transport sweeps on unstructured, 3D, tetrahedral meshes. It has been designed to be small and require minimal dependencies to make collaboration and experimentation as easy as possible. Tycho 2 has been released as open source software. The software is currently in a beta release with plans for a stable release (version 1.0) before the end of the year. The code is parallelized via MPI across spatial cells and OpenMP across angles. Currently, several parallelization algorithms are implemented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SENSMG: First-Order Sensitivities of Neutron Reaction Rates, Reaction-Rate Ratios, Leakage, k eff , and α Using PARTISN

SENSMG is a tool for computing first-order sensitivities of neutron reaction rates, reaction-rate ratios, leakage, k eff , and α using the PARTISN multigroup discrete-ordinates code. SENSMG computes sensitivities to all of the transport cross sections and data (total, fission, nu, chi, and all scattering moments), two edit cross sections (absorption and capture), and the density for every isotope and energy group. It also computes sensitivities to the mass density for every material and derivatives with respect to all interface locations. The tool can be used for one-dimensional spherical (r) and two-dimensional cylindrical (r-z) geometries. The tool can be used for fixed-source and eigenvalue problems. The tool implements Generalized Perturbation Theory (GPT) as discussed by Williams and Stacey. Section II of this report describes the theory behind adjoint-based sensitivities, gives the equations that SENSMG solves, and defines the sensitivities that are output. Section III describes the user interface, including the input file and command line options. Section IV describes the output. Section V gives some notes about the coding that may be of interest. Section VI discusses verification, which is ongoing. Section VII lists needs and ideas for future work. Appendix A lists all of the input files whose results are presented in Sec. VI.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

MFIX-DEM Phi: Performance and Capability Improvements Towards Industrial Grade Open-source DEM Framework with Integrated Uncertainty Quantification

Two major challenges hinder the effective use and adoption of multiphase computational fluid dynamics tools by the industry. The first is the need for significant computational resources, which is inversely proportional to the accuracy of solutions due to computational intensity of the algorithms. The second barrier is assessing the prediction credibility and confidence in the simulation results. In this project, a multi-tiered approach has been proposed under four broad activities to overcome these challenges while addressing all of the objectives outlined in FOA-0001238 through Phases 1 and 2 of the project. The present report consists of the results for only Phase 1, which was the funded performance period. From the start the project, all of the objectives outlined in FOA were addressed through four major activity tasks in an integrated and balanced fashion to improve adoption of MFIX suite of solvers for industrial use. The first task aimed to improve the performance of MFIX-DEM specifically targeting to acquire the peak performance on Intel Xeon and Xeon Phi based systems, which are expected to be one of the primary high-performance computing platforms both affordable and available for the industrial users in the next two to five years. However, due to a number of changes in course of the project, the scope of the performance improvements related task was significantly reduced to avoid duplicate work. Hence, more emphasis was placed on the other three tasks as discussed below. The second task aimed at physical modeling enhancements through implementation of polydispersity capability and validation of heat transfer models in MFIX. An extended verification and validation (V&V) study was performed for the new polydispersity feature implemented in MFIX-DEM both for granular and coupled gas-solid flows. The features of the polydispersity capability and results for an industrially relevant problem were disseminated through journal papers (one published and one under review at the time of writing of the final technical report). As part of the validation efforts, another industrially relevant problem of interest based on rotary drums was studied for several modes of heat transfer and results were presented in conferences. Third task was aimed towards an important and unique contribution of the project, which was to develop a unified uncertainty quantification framework by integrating MFIX-DEM with a graphical user interface (GUI) driven uncertainty quantification (UQ) engine, i.e., MFIX-GUI and PSUADE. The goal was to enable a user with only modest knowledge of statistics to effectively utilize the UQ framework offered with MFIX-DEM Phi to perform UQ analysis routinely. For Phase 1, a proof-of-concept demonstration of the proposed framework was completed and shared. Direct industry involvement was one of the key virtues of this project, which was performed through forth task. For this purpose, even at the proposal stage, the project team received strong interest in the proposed capabilities from two major corporations, which were further expanded throughout Phase 1 and a new collaboration with another major corporation from chemical industry was also initiated. The level of interest received and continued collaboration for the project during Phase 1 clearly shows the relevance and potential impact of the project for the industrial users.

42 ENGINEERING↗
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
DOE OSTI2020

Model Development for Multi-Component Fuel Vaporization and Flash Boiling

The objectives of this project are to improve the multi-component fuel droplet and film vaporization models used in internal combustion engine simulation, and to develop a comprehensive model to predict the characteristics of multi-component flash boiling spray. This work explores two approaches to fuel composition treatment for modeling multi-component fuel vaporization: one based on discretization and surrogates, and the other based on continuous thermodynamic distribution of fuel properties. The experimental data collected for model validation are done under three fuel form factors: droplet, spray, and thin film. The study on sprays also include experimentation under non-flash and flash boiling conditions, a phenomenon that enhances fuel vaporization. The main goals of this work are: Design and develop a multi-component fuel droplet and wall film vaporization model using both discrete and continuous thermodynamics methods. Design and develop an analytical model for multi-component flash boiling. Integrate the multi-component droplet and film model into multi-dimensional engine calculations to predict the fuel vaporization process under engine operation condition. Conduct multi-component droplet and fuel film vaporization experiments in a non-combusting chamber to verify the proposed vaporization models. Characterize flash boiling phenomena of multi-component fuel sprays by optical and laser diagnostic techniques. This report will detail the experimental setup and the numerical basis for developing a model to achieve the main goals listed above. Key features of observed multi-component fuel vaporization will be summarized at the end of each experimental sections, and corresponding model performance evaluation will be presented at the end of each model development section.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

NEAMS Technical Area Support in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work transcribed in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include a factor of 10^4 improvement in dependency resolution speed, sorting of user objects, ability to compute residuals and Jacobians together, transfer fixes, support for the mortar method in finite volume discretizations, addition of generalized advection schemes for fluid simulations, 10^2 speedup in some Griffin simulations due to a new matrix-only solve type, and much more.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

High-fidelity multiscale model development for accelerated fuel qualification for high discharge burnup [Poster]

Reactor vendors seek economic benefits associated with increasing nuclear fuel service lifetime in the existing light-water reactors fleet. FFRD phenomena represent a major safety concern which still needs to be addressed. Formation of high burnup structure (HBS) in conjunction with LOCA can lead to relocating fuel to escape the fuel pin and get dispersed into the primary coolant system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Onset of Fluidization in MP-PIC Simulations using MFIX-Exa

Fluidized bed reactors are used across a variety of industries, including for energy processes like pyrolysis that result in low-cost energy products. Design and scale-up of fluidized beds is de-risked by modeling and simulation, utilizing tools like NETL’s MFIX-Exa High-Performance Computing (HPC) code for reacting multiphase flow. This report summarizes an investigation into the breadth of problems to which MFIX-Exa may be applied, specifically with regard to low fluid velocities and the onset of fluidization. A simple fluidization study is conducted both experimentally and numerically for particles of interest, then reactor simulations are compared to cold flow experiments for uniform distributor plates. Approaches for modeling bubble caps are also presented.

discrete particle method↗
DOE OSTI2012

Scalar-Filtered Mass-Density-Function Simulation of Swirling Reacting Flows on Unstructured Grids

The scalar-filtered mass-density-function methodology is employed for large-eddy simulation of two swirl-stabilized nonpremixed flames. These are low-swirl (SM1) and high-swirl (SMA2) Sydney methane flames, both of which have been the subject of detailed laboratory measurements. Combustion chemistry is modeled via a flamelet model for the low-swirl flame, and a detailed finite-rate kinetics model for the high-swirl flame. The scalar-filtered mass-density-function is simulated by a Lagrangian Monte Carlo method on a domain discretized by unstructured grids. The simulated results are assessed via comparison with experimental data and show very good agreement. This demonstrates the capability of scalar-filtered mass-density-function for large-eddy simulation of complex flows and warrants future applications of the methodology for large-eddy simulation of practical combustor configurations.

42 ENGINEERING↗
DOE OSTI2013

Mathematical and Numerical Aspects of the Adaptive Fast Multipole Poisson-Boltzmann Solver

This paper summarizes the mathematical and numerical theories and computational elements of the adaptive fast multipole Poisson-Boltzmann (AFMPB) solver. We introduce and discuss the following components in order: the Poisson-Boltzmann model, boundary integral equation reformulation, surface mesh generation, the nodepatch discretization approach, Krylov iterative methods, the new version of fast multipole methods (FMMs), and a dynamic prioritization technique for scheduling parallel operations. For each component, we also remark on feasible approaches for further improvements in efficiency, accuracy and applicability of the AFMPB solver to large-scale long-time molecular dynamics simulations. Lastly, the potential of the solver is demonstrated with preliminary numerical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

When ancient numerical demons meet physics-informed machine learning: adjoint-based gradients for implicit differentiable modeling

Recent advances in differentiable modeling, a genre of physics-informed machine learning that trains neural networks (NNs) together with process-based equations, have shown promise in enhancing hydrological models' accuracy, interpretability, and knowledge-discovery potential. Current differentiable models are efficient for NN-based parameter regionalization, but the simple explicit numerical schemes paired with sequential calculations (operator splitting) can incur numerical errors whose impacts on models' representation power and learned parameters are not clear. Implicit schemes, however, cannot rely on automatic differentiation to calculate gradients due to potential issues of gradient vanishing and memory demand. Here we propose a “discretize-then-optimize” adjoint method to enable differentiable implicit numerical schemes for the first time for large-scale hydrological modeling. The adjoint model demonstrates comprehensively improved performance, with Kling–Gupta efficiency coefficients, peak-flow and low-flow metrics, and evapotranspiration that moderately surpass the already-competitive explicit model. Therefore, the previous sequential-calculation approach had a detrimental impact on the model's ability to represent hydrological dynamics. Furthermore, with a structural update that describes capillary rise, the adjoint model can better describe baseflow in arid regions and also produce low flows that outperform even pure machine learning methods such as long short-term memory networks. The adjoint model rectified some parameter distortions but did not alter spatial parameter distributions, demonstrating the robustness of regionalized parameterization. Despite higher computational expenses and modest improvements, the adjoint model's success removes the barrier for complex implicit schemes to enrich differentiable modeling in hydrology.

58 GEOSCIENCES↗
NASA NTRS1974

Optimal aeroelastic design of an oblique wing structure

A procedure is presented for determining the optimal cover panel thickness of a wing structure to meet specified strength and static aeroelastic divergence requirements for minimum weight. Efficient reanalysis techniques using discrete structural and aerodynamic methods are used in conjunction with redesign algorithms driven by optimality criteria. The optimality conditions for the divergence constraint are established, and expressions are obtained for derivatives of the dynamic pressure at divergence with respect to design variables. The procedure is applied to an oblique wing aircraft where strength and stiffness are critical design considerations for sizing the cover thickness of the wing structure.

Gwin, L. B.↗