Search NASA⌕ Search

SEARCH · Search NASA

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Insights into cellulase-lignin non-specific binding revealed by computational redesign of the surface of green fluorescent protein

Biological-mediated conversion of pretreated lignocellulosic biomass to biofuels and biochemicals is a promising avenue towards energy sustainability. However, a critical impediment to the commercialization of cellulosic biofuel production is the high cost of cellulase enzymes needed to deconstruct biomass into fermentable sugars. One major factor driving cost is cellulase adsorption and inactivation in the presence of lignin, yet we currently have a poor understanding of the protein structure-function relationships driving this adsorption. In this work, we have systematically investigated the role of protein surface potential on lignin adsorption using a model monomeric fluorescent protein. We have designed and experimentally characterized 16 model protein variants spanning the physiological range of net charge (-24 to +16 total charges) and total charge density (0.28 to 0.40 charges per sequence length) typical for natural proteins. Protein designs were expressed, purified, and subjected to in silico and in vitro biophysical measurements to evaluate the relationship between protein surface potential and lignin adsorption properties. The designs were comparable to model fluorescent protein in terms of thermostability and heterologous expression yield, although the majority of the designs unexpectedly formed homodimers. Protein adsorption to lignin was studied at two different temperatures using Quartz Crystal Microbalance with Dissipation Monitoring and a subtractive mass balance assay. We found a weak correlation between protein net charge and protein-binding capacity to lignin. No other single characteristic, including apparent melting temperature and 2nd virial coefficient, showed correlation with lignin binding. Analysis of an unrelated cellulase dataset with mutations localized to a family I carbohydrate-binding module showed a similar correlation between net charge and lignin binding capacity. Altogether, our study provides strategies to identify highly active, low lignin-binding cellulases by either rational design or by computational screening genomic databases.

09 BIOMASS FUELS↗
DOE OSTI2019

Computational screening, synthesis and testing of metal–organic frameworks with a bithiazole linker for carbon dioxide capture and its green conversion into cyclic carbonates

Computational crystal construction algorithms were used to create twelve metal–organic frameworks containing a newly synthesized [2,2'-bithiazole]-5,5'-dicarboxylic acid (H 2 TzTz) spacer and assorted transition metal nodes. Among the twelve structures, the zirconium-based MOF of general formula [Zr 6 O 4 (OH) 4 (TzTz) 6 ] (1) was found to be the best candidate for carbon dioxide uptake, as judged from the results of the grand canonical Monte Carlo (GCMC) simulations of CO 2 adsorption isotherms. Guided by the simulation results, 1 was synthesized in the laboratory and thoroughly characterized. 1 is isoreticular to its bithiophene and bis(benzene) (UiO-67) analogues; it crystallizes in the cubic Pn$\bar{3}$ with combining macron] space group with fcu topology, and it features octahedral [Zr 6 ] nodes connected by twelve carboxylate groups from six bridging TzTz 2- spacers. It is a predominantly microporous material (micropore volume = 84% of the total pore volume), with a BET area of 840 m 2 g -1 and a maximum CO 2 uptake at ambient pressure of 2.3 mmol g -1 (10.0 wt%) or 1.7 mmol g -1 (7.5 wt%) at 273 or 298 K, respectively. The CO 2 affinity (isosteric heat of adsorption Qst = 18.7 kJ mol -1 ; CO 2 /N 2 Henry selectivity = 10; CO 2 /N 2 IAST selectivity = 8.4) is similar to that of its bithiophene analogue. After partial removal of solvent (activation), 1 was tested as a heterogeneous catalyst in the reaction of CO 2 with epoxides bearing a –CH 2 X pendant arm (X = Cl: epichlorohydrin; X = Br: epibromohydrin) to give the corresponding cyclic carbonates at T = 393 K and p CO 2 = 1 bar under green (solvent- and co-catalyst-free) conditions. A good conversion of 74% and a turnover frequency of 12.3 mmol (cyclic carbonate) per mmol Zr per h have been recorded with epibromohydrin as a substrate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2015

A WENO-solver combined with adaptive momentum discretization for the Wigner transport equation and its application to resonant tunneling diodes

We present a novel numerical scheme for the deterministic solution of the Wigner transport equation, especially suited to deal with situations in which strong quantum effects are present. The unique feature of the algorithm is the expansion of the Wigner function in local basis functions, similar to finite element or finite volume methods. This procedure yields a discretization of the pseudo-differential operator that conserves the particle density on arbitrarily chosen grids. The high flexibility in refining the grid spacing together with the weighted essentially non-oscillatory (WENO) scheme for the advection term allows for an accurate and well-resolved simulation of the phase space dynamics. A resonant tunneling diode is considered as test case and a detailed convergence study is given by comparing the results to a non-equilibrium Green's functions calculation. The impact of the considered domain size and of the grid spacing is analyzed. The obtained convergence of the results towards a quasi-exact agreement of the steady state Wigner and Green's functions computations demonstrates the accuracy of the scheme, as well as the high flexibility to adjust to different physical situations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Electromagnetic scaling functions within the Green's function Monte Carlo approach

We have studied the scaling properties of the electromagnetic response functions of 4 He and 12 C nuclei computed by the Green's function Monte Carlo approach, retaining only the one-body current contribution. Longitudinal and transverse scaling functions have been obtained in the relativistic and nonrelativistic cases and compared to experiment for various kinematics. The characteristic asymmetric shape of the scaling function exhibited by data emerges in the calculations in spite of the nonrelativistic nature of the model. The results are mostly consistent with scaling of zeroth, first, and second kinds. Our analysis reveals a direct correspondence between the scaling and the nucleon-density response functions. In conclusion, the scaling function obtained from the proton-density response displays scaling of the first kind, even more evidently than the longitudinal and transverse scaling functions

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2016

Characterization of the porous structures of the green body and sintered biomedical titanium scaffolds with micro-computed tomography

The present research was aimed at gaining an understanding of the porous structure changes from the green body through water leaching and sintering to titanium scaffolds. Micro-computed tomography (micro-CT) was performed to generate 3D models of titanium scaffold preforms containing carbamide space-holding particles and sintered scaffolds containing macro- and micro-pores. The porosity values and structural parameters were determined by means of image analysis. The result showed that the porosity values, macro-pore sizes, connectivity densities and specific surface areas of the titanium scaffolds sintered at 1200 °C for 3 h did not significantly deviate from those of the green structures with various volume fractions of the space holder. Titanium scaffolds with a maximum specific surface area could be produced with an addition of 60–65 vol% carbamide particles to the matrix powder. The connectivity of pores inside the scaffold increased with rising volume fraction of the space holder. The shrinkage of the scaffolds prepared with > 50 vol% carbamide space holder, occurring during sintering, was caused by the reductions of macro-pore sizes and micro-pore sizes as well as the thickness of struts. In conclusion, the final porous structural characteristics of titanium scaffolds could be estimated from those of the green body. - Highlights: •Porous structures of green body and sintered titanium scaffolds was studied. •Porous structures of both samples were quantitatively characterized with micro-CT. •Porous structures of scaffolds could be controlled from the green body. •Shrinkage mechanisms of titanium scaffolds during sintering was established.

36 MATERIALS SCIENCE↗
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↗