Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Implementation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗
DOE OSTI2016

Using Hybrid Model OpenSHMEM + CUDA to Implement the SHOC Benchmark Suite

This work describes the process of porting the Scalable HeterOgeneous Computing (SHOC) benchmark suite from the hybrid MPI + CUDA implementation to OpenSHMEM + CUDA. SHOC includes a wide variety of benchmark kernels used to measure accelerator performance in both single node and cluster configurations. The hybrid model implementation attempts to place all major computation on accelerator devices, and uses MPI to synchronize and aggregate results. In some cases, MPI Groups are used to gradually reduce the number of accelerators used for computation as the problem size drops. Porting this behavior to OpenSHMEM required implementing several synchronizing collective operations, and using SHMEM teams to replace MPI Group functionality. Benchmark results on a Cray XK7 system with one GPU per compute node show that SHMEM performance is equal to MPI performance in these hybrid tasks. These results and porting experience show that using OpenSHMEM for accelerator devices benefits from adding functionality for synchronization and teams, and would further benefit from adding support for communication within accelerator kernels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Implementing the three-particle quantization condition for π + π + K + and related systems

Recently, the formalism needed to relate the finite-volume spectrum of systems of nondegenerate spinless particles has been derived. In this work we discuss a range of issues that arise when implementing this formalism in practice, provide further theoretical results that can be used to check the implementation, and make available codes for implementing the three-particle quantization condition. Specifically, we discuss the need to modify the upper limit of the cutoff function due to the fact that the left-hand cut in the scattering amplitudes for two nondegenerate particles moves closer to threshold; we describe the decomposition of the three-particle amplitude K df,3 into the matrix basis used in the quantization condition, including both s and p waves, with the latter arising in the amplitude for two nondegenerate particles; we derive the threshold expansion for the lightest three-particle state in the rest frame up to O(1/L 5 ); and we calculate the leading-order predictions in chiral perturbation theory for K df,3 in the π + π + K + and π + K + K + systems. We focus mainly on systems with two identical particles plus a third that is different (“2+1” systems). We describe the formalism in full detail, and present numerical explorations in toy models, in particular checking that the results agree with the threshold expansion, and making a prediction for the spectrum of π + π + K + levels using the two- and three-particle interactions predicted by chiral perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Implementing the three-particle quantization condition including higher partial waves

We present an implementation of the relativistic three-particle quantization condition including both s - and d -wave two-particle channels. For this, we develop a systematic expansion of the three-particle K matrix, \( \mathcal{K} \) df,3 , about threshold, which is the generalization of the effective range expansion of the two-particle K matrix, \( \mathcal{K} \) 2 . Relativistic invariance plays an important role in this expansion. We find that d -wave two-particle channels enter first at quadratic order. We explain how to implement the resulting multichannel quantization condition, and present several examples of its application. We derive the leading dependence of the threshold three-particle state on the two-particle d -wave scattering amplitude, and use this to test our implementation. We show how strong two-particle d -wave interactions can lead to significant effects on the finite-volume three-particle spectrum, including the possibility of a generalized three-particle Efimov-like bound state. We also explore the application to the 3 π + system, which is accessible to lattice QCD simulations, where we study the sensitivity of the spectrum to the components of \( \mathcal{K} \) df,3 . Finally, we investigate the circumstances under which the quantization condition has unphysical solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Numerical implementation, verification and validation of two-phase flow four-equation drift flux model with Jacobian-free Newton–Krylov method

This study presents a numerical investigation on using the Jacobian-free Newton–Krylov (JFNK) method to solve the two-phase flow four-equation drift flux model with realistic constitutive correlations (‘closure models’). The drift flux model is based on Isshi and his collaborators’ work. Additional constitutive correlations for vertical channel flow, such as two-phase flow pressure drop, flow regime map, wall boiling and interfacial heat transfer models, were taken from the RELAP5-3D Code Manual and included to complete the model. The staggered grid finite volume method and fully implicit backward Euler method was used for the spatial discretization and time integration schemes, respectively. The Jacobian-free Newton–Krylov method shows no difficulty in solving the two-phase flow drift flux model with a discrete flow regime map. In addition to the Jacobian-free approach, the preconditioning matrix is obtained by using the default finite differencing method provided in the PETSc package, and consequently the labor-intensive implementation of complex analytical Jacobian matrix is avoided. Extensive and successful numerical verification and validation have been performed to prove the correct implementation of the models and methods. Code-to-code comparison with RELAP5-3D has further demonstrated the successful implementation of the drift flux model.

42 ENGINEERING↗
DOE OSTI2020

Efficient discontinuous Galerkin implementations and preconditioners for implicit unsteady compressible flow simulations

This work presents and compares efficient implementations of high-order discontinuous Galerkin methods: a modal matrix-free discontinuous Galerkin (DG) method, a hybridizable discontinuous Galerkin (HDG) method, and a primal formulation of HDG, applied to the implicit solution of unsteady compressible flows. The matrix-free implementation allows for a reduction of the memory footprint of the solver when dealing with implicit time-accurate discretizations. HDG reduces the number of globally-coupled degrees of freedom relative to DG, at high order, by statically condensing element-interior degrees of freedom from the system in favor of face unknowns. The primal formulation further reduces the element-interior degrees of freedom by eliminating the gradient as a separate unknown. This paper introduces a p-multigrid preconditioner implementation for these discretizations and presents results for various flow problems. Benefits of the p-multigrid strategy relative to simpler, less expensive, preconditioners are observed for stiff systems, such as those arising from low-Mach number flows at high-order approximation. The p-multigrid preconditioner also shows excellent scalability for parallel computations. Finally, additional savings in both speed and memory occur with a matrix-free/reduced version of the preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Implementation of replica-exchange umbrella sampling in GAMESS

In GAMESS-US quantum chemistry package, both replica-exchange molecular dynamics (REMD) and the umbrella sampling (US) methods are available. However, prior to this work, the replica-exchange umbrella sampling (REUS) method, combining REMD and US, was not implemented. Here, REUS affords a more efficient conformational sampling than the conventional US method and is often used in MD simulations. We implemented REUS in GAMESS-US and improved some routines of the existing REMD code. A density-functional tight-binding REMD (DFTB-REMD) simulation of alanine deca-peptide and a DFTB-MD REUS simulation of proton transfer in malonaldehyde were performed. The results of DFTB-MD REUS simulations performed with three different simulation packages were compared to validate our GAMESS-US implementation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

we report a multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

FutureGen 2.0 Monitoring Program: An Overview of the Monitoring Approach and Technologies Selected for Implementation

The FutureGen 2.0 Project will design and build a first-of-its-kind, near-zero emissions coal-fueled power plant with carbon capture and storage (CCS). To assess storage site performance and meet the regulatory requirements of the Class VI Underground Injection Control (UIC) Program for CO2 Geologic Sequestration, the FutureGen 2.0 project will implement a suite of monitoring technologies designed to 1) evaluate CO2 mass balance and 2) detect any unforeseen loss in CO2 containment. The monitoring program will include direct monitoring of the injection stream and reservoir, and early-leak-detection monitoring directly above the primary confining zone. It will also implement an adaptive monitoring strategy whereby monitoring results are continually evaluated and the monitoring network is modified as required, including the option to drill additional wells in out-years. Wells will be monitored for changes in CO2 concentration and formation pressure, and other geochemical/isotopic signatures that provide indication of CO2 or brine leakage. Indirect geophysical monitoring technologies that were selected for implementation include passive seismic, integrated surface deformation, time-lapse gravity, and pulsed neutron capture logging. Near-surface monitoring approaches that have been initiated include surficial aquifer and surface- water monitoring, soil-gas monitoring, atmospheric monitoring, and hyperspectral data acquisition for assessment of vegetation conditions. Initially, only the collection of baseline data sets is planned; the need for additional near- surface monitoring will be continually evaluated throughout the design and operational phases of the project, and selected approaches may be reinstituted if conditions warrant. Given the current conceptual understanding of the subsurface environment, early and appreciable impacts to near-surface environments are not expected.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

LuGo: An enhanced quantum phase estimation implementation

Quantum Phase Estimation (QPE) is a cardinal algorithm in quantum computing that plays a crucial role in various applications, including cryptography, molecular simulation, and solving systems of linear equations. However, the standard implementation of QPE faces challenges related to time complexity and circuit depth, which limit its practicality for large-scale computations. We introduce LuGo, a novel framework designed to enhance the performance of QPE by reducing circuit duplication, as well as using parallelization techniques to achieve faster generation of the QPE circuit and gate reduction. We validate the effectiveness of our framework by generating quantum linear solver circuits, which require both QPE and inverse QPE, to solve linear systems of equations. LuGo achieves significant improvements in both computational efficiency and hardware requirements without compromising on accuracy. Compared to a standard QPE implementation, LuGo reduces time consumption to generate a circuit that solves a 2 6 × 2 6 system matrix by a factor of 50.68 and over 31× reduction of quantum gates and circuit depth, with no fidelity loss on an ideal quantum simulator. Furthermore, we demonstrated the versatility and scalability of LuGo enabled HHL algorithm by simulating a canonical Hele-Shaw fluid problem using a quantum simulator. With these advantages, LuGo paves the way for more efficient implementations of QPE, enabling broader applications across several quantum computing domains.

Quantum algorithm↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
DOE OSTI2014

Clinical Implementation of an Online Adaptive Plan-of-the-Day Protocol for Nonrigid Motion Management in Locally Advanced Cervical Cancer IMRT

Purpose: To evaluate the clinical implementation of an online adaptive plan-of-the-day protocol for nonrigid target motion management in locally advanced cervical cancer intensity modulated radiation therapy (IMRT). Methods and Materials: Each of the 64 patients had four markers implanted in the vaginal fornix to verify the position of the cervix during treatment. Full and empty bladder computed tomography (CT) scans were acquired prior to treatment to build a bladder volume-dependent cervix-uterus motion model for establishment of the plan library. In the first phase of clinical implementation, the library consisted of one IMRT plan based on a single model-predicted internal target volume (mpITV), covering the target for the whole pretreatment observed bladder volume range, and a 3D conformal radiation therapy (3DCRT) motion-robust backup plan based on the same mpITV. The planning target volume (PTV) combined the ITV and nodal clinical target volume (CTV), expanded with a 1-cm margin. In the second phase, for patients showing >2.5-cm bladder-induced cervix-uterus motion during planning, two IMRT plans were constructed, based on mpITVs for empty-to-half-full and half-full-to-full bladder. In both phases, a daily cone beam CT (CBCT) scan was acquired to first position the patient based on bony anatomy and nodal targets and then select the appropriate plan. Daily post-treatment CBCT was used to verify plan selection. Results: Twenty-four and 40 patients were included in the first and second phase, respectively. In the second phase, 11 patients had two IMRT plans. Overall, an IMRT plan was used in 82.4% of fractions. The main reasons for selecting the motion-robust backup plan were uterus outside the PTV (27.5%) and markers outside their margin (21.3%). In patients with two IMRT plans, the half-full-to-full bladder plan was selected on average in 45% of the first 12 fractions, which was reduced to 35% in the last treatment fractions. Conclusions: The implemented online adaptive plan-of-the-day protocol for locally advanced cervical cancer enables (almost) daily tissue-sparing IMRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Implementation of dietary methionine restriction using casein after selective, oxidative deletion of methionine

Dietary methionine restriction (MR) is normally implemented using diets formulated from elemental amino acids (AA) that reduce methionine content to 0.17%. However, translational implementation of MR with elemental AA-based diets is intractable due to poor palatability. To solve this problem and restrict methionine using intact proteins, casein was subjected to mild oxidation to selectively reduce methionine. Diets were then formulated using oxidized casein, adding back methionine to produce a final concentration of 0.17%. The biological efficacy of dietary MR using the oxidized casein (Ox Cas) diet was compared with the standard elemental MR diet in terms of the behavioral, metabolic, endocrine, and transcriptional responses to the four diets. The Ox Cas MR diet faithfully reproduced the expected physiological, biochemical, and transcriptional responses in liver and inguinal white adipose tissue. Collectively, these findings demonstrate that dietary MR can be effectively implemented using casein after selective oxidative reduction of methionine.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Efficient implementation of the many-body Reactive Bond Order (REBO) potential on GPU

The second generation Reactive Bond Order (REBO) empirical potential is commonly used to accurately model a wide range hydrocarbon materials. It is also extensible to other atom types and interactions. REBO potential assumes complex multi-body interaction model, that is difficult to represent efficiently in the SIMD or SIMT programming model. Hence, despite its importance, no efficient GPGPU implementation has been developed for this potential. Here we present a detailed description of a highly efficient GPGPU implementation of molecular dynamics algorithm using REBO potential. The presented algorithm takes advantage of rarely used properties of the SIMT architecture of a modern GPU to solve difficult synchronizations issues that arise in computations of multi-body potential. Techniques developed for this problem may be also used to achieve efficient solutions of different problems. The performance of proposed algorithm is assessed using a range of model systems. It is compared to highly optimized CPU implementation (both single core and OpenMP) available in LAMMPS package. These experiments show up to 6x improvement in forces computation time using single processor of the NVIDIA Tesla K80 compared to high end 16-core Intel Xeon processor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Multiresolution molecular mechanics: Implementation and efficiency

Atomistic/continuum coupling methods combine accurate atomistic methods and efficient continuum methods to simulate the behavior of highly ordered crystalline systems. Coupled methods utilize the advantages of both approaches to simulate systems at a lower computational cost, while retaining the accuracy associated with atomistic methods. Many concurrent atomistic/continuum coupling methods have been proposed in the past; however, their true computational efficiency has not been demonstrated. The present work presents an efficient implementation of a concurrent coupling method called the Multiresolution Molecular Mechanics (MMM) for serial, parallel, and adaptive analysis. First, we present the features of the software implemented along with the associated technologies. The scalability of the software implementation is demonstrated, and the competing effects of multiscale modeling and parallelization are discussed. Then, the algorithms contributing to the efficiency of the software are presented. These include algorithms for eliminating latent ghost atoms from calculations and measurement-based dynamic balancing of parallel workload. The efficiency improvements made by these algorithms are demonstrated by benchmark tests. The efficiency of the software is found to be on par with LAMMPS, a state-of-the-art Molecular Dynamics (MD) simulation code, when performing full atomistic simulations. Speed-up of the MMM method is shown to be directly proportional to the reduction of the number of the atoms visited in force computation. Finally, an adaptive MMM analysis on a nanoindentation problem, containing over a million atoms, is performed, yielding an improvement of 6.3–8.5 times in efficiency, over the full atomistic MD method. For the first time, the efficiency of a concurrent atomistic/continuum coupling method is comprehensively investigated and demonstrated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Efficient algorithms and implementations of entropy-based moment closures for rarefied gases

We present efficient algorithms and implementations of the 35-moment system equipped with the maximum-entropy closure in the context of rarefied gases. While closures based on the principle of entropy maximization have been shown to yield very promising results for moderately rarefied gas flows, the computational cost of these closures is in general much higher than for closure theories with explicit closed-form expressions of the closing fluxes, such as Grad's classical closure. Following a similar approach as Garrett et al. (2015) , we investigate efficient implementations of the computationally expensive numerical quadrature method used for the moment evaluations of the maximum-entropy distribution by exploiting its inherent fine-grained parallelism with the parallelism offered by multi-core processors and graphics cards. We show that using a single graphics card as an accelerator allows speed-ups of two orders of magnitude when compared to a serial CPU implementation. To accelerate the time-to-solution for steady-state problems, we propose a new semi-implicit time discretization scheme. The resulting nonlinear system of equations is solved with a Newton type method in the Lagrange multipliers of the dual optimization problem in order to reduce the computational cost. Additionally, fully explicit time-stepping schemes of first and second order accuracy are presented. We investigate the accuracy and efficiency of the numerical schemes for several numerical test cases, including a steady-state shock-structure problem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Implementation of metal-friendly EAM/FS-type semi-empirical potentials in HOOMD-blue: A GPU-accelerated molecular dynamics software

We present an implementation of EAM and FS interatomic potentials, which are widely used in simulating metallic systems, in HOOMD-blue, a software designed to perform classical molecular dynamics simulations using GPU accelerations. We first discuss the details of our implementation and then report extensive benchmark tests. We demonstrate that single-precision floating point operations efficiently implemented on GPUs can produce sufficient accuracy when compared against double-precision codes, as demonstrated in test simulations of calculations of the glass-transition temperature of Cu 64.5 Zr 35.5 , and pair correlation function of liquid Ni 3 Al. Our code scales well with the size of the simulating system on NVIDIA Tesla M40 and P100 GPUs. Compared with another popular software LAMMPS running on 32 cores of AMD Opteron 6220 processors, the GPU/CPU performance ratio can reach as high as 4.6. In conclusion, the source code can be accessed through the HOOMD-blue web page for free by any interested user.

36 MATERIALS SCIENCE↗