SEARCH · Engineering Papers

Results for “reference 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83

All sources NASA NTRS DOE OSTI

DOE OSTI2015

3D RISM theory with fast reciprocal-space electrostatics

The calculation of electrostatic solute-solvent interactions in 3D RISM (“three-dimensional reference interaction site model”) integral equation theory is recast in a form that allows for a computational treatment analogous to the “particle-mesh Ewald” formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Exchange repulsive potential adaptable for electronic structure changes during chemical reactions

Hybrid methods combining quantum mechanical (QM) and classical calculations are becoming important tools in chemistry. The popular approach to calculate the interaction between QM and classical calculations employs interatomic potentials. In most cases, the interatomic potential is constructed of an electrostatic (ES) potential and a non-ES potential. Because QM treatment is employed in the calculation of the ES potential, the electronic change can be considered in this ES potential. However, QM treatment of the non-ES potential is difficult because of high computational cost. To overcome this difficulty of evaluating the non-ES potential, we proposed an exchange repulsive potential as the main part of the non-ES potential on the basis of a QM approach. This potential is independent of empirical parameters and adaptable for electronic structure. We combined this potential with the reference interaction site model self-consistent field explicitly including spatial electron density distribution and successfully applied it to the chemical reactions in aqueous phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

A hybrid framework of first principles molecular orbital calculations and a three-dimensional integral equation theory for molecular liquids: Multi-center molecular Ornstein–Zernike self-consistent field approach

In this study, we reported the development of a new quantum mechanics/molecular mechanics (QM/MM)-type framework to describe chemical processes in solution by combining standard molecular-orbital calculations with a three-dimensional formalism of integral equation theory for molecular liquids (multi-center molecular Ornstein–Zernike (MC-MOZ) method). The theoretical procedure is very similar to the 3D-reference interaction site model self-consistent field (RISM-SCF) approach. Since the MC-MOZ method is highly parallelized for computation, the present approach has the potential to be one of the most efficient procedures to treat chemical processes in solution. Benchmark tests to check the validity of this approach were performed for two solute (solute water and formaldehyde) systems and a simple S{sub N}2 reaction (Cl{sup −} + CH{sub 3}Cl → ClCH{sub 3} + Cl{sup −}) in aqueous solution. The results for solute molecular properties and solvation structures obtained by the present approach were in reasonable agreement with those obtained by other hybrid frameworks and experiments. In particular, the results of the proposed approach are in excellent agreements with those of 3D-RISM-SCF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Theoretical study of chemical reactions in solution

Quantum chemical calculations in solution are becoming more and more important in chemistry. Reference interaction site model self-consistent field (RISM-SCF) is one of the powerful approaches to perform quantum chemical calculations in solution. In this work, we developed a new generation of RISM-SCF, where a robust fitting method was newly introduced. We applied the new method to tautomerization reaction of cytosine in aqueous phase. Our calculation reproduced experimentally obtained relative stabilities and relative free energies correctly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Design principles for high–pressure force fields: Aqueous TMAO solutions from ambient to kilobar pressures

Accurate force fields are one of the major pillars on which successful molecular dynamics simulations of complex biomolecular processes rest. They have been optimized for ambient conditions, whereas high-pressure simulations become increasingly important in pressure perturbation studies, using pressure as an independent thermodynamic variable. Here, we explore the design of non-polarizable force fields tailored to work well in the realm of kilobar pressures – while avoiding complete reparameterization. Our key is to first compute the pressure-induced electronic and structural response of a solute by combining an integral equation approach to include pressure effects on solvent structure with a quantum-chemical treatment of the solute within the embedded cluster reference interaction site model (EC-RISM) framework. Next, the solute’s response to compression is taken into account by introducing pressure-dependence into selected parameters of a well-established force field. In our proof-of-principle study, the full machinery is applied to N,N,N-trimethylamine-N-oxide (TMAO) in water being a potent osmolyte that counteracts pressure denaturation. EC-RISM theory is shown to describe well the charge redistribution upon compression of TMAO(aq) to 10 kbar, which is then embodied in force field molecular dynamics by pressure-dependent partial charges. The performance of the high pressure force field is assessed by comparing to experimental and ab initio molecular dynamics data. Beyond its broad usefulness for designing non-polarizable force fields for extreme thermodynamic conditions, a good description of the pressure-response of solutions is highly recommended when constructing and validating polarizable force fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

On the theory of electric double layer with explicit account of a polarizable co-solvent

We present a continuation of our theoretical research into the influence of co-solvent polarizability on a differential capacitance of the electric double layer. We formulate a modified Poisson-Boltzmann theory, using the formalism of density functional approach on the level of local density approximation taking into account the electrostatic interactions of ions and co-solvent molecules as well as their excluded volume. We derive the modified Poisson-Boltzmann equation, considering the three-component symmetric lattice gas model as a reference system and minimizing the grand thermodynamic potential with respect to the electrostatic potential. We apply present modified Poisson-Boltzmann equation to the electric double layer theory, showing that accounting for the excluded volume of co-solvent molecules and ions slightly changes the main result of our previous simplified theory. Namely, in the case of small co-solvent polarizability with its increase under the enough small surface potentials of electrode, the differential capacitance undergoes the significant growth. Oppositely, when the surface potential exceeds some threshold value (which is slightly smaller than the saturation potential), the increase in the co-solvent polarizability results in a differential capacitance decrease. However, when the co-solvent polarizability exceeds some threshold value, its increase generates a considerable enhancement of the differential capacitance in a wide range of surface potentials. We demonstrate that two qualitatively different behaviors of the differential capacitance are related to the depletion and adsorption of co-solvent molecules at the charged electrode. We show that an additive of the strongly polarizable co-solvent to an electrolyte solution can shift significantly the saturation potential in two qualitatively different manners. Namely, a small additive of strongly polarizable co-solvent results in a shift of saturation potential to higher surface potentials. On the contrary, a sufficiently large additive of co-solvent shifts the saturation potential to lower surface potentials. We obtain that an increase in the co-solvent polarizability makes the electrostatic potential profile longer-ranged. However, increase in the co-solvent concentration in the bulk leads to non-monotonic behavior of the electrostatic potential profile. An increase in the co-solvent concentration in the bulk at its sufficiently small values makes the electrostatic potential profile longer-ranged. Oppositely, when the co-solvent concentration in the bulk exceeds some threshold value, its further increase leads to decrease in electrostatic potential at all distances from the electrode.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

How strongly do hydrogen and water molecules stick to carbon nanomaterials?

The interaction strength of molecular hydrogen and water to carbon nanomaterials is relevant to, among many applications, hydrogen storage, water treatment, and water flow. Yet, accurate interaction energies for hydrogen and water with carbon nanotubes (CNTs) remain scarce despite the importance of having reliable benchmark data to inform experiments and to validate computational models. Here, benchmark fixed-node diffusion Monte Carlo (DMC) interaction energies are provided for hydrogen and water monomers inside and outside a typical zigzag CNT. The DMC interaction energies provide valuable insight into molecular interactions with CNTs in general and are also expected to be particularly relevant to gas uptake studies on CNTs. In addition, a selection of density functional theory (DFT) exchange-correlation (xc) functionals and force field potentials that ought to be suitable for these systems is compared. An unexpected variation is found in the performance of DFT van der Waals (vdW) models in particular. An analysis of the peculiar discrepancy between different vdW models indicates that medium-range correlation (at circa 3 to 5 Å) plays a key role inside CNTs and is poorly predicted by some vdW models. Using accurate reference information, this work reveals which xc functionals and force fields perform well for molecules interacting with CNTs. The results will be valuable to future work on these and related systems that involve molecules interacting with low-dimensional systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Accounting for Debye sheath expansion for proud Langmuir probes in magnetic confinement fusion plasmas

A Child-Langmuir law-based method for accounting for Debye sheath expansion while fitting the current-voltage I-V characteristic of proud Langmuir probes (electrodes that extend into the volume of the plasma) is described. For Langmuir probes of a typical size used in tokamak plasmas, these new estimates of electron temperature and ion saturation current density values decreased by up to 60% compared to methods that did not account for sheath expansion. Changes to the collection area are modeled using the Child-Langmuir law and effective expansion perimeter l p , and the model is thus referred to as the “perimeter sheath expansion method.” l p is determined solely from electrode geometry, so the method may be employed without prior measurement of the magnitude of the sheath expansion effects for a given Langmuir probe and can be used for electrodes of different geometries. This method correctly predicts the non-saturating ΔI/ΔV slope for cold, low-density plasmas where sheath-expansion effects are strong, as well as for hot plasmas where ΔI/ΔV ~ 0, though it is shown that the sheath can still significantly affect the collection area in these hot conditions. The perimeter sheath expansion method has several advantages compared to methods where the non-saturating current is fitted: (1) It is more resilient to scatter in the I-V characteristics observed in turbulent plasmas. (2) It is able to separate the contributions to the ΔI/ΔV slope from sheath expansion to that of the high energy electron tail in high Te conditions. (3) It calculates the change in the collection area due to the Debye sheath for conditions where ΔI/ΔV ~ 0 and for V = V f .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Chemical bond effects in classical site density functional theory of inhomogeneous molecular liquids

Intra-molecular interactions or chemical bonds represent one of the main distinguishing characteristic of molecular fluids. Development of accurate and practical methods to treat these effects is one of the long standing problems in classical site density functional theory (SDFT). One particular instance when these issues become particularly severe is the case of classical interactions potentials with auxiliary sites or dummy atoms. In this situation current SDFT implementations, such as three-dimensional reference interaction site model (3D-RISM), lead to nonphysical results. We re-examine this issue in this work using our recent reformulation of SDFT. We put forward a simple practical solution to this problem, and illustrate its utility for the case of spherical solutes in diatomic liquids.

Chuev, Gennady N.↗

DOE OSTI2020

Exact sum rules with approximate ground states

Electromagnetic and weak transitions tell us a great deal about the structure of atomic nuclei. Yet modeling transitions can be difficult: it is often easier to compute the ground state, if only as an approximation, than excited states. One alternative is through transition sum rules, in particular the non-energy-weighted and energy-weighted sum rules, which can be computed as expectation values of operators. We investigate by computing sum rules for a variety of nuclei, comparing the numerically exact full configuration-interaction shell model, as a reference, to Hartree-Fock, projected Hartree-Fock, and the nucleon pair approximation. Furthermore, these approximations yield reasonable agreement, which we explain by prior work on the systematics of transition moments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2016

Cluster abundance in chameleon f ( R ) gravity I: toward an accurate halo mass function prediction

We refine the mass and environment dependent spherical collapse model of chameleon f ( R ) gravity by calibrating a phenomenological correction inspired by the parameterized post-Friedmann framework against high-resolution N -body simulations. We employ our method to predict the corresponding modified halo mass function, and provide fitting formulas to calculate the enhancement of the f ( R ) halo abundance with respect to that of General Relativity (GR) within a precision of ∼< 5% from the results obtained in the simulations. Similar accuracy can be achieved for the full f ( R ) mass function on the condition that the modeling of the reference GR abundance of halos is accurate at the percent level. We use our fits to forecast constraints on the additional scalar degree of freedom of the theory, finding that upper bounds competitive with current Solar System tests are within reach of cluster number count analyses from ongoing and upcoming surveys at much larger scales. Importantly, the flexibility of our method allows also for this to be applied to other scalar-tensor theories characterized by a mass and environment dependent spherical collapse.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Common origin of warm dark matter and dark radiation

We consider a cosmological scenario where a relativistic particle and a stable massive particle are simultaneously produced from the decay of a late-decaying particle after Big-Bang Nucleosynthesis but before matter-radiation equality. The relativistic and massive particles behave as dark radiation and warm dark matter, respectively. Due to a common origin, the warmness and abundances are closely related. We refer to the models that lead to such a scenario as Common Origin of Warm and Relativistic Decay Products (COWaRD). We show that COWaRD predicts a correlation between the amount of dark radiation and suppression of the large scale structure, which can be tested in future precision cosmology observations. Furthermore, we demonstrate that COWaRD is realized, as an example, in a class of supersymmetric axion models and that future observations by the next generation Cosmic Microwave Background, Large Scale Structure, and 21-cm surveys can reveal the structure of the theory.

cosmology of theories beyond the SM↗

DOE OSTI2018

Automated map sharpening by maximization of detail and connectivity

An algorithm for automatic map sharpening is presented that is based on optimization of the detail and connectivity of the sharpened map. The detail in the map is reflected in the surface area of an iso-contour surface that contains a fixed fraction of the volume of the map, where a map with high level of detail has a high surface area. The connectivity of the sharpened map is reflected in the number of connected regions defined by the same iso-contour surfaces, where a map with high connectivity has a small number of connected regions. By combining these two measures in a metric termed the `adjusted surface area', map quality can be evaluated in an automated fashion. This metric was used to choose optimal map-sharpening parameters without reference to a model or other interpretations of the map. Map sharpening by optimization of the adjusted surface area can be carried out for a map as a whole or it can be carried out locally, yielding a locally sharpened map. To evaluate the performance of various approaches, a simple metric based on map–model correlation that can reproduce visual choices of optimally sharpened maps was used. The map–model correlation is calculated using a model with B factors (atomic displacement factors; ADPs) set to zero. This model-based metric was used to evaluate map sharpening and to evaluate map-sharpening approaches, and it was found that optimization of the adjusted surface area can be an effective tool for map sharpening.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Evaluation of Composite Structural Materials for Heliostat Cost Reduction

Structures manufactured from steel comprise up to 40% of a CSP heliostat's cost. Composite structures represent a potential opportunity to reduce this cost. A reference heliostat structural model has been created with a reflector area of 25m 2 . The design, constructed of low-carbon steel, pro-vides baseline deflection and stiffness under a 21 m/s operating wind speeds. Wind loads on the tracker structure are determined for both operating and stow conditions. An established roster of suitable metal alternative materials is considered including: glass, basalt, and carbon reinforced polymer (GFRP, BFRP, and CFRP respectively). Three heliostat components are investigated: the pylon, torque tube, and the purlin-strut assembly. Composite material properties are substituted for those of steel, and the beams are re-sized to match the original steel components' deflection under given wind loads. Weight and cost changes resulting from this resizing are evaluated. It is found that GFRP and BFRP represent a 3X–6X cost premium for the same operating deflection character-istics as steel across all three investigated component classes; with weight reduction only achieved for the purlin-strut assembly. While CFRP components can achieve approximately 25–75% weight savings depending on the application, this comes with a 9X–14X cost increase over the steel base-line for tube-type structures and roughly 5X cost increase when replacing c-channel structures. This work does not rule out the possibility of cost savings when the heliostat design and kinematics to take advantage of composites' specific properties.

14 SOLAR ENERGY↗

DOE OSTI2014

TH-C-12A-11: Target Correlation of a 3D Surface Surrogate for Left Breast Irradiation Using the Respiratory-Gated Deep-Inspiration Breath-Hold Technique

Purpose: To evaluate the use of 3D optical surface imaging as a new surrogate for respiratory motion gated deep-inspiration breath-hold (DIBH) technique for left breast cancer patients. Methods: Patients with left-sided breast cancer after lumpectomy or mastectomy were selected as candidates for DIBH technique for their external beam radiation therapy. Treatment plans were created on both free breathing (FB) and DIBH CTs to determine whether DIBH was beneficial in reducing heart doses. The Real-time Position Management (RPM) system was used to acquire patient's breathing trace during DIBH CT acquisition and treatment delivery. The reference 3D surface models from FB and DIBH CTs were generated and transferred to the “AlignRT” system for patient positioning and real-time treatment monitoring. MV Cine images were acquired for each beam as quality assurance for intra-fractional position verification. The chest wall excursions measured on these images were used to define the actual target position during treatment, and to investigate the accuracy and reproducibility of RPM and AlignRT. Results: Reduction in heart dose can be achieved for left-sided breast patients using DIBH. Results showed that RPM has poor correlation with target position, as determined by the MV Cine imaging. This indicates that RPM may not be an adequate surrogate in defining the breath-hold level when used alone. Alternatively, the AlignRT surface imaging demonstrated a better correlation with the actual CW excursion during DIBH. Both the vertical and magnitude real-time deltas (RTDs) reported by AlignRT can be used as the gating parameter, with a recommend threshold of ±3 mm and 5 mm, respectively. Conclusion: 3D optical surface imaging serves as a superior target surrogate for the left breast treatment when compared to RPM. Working together with the realtime MV Cine imaging, they ensure accurate patient setup and dose delivery, while minimizing the imaging dose to patients.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

pnnl/90.1-PRM_ECB-ReviewManual

The 90.1 Energy Cost Budget and Performance Rating Method Submittal Review Manual (the Manual) is a comprehensive reference for reviewing modeling based energy code compliance submittals.

Xu, Weili [PNNL]↗

DOE OSTI2020

SMART Mobility. Advanced Fueling Infrastructure Capstone Report

The U.S. Department of Energy’s Systems and Modeling for Accelerated Research in Transportation (SMART) Mobility Consortium is a multiyear, multi-laboratory collaborative, managed by the Energy Efficient Mobility Systems Program of the Office of Energy Efficiency and Renewable Energy, Vehicle Technologies Office, dedicated to further understanding the energy implications and opportunities of advanced mobility technologies and services. The first three-year research phase of SMART Mobility occurred from 2017 through 2019 and included five research pillars: Connected and Automated Vehicles, Mobility Decision Science, Multi-Modal Freight, Urban Science, and Advanced Fueling Infrastructure. A sixth research thrust integrated aspects of all five pillars to develop a SMART Mobility Modeling Workflow to evaluate new transportation technologies and services at scale. This report summarizes the work of the Advanced Fueling Infrastructure Pillar. This Pillar investigated the charging infrastructure needs of electric ride-hailing and car-sharing vehicles, automated shuttle buses, and freight-delivery truck fleets. For information about the other Pillars and about the SMART Mobility Modeling Workflow, please refer to the relevant Pillar’s Capstone Report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2018

A Composition-Based Approach for Predicting Hanford Tank Waste Liquid Viscosity

Liquid viscosity is an important Hanford waste physical property that governs, in part, the performance of key Hanford flowsheet unit operations and processes, including waste filtration, ion exchange, pumping and transport, and mixing. Evaluation and optimization of these processes rely on liquid viscosity correlations subject to substantial error (at least ±60%) that can lead to overly conservative, and sometimes physically unachievable, limits for validation and operation of these unit operations. The current liquid viscosity correlations used for the Hanford flowsheet attempt to capture the effects of specific waste chemistry (namely, the speciation of various salts) on liquid viscosity using liquid density alone. It is postulated that the ±60% or greater error associated with these “density-based” liquid viscosity predictions is driven in part by the inability of waste density to capture the full impact of waste chemistry on liquid viscosity, and that improved prediction accuracy could be achieved by implementing a composition-based approach to evaluating liquid viscosity for Hanford wastes. Recent publications, namely those of Lencka et al. (1998), Hu (2004), and Laliberté (2007a), provide correlations and comprehensive, consolidated databases of single-salt solution viscosity that enable composition-based predictions for Hanford-type liquids (i.e., concentrated multi-component aqueous electrolyte solutions). In the current study, the viscosities of Hanford wastes and waste simulants were predicted using a composition-based model developed by Laliberté (2007a) (hereafter referred to as the “Laliberté viscosity model”) to evaluate its accuracy relative to the current density-based models used to predict waste viscosity for engineering applications at the Hanford Site. A viscosity database that captured the composition and viscosity of wastes and waste simulants was developed to support this evaluation. To facilitate application of the Laliberté viscosity model, waste and waste simulant compositions were converted to an equivalent salt recipe composed of salts with known viscosity coefficients using a number of simplifying assumptions that have the potential to impact the fidelity of the model prediction. Despite these simplifications, the Laliberté viscosity model was able to predict viscosity more accurately than the density-based correlations currently available for Hanford flowsheet use, such as that proposed in RPP-RPT- 51652. The use of composition-based viscosity models, like that proposed by Laliberté (2007a), shows promise for leading to improved prediction of liquid waste viscosity. As implemented in this study, application of the Laliberté viscosity model could lead to a decrease in the uncertainty associated with liquid viscosity prediction by as much as a factor of two, allowing for reduction in unit operation design conservatism and for the more realistic selection of test simulants for flowsheet process evaluation and optimization, all of which could potentially reduce flowsheet design margins, lower costs and reduce duration of the Hanford mission, and reduce Hanford flowsheet risks.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗