Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Comparing Emulation Methods for Computer Models with High Dimensional Output

This Master’s Capstone project will present a comparison of statistical models for computer simulation studies. The four models included in this comparison study were chosen for both their proven and diverse methodologies. We will present a case study on hurricane flood data in the Delaware bay which highlights the strengths and weaknesses of each model when applied to a very large spatial field. As computers have gotten faster, we have become interested in modeling increasingly large spatial fields in both size and resolution. Statistical algorithms that are able to efficiently handle these fields have never been more important. We therefore find this comparison to be extremely topical.

42 ENGINEERING↗
DOE OSTI2022

System and method for computed axial lithography (CAL) for 3D additive manufacturing

In one aspect the present disclosure relates to a system for forming a three dimensional (3D) object from a volume of photo-responsive material contained within a material container. A plurality of optical projection subsystems are adapted to be arranged circumferentially around a material container containing the volume of photo-responsive material, and are controlled by a controller. The optical projection subsystems direct optical projections at a plurality of angles θ through the volume of photo-responsive material. The optical projections are further directed about a z axis extending through the volume of photo-responsive material. The plurality of optical subsystems are further controlled by the controller to provide each of the optical projections with a calculated three-dimensional intensity distribution acting over a fixed temporal exposure period, which is sufficient to at least one of cure or remove selected portions of the volume of photo-responsive material, and leave other portions unmodified, to form a desired 3D part.

Kelly, Brett↗
DOE OSTI2022

System, method, and computer program for creating an internal conforming structure

A system for creating an internal formation of a tubular structure having an inner surface via additive manufacturing. The system broadly includes a computer modeling system and an additive manufacturing system. The computer modeling system may include a processor for generating a lattice cellular component via computer-aided design software according to inputs received from a user. The processor may also generate an internal formation lattice structure based on the lattice cellular component and modify the lattice structure to follow and/or conform to the curvature of the inner surface of the outer wall of the tubular structure. The additive manufacturing system may be configured to produce the lattice structure and the tubular structure via additive manufacturing material deposited layer by layer according to the lattice structure.

42 ENGINEERING↗
DOE OSTI2024

System and method for computed axial lithography (CAL) for 3D additive manufacturing

In one aspect the present disclosure relates to a system for forming a three dimensional (3D) object from a volume of photo-responsive material contained within a material container. A plurality of optical projection subsystems are adapted to be arranged at least partially circumferentially around a material container containing the volume of photo-responsive material, and are controlled by a controller. The optical projection subsystems direct optical projections at a plurality of angles θ through the volume of photo-responsive material. The optical projections are further directed about a z axis extending through the volume of photo-responsive material. The plurality of optical subsystems are further controlled by the controller to provide each of the optical projections with a calculated three-dimensional intensity distribution acting over a fixed temporal exposure period, which is sufficient to at least one of cure or remove selected portions of the volume of photo-responsive material, and leave other portions unmodified, to form a desired 3D part.

Kelly, Brett↗
DOE OSTI2021

Gas-phase ion-molecule interactions in a collision reaction cell with triple quadrupole-inductively coupled plasma mass spectrometry: Investigations with N 2 O as the reaction gas

Nitrous oxide (N 2 O) was used as a reaction gas to investigate the gas phase ion-molecule interactions using the Agilent 8900 QQQ-ICP-MS. A multi-element standard containing 45 elements with masses ranging from 9 to 208 u was measured in the presence and absence of N 2 O. The main product ion species observed were oxides and nitrides. Comparison of the N 2 O reaction results with similar measurements conducted with O 2 revealed that N 2 O was more effective at forming oxides in general: the elements Cd and Pb were shown to produce oxides with N 2 O where the reaction did not occur with O 2 . Nitrous oxide was also shown to produce a significant amount of nitride species in a few cases. The general reactivity was shown to be consistent with density functional theory (DFT)-predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels at least 1% of the unreacted ion. Our results show that reaction enthalpy is a reasonable predictor of reactivity with N2O on the timescales of the interactions in non-thermal ICP-MS/MS systems. Our work demonstrates the utility of two relatively new platforms (commercial elemental ICP-MS/MS and EMSL Arrows interface to the NWChem program suite), which allows for the study of a large number of elements within a short period. While DFT with the basis sets utilized here is not the most accurate computational method, it is also not computationally expensive and is shown to be suitable for predicting gas phase reactivity in the QQQ-ICP-MS for the majority of ions studied. Here, the ease and rapidity of data collection and DFT calculations has the potential to be very impactful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

pnnl/neuromancer

Dynamics-based deep learning methods to modernize current scientific computing methods. Neuromancer is currently capable of solving inverse problems for a system of ordinary differential equations. The functionality includes system identification and constrained optimal control of unknown or partially known ODEs.

Skomski, Elliott↗
DOE OSTI2024

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. Here, the ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Method for enhancing a computer to estimate an uncertainty of an onset of a signal of interest in time-series noisy data

A computer-implemented method of enhancing a computer to estimate an uncertainty of an onset of a signal of interest in time-series noisy data. A first mathematical model of first time series data that contains only noise is calculated. A second mathematical model of second time series data that contains the noise and an onset of a signal of interest in the second time series data is calculated. A difference is evaluated between a first combination, being the first mathematical model and the second mathematical model, and a second combination, being the first time series data and the second time series data, wherein evaluating is performed using a generalized entropy metric. A specific time when an onset of the signal of interest occurs is estimated from the difference. An “a posteriori” distribution is derived for an uncertainty of the specific time at which the onset occurs.

Stracuzzi, David John↗
DOE OSTI2025

Reduced-basis method for few-body bound-state emulation

Recent advances in both theoretical and computational methods have enabled large-scale, precision calculations of the properties of atomic nuclei. With the growing complexity of modern nuclear theory, however, also comes the need for novel methods to perform systematic studies and quantify the uncertainties of models when confronted with experimental data. Here, this study presents an application of such an approach, the reduced basis method, to substantially lower computational costs by constructing a significantly smaller Hamiltonian subspace informed by previous solutions. Our method shows comparable efficiency and accuracy to other dimensionality reduction techniques on an artificial three-body bound system while providing a richer representation of physical information in its projection and training subspace. This methodological advancement can be applied in other contexts and has the potential to greatly improve our ability to systematically explore theoretical models and thus enhance our understanding of the fundamental properties of nuclear systems.

cluster models↗
DOE OSTI2021

From classical thermodynamics to phase-field method

Phase-field method is a density-based computational method at the mesoscale for modeling and predicting the temporal microstructure and property evolution during materials processes. The focus of this article is on connecting the most common phase-field equations to the very basic first and second laws of classical thermodynamics through rudimentary irreversible thermodynamics. It briefly discusses the relations of the continuum phase-field equations to their counter parts at the microscopic and atomic levels. It attempts to clarify the contributions of long-range elastic, electrostatic, and magnetic interactions to domain structure evolution during structural, ferroelectric, and ferromagnetic phase transformations by separating order parameter changes due to the presence of quasi-static fields and those arising from phase transformations. A few examples are presented to demonstrate the possibility of employing the phase-field method to provide guidance to designing materials for optimum properties or discovering novel mesoscale phenomena or new materials functionalities. Here, the article ends with a brief perspective on a number of potential future directions on the development and applications of phase-field method beyond its traditional applications to structural alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Massively parallel transport sweeps on meshes with cyclic dependencies

When solving the first-order form of the linear Boltzmann equation, a common misconception is that the matrix-free computational method of “sweeping the mesh”, used in conjunction with the Discrete Ordinates method, is too complex or does not scale well enough to be implemented in modern high performance computing codes. This has led to considerable efforts in the development of matrix-based methods that are computationally expensive and is partly driven by the requirements placed on modern spatial discretizations. In particular, modern transport codes are required to support higher order elements, a concept that invariably adds a lot of complexity to sweeps because of the introduction of cyclic dependencies with curved mesh cells. In this article we will present a comprehensive implementation of sweeping, to a piecewise-linear DFEM spatial discretization with particular focus on handling cyclic dependencies and possible extensions to higher order spatial discretizations. We find that these methods are implemented in a new C++ simulation framework called Chi-Tech (). We present some typical simulation results with some performance aspects that one can expect during real world simulations, we also present a scaling study to >100k processes where Chi-Tech maintains greater than 80% efficiency solving a total of 87.7 trillion angular flux unknowns for a 116 group simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Methods and devices for preventing computationally explosive calculations in a computer for model parameters distributed on a hierarchy of geometric simplices

A computer-implemented method of preventing computationally explosive calculations. The method includes obtaining, by a processor of the computer, measured data of one of a physical process or a physical object; performing hierarchical numerical modeling of a physical process inclusive of an Earth model containing at least one of (a) infrastructure in the ground and (b) a formation feature in the ground, wherein predicted data is generated; comparing the measured data to the predicted data to calculate an estimated error; analyzing the estimated error via an inversion process to update the at least one of the Earth model and infrastructure model so as to reduce the estimated error and to determine a final composite Earth model of at least one of the infrastructure and the feature; and using the final composite Earth model to characterize at least one of the process and the physical object.

Source record↗
DOE OSTI2021

Validation of simulated wave energy converter responses to focused waves

Focused waves are increasingly used as a design method to predict extreme loads for offshore structures, such as wave energy converters (WECs). However, the ability of computational methods to accurately simulate the prescribed focused waves and resulting wave structure interactions is not well validated. This study uses two computational fluid dynamic methods, WEC-Sim and STAR-CCM+, to evaluate the responses of two taut-moored WEC-like buoys to three different focused waves. WEC-Sim is a computationally efficient, mid-fidelity model, based on linearised potential flow theory. In contrast, STAR-CCM+ is a high-fidelity, three-dimensional, unsteady, Reynolds-averaged Navier–Stokes-based model. The two buoy geometries considered are a hemispherical-bottom cylinder and a truncated cylinder with a cylindrical moon pool. Three prescribed focused waves are generated in both codes, and the simulated displacements and mooring loads are obtained for each of the hydrodynamic bodies. The resulting STAR-CCM+ generated focused waves have approximately the same accuracy, in comparison with the analytic solution, as experimentally generated focused waves. The WEC-Sim simulated responses, when using viscous coefficients tuned to the measured displacements, are, on average, within 11.5% of experimentally measured values, whereas the STAR-CCM+ simulated responses are, on average, within 14.2% of the experimental values.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Determination of excited state molecular structures from time-resolved gas-phase X-ray scattering

We present a comprehensive investigation of a recently introduced method to determine transient structures of molecules in excited electronic states with sub-ångstrom resolution from time-resolved gas-phase scattering signals. The method, which is examined using time-resolved X-ray scattering data measured on the molecule N-methylmorpholine (NMM) at the Linac Coherent Light Source (LCLS), compares the experimentally measured scattering patterns against the simulated patterns corresponding to a large pool of molecular structures to determine the full set of structural parameters. In addition, we examine the influence of vibrational state distributions and find the effect negligible within the current experimental detection limits, despite that the molecules have a comparatively high internal vibrational energy. Here, the excited state structures determined using three structure pools generated using three different computational methods are in good agreement, demonstrating that the procedure is largely independent of the computational chemistry method employed as long as the pool is sufficiently expansive in the vicinity of the sought structure and dense enough to yield good matches to the experimental patterns.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Experimental equation of state of 11 lanthanide nitrides (NdN to LuN) and pressure induced phase transitions in NdN, SmN, EuN, and GdN

Through an extensive data analysis of powder X-ray diffraction data obtained at pressures up to at least 78 GPa, we report the experimental equations of state for all lanthanide nitrides between NdN and LuN, excluding the radioactive Pm. By fitting the obtained unit cell volumes as a function of pressure with the third order Birch–Murnaghan equation of state, we find that the bulk modulus increases with an increasing lanthanide number from K 0 = 146(12) GPa for NdN to 182(7) GPa in EuN. Hereafter, the bulk modulus reaches a plateau for the rest of the series except for TmN which has a lower bulk modulus. We find that the first derivative of the bulk modulus is around 4 for all compounds except TbN, which displays a significantly different compression behavior. In addition, we find a B1 to B10 pressure-induced phase transition in NdN, SmN, EuN, and GdN at increasingly higher pressures. In fact, we observe that the onset pressure of the phase transition increases linearly with Ln atomic number. From TbN and onwards, we do not observe any sign of a B1 to B10 transition indicating that the transition pressure exceeds the studied pressure range. Therefore, we believe that, for the heavier lanthanides, the linear relationship between the onset pressure and the lanthanide number does not hold and even higher pressures are needed to observe the transition. This coherent study of the series of lanthanide nitrides offers a unique opportunity for benchmark studies of computational methods applied to compounds with 4f electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗