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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Simulated annealing in networks for computing possible arrangements for red and green cones

Attention is given to network models in which each of the cones of the retina is given a provisional color at random, and then the cones are allowed to determine the colors of their neighbors through an iterative process. A symmetric-structure spin-glass model has allowed arrays to be generated from completely random arrangements of red and green to arrays with approximately as much disorder as the parafoveal cones. Simulated annealing has also been added to the process in an attempt to generate color arrangements with greater regularity and hence more revealing moirepatterns than than the arrangements yielded by quenched spin-glass processes. Attention is given to the perceptual implications of these results.

Ahumada, Albert J., Jr.↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
DOE OSTI2024

Computational Advances in Ionic Liquid Applications for Green Chemistry: A Critical Review of Lignin Processing and Machine Learning Approaches

The valorization and dissolution of lignin using ionic liquids (ILs) is critical for developing sustainable biorefineries and a circular bioeconomy. This review aims to critically assess the current state of computational and machine learning methods for understanding and optimizing IL-based lignin dissolution and valorization processes reported since 2022. The paper examines various computational approaches, from quantum chemistry to machine learning, highlighting their strengths, limitations, and recent advances in predicting and optimizing lignin-IL interactions. Key themes include the challenges in accurately modeling lignin’s complex structure, the development of efficient screening methodologies for ionic liquids to enhance lignin dissolution and valorization processes, and the integration of machine learning with quantum calculations. These computational advances will drive progress in IL-based lignin valorization by providing deeper molecular-level insights and facilitating the rapid screening of novel IL-lignin systems.

09 BIOMASS FUELS↗
DOE OSTI2019

Locally self-consistent embedding approach for disordered electronic systems

We introduce an embedding scheme for the locally self-consistent method to study disordered electron systems. We test this approach in a tight-binding basis and apply it to the single band Anderson model. The local interaction zone is used to efficiently compute the local Green's function of a supercell embedded into a local typical medium. We find a quick convergence as the size of the local interaction zone which reduces the computational costs as expected. This method captures the Anderson localization transition and accurately predicts the critical disorder strength. The current work opens the path towards the development of a typical medium embedding scheme for the O(N) multiple scattering methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the fractional Laplacian of variable order

We present a novel definition of variable-order fractional Laplacian on R n based on a natural generalization of the standard Riesz potential. Our definition holds for values of the fractional parameter spanning the entire open set (0, n/2). We then discuss some properties of the fractional Poisson’s equation involving this operator and we compute the corresponding Green’s function, for which we provide some instructive examples for specific problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗
DOE OSTI2021

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Improvement of Numerical Methods in Pronghorn

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Geomechanical response due to nonisothermal fluid injection into a reservoir

The geomechanical response of a porous reservoir due to injection of fluid can result from a complex interplay between the changes in porepressure and temperature near the wellbore. As a result, predictions are usually made using either simplified analytical models, which may apply unrealistic assumptions in order to produce a tractable model, or detailed numerical simulations that can be computationally expensive.

58 GEOSCIENCES↗
DOE OSTI2019

Volume Flow Rate Estimation for Small Explosions at Mt. Etna, Italy, From Acoustic Waveform Inversion

Rapid assessment of the volume and the rate at which gas and pyroclasts are injected into the atmosphere during volcanic explosions is key to effective eruption hazard mitigation. Here, we use data from a dense infrasound network deployed in 2017 on Mt. Etna, Italy, to estimate eruptive volume flow rates (VFRs) during small gas-and-ash explosions. We use a finite-difference time-domain approximation to compute the acoustic Green's functions and perform a full waveform inversion for a multipole source, combining monopole and horizontal dipole terms. The inversion produces realistic estimates of VFR, on the order of 4 × 10 4 m 3 /s and well-defined patterns of source directivity. This is the first application of acoustic waveform inversion at Mt. Etna. Our results demonstrate that acoustic waveform inversion is a mature and robust tool for assessment of source parameters and holds potential as a tool to provide rapid estimates of VFR in near real time.

58 GEOSCIENCES↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2018

Thermal conductivity of InN with point defects from first principles

We present ab initio calculations of thermal conductivity of InN with vacancies and substitutional defects using a full solution of the Peierls-Boltzmann transport equation. Our parameter-free calculations are in good agreement with experimental measurements demonstrating the predictive power of this approach. Phonon-defect scattering rates are computed from a Green's function methodology that is nonperturbative and includes interatomic force constant variance induced near the defects. Restricting calculations to first-order perturbation approaches can overestimate optic phonon scattering rates by nearly three orders of magnitude. On the other hand, neglecting the force variance weakens the scattering rates by about an order of magnitude, mostly in the low-frequency region below 2 THz. Here, this work elucidates important properties of phonon-defect scattering in thermal transport and demonstrates the predictive power of the coupling of Peierls-Boltzmann transport, Green's function methods, and density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1978

Information extraction from digital images of the earth and the planets

An overview is presented of recent developments at the JPL Image Processing Laboratory which emphasize the utilization of a digital computer to automate the process of information extraction from digital imagery. Consideration is given to: (1) the analysis of Viking Orbiter stereo imagery to determine elevation profiles of the Mars surface, (2) the use of Viking Lander stereo imagery to determine nonhazardous surface sample acquisition strategies, (3) the correlation of Landsat imagery with geographically referenced cultural data to determine land use trends, and (4) the generation of mosaics using digital computer techniques.

Green, W. B.↗
NASA NTRS1989

A 3D visualization system for molecular structures

The properties of molecules derive in part from their structures. Because of the importance of understanding molecular structures various methodologies, ranging from first principles to empirical technique, were developed for computing the structure of molecules. For large molecules such as polymer model compounds, the structural information is difficult to comprehend by examining tabulated data. Therefore, a molecular graphics display system, called MOLDS, was developed to help interpret the data. MOLDS is a menu-driven program developed to run on the LADC SNS computer systems. This program can read a data file generated by the modeling programs or data can be entered using the keyboard. MOLDS has the following capabilities: draws the 3-D representation of a molecule using stick, ball and ball, or space filled model from Cartesian coordinates, draws different perspective views of the molecule; rotates the molecule on the X, Y, Z axis or about some arbitrary line in space, zooms in on a small area of the molecule in order to obtain a better view of a specific region; and makes hard copy representation of molecules on a graphic printer. In addition, MOLDS can be easily updated and readily adapted to run on most computer systems.

Green, Terry J.↗
NASA NTRS1990

The cleanroom case study in the Software Engineering Laboratory: Project description and early analysis

This case study analyzes the application of the cleanroom software development methodology to the development of production software at the NASA/Goddard Space Flight Center. The cleanroom methodology emphasizes human discipline in program verification to produce reliable software products that are right the first time. Preliminary analysis of the cleanroom case study shows that the method can be applied successfully in the FDD environment and may increase staff productivity and product quality. Compared to typical Software Engineering Laboratory (SEL) activities, there is evidence of lower failure rates, a more complete and consistent set of inline code documentation, a different distribution of phase effort activity, and a different growth profile in terms of lines of code developed. The major goals of the study were to: (1) assess the process used in the SEL cleanroom model with respect to team structure, team activities, and effort distribution; (2) analyze the products of the SEL cleanroom model and determine the impact on measures of interest, including reliability, productivity, overall life-cycle cost, and software quality; and (3) analyze the residual products in the application of the SEL cleanroom model, such as fault distribution, error characteristics, system growth, and computer usage.

Green, Scott↗
NASA NTRS1993

Application of the multigrid solution technique to hypersonic entry vehicles

A multigrid solution procedure has been incorporated in a version of the Langley Aerothermodynamic Upwind Relaxation Algorithm. The multigrid scheme is based on the Full Approximation Storage approach and uses Full Multigrid to obtain a well defined fine mesh starting solution. Predictions were obtained using standard transfer operators and a 'V-cycle' was used to control grid sequencing. Computed hypersonic flow solutions compared with experimental data for a 15 degree sphere cone, blended-wing body, and shuttle-like geometries are presented. It is shown that the algorithm accurately predicts heating rates, and when compared with the single grid algorithm computes solutions in one-third the computational time.

Greene, Francis A.↗
NASA NTRS1994

LET spectra of fragmented 600-MeV/A iron beam on aluminum and polyethylene: A preliminary analysis

An iron beam experiment recently conducted at the Lawrence Berkeley Laboratory Bevalac by Benton et al. provides an opportunity for verifying the new Green's function computer code (GRNTRN) and assessing the related nuclear database. The iron beam with 600 MeV/A at extraction traversed a series of beam transport elements, lead foil, and several triggering devices before impacting the target. Of these, only the 2.24 g/sq cm lead foil and target are considered in the transport analysis with an assumed (inferred from calibration) beam energy of 557 MeV/A. A thin layer of CR-39 plastic nuclear track detector (PNTD) was placed in front of the target to monitor the incident flux and a stack of four PNTD's placed behind to measure the linear energy transfer (LET) distribution of the transported beam. Test data are analyzed for three separate targets: 2 g/sq cm aluminum, 5-cm polyethylene, and 8-cm polyethylene. The two-layer GRNTRN results were mapped into the detector response function for comparison with the measured LET spectra. Reasonable agreement is obtained. Future research and analysis can be improved by using a more accurate isotope set or including other important media which significantly alter the beam. The assumption that the fragmentation cross sections are too small for aluminum is inconclusive because the production of fragments by 2 g/sq cm aluminum can be more significantly affected by the fragment contribution originating from the media in front of the target than that currently considered in the analysis.

Shinn, Judy L.↗