Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HUBBARD MODEL”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Optimized Jupiter, Saturn, and Uranus interior models

A method for inverting gravitational data is the basis for the present models of Saturn, Jupiter, and Uranus, which have been computed to the third-order in rotational disturbance to the total potential and yield exact matches for recent, highly-accurate determinations of gravitational harmonics for these planets. These results indicate that a major He depletion has either not occurred in the outer layers of Jupiter or of Saturn, or that such a depletion's effect on the density profile has been masked by the presence of denser components. Considerable uncertainty characterizes Uranian models; the one presently preferred has substantial heavy-element enhancement at depth and little ice-rock separation.

Hubbard, W. B.↗
DOE OSTI2026

Band-mixing effects in one-dimensional charge transfer insulators

The low-energy properties of transition metal oxides (TMOs) are governed by the electrons occupying strongly correlated $d$-orbitals that are hybridized with surrounding ligand oxygen $p$ orbitals to varying degrees. Their physics is thus established by a complex interplay between the transition-metal (TM)-ligand hopping $t$, charge transfer energy $\Delta_\mathrm{CT}$, and on-site TM Hubbard repulsion $U$. Here, we study the spectral properties of a one-dimensional (1D) analog of such a $pd$ system, with alternating TM $d$ and ligand anion $p$ orbitals situated along a chain. Using the density matrix renormalization group method, we study the model's single-particle spectral function, x-ray absorption spectrum, and dynamical spin structure factor as a function of $\Delta_\mathrm{CT}$ and $U$. In particular, we present results spanning from the Mott insulating ($\Delta_\mathrm{CT} > U$) to negative charge transfer regime $\Delta_\mathrm{CT} < 0$ to gain a better understanding of the ground and momentum-resolved excited state properties of these different regimes. Our results can guide new studies on TMOs that seek to situate them within the Mott-Hubbard/charge transfer insulator classification scheme.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

BASIN-3D: A brokering framework to integrate diverse environmental data

Diverse observational and simulation datasets are needed to understand and predict complex ecosystem behavior over seasonal to decadal and century time-scales. Integration of these datasets poses a major barrier towards advancing environmental science, particularly due to differences in the structure and formats of data provided by various sources. Here, we describe BASIN-3D (Broker for Assimilation, Synthesis and Integration of eNvironmental Diverse, Distributed Datasets), a data integration framework designed to dynamically retrieve and transform heterogeneous data from different sources into a common format to provide an integrated view. BASIN-3D enables users to adopt a standardized approach for data retrieval and avoid customizations for the data type or source. We demonstrate the value of BASIN-3D with two use cases that require integration of data from regional to watershed spatial scales. The first application uses the BASIN-3D Python library to integrate time-series hydrological and meteorological data to provide standardized inputs to analytical and machine learning codes in order to predict the impacts of hydrological disturbances on large river corridors of the United States. The second application uses the BASIN-3D Django framework to integrate diverse time-series data in a mountainous watershed in East River, Colorado, United States to enable scientific researchers to explore and download data through an interactive web portal. Thus, BASIN-3D can be used to support data integration for both web-based tools, as well as data analytics using Python scripting and extensions like Jupyter notebooks. The framework is expected to be transferable to and useful for many other field and modeling studies.

Varadharajan, C↗
NASA NTRS1985

Theoretical quasar emission-line ratios. VII - Energy-balance models for finite hydrogen slabs

The present energy balance calculations for finite, isobaric, hydrogen-slab quasar emission line clouds incorporate probabilistic radiative transfer (RT) in all lines and bound-free continua of a five-level continuum model hydrogen atom. Attention is given to the line ratios, line formation regions, level populations and model applicability results obtained. H lines and a variety of other considerations suggest the possibility of emission line cloud densities in excess of 10 to the 10th/cu cm. Lyman-beta/Lyman-alpha line ratios that are in agreement with observed values are obtained by the models. The observed Lyman/Balmer ratios can be achieved with clouds whose column depths are about 10 to the 22nd/sq cm.

Hubbard, E. N.↗
DOE OSTI2020

Calculating and Characterizing the Charge Distributions in Solids

Accurate estimation of the partial atomic charges on metal centers is useful for understanding electronic and catalytic properties of materials. However, different methods of calculating these charges may give quite different results; this issue has been more widely studied for molecules than for solids. Here we study the charges on the metal centers of a test set of 18 solids containing transition metals by using density functional theory with several density functionals (PBE, PBE+U, TPSS, revTPSS, HLE17, revM06-L, B3LYP, B3LYP*, and other exchange-modified B3LYP functionals) and four charge models (Bader, Hirshfeld, CM5, and DDEC6). The test set contains 12 systems with nonmagnetic metal centers (eight metal oxides (MO 2 ), two metal sulfides (MS 2 ), and two metal selenides (MSe 2 )) and six ferromagnetic transition metal complexes. Our study shows that, among the four types of charges, Bader charges are the highest and Hirshfeld charges are the lowest for all the systems, regardless of the functional being used. The CM5 charges are bigger than DDEC6 charges for MX 2 with M = Ti or Mo and X = S or Se, but for the other 14 cases they are lower. We found that the most of the systems are sensitive to the Hubbard U parameters in PBE+U and to the percentage X of Hartree–Fock exchange in exchange-modified B3LYP; as we increase U or X, the charges on the metal atoms in MX 2 increase steadily. Testing different density functionals shows charges calculated with higher Hubbard U parameters in PBE+U are comparable to B3LYP (with 20% Hartree–Fock exchange). Among four meta-GGA functionals studied, the charges with HLE17 have the closest agreement with B3LYP. In conclusion, the variation of charges with choice of charge model is greater than the variation with choice of density functional.

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

Local-scale heterogeneity of soil thermal dynamics and controlling factors in a discontinuous permafrost region

In permafrost regions, the strong spatial and temporal variability in soil temperature cannot be explained by the weather forcing only. Understanding the local heterogeneity of soil thermal dynamics and their controls is essential to understand how permafrost systems respond to climate change and to develop process-based models or remote sensing products for predicting soil temperature. In this study, we analyzed soil temperature dynamics and their controls in a discontinuous permafrost region on the Seward Peninsula, Alaska. We acquired one-year temperature time series at multiple depths (at 5 or 10 cm intervals up to 85 cm depth) at 45 discrete locations across a 2.3 km 2 watershed. We observed a larger spatial variability in winter temperatures than that in summer temperatures at all depths, with the former controlling most of the spatial variability in mean annual temperatures. We also observed a strong correlation between mean annual ground temperature at a depth of 85 cm and mean annual or winter season ground surface temperature across the 45 locations. We demonstrate that soils classified as cold, intermediate, or warm using hierarchical clustering of full-year temperature data closely match their co-located vegetation (graminoid tundra, dwarf shrub tundra, and tall shrub tundra, respectively). We show that the spatial heterogeneity in soil temperature is primarily driven by spatial heterogeneity in snow cover, which induces variable winter insulation and soil thermal diffusivity. These effects further extend to the subsequent summer by causing variable latent heat exchanges. Finally, we discuss the challenges of predicting soil temperatures from snow depth and vegetation height alone by considering the complexity observed in the field data and reproduced in a model sensitivity analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Complementary Study of Radiative Heat Transfer and Flow Physics from Moderate-scale Hydrocarbon Pool Fire Simulations

As part of the Advanced Simulation and Computing Verification and Validation (ASCVV) program, a 0.3-m diameter hydrocarbon pool fire with multiple fuels was modeled and simulated. In the study described in this report, systematic examination was performed on the radiation model used in a series of coupled Fuego/Nalu simulations. A calibration study was done with a medium-scale methanol pool fire and the effect of calibration traced throughout the radiation model. This analysis provided a more detailed understanding of the effect of radiation model parameters on each other and on other quantities in the simulations. Heptane simulation results were also examined using this approach and possible areas for further improvement of the models were identified. The effect of soot on radiative losses was examined by comparing heptane and methanol results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1985

Theoretical Thermodynamics of Mixtures at High Pressures

The development of an understanding of the chemistry of mixtures of metallic hydrogen and abundant, higher-z material such as oxygen, carbon, etc., is important for understanding of fundamental processes of energy release, differentiation, and development of atmospheric abundances in the Jovian planets. It provides a significant theoretical base for the interpretation of atmospheric elemental abundances to be provided by atmospheric entry probes in coming years. Significant differences are found when non-perturbative approaches such as Thomas-Fermi-Dirac (TFD) theory are used. Mapping of the phase diagrams of such binary mixtures in the pressure range from approx. 10 Mbar to approx. 1000 Mbar, using results from three-dimensional TFD calculations is undertaken. Derivation of a general and flexible thermodynamic model for such binary mixtures in the relevant pressure range was facilitated by the following breakthrough: there exists an accurate nd fairly simple thermodynamic representation of a liquid two-component plasma (TCP) in which the Helmholtz free energy is represented as a suitable linear combination of terms dependent only on density and terms which depend only on the ion coupling parameter. It is found that the crystal energies of mixtures of H-He, H-C, and H-O can be satisfactorily reproduced by the same type of model, except that an effective, density-dependent ionic charge must be used in place of the actual total ionic charge.

Hubbard, W. B.↗
DOE OSTI2023

Ab initio modeling and thermodynamics of hydrated plutonium oxalates

An ab initio study on the plutonium oxalate hydrates: Pu 2 (C 2 O 4 ) 3 ·10H 2 O and Pu(C 2 O 4 ) 2 ·6H 2 O, using PBE exchange-correlation with D3 dispersion correction and Hubbard correction for the plutonium atoms, was performed and compared to experimental vibrational spectral and thermodynamic property values. Here we demonstrated that this technique can accurately predict the experimental infrared spectra Pu(III) oxalate hydrate, as well as the Raman peak of PuO 2 (used to calculate the thermodynamic properties of the oxalates). For Pu(IV) oxalate hydrate we found that our predicted structure agreed qualitatively with PXRD measurements, the only available experimental determination of the structure. Using this method at standard temperature and pressure, we predicted standard enthalpies of formation of -6,755 kJ mol -1 and -3,923 kJ mol -1 and standard Gibbs free energy of formation of -5,899 kJ mol -1 and -3,386 kJ mol -1 for Pu 2 (C 2 O 4 ) 3 ·10H 2 O and Pu(C 2 O 4 ) 2 ·6H 2 O, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic Structure of InAs and InSb Surfaces: Density Functional Theory and Angle-Resolved Photoemission Spectroscopy

The electronic structure of surfaces plays a key role in the properties of quantum devices. However, surfaces are also the most challenging to simulate and engineer. Here, in this work, the electronic structure of InAs(001), InAs(111), and InSb(110) surfaces is studied using a combination of density functional theory (DFT) and angle-resolved photoemission spectroscopy (ARPES). Large-scale first principles simulations are enabled by using DFT calculations with a machine-learned Hubbard U correction [npj Comput. Mater. 6, 180 (2020)]. To facilitate direct comparison with ARPES results, a “bulk unfolding” scheme is implemented by projecting the calculated band structure of a supercell surface slab model onto the bulk primitive cell. For all three surfaces, a good agreement is found between DFT calculations and ARPES. For InAs(001), the simulations clarify the effect of the surface reconstruction. Different reconstructions are found to produce distinctive surface states, which may be detected by ARPES with low photon energies. For InAs(111) and InSb(110), the simulations help elucidate the effect of oxidation. Owing to larger charge transfer from As to O than from Sb to O, oxidation of InAs(111) leads to significant band bending and produces an electron pocket, whereas oxidation of InSb(110) does not. The combined theoretical and experimental results may inform the design of quantum devices based on InAs and InSb semiconductors, for example, topological qubits utilizing the Majorana zero modes.

42 ENGINEERING↗
DOE OSTI2021

Revolutionizing observations and predictability of Arctic system dynamics through next-generation dense, heterogeneous and intelligent wireless sensor networks with embedded AI

The Arctic region is a complex and dynamic system, where changes in temperature and partitioning between water in solid, liquid and vapor phase are fundamentally reshaping hydro-biogeochemical fluxes from the bedrock to the atmosphere. The increasing temperature and frequency of extreme heat, precipitation and fire events (Wang et al., 2017; Meredith et al., 2019) have critical consequences for the Arctic ecosystems and people, as well as for the global climate system. Although many observational and remote sensing activities have been on-going, detailed ecosystem processes – particularly related to controls on and impacts of permafrost thaw – are difficult to detect or measure. To improve the understanding and prediction of these extreme events and their consequences, we need a revolution in the way we collect ground-based observations from hillslope to pan-Arctic scales and couple them with satellite imagery and models to deliver rich, actionable, and scale-appropriate data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Transfer Tube Material Properties Obtained from Test Data

A series of tension and compression tests were performed on typical transfer tubes made of 304L corrosion-resistant steel. The results were compiled and converted to true stress and true plastic strain for comparison to numerical models. The parameters of a Johnson-Cook constitutive model and fracture criterion were calibrated to mimic the test data. Comparisons were made between the converted data and the empirical formulas in the Johnson-Cook model. Simulations were performed to evaluate the implementation of the formula in finite element analysis, and the parameters were adjusted to improve the correlation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A hybrid data–model approach to map soil thickness in mountain hillslopes

Abstract. Soil thickness plays a central role in the interactions between vegetation, soils, and topography, where it controls the retention and release of water, carbon, nitrogen, and metals. However, mapping soil thickness, here defined as the mobile regolith layer, at high spatial resolution remains challenging. Here, we develop a hybrid model that combines a process-based model and empirical relationships to estimate the spatial heterogeneity of soil thickness with fine spatial resolution (0.5 m). We apply this model to two aspects of hillslopes (southwest- and northeast-facing, respectively) in the East River watershed in Colorado. Two independent measurement methods – auger and cone penetrometer – are used to sample soil thickness at 78 locations to calibrate the local value of unconstrained parameters within the hybrid model. Sensitivity analysis using the hybrid model reveals that the diffusion coefficient used in hillslope diffusion modeling has the largest sensitivity among all input parameters. In addition, our results from both sampling and modeling show that, in general, the northeast-facing hillslope has a deeper soil layer than the southwest-facing hillslope. By comparing the soil thickness estimated between a machine-learning approach and this hybrid model, the hybrid model provides higher accuracy and requires less sampling data. Modeling results further reveal that the southwest-facing hillslope has a slightly faster surface soil erosion rate and soil production rate than the northeast-facing hillslope, which suggests that the relatively less dense vegetation cover and drier surface soils on the southwest-facing slopes influence soil properties. With seven parameters in total for calibration, this hybrid model can provide a realistic soil thickness map with a relatively small amount of sampling dataset comparing to machine-learning approach. Integrating process-based modeling and statistical analysis not only provides a thorough understanding of the fundamental mechanisms for soil thickness prediction but also integrates the strengths of both statistical approaches and process-based modeling approaches.

58 GEOSCIENCES↗
NASA NTRS1981

Observations of electron gyroharmonic waves and the structure of the Io torus

Narrow-banded emissions observed by the planetary radio astronomy experiment on the Voyager 1 spacecraft as it traversed the Io plasma torus are discussed. It is found that the waves occur between harmonics of the electron gyrofrequency; they are Jovian analogue of electrostatic emissions observed and theoretically studied for the terrestrial magnetosphere. It is noted that the observed frequencies always include the component near the upper hybrid resonant frequency but that the distribution of the other observed emissions varies in a systematic way with position in the torus. A detailed discussion of the observations is given. Also included is a refined model of the electron density variation, based on identification of the upper hybrid resonant frequency line.

Birmingham, T. J.↗
DOE OSTI2024

Magnetism and finite-temperature effects in UZr 2 : A density functional theory analysis

The structure and the thermophysical properties of d-UZr 2 are investigated using 0 K density functional theory and ab initio molecular dynamics (AIMD). Modeling the true paramagnetic state of this intermetallic compound has been challenging using first-principles calculations. For the first time, we find that the generalized gradient approximation method without applying an on-site Coulomb interaction term (Hubbard U) can result in a ground state that is antiferromagnetic (AFM). We believe that this weak AFM ground state is the closest to the real paramagnetic state. We found that structure optimization at finite temperatures using AIMD is necessary to achieve this ground state instead of the non-magnetic and ferromagnetic states previously reported in the literature. Our findings indicate that applying the Hubbard U on uranium f-orbitals in this metallic system is unnecessary and not recommended, as it leads to a large overestimation of the volume and introduces an unphysical strong spin polarization. Finally, our approach results in atomic volume, thermal expansion, and heat capacities that have strong agreement with experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elucidating the impact of fiber source on polypropylene/hemp composite performance for the automotive industry

Given their high strength-to-weight ratio, there is an ever-increasing volume of plastics being used in the automotive industry as plastics aid in the charge to lightweight vehicles for improved fuel efficiency. However, these plastics are often landfilled at their end-of-life, which has given rise to the demand for sustainable materials and waste management alternatives compared to purely synthetic systems. Natural fiber composites have been explored as a viable material option to reduce the environmental impact of plastic use in automobiles while simultaneously ensuring the part performance is not sacrificed. Herein, we explored the use of hemp/polypropylene (PP) composites in which US-sourced hemp is compared to internationally sourced and industrially available hemp. There appears to be a minimal impact on the composite properties regardless of fiber sourcing, and the addition of a natural filler to the PP matrix results in up to a 367% increase in Young's modulus, 126% increase in heat deflection temperature, and comparable water uptake performance. It should be noted that the natural filler addition does increase density by up to 13% due to the higher density of natural fibers compared to the low-density PP matrix. A modified rule of mixtures calculation revealed that the composite materials produced in this study demonstrated good agreement with analytical modeling. Finally, a screening analysis was performed exploring the transportation of hemp fibers, and the results build a strong case for regionalized manufacturing of automotive parts.

Hubbard, Amber M. [Oak Ridge National Laboratory (↗
DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗
NASA NTRS1973

Lunar rock types - The role of plagioclase in non-mare and highland rock types

Some nonmare and highland rock types (14310 type KREEP and very high Al2O3 basalts) have the internal chemical variations expected for a plagioclase-liquid system. The observed Eu variations in these rock types suggest a D(Eu) 1/p of 0.6 to 0.7. The Sr variations suggest a D(Sr) 1/p of about 0.6, with values as low as 0.35 suggested for some materials from sample 14063. Common Apollo 14 KREEP and Apollo 15 KREEP do not show internal Sr, Eu, Al2O3 variations consistent with the D(Eu, Sr) 1/p values derived for 14310 type KREEP. Major element and experimental data indicate that olivine or pyroxene is a large, perhaps dominant, controller of chemical variations within common Apollo 14 KREEP. The application of these distribution coefficients to pure anorthosites like 15415 yields the model dependent conclusion that the silicate liquids with which such anorthosites may have been chemically equilibrated have not yet been analyzed and perhaps not directly sampled.

Hubbard, N. J.↗