Search NASA⌕ Search

SEARCH · Search NASA

Results for “Basis function”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fourier Analysis Method for Wireless Power Transfer Coil Design

Wireless power transfer (WPT) systems for Electric Vehicles (EVs) are designed to meet specifications such as stray field, power transfer, efficiency, and ground clearance. Typical design approaches include iterative analysis of predetermined coil geometries to identify candidates that meet these constraints. This work instead directly generates WPT coil shapes and magnetic fields to meet specifications and constraints through the optimization of Fourier basis function coefficients. The proposed Fourier Analysis Method (FAM) applies to arbitrary planar coil geometries and does not rely on iterative finite-element analysis (FEA) simulations. This flexibility allows for rapid design evaluation across a larger range of coil geometries and design specifications. A prototype coil is built to compare FAM outputs to experimental measurements and FEA simulations. The FAM is then used to illustrate the tradeoff of coil current and stray field for a given power level showing that the method is capable of generating optimized coil shapes to meet arbitrary field constraints.

Foote, Andrew↗
DOE OSTI2012

Pseudo spectral Chebyshev representation of few-group cross sections on sparse grids

This paper presents a pseudo spectral method for representing few-group homogenised cross sections, based on hierarchical polynomial interpolation. The interpolation is performed on a multi-dimensional sparse grid built from Chebyshev nodes. The representation is assembled directly from the samples using basis functions that are constructed as tensor products of the classical one-dimensional Lagrangian interpolation functions. The advantage of this representation is that it combines the accuracy of Chebyshev interpolation with the efficiency of sparse grid methods. As an initial test, this interpolation method was used to construct a representation for the two-group macroscopic cross sections of a VVER pin cell. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

A non-linear discontinuous Petrov-Galerkin method for removing oscillations in the solution of the time-dependent transport equation

This paper describes a Non-Linear Discontinuous Petrov-Galerkin method and its application to the one-speed Boltzmann Transport Equation (BTE) for space-time problems. The purpose of the method is to remove unwanted oscillations in the transport solution which occur in the vicinity of sharp flux gradients, while improving computational efficiency and numerical accuracy. This is achieved by applying artificial dissipation in the solution gradient direction, internal to an element using a novel finite element (FE) Riemann approach. The added dissipation is calculated at each node of the finite element mesh based on local behaviour of the transport solution on both the spatial and temporal axes of the problem. Thus a different dissipation is used in different elements. The magnitude of dissipation that is used is obtained from a gradient-informed scaling of the advection velocities in the stabilisation term. This makes the method in its most general form non-linear. The method is implemented within a very general finite element Riemann framework. This makes it completely independent of choice of angular basis function allowing one to use different descriptions of the angular variation. Results show the non-linear scheme performs consistently well in demanding time-dependent multi-dimensional neutron transport problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Application of adaptive hierarchical sparse grid collocation to the uncertainty quantification of nuclear reactor simulators

Recent efforts in the application of uncertainty quantification to nuclear systems have utilized methods based on generalized perturbation theory and stochastic sampling. While these methods have proven to be effective they both have major drawbacks that may impede further progress. A relatively new approach based on spectral elements for uncertainty quantification is applied in this paper to several problems in reactor simulation. Spectral methods based on collocation attempt to couple the approximation free nature of stochastic sampling methods with the determinism of generalized perturbation theory. The specific spectral method used in this paper employs both the Smolyak algorithm and adaptivity by using Newton-Cotes collocation points along with linear hat basis functions. Using this approach, a surrogate model for the outputs of a computer code is constructed hierarchically by adaptively refining the collocation grid until the interpolant is converged to a user-defined threshold. The method inherently fits into the framework of parallel computing and allows for the extraction of meaningful statistics and data that are not within reach of stochastic sampling and generalized perturbation theory. This paper aims to demonstrate the advantages of spectral methods-especially when compared to current methods used in reactor physics for uncertainty quantification-and to illustrate their full potential. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Effective Elastic Coefficients for Composites with Reinforcement Under Shear External Loading - 18096

The elastic coefficients of composite material with embedded reinforcement bar that has both longitudinal and transverse ribs, under shear loading, are calculated by solving an elasto-static boundary-value problem (bvp) defined in a unit cell on the microscale. The solutions to the bvp's are obtained numerically by using finite elements of quadratic 3-D basis functions in ABAQUS for all of the three macro-scale shear strains imposed externally. The effective elastic coefficients are shown to be a bit larger for external shear deformations on the macro-scale when it encounters longitudinal rib during the shear. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2020

Determination of Elastic Coefficients for a Composite with Crossing Steel - 20086

For a composite material with an embedding of crossing reinforcement bars, the elastic coefficients are calculated by solving a boundary-value problem (elastostatic problem). The background material is concrete whereas the bars are made of steel. The bars cross each other horizontally at a right angle. The boundary-value problem is solved numerically by using 3-D quadratic basis functions in ABAQUS{sup R}. Under external shear loading conditions, the effective elastic coefficients on the macroscale are quite larger for a cross-section in which the bars are immediately below the ones in the cross-section normal to it. The shear modulus in the off-plane direction is about one and half times larger than that in the in-plane direction. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Improved Variable-Transformation Algorithm for the Hexagonal Variational Nodal Method

Hexagonal Variational Nodal Method (VNM) has been widely employed in fast reactor core calculation codes. It provides the weak form solution of the second order transport equation. It employs a direct three-dimensional expansion instead of the transverse integration technique, with the singular term caused by transverse integration over the hexagonal nodal avoided. However, as the increasing Pn expansion order in angle, the number of moments increases dramatically, which significantly increases the requirements in both memory and computational time. As a consequence, Yang proposed the Variable-Transformation Algorithm. It firstly used the reverse Cuthill-McKee (RCMK) reordering method to transform the response matrices into block-diagonal. Then, for the unreducible matrices, orthogonal transformation was applied, after which the matrices were sparse enough to be reduced by using the RCMK reordering once more. In this summary, based on a newly developed hexagonal VNM code NECP-VioletHex, the Variable Transformation Algorithm was investigated and improved by adopting more constraints to the basis functions to further take advantage of the symmetric properties of the hexagonal-z nodes. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Verification of the 2D Finite Element SP3 Solutions

Method(FEM) is the free geometry handling. However many of the geometries of the nuclear engineering problems are not unstructured but rather structured, such as a hexahedral, or triangular prism. Thus, the need for the use of FEM has less portion in nuclear engineering. However it is also certain that it is very hard to exclude the possibility for use of an unstructured geometry such as a PBMR core problem, fusion blanket analysis of the source problem, another design issue with an unstructured geometry. In addition, research regarding the acceleration scheme such as Coarse Mesh Finite Difference (CMFD) showed up in nuclear engineering in recent years. By considering these facts, the FEM cannot be totally exclusive. To utilize the FEM for nuclear engineering problems, several types of assistance are needed. The first is parallelization because of its slow calculation speed. The second is an acceleration scheme. With these two types of assistance the FEM can be utilized freely in the nuclear engineering field. In this research, before conducting parallelization and acceleration, the final verification for the implementation of the up-scattering treatment with an arbitrary geometry conducted through a comparison with the McCARD solution through the well-known C5G7 2D benchmark problem. Note that the earlier work for verification of the FEM for the multi-group problem without up-scattering and with a structured geometry including a mathematical formulation has already been done. The main points that will be tested here are as follows. The first is a check of the right implementation for the multi-group problem with an up-scattering cross section from the perspectives of the multiplication factor and power. The second is the transport effect for this problem by comparing the solutions from P{sub 1} and Simplified P{sub 3}(SP{sub 3}). The third is the effectiveness of the higher order basis function for the same number of unknowns. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A Domain-Decomposed A-ϕ Formulation Based on Lagrange Multipliers for Low-Frequency Problems

A domain-decomposed A-ϕ formulation based on Lagrange multipliers is proposed to simulate low-frequency elec- tromagnetic problems. This method partitions the computational domain into smaller subdomains, allowing each subdomain to be independently formulated using Lagrange multipliers as Dirichlet boundary conditions, while ensuring continuity of the fields across the interfaces. A mixed finite element method, utilizing both vector and scalar basis functions, is employed to discretize the formulation, resulting in a global system to be solved. The proposed method is validated using TEAM Problem 7 at 50 Hz, demonstrating its effectiveness in handling complex geometries and addressing the low-frequency breakdown issues commonly encountered in traditional finite element methods.

Hossain, Amzad↗
DOE OSTI2015

Accuracy of Lagrange-sinc functions as a basis set for electronic structure calculations of atoms and molecules

We developed a self-consistent field program based on Kohn-Sham density functional theory using Lagrange-sinc functions as a basis set and examined its numerical accuracy for atoms and molecules through comparison with the results of Gaussian basis sets. The result of the Kohn-Sham inversion formula from the Lagrange-sinc basis set manifests that the pseudopotential method is essential for cost-effective calculations. The Lagrange-sinc basis set shows faster convergence of the kinetic and correlation energies of benzene as its size increases than the finite difference method does, though both share the same uniform grid. Using a scaling factor smaller than or equal to 0.226 bohr and pseudopotentials with nonlinear core correction, its accuracy for the atomization energies of the G2-1 set is comparable to all-electron complete basis set limits (mean absolute deviation ≤1 kcal/mol). The same basis set also shows small mean absolute deviations in the ionization energies, electron affinities, and static polarizabilities of atoms in the G2-1 set. In particular, the Lagrange-sinc basis set shows high accuracy with rapid convergence in describing density or orbital changes by an external electric field. Moreover, the Lagrange-sinc basis set can readily improve its accuracy toward a complete basis set limit by simply decreasing the scaling factor regardless of systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1975

Determination of discontinuities in the derivative of a piecewise-smooth function on the basis of the results of measurements

A piecewise-smooth function with discontinuity in the first derivative on a given interval is considered. The values of the function are measured at a sequence of points in the interval and a random error is included in the measurements. A method is proposed to estimate the position of the discontinuity in the derivative. Regression lines are associated with each measurement point and account for k - 1 points preceding or following the point. The estimate for the position of the discontinuity is the measurement point with the largest angle between the regression lines. The error in the estimate is analyzed and the results are verified.

Gurin, L. S.↗
DOE OSTI2021

Achieving performance portability in Gaussian basis set density functional theory on accelerator based architectures in NWChemEx

The numerical integration of the exchange–correlation (XC) potential is one of the primary computational bottlenecks in Gaussian basis set Kohn–Sham density functional theory (KS-DFT). To achieve optimal performance and accuracy, care must be taken in this numerical integration to preserve local sparsity as to allow for near linear weak scaling with system size. This leads to an integration scheme with several performance critical kernels which must be hand optimized for each architecture of interest. As the set of available accelerator hardware goes more diverse, a key challenge for developers of KS-DFT software is to maintain performance portability across a wide range of computational architectures. In this article, we examine a modular software design pattern which decouples the implementation details of performance critical kernels from the expression of high-level algorithmic workflows in a device-agnostic language such as C++; thus allowing for developers to target existing and emerging accelerator hardware within a single code base. We consider the efficacy of such a design pattern in the numerical integration of the XC potential by demonstrating its ability to achieve performance portability across a set of accelerator architectures which are representative of those on current and future U.S. Department of Energy Leadership Computing Facilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Calculation of the Energy of the Interelectron Interaction in Molecules in a Basis of Slater Functions

A calculation of the electron interaction energy for the CH molecule with an open electron shell is performed using a basis of Slater atomic orbitals. The molecular orbitals are represented in the form of a linear combination of the 1s, 2s, 2p{sub x}, 2p{sub y}, and 2p{sub z} atomic orbitals of the C atom and the 1s orbital of the H atom. As a result of solving the Hartree–Fock–Roothaan equation, numerical values of the coefficients in this linear combination have been found.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Adaptive local basis set for Kohn–Sham density functional theory in a discontinuous Galerkin framework II: Force, vibration, and molecular dynamics calculations

Recently, we have proposed the adaptive local basis set for electronic structure calculations based on Kohn–Sham density functional theory in a pseudopotential framework. The adaptive local basis set is efficient and systematically improvable for total energy calculations. In this paper, we present the calculation of atomic forces, which can be used for a range of applications such as geometry optimization and molecular dynamics simulation. We demonstrate that, under mild assumptions, the computation of atomic forces can scale nearly linearly with the number of atoms in the system using the adaptive local basis set. We quantify the accuracy of the Hellmann–Feynman forces for a range of physical systems, benchmarked against converged planewave calculations, and find that the adaptive local basis set is efficient for both force and energy calculations, requiring at most a few tens of basis functions per atom to attain accuracies required in practice. Since the adaptive local basis set has implicit dependence on atomic positions, Pulay forces are in general nonzero. However, we find that the Pulay force is numerically small and systematically decreasing with increasing basis completeness, so that the Hellmann–Feynman force is sufficient for basis sizes of a few tens of basis functions per atom. We verify the accuracy of the computed forces in static calculations of quasi-1D and 3D disordered Si systems, vibration calculation of a quasi-1D Si system, and molecular dynamics calculations of H{sub 2} and liquid Al–Si alloy systems, where we show systematic convergence to benchmark planewave results and results from the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Functional and structural basis of E. coli enolase inhibition by SF2312: a mimic of the carbanion intermediate

Many years ago, the natural secondary metabolite SF2312, produced by the actinomycete Micromonospora, was reported to display broad spectrum antibacterial properties against both Gram-positive and Gram-negative bacteria. Recent studies have revealed that SF2312, a natural phosphonic acid, functions as a potent inhibitor of human enolase. The mechanism of SF2312 inhibition of bacterial enolase and its role in bacterial growth and reproduction, however, have remained elusive. In this work, we detail a structural analysis of E. coli enolase bound to both SF2312 and its oxidized imide-form. Our studies support a model in which SF2312 acts as an analog of a high energy intermediate formed during the catalytic process. Biochemical, biophysical, computational and kinetic characterization of these compounds confirm that altering features characteristic of a putative carbanion (enolate) intermediate significantly reduces the potency of enzyme inhibition. When SF2312 is combined with fosfomycin in the presence of glucose-6 phosphate, significant synergy is observed. This suggests the two agents could be used as a potent combination, targeting distinct cellular mechanism for the treatment of bacterial infections. Together, our studies rationalize the structure-activity relationships for these phosphonates and validate enolase as a promising target for antibiotic discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Universal reduced basis for the calibration of covariant energy density functionals

The reduced basis method is used to construct a “universal” basis of Dirac orbitals that may be applicable throughout the nuclear chart to calibrate covariant energy density functionals. Relative to the successful development of a reduced basis emulator for the nonrelativistic Schrödinger equation, the Dirac equation adds an extra layer of complexity due to the existence of negative energy states, which complicates building an efficient reduced basis. However, once this problem is mitigated, the resulting reduced basis is able to accurately and efficiently reproduce the high-fidelity model at a fraction of the computational cost. We are confident that the resulting reduced basis will serve as a foundational element in developing rapid and accurate emulators. In turn, these emulators will play a critical role in the Bayesian optimization of covariant energy density functionals.

Bayesian methods↗