Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Green computing”

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 901 records · Page 50

Navier-Stokes simulations of Orbiter aerodynamic characteristics including pitch trim and bodyflap

An analysis of the longitudinal aerodynamics of the shuttle orbiter in the hypersonic flight regime is made through the use of computational fluid dynamics. Particular attention is given to establishing the cause of the 'pitching moment anomaly,' which occurred on the orbiter's first flight, and to computing the aerodynamics of a complete orbiter configuration at flight conditions. Data from ground-based facilities as well as orbiter flight data are used to validate the computed results. Analysis shows that the pitching moment anomaly is a real-gas chemistry effect that was not simulated in ground-based facilities, which used air as a test gas. Computed flight aerodynamics for the orbiter are within 5% of the measured flight values and trim bodyflap deflections are predicted to within 10%.

Weilmuenster, K. James↗

The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation

We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.

42 ENGINEERING↗

Sum-of-Fractions Method

Sum-of-fractions is a method intended to make sure a subcritical margin for aqueous solutions and slurries of fissionable isotopes exists. The method indicates that a system is subcritical if the sum of the ratios of the mass of each isotope (in a mixture) to its individual minimum subcritical mass limit is less than or equal to one. Historically, the basis of the sum-of-fractions has been derived from allowances given in the American National Standards Institute (ANSI)/ American Nuclear Society (ANS)-8.15-1981. However, the allowance was removed in ANSI/ANS-8.15-2014 due to a lack of technical basis. A methodology was developed to assess the validity of using the sum-of-fractions for water- or polyethylene-moderated systems for the following nuclides: 232U, 233U, 234U, 235U, 237Np, 236Pu, 238Pu, 239Pu, 240Pu, 241Pu, 242Pu, 241Am, 242mAm, 243Am, 242Cm, 243Cm, 244Cm, 245Cm, 246Cm, 247Cm, 249Cf, and 251Cf. The methodology uses available benchmark data for mixtures of 233U, 235U, and 239Pu to establish the calculational margin, and a mass limit reduction to establish the margin of subcriticality. Water- or polyethylene-moderated and -reflected mixtures containing the nuclides are evaluated with the code system, SCALE 6.2.4. Including the calculational margin, subcritical mass limits for each nuclide were computed for optimally water- or polyethylene-moderated and fully reflected systems. These masses were used to create nuclide mixtures in which the sum of the mass to subcritical mass limit ratios is one. The various nuclide mixtures were modeled over a range of moderation and demonstrate the keff does not exceed the calculational margin. For additional assurance of subcriticality, a significant mass reduction is applied to each computed minimum critical mass of the nuclides without adequate benchmark data consistent with the method in ANSI/ANS-8.15-2014.

criticality safety, Actinide↗

A systematic benchmark of the ab initio Bethe-Salpeter equation approach for low-lying optical excitations of small organic molecules

The predictive power of the ab initio Bethe-Salpeter equation (BSE) approach, rigorously based on many-body Green’s function theory but incorporating information from density functional theory, has already been demonstrated for the optical gaps and spectra of solid-state systems. Interest in photoactive hybrid organic/inorganic systems has recently increased and so has the use of the BSE for computing neutral excitations of organic molecules. However, no systematic benchmarks of the BSE for neutral electronic excitations of organic molecules exist. Here, we study the performance of the BSE for the 28 small molecules in Thiel’s widely used time-dependent density functional theory benchmark set [Schreiber et al., J. Chem. Phys. 128, 134110 (2008)]. We observe that the BSE produces results that depend critically on the mean-field starting point employed in the perturbative approach. We find that this starting point dependence is mainly introduced through the quasiparticle energies obtained at the intermediate GW step and that with a judicious choice of starting mean-field, singlet excitation energies obtained from BSE are in excellent quantitative agreement with higher-level wavefunction methods. The quality of the triplet excitations is slightly less satisfactory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relative sensitivity of plastic scintillator: A comparative analysis with 60 Co gamma rays, deuterium–deuterium, and deuterium–tritium neutrons

A plastic scintillator has found extensive application in the realm of high-energy physics and national security science. Many applications in those fields often involve the simultaneous production of photons, neutrons, and charged particles, which makes the relative sensitivity information for these different radiation types important. In this study, we have adopted a multi-head detector comprised of a plastic scintillator and high gain phototubes, which provides a large dynamic range and linearity. A comparative study on the relative sensitivities of plastic scintillators was facilitated by adopting three distinct radiation calibration sources (i.e., 60 Co γ rays, DD neutrons, and DT neutrons). Neutrons from a DD source generate a comparable level of scintillation to gamma rays emitted by 60 Co (i.e., 60 Co-γ/DD-n = 0.92 ± 16%). DT neutrons induce ~3.5 times the scintillation observed with DD neutrons (i.e., DT-n/DD-n = 3.5 ± 28%). In addition, the Geant4 simulation granted us valuable insights into the relative sensitivity of the scintillator. This comparative study will provide a useful database for users in diverse applications.

47 OTHER INSTRUMENTATION↗

TU-AB-202-06: Quantitative Evaluation of Deformable Image Registration in MRI-Guided Adaptive Radiation Therapy

Purpose: To assess the performance of the deformable image registration algorithm used for MRI-guided adaptive radiation therapy using image feature analysis. Methods: MR images were collected from five patients treated on the MRIdian (ViewRay, Inc., Oakwood Village, OH), a three head Cobalt-60 therapy machine with an 0.35 T MR system. The images were acquired immediately prior to treatment with a uniform 1.5 mm resolution. Treatment sites were as follows: head/neck, lung, breast, stomach, and bladder. Deformable image registration was performed using the ViewRay software between the first fraction MRI and the final fraction MRI, and the DICE similarity coefficient (DSC) for the skin contours was reported. The SIFT and Harris feature detection and matching algorithms identified point features in each image separately, then found matching features in the other image. The target registration error (TRE) was defined as the vector distance between matched features on the two image sets. Each deformation was evaluated based on comparison of average TRE and DSC. Results: Image feature analysis produced between 2000–9500 points for evaluation on the patient images. The average (± standard deviation) TRE for all patients was 3.3 mm (±3.1 mm), and the passing rate of TRE<3 mm was 60% on the images. The head/neck patient had the best average TRE (1.9 mm±2.3 mm) and the best passing rate (80%). The lung patient had the worst average TRE (4.8 mm±3.3 mm) and the worst passing rate (37.2%). DSC was not significantly correlated with either TRE (p=0.63) or passing rate (p=0.55). Conclusions: Feature matching provides a quantitative assessment of deformable image registration, with a large number of data points for analysis. The TRE of matched features can be used to evaluate the registration of many objects throughout the volume, whereas DSC mainly provides a measure of gross overlap. We have a research agreement with ViewRay Inc.

60 APPLIED LIFE SCIENCES↗

Scalability of OpenFOAM Density-Based Solver with Runge–Kutta Temporal Discretization Scheme

Compressible density-based solvers are widely used in OpenFOAM, and the parallel scalability of these solvers is crucial for large-scale simulations. In this paper, we report our experiences with the scalability of OpenFOAM’s native rhoCentralFoam solver, and by making a small number of modifications to it, we show the degree to which the scalability of the solver can be improved. The main modification made is to replace the first-order accurate Euler scheme in rhoCentralFoam with a third-order accurate, four-stage Runge-Kutta or RK4 scheme for the time integration. The scaling test we used is the transonic flow over the ONERA M6 wing. This is a common validation test for compressible flows solvers in aerospace and other engineering applications. Numerical experiments show that our modified solver, referred to as rhoCentralRK4Foam, for the same spatial discretization, achieves as much as a 123.2% improvement in scalability over the rhoCentralFoam solver. As expected, the better time resolution of the Runge–Kutta scheme makes it more suitable for unsteady problems such as the Taylor–Green vortex decay where the new solver showed a 50% decrease in the overall time-to-solution compared to rhoCentralFoam to get to the final solution with the same numerical accuracy. Finally, the improved scalability can be traced to the improvement of the computation to communication ratio obtained by substituting the RK4 scheme in place of the Euler scheme. All numerical tests were conducted on a Cray XC40 parallel system, Theta, at Argonne National Laboratory.

Li, Sibo↗

Continental Spatio-Temporal Data Analysis with Linear Spectral Mixture Model Using FOSS

This work demonstrates the development and implementation of a Fully Constrained Least Squares (FCLS) unmixing model developed in C++ programming language with OpenCV package and boost C++ libraries in the NASA Earth Exchange (NEX). Visualization of the results is supported by GRASS GIS and statistical analysis is carried in R in a Linux system environment. FCLS was first tested on computer simulated data with Gaussian noise of various signal-to-noise ratio, and Landsat data of an agricultural scenario and an urban environment using a set of global end members of substrate (soils, sediments, rocks, and non-photosynthetic vegetation), vegetation that includes green photosynthetic plants and dark objects which encompasses absorptive substrate materials, clear water, deep shadows, etc. For the agricultural scenario, a spectrally diverse collection of 11 scenes of Level 1 terrain corrected, cloud free Landsat-5 TM data of Fresno, California, USA were unmixed and the results were validated with the corresponding ground data. To study an urbanized landscape, a clear sky Landsat-5 TM data were unmixed and validated with coincident World View-2 abundance maps (of 2 m spatial resolution) for an area of San Francisco, California, USA. The results were evaluated using descriptive statistics, correlation coefficient, RMSE, probability of success, boxplot and bivariate distribution function. Finally, FCLS was used for sub-pixel land cover analysis of the monthly WELD (Wen-enabled Landsat data) repository from 2008 to 2011 of North America. The abundance maps in conjunction with DMSP-OLS nighttime lights data were used to extract the urban land cover features and analyze their spatial-temporal growth.

Landsat Satellites↗

Compton polarimeter during the Qweak experiment

Compton Polarimeter was successfully used during the Q$_{weak}$ experiment from Nov 2010 to May 2012 in Hall C, Jefferson Lab. In this newly installed Compton polarimeter the electron beam collides with green laser stored in a low gain Fabry-Perot Cavity; the scattered electrons are detected in 4 planes of a novel diamond micro strip detector while the back scattered photons are detected in four crystals of PbWO$_4$ detector. The analyzing power is computed using a Geant4 Monte Carlo simulation to match the real detector energy resolution and non-linearity in the response of the detector are folded in. The Compton polarimeter has achieved the design goals of 1\% statistical uncertainty per hour and we expect to achieve 1\% systematic uncertainty. We have done specific tests to understand various properties of the electron and photon detectors. We will share our study of the contributions to the electron and photon detector's systematic uncertainty along with the preliminary polarization results.

Narayan, Amrendra↗

Two-color pyrometry based flame to fuel surface radiative heat flux diagnostic using flamelets

In this study, a non-intrusive two-color pyrometry based radiative heat flux diagnostic is presented that accounts for the self-absorption effects. This diagnostic relies on looking up the flame temperature and soot volume fraction using the ratio of color intensities obtained from digital cameras and numerical solutions of 1D steady-state diffusion flames. Virtual two color pyrometry of computed flame show a unique one-to-one mapping of strained flames with intensity ratio; thereby allowing flames to be indexed and creation of a two-color flamelet manifold (TCFM). The TCFM is a function of two variables; the ratio of red to green intensities (flame index) and physical distance to the stoichiometric surface defined using a level-set function. The developed diagnostic is used in upward flame spread experiments where an approximate flame hull reconstruction method is applied to provide the 3D level-set function and intensity ratio for use of TCFM to provide radiative absorption properties for 3D ray tracing. Accounting for self- absorption effects provide improved soot volume fraction estimates. In conclusion, the radiative heat flux estimates are found to be in very good agreement with literature and improved results are obtained near the flame base.

42 ENGINEERING↗

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Strategies for designing biocatalysts with new functions

Enzymes can be optimized to accelerate chemical transformations via a range of methods. In this review, we showcase how protein engineering and computational design techniques can be interfaced to develop highly efficient and selective biocatalysts.

09 BIOMASS FUELS↗

Approximating accelerator impedances with resonator networks

It is common in the accelerator community to use the impedance of accelerator components to describe wake interactions in the frequency domain. However, it is often desirable to understand such wake interactions in the time domain in a general manner for excitations that are not necessarily Gaussian in nature. Here, the conventional method for doing this involves taking the inverse Fourier Transform of the component impedance, obtaining the Green's Function, and then convolving it with the desired excitation distribution. This method can prove numerically cumbersome, for a convolution integral must be evaluated for each individual point in time when the wake function is desired. An alternative to this method would be to compute the wake function analytically, which would sidestep the need for repetitive integration. Only a handful of cases, however, are simple enough for this method to be tenable. One of these cases is the case where the component in question is an RLC resonator, which has a closed-form analytical wake function solution. This means that a component which can be represented in terms of resonators can leverage this solution. As it happens, common network synthesis techniques may be used to map arbitrary impedance profiles to RLC resonator networks in a manner the accelerator community has yet to take advantage of. In this work, we will use Foster Canonical Resonator Networks and partial derivative descent optimization to develop a technique for synthesizing resonator networks that well approximate the impedances of real-world accelerator components. We will link this synthesis to the closed-form resonator wake function solution, giving rise to a powerful workflow that may be used to streamline beam dynamics simulations.

43 PARTICLE ACCELERATORS↗

Multiple scattering theory for dense plasmas

Dense plasmas occur in stars, giant planets, and in inertial fusion experiments. Accurate modeling of the electronic structure of these plasmas allows for prediction of material properties that can in turn be used to simulate these astrophysical objects and terrestrial experiments. But modeling them remains a challenge. Here we explore the Korringa-Kohn-Rostoker Green's function (KKR-GF) method for this purpose. We find that it is able to predict equation of state in good agreement with other state-of-the-art methods, where they are accurate and viable. In addition, it is shown that the computational cost does not significantly change with temperature, in contrast with other approaches. Moreover, the method does not use pseudopotentials—core states are calculated self consistently. We conclude that KKR-GF is a very promising method for dense plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NMSBA Sustainable Engineering (Final Report)

This report summarizes the guidance provided to Sustainable Engineering to help them learn about equation-oriented optimization and the Sandia-developed software packages Pyomo and IDAESPSE. This was a short 10-week project (October 2021 – December 2021) and the goal was to help the company learn about the IDAES framework and how it could be used for their future projects. The company submitted an SBIR proposal related to developing a green ammonia process model with IDAES and if that proposal is successful this NMSBA project could lead to future collaboration opportunities.

42 ENGINEERING↗

Holographic Helmet-Mounted Display Unit

Helmet-mounted display unit designed for use in testing innovative concepts for display of information to aircraft pilots. Operates in conjunction with computers generating graphical displays. Includes two ocular subunits containing miniature cathoderay tubes and optics providing 40 degrees vertical, 50 degrees horizontal field of view to each eye, with or without stereopsis. In future color application, each ocular subunit includes trichromatic holographic combiner tuned to red, green, and blue wavelengths of phosphors used in development of miniature color display devices.

Burley, James R., II↗

Code verification for practically singular equations

We report the method-of-moments implementation of the electric-field integral equation (EFIE) yields many code-verification challenges due to the various sources of numerical error and their possible interactions. Matters are further complicated by singular integrals, which arise from the presence of a Green's function. To address these singular integrals, an approach is presented in wherein both the solution and Green's function are manufactured. Because the arising equations are poorly conditioned, they are reformulated as a set of constraints for an optimization problem that selects the solution closest to the manufactured solution. In this paper, we demonstrate how, for such practically singular systems of equations, computing the truncation error by inserting the exact solution into the discretized equations cannot detect certain orders of coding errors. On the other hand, the discretization error from the optimal solution is a more sensitive metric that can detect orders less than those of the expected convergence rate.

97 MATHEMATICS AND COMPUTING↗

Switching limits of top-gated carbon nanotube field-effect transistors

The performance limits of carbon nanotube field-effect transistors (CNFETs) based on a recently reported process are studied by computational techniques at temperatures between 300 K and 4 K. Here, the impact of band-to-band tunneling (BTBT) and source-to-drain tunneling (SDT) is examined for devices with varying gate length through the use of simulations based on the non-equilibrium Green’s function (NEGF) formalism, and calibrated to measurements. Additionally, the case of junctionless chemical doping profiles is analyzed in contrast to electrostatic doping recently reported for test structures. The switching limits of CNFETs are further explored for devices based on carbon nanotubes (CNTs) with more favorable electronic structures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗